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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Analytic Modeling of the Hydrodynamic, Thermal, and Structural Behavior of Foil Thrust Bearings

A simulation and modeling effort is conducted on gas foil thrust bearings. A foil bearing is a self acting hydrodynamic device capable of separating stationary and rotating components of rotating machinery by a film of air or other gaseous lubricant. Although simple in appearance these bearings have proven to be complicated devices in analysis. They are sensitive to fluid structure interaction, use a compressible gas as a lubricant, may not be in the fully continuum range of fluid mechanics, and operate in the range where viscous heat generation is significant. These factors provide a challenge to the simulation and modeling task. The Reynolds equation with the addition of Knudsen number effects due to thin film thicknesses is used to simulate the hydrodynamics. The energy equation is manipulated to simulate the temperature field of the lubricant film and combined with the ideal gas relationship, provides density field input to the Reynolds equation. Heat transfer between the lubricant and the surroundings is also modeled. The structural deformations of the bearing are modeled with a single partial differential equation. The equation models the top foil as a thin, bending dominated membrane whose deflections are governed by the biharmonic equation. A linear superposition of hydrodynamic load and compliant foundation reaction is included. The stiffness of the compliant foundation is modeled as a distributed stiffness that supports the top foil. The system of governing equations is solved numerically by a computer program written in the Mathematica computing environment. Representative calculations and comparisons with experimental results are included for a generation I gas foil thrust bearing.

Bruckner, Robert J.↗
NASA NTRS1993

Comparisons between DSMC and the Navier-Stokes equations for reentry flows

A detailed comparison is made between Navier-Stokes and DSMC calculations for flows near the continuum limit to assess the accuracy of the continuum equations in this regime. Meaningful comparisons require the use of similar physical models. This necessitates the inclusion of a separate rotational energy equation and use of slip boundary conditions. Inclusion of slip boundary conditions resulted in improved agreement between surface properties. Moreover, good agreement was obtained for the various temperatures in the nonequilibrium portion of the flow field that does not contain the shock region. Departures are noted in the shock region and in regions where thermal diffusion effects are important.

Olynick, David R.↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

A study of the interactions between large-scale coherent structures and fine-grained turbulence in a round jet

An approximate energy integral description is used to study the development of large-scale coherent structures in the technologically important problem of the round turbulent jet. The analysis begins from the radially integrated form of the kinetic energy equations of the mean flow, the large-scale structure and the fine-grained turbulence, which are obtained by using the usual Reynolds time average and a conditional average with reference to the frequency of the idealized monochromatic component of the large-scale wavelike structure. This is the basis for obtaining the amplitude equations for the three components of the flow in terms of the mean flow momentum thickness, the large-scale structure kinetic energy and the fine-grained turbulence kinetic energy across the jet.

Mankbadi, R.↗
NASA NTRS1996

Two-Flux Green's Function Analysis for Transient Spectral Radiation in a Composite

An analysis is developed for obtaining transient temperatures in a two-layer semitransparent composite with spectrally dependent properties. Each external boundary of the composite is subjected to radiation and convection. The two-flux radiative transfer equations are solved by deriving a Green's function. This yields the local radiative heat source needed to numerically solve the transient energy equation. An advantage of the two-flux method is that isotropic scattering is included without added complexity. The layer refractive indices are larger than one. This produces internal reflections at the boundaries and the internal interface; the reflections are assumed diffuse. Spectral results using the Green's function method are verified by comparing with numerical solutions using the exact radiative transfer equations. Transient temperature distributions are given to illustrate the effect of radiative heating on one side of a composite with external convective cooling. The protection of a material from incident radiation is illustrated by adding scattering to the layer adjacent to the radiative source.

Siegel, Robert↗
NASA NTRS1999

Steady Secondary Flows Generated by Periodic Compression and Expansion of an Ideal Gas in a Pulse Tube

This study establishes a consistent set of differential equations for use in describing the steady secondary flows generated by periodic compression and expansion of an ideal gas in pulse tubes. Also considered is heat transfer between the gas and the tube wall of finite thickness. A small-amplitude series expansion solution in the inverse Strouhal number is proposed for the two-dimensional axisymmetric mass, momentum and energy equations. The anelastic approach applies when shock and acoustic energies are small compared with the energy needed to compress and expand the gas. An analytic solution to the ordered series is obtained in the strong temperature limit where the zeroth-order temperature is constant. The solution shows steady velocities increase linearly for small Valensi number and can be of order I for large Valensi number. A conversion of steady work flow to heat flow occurs whenever temperature, velocity or phase angle gradients are present. Steady enthalpy flow is reduced by heat transfer and is scaled by the Prandtl times Valensi numbers. Particle velocities from a smoke-wire experiment were compared with predictions for the basic and orifice pulse tube configurations. The theory accurately predicted the observed steady streaming.

Lee, Jeffrey M.↗
NASA NTRS1973

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction, volume 1

The turbulent kinetic energy equation is coupled with boundary layer equations to solve the characteristics of compressible turbulent boundary layers with mass injection and combustion. The Reynolds stress is related to the turbulent kinetic energy using the Prandtl-Wieghardt formulation. When a lean mixture of hydrogen and nitrogen is injected through a porous plate into the subsonic turbulent boundary layer of air flow and ignited by external means, the turbulent kinetic energy increases twice as much as that of noncombusting flow with the same mass injection rate of nitrogen. The magnitudes of eddy viscosity between combusting and noncombusting flows with injection, however, are almost the same due to temperature effects, while the distributions are different. The velocity profiles are significantly affected by combustion; that is, combustion alters the velocity profile as if the mass injection rate is increased, reducing the skin-friction as a result of a smaller velocity gradient at the wall. If pure hydrogen as a transpiration coolant is injected into a rocket nozzle boundary layer flow of combustion products, the temperature drops significantly across the boundary layer due to the high heat capacity of hydrogen. At a certain distance from the wall, hydrogen reacts with the combustion products, liberating an extensive amount of heat. The resulting large increase in temperature reduces the eddy viscosity in this region.

Omori, S.↗
NASA NTRS1973

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction

The turbulent kinetic energy equation is coupled with boundary layer equations to solve the characteristics of compressible turbulent boundary layers with mass injection and combustion. The Reynolds stress is related to the turbulent kinetic energy using the Prandtl-Wieghardt formulation. When a lean mixture of hydrogen and nitrogen is injected through a porous plate into the subsonic turbulent boundary layer of air flow and ignited by external means, the turbulent kinetic energy increases twice as much as that of noncombusting flow with the same mass injection rate of nitrogen. The magnitudes of eddy viscosity between combusting and noncombusting flows with injection, however, are almost the same due to temperature effects, while the distributions are different. The velocity profiles are significantly affected by combustion. If pure hydrogen as a transpiration coolant is injected into a rocket nozzle boundary layer flow of combustion products, the temperature drops significantly across the boundary layer due to the high heat capacity of hydrogen. At a certain distance from the wall hydrogen reacts with the combustion products, liberating an extensive amount of heat.

Omori, S.↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗
NASA NTRS1976

Time-asymptotic solutions of the Navier-Stokes equation for free shear flows using an alternating-direction implicit method

An uncoupled time asymptotic alternating direction implicit method for solving the Navier-Stokes equations was tested on two laminar parallel mixing flows. A constant total temperature was assumed in order to eliminate the need to solve the full energy equation; consequently, static temperature was evaluated by using algebraic relationship. For the mixing of two supersonic streams at a Reynolds number of 1,000, convergent solutions were obtained for a time step 5 times the maximum allowable size for an explicit method. The solution diverged for a time step 10 times the explicit limit. Improved convergence was obtained when upwind differencing was used for convective terms. Larger time steps were not possible with either upwind differencing or the diagonally dominant scheme. Artificial viscosity was added to the continuity equation in order to eliminate divergence for the mixing of a subsonic stream with a supersonic stream at a Reynolds number of 1,000.

Rudy, D. H.↗
NASA NTRS1992

Electrons in the ionosphere of Titan

A theoretical model of the spatial and energy distribution of electrons in the ionosphere of Titan has been constructed using the two-stream electron transport method and the electron energy equation. The calculated electron spectra show abrupt decreases that can contribute to the 'bite-out' signature observed by the Voyager 1 plasma science instrument. Energy deposition rates in the exosphere of Titan by photoelectrons and magnetospheric electrons are calculated. Calculated N2 EUV airglow emission rates lead to the conclusion that airglow emission due to photoelectron impact is much more important than airglow emission due to magnetospheric electron interactions. Thermal electron temperatures at Titan are calculated for the first time. The electron gas remains well thermalized with the neutral atmosphere for radial distances from the center of Titan less than 3500 km. For radial distances beyond about 4000 km, energy transport dominates the energetics, and the electrons are almost isothermal along magnetic field lines.

Gan, L.↗
NASA NTRS1995

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Blaisdell, Gregory A.↗
NASA NTRS1975

Chemically Reacting One-Dimensional Gas-Particle Flows

The governing equations for the one-dimensional flow of a gas-particle system are discussed. Gas-particle effects are coupled via the system momentum and energy equations with the gas assumed to be chemically frozen or in chemical equilibrium. A computer code for calculating the one-dimensional flow of a gas-particle system is discussed and a user's input guide presented. The computer code provides for the expansion of the gas-particle system from a specified starting velocity and nozzle inlet geometry. Though general in nature, the final output of the code is a startline for initiating the solution of a supersonic gas-particle system in rocket nozzles. The startline includes gasdynamic data defining gaseous startline points from the nozzle centerline to the nozzle wall and particle properties at points along the gaseous startline.

Tevepaugh, J. A.↗
NASA NTRS1981

An improved model for the combustion of AP composite propellants

This paper presents several improvements to the BDP model of steady-state burning of AP composite solid propellants. The Price-Boggs-Derr model of AP monopropellant burning is incorporated to represent the AP. A separate energy equation is written for the binder to permit a different surface temperature from the AP; this includes an analysis of the sharing of primary diffusion flame energy, and correction of a BDP model inconsistency in treating the binder regression rate. A method for assembling component contributions to calculate the burning rates of multimodal propellants is also presented. Results are shown in the form of representative burning rate curves, comparisons with data, and calculated internal details of interest. Ideas for future work are discussed in an Appendix.

Cohen, N. S.↗
NASA NTRS1993

Shear and compression buckling analysis for anisotropic panels with centrally located elliptical cutouts

An approximate analysis for buckling of biaxial- and shear-loaded anisotropic panels with centrally located elliptical cutouts is presented in the present paper. The analysis is composed of two parts, a prebuckling analysis and a buckling analysis. The prebuckling solution is determined using Lekhnitskii's complex variable equations of plane elastostatics combined with a Laurent series approximation and a boundary collocation method. The buckling solution is obtained using the principle of minimum potential energy. A by-product of the minimum potential energy equation is an integral equation which is solved using Gaussian quadrature. Comparisons with documented experimental results and finite element analyses indicate that the approximate analysis accurately predicts the buckling loads of square biaxial- and shear-loaded panels having elliptical cutouts with major axes up to sixty percent of the panel width. Results of a parametric study are presented for shear- and compression-loaded rectangular anisotropic panels with elliptical cutouts. The effects of panel aspect ratio, cutout shape, cutout size, cutout orientation, laminate anisotropy, and combined loading on the buckling load are examined.

Britt, V. O.↗
DOE OSTI2023

Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2

Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.

Meteorology & Atmospheric Sciences↗
NASA NTRS1986

Three-dimensional computation of mixing of transverse injector in a ducted supersonic airstream

Numerical solutions of the three-dimensional mass-averaged elliptic Navier-Stokes equations, including species transport, are obtained for nonreacting, turbulent, mixing flow fields for the case of transverse sonic injection of a secondary gas into a supersonic airstream through a circular orifice injector. Results are presented for flow through a constant area duct and through a duct with a rearward-facing step upstream of the injector. The equations are numerically integrated using MacCormack's explicit method and turbulence is included using the Baldwin-Lomax algebraic eddy viscosity model. In the species transport and energy equations, diffusion coefficients based on Fick's law and an assumption of unit Lewis number are applied. The computations were performed on a CDC-VPS-32 (extended version of Cyber-205) using a grid consisting of approximately 200,000 points. The computed results are compared with experimentally observed penetration and spreading boundaries for an injected gas at two dynamic pressure ratios. Three-dimensional flow field structures are dipicted in terms of static pressure, mass fractions of species and velocity vectors.

Uenishi, K.↗