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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A near-wall turbulence model and its application to fully developed turbulent channel and pipe flows

A near wall turbulence model and its incorporation into a multiple-time-scale turbulence model are presented. In the method, the conservation of mass, momentum, and the turbulent kinetic energy equations are integrated up to the wall; and the energy transfer rate and the dissipation rate inside the near wall layer are obtained from algebraic equations. The algebraic equations for the energy transfer rate and the dissipation rate inside the near wall layer were obtained from a k-equation turbulence model and the near wall analysis. A fully developed turbulent channel flow and fully developed turbulent pipe flows were solved using a finite element method to test the predictive capability of the turbulence model. The computational results compared favorably with experimental data. It is also shown that the present turbulence model could resolve the over shoot phenomena of the turbulent kinetic energy and the dissipation rate in the region very close to the wall.

Kim, S.-W.↗
NASA NTRS2002

Modeling Disturbance Dynamics in Transitional and Turbulent Boundary Layers

The dynamics of an ensemble of linear disturbances in boundary-layer flows at various Reynolds numbers is studied through an analysis of the transport equations for the mean disturbance kinetic energy and energy dissipation rate. Effects of adverse and favorable pressure-gradients on the disturbance dynamics are also included in the analysis. Unlike the fully turbulent regime where nonlinear phase scrambling of the fluctuations affects the flow field even in proximity to the wall, the early stage transition regime fluctuations studied here are influenced across the boundary layer by the solid boundary. In addition, the dominating dynamics in the disturbance kinetic energy equation is governed by the energy production, pressure-transport and viscous diffusion - also in contrast to the fully turbulent regime. For the disturbance dissipation rate, a dynamic balance exists between the destruction and diffusion of dissipation.

Grosch, C. E.↗
NASA NTRS1980

Mass transfer between binary stars

The transfer of mass from one component of a binary system to another by mass ejection is analyzed through a stellar wind mechanism, using a model which integrates the equations of motion, including the energy equation, with an initial static atmosphere and various temperature fluctuations imposed at the base of the star's corona. The model is applied to several situations and the energy flow is calculated along the line of centers between the two binary components, in the rotating frame of the system, thereby incorporating the centrifugal force. It is shown that relatively small disturbances in the lower chromosphere or photosphere can produce mass loss through a stellar wind mechanism, due to the amplification of the disturbance propagating into the thinner atmosphere. Since there are many possible sources of the disturbance, the model can be used to explain many mass ejection phenomena.

Modisette, J. L.↗
NASA NTRS2023

Simulations of Zero Boiloff, Densification, and Solidification in a Large-Scale Liquid Hydrogen Tank

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Colin Philip Mahony↗
DOE OSTI2026

Dark Energy Survey: Implications for cosmological expansion models from the final DES baryon acoustic oscillation and supernova data

The Dark Energy Survey (DES) recently released the final results of its two principal probes of the expansion history: Type Ia supernovae (SNe) and baryonic acoustic oscillations (BAO). In this paper, we explore the cosmological implications of these data in combination with external cosmic microwave background (CMB), big bang nucleosynthesis (BBN), and age-of-the-Universe information. The BAO measurement, which is ∼ 2 σ away from Planck ’s Λ CDM predictions, pushes for low values of Ω m compared to Planck, in contrast to SN which prefers a higher value than Planck. We identify several tensions among datasets in the Λ CDM model that cannot be resolved by including either curvature ( k Λ CDM ) or a constant dark energy equation of state ( w CDM ). By combining BAO + SN + CMB despite these mild tensions, we obtain Ω k = - 5.5 - 4.2 + 4.6 × 10 - 3 in k Λ CDM , and w = - 0.94 8 - 0.027 + 0.028 in w CDM . In w CDM , BAO and SN push again in different directions of parameter space, favoring, respectively, w < - 1 and w > - 1 . If we open the parameter space to w 0 w a CDM [where the equation of state of dark energy varies as w ( a ) = w 0 + ( 1 - a ) w a ], all the datasets are mutually more compatible, and we find concordance in the [ w 0 > - 1 , w a < 0 ] quadrant, with BAO pushing for w a < 0 and SN for [ w 0 > - 1 , w a < 0 ] . For DES BAO and SN in combination with Planck -CMB, we find a 3.2 σ deviation from Λ CDM , with w 0 = - 0.67 3 - 0.097 + 0.098 , w a = - 1.3 7 - 0.50 + 0.51 , a Hubble constant of H 0 = 67.8 1 - 0.86 + 0.96 km s - 1 Mpc - 1 , and an abundance of matter of Ω m = 0.310 9 - 0.0099 + 0.0086 . For the combination of all the background cosmological probes considered (including CMB’s angular acoustic scale θ ⋆ ), we still find a deviation of 2.8 σ from Λ CDM in the w 0 - w a plane. Assuming a minimal neutrino mass, this work provides tentative evidence for non- Λ CDM physics, which is consistent with recent claims in support of evolving dark energy, or a source of unknown systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1975

Shock wave relations in lunar ash flow

A detailed analysis of steady normal shock waves in a layer of lunar ash by the theory of two phase flow of a mixture of a gas and small solid particles is presented. New terms of pressure gradient and virtual mass forces in the particle momentum equation and the particle kinetic energy and work done in the particle energy equation are included in the fundamental equations and their influences are investigated. The flow variables as functions of the free stream Mach number, initial particle volume fraction and the density ratio G are presented. The thickness of the relaxation zone is found to increase with decreasing initial particle volume fraction and almost independent of G for any given values of free stream Mach number. A new empirical relation of the drag coefficient for the spherical particles in the mixture is proposed to fit experimental data.

Pai, S. I.↗
NASA NTRS1978

Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities

An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.

Anderson, O. L.↗
NASA NTRS1979

Computational simulation of turbulent vortex merger and decay

The interaction and eventual merger of corotational vortices and the decay of a single vortex have been studied by employing zero-, one- and two-equation turbulent-flow models in order to gain a better understanding of the role of turbulence. An implicit finite-difference procedure is used to integrate the unsteady, two-dimensional equations in a cross-plane. The zero- and one-equation formulations utilize a mixing-length model, which incorporates the streamline curvature effect by prescribing a spatially-varying mixing-length. In the two-equation model, the turbulence kinetic energy equation and a modified rate of dissipation equation which includes a streamline curvature correction are solved. Computational results of different models applied to various flow-configurations are presented and compared with available experimental data whenever possible.

Raj, P.↗
NASA NTRS1983

Thermodynamic analysis of journal bearings operating under steady state loading in laminar regime

Thermohydrodynamic effects in journal bearings operating under steady load in laminar regime are investigated. An analytical model for the finite and infinitely long journal bearings is formulated. The model includes correction factors for the cavitation effects in the unloaded region of the bearing and the mixing of the recirculating oil and supply oil at the oil inlet. A finite difference computer program is developed to numerically solve the governing equations of the continuity, Reynolds, energy, Laplace heat conduction, and a viscosity-temperature relation simultaneously. The program includes a numerical technique for obtaining an isothermal shaft temperature. The numerical results of temperature distribution and the heat effects on the bearing load carrying capacity agree closely with those of experimental findings. Several different sets of simpler boundary conditions for the energy equation are studied.

Khonsari, M. M.↗
NASA NTRS1990

Modeling of near-wall turbulence

An improved k-epsilon model and a second order closure model is presented for low Reynolds number turbulence near a wall. For the k-epsilon model, a modified form of the eddy viscosity having correct asymptotic near wall behavior is suggested, and a model for the pressure diffusion term in the turbulent kinetic energy equation is proposed. For the second order closure model, the existing models are modified for the Reynolds stress equations to have proper near wall behavior. A dissipation rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The calculations are compared with direct numerical simulations. It is shown that the present models, both the k-epsilon model and the second order closure model, perform well in predicting the behavior of the near wall turbulence. Significant improvements over previous models are obtained.

Shih, T. H.↗
NASA NTRS1992

The Proteus Navier-Stokes code

An effort is currently underway at NASA Lewis to develop two and three dimensional Navier-Stokes codes, called Proteus, for aerospace propulsion applications. Proteus solves the Reynolds-averaged, unsteady, compressible Navier-Stokes equations in strong conservation law form. Turbulence is modeled using a Baldwin-Lomax based algebraic eddy viscosity model. In addition, options are available to solve thin layer or Euler equations, and to eliminate the energy equation by assuming constant stagnation enthalpy. An extensive series of validation cases have been run, primarily using the two dimensional planar/axisymmetric version of the code. Several flows were computed that have exact solution such as: fully developed channel and pipe flow; Couette flow with and without pressure gradients; unsteady Couette flow formation; flow near a suddenly accelerated flat plate; flow between concentric rotating cylinders; and flow near a rotating disk. The two dimensional version of the Proteus code has been released, and the three dimensional code is scheduled for release in late 1991.

Towne, Charles E.↗
NASA NTRS2004

Procedure for Determining One-Dimensional Flow Distributions in Arbitrarily Connected Passages Without the Influence of Pumping

A calculation procedure is presented which allows the one-dimensional determination of flow distributions in arbitrarily connected (branching) flow passages having multiple inlets and exits. The procedure uses an adaptation of the finite element technique, iteratively coupled with an accurate one-dimensional flow solver. The procedure eliminates the usual restrictions inherent with finite element flow calculations. Unlike existing one-dimensional methods, which require simplifications to the flow equations (uncoupling the momentum and energy equations), to allow for arbitrary branching and multiple inlets and exits, the only limitation of the described methodology is that, at present, it can only accommodate non-rotating configurations (no pumping effects). The calculation procedure is robust, and will always converge for physically possible flow. The procedure is described, and its use is illustrated by an example.

Meitner, Peter L.↗
NASA NTRS2017

The Root Cause of the Overheating Problem

Previously we identified the receding flow, where two fluid streams recede from each other, as an open numerical problem, because all well-known numerical fluxes give an anomalous temperature rise, thus called the overheating problem. This phenomenon, although presented in several textbooks, and many previous publications, has scarcely been satisfactorily addressed and the root cause of the overheating problem not well understood. We found that this temperature rise was solely connected to entropy rise and proposed to use the method of characteristics to eradicate the problem. However, the root cause of the entropy production was still unclear. In the present study, we identify the cause of this problem: the entropy rise is rooted in the pressure flux in a finite volume formulation and is implanted at the first time step. It is found theoretically inevitable for all existing numerical flux schemes used in the finite volume setting, as confirmed by numerical tests. This difficulty cannot be eliminated by manipulating time step, grid size, spatial accuracy, etc, although the rate of overheating depends on the flux scheme used. Finally, we incorporate the entropy transport equation, in place of the energy equation, to ensure preservation of entropy, thus correcting this temperature anomaly. Its applicability is demonstrated for some relevant 1D and 2D problems. Thus, the present study validates that the entropy generated ab initio is the genesis of the overheating problem.

Liou, Meng-Sing↗
NASA NTRSDate not supplied

Energy of Fluid (EOF) Simulation of Large-Scale Densification and Solidification of Hydrogen

The current research develops a finite volume-based Computation Fluid Dynamic (CFD) model utilizing an Energy of Fluid (EOF) approach to generate a simulation of the densification of liquid hydrogen (LH2) in the Integrated Refrigeration and Storage (IRAS) experimental tank for the Ground Operations Demonstration Unit for Liquid Hydrogen Project (GODU-LH2) at the Kennedy Space Center (KSC). The computational code will integrate a commercial software pressure-based model with User Defined Functions (UDF) generated and added to the model. The modified model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. This method will decrease computational time when compared to the built-in commercially available Enthalpy solvers. Currently the research focused on utilizing User Defined Scalars (UDS) compared to the enthalpy model and a lumped node analysis to show the written UDS is as accurate as the built-in solvers, and the models will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

C. Mahony↗