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,549 records · Page 86

Application of Navier-Stokes analysis to stall flutter

A solution procedure was developed to investigate the two-dimensional, one- or two-dimensional flutter characteristics of arbitrary airfoils. This procedure requires a simultaneous integration in time of the solid and fluid equations of motion. The fluid equations of motion are the unsteady compressible Navier-Stokes equations, solved in a body-fitted moving coordinate system using an approximate factorization scheme. The solid equations of motion are integrated in time using an Euler implicit scheme. Flutter is said to occur if small disturbances imposed on the airfoil attitude lead to divergent oscillatory motions at subsequent times. The flutter characteristics of airfoils in subsonic speed at high angles of attack and airfoils in high subsonic and transonic speeds at low angles of attack are investigated. The stall flutter characteristics are also predicted using the same procedure.

Wu, J. C.↗

A computational study of heat transfer for laminar hypervelocity flows

Computational Fluid Dynamics (CFD) codes are routinely used to predict the flowfield and the heating environment around complex reentry configurations. At hypervelocities, where the velocity is greater than 3 km/sec, the AFWAL version of the blunt body code predicts the correct surface pressure distributions but underpredicts laminar wall heat fluxes. This study was performed to determine the reasons for the underprediction. The computer code chosen solves thin-layer Navier-Stokes equations in a time-asymptotic manner and assumes a constant isentropic exponent. Flowfields around a spherical configuration at various entry velocities are computed. The computed pressure distributions agree well with the tabulated, inviscid results of Lyubimov and Rusanov for entry velocities ranging from 0.6 to 5.92 km/sec. At hypervelocities, the calculated stagnation point heat transfer rates were lower by roughly fifty percent when compared to engineering correlations available in the literature. Good comparisons between heat transfer rates are obtained at hypervelocity entry conditions provided the CFD code is modified to include equilibrium air properties.

Balakrishnan, A.↗

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.↗

Unsteady stagnation-point heat transfer during passage of a concentrated vortex

The unsteady boundary layer due to a single rectilinear vortex filament approaching a 2-D stagnation point is investigated. Assuming the vortex remains far from the surface, incompressible potential flow theory is used to determine the time dependent inviscid flow field. The unsteady boundary layer equations are solved by an alternating-direction-implicit finite-difference method. Two mechanisms which cause fluctuations in heat transfer are the unsteady velocity field in the boundary layer and secondly, the unsteady total temperature at the edge of the boundary layer. The relative importance of these mechanisms is dependent upon the total temperature fluctuations relative to the imposed temperature difference. As a vortex approaches a stagnation point it may be forced to one side of the stagnation line or the other, depending on its initial position. Results are presented for both of these cases.

Rigby, David L.↗

Solar coronal loop heating by cross-field wave transport

Solar coronal arches heated by turbulent ion-cyclotron waves may suffer significant cross-field transport by these waves. Nonlinear processes fix the wave-propagation speed at about a tenth of the ion thermal velocity, which seems sufficient to spread heat from a central core into a large cool surrounding cocoon. Waves heat cocoon ions both through classical ion-electron collisions and by turbulent stochastic ion motions. Plausible cocoon sizes set by wave damping are in roughly kilometers, although the wave-emitting core may be only 100 m wide. Detailed study of nonlinear stabilization and energy-deposition rates predicts that nearby regions can heat to values intermediate between the roughly electron volt foot-point temperatures and the about 100 eV core, which is heated by anomalous Ohmic losses. A volume of 100 times the core volume may be affected. This qualitative result may solve a persistent problem with current-driven coronal heating; that it affects only small volumes and provides no way to produce the extended warm structures perceptible to existing instruments.

Amendt, Peter↗

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.↗