SEARCH · Engineering Papers

Results for “Conservation equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A hybrid multigrid technique for computing steady-state solutions to supersonic flows

Recently, Li and Sanders have introduced a class of finite difference schemes to approximate generally discontinuous solutions to hyperbolic systems of conservation laws. These equations have the form together with relevant boundary conditions. When modelling hypersonic spacecraft reentry, the differential equations above are frequently given by the compressible Euler equations coupled with a nonequilibrium chemistry model. For these applications, steady state solutions are often sought. Many tens (to hundreds) of super computer hours can be devoted to a single three space dimensional simulation. The primary difficulty is the inability to rapidly and reliably capture the steady state. In these notes, we demonstrate that a particular variant from the schemes presented can be combined with a particular multigrid approach to capture steady state solutions to the compressible Euler equations in one space dimension. We show that the rate of convergence to steady state coming from this multigrid implementation is vastly superior to the traditional approach of artificial time relaxation. Moreover, we demonstrate virtual grid independence. That is, the rate of convergence does not depend on the degree of spatial grid refinement.

Sanders, Richard↗

NASA NTRS1993

The reionization of the universe: The feedback of galaxy formation on the intergalactic medium

The thermal and ionization evolution of a uniform intergalactic medium (IGM) composed of H and He, undergoing reionization, including the mean effect of gas clumps embedded in a smoothly distributed ambient gas were calculated. The rate equations for ionization and recombination were solved together with the equations of energy conservation, including the effects of cosmological expansion, radiative and Compton cooling, and the diffuse flux emitted by the gas, and radiative transfer. The contribution to the continuum opacity of the universe due to the observed quasar absorption line clouds (QALC'S) were included. A variety of sources of photoionization, including quasars and primeval galaxies, as well as the possibility that hydrodynamical processes deposit thermal energy in the IGM were considered. Applications of these calculations including the evolution of the Ly-alpha forest clouds are described. A self-consistent treatment of the thermal and ionization history of the intergalactic medium (IGM) must take account of the growth of structure in the universe, since the mean density of the IGM corresponds primarily to the time-varying uncollapsed fraction of the baryon-electron component of the matter, and the collapsed fraction, in turn, can have a feedback effect on this uncollapsed fraction by releasing ionizing radiation and thermal energy and by contributing to the opacity of the universe. The coupled evolution of the IGM and the emerging structure with a special focus on the reionization of the IGM, which is believed to have been completed by some redshift z is approximately greater than 4, as inferred from the absence of the Gunn-Peterson effect in the spectra of high z quasars, are studied. The results and implications of detailed, numerical calculations of the thermal and ionization balance and radiative transfer in a uniform IGM of H and He, including the mean effect of an evolving distribution of gas clumps embedded in a smoothly distributed ambient gas is described.

Shapiro, Paul R.↗

NASA NTRS1993

Global fields of soil moisture and land surface evapotranspiration derived from observed precipitation and surface air temperature

The global fields of normal monthly soil moisture and land surface evapotranspiration are derived with a simple water budget model that has precipitation and potential evapotranspiration as inputs. The precipitation is observed and the potential evapotranspiration is derived from the observed surface air temperature with the empirical regression equation of Thornthwaite (1954). It is shown that at locations where the net surface radiation flux has been measured, the potential evapotranspiration given by the Thornthwaite equation is in good agreement with those obtained with the radiation-based formulations of Priestley and Taylor (1972), Penman (1948), and Budyko (1956-1974), and this provides the justification for the use of the Thornthwaite equation. After deriving the global fields of soil moisture and evapotranspiration, the assumption is made that the potential evapotranspiration given by the Thornthwaite equation and by the Priestley-Taylor equation will everywhere be about the same; the inverse of the Priestley-Taylor equation is used to obtain the normal monthly global fields of net surface radiation flux minus ground heat storage. This and the derived evapotranspiration are then used in the equation for energy conservation at the surface of the earth to obtain the global fields of normal monthly sensible heat flux from the land surface to the atmosphere.

Mintz, Y.↗

NASA NTRS1994

Theoretical studies of a molecular beam generator

We consider the problem of modeling a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem is modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann population distribution which means that each degree of freedom has its own temperature and consequently each system state can be characterized by these temperatures. This suggests the formulation of a second model with a vibrational degree of freedom as having its own temperature along with a rotational-translational degree of freedom with its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas expands through the nozzle throat there is a sudden drop in the rotational-translational temperature along with a finite relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translational degree of freedom. That is, the temperature change that occurs when the N2 gas passes through the nozzle throat is so great that the changes in the vibrational degree of freedom lags behind the rotational-translational energy changes. The resulting relaxation time is finite. It is in this context that the term nonequilibrium is used. That is, the term nonequilibrium denotes the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions resulting from the high temperatures which could also add nonequilibrium effects. We develop the basic equations for the two models in various forms in order to check the derivations with other sources. The final form which is solved numerically consists of the scaled equations in a conservative dimensionless form.

Heinbockel, John H.↗

NASA NTRS1995

Computation of three-dimensional mixed convective boundary layer flow

The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.

Gadepalli, Prashandt↗

NASA NTRS1996

Computation Of Flow About A Helicopter Rotor

Report describes computations of flow about hovering helicopter rotor. Thin-layer Navier-Stokes equations in conservation-law form solved numerically on embedded grids.

Duque, Earl P.↗

NASA NTRS1996

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗

NASA NTRS1997

Shapes of Buoyant and Nonbuoyant Methane Laminar Jet Diffusion Flames

Laminar gas jet diffusion flames represent a fundamental combustion configuration. Their study has contributed to numerous advances in combustion, including the development of analytical and computational combustion tools. Laminar jet flames are pertinent also to turbulent flames by use of the laminar flamelet concept. Investigations into the shapes of noncoflowing microgravity laminar jet diffusion flames have primarily been pursued in the NASA Lewis 2.2-second drop tower, by Cochran and coworkers and by Bahadori and coworkers. These studies were generally conducted at atmospheric pressure; they involved soot-containing flames and reported luminosity lengths and widths instead of the flame-sheet dimensions which are of Greater value to theory evaluation and development. The seminal model of laminar diffusion flames is that of Burke and Schumann, who solved the conservation of momentum equation for a jet flame in a coflowing ambient by assuming the velocity of fuel, oxidizer and products to be constant throughout. Roper and coworkers improved upon this model by allowing for axial variations of velocity and found flame shape to be independent of coflow velocity. Roper's suggestion that flame height should be independent of gravity level is not supported by past or present observations. Other models have been presented by Klajn and Oppenheim, Markstein and De Ris, Villermaux and Durox, and Li et al. The common result of all these models (except in the buoyant regime) is that flame height is proportional to fuel mass flowrate, with flame width proving much more difficult to predict. Most existing flame models have been compared with shapes of flames containing soot, which is known to obscure the weak blue emission of flame sheets. The present work involves measurements of laminar gas jet diffusion flame shapes. Flame images have been obtained for buoyant and nonbuoyant methane flames burning in quiescent air at various fuel flow-rates, burner diameters and ambient pressures. Soot concentrations were minimized by selecting conditions at low flowrates and low ambient pressures; this allows identification of actual flame sheets associated with blue emissions of CH and CO2. The present modeling effort follows that of Roper and is useful in explaining many of the trends observed.

Sunderland, Peter B.↗

NASA NTRS1998

A 2.5D Single Passage CFD Model for Centrifugal Pumps

This paper describes the single passage model based on CFD to analyze the flow in blade passages of a centrifugal pump. The model consists of the flow passage between two impeller blades and the spaces in the inlet eye as well as in the volute. The incompressible Navier-Stokes equations in the conservation form are solved by a finite difference method. The code is designed to investigate the velocity and pressure distributions and intended to investigate how the pump design affects fluid flow through the rotor as well as the pump performance. An early part of the paper investigates the behavior of the model as well as validity of the assumptions made in the model. Then, applications to a rotodynamic heart pump are presented.

Nakamura S.↗

NASA NTRS1999

A One Dimensional, Time Dependent Inlet/Engine Numerical Simulation for Aircraft Propulsion Systems

The NASA Lewis Research Center (LeRC) and the Arnold Engineering Development Center (AEDC) have developed a closely coupled computer simulation system that provides a one dimensional, high frequency inlet/engine numerical simulation for aircraft propulsion systems. The simulation system, operating under the LeRC-developed Application Portable Parallel Library (APPL), closely coupled a supersonic inlet with a gas turbine engine. The supersonic inlet was modeled using the Large Perturbation Inlet (LAPIN) computer code, and the gas turbine engine was modeled using the Aerodynamic Turbine Engine Code (ATEC). Both LAPIN and ATEC provide a one dimensional, compressible, time dependent flow solution by solving the one dimensional Euler equations for the conservation of mass, momentum, and energy. Source terms are used to model features such as bleed flows, turbomachinery component characteristics, and inlet subsonic spillage while unstarted. High frequency events, such as compressor surge and inlet unstart, can be simulated with a high degree of fidelity. The simulation system was exercised using a supersonic inlet with sixty percent of the supersonic area contraction occurring internally, and a GE J85-13 turbojet engine.

Doug Garrard↗

NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗

NASA NTRS2002

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one- dimensional convection-diffusion-reaction (CDR) equations. Accuracy, stability, conservation, and dense-output are first considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, (N = 2), additive Runge-Kutta (ARK(sub 2)) methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms of the partitioned method are of equal order to those of the elemental methods. Derived ARK(sub 2) methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z['] yields -infinity, and retain high stability efficiency in the absence of stiffness, z['] yield 0. Extrapolation-type stage- value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK(sub 2) error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK(sub 2) methods.

Kennedy, Christopher A.↗

NASA NTRS2001

An Implicit LU/AF FDTD Method

There has been some recent work to develop two and three-dimensional alternating direction implicit (ADI) FDTD schemes. These ADI schemes are based upon the original ADI concept developed by Peaceman and Rachford and Douglas and Gunn, which is a popular solution method in Computational Fluid Dynamics (CFD). These ADI schemes work well and they require solution of a tridiagonal system of equations. A new approach proposed in this paper applies a LU/AF approximate factorization technique from CFD to Maxwell s equations in flux conservative form for one space dimension. The result is a scheme that will retain its unconditional stability in three space dimensions, but does not require the solution of tridiagonal systems. The theory for this new algorithm is outlined in a one-dimensional context for clarity. An extension to two and threedimensional cases is discussed. Results of Fourier analysis are discussed for both stability and dispersion/damping properties of the algorithm. Results are presented for a one-dimensional model problem, and the explicit FDTD algorithm is chosen as a convenient reference for comparison.

Beggs, John H.↗

NASA NTRS2009

Adiabatic Shock Capturing in Perfect Gas Hypersonic Flows

This paper considers the streamline-upwind Petrov/Galerkin (SUPG) method applied to the compressible Euler and Navier-Stokes equations in conservation-variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, is briefly reviewed. Of particular interest is the behavior of the shock capturing operator, which is required to regularize the scheme in the presence of strong, shock-induced gradients. A standard shock capturing operator which has been widely used in previous studies by several authors is presented and discussed. Specific modifications are then made to this standard operator which are designed to produce a more physically consistent discretization in the presence of strong shock waves. The actual implementation of the term in a finite dimensional approximation is also discussed. The behavior of the standard and modified scheme is then compared for several supersonic/hypersonic flows. The modified shock capturing operator is found to preserve enthalpy in the inviscid portion of the flowfield substantially better than the standard operator.

Kirk, Benjamin S.↗

NASA NTRS2014

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

NASA NTRS1984

Entropy condition satisfying approximations for the full potential equations of transonic flow

A class of conservative difference approximations for the steady full potential equation was presented. They are, in general, easier to program than the usual density biasing algorithms, and in fact, differ only slightly from them. Rigorous proof indicated that these new schemes satisfied a new discrete entropy inequality, which ruled out expansion shocks, and that they have sharp, steady, discrete shocks. A key tool in the analysis is the construction of a new entropy inequality for the full potential equation itself. Results of some numerical experiments using the new schemes are presented.

Osher, S.↗

NASA NTRS1985

Entropy condition satisfying approximations for the full potential equation of transonic flow

Osher, S.↗

NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗