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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A numerical study of the nonlinear thermal stability of solar loops

A time-dependent numerical model is used to investigate the nonlinear thermal stability of static loops of various heights. Simulations show that the instability of a hot state with loop heights of less than about 1000 km is physically significant, with an initially hot atmosphere in low-lying compact loops evolving to an extended atmosphere with temperatures far below 100,000 K. Results also show that high-lying loops are stable to all reasonable perturbations, including those of large initial amplitude and long wavelength. The simulation results suggest that low-lying compact loops should not be common to the sun, and that cool loops with temperatures near 100,000 K must be formed in the cool state initially and cannot evolve from preexisiting loops.

Klimchuk, J. A.↗
NASA NTRS1992

Thermal stability during rotation in space - A scaling and numerical analysis

An analysis of thermal stability in a fluid during rotation is presented. The conditions resulting in thermal stability are examined and numerically verified. It is shown that at least one thermal field exists for any rotating fluid, with or without a gravitational background field, where convection does not occur. The numerical model used is fully nonlinear and 3D. Results are presented on the basis of both steady-state and fully transient models. The thermally stable condition is further investigated in the context of the average net acceleration and the inhomogeneous part of the acceleration field (gradient of acceleration). It is shown that in the thermally stable condition the flow that would result from the average net acceleration is equal and opposite to the flow that would result from the inhomogeneous part of the acceleration field. It is concluded that the thermally stable condition is actually a flow mode transition at Ad = 1.

Arnold, William A.↗
NASA NTRS1993

Numerical solution of wave equations for the stability of the inner cometo-sheath

Numerical solution of the MHD wave equations for stability of the cometary sheath determined by the balance between the inward Lorentz body force and the outward ion-neutral drag force is obtained by using a two-point boundary value method. The eigenvalues and the eigenfunctions are obtained numerically by treating the cometary inner sheath as a layer of finite thickness, bounded by the contact surface, i.e., the diamagnetic cavity boundary. The magnetic field structure discovered in the ionosphere of Comets Halley and Giacobini-Zinner is found to be unstable. The effects of finite plasma pressure, dissociative recombination, and mass loading due to photoionization are found to be stabilizing but are unable to quench the instability completely. It is also found that the higher the neutral production rate the lesser is the growth rate for the instability.

Srivastava, Krishna M.↗
NASA NTRS1986

Impact of turbulence modeling on numerical accuracy and efficiency of compressible flow simulations

Discussed is the numerical implementation of turbulence models used in viscous compressible flow simulations and their performance described with respect to numerical accuracy, efficiency and stability. The basic approach utilizes the Reynolds-averaged compressible Navier-Stokes equations in which the Reynolds stresses and heat fluxes are mathematically modeled by suitable turbulence models. The turbulence models investigated include zero-, one-, and two-equation eddy viscosity models. The flow fields investigated include theoretical and supersonic flows about two-dimensional and axisymmetric bodies. Discussions concerning the numerical implementation of models include differencing procedures and boundary conditions used to assume numerical stability and accuracy. Numerical performance is also evaluated by comparing computations with experimental results.

Coakley, Thomas J.↗
NASA NTRS1980

A numerical study of a marine subtropical stratus cloud layer and its stability

A model for numerical simulation of stratus cloud layers is constructed by combining a second-order closure, turbulent transfer model with a thermal radiative transfer model. The turbulent transfer model allows water vapor saturation. The combined turbulence-radiation model is applied to both a horizontally uniform one-dimensional case and a horizontally nonuniform two-dimensional case. In the latter, the dynamics of mesoscale circulations are also incorporated. Results of the two-dimensional simulation show that the layer cloud instability occurs where the sea surface temperature is high and the large-scale subsidence is weak. The simulated instability is analyzed in view of an instability criterion, the eddy kinetic energy budget, and evaporative cooling near the cloud top.

Moeng, C.-H.↗
NASA NTRS1992

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗
NASA NTRS1981

A Galerkin type finite element method for rotary-wing aeroelasticity in hover and forward flight

A Galerkin finite element method for the spatial discretization of the nonlinear, nonselfadjoint, partial differential equations governing rotary-wing aeroelasticity is presented. This method reduces algebraic manipulative labor significantly when compared to the global Galerkin method based on assumed modes. Furthermore, the Galerkin finite element method is ideally suited to treat rotor blades with discontinuous mass and stiffness distribution and structurally redundant configurations as they appear in bearingless rotors. Implementation of the method is illustrated for the coupled flap-lag aeroelastic problem of hingeless rotor blades in hover and forward flight. Numerical results for stability and response illustrate the numerical properties and convergence behavior of the method. It is concluded that the Galerkin finite element method is a practical tool for solving rotary-wing aeroelastic stability and response problems.

Straub, F. K.↗
NASA NTRS1993

Analytical and Numerical Modeling of Flame-Balls in Hydrogen-Air Mixtures

Flame-balls are stationary spherical premixed flames observed in certain near-limit mixtures. It is believed that radiative heat losses are an important stabilizing influence. Numerical solutions of flame balls are constructed for hydrogen-air mixtures using an accurate description of the chemical kinetics, diffusive transport, and radiation losses. A lean limit equivalence ratio of 0.0866 is predicted and a rich limit of 2,828. For any equivalence ratio between the two limits there are two solutions. One is characterized by a small flame, incomplete reactant consumption, and negligible radiation losses. The other by a large flame, complete consumption of one of the reactants, and significant radiation losses. The maximum temperature varies between 1200 and 900 K as the two solution branches are traversed. Much of our discussion is a reprise and modification of previously published analytical results, for these provide physical insight into the nature of the solutions, and suggest that a portion of the large flame branch near the lean limit is stable and so corresponds to observable flames.

Buckmaster, J.↗
NASA NTRS1993

Simulations of Multidimensional Burner-stabilized Flames

Detailed numerical simulations have been performed to study the structure and dynamics of downward propagating burner-stabilized flames in lean hydrogen-air mixtures. These simulations include the effects of fluid convection, detailed hydrogen-oxygen chemistry, multispecies diffusion, thermal conduction, viscosity, and heat losses to the burner. One-dimensional calculations have been carried out to investigate burner boundary conditions. Well known features of a burner-stabilized flame have been reproduced. Two-dimensional calculations show the presence of cellular structures at the burner surface. At low inlet velocities, these cellular structures are suppressed by the increased heat loss to the burner. These simulations have been examined in detail to gain understanding of the similarities and differences in structure and stability of freely propagating and burner-stabilized flames.

Patnaik, G.↗
NASA NTRS1986

A numerical study of the thermal stability of low-lying coronal loops

The nonlinear evolution of loops that are subjected to a variety of small but finite perturbations was studied. Only the low-lying loops are considered. The analysis was performed numerically using a one-dimensional hydrodynamical model developed at the Naval Research Laboratory. The computer codes solve the time-dependent equations for mass, momentum, and energy transport. The primary interest is the active region filaments, hence a geometry appropriate to those structures was considered. The static solutions were subjected to a moderate sized perturbation and allowed to evolve. The results suggest that both hot and cool loops of the geometry considered are thermally stable against amplitude perturbations of all kinds.

Klimchuk, J. A.↗
NASA NTRS1993

Linear stability of streamwise developing longitudinal vortices

The linear stability of numerical solutions to the quasicylindrical equations of motion for swirling flows is investigated. Initial conditions are derived from Batchelor's similarity solution for a trailing line vortex. The stability calculations are performed using a 2nd-order accurate finitedifference scheme on a staggered grid with the accuracy of the computed eigenvalues enhanced through Richardson extrapolation. The streamwise development of both viscous and inviscid instability modes is discussed for low and high swirl cases. The effect of adverse and favorable pressure gradients on the viscous stability of a longitudinal vortex is also investigated.

Spall, Robert E.↗
NASA NTRS1988

Helicopter stability and control modeling improvements and verification on two helicopters

A linearized model of helicopter flight dynamics is developed which includes the flapping, lead-lag, and dynamic inflow degrees of freedom (DOF). The model is a combination of analytical terms and numerically determined stability derivatives, and is used to investigate the importance of the rotor DOF to stability and control modeling. The results show that the rotor DOF can have a significant impact on some of the natural modes in a linear model. The flap and dynamic inflow DOF show the greatest influence. Flapping exhibits strong coupling to the body, dynamic inflow, and to lead-lag to a lesser extent. Dynamic inflow tends to damp the high-frequency flapping modes, and reduces the damping on coupled body-flap motion. Dynamic inflow also couples to the flapping motion to produce complex roots. With body-flap and lag regressing modes as exceptions, the results show essentially similar behavior for most modes of articulated and hingeless rotor helicopters.

Schrage, D. P.↗
NASA NTRS1982

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗
NASA NTRS1989

Hyperbolic heat conduction problems involving non-Fourier effects - Numerical simulations via explicit Lax-Wendroff/Taylor-Galerkin finite element formulations

Numerical simulations are presented for hyperbolic heat-conduction problems that involve non-Fourier effects, using explicit, Lax-Wendroff/Taylor-Galerkin FEM formulations as the principal computational tool. Also employed are smoothing techniques which stabilize the numerical noise and accurately predict the propagating thermal disturbances. The accurate capture of propagating thermal disturbances at characteristic time-step values is achieved; numerical test cases are presented which validate the proposed hyperbolic heat-conduction problem concepts.

Tamma, Kumar K.↗
NASA NTRS1989

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗
NASA NTRS1991

The stability to two-dimensional wakes and shear layers at high Mach numbers

This study is concerned with the stability properties of laminar free-shear-layer flows, and in particular symmetric two-dimensional wakes, for the supersonic through the hypersonic regimes. Emphasis is given to the use of proper wake profiles that satisfy the equations of motion at high Reynolds numbers. In particular the inviscid stability of a developing two-dimensional wake is studied as it accelerates at the trailing edge of a splitter plate. The nonparallelism of the flow is a leading-order effect in the calculation of the basic state, which is obtained numerically. Neutral stability characteristics are computed and the hypersonic stability is obtained by increasing the Mach number. It is found that the stability characteristics are altered significantly as the wake develops. Multiple modes (secondary modes) are found in the near wake that are closely related to the corresponding Blasius ones, but as the wake develops mode multiplicity is delayed to higher and higher Mach numbers. At a distance of about one plate length from the trailing edge, there is only one mode in a Mach number range of 0-20. The dominant mode emerging at all wake stations, and for high enough Mach numbers, is the so-called vorticity mode that is centered around the generalized inflection point layer. The structure of the dominant mode is also obtained analytically for all streamwise wake locations and it is shown how the far-wake limit is approached. Asymptotic results for the hypersonic mixing layer given by a tanh and a Lock distribution are also given.

Papageorgiou, Demetrios T.↗