Search NASA⌕ Search

SEARCH · Search NASA

Results for “energy 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Navier-Stokes Entropy Controlled Combustion Instability Analysis for Liquid Propellants

Navier-Stokes solutions are used to calculate oscillatory components of pressure, velocity, and density, which in turn provide necessary data to compute energy growth factors to determine combustion instability. It is shown that wave instabilities are associated with changes in entropy and the space and time averages of oscillatory components of pressure, velocity and density, together with the mean flow field in the energy equation. Compressible laminar and turbulent flows and reacting flows with hydrogen/oxygen combustion are considered. The SSME combustion/thrust chamber is used for illustration of the theory. The analysis shows that the increase of mean pressure and disturbances consistently results in the increase of instability. It is shown that adequate combustion instability analysis requires at least third order nonlinearity in energy growth or decay.

Chung, T. J.↗
NASA NTRS1975

Application of a transitional boundary-layer theory in the low hypersonic Mach number regime

An investigation is made to assess the capability of a finite-difference boundary-layer procedure to predict the mean profile development across a transition from laminar to turbulent flow in the low hypersonic Mach-number regime. The boundary-layer procedure uses an integral form of the turbulence kinetic-energy equation to govern the development of the Reynolds apparent shear stress. The present investigation shows the ability of this procedure to predict Stanton number, velocity profiles, and density profiles through the transition region and, in addition, to predict the effect of wall cooling and Mach number on transition Reynolds number. The contribution of the pressure-dilatation term to the energy balance is examined and it is suggested that transition can be initiated by the direct absorption of acoustic energy even if only a small amount (1 per cent) of the incident acoustic energy is absorbed.

Shamroth, S. J.↗
NASA NTRS1984

A global low order spectral model designed for climate sensitivity studies

A two level, global, spectral model using pressure as a vertical coordinate is developed. The system of equations describing the model is nonlinear and quasi-geostrophic. A moisture budget is calculated in the lower layer only with moist convective adjustment between the two layers. The mechanical forcing of topography is introduced as a lower boundary vertical velocity. Solar forcing is specified assuming a daily mean zenith angle. On land and sea ice surfaces a steady state thermal energy equation is solved to calculate the surface temperature. Over the oceans the sea surface temperatures are prescribed from the climatological average of January. The model is integrated to simulate the January climate.

Hanna, A. F.↗
NASA NTRS1991

Combustion wave instabilities in turbulent flows

Navier-Stokes solutions are used to calculate oscillatory components of pressure, velocity, and density, which in turn provide necessary data to compute energy growth factors to determine combustion instability. It is shown that wave instabilities are associated with changes in entropy and the space and time averages of oscillatory components of pressure, velocity and density, together with the mean flow field in the energy equation. Compressible laminar and turbulent flows and reacting flows with hydrogen/oxygen combustion are considered in this paper. The space shuttle main engine combustion/thrust chamber is used for illustration of the theory. The analysis shows that the increase of mean pressure and disturbances consistently results in the increase in instability. It is shown that adequate combustion instability analysis requires at least third order nonlinearity in energy growth or decay.

Yoon, W. S.↗
NASA NTRS1988

A mathematical model of the middle and high latitude ionosphere

A time-dependent three-dimensional model of the middle and high latitude ionosphere is described. The density distributions of six ion species NO(+), O(2+), N(2+), O(+), N(+), He(+), and the electron and ion temperatures are obtained from a numerical solution of the appropriate continuity, momentum, and energy equations. The equations are solved as a function of height for an inclined magnetic field at E and F region altitudes. The three-dimensional nature of the model is obtained by following flux tubes of plasma as they convect or corotate through a moving neutral atmosphere. The model takes account of field-aligned diffusion, cross-field electrodynamic drifts, thermospheric winds, polar wind escape, energy-dependent chemical reactions, neutral composition changes, ion production due to solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow and local heating and cooling processes. The model also takes account of the offset between the geomagnetic and geographic poles.

Schunk, R. W.↗
NASA NTRS1989

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.↗
DOE OSTI2023

Connecting primordial gravitational waves and dark energy

Cosmic acceleration manifested in the early universe as inflation, generating primordial gravitational waves detectable in the cosmic microwave background (CMB) radiation. Cosmic acceleration is occurring again at present as dark energy, detectable in cosmic distance and structure surveys. We explore the intriguing idea of connecting the two occurrences through quintessential inflation by an α-attractor potential without a cosmological constant. For this model we demonstrate robustness of the connection 1 + w 0 ≈ 4/(3N 2 r) between the present day dark energy equation of state parameter w 0 and the primordial tensor to scalar ratio r for a wide range of initial conditions. Analytic and numerical solutions produce current thawing behavior, resulting in a tight relation w a ≈ -1.53(1 + w 0 )≈ -0.2 (4 × 10 -3 /r). Upcoming CMB and galaxy redshift surveys can test this consistency condition. In conclusion, within this model, lack of detection of a dark energy deviation from Λ predicts a higher r, and lack of detection of r predicts greater dark energy dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Self-similar Lagrangian hydrodynamics of beam-heated solar flare atmospheres

The one-dimensional hydrodynamic problem in Lagrangian coordinates (Y, t) is considered for which the specific energy input Q has a power-law dependence on both Y and t, and the initial density distribution is rho(0) which is directly proportional to Y exp gamma. In regimes where the contributions of radiation, conduction, quiescent heating, and gravitational terms in the energy equation are negligible compared to those arising from Q, the problem has a self-similar solution, with the hydrodynamic variables depending only on a single independent variable which is a combination of Y, t, and the dimensional constants of the problem. It is then shown that the problem of solar flare chromospheric heating due to collisional interaction of a beam of electrons (or protons) with a power-law energy spectrum can be approximated by such forms of Q(Y, t) and rho(0)(Y), and that other terms are negligible compared to Q over a restricted regime early in the flare.

Brown, John C.↗
NASA NTRS1986

An introduction to three-dimensional climate modeling

The development and use of three-dimensional computer models of the earth's climate are discussed. The processes and interactions of the atmosphere, oceans, and sea ice are examined. The basic theory of climate simulation which includes the fundamental equations, models, and numerical techniques for simulating the atmosphere, oceans, and sea ice is described. Simulated wind, temperature, precipitation, ocean current, and sea ice distribution data are presented and compared to observational data. The responses of the climate to various environmental changes, such as variations in solar output or increases in atmospheric carbon dioxide, are modeled. Future developments in climate modeling are considered. Information is also provided on the derivation of the energy equation, the finite difference barotropic forecast model, the spectral transform technique, and the finite difference shallow water waved equation model.

Washington, W. M.↗
NASA NTRS1971

Finite difference solution to the three-dimensional, incompressible thermal energy boundary-layer equation

An implicit numerical method has been adopted for the solution to the three-dimensional, energy boundary-layer equation. The energy equation is written in terms of a dimensionless temperature function, the relative stagnation-enthalpy difference, and transformed by the introduction of a Blasius-type transformation of coordinates as well as dimensionless stream functions. The method is applied to the problem consisting of an infinite cylinder joined with its axis perpendicular to a thin, flat, heated plate. A Prandtl number equal to one is simply considered.

Fillo, J. A.↗
NASA NTRS2018

Radiative Heating of Large Meteoroids During Atmospheric Entry

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries was developed, which allows the commonly assumed heat transfer coefficient of 0.1 to be assessed. This model uses chemically reacting computational fluid dynamics (CFD), coupled with radiation transport and surface ablation. Coupled radiation accounts for the impact of radiation on the flowfield energy equations, while coupled ablation explicitly models the injection of ablation products within the flowfield and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation is shown to reduce the stagnation point radiative heating by over 60%. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to cases with only coupled radiation. This large reduction is partially the result of the low ionization energies of meteoric ablation products relative to air species. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock-layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients, the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases, for both laminar and turbulent, which is significantly lower than the canonical value of 0:1. When the new heat transfer model is applied to a Tunguska-like 15 Mt entry, the effect of the new model is to lower the height of burst by up to 2 km, depending on assumed entry angle. This, in turn, results in a significantly larger ground damage footprint than when the canonical heating assumption is used.

Christopher O Johnston↗
NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was 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 distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗
DOE OSTI2025

Comparison of Artificial Compressibility Methods for Coupled Laminar Fluid Flow and Heat Transfer

Artificial compressibility (AC) methods for incompressible flows perturb the mass continuity equation by the time derivative of the pressure field to achieve direct coupling between pressure and velocity without the need for a projection method. Several variations of the original AC method have been proposed (e.g., the entropically damped artificial compressibility [EDAC] method) that can, in some scenarios, improve on the original model. In this paper, we methodically compare the results of the original AC model and the EDAC model for both steady-state and transient problems with coupling between the momentum and energy equations. We demonstrate that the EDAC model and the AC model yield similar numerical results, and the EDAC model is 10%-20% faster than the AC model for a given problem.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
NASA NTRS1982

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is paramterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations.

Hanna, A. F.↗
NASA NTRS1984

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is parameterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations. Previously announced in STAR as N83-14822

Hanna, A. F.↗
NASA NTRS1983

Short-term climatic fluctuations forced by thermal anomalies

A two level, global, spectral model using pressure as a vertical coordinate was developed. The system of equations describing the model is nonlinear and quasi-geostrophic (linear balance). Static stability is variable in the model. A moisture budget is calculated in the lower layer only. Convective adjustment is used to avoid supercritical temperature lapse rates. The mechanical forcing of topography is introduced as a vertical velocity at the lower boundary. Solar forcing is specified assuming a daily mean zenith angle. The differential diabatic heating between land and sea is parameterized. On land and sea ice surfaces, a steady state thermal energy equation is solved to calculate the surface temperature. On the oceans, the sea surface temperature is specified as the climatological average for January. The model is used to simulate the January, February and March circulations. Previously announced in STAR as N83-14822

Hanna, A. F.↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
DOE OSTI2023

Dark energy with a triplet of classical U(1) fields

Here, we present a new mechanism for cosmic acceleration consisting of a scalar field coupled to a triplet of classical U(1) gauge fields. The gauge fields are arranged in a homogeneous, isotropic configuration, with both electric- and magneticlike vacuum expectation values. The gauge fields provide a masslike term via a Chern-Simons interaction that suspends the scalar away from its potential minimum, thereby enabling potential-dominated evolution. We show this mechanism can drive a brief period of acceleration, such as dark energy, without the need for fine-tunings. We obtain simple analytic results for the dark energy equation of state and dependence on model parameters. In this model, the presence of the gauge field generically leads to a suppression of long-wavelength gravitational waves, with implications for the experimental search for cosmic microwave background B modes and direct detection of a stochastic gravitational wave background.

79 ASTRONOMY AND ASTROPHYSICS↗