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 1,153 records · Page 64

Computational complexities and storage requirements of some Riccati equation solvers

The linear optimal control problem of an nth-order time-invariant dynamic system with a quadratic performance functional is usually solved by the Hamilton-Jacobi approach. This leads to the solution of the differential matrix Riccati equation with a terminal condition. The bulk of the computation for the optimal control problem is related to the solution of this equation. There are various algorithms in the literature for solving the matrix Riccati equation. However, computational complexities and storage requirements as a function of numbers of state variables, control variables, and sensors are not available for all these algorithms. In this work, the computational complexities and storage requirements for some of these algorithms are given. These expressions show the immensity of the computational requirements of the algorithms in solving the Riccati equation for large-order systems such as the control of highly flexible space structures. The expressions are also needed to compute the speedup and efficiency of any implementation of these algorithms on concurrent machines.

Utku, Senol↗

On the parallel solution of parabolic equations

Parallel algorithms for the solution of linear parabolic problems are proposed. The first of these methods is based on using polynomial approximation to the exponential. It does not require solving any linear systems and is highly parallelizable. The two other methods proposed are based on Pade and Chebyshev approximations to the matrix exponential. The parallelization of these methods is achieved by using partial fraction decomposition techniques to solve the resulting systems and thus offers the potential for increased time parallelism in time dependent problems. Experimental results from the Alliant FX/8 and the Cray Y-MP/832 vector multiprocessors are also presented.

Gallopoulos, E.↗

Cracking of coated materials under transient thermal stresses

The crack problem for a relatively thin layer bonded to a very thick substrate under thermal shock conditions is considered. The effect of surface cooling rate is studied by assuming the temperature boundary condition to be a ramp function. Among the crack geometries considered are the edge crack in the coating layer, the broken layer, the edge crack going through the interface, the undercoat crack in the substrate and the embedded crack crossing the interface. The primary calculated quantity is the stress intensity factor at various singular points and the main variables are the relative sizes and locations of cracks, the time, and the duration of the cooling ramp. The problem is solved and rather extensive results are given for two material pairs, namely a stainless steel layer welded on a ferritic medium and a ceramic coating on a steel substrate.

Rizk, A. A.↗

Parallel solution of pentadiagonal systems using generalized odd-even elimination

A method for the solution of pentadiagonal systems of linear equations is presented. The method is a generalization of ordinary odd-even elimination used for tridiagonal systems. Using n processors, an n x n pentadiagonal system can be solved using the new method (generalized odd-even elimination) in time proportional to log(2) n.

Levit, Creon↗

Rapid in-plane maneuvering of the flexible orbiting SCOLE

The rapid two-dimensional slewing and vibrational control of the unsymmetrical flexible SCOLE (Spacecraft Control Laboratory Experiment) with multi-bounded controls has been considered. Pontryagin's Maximum Principle has been applied to the nonlinear equations of the system to derive the necessary conditions for the optimal control. The resulting two-point boundary-value problem is then solved by using the quasilinearization technique, and the near-minimum time is obtained by sequentially shortening the slewing time until the controls are near the bang-bang type. The trade-off between the minimum time and the minimum flexible amplitude requirements has been discussed. The numerical results show that the responses of the nonlinear system are significantly different from those of the linearized system for rapid slewing.

Bainum, Peter M.↗

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗

Progress in the development of advanced computational models for the analysis of generalized supersonic jet flowfields

This paper describes two three-dimensional computational codes (PARCH and CRAFT) for computing a generalized supersonic jet flowfields with subsonic external streams, which solve the Reynolds-averaged Navier-Stokes equations in a time asymptotic manner. The application of these codes to an imperfectly expanded jet flowfield is demonstrated. Special attention is given to numerical issues and to turbulence modeling. The application of the PARCH and CRAFT codes to other related flowfield problems, such as gas turbine and nozzle flows, and to generalized missile-plume flowfields is discussed.

Dash, S. M.↗

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗

The numerical simulation of a high-speed axial flow compressor

The advancement of high-speed axial-flow multistage compressors is impeded by a lack of detailed flow-field information. Recent development in compressor flow modeling and numerical simulation have the potential to provide needed information in a timely manner. The development of a computer program is described to solve the viscous form of the average-passage equation system for multistage turbomachinery. Programming issues such as in-core versus out-of-core data storage and CPU utilization (parallelization, vectorization, and chaining) are addressed. Code performance is evaluated through the simulation of the first four stages of a five-stage, high-speed, axial-flow compressor. The second part addresses the flow physics which can be obtained from the numerical simulation. In particular, an examination of the endwall flow structure is made, and its impact on blockage distribution assessed.

Mulac, Richard A.↗

Advanced design for orbital debris removal in support of solar system exploration

The development of an Autonomous Space Processor for Orbital Debris (ASPOD) is the ultimate goal. The craft will process, in situ, orbital debris using resources available in low Earth orbit (LEO). The serious problem of orbital debris is briefly described and the nature of the large debris population is outlined. This year, focus was on development of a versatile robotic manipulator to augment an existing robotic arm; incorporation of remote operation of robotic arms; and formulation of optimal (time and energy) trajectory planning algorithms for coordinating robotic arms. The mechanical design of the new arm is described in detail. The versatile work envelope is explained showing the flexibility of the new design. Several telemetry communication systems are described which will enable the remote operation of the robotic arms. The trajectory planning algorithms are fully developed for both the time-optimal and energy-optimal problem. The optimal problem is solved using phase plane techniques while the energy optimal problem is solved using dynamics programming.

Source record↗

Massively parallel support for a case-based planning system

Case-based planning (CBP), a kind of case-based reasoning, is a technique in which previously generated plans (cases) are stored in memory and can be reused to solve similar planning problems in the future. CBP can save considerable time over generative planning, in which a new plan is produced from scratch. CBP thus offers a potential (heuristic) mechanism for handling intractable problems. One drawback of CBP systems has been the need for a highly structured memory to reduce retrieval times. This approach requires significant domain engineering and complex memory indexing schemes to make these planners efficient. In contrast, our CBP system, CaPER, uses a massively parallel frame-based AI language (PARKA) and can do extremely fast retrieval of complex cases from a large, unindexed memory. The ability to do fast, frequent retrievals has many advantages: indexing is unnecessary; very large case bases can be used; memory can be probed in numerous alternate ways; and queries can be made at several levels, allowing more specific retrieval of stored plans that better fit the target problem with less adaptation. In this paper we describe CaPER's case retrieval techniques and some experimental results showing its good performance, even on large case bases.

Kettler, Brian P.↗

An effective multigrid method for high-speed flows

The use is considered of a multigrid method with central differencing to solve the Navier-Stokes equations for high speed flows. The time dependent form of the equations is integrated with a Runge-Kutta scheme accelerated by local time stepping and variable coefficient implicit residual smoothing. Of particular importance are the details of the numerical dissipation formulation, especially the switch between the second and fourth difference terms. Solutions are given for 2-D laminar flow over a circular cylinder and a 15 deg compression ramp.

Swanson, R. C.↗

Block-structured grids for complex aerodynamic configurations: Current status

The status of CFD methods based on the use of block-structured grids for analyzing viscous flows over complex configurations is examined. The objective of the present study is to make a realistic assessment of the usability of such grids for routine computations typically encountered in the aerospace industry. It is recognized at the very outset that the total turnaround time, from the moment the configuration is identified until the computational results have been obtained and postprocessed, is more important than just the computational time. Pertinent examples will be cited to demonstrate the feasibility of solving flow over practical configurations of current interest on block-structured grids.

Vatsa, Veer N.↗

Dominant takeover regimes for genetic algorithms

The genetic algorithm (GA) is a machine-based optimization routine which connects evolutionary learning to natural genetic laws. The present work addresses the problem of obtaining the dominant takeover regimes in the GA dynamics. Estimated GA run times are computed for slow and fast convergence in the limits of high and low fitness ratios. Using Euler's device for obtaining partial sums in closed forms, the result relaxes the previously held requirements for long time limits. Analytical solution reveal that appropriately accelerated regimes can mark the ascendancy of the most fit solution. In virtually all cases, the weak (logarithmic) dependence of convergence time on problem size demonstrates the potential for the GA to solve large N-P complete problems.

Noever, David↗

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗

CFD Analysis of Hypersonic Flowfields With Surface Thermochemistry and Ablation

In the past forty years much progress has been made in computational methods applied to the solution of problems in spacecraft hypervelocity flow and heat transfer. Although the basic thermochemical and physical modeling techniques have changed little in this time, several orders of magnitude increase in the speed of numerically solving the Navier-Stokes and associated energy equations have been achieved. The extent to which this computational power can be applied to the design of spacecraft heat shields is dependent on the proper coupling of the external flow equations to the boundary conditions and governing equations representing the thermal protection system in-depth conduction, pyrolysis and surface ablation phenomena. A discussion of the techniques used to do this in past problems as well as the current state-of-art is provided. Specific examples, including past missions such as Galileo, together with the more recent case studies of ESA/Rosetta Sample Comet Return, Mars Pathfinder and X-33 will be discussed. Modeling assumptions, design approach and computational methods and results are presented.

Henline, W. D.↗

Numerical Solutions of Supersonic and Hypersonic Laminar Compression Corner Flows

An efficient time-splitting, second-order accurate, numerical scheme is used to solve the complete Navier-Stokes equations for supersonic and hypersonic laminar flow over a two-dimensional compression corner. A fine, exponentially stretched mesh spacing is used in the region near the wall for resolving the viscous layer. Good agreement is obtained between the present computed results and experimental measurement for a Mach number of 14.1 and a Reynolds number of 1.04 x 10(exp 5) with wedge angles of 15 deg, 18 deg, and 24 deg. The details of the pressure variation across the boundary layer are given, and a correlation between the leading edge shock and the peaks in surface pressure and heat transfer is observed.

Hung, C. M.↗

Simulation of Unsteady Hypersonic Combustion Around Projectiles in an Expansion Tube

The temporal evolution of combustion flowfields established by the interaction between wedge-shaped bodies and explosive hydrogen-oxygen-nitrogen mixtures accelerated to hypersonic speeds in an expansion tube is investigated. The analysis is carried out using a fully implicit, time-accurate, computational fluid dynamics code that we developed recently for solving the Navier-Stokes equations for a chemically reacting gas mixture. The numerical results are compared with experimental data from the Stanford University expansion tube for two different gas mixtures at Mach numbers of 4.2 and 5.2. The experimental work showed that flow unstart occurred for the Mach 4.2 cases. These results are reproduced by our numerical simulations and, more significantly, the causes for unstart are explained. For the Mach 5.2 mixtures, the experiments and numerical simulations both produced stable combustion. However, the computations indicate that in one case the experimental data were obtained during the transient phase of the flow; that is, before steady state had been attained.

Yungster, S.↗