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

The Application of a Genetic Algorithm to the Optimization of a Mesoscale Model for Emergency Response

Besides solving the equations of momentum, heat, and moisture transport on the model grid, mesoscale weather models must account for subgrid-scale processes that affect the resolved model variables. These are simulated with model parameterizations, which often rely on values preset by the user. Such “free” model parameters, along with others set to initialize the model, are often poorly constrained, requiring that a user select each from a range of plausible values. Finding the values to optimize any forecasting tool can be accomplished with a search algorithm, and one such process—the genetic algorithm (GA)—has become especially popular. As applied to modeling, GAs represent a Darwinian process: an ensemble of simulations is run with a different set of parameter values for each member, and the members subsequently judged to be most accurate are selected as “parents” who pass their parameters onto a new generation. At the U.S. Department of Energy’s Savannah River Site in South Carolina, we are applying a GA to the Regional Atmospheric Modeling System (RAMS) mesoscale weather model, which supplies input to a model to simulate the dispersion of an airborne contaminant as part of the site’s emergency response preparations. An ensemble of forecasts is run each day, weather data are used to “score” the individual members of the ensemble, and the parameters from the best members are used for the next day’s forecasts. As meteorological conditions change, the parameters change as well, maintaining a model configuration that is best adapted to atmospheric conditions. Significance Statement We wanted to develop a forecasting system by which a weather model is run over the Savannah River Site each day and repeatedly adjusted according to how well it performed the previous day. To run the model, a series of values (parameters) must be set to control how the model will calculate winds, temperatures, and other desired variables. Each day the model was run several times using different combinations of these parameters and later compared with observed meteorological conditions. Parameters that produced the most accurate forecasts were preferentially reused to create the forecasts for the next day. The process was tested for the summer of 2020 and exhibited lower errors than forecasts produced by the model using default values of the parameters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Thermal transport exceeding bulk heat conduction due to nonthermal micro/nanoscale phonon populations

While classical size effects usually lead to a reduced effective thermal conductivity, we report here that nonthermal phonon populations produced by a micro/nanoscale heat source can lead to enhanced heat conduction, exceeding the prediction from Fourier's law. We study nondiffusive thermal transport by phonons at small distances within the framework of the Boltzmann transport equation (BTE) and demonstrate that the transport is significantly affected by the distribution of phonons emitted by the source. We discuss analytical solutions of the steady-state BTE for a source with a sinusoidal spatial profile, as well as for a three-dimensional Gaussian “hot spot,” and provide numerical results for single crystal silicon at room temperature. If a micro/nanoscale heat source produces a thermal phonon distribution, it gets hotter than that predicted by the heat diffusion equation; however, if the source predominantly produces low-frequency acoustic phonons with long mean free paths, it may get significantly cooler than that predicted by the heat equation, yielding an enhanced heat transport beyond bulk heat conduction.

42 ENGINEERING↗
DOE OSTI2022

Sheath formation around a dielectric droplet in a He atmospheric pressure plasma

Interactions at the interface between atmospheric pressure plasmas and liquids are being investigated to address applications ranging from nanoparticle synthesis to decontamination and fertilizer production. Furthermore, many of these applications involve activation of droplets wherein the droplet is fully immersed in the plasma and synergistically interacts with the plasma. To better understand these interactions, two-dimensional modeling of radio frequency (RF) glow discharges at atmospheric pressure operated in He with an embedded lossy dielectric droplet (tens of microns in size) was performed. The properties of the sheath that forms around the droplet were investigated over the RF cycle. The electric field in the bulk plasma polarizes the dielectric droplet while the electron drift in the external electric field is shadowed by the droplet. The interaction between the bulk and sheath electric fields produces a maximum in E/N (electric field/gas number density) at the equator on one side of the droplet where the bulk and sheath fields are aligned in the same direction and a minimum along the opposite equator. Due to resistive heating, the electron temperature T e is maximum 45° above and below the equator of the droplet where power deposition per electron is the highest. Although the droplet is, on the average, negatively charged, the charge density on the droplet is positive on the poles and negative on the equator, as the electron motion is primarily due to diffusion at the poles but due to drift at the equator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗
DOE OSTI2020

Conditional Generative Adversarial Networks for Solving Heat Transfer Problems

Generative Adversarial Networks (GANs) have been used as a deep learning approach to solving physics and engineering problems. Using deep learning for these problems is attractive in that reasonably accurate models can be inferred from only raw data, eliminating the need to define the exact physical equations governing a problem. We expand on previous work using GANs to generate steady-state solutions to the two-dimensional heat equation. Using a basic conditional GAN (cGAN), we generate accurate solutions for rectangular domains conditioned on four edge boundary conditions (MAE < 0.5%). For finding steady-state solutions over arbitrary two-dimensional domains (not constrained to rectangles), we use a cGAN designed for image-to-image translation. We train this GAN on various types of geometric domains (circles, squares, triangles, shapes with one circular or rectangular hole), achieving accurate results on test data made up of geometries similar to those in training (MAE < 1%). For both of these GANs, we experiment with different loss function terms, showing that a term using the gradients of solution images significantly improves the basic cGAN but not the image-to-image GAN. Lastly, we show that the image-to-image GAN performs poorly when applied to two-dimensional geometries that vary in structure from training data (MAE < 8% for shapes with multiple holes or different shaped holes). This demonstrates the cGAN's lack of generalizability. While the cGAN is an accurate and computationally efficient method when trained and tested on similarly structured data, it is a much less reliable method when applied to data that is slightly different in structure from the training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Reactor Scale-Up Methodology from Lab-Scale to Pilot-Scale Operations: Numerical Modeling of THFA Dehydration to DHP in Packed-Bed Reactors

This manuscript discusses developing a model-based scale-up methodology for a successful technology transfer of gas-phase catalytic reactors from a lab-scale to a pilot-scale operation. The manuscript demonstrates the methodology for gas-phase dehydration of tetrahydrofurfuryl alcohol (THFA) to dihydropyran (DHP) process over commercial Al 2 O 3 catalysts. A two-dimensional reactor model was developed using COMSOL Multiphysics 6.1 software. The model solves heat and mass transport equations in bed-scale and particle scales simultaneously. This powerful feature enables accurate prediction of the heat and mass transfer limitations in pilot-scale reactors, if any exists. Further, the model uses isothermal lab-scale experimental data to derive and validate the reaction chemistry, flow fields and boundary conditions. The model was then scaled-up to project conversion, selectivity, yield and formation rate of DHP in a pilot-scale reactor. The results highlight the complex nature of chemistry, heat, and mass transfer effects in lab-scale and pilot-scale reactors. The model results inform the possible operational limitations of the pilot-scale reactor and design strategies to improve process efficiency. Although the scale-up approach is explained through the THFA dehydration process, the methodology is applicable to any catalytic packed-bed reactor models for a successful process scale-up.

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

The interpretation of temperature and salinity variables in numerical ocean model output and the calculation of heat fluxes and heat content

Abstract. The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρcp0Θ, where cp0 is a constant “isobaric heat capacity”. However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model. We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A two-phase three-field modeling framework for heat pipe application in nuclear reactors

Heat pipes and two-phase thermosyphons are highly efficient heat transfer devices utilizing continuous evaporation and condensation of working fluid for two-phase heat transport in closed systems. Because of the nearly isothermal and fully passive phase-change heat transfer mechanism, heat pipes and thermosyphons have found many applications in nuclear engineering, space technologies, and other energy systems. High-temperature heat pipes are used in nuclear microreactors to remove fission power from the primary system and are coupled with power conversion systems or process heat applications. Modeling of the two-phase flow phenomena inside a heat pipe is essential to its design and safety analysis. In this study, a comprehensive one-dimensional two-phase three-field flow model has been developed for the analysis of heat pipes in normal operation conditions and transients. The conservation or field equations of mass, momentum, and energy were developed for the liquid film, vapor, and droplet. In addition, constitutive models or correlations were reviewed thoroughly and provided for the closure of the three-field equations. Specific constitutive equations regarding interfacial mass and heat transfer at two interfaces, namely film-gas interface and gas-droplet interface, were reviewed for droplet entrainment and deposition rates as well as film and droplet evaporation rates. Furthermore, mechanistic correlations of annular flow film thickness were recommended for the modeling of the thermosyphons without a wick as a critical constitutive correlation. Furthermore, experimental data needs from new experiments using a prototype working fluid or surrogate fluids for the model validation of high-temperature heat pipes in microreactors were recommended for future research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Depletion Chain Simplification With Pseudo-Nuclides to Model Decay Effects

This work introduces the novel usage of pseudo-nuclides to model decay effects which are otherwise absent from simplified depletion libraries. Pseudo-nuclides are artificial nuclides which, when added to a simplified depletion library, can preserve quantities of interest such as decay energy release or decay photon activity which are otherwise significantly under-predicted by simplified libraries. Several dozen pseudo-nuclides are capable of accurately modeling decay effects of hundreds of short-lived radionuclides. The computational resources needed to model the decay effects in depletion systems with simplified depletion libraries with pseudo-nuclides are significantly less than those needed when using simplified depletion libraries without pseudo-nuclides. When compared to the similar method of decay heat precursors for a given system, pseudo-nuclides allow for smaller depletion libraries to be generated because of their ability to more accurately model the fission product irradiation effect on decay heat.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A kinetic model for pure electron plasma compression via the rotating wall technique

An approximate model for pure electron plasma compression is developed for the case where the rotating wall (RW) electric field couples to the ${E}\times{B}$rotation and axial bounce motion of the electrons. The key assumption in the model is that, throughout the compression, the plasma remains in a slowly evolving thermal equilibrium defined by the plasma temperature and angular momentum. Linearized drift kinetic theory is employed to derive an expression for torque exerted by the RW field on the plasma through coupling to the resonant plasma particles, and averaging is used to find the torque that both compresses and heats the plasma. The evolution equations for the angular velocity and temperature of the plasma include the compression and heating from the torque and cooling from cyclotron radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Anisotropic Electron Fluid Closure Validated by in Situ Spacecraft Observations in the far Exhaust of Guide-field Reconnection

Anisotropic electron heating predictions of the equations of state of Le et al. (2009); https://doi.org/10.1103/PhysRevLett.102.085001 (“the EoS”) are applied to reconnection in the Earth's magnetosheath. The model is applicable to open systems where electrons are streaming along magnetic field lines into the reconnection region, sourced by fixed external reservoirs of plasma ambient to the reconnection region. While spacecraft observations have previously shown the EoS to hold in the region near the X-line, we find that for an event observed far downstream (~100 d i ) from the X-line both inflows and the exhaust follow the predictions of the EoS. Furthermore, the model underlying the EoS is extended to include additional perpendicular heating terms relevant to the considered event undergoing active magnetic compression with local electron trapping. Finally, the agreement between the spacecraft observations and the EoS emphasizes the roles of the ambient plasma sources and the dynamics of trapped and passing electrons in setting and controlling the anisotropic electron heating at large scale in naturally occurring plasma configurations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Can Biomass Burning Aerosol Induced Surface Cooling Be Amplified Through Sea Surface Temperature‐Cloud Feedback Over the Southeast Atlantic?

Biomass burning (BB) aerosols exert a strong surface cooling effect over the southeast Atlantic (SEA) via aerosol-radiation and aerosol-cloud interactions. The reduction of the sea surface temperature (SST) can trigger the SST-low cloud feedback. Whether this feedback can amplify the surface cooling effect is examined. The modeling results from the Community Earth System Model version 2 (CESM2) demonstrate that counterintuitively the cloud radiative effect (CRE) caused by the BB aerosols is weaker if SST-low cloud feedback is considered compared to fixed-SST simulation (-2.99 W m -2 vs. -4.79 W m -2 ). This is caused by (a) stronger sea breeze due to larger sea-land temperature contrast causing less smoke transport over SEA and (b) less moisture supply from surface due to colder SST. Changes in SST also lead to counterclockwise rotation of ocean circulation anomalies. Consequently, the excess heat transport from the equator reverses the direction of SST-cloud feedback in this region.

0305↗
DOE OSTI2021

Ab initio anharmonic thermodynamic properties of cubic Ca Si O 3 perovskite

We present an ab initio study of the thermodynamic properties of cubic CaSiO 3 perovskite (CaPv) over the pressure and temperature range of the Earth's lower mantle. We compute the anharmonic phonon dispersions throughout the Brillouin zone by utilizing the phonon quasiparticle approach, which characterizes the intrinsic temperature dependence of phonon frequencies and, in principle, captures full anharmonicity. Such temperature-dependent phonon dispersions are used to calculate ab initio free energy in the thermodynamic limit (N → ∞) within the framework of the phonon gas model. Accurate free energy calculations enable us to investigate cubic CaPv's thermodynamic properties, e.g., thermal expansivity, Grüneisen parameter, bulk modulus, heat capacity and thermal equation of state, where anharmonic effects are demonstrated. Furthermore, the present methodology provides an important theoretical approach for exploring phase boundaries, thermodynamic, and thermoelastic properties of strongly anharmonic materials at high pressures and temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗