Engineering Papers⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Energy-momentum-conserving stochastic differential equations and algorithms for the nonlinear Landau-Fokker-Planck equation

Coulomb collision is a fundamental diffusion process in plasmas that can be described by the Landau-Fokker-Planck (LFP) equation or the stochastic differential equation (SDE). While energy and momentum are conserved exactly in the LFP equation, they are conserved only on average by the conventional corresponding SDEs, suggesting that the underlying stochastic process may not be well defined by such SDEs. Here, in this study, we derive new SDEs with exact energy-momentum conservation for the Coulomb collision by factorizing the collective effect of field particles into individual particles and enforcing Newton's third law. These SDEs, when interpreted in the Stratonovich sense, have a particularly simple form that represents pure diffusion between particles without drag. To demonstrate that the new SDEs correspond to the LFP equation, we develop numerical algorithms that converge to the SDEs and preserve discrete conservation laws. Simulation results are presented in a benchmark of various relaxation processes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Multiple-compartment venting program

Computer program solves time-dependent energy and state equations for gas reservoirs, using the solutions of conductor conservation equations as mass and energy rate changes to reservoirs.

Leblanc, L. P.↗
NASA NTRS1995

Eight-moment approximation solar wind models

Heat conduction from the corona is important in the solar wind energy budget. Until now all hydrodynamic solar wind models have been using the collisionally dominated gas approximation for the heat conductive flux. Observations of the solar wind show particle distribution functions which deviate significantly from a Maxwellian, and it is clear that the solar wind plasma is far from collisionally dominated. We have developed a numerical model for the solar wind which solves the full equation for the heat conductive flux together with the conservation equations for mass, momentum, and energy. The equations are obtained by taking moments of the Boltzmann equation, using an 8-moment approximation for the distribution function. For low-density solar winds the 8-moment approximation models give results which differ significantly from the results obtained in models assuming the gas to be collisionally dominated. The two models give more or less the same results in high density solar winds.

Olsen, Espen Lyngdal↗
NASA NTRS1996

Classical-Quantum Correspondence by Means of Probability Densities

Within the frame of the recently introduced phase space representation of non relativistic quantum mechanics, we propose a Lagrangian from which the phase space Schrodinger equation can be derived. From that Lagrangian, the associated conservation equations, according to Noether's theorem, are obtained. This shows that one can analyze quantum systems completely in phase space as it is done in coordinate space, without additional complications.

Vegas, Gabino Torres↗
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗
NASA NTRS1974

Thermoacoustic convection of fluids in low gravity

The heat flow in a confined perfect gas in low gravity is investigated, including the effects of conduction and thermal convection. Buoyancy-driven flow is neglected, due to the low-gravity environment, but the effect of thermoacoustic motion due to fluid compressibility is included. One-dimensional mathematical models are constructed from the conservation equations for a compressible, viscous, heat-conducting fluid. A conservative, time-dependent finite-difference method is used to generate numerical solutions on a digital computer. Problems for flat plates and cylindrical segments are solved for specified thermal boundary conditions. Numerical results are given which indicate that thermoacoustic convection can significantly increase the transient heat flow over conduction model predictions for cases where a confined gas is rapidly heated.

Spradley, L. W.↗
DOE OSTI2024

Null Raychaudhuri: canonical structure and the dressing time

Abstract We initiate a study of gravity focusing on generic null hypersurfaces, non-perturbatively in the Newton coupling. We present an off-shell account of the extended phase space of the theory, which includes the expected spin-2 data as well as spin-0, spin-1 and arbitrary matter degrees of freedom. We construct the charges and the corresponding kinematic Poisson brackets, employing a Beltrami parameterization of the spin-2 modes. We explicitly show that the constraint algebra closes, the details of which depend on the non-perturbative mixing between spin-0 and spin-2 modes. Finally we show that the spin zero sector encodes a notion of a clock, called dressing time, which is dynamical and conjugate to the constraint. It is well-known that the null Raychaudhuri equation describes how the geometric data of a null hypersurface evolve in null time in response to gravitational radiation and external matter. Our analysis leads to three complementary viewpoints on this equation. First, it can be understood as a Carrollian stress tensor conservation equation. Second, we construct spin-0, spin-2 and matter stress tensors that act as generators of null time reparametrizations for each sector. This leads to the perspective that the null Raychaudhuri equation can be understood as imposing that the sum of CFT-like stress tensors vanishes. Third, we solve the Raychaudhuri constraint non-perturbatively. The solution relates the dressing time to the spin-2 and matter boost charge operators. Finally we establish that the corner charge corresponding to the boost operator in the dressing time frame is monotonic. These results show that the notion of an observer can be thought of as emerging from the gravitational degrees of freedom themselves. We briefly mention that the construction offers new insights into focusing conjectures.

Physics↗
NASA NTRS1989

Variational objective analysis for cyclone studies

Significant accomplishments during 1987 to 1988 are summarized with regard to each of the major project components. Model 1 requires satisfaction of two nonlinear horizontal momentum equations, the integrated continuity equation, and the hydrostatic equation. Model 2 requires satisfaction of model 1 plus the thermodynamic equation for a dry atmosphere. Model 3 requires satisfaction of model 2 plus the radiative transfer equation. Model 4 requires satisfaction of model 3 plus a moisture conservation equation and a parameterization for moist processes.

Achtemeier, Gary L.↗
NASA NTRS2002

A Novel Approach for Modeling Chemical Reaction in Generalized Fluid System Simulation Program

The Generalized Fluid System Simulation Program (GFSSP) is a computer code developed at NASA Marshall Space Flight Center for analyzing steady state and transient flow rates, pressures, temperatures, and concentrations in a complex flow network. The code, which performs system level simulation, can handle compressible and incompressible flows as well as phase change and mixture thermodynamics. Thermodynamic and thermophysical property programs, GASP, WASP and GASPAK provide the necessary data for fluids such as helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, water, a hydrogen, isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, several refrigerants, nitrogen trifluoride and ammonia. The program which was developed out of need for an easy to use system level simulation tool for complex flow networks, has been used for the following purposes to name a few: Space Shuttle Main Engine (SSME) High Pressure Oxidizer Turbopump Secondary Flow Circuits, Axial Thrust Balance of the Fastrac Engine Turbopump, Pressurized Propellant Feed System for the Propulsion Test Article at Stennis Space Center, X-34 Main Propulsion System, X-33 Reaction Control System and Thermal Protection System, and International Space Station Environmental Control and Life Support System design. There has been an increasing demand for implementing a combustion simulation capability into GFSSP in order to increase its system level simulation capability of a liquid rocket propulsion system starting from the propellant tanks up to the thruster nozzle for spacecraft as well as launch vehicles. The present work was undertaken for addressing this need. The chemical equilibrium equations derived from the second law of thermodynamics and the energy conservation equation derived from the first law of thermodynamics are solved simultaneously by a Newton-Raphson method. The numerical scheme was implemented as a User Subroutine in GFSSP.

Sozen, Mehmet↗
NASA NTRS2003

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗
NASA NTRS2004

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗
NASA NTRS2003

Effects of Mean Flow Profiles on Instability of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to investigate the effects of the mean flow profiles on the instability characteristics in the near-injector region of low-density gas jets injected into high-density ambient gas mediums. To achieve this, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round gas jet injected vertically upwards into a high-density ambient gas were performed by assuming three different sets of mean velocity and density profiles. The flow was assumed to be isothermal and locally parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The first set of mean velocity and density profiles assumed were those used by Monkewitz and Sohn for investigating absolute instability in hot jets. The second set of velocity and density profiles assumed for this study were the ones used by Lawson. And the third set of mean profiles included a parabolic velocity profile and a hyperbolic tangent density profile. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results for each set of mean profiles were delineated. Additional information is included in the original extended abstract.

Vedantam, Nanda Kishore↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

The efficiency of several algorithms used for numerical integration of stiff ordinary differential equations was compared. The methods examined included two general purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes were applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code available for the integration of combustion kinetic rate equations. It is shown that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient then evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗
NASA NTRS1993

A new approach to the study of the turbulence layer

Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.

Levenson, Ruth A.↗
NASA NTRS1984

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗
NASA NTRS1984

A comparison of the efficiency of numerical methods for integrating chemical kinetic rate equations

A comparison of the efficiency of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations is presented. The methods examined include two general-purpose codes EPISODE and LSODE and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an iterative solution of the algebraic energy conservation equation to compute the temperature can be more efficient than evaluating the temperature by integrating its time-derivative.

Radhakrishnan, K.↗