Search NASA⌕ Search

SEARCH · Search NASA

Results for “heat equation”

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
DOE OSTI2022

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

AEROICESTUDY-Colorado State University Ice Spectrometer

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0 °C down to -26 to -29 °C) of atmospheric ice nucleating particle (INP) concentrations taken at ARM's Southern Great Plains (SGP) Site Guest Instrument Facility (GIF), in Lamont, Oklahoma, during the Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY). A filter sampler was mounted approximately 6 m above ground level atop the GIF rooftop platform. Single-use filter units open to the atmosphere fitted with pre-cleaned and pre-sterilized, 47-mm diameter Nuclepore polycarbonate filters were used. Filters were typically drawn for either a full day period (9AM-5PM CDT) or broken up into morning (9AM-12PM CDT) and afternoon (1PM-5PM CDT) periods. Additionally, 4 overnight filters were taken (6PM-8AM CDT). The volume of air filtered was typically between 3000 and 8000 standard liters for daily samples. After collection, filters were stored at -20 °C until processed. Initial processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurement of INPs, collected aerosol particles were re-suspended in 8 mL of 0.1 µm-filtered deionized (DI) water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples were cooled at 0.33 °C per min and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Two-tailed, 95% confidence intervals for binomial sampling were based on Equation 2 from Agresti and Coull (1998). Aliquots of suspensions from selected samples were also heat treated (95 °C for 20 min) to denature and deactiv

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Application of Recursive Estimation to Heat Tracing for Groundwater/Surface‐Water Exchange

Abstract We present and demonstrate a recursive‐estimation framework to infer groundwater/surface‐water exchange based on temperature time series collected at different vertical depths below the sediment/water interface. We formulate the heat‐transport problem as a state‐space model (SSM), in which the spatial derivatives in the convection/conduction equation are approximated using finite differences. The SSM is calibrated to estimate time‐varying specific discharge using the Extended Kalman Filter (EKF) and Extended Rauch‐Tung‐Striebel Smoother (ERTSS). Whereas the EKF is suited to real‐time (“online”) applications and uses only the past and current measurements for estimation (filtering), the ERTSS is intended for near‐real time or batch‐processing (“offline”) applications and uses a window of data for batch estimation (smoothing). The two algorithms are demonstrated with synthetic and field‐experimental data and are shown to be efficient and rapid for the estimation of time‐varying flux over seasonal periods; further, the recursive approaches are effective in the presence of rapidly changing flux and (or) nonperiodic thermal boundary conditions, both of which are problematic for existing approaches to heat tracing of time‐varying groundwater/surface‐water exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗
DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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↗
DOE OSTI2026

Effect of soot concentration on radiative transport in a 2 m JP-8 pool fire

Hubbard et al. (2024) investigated conjugate heat transfer in a 7.9 m JP-8 pool fire, utilizing both heavily sooting and non-sooting fuel models. Their findings indicated that all models produced pool radiative heat fluxes within the measurement uncertainty for large-scale outdoor experiments. This study aims to further investigate the influence of soot concentration on radiative transport. To this end, Large Eddy Simulations (LES) of a 2 m diameter JP-8 pool fire were conducted, employing four turbulent combustion chemistry models: the Eddy Dissipation Concept (EDC) and three variants of flamelet model. Additionally, three variations of the one-equation subgrid scale turbulence kinetic energy model, and multiple mesh resolutions were utilized. The primary validation metrics were pool radiative heat flux and plume radiative fraction, with simulation data compared to experimental results published by Blanchat et al. (2010). The EDC model produced high soot concentrations near the pool surface, resulting in a significant shielding effect. In contrast, the flamelet combustion model has different kinetic limiters for soot production, leading to soot formation higher in the flame, where it can be more readily oxidized. Flamelet models incorporating both acetylene and benzene contributions to soot growth demonstrated better agreement with experimental data. To assess the sensitivity of pool radiative heat flux to soot mass concentration, we used a one-dimensional approximation over a wide range of soot concentrations and compared it to three-dimensional solutions obtained using Discrete Ordinates. Furthermore, this analysis illustrates the impacts of soot shielding and Turbulence Radiation Interactions (TRI), with TRI emission enhancement factors predicted to be in the range of 3–5 near the edges of the plume.

Computational fluid dynamics↗
DOE OSTI2022

Characterization of the high-pressure and high-temperature phase diagram and equation of state of chromium

The high-pressure and high-temperature phase diagram of chromium has been investigated both experimentally (in situ), using a laser-heated diamond-anvil cell technique coupled with synchrotron powder X-ray diffraction, and theoretically, using ab initio density-functional theory simulations. In the pressure–temperature range covered experimentally (up to 90 GPa and 4500 K, respectively) only the solid body-centred-cubic and liquid phases of chromium have been observed. Experiments and computer calculations give melting curves in agreement with each other that can both be described by the Simon–Glatzel equation $T_{m}(P) = 2136K (1 + P/25.9)^{0.41}$ In addition, a quasi-hydrostatic equation of state at ambient temperature has been experimentally characterized up to 131 GPa and compared with the present simulations. Both methods give very similar third-order Birch–Murnaghan equations of state with bulk moduli of 182–185 GPa and respective pressure derivatives of 4.74–5.15. According to the present calculations, the obtained melting curve and equation of state are valid up to at least 815 GPa, at which pressure the melting temperature is 9310 K. Finally, from the obtained results, it was possible to determine a thermal equation of state of chromium valid up to 65 GPa and 2100 K.

36 MATERIALS SCIENCE↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2023

Temperature separation under compression of moderately coupled plasma

In moderately coupled plasmas, a significant fraction of the internal energy resides in electric fields. As these plasmas are heated or compressed, the shifting partition of energy between particles and fields leads to surprising effects, particularly when ions and electrons have different temperatures. In this work, quasi-equations of state (quasi-EOS) are derived for two-temperature moderately coupled plasma in a thermodynamic framework and expressed in a simple form. These quasi-EOS readily yield expressions for correlation heating, in which heating of the electrons causes a rapid increase in ion temperature even in the absence of collisional energy exchange between species. It is also shown that, remarkably, compression of moderately coupled plasma drives a temperature difference between electrons and ions, even when the species start at equal temperatures. These additional channels for ion heating may be relevant in designing ignition schemes for inertial confinement fusion.

Fetsch, H. (ORCID:0000000256133975)↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure

The objective of this paper is to develop an analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure. This expression for the temperature is obtained for both constant and variable specific heat models and can be used in common models for the material yield strength that include thermal softening due to increase in temperature. Additionally, the material parameters in the model can easily be calibrated to match known experimental data and the model has been used to determine the value of the shock pressure at the onset of melting during shock loading.

42 ENGINEERING↗
DOE OSTI2020

Equilibrium Thermodynamic Properties of Aqueous Solutions of Ionic Liquid 1-Ethyl-3-Methylimidazolium Methanesulfonate [EMIM][MeSO3]

Abstract The ionic liquid 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO 3 ]) has been considered as a promising alternative desiccant to triethylene glycol and lithium bromide commonly used in the industry. In this paper, the water activity coefficient of this binary system was measured from 303 K to 363 K with water concentration from 18% to 92%. The interaction energies between the ionic liquid molecules ( $${g}_{22}$$ g 22 ) and between the ionic liquid and water molecules ( $${g}_{12}$$ g 12 ) for the [EMIM][MeSO 3 ]/water binary system were determined from the water activity coefficient data using the Non-Random Two-Liquid (NRTL) model. The magnitude of the interaction energy between the [EMIM][MeSO 3 ] and water molecules ( $${g}_{12}$$ g 12 ) was found to be in the range of 45~49 kJ/mol, which was about 20% larger than that between the water molecules ( $${g}_{11}$$ g 11 ) in the [EMIM][MeSO 3 ]/water system. The large ( $${g}_{12}$$ g 12 ) can explain many observed macroscopic thermodynamic properties such as strong hygroscopicity in the ionic liquid [EMIM][MeSO 3 ]. These interaction energies were used to determine the heat of desorption of the [EMIM][MeSO 3 ]/water system, and the obtained heat of desorption was in good agreement with that calculated from the conventional Clausius-Clapeyron Equation.

42 ENGINEERING↗
DOE OSTI2021

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Post-detonation fireball modeling: Validation of freeze out approximations

A numerical simulation study was performed to examine the post-detonation reaction processes produced by the detonation of a 12 mm diameter hemispherical pentaerythritol tetranitrate (PETN) explosive charge. The simulations used a finite rate detailed chemical reaction model consisting of 59 species and 368 reactions to capture post-detonation reaction processes including air dissociation from Mach 19+ shock waves that initially break out of the PETN charge, reactions within the detonation products during expansion, and afterburning when the detonation products mix with the shock heated air. The multi-species and thermodynamically complete Becker-Kistiakowsky-Wilson real-gas equation of state is used for the gaseous phase to allow for the mixing of reactive species. A recent simplified reactive burn model is used to propagate the detonation through the charge and allow for detailed post-detonation reaction processes. The computed blast, shock structures, and mole fractions of species within the detonation products agree well with experimental measurements. A comparison of the simulation results to equilibrium calculations indicates that the assumption of a local equilibrium is fairly accurate until the detonation products rapidly cool to temperatures in the range of 1500-1900 K by expansion waves. Below this range, the computed results show mole fractions that are nearly chemically frozen within the detonation products for a significant portion of expansion. In conclusion, these results are consistent with the freeze out approximation used in the blast modeling community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Verification of a specialized hydrodynamic simulation code for modeling deflagration and detonation of high explosives

A specialized hydrodynamic simulation code has been developed and verified for the simulation of one-dimensional unsteady problems involving the detonation and deflagration of high explosives. To model all the relevant physical processes in these problems, a code is required to simulate compressible hydrodynamics, unsteady thermal conduction, and chemical reactions with complex rate laws. Several verification exercises are presented which test the implementation of these capabilities. The code also requires models for physics processes such as equations of state and conductivity for pure materials and mixtures as well as rate laws for chemical reactions. Additional verification tests are required to ensure that these models are implemented correctly. Though this code is limited in the types of problems it can simulate, its computationally efficient formulation allows it to be used in calibration studies for reactive burn models for high explosives. Furthermore, this study demonstrates how a series of verification tests can be used to ensure that the various physics processes needed to simulate complex phenomenon can be tested to ensure that they are correctly implemented.

97 MATHEMATICS AND COMPUTING↗