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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Magnetic-field generation and its effect on ablative Rayleigh–Taylor instability in diffusive ablation fronts

The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A subgrid model for electron-scale turbulent transport in global ion-scale gyrokinetic simulations of tokamak plasmas

A subgrid electron–temperature–gradient (ETG) model is demonstrated here which averages local electron-scale turbulence over intermediate scales in space and time to include in global ion–temperature–gradient (ITG) simulations. This approach results in ion-scale equations, which incorporate the electron heat transport from ETG turbulence and the effects of electron-scale turbulence on the ion scale. Flux-tube ETG Cyclone Base Case simulations are performed at different radial locations and a kinetic form of the flux is added to global ion-scale (IS) simulations as a source term. Analytic radial profiles of ETG heat flux are constructed and compared to flux-tube results at multiple radial locations. Different ratios of ITG to ETG heat flux levels are considered and the results of capturing ETG heat transport in global ion-scale simulations are discussed. The potential coupling of the ETG streamer potential and intermediate-scale zonal flows to the ion scale is further addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

MBMS1.0: An Open-Source Code for Modeling and Simulation of Membrane-Based Dehumidification and Energy Recovery

Membrane-based dehumidification is currently being considered as a promising solution for the building application due to its low cost and very limited energy consumption. Developing a simple and efficient open-source code simulation tool is important for boosting the optimization and evaluation of such device in HVAC community. This paper reports a first-order physics based model which accounts for the fundamental heat and mass transfer of humid-air vapor at feed side to flow stream at permeate side. The current model comprises two membrane mass transfer submodels (i.e. microstructure model and performance map model); and it adopts a segment-by-segment methodology for discretizing heat and mass transfer governing equations. The model is capable of simulating both dehumidifiers and energy recovery ventilators with parallel-flow cross-flow, and counter-flow configurations. The model was validated with the measurements at appropriate device. The practices in dehumidification and energy recovery exchangers are also discussed. The model and open-source codes are expected to become a solid fundament for developing a more comprehensive and accurate membrane-based dehumidification in the future.

Gao, Zhiming↗
DOE OSTI2020

FY20 Verification of BISON Using Analytic and Manufactured Solutions

In 2010, the U.S. Department of Energy created its first Energy Innovation Hub, which is focused on developing high-fidelity and high-resolution Modeling and Simulation (M&S) tools for modeling of Light Water Reactors (LWRs). This hub, Consortium for Advanced Simulation of LWRs (CASL), has developed an LWR simulation tool called Virtual Environment for Reactor Applications (VERA). The multi-physics capability of VERA is achieved through the coupling of single-physics codes, including BISON, CTF, MPACT, and MAMBA. BISON is a fuel performance code which models the thermo-mechanical behavior of nuclear fuel using high performance M&S. It is capable of modeling traditional LWR fuel rods, fuel plates, and TRi-structural ISOtropic (TRISO) fuel particles. It can employ three-dimensional Cartesian, two-dimensional axisymmetric cylindrical, or one-dimensional radial spherical geometry. It includes empirical models for a large variety of fuel physics: temperature- and burnup-dependent thermal properties, fuel swelling and densification, fission gas production, cladding creep, fracture, cladding plasticity, and gap/plenum models. This document details a series of code verification test problems that are used to test BISON. These problems add confidence that the BISON code is a faithful representation of its underlying mathematical model. The suite of verification tests are mapped to the underlying conservation equations solved by the code: heat conduction, mechanics, and species conservation. Twenty-two problems are added for the heat conduction solution, two for the mechanics solution, and none for species conservation. Method of Manufactured Solutions (MMS) capability is demonstrated with three problems, and temperature drops across the fuel gap are tested.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Predictive Models and Novel Accelerated Tests for the Reliability of Cell Metallization in Photovoltaic Modules (Final Report)

Studies of metallization corrosion m photovoltaics have mainly been limited to comparisons of modules placed in accelerated chambers to fielded modules [1]. Damp Heat accelerated tests and phenomenological equations [2] are used to assess metallization corrosion without understanding the effect of UV light and temperature and humidity cycles on encapsulant adhesion degradation. The roles of encapsulant in-and out-diffusions of moisture and encapsulant impurities are important. Furthermore, few photovoltaic metallization corrosion studies included the role of bias and leakage currents, which are crucial in the electrochemical reaction of metallization. Leakage currents can highly accelerate the corrosion mechanism and are important to include in the studies of corrosion. In our research plan, we will address the following gaps in the PV community's understanding of metallization corrosion: (1) metallization corrosion with bias, humidity, and impurities in the encapsulant or metallization; (2) humidity diffusion through fresh and degraded encapsulants; and (3) comparison of model predictions with outdoor field modules and SunPower's extensive data for its back-contact and front-contact fleets [2] along with NREL's store of >20 year old modules. Our goal is to build models and accelerated tests to predict long term degradation of metallization corrosion of photovoltaic modules in the field. Our studies will include metals used in c-Si solar cells (Cu, Ag, and Al) and commonly used encapsulants (EV A (ethylene vinyl acetate), TPO (thermoplastic olefin), and silicone).

14 SOLAR ENERGY↗
DOE OSTI2025

Theoretical Modeling of Reactor Relevant Conditions for Plasma Jet Driven Magneto-Inertial Fusion

The Charger Advanced Power and Propulsion Laboratory (CAPP), a laboratory within the Propulsion Research Center (PRC) at the University of Alabama in Huntsville (UAH) is working with Los Alamos National Laboratory (LANL). to develop models and inform on promising paths for high gain magneto-inertial fusion (MIF) conditions. This report provides a framework for identifying promising conditions for achieving ignition in plasma-jet-driven magneto-inertial fusion (PJMIF)[1]. As proposed, for the first part of the contract, UAH proposes to develop a gain over unity set of stagnation conditions to provide a state of plasma conditions to achieve to set long terms goals for the PJMIF program. Specifically, UAH will model PJMIF stagnation conditions to include radiation, heat transfer, two temperature energy equations, fusion reactivity and nonlocal fusion product deposition, but no hydrodynamics for these purposes. These calculations will use a stationary plasma model to reduce simulation complexity—focusing on a DT target at 10 keV. Subsequent work will include a DD plasma layer acting as an afterburner. UAH will assume an initial magnetic field without any consideration of the topology, just assume a field strength, most likely scaled with consideration of the local hall parameter. This effort will inform the team on the tradeoff between mass, peak target field, etc and the achievable gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A Comprehensive Review of Latent Space Dynamics Identification Algorithms for Intrusive and Non-Intrusive Reduced-Order-Modeling

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

Computational Engineering, Finance, and Science (c↗
DOE OSTI2020

Three-Phase Modeling in Sockeye

The heat pipe simulation tool Sockeye is given the capability to allow the presence of the solid phase of the heat pipe working fluid, which is necessary for modeling transients to and from room temperature, since microreactor designs feature high-temperature heat pipes employing metal working fluids. The liquid phase equations are converted to equations for the mixture of liquid and solid phases, with a simple latent heat relation determining the fractions of these phases. This capability is introduced while preserving equivalent solutions for the liquid-vapor regime. Theory is presented, along with numerical examples to illustrate the propagation of a melt front and startup from the solid state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Photochemically induced acousto-optics in gases

Acousto-optics consists of launching acoustic waves in a medium (usually a crystal) in order to modulate its refractive index and create a tunable optical grating. Here, in this article, we present the theoretical basis of an alternative scheme to generate acousto-optics in a gas, where the acoustic waves are initiated by the localized absorption (and thus gas heating) of spatially modulated UV light, as was demonstrated by Michine and Yoneda [Commun. Phys. 3, 24 (2020)]. We identify the chemical reactions initiated by the absorption of UV light via the photodissociation of ozone molecules present in the gas, and calculate the resulting temperature increase in the gas as a function of space and time. Solving the Euler fluid equations shows that the modulated, isochoric heating initiates a mixed acoustic-entropy wave in the gas, whose high-amplitude density (and thus refractive index) modulation can be used to manipulate a high-power laser. We calculate that diffraction efficiencies near 100% can be obtained using only a few millimeters of gas containing a few percent ozone fraction at room temperature, with UV fluences of less than 100 mJ/cm 2 —consistent with the experimental measurements. Our analysis suggests possible ways to optimize the diffraction efficiency by changing the buffer gas composition. Gases have optics damage thresholds 2–3 orders of magnitude beyond those of solids; these optical elements should therefore be able to manipulate kilojoule-class lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Physics-Informed Neural Network (PINN) Prediction of Mixed Mass-Heat-Crystallization Limited Methane Hydrate Formation and Dissociation in Micro-Confinement

The creation and use of Physics-Informed Neural Networks (PINNs) for simulating the dynamics of methane hydrate formation and dissociation will be presented. The PINN framework's main benefit is its capacity to impose physical consistency with only a partial comprehension of the governing equations. This makes the algorithm especially useful for systems with little experimental evidence or a lack of theoretical knowledge. A strong basis for forecasting methane hydrate behavior over the verified operating ranges of 30.0-80.9 bar pressure and 1.0-4.0 K sub-cooling conditions is provided by the combination of conductive heat transfer equations and mixed mass-transfer–crystallization kinetics. PINNs were more accurate at predicting the mixed mass-heat-crystallization limited kinetics than conventional Artificial Neural Networks (ANNs), demonstrating remarkable predictive accuracy for methane hydrate production over the ANN model. The efficiency of incorporating physical limitations from first principles into machine learning frameworks for methane hydrate crystallizations is reinforced by these findings. For hydrate-related applications in energy generation, carbon sequestration, and climate modelling, our study establishes PINNs as a computational tool that is both scalable and efficient. The proven capacity to close the gap between conventional physics-based simulations and solely data-driven models creates new opportunities for expedited hydrate research and practical applications.

Hartman, Ryan L [NYU Tandon School of Engineering]↗
DOE OSTI2021

Evaluation of pressure reconstruction techniques for Model Order Reduction in incompressible convective heat transfer

This paper compares pressure reconstruction strategies in Model Order Reduction for incompressible flows with convective heat transfer. The Navier-Stokes equation are reduced along with the passive scalar transport equation for the temperature using the POD-Galerkin technique. Six different pressure reconstruction methods are evaluated, two of which are novel to the best of the authors’ knowledge. Accurate pressure reconstruction is key to avoid error buildup in when solving for the conservation of linear momentum at the reduced level. The six approaches are compared using Direct Numerical Simulations of convective heat exchange processes in a 3D Backward Facing Step with a heated cylinder. Additionally, when comparing time-averaged metrics, we observe that the reconstruction methods that approximate the reduced pressure field using techniques borrowed from full order models (mechanical analogy, pressure Poisson, and velocity supremizers) yield higher errors than the methods that seek to stabilize the reduced systems (reduced residual stabilization, artificial divergence, and Uzawa operator).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Self-powered Heating: Efficiency Analysis

Conventional fuel-fired heating devices such as furnaces, boilers, and water heaters have fuel efficiency less than 100% on the basis of higher heating value. They also require electricity from the electric grid to power parasitic loads such as blowers, pumps, fans, and ignitors. The primary energy efficiency of the device accounts for both fuel used on-site and primary energy used off-site to produce electric power used by the device. This work compares conventional fuel-fired heating devices to two types of self-powered devices. A self-powered device (SPD) integrates a power cycle onboard to eliminate consumption of grid electricity. We assume that all heat rejected by the onboard power cycle is added to the process fluid, so that, compared with a conventional device, the same amount of heat is provided to the process fluid and the same amount of fuel is consumed, but grid electricity consumption is eliminated. The first SPD type is the basic one: exactly the electricity required is generated. The second type considered is the SPD with heat pump (SPD-HP), in which the power cycle generates more electricity than needed for parasitic loads, and the excess electricity is used to power a heat pump. The heat pump extracts additional heat from the ambient to boost efficiency. Both SPD and SPD-HP self-consume all the generated electricity, in contrast to combined heat and power (CHP) systems that export electricity. In this work, equations are derived to express the efficiency of three classes of heating devices: conventional (consuming grid electricity), self-powered (consuming no grid electricity), and self-powered with heat pump. The efficiency of each is derived as a function of up to six factors: (1) the fraction of combustion heat captured, (2) the rate of parasitic power consumption, (3) the fraction of electric energy dissipated as useful heat, (4) the power cycle conversion efficiency, (5) the grid efficiency, when applicable, and (6) the heat pump COP, when applicable. Scenarios are identified in which it is possible to achieve efficiency greater than 100% on a higher heating value basis. Plausible configurations using existing technology options are outlined.

Gluesenkamp, Kyle↗
DOE OSTI2023

Motorette Thermal Performance Testing and Modeling for an Electric Motor with Additively Manufactured Hollow Conductors with Integrated Heat Pipes

This paper discusses the design, build and test of a motorette to characterize the thermal performance of an additively manufactured coil integrated with heat pipes. The motorette is powered using a variable frequency AC power supply and cooled through two independent cooling systems. The first cooling system cools the stator core using forced convection of air flowing over a finned heat sink mounted on the stator outer diameter. The second cooling system cools the heat pipes using forced convection of 50-50 water ethylene glycol (WEG) mixture flowing through condenser chambers. A 3D thermal FEA model of the setup is built and heat transfer coefficients (HTC) of convective boundaries are computed using coolant flow rates from test data and empirical equations. Temperature at locations on the coil and heat sink are recorded and compared to thermal model predicted values. A maximum temperature error of 22.7% occurs at 180ARMS and 800 Hz operating point for the AM dual coil. The thermal model provides higher temperatures compared to test results and hence it is on the conservative side but in general, there is good correlation between test and model results. Lastly, opportunities for improvement to both test data measurement and 3D thermal FEA modeling are discussed.

additive manufacturing↗
DOE OSTI2022

Boiling and cavitation caused by transient heat transfer in superfluid helium-4

Superfluid helium-4 (He II) has been widely utilized as a coolant in various scientific and engineering applications due to its superior heat transfer capability. An important parameter required in the design of many He II based cooling systems is the peak heat flux $q^*$, which refers to the threshold heat flux above which boiling spontaneously occurs in He II. Past experimental and numerical studies showed that $q^*$ increases when the heating time $t_h$ is reduced, which leads to an intuitive expectation that very high $q^*$ may be achievable at sufficiently small $t_h$. Knowledge on how $q^*$ actually behaves at small $t_h$ is important for applications such as laser ablation in He II. Furthermore we present a numerical study on the evolution of the thermodynamic state of the He II in front of a planar heater by solving the He II two-fluid equations of motion. For an applied heat flux, we determine the heating time beyond which the He II near the heater transits to the vapor phase. As such, a curve correlating $q^*$ and $t_h$ can be obtained, which nicely reproduces some relevant experimental data. Surprisingly, we find that there exists a critical peak heat flux $q^*_c$, above which boiling occurs nearly instantaneously regardless of $t_h$. We reveal that the boiling in this regime is essentially cavitation caused by the combined effects of the first-sound and the second-sound waves in He II. Based on this physical picture, an analytical model for $q^*_c$ is developed, which reproduces the simulated $q^*_c$ values at various He II bath temperatures and hydrostatic head pressures. This work represents a major progress in our understanding of transient heat transfer in He II.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermophysical Properties Experimentally Tested for NaCl-KCl-MgCl 2 Eutectic Molten Salt as a Next-Generation High-Temperature Heat Transfer Fluids in Concentrated Solar Power Systems

A new eutectic chloride molten salt, MgCl 2 -KCl-NaCl (wt% 45.98–38.91–15.11), has been recognized as one of the most promising high-temperature heat transfer fluids (HTF) for both heat transfer and thermal storage for the third-generation concentrated solar power (CSP) systems. For the first time, some essential thermophysical properties of this eutectic chloride molten salt needed for basic heat transfer and energy storage analysis in the application of concentrating solar power systems have been experimentally tested and provided as functions of temperature in the range from 450 °C to 700 °C. The studied properties include heat capacity, melting point, heat of fusion, viscosity, vapor pressure, density, and thermal conductivity. The property equations provide essential database for engineers to use to calculate convective heat transfer in concentrated solar receivers, heat exchangers, and thermal storage for concentrated solar power plants.

14 SOLAR ENERGY↗
DOE OSTI2023

The Effect of Tropical Pacific Air‐Sea Coupling on the Rainfall Response to Quadrupled CO 2 Forcing

We perform quadrupled CO 2 climate simulations with the Community Earth System Model version 1 (CESM1) to study how air-sea coupling affects the response of tropical rainfall under global warming. We use a hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows directed perpendicular to the wind, and the near-equator frictional flows directed in the same direction as the wind. We show that the Pacific Ocean's enhanced equatorial warming pattern (EEW) and equatorward ITCZ contraction observed in previous climate simulations emerge when the ocean model includes wind-driven Ekman and frictional flows. Furthermore, the near-equator frictional flow contributes more than half of the heat convergence in the equatorial Pacific Ocean. Finally, we show that although Ekman flow and near-equator frictional flow can both result in EEW, their coupled interactions with the Hadley circulation lead to opposite feedbacks on EEW's strength.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗
DOE OSTI2022

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗