IMPLEMENTATION OF THE ENERGY EQUATION SOLVER TO THE LATTICE BOLTZMANN METHOD-BASED CODE PRATHAM
Explore the source record for details and available documents.
Jarrah, Ibrahim↗
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.
Analytical study of the liquid phase transient behavior of a high temperature heat pipe
The transient operation of the liquid phase of a high temperature heat pipe is studied. The study was conducted in support of advanced heat pipe applications that require reliable transport of high temperature drops and significant distances under a broad spectrum of operating conditions. The heat pipe configuration studied consists of a sealed cylindrical enclosure containing a capillary wick structure and sodium working fluid. The wick is an annular flow channel configuration formed between the enclosure interior wall and a concentric cylindrical tube of fine pore screen. The study approach is analytical through the solution of the governing equations. The energy equation is solved over the pipe wall and liquid region using the finite difference Peaceman-Rachford alternating direction implicit numerical method. The continuity and momentum equations are solved over the liquid region by the integral method. The energy equation and liquid dynamics equation are tightly coupled due to the phase change process at the liquid-vapor interface. A kinetic theory model is used to define the phase change process in terms of the temperature jump between the liquid-vapor surface and the bulk vapor. Extensive auxiliary relations, including sodium properties as functions of temperature, are used to close the analytical system. The solution procedure is implemented in a FORTRAN algorithm with some optimization features to take advantage of the IBM System/370 Model 3090 vectorization facility. The code was intended for coupling to a vapor phase algorithm so that the entire heat pipe problem could be solved. As a test of code capabilities, the vapor phase was approximated in a simple manner.
Roche, Gregory Lawrence↗
Creeping flow solution of the Leidenfrost phenomenon
Creeping flow solution of Leidenfrost phenomenon by use of Navier-Stokes, continuity, and energy equations
NAVIER-STOKES EQUATION↗
On the Lagrangian description of unsteady boundary layer separation. Part 1: General theory
Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.
Vandommelen, Leon L.↗
On the Lagrangian description of unsteady boundary-layer separation. I - General theory
Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.
Van Dommelen, Leon L.↗
Thermal shaft effects on load-carrying capacity of a fully coupled, variable-properties cryogenic journal bearing
The purpose of this work was to perform a rather complete analysis for a cryogenic (oxygen) journal bearing. The Reynolds equation required coupling and simultaneous solution with the fluid energy equation. To correctly account for the changes in the fluid viscosity, the fluid energy equation was coupled with the shaft and bearing heat conduction energy equations. The effects of pressure and temperature on the density, viscosity, and load-carrying capacity were further discussed as analysis parameters, with respect to relative eccentricity and the angular velocity. The isothermal fluid case and the adiabatic fluid case represented the limiting boundaries. The discussion was further extrapolated to study the Sommerfeld number dependency on the fluid Nusselt number and its consequence on possible total loss of load-carrying capacity and/or seizure (catastrophic failure).
Braun, M. J.↗
Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++
Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Simulation of the impact of particle recycling on the plasma in MPS-LD device based on the BOUT++ LPD module
A linear plasma device (LPD) module has been developed under the BOUT++ framework to simulate plasma transport in the MPS-LD. However, previously, the LPD module used a simplistic neutral particle model that only includes particle density and velocity, which prevents the full understanding of the plasma-neutrals interactions. In this work, we further optimize the neutral model by using a more complete neutral fluid model containing the continuity equation, momentum equation, and energy equation. The reactions such as charge exchange, excitation, and radiation collisions are included. Since the neutral particle source is mainly provided by particle recycling from the target, a particle recycling model is employed, which includes both fast reflection and slow thermal release. The upgraded LPD module is applied to simulate the argon (Ar) discharge experiment of MPS-LD, and the benchmark against experiment measurement and SOLPS-ITER simulation results are presented. Good agreements are obtained, showing the validation of the upgraded module. After that, the impact of particle recycling on Ar plasma is investigated. It is found that a higher recycling coefficient (R) promotes the achievement of high-density plasma at the target. In conclusion, the recycled Ar atoms change target plasma pressure as well as plasma-neutral collisions, which both contribute to plasma momentum loss, thus promoting the rollover of ion flux to the target.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Application of a general boundary layer analysis to turbulent boundary layers subjected to strong favorable pressure gradients
Theoretical predictions of turbulent boundary layer development under the influence of strong favorable pressure gradients made using a finite-difference calculation procedure are compared to experimental data. Comparisons are presented for low speed flows with and without wall heat transfer as well as for supersonic flows with adiabatic walls. The turbulence model used is governed by an integral form of the turbulence kinetic energy equation and the results are compared with predictions made using a conventional equilibrium turbulence model based upon Prandtl's mixing length, a Clauser-type eddy viscosity model used by Cebecci and Mosinskis, and a two-equation turbulence energy model of Launder and Jones.
Kreskovsky, J. P.↗
Thermodynamic energy balance equations for Space Shuttle Orbiter gas compartment during ascent and re-entry
Thermodynamic energy balance equations are derived and applied to midsection Orbiter-payload atmospheric thermal math models (TMMs) to predict Orbiter component, element, compartment, internal insolation and structure temperatures in support of NASA/JSC mission planning, postflight thermal analysis and payload thermal integration planning. The equations are extended and applied to the forward section, midsection, and aft section of the TMMs for five Orbiter mission phases: prelaunch on pad with purge, lift-off to ascent, re-entry to touchdown, post landing without purge, and post-landing with purge. Predicted results from the 390 node/DFI atmospheric TMM are in good agreement with STS-1 flight measurement data.
Ting, P. C.↗
Nozzle Flow with Vibrational Nonequilibrium
Flow of nitrogen gas through a converging-diverging nozzle is simulated. The flow is modeled using the Navier-Stokes equations that have been modified for vibrational nonequilibrium. The energy equation is replaced by two equations. One equation accounts for energy effects due to the translational and rotational degrees of freedom, and the other accounts for the affects due to the vibrational degree of freedom. The energy equations are coupled by a relaxation time which measures the time required for the vibrational energy component to equilibrate with the translational and rotational energy components. An improved relaxation time is used in this thesis. The equations are solved numerically using the Steger-Warming flux vector splitting method and the Implicit MacCormack method. The results show that uniform flow is produced outside of the boundary layer. Nonequilibrium exists in both the converging and diverging nozzle sections. The boundary layer region is characterized by a marked increase in translational-rotational temperature. The vibrational temperature remains frozen downstream of the nozzle, except in the boundary layer.
Landry, John Gary↗
Analytical study of laser supported combustion waves in hydrogen
A one-dimensional energy equation, with constant pressure and area, was used to model the LSC wave. This equation balances convection, conduction, laser energy absorption, radiation energy loss and radiation energy transport. Solutions of this energy equation were obtained to give profiles of temperature and other properties, as well as the relation between laser intensity and mass flux through the wave. The flow through the LSC wave was then conducted through a variable pressure, variable area streamtube to accelerate it to high speed, with the propulsion application in mind. A numerical method for coupling the LSC wave model to the streamtube flow was developed, and a sample calculation was performed. The result shows that 42% of the laser power has been radiated away by the time the gas reaches the throat. It was concluded that in the radially confined flows of interest for propulsion applications, transverse velocities would be less important than in the unconfined flows where air experiments have been conducted.
Kemp, N. H.↗
The Numerical Analysis of a Turbulent Compressible Jet
A numerical method to simulate high Reynolds number jet flows was formulated and applied to gain a better understanding of the flow physics. Large-eddy simulation was chosen as the most promising approach to model the turbulent structures due to its compromise between accuracy and computational expense. The filtered Navier-Stokes equations were developed including a total energy form of the energy equation. Subgrid scale models for the momentum and energy equations were adapted from compressible forms of Smagorinsky's original model. The effect of using disparate temporal and spatial accuracy in a numerical scheme was discovered through one-dimensional model problems and a new uniformly fourth-order accurate numerical method was developed. Results from two- and three-dimensional validation exercises show that the code accurately reproduces both viscous and inviscid flows. Numerous axisymmetric jet simulations were performed to investigate the effect of grid resolution, numerical scheme, exit boundary conditions and subgrid scale modeling on the solution and the results were used to guide the three-dimensional calculations. Three-dimensional calculations of a Mach 1.4 jet showed that this LES simulation accurately captures the physics of the turbulent flow. The agreement with experimental data was relatively good and is much better than results in the current literature. Turbulent intensities indicate that the turbulent structures at this level of modeling are not isotropic and this information could lend itself to the development of improved subgrid scale models for LES and turbulence models for RANS simulations. A two point correlation technique was used to quantify the turbulent structures. Two point space correlations were used to obtain a measure of the integral length scale, which proved to be approximately 1/2 D(sub j). Two point space-time correlations were used to obtain the convection velocity for the turbulent structures. This velocity ranged from 0.57 to 0.71 U(sub j).
DeBonis, James R.↗
Cosmological measurements from void-galaxy and galaxy-galaxy clustering in the Sloan Digital Sky Survey
We present the cosmological implications of measurements of void-galaxy and galaxy-galaxy clustering from the Sloan Digital Sky Survey (SDSS) Main Galaxy Sample (MGS), Baryon Oscillation Spectroscopic Survey (BOSS), and extended BOSS (eBOSS) luminous red galaxy catalogues from SDSS Data Release 7, 12, and 16, covering the redshift range 0.07 < $z$ < 1.0. We fit a standard ΛCDM cosmological model as well as various extensions, including a constant dark energy equation of state not equal to −1, a time-varying dark energy equation of state, and these same models allowing for spatial curvature. Results on key parameters of these models are reported for void-galaxy and galaxy-galaxy clustering alone, both of these combined, and all these combined with measurements from the cosmic microwave background (CMB) and supernovae (SN). For the combination of void-galaxy and galaxy-galaxy clustering, we find tight constraints of Ω m = 0.356 ± 0.024 for a base ΛCDM cosmology, Ω m = $0.391^{+0.028}_{-0.021}$, w = $-1.50^{+0.43}_{-0.28}$ additionally allowing the dark energy equation of state $w$ to vary, and Ω m = $0.331^{+0.067}_{-0.094}$, $w=-1.41^{+0.70}_{-0.31}$, and Ω k = $0.06^{+0.18}_{-0.13}$ further extending to non-flat models. The combined SDSS results from void-galaxy and galaxy-galaxy clustering in combination with CMB+SN provide a 30 per cent improvement in parameter Ω m over CMB+SN for ΛCDM, a 5 per cent improvement in parameter Ω m when $w$ is allowed to vary, and a 32 per cent and 68 per cent improvement in parameters Ω m and Ω k when allowing for spatial curvature.
79 ASTRONOMY AND ASTROPHYSICS↗
Pressure updating methods for the steady-state fluid equations
We consider the steady state equations for a compressible fluid. Since we wish to solve for a range of speeds we must consider the equations in conservation form. For transonic speeds these equations are of mixed type. Hence, the usual approach is to add time derivatives to the steady state equations and then march these equations in time. One then adds a time derivative of the density to the continuity equation, a derivative of the momentum to the momentum equation and a derivative of the total energy to the energy equation. This choice is dictated by the time consistent equations. However, since we are only interested in the steady state this is not necessary. Thus we shall consider the possibility of adding a time derivative of the pressure to the continuity equation and similar modifications for the energy equation. This can then be generalized to adding combinations of time derivatives to each equation since these vanish in the steady state. When using acceleration techniques such as residual smoothing, multigrid, etc. these are applied to the pressure rather than the density. Hence, the code duplicates the behavior of the incompressible equations for low speeds.
Fiterman, A.↗
Multi-stage preconditioners for thermal–compositional–reactive flow in porous media
We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.
42 ENGINEERING↗
Turbulence models and Reynolds analogy for two-dimensional supersonic compression ramp flow
Results of the application of turbulence models and the Reynolds analogy to the Navier-Stokes computations of Mach 2.9 two-dimensional compression ramp flows are presented. The Baldwin-Lomax eddy viscosity model and the kappa-epsilon turbulence transport equations for the turbulent momentum flux modeling in the Navier-Stokes equations are studied. The Reynolds analogy for the turbulent heat flux modeling in the energy equation was also studied. The Navier-Stokes equations and the energy equation were numerically solved for the flow properties. The Reynolds shear stress, the skin friction factor, and the surface heat transfer rate were calculated and compared with their measurements. It was concluded that with a hybrid kappa-epsilon turbulence model for turbulence modeling, the present computations predicted the skin friction factors of the 8 deg and 16 deg compression ramp flows and with the turbulent Prandtl number Pr(sub t) = 0.93 and the ratio of the turbulent thermal and momentum transport coefficients mu(sub q)/mu(sub t) = 2/Prt, the present computations also predicted the surface heat transfer rates beneath the boundary layer flow of the 16 compression ramp.
Wang, Chi R.↗
Thermodynamic evaluation of transonic compressor rotors using the finite volume approach
A method was developed which calculates two-dimensional, transonic, viscous flow in ducts. The finite volume, time marching formulation is used to obtain steady flow solutions of the Reynolds-averaged form of the Navier Stokes equations. The entire calculation is performed in the physical domain. The method is currently limited to the calculation of attached flows. The features of the current method can be summarized as follows. Control volumes are chosen so that smoothing of flow properties, typically required for stability, is now needed. Different time steps are used in the different governing equations to improve the convergence speed of the viscous calculations. A new pressure interpolation scheme is introduced which improves the shock capturing ability of the method. A multi-volume method for pressure changes in the boundary layer allows calculations which use very long and thin control volumes. A special discretization technique is also used to stabilize these calculations. A special formulation of the energy equation is used to provide improved transient behavior of solutions which use the full energy equation. The method is then compared with a wide variety of test cases. The freestream Mach numbers range from 0.075 to 2.8 in the calculations. Transonic viscous flow in a converging diverging nozzle is calculated with the method; the Mach number upstream of the shock is approximately 1.25. The agreement between the calculated and measured shock strength and total pressure losses is good. Essentially incompressible turbulent boundary layer flow in a adverse pressure gradient is calculated and the computed distribution of mean velocity and shear stress are in good agreement with the measurements. At the other end of the Mach number range, a flat plate turbulent boundary layer with a freestream Mach number of 2.8 is calculated using the full energy equation; the computed total temperature distribution and recovery factor agree well with the measurements when a variable Prandtl number is used through the boundary layer.
Nicholson, S.↗