Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stability”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

The Apollo Capsule Optimization for Improved Stability and Computational/Experimental Data Comparisons

Numerical optimization was employed on the Apollo Command Module to modify its external shape. The Apollo Command Module (CM) that was used on all NASA human space flights during the Apollo Space Program is stable and trimmed in an apex forward (alpha of approximately 40 to 80 degrees) position. This poses a safety risk if the CM separates from the launch tower during abort. Optimization was employed on the Apollo CM to remedy the undesirable stability characteristics of the configuration. Geometric shape changes were limited to axisymmetric modifications that altered the radius of the apex (R(sub A)), base radius (R(sub O)), corner radius (R(sub C)), and the cone half angle (theta), while the maximum diameter of the CM was held constant. The results of multipoint optimization on the CM indicated that the cross-range performance can be improved while maintaining robust apex-aft stability with a single trim point. Navier-Stokes computations were performed on the baseline and optimized configurations and confirmed the Euler-based optimization results. Euler Analysis of ten alternative CM vehicles with different values of the above four parameters are compared with the published experimental results of numerous wind tunnel tests during the late 1960's. These comparisons cover a wide Mach number range and a full 180-degree pitch range and show that the Euler methods are capable of fairly accurate force and moment computations and can separate the vehicle characteristics of these ten alternative configurations.

Cliff, Susan E.↗
NASA NTRS1993

Effects of gas on the global stability of galactic disks - Radial flows

We study numerically the effect of gas on the global stability of a two-component self-gravitating galactic disk embedded in a live halo. The stars are evolved by using a 3D collisionless N-body code, and the gas is represented by an ensemble of finite size inelastic particles. The gravitational interaction of stars and gas is calculated using a TREE method. We find that the evolution of the gaseous distribution in the globally unstable disks can be described by two different regimes. When the gas mass fraction is less than about 10 percent, the gas is channeled toward the galactic center by a growing stellar bar. For higher gas fractions, the gas becomes highly inhomogeneous, and the bar instability in the disk is heavily damped. The gas falls toward the inner kpc due to dynamical friction. Domains of both regimes depend on the efficiency of dissipation in the gas. We also discuss the relevance of the Jeans instability and give an empirical criterion for the global bar instability in a two-component self-gravitating disk.

Shlosman, Isaac↗
NASA NTRS1984

Stability of a laminar premixed supersonic free shear layer with chemical reactions

The stability of a two-dimensional compressible supersonic flow in the wake of a flat plate is discussed. The fluid is a multi-species mixture which is undergoing finite rate chemical reactions. The spatial stability of an infinitesimal disturbance in the fluid is considered. Numerical solutions of the eigenvalue stability equations for both reactive and nonreactive supersonic flows are presented and discussed. The chemical reactions have significant influence on the stability behavior. For instance, a neutral eigenvalue is observed near the freestream Mach number of 2.375 for the nonreactive case, but disappears when the reaction is turned on. For reactive flows, the eigenvalues are not very dependent on the free stream Mach number.

Menon, S.↗
NASA NTRS1988

Cage stability analysis for SSME HPOTP bearings

A numerical model of cage motion (CAGEDYN) was used to analyze the stability of bearing cages in the Space Shuttle main engine (SSME) high pressure oxygen turbopump (HPOTP). The stability of existing bearing geometries, as well as perturbations of these geometries, was analyzed for various operating conditions. Results of the analyses show that some combinations of operating parameters, exacerbated by the sparse lubrication that exist in the HPOTP bearings, can cause unstable cage oscillations. Frequencies of cage oscillations were predicted by the CAGEDYN numerical model by Fourier analysis of predicted cage motions. Under conditions that cause unstable cage motion, high frequency oscillations were predicted that could cause premature cage failures.

Merriman, T. L.↗
NASA NTRS1983

Finite element formulation for transient heat treat problems

The macrothermomechanical behavior of materials subjected to rapid thermal or mechanical loading such as occurs in most heat treatments is described. The equations are developed for Lagrangian, Eulerian, and intermediary kinematic descriptions and are independent of the constitutive laws and the equation of state; they can be solved numerically for a specified material and boundary conditions. The coupled transport effects between dissipation and energy are included. The conventional linearized stability approach indicates the numerical procedure to be stable, with certain restriction on the time step size.

Mullen, R. L.↗
NASA NTRS1983

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗
NASA NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗
NASA NTRS1993

Planetary circulations in the presence of transient and self-induced heating

The research program focuses on large-scale circulations and their interaction with the global convective pattern. An 11-year record of global cloud imagery and contemporaneous fields of motion and temperature have been used to investigate organized convection and coherent variability of the tropical circulation operating on intraseasonal time scales. This study provides a detailed portrait of tropical variability associated with the so-called Madden-Julian Oscillation (MJO). It reveals the nature, geographical distribution, and seasonality of discrete convective signal, which is a measure of feedback between the circulation and the convective pattern. That discrete spectral behavior has been evaluated in light of natural variability of the ITCZ associated with climatological convection. A composite signature of the MJO, based on cross-covariance statistics of cloud cover, motion, and temperature, has been constructed to characterize the lifecycle of the disturbance in terms of these properties. The composite behavior has also been used to investigate the influence the MJO exerts on the zonal-mean circulation and the involvement of the MJO in transfers of momentum between the atmosphere and the solid Earth. The aforementioned observational studies have led to the production of two animations. One reveals the convective signal in band-pass filtered OLR and compares it to climatological convection. The other is a 3-dimensional visualization of the composite lifecycle of the MJO. With a clear picture of the MJO in hand, feedback between the circulation and the convective pattern can be diagnosed meaningfully in numerical simulations. This process is being explored in calculations with the linearized primitive equations on the sphere in the presence of realistic stability and shear. The numerical framework represents climatological convection as a space-time stochastic process and wave-induced convection in terms of the vertically-integrated moisture flux convergence. In these calculations, frictional convergence near the equator emerges as a key to feedback between the circulation and the convective pattern. At low latitudes, nearly geostrophic balance in the boundary layer gives way to frictional balance. This shifts the wave-induced convection into phase with the temperature anomaly and allows the attending heating to feed back positively onto the circulation. The calculations successfully reproduce the salient features of the MJO. They are being used to understand the growth and decay phases of the composite lifecycle and the conditions that favor amplification of the MJO.

Salby, Murry L.↗
NASA NTRS1975

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method of stabilization of the differential equations of the Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as an independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Liapunov-stable set of differential equations. In contrast to other publications, here the equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗
NASA NTRS1976

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method for stabilizing the differential equations of Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. The equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗
NASA NTRS2005

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible. (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred. (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used. (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for non-linear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows.

Yee, H. C.↗
NASA NTRS1993

Nonlinear evolution of the first mode supersonic oblique waves in compressible boundary layers. Part 1: Heated/cooled walls

The nonlinear stability of an oblique mode propagating in a two-dimensional compressible boundary layer is considered under the long wave-length approximation. The growth rate of the wave is assumed to be small so that the concept of unsteady nonlinear critical layers can be used. It is shown that the spatial/temporal evolution of the mode is governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity and density. Expressions for the linear growth rate show clearly the effects of wall heating and cooling and in particular how heating destabilizes the boundary layer for these long wavelength inviscid modes at O(1) Mach numbers. A generalized expression for the linear growth rate is obtained and is shown to compare very well for a range of frequencies and wave-angles at moderate Mach numbers with full numerical solutions of the linear stability problem. The numerical solution of the nonlinear unsteady critical layer problem using a novel method based on Fourier decomposition and Chebychev collocation is discussed and some results are presented.

Gajjar, J. S. B.↗
NASA NTRS1976

Elliptic systems and numerical transformations

Properties of a transformation method, which was developed for solving fluid dynamic problems on general two dimensional regions, are discussed. These include construction error of the transformation and applications to mesh generation. An error and stability analysis for the numerical solution of a model parabolic problem is also presented.

Mastin, C. W.↗
NASA NTRS1990

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗
NASA NTRS1990

The inviscid stability of supersonic flow past axisymmetric bodies

The supersonic flow past a sharp cone is studied. The associated boundary layer flow (i.e., the velocity and temperature field) is computed. The inviscid linear temporal stability of axisymmetric boundary layers in general is considered, and in particular, a so-called 'triply generalized' inflection condition for 'subsonic' nonaxisymmetric neutral modes is presented. Preliminary numerical results for the stability of the cone boundary layer are presented for a freestream Mach number of 3.8. In particular, a new inviscid mode of instability is seen to occur in certain regimes, and this is shown to be related to a viscous mode found by Duck and Hall (1988).

Duck, Peter W.↗
NASA NTRS1952

Calculation of the lateral-dynamic stability of aircraft

Graphs and formulas are given with the aid of which all the aerodynamic coefficients required for computing the lateral dynamic stability can be determined. A number of numerical examples are given for obtaining the stability derivatives and solving the characteristic-stability equation. Approximate formulas are derived with the aid of which rapid preliminary computations may be made and the stability coefficients corrected for certain modifications of the airplane. A derivation of the lateral-dynamic-stability equations is included.

STABILITY, LATERAL - DYNAMIC↗
NASA NTRS2016

Calculated Low-Speed Steady and Time-Dependent Aerodynamic Derivatives for Several Different Wings Using a Discrete Vortex Method

Calculated numerical values for some aerodynamic terms and stability Derivatives for several different wings in unseparated inviscid incompressible flow were made using a discrete vortex method involving a limited number of horseshoe vortices. Both longitudinal and lateral-directional derivatives were calculated for steady conditions as well as for sinusoidal oscillatory motions. Variables included the number of vortices used and the rotation axis/moment center chordwise location. Frequencies considered were limited to the range of interest to vehicle dynamic stability (kb <.24 ). Comparisons of some calculated numerical results with experimental wind-tunnel measurements were in reasonable agreement in the low angle-of-attack range considering the differences existing between the mathematical representation and experimental wind-tunnel models tested. Of particular interest was the presence of induced drag for the oscillatory condition.

Riley, Donald R.↗