Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

Adaptive Optimization for System Performance and Combined Bernstein Polynomial, Optimal Reciprocal Collision Avoidance, Differential Dynamic Programming for Trajectory Replanning and Collision Avoidance for UAM Vehicles

The emerging urban air mobility (UAM) sector in aerospace is driving development of unconventional multi-modal vehicle configurations and autonomous flight. The combination of multi-modal vehicle dynamics, complex environment, requirements to deal with flight contingencies in an efficient and safe manner, as well as necessity for precise trajectory following and performance, are the driving influence behind adaptive optimization for system performance. We are interested in trajectory optimization algorithm that would system parameter estimation and identifying the optimal switching time between modes of hybrid dynamical systems. This presentation discusses a parameterized optimal control trajectory optimization algorithm that is an extended and generalized version of Differential Dynamic Programming (DDP), titled Parameterized Differential Dynamic Programming (PDDP). DDP is an efficient trajectory optimization algorithm relying on second order approximations of a system’s dynamics and cost function and has recently been applied to optimize systems with time invariant parameters. Experiments are presented applying PDDP to solve model predictive control (MPC) and moving horizon estimation (MHE) tasks simultaneously. In particular, PDDP is used to determine the optimal transition point between flight regimes of a complex urban air mobility (UAM) class vehicle exhibiting multiple phases of flight and to identify and compensate for actuation faults.

optimization↗

Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center has piloted an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project was the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the NASA Lunar Architecture Team and may allow for increased science productivity and return on investment than with only the Surface Habitat. The Forge team developed two concepts, a stationary laboratory and a mobile laboratory. Human factors evaluations have been conducted to determine the acceptability of the resulting design concepts. The human factors evaluations conducted include a card sort evaluation of the laboratories by science discipline and the System Usability Scale (SUS). The card sort is focused subject matter expert expectations for the relative positioning of the laboratory functional areas, determining whether a given lab area should be adjacent to, or separated from, another. The SUS will focus on layout efficiency. Both will use human factors personnel and both life and physical scientists as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of the card sort and SUS evaluations and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of human factors rapid assessment processes in early design trades will also be provided.

Human Factors↗

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗

Design of a Thrust Vector Servoelastic Filter for SLS Block 2 Booster Thrust Vector Control System

This paper presents a method to improve the thrust vector servoelastic (TVSE) response of a launch vehicle. TVSE is the coupling between the Thrust Vector Control (TVC) system nozzle dynamics and the vehicle flexible modes. A poorly phased TVC response, i.e. in phase with the structural velocity, can lead to instability of the flexible modes due to the energy injection from the thrust. For large boosters with heavy nozzles, TVSE is more likely to cause instability. This instability is part of the plant dynamics and may not be readily detected by the flight control system. In this paper, TVSE is treated as a disturbance rejection problem with the inertial force of the flexible motion viewed as an external disturbance to the TVC actuator control system. By effectively rejecting this disturbance, the engine can be decoupled from the flexible dynamics thus improving TVSE stability. The actuator disturbance transfer function is modified by inserting a filter into the servo-loop control system. The filter aims to reduce the magnitude and introduce phase lead of the disturbance transfer function over the TVSE sensitive frequency range while maintaining adequate stability margins of the actuator control system. A prefilter is also designed to improve command tracking performance. The filter and prefilter resulting from this method are usually of low order and can be included to an existing actuator controller if the controller is programmable. This method was demonstrated for a hydraulic Simplex actuator model. Its applications to the Space Launch System (SLS) Block 1 and Block 2 TVC actuator are also described. The Block 2 actuator controller was modified based on the results of this method. The tradeoff between TVSE stability and the servo-loop stability can be balanced by adjusting the amount of phase lead introduced by the TVSE filter. This method provides a useful way to solve the TVSE instability issue without costly and time-consuming hardware solutions. This method is based on the basic control principles of loop shaping and disturbance rejection and has been proven to be effective under proper circumstances.

Thrust vector control↗

Design of a Thruster Vector Servo-Elastic Filter

This presents a method to improve the thruster vector servo elasticity (TVSE) stability. TVSE is the coupling between the TVC system engine dynamics and the vehicle flexible modes. It can lead to instability of the flexible modes due to the energy injection from the thrust. For large boosters with heavy TVC engine, TVSE is more likely to cause instability. This instability is part of the plant dynamics and is not detected by the flight control. In this paper, the TVSE is treated as a disturbance rejection problem with the inertial force of the flexible motion viewed as an external disturbance to the TVC actuator control system. By effectively rejecting this disturbance, the engine can be decoupled from the flexible dynamics thus improving TVSE stability. The actuator disturbance transfer function is modified by inserting a TVSE filter into the control system. The filter aims to reduce the magnitude and introduce phase lead of the disturbance transfer function over the TVSE sensitive frequency range while maintaining adequate stability margins of the actuator control system. A prefilter is also designed for good command tracking. The filter and prefilter resulting from this method are usually of low order and can be included to the existing actuator controller if the controller is programmable. This method was demonstrated for a hydraulic Simplex actuator model. Its applications to SLS B1 and B2 TVC actuator are also described. B2 actuator controller is modified based on this method. The tradeoff between TVSE stability and the control system stability can be balanced by adjusting the amount of phase lead introduced by the TVSE filter. This method provides a useful way to solve the TVSE instability issue without costly and time-consuming a TVC system redesign. This method is based on the basic control principles of loop shaping and disturbance rejection and has been proven to be effective under proper circum-stances.

Thrust vector control↗

On recursive least-squares filtering algorithms and implementations

In many real-time signal processing applications, fast and numerically stable algorithms for solving least-squares problems are necessary and important. In particular, under non-stationary conditions, these algorithms must be able to adapt themselves to reflect the changes in the system and take appropriate adjustments to achieve optimum performances. Among existing algorithms, the QR-decomposition (QRD)-based recursive least-squares (RLS) methods have been shown to be useful and effective for adaptive signal processing. In order to increase the speed of processing and achieve high throughput rate, many algorithms are being vectorized and/or pipelined to facilitate high degrees of parallelism. A time-recursive formulation of RLS filtering employing block QRD will be considered first. Several methods, including a new non-continuous windowing scheme based on selectively rejecting contaminated data, were investigated for adaptive processing. Based on systolic triarrays, many other forms of systolic arrays are shown to be capable of implementing different algorithms. Various updating and downdating systolic algorithms and architectures for RLS filtering are examined and compared in details, which include Householder reflector, Gram-Schmidt procedure, and Givens rotation. A unified approach encompassing existing square-root-free algorithms is also proposed. For the sinusoidal spectrum estimation problem, a judicious method of separating the noise from the signal is of great interest. Various truncated QR methods are proposed for this purpose and compared to the truncated SVD method. Computer simulations provided for detailed comparisons show the effectiveness of these methods. This thesis deals with fundamental issues of numerical stability, computational efficiency, adaptivity, and VLSI implementation for the RLS filtering problems. In all, various new and modified algorithms and architectures are proposed and analyzed; the significance of any of the new method depends crucially on specific application.

Hsieh, Shih-Fu↗

Evaluation of a transfinite element numerical solution method for nonlinear heat transfer problems

Laplace transform techniques have been widely used to solve linear, transient field problems. A transform-based algorithm enables calculation of the response at selected times of interest without the need for stepping in time as required by conventional time integration schemes. The elimination of time stepping can substantially reduce computer time when transform techniques are implemented in a numerical finite element program. The coupling of transform techniques with spatial discretization techniques such as the finite element method has resulted in what are known as transfinite element methods. Recently attempts have been made to extend the transfinite element method to solve nonlinear, transient field problems. This paper examines the theoretical basis and numerical implementation of one such algorithm, applied to nonlinear heat transfer problems. The problem is linearized and solved by requiring a numerical iteration at selected times of interest. While shown to be acceptable for weakly nonlinear problems, this algorithm is ineffective as a general nonlinear solution method.

Cerro, J. A.↗

Diffusive and radiative effects on vaporization times of drops in film boiling

Diffusive and radiative effects are incorporated into an analysis for the vaporization time of drops in film boiling. The momentum, energy, and continuity equations are solved with some appropriate simplifications so as to obtain a simple closed form solution for the overall film boiling heat transfer coefficient. Next, a theoretical expression for the droplet vaporization time is developed and compared to the measured vaporization times of water droplets vaporizing into air, argon, nitrogen, and helium. The agreement between experiment and theory is good. Under the helium blanket, the diffusive evaporative component is significant in comparison to the film boiling component.

Baumeister, K. J.↗

Multi-Robot Assembly Scheduling for the Lunar Crater Radio Telescope on the Far-Side of the Moon

The Lunar Crater Radio Telescope (LCRT) is a pro- posed ultra-long-wavelength radio telescope to be constructed on the far side of the moon. The proposed telescope will be constructed by deploying a 1km wire mesh in a 3-5km crater using a team of wall-climbing DuAxel robots. In this work, we consider the problem of generating minimum-time assembly sequences for LCRT, using realistic models of travel speed and lighting. Specifically, we pose the assembly sequencing problem as a mixed-integer linear program (MILP), which we solve to global optimality using commercial solvers. We present methods for modeling time-varying travel and assembly times, based on variable lighting conditions (including crater shadowing), and show how such time-varying parameters can be incorporated into the MILP. Finally, we present numerical studies of our method, showing how makespan varies with the number of assembly robots.

Schwager, Mac↗

Application of the space-time conservation element and solution element method to two-dimensional advection-diffusion problems

The existing 2-D alpha-mu scheme and alpha-epsilon scheme based on the method of space-time conservation element and solution element, which were constructed for solving the linear 2-D unsteady advection-diffusion equation and unsteady advection equation, respectively, are tested. Also, the alpha-epsilon scheme is modified to become the V-E scheme for solving the nonlinear 2-D inviscid Burgers equation. Numerical solutions of six test problems are presented in comparison with their exact solutions or numerical solutions obtained by traditional finite-difference or finite-element methods. It is demonstrated that the 2-D alpha-mu, alpha-epsilon, and nu-epsilon schemes can be used to obtain numerical results which are more accurate than those based on some of the traditional methods but without using any artificial tuning in the computation. Similar to the previous 1-D test problems, the high accuracy and simplicity features of the space-time conservation element and solution element method have been revealed again in the present 2-D test results.

Wang, Xiao-Yen↗

WHAMII - An enumeration and insertion procedure with binomial bounds for the stochastic time-constrained traveling salesman problem

This paper presents an algorithm (WHAMII) designed to solve the Artificial Intelligence Design Challenge at the 1987 AIAA Guidance, Navigation and Control Conference. The problem under consideration is a stochastic generalization of the traveling salesman problem in which travel costs can incur a penalty with a given probability. The variability in travel costs leads to a probability constraint with respect to violating the budget allocation. Given the small size of the problem (eleven cities), an approach is considered that combines partial tour enumeration with a heuristic city insertion procedure. For computational efficiency during both the enumeration and insertion procedures, precalculated binomial probabilities are used to determine an upper bound on the actual probability of violating the budget constraint for each tour. The actual probability is calculated for the final best tour, and additional insertions are attempted until the actual probability exceeds the bound.

Dahl, Roy W.↗

Acceleration of FDNS flow simulations using initial flowfields generated with a parabolized Navier-Stokes method

The capability of computational fluid dynamics (CFD) analysis to predict complex flowfields has been greatly advanced by the widespread use of flow simulation programs based on the Navier-Stokes (NS) equations. The flow physics are theoretically well represented, and with proper care the numerical solution should not introduce appreciable uncertainties. However, the computational cost of a typical simulation in terms of both memory and execution time are large by current standards. Therefore, the efficiency of the numerical algorithm in solving the set of model equations is one factor determining the usefulness of CFD tools. The objective of the present study is to reduce the execution time of the FDNS flow simulation code by using the parabolized Navier-Stokes (PNS) equations to provide a 'good' starting condition. The technique is not universal, however the PNS model can be applied to convection dominated flows with moderate deflection in geometries that are free of large obstructions.

Tenpas, Peter W.↗

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations was solved for in the steady state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time step, while the turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positively. Low Reynolds number modifications to the original two equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗

A time-dependent theoretical model of the polar wind Preliminary results

The coupled, time-dependent continuity, momentum, and energy equations of a two-ion O(+) and H(+) quasi-neutral plasma were solved in order to extend understanding of polar wind behavior. This numerical code allows studies of the time dependent behavior of polar wind-type flows into and out of the ionosphere. Initial studies indicate that the typical time constants for electron and ion temperature changes are of the order of minutes and tens of minutes, respectively. The response time of the minor high altitude ion O(+) is less than an hour, whereas that of the major ion, H(+), is many hours. The initial test runs also demonstrate the fact that temporary supersonic flows of both O(+) and H(+) are possible, especially in the presence of significant ion heating.

Gombosi, T. I.↗

Acquisition time for sideband modulated phase lock loops.

In this paper a definition for acquisition of a sideband modulated carrier is proposed for the phase lock loop, and the acquisition differential equation solved for a region of the modulation parameter space to determine the expected value of acquisition time. Principally the case of modulation by a single unmodulated subcarrier is considered. The numerical analysis is made for the noiseless case.

Gilbert, A. L.↗

Constraints on the Martian cratering rate imposed by the SNC meteorites and Vallis Marineris layered deposits

Following two independent lines of evidence -- estimates of the age and formation time of a portion of the Martian geologic column exposed in the layered deposits and the crystallization and ejection ages of the SNC meteorites -- it appears that the Martian cratering rate must be double the lunar rate or even higher. This means models such as NHII or NHIII (Neukum and Hiller models II and III), which estimate the Martian cratering rate as being several times lunar are probably far closer to reality on Mars than lunar rates. The effect of such a shift is profound: Mars is transformed from a rather Moon-like place into a planet with vigorous dynamics, multiple large impacts, erosion, floods, and volcanism throughout its history. A strong shift upward in cratering rates on Mars apparently solves some glaring problems; however, it creates others. The period of time during which Earth-like atmospheric conditions existed, the liquid water era on Mars, persists in NHIII up to only 0.5 b.y. ago. Scenarios of extended Earth-like conditions on Mars have been discounted in the past because they would have removed many of the craters from the early bombardment era found in the south. It does appear that some process of crater removal was quite vigorous in the north during Mars' past. Evidence exists that the northern plains may have been the home of long-lived seas or perhaps even a paleo-ocean, so models exist for highly localized destruction of craters in the north. However, the question of how the ancient crater population could be preserved in the south under a long liquid-water era found in any high-cratering-rate models is a serious question that must be addressed. It does appear to be a higher-order problem because it involves low-energy dynamics acting in localized areas, i.e., erosion of craters in the south of Mars, whereas the two problems with the low-cratering-rate models involve high-energy events acting over large areas: the formation of the Vallis Marineris, the SNC ejecting impacts, and the global atmospheric pressure and temperature conditions that allow liquid water to exist as a robust entity anywhere on the Martian surface. In any case, it appears Mars is a more complex and dynamic planet than previously supposed. It has canyons dating from the middle to late period of its history that contain apparent lake sediments bedded deeper than most sediments on Earth. Recent multiple, violent impacts on Mars have apparently provided us with multiple random samples of its surface that all crystallized less than 1.5 b.y. ago. These things cannot be accommodated in our present cratering chronologies of Mars, based on 1x lunar cratering rates, without great difficulties. These difficulties suggest that a new chronology, probably based on NHII or even NHIII, should be adopted; this new chronology will provide us with a new view of Mars as a dynamic planet of rich history.

Brandenburg, J. E.↗

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been solved for in the steady-state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time-step, while the turbulence equations are advanced in a point-implicit scheme with a time-step which guarantees stability and positivity. Low Reynolds number modifications to the original two-equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗