Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

The propagation and stability of time-dependent galactodetonation waves

A new analytic approach is applied to the problem of star formation fronts in differentially rotating disks recently studied by Cowie and Rybicki. A quasi-linear equation for the inclination angle of the front is derived, and its characteristic equations are solved. The stability of these fronts found numerically in the previous study is revealed to be a direct consequence of the structure of the front characteristics. The time-dependent solutions found here are potentially useful in describing the evolution of local wave fronts shearing in a galactic velocity field, as observed in spirals lacking 'grand-design' structure.

Balbus, S. A.↗
NASA NTRS1985

Convenient stability criteria for difference approximations of hyperbolic initial-boundary value problems

New convenient stability criteria are provided in this paper for a large class of finite difference approximations to initial-boundary value problems associated with the hyperbolic system u sub t = au sub x + Bu + f in the quarter plane x or = 0, t or = 0. Using the new criteria, stability is easily established for numerous combinations of well known basic schemes and boundary conditins, thus generalizing many special cases studied in recent literature.

Goldberg, M.↗
NASA NTRS1985

Magnetohydrodynamic thermal instabilities in cool inhomogeneous atmospheres

The stability of magnetic loops to current-driven filamentation instabilities is investigated. The unperturbed atmosphere is assumed to be composed of an (upper) isothermal optically thin low-density portion and a (lower) higher-density portion which is in radiative equilibrium; in both cases, the atmosphere is in hydrostatic equilibrium, so that gravitational stratification is taken into account. In order to provide specific equilibrium conditions for evaluation of the dispersion relation, conditions appropriate for the surface of a solar-type star are adopted; i.e., a fairly low temperature (T = 5000 K) appropriate for a 'precoronal' state associated, for example, with magnetic flux emerging from photospheric levels under the action of magnetic buoyancy. A linear stability analysis is performed, and numerical results show that physically plausible current densities, which would be generated by typical loop-footpoint motions, are effective in driving MHD instabilities in such a plasma. The instability growth rates are strongly dependent on the assumed current density distribution and on the density scale height.

Bodo, G.↗
NASA NTRS1985

Dynamics of liquid sloshing in upright and inverted bladdered tanks

A system identification methodology was used to examine the dynamics of liquid sloshing in the upright and inverted bladdered hydrazine tanks of the Tracking and Data Relay Satellite, (TDRS) and to evaluate the effects of bladder stiffness on the sloshing parameters. Mechanical models of the two systems were developed using the numerical values derived from static stability tests and from slosh frequency response tests of a full-size model tank fitted with a prototype bladder. For the upright tank (liquid below the bladder) a modified conventional pendulum was used. In the inverted tank (liquid above the bladder) where sloshing is unconventional due to the highly nonsymmetrical orientation of the liquid held by the bladder, a mechanical model using an inverted pendulum which is able to undergo small oscillations as well as large reorientations was necessary. Both thrusting and low-gravity conditions are considered.

Dodge, F. T.↗
NASA NTRS1987

Aeroelastic analysis of launch vehicles in transonic flight

A numerical study of the aeroelastic stability of launch vehicles in transonic flight is performed. The unsteady aerodynanic flowfield about the vehicle is simulated using the three dimensional Reynolds-Averaged Navier-Stokes equations with the account for turbulent mixing done by a two-layer algebraic eddy viscosity model. The structural-dynamic formulation considers free-free flexural vibration of an elongated beam with variable properties and the equations are cast in modal form. The two sets of equations are integrated simultaneously in time, after a initial perturbation, to ascertain the aeroelastic stability of a given configuration. The method propossed is described in detail, and results that illustrate its application are presented.

Azevedo, Joao Luiz F.↗
NASA NTRS1988

A study of numerical methods for hyperbolic conservation laws with stiff source terms

The proper modeling of nonequilibrium gas dynamics is required in certain regimes of hypersonic flow. For inviscid flow this gives a system of conservation laws coupled with source terms representing the chemistry. Often a wide range of time scales is present in the problem, leading to numerical difficulties as in stiff systems of ordinary differential equations. Stability can be achieved by using implicit methods, but other numerical difficulties are observed. The behavior of typical numerical methods on a simple advection equation with a parameter-dependent source term was studied. Two approaches to incorporate the source term were utilized: MacCormack type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. Various comparisons over a wide range of parameter values were made. In the stiff case where the solution contains discontinuities, incorrect numerical propagation speeds are observed with all of the methods considered. This phenomenon is studied and explained.

Leveque, R. J.↗
NASA NTRS1988

Improved 3-D turbomachinery CFD algorithm

The building blocks of a computer algorithm developed for the time-accurate flow analysis of rotating machines are described. The flow model is a finite volume method utilizing a high resolution approximate Riemann solver for interface flux definitions. This block LU implicit numerical scheme possesses apparent unconditional stability. Multi-block composite gridding is used to orderly partition the field into a specified arrangement. Block interfaces, including dynamic interfaces, are treated such as to mimic interior block communication. Special attention is given to the reduction of in-core memory requirements by placing the burden on secondary storage media. Broad applicability is implied, although the results presented are restricted to that of an even blade count configuration. Several other configurations are presently under investigation, the results of which will appear in subsequent publications.

Janus, J. Mark↗
NASA NTRS1988

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗
NASA NTRS1989

Application of stochastic robustness to aircraft control systems

Stochastic robustness, a simple numerical procedure for estimating the stability robustness of linear, time-invariant systems, is applied to a forward-swept-wing aircraft control system. Based on Monte Carlo evaluation of the system's closed-loop eignevalues, this analysis approach introduces the probability of instability as a scalar stability robustness measure. The related stochastic root locus provides insight into robustness characteristics of the closed-loop system. Three Linear Quadratic controllers of decreasing robustness are chosen to demonstrate the use of stochastic robustness to analyze and compare control designs. Examples are presented illustrating the use of stochastic robustness analysis to address the effects of actuator dynamics and unmodeled dynamics on the stability robustness of the forward-swept-wing aircraft.

Stengel, Robert F.↗
NASA NTRS1990

Numerical investigation of laminar-turbulent transition in a flat plate wake

Lamina-turbulent transition of high-deficit flat plate wakes is investigated by direct numerical simulations using the complete Navier Stokes equations. The simulations are based on a spatial model so that both the base flow and the disturbance flow can develop in the downstream direction. The Navier Stokes equations are used in a vorticity-velocity form and are solved using a combination of finite difference and spectral approximations. Fourier series are used in the spanwise direction. Second-order finite-differences are used to approximate the spatial derivatives in the streamwise and transverse directions. For the temporal discretion, a combination of ADI, Crank-Nicolson, and Adams-Bashforth methods is employed. The discretized velocity equations are solved using fast Helmholtz solvers. Code validation is accomplished by comparison of the numerical results to both linear stability and to experiments. Calculations of two- and/or three-dimensional sinuous and mode disturbances in the wake of flat plate are undertaken. For calculations of two-dimensional disturbances, the wake is forced at an amplitude level so that nonlinear disturbance development may be observed. In addition, the forcing amplitude is varied in order to determine its effect on the disturbance behavior. To investigate the onset of three-dimensionality, the wake is forced with a small-amplitude three-dimensional disturbance and a larger amplitude two-dimensional disturbance. The two-dimensional forcing amplitude is varied in order to determine its influence on the three-dimensional flow field.

Fasel, Hermann F.↗
NASA NTRS1990

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions, are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1990

Instability during directional solidification - Gravitational effects

After an elementary introduction to solidification, recent research that attempts to understand the interaction of fluid flow with a crystal-melt interface under gravitational influence using linear stability analysis is reviewed. Recent numerical calculations of fluid flow and interface morphology during alloy solidification are briefly discussed.

Coriell, S. R.↗
NASA NTRS1992

The inviscid stability of supersonic flow past heated or cooled axisymmetric bodies

The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.

Shaw, Stephen J.↗
NASA NTRS1995

Thermally Stable, Piezoelectric and Pyroelectric Polymeric Substrates and Method Relating Thereto

Production of an electric voltage in response to mechanical excitation (piezoelectricity) or thermal excitation (pyroelectricity) requires a material to have a preferred dipole orientation in its structure. This preferred orientation or polarization occurs naturally in some crystals such as quartz and can be induced into some ceramic and polymeric materials by application of strong electric or mechanical fields. For some materials, a combination of mechanical and electrical orientation is necessary to completely polarize the material. The only commercially available piezoelectric polymer is poly(vinylidene fluoride) (PVF2). However, this polymer has material and process limitations which prohibit its use in numerous device applications where thermal stability is a requirement. By the present invention, thermally stable, piezoelectric and pyroelectric polymeric substrates were prepared from polymers having a softening temperature greater than 1000C. A metal electrode material is deposited onto the polymer substrate and several electrical leads are attached to it. The polymer substrate is heated in a low dielectric medium to enhance molecular mobility of the polymer chains. A voltage is then applied to the polymer substrate inducing polarization. The voltage is then maintained while the polymer substrate is cooled 'freezing in' the molecular orientation. The novelty of the invention resides in the process of preparing the piezoelectric and pyroelectric polymeric substrate. The nonobviousness of the invention is found in heating the polymeric substrate in a low dielectric medium while applying a voltage.

Simpson, Joycelyn O.↗
NASA NTRS2001

Panel Flutter and Sonic Fatigue Analysis for RLV

A methodology is presented for the flutter analysis of the seal of thermal protection system (TPS) panel of X-33 Advanced Technology Demonstrator test vehicle. The seal is simulated as a two-dimensional cantilevered panel with an elastic stopper, which is modeled as an equivalent spring. This cantilever beam-spring model under the aerodynamic pressure at supersonic speeds turns out to be an impact nonlinear dynamic system. The flutter analysis of the seal is thus carried out using, time domain numerical simulation with a displacement stability criterion. The flutter boundary of the seal is further verified with a family of three traditional and one nontraditional panel flutter models. The frequency domain method that applies eigenanalysis on the traditional panel flutter problem was used. The results showed that the critical dynamic pressure could be more than doubled with properly chosen material for the base stopper. The proposed methodology can be easily extended to three-dimensional panel seals with flow angularity.

Mei, Chuh↗
NASA NTRS2004

On the Total Variation of High-Order Semi-Discrete Central Schemes for Conservation Laws

We discuss a new fifth-order, semi-discrete, central-upwind scheme for solving one-dimensional systems of conservation laws. This scheme combines a fifth-order WENO reconstruction, a semi-discrete central-upwind numerical flux, and a strong stability preserving Runge-Kutta method. We test our method with various examples, and give particular attention to the evolution of the total variation of the approximations.

Bryson, Steve↗
NASA NTRS2010

Detailed Comparison of DNS to PSE for Oblique Breakdown at Mach 3

A pair of oblique waves at low amplitudes is introduced in a supersonic flat-plate boundary layer. Their downstream development and the concomitant process of laminar to turbulent transition is then investigated numerically using Direct Numerical Simulations (DNS) and Parabolized Stability Equations (PSE). This abstract is the last part of an extensive study of the complete transition process initiated by oblique breakdown at Mach 3. In contrast to the previous simulations, the symmetry condition in the spanwise direction is removed for the simulation presented in this abstract. By removing the symmetry condition, we are able to confirm that the flow is indeed symmetric over the entire computational domain. Asymmetric modes grow in the streamwise direction but reach only small amplitude values at the outflow. Furthermore, this abstract discusses new time-averaged data from our previous simulation CASE 3 and compares PSE data obtained from NASA's LASTRAC code to DNS results.

Mayer, Christian S. J.↗