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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Numerical algorithm conjugating steady and transient, separated, compressible flow and a solid body having arbitrarily distributed heat sources

A numerical algorithm that is second-order accurate in time has been developed for the conjugated problem of a separated, compressible flow field and a conductive solid body. The full two-dimensional time-dependent Navier-Stokes equations are coupled with the time-dependent energy equation for the solid body and are solved simultaneously. using implicit algorithms. The energy equation for the solid body may include arbitrarily distributed heat sources. The algorithm has been exmined for the case of two-dimensional supersonic compression-corner interaction, with a heat source embedded in the wall in the vicinity of the separation bubble and the attached boundary layer. The effect of the heat source on the flow field is studied for steady and transient cases.

Degani, D.↗
NASA NTRS1992

Convection treatment for high Rayleigh number, laminar, natural convection calculation

The problem of natural convection in a square cavity of a Boussinesq incompressible fluid is computed using a consistent central-difference approximation, a consistent second-order upwind scheme, and the hybrid scheme for the energy equation. It was found that, for this problem, the second-order upwind scheme yields satisfactory solutions and that an apparently low-order combination of schemes, with the second-order upwind scheme for the momentum equations and the hybrid scheme for the energy equation, can produce satisfactory predictions for this problem.

Shyy, Wei↗
NASA NTRS1991

Benchmark solutions for the galactic heavy-ion transport equations with energy and spatial coupling

Nontrivial benchmark solutions are developed for the galactic heavy ion transport equations in the straightahead approximation with energy and spatial coupling. Analytical representations of the ion fluxes are obtained for a variety of sources with the assumption that the nuclear interaction parameters are energy independent. The method utilizes an analytical LaPlace transform inversion to yield a closed form representation that is computationally efficient. The flux profiles are then used to predict ion dose profiles, which are important for shield design studies.

Ganapol, Barry D.↗
NASA NTRS1991

Time-accurate finite difference solutions to the incompressible Navier-Stokes/energy equations

Two new algorithms for solving coupled Navier-Stokes and energy equations are presented. The algorithms are compared with available Navier-Stokes solutions for the model problem of a driven-cavity-flow over a range of Reynolds numbers. It is noted that the algorithms represent two different implementations of the fractional step approach to the solution of the coupled Navier-Stokes and energy equations using the Boussinesq approximation.

Mirfakhraee, Ali↗
NASA NTRS1991

Radiation-Controlled, Opposed-Flow Flame Spread in a Microgravity Environment

The effects of surface and gas-phase radiation on the rate and the structure of laminar flame spread over thin fuels are investigated using a flame-spread model which consists of the continuity, momentum, species, and energy equations in the gas and the continuity and energy equations in the solid. Numerical calculations, complemented by scaling arguments, show that, at high velocities of the oxidizer flow, radiation effects are unimportant; the spread rate decreases with increasing opposing velocity due to finite-rate gas-phase kinetics. However, radiation becomes progressively important when the opposing velocity is below a certain value: the flame cools, shrinks in size, and its spread rate falls sharply with decreasing opposing velocity.

Bhattacharjee, Subrata↗
NASA NTRS1986

A fully coupled variable properties thermohydraulic model for a cryogenic hydrostatic journal bearing

The goal set forth here is to continue the work started by Braun et al. (1984-1985) and present an integrated analysis of the behavior of the two row, 20 staggered pockets, hydrostatic cryogenic bearing used by the turbopumps of the Space Shuttle main engine. The variable properties Reynolds equation is fully coupled with the two-dimensional fluid film energy equation. The three-dimensional equations of the shaft and bushing model the boundary conditions of the fluid film energy equation. The effects of shaft eccentricity, angular velocity, and inertia pressure drops at pocket edge are incorporated in the model. Their effects on the bearing fluid properties, load carrying capacity, mass flow, pressure, velocity, and temperature form the ultimate object of this paper.

Braun, M. J.↗
NASA NTRS1990

Numerical simulation of 3-D Benard convection with gravitational modulation

In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.

Biringen, S.↗
DOE OSTI2023

Drift-ordered fluid vorticity equation with energy consistency

Although drift-ordered fluid models are widely applied in tokamak edge turbulence simulations, the models used are acknowledged not to conserve energy or even electrical charge. The present paper aims to remove many of the existing pitfalls in drift-fluid models, however, with the objective of finding a solution simple enough to be implemented in numerical applications. Our main result is an improved version of the drift-Braginskii equations involving a generalized vorticity function. In the new drift-Braginskii system, the quasi-neutrality condition translates into a transport equation for a generalized vorticity, expressed in conservation form, and related to the total mass-weighted circulation. It is found that kinetic energy conservation can be achieved if the polarization flow is defined recursively. In conclusion, the resulting model conserves the kinetic energy associated with E×B and diamagnetic flows and retains the associated perpendicular kinetic energy flux.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Basic lubrication equations

Lubricants, usually Newtonian fluids, are assumed to experience laminar flow. The basic equations used to describe the flow are the Navier-Stokes equation of motion. The study of hydrodynamic lubrication is, from a mathematical standpoint, the application of a reduced form of these Navier-Stokes equations in association with the continuity equation. The Reynolds equation can also be derived from first principles, provided of course that the same basic assumptions are adopted in each case. Both methods are used in deriving the Reynolds equation, and the assumptions inherent in reducing the Navier-Stokes equations are specified. Because the Reynolds equation contains viscosity and density terms and these properties depend on temperature and pressure, it is often necessary to couple the Reynolds with energy equation. The lubricant properties and the energy equation are presented. Film thickness, a parameter of the Reynolds equation, is a function of the elastic behavior of the bearing surface. The governing elasticity equation is therefore presented.

Hamrock, B. J.↗
DOE OSTI2022

Toward full simulations for a liquid metal blanket: part 2. Computations of MHD flows with volumetric heating for a PbLi blanket prototype at Ha ~10 4 and Gr ~10 12

On the pathway toward full simulations for a liquid metal (LM) blanket, this part 2 extends a previous study of purely magnetohydrodynamic (MHD) flows in a DCLL blanket in reference Chen et al (2020 Nucl. Fusion 60 076003) to more general conditions when the MHD flow is coupled with heat transfer. The simulated prototypic blanket module includes all components of a real LM blanket system, such as supply ducts, inlet and outlet manifolds, multiple poloidal ducts and a U-turn zone. Volumetric heating generated by fusion neutrons is added to simulate thermal effects in the flowing lead–lithium (PbLi) breeder. The MHD flow equations and the energy equation are solved with a DNS-type finite-volume code ‘MHD-UCAS’ on a very fine mesh of 470 × 10 6 cells. The applied magnetic field is 5 T (Hartmann number Ha ~ 10 4 ), the PbLi velocity in the poloidal ducts is 10 cm s –1 (Reynolds number Re ~10 5 ), whereas the maximum volumetric heating is 30 MW m –3 (Grashof number Gr ~ 10 12 ). Four cases have been simulated, including forced- and mixed-convection flows, and either an electrically conducting or insulating blanket structure. Various comparisons are made between the four computed cases and also against the purely MHD flows computed earlier in reference Chen et al (2020 Nucl. Fusion 60 076003) with regards to the (1) MHD pressure drop, (2) flow balancing, (3) temperature field, (4) flows in particular blanket components, and (5) 3D and turbulent flow effects. The strongest buoyancy effects were found in the poloidal ducts. In the electrically non-conducting blanket, the buoyancy forces lead to significant modifications of the flow structure, such as formation of reverse flows, whereas their effect on the MHD pressure drop is relatively small. In the electrically conducting blanket case, the buoyancy effects on the flow and MHD pressure drop are almost negligible.

Physics↗
NASA NTRS1983

Thermal elastohydrodynamic lubrication of line contacts

A numerical solution to the problem of thermal elastohydrodynamic lubrication of line contacts was obtained by using a finite difference formulation. The solution procedure consists of simultaneous solution of the thermal Reynolds equation, the elasticity equation, and the energy equation subject to appropriate boundary conditions. Pressure distribution, film shape, and temperature distribution were obtained for fully flooded conjunctions, a paraffinic lubricant, and various dimensionless speed parameters while the dimensionless load and materials parameters were held constant. Reduction in the minimum film thickness due to thermal effects (as a ratio of thermal to isothermal minimum film thickness) is given by a simple formula as a function of the thermal loading parameter Q: H(min)/H(min,I) = 10/10+ Q(0.4). Plots of pressure distribution, film shape, temperature distribution, and flow are shown for some representative cases.

Ghosh, M. K.↗
NASA NTRS1985

Thermal elastohydrodynamic lubrication of line contacts

A numerical solution to the problem of thermal elastohydrodynamic lubrication of line contacts was obtained by using a finite difference formulation. The solution procedure consists of simultaneous solution of the thermal Reynolds equation, the elasticity equation, and the energy equation subject to appropriate boundary conditions. Pressure distribution, film shape, and temperature distribution were obtained for fully flooded conjunctions, a paraffinic lubricant, and various dimensionless speed parameters while the dimensionless load and materials parameters were held constant. Reduction in the minimum film thickness due to thermal effects (as a ratio of thermal to isothermal minimum film thickness) is given by a simple formula as a function of the thermal loading parameter Q: H(min)/H(min, I) = 10/10 + Q(0.4). Plots of pressure distribution, film shape, temperature distribution, and flow are shown for some representative cases.

Ghosh, M. K.↗
NASA NTRS1989

Computational study of 3-D Benard convection with gravitational modulation

In this numerical study the effects of a modulated gravitational field on three-dimensional Rayleigh-Benard convection with heating from above or from below is investigated. The full, nonlinear, time-dependent, Boussinesq Navier-Stokes equations and the energy equation are solved by a semiimplicit, pseudo-spectral procedure. This study has been motivated by the need to better understand the effects of vibration (G-Jitter) on fluids systems especially in the low gravity environment.

Biringen, S.↗
NASA NTRS1990

Combustion instability analysis

A theory and computer program for combustion instability analysis are presented. The basic theoretical foundation resides in the concept of entropy-controlled energy growth or decay. Third order perturbation expansion is performed on the entropy-controlled acoustic energy equation to obtain the first order integrodifferential equation for the energy growth factor in terms of the linear, second, and third order energy growth parameters. These parameters are calculated from Navier-Stokes solutions with time averages performed on as many Navier-Stokes time steps as required to cover at least one peak wave period. Applications are made for a 1-D Navier-Stokes solution for the Space Shuttle Main Engine (SSME) thrust chamber with cross section area variations taken into account. It is shown that instability occurs when the mean pressure is set at 2000 psi with 30 percent disturbances. Instability also arises when the mean pressure is set at 2935 psi with 20 percent disturbances. The system with mean pressures and disturbances more adverse that these cases were shown to be unstable.

Chung, T. J.↗
NASA NTRS1993

Time-dependent thermocapillary convection in a Cartesian cavity - Numerical results for a moderate Prandtl number fluid

The present numerical simulation explores a thermal-convective mechanism for oscillatory thermocapillary convection in a shallow Cartesian cavity for a Prandtl number 6.78 fluid. The computer program developed for this simulation integrates the two-dimensional, time-dependent Navier-Stokes equations and the energy equation by a time-accurate method on a stretched, staggered mesh. Flat free surfaces are assumed. The instability is shown to depend upon temporal coupling between large scale thermal structures within the flow field and the temperature sensitive free surface. A primary result of this study is the development of a stability diagram presenting the critical Marangoni number separating steady from the time-dependent flow states as a function of aspect ratio for the range of values between 2.3 and 3.8. Within this range, a minimum critical aspect ratio near 2.3 and a minimum critical Marangoni number near 20,000 are predicted below which steady convection is found.

Peltier, L. J.↗
NASA NTRS1996

Three-Dimensional Simulations of Marangoni-Benard Convection in Small Containers by the Least-Squares Finite Element Method

This paper reports a numerical study of the Marangoni-Benard (MB) convection in a planar fluid layer. The least-squares finite element method (LSFEM) is employed to solve the three-dimensional Stokes equations and the energy equation. First, the governing equations are reduced to be first-order by introducing variables such as vorticity and heat fluxes. The resultant first-order system is then cast into a div-curl-grad formulation, and its ellipticity and permissible boundary conditions are readily proved. This numerical approach provides an equal-order discretization for velocity, pressure, vorticity, temperature, and heat conduction fluxes, and therefore can provide high fidelity solutions for the complex flow physics of the MB convection. Numerical results reported include the critical Marangoni numbers (M(sub ac)) for the onset of the convection in containers with various aspect ratios, and the planforms of supercritical MB flows. The numerical solutions compared favorably with the experimental results reported by Koschmieder et al..

Yu, Sheng-Tao↗
NASA NTRS1972

Analysis of heat transfer in a porous cooled wall with variable pressure and temperature along the coolant exit boundary

Fluid from a reservior at constant pressure and temperature is forced through a porous wall of uniform thickness. The boundary through which the fluid exits has specified variations in pressure and temperature along it in one direction so that the flow and heat transfer are two-dimensional. The local fluid and matrix temperatures are assumed to be equal and therefore a single energy equation governs the temperature distribution within the wall. The solution is obtained by transforming this energy equation into potential plane coordinates, which results in a separable equation. A technique yielding an integral equation is used to adapt the general solution so that it satisfies the variable-pressure boundary condition. Analytical expressions are given for the normal exit velocity and heat flux along the exit boundary. Illustrative examples are carried out which indicate to what extent the solution is locally one-dimensional.

Siegel, R.↗
NASA NTRS1982

'Scaling' analysis of the ice accretion process on aircraft surfaces

A comprehensive set of scaling parameters is developed for the ice accretion process by analyzing the energy equations of the dynamic freezing zone and the already frozen ice layer, the continuity equation associated with supercooled liquid droplets entering into and impacting within the dynamic freezing zone, and energy equation of the ice layer. No initial arbitrary judgments are made regarding the relative magnitudes of each of the terms. The method of intrinsic reference variables in employed in order to develop the appropriate scaling parameters and their relative significance in rime icing conditions in an orderly process, rather than utilizing empiricism. The significance of these parameters is examined and the parameters are combined with scaling criteria related to droplet trajectory similitude.

Keshock, E. G.↗