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 253 records · Page 14

Equation of State of TiN at High Pressures and Temperatures: A Possible Host for Nitrogen in Planetary Mantles

Nitrogen, the most abundant element in Earth's atmosphere, is also a primary component of solid nitride minerals found in meteorites and on Earth's surface. If they remain stable to high pressures and temperatures, these nitrides may also be important reservoirs of nitrogen in planetary interiors. We used synchrotron X-ray diffraction to measure the thermal equation of state and phase stability of titanium nitride (TiN) in a laser-heated diamond anvil cell at pressures up to ~70 GPa and temperatures up to ~2500 K. TiN maintains the cubic B1 (NaCl-type) crystal structure over the entire pressure and temperature range explored. It has K 0 = 274(4) GPa, K 0 ' = 3.9(2), and γ 0 = 1.39(4) for a fixed V 0 = 76.516(30) Å 3 (based on experimental measurements), q = 1, and θ 0 = 579 K. Additionally, we collected Raman spectra of TiN up to ~60 GPa, where we find that the transverse acoustic (TA), longitudinal acoustic (LA), and transverse optical (TO) phonon modes exhibit mode Grüneisen parameters of 1.66(17), 0.54(15), and 0.93(4), respectively. Based on our equation of state, TiN has a density of ~5.6–6.4 g/cm 3 at Earth's lower mantle conditions, significantly more dense than both the mantle of the Earth and the estimated densities of the mantles of other terrestrial planets, but less dense than planetary cores. Furthermore, we find that TiN remains stable against physical decomposition at the pressures and temperatures found within Earth's mantle, making it a plausible reservoir for deep planetary nitrogen if chemical conditions allow its formation.

58 GEOSCIENCES↗

Evaluating Diesel/Biofuel Blends Using Artificial Neural Networks and Linear/Nonlinear Equations

Abstract The use of biomass-derived additives in diesel fuel mixtures has the potential to increase the fuel’s efficiency, decrease the formation of particulate matter during its combustion, and retain the fuel’s behavior in cold weather. To this end, identifying compounds that enable these behaviors is paramount. The present work utilizes a series of linear and non-linear equations in series with artificial neural networks to predict the cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value of multi-component blends. Property values of pure components are predicted using artificial neural networks trained with existing experimental data, and these predictions and their expected errors are propagated through linear and non-linear equations to obtain property predictions for multi-component blends. Individual component property prediction errors, defined by blind prediction median absolute error, are 4.91 units, 7.84 units, 0.06 cSt, 4.00 °C, and 0.55 MJ/kg for cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value respectively. On average, property predictions for blends are shown to be accurate to within 6% of the blends’ experimental values. Further, a multitude of compounds expected to be produced from catalytically upgrading products of fast pyrolysis are evaluated with respect to their behavior in diesel fuel blends.

09 BIOMASS FUELS↗

Thermal Management System for an Electric Machine With Additively Manufactured Hollow Conductors With Integrated Heat Pipes

Here, this paper discusses steps taken to develop a novel thermal management system for an aircraft propulsion electric machine containing additively manufactured coils integrated with heat pipes aimed at boosting its specific power. Experimental setups are used to size and characterize heat pipes for the application and 3D thermal finite element analysis is used to determine optimum heat transfer coefficient of the convective boundaries. For some of the convective boundaries, fin-based surface area enhancement is required to reach a target combined overall heat transfer coefficient and surface area performance (UA). This enhancement is worked out using a combination of the Engineering Equation Solver tool and 3D thermal FEA. The thermal management system's UA, and by extension its specific power, sensitivity to coolant temperature is explored. Temperature distribution plots of optimized machine components are also presented and discussed. Lastly, additional heat pipe testing is carried out to study its maximum heat transfer capability's sensitivity to condenser coolant temperature and configuration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Modeling Tool Development and Validation for Solar Industry Process Heat Using Particle Thermal Energy Storage

U.S. industry sectors used 26.2 quadrillion Btu and accounted for 33% of total energy consumption in 2021 according to the Energy Information Agency. Industrial process heat accounts for 70% of industrial energy use with application temperatures ranging from 60 degrees -1100 degrees C. Industry processes, heavily relying on fossil fuels of cheap coal or natural gas, differ widely in operating conditions and load requirements which makes them difficult to standardize and imposes great challenges in decarbonization. Industry processes require reliable energy supply and vary widely in temperature ranges. Storing energy from renewable sources is necessary to improve reliability and to mitigate renewable intermittency when replacing carbon fuel-based heat supplies to achieve energy savings and reduce emissions. To this end, we have developed a particle-based thermal energy storage (TES) technology using low-cost and highly stable silica sand as a storage medium. The economic and performance-based analysis is key for renewable energy sources to reliably supply industry process heat and ultimately displace fossil fuels for decarbonization. The diversified industrial processes need case-by-case analysis and design. Therefore, an adaptive modeling tool is key for renewable power with energy storage to meet industry demands. Thus, a modeling tool to simulate a solar industry process heat system using the particle TES has been developed using the object-oriented equation-based language Modelica and the commercial platform of Modelon Impact. The Modelica-based software tool provides a general simulation environment for the design of reliable solar energy sources integrated with TES for various industrial process applications at different temperatures for economic competence with fossil fuels such as coal and natural gases. It uses both customized and standard component modeling modules in Modelon libraries for the flexibility to be adapted to a specific energy demand application. The particle TES system establishes a uniform energy supply platform with an efficient heat exchanger and particle thermal energy reservoir integrated with renewable powers. The particle TES system can provide a wide temperature range and can have a large storage temperature difference that increases storage energy density; therefore, it can be an adaptable energy storage system integrated with renewable power to supply 24/7 heat for industry decarbonization.

concentrated solar thermal↗

Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. In this paper, we present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method. The numerical scheme is benchmarked by simulating the pitch angle scattering of a particle beam, and comparing with the analytical solution. Benchmark results show excellent agreement with theoretical predictions, showcasing the remarkable long time accuracy of the proposed algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

36 MATERIALS SCIENCE↗

The Rankine-Kirchhoff approximations for moist thermodynamics

The Rankine-Kirchhoff equation for the saturation vapor pressure has been derived independently many times, sometimes with the presumption of novelty. As a result, the equation is referred to in the literature by a bewildering number of different names and, often, by no name at all. In this work, the history of the Rankine-Kirchhoff equation is reviewed. The three approximations used to derive it (ideal gas, constant heat capacities, and zero-volume condensates) can be used to derive a whole system of accurate and analytic expressions of relevance to the atmosphere. Consequently, it is proposed that this set of approximations be elevated to a named entity, e.g., the Rankine-Kirchhoff approximations.

54 ENVIRONMENTAL SCIENCES↗

Numerical Modeling of Atmospheric Rime Ice Accretion on an Airfoil Using an Eulerian Approach

This paper presents an approach to numerically simulate the inherently unsteady rime ice accretion problem on a two-dimensional airfoil and elucidate the associated variations under different icing conditions. The airflow field and the water impingement on the airfoil are obtained based on an Eulerian two-phase model. A dynamic mesh strategy is employed to unsteadily account for the changes in the ice profile and its impact on the air and droplet flow by continuously reconstructing the computational grid at each time-step through smoothing and layering mechanisms. All main icing modules including the airflow field, droplet trajectory, icing thickness profile, and mesh management are fully coupled within the same computational framework without resorting to any external tools. Classical icing theory is employed to model the rime ice roughness, and it is assumed that the ice accretes in a direction normal to the airfoil surface. The governing Reynolds-averaged Navier–Stokes (RANS) conservation equations along with the energy and continuity equations are solved to produce the velocity and temperature fields. A convective film heat transfer coefficient is computed based on the surface heat flux and a recovery temperature which takes into account the dissipative heat release in the boundary layer in the vicinity of the airfoil surface. With the implemented strategy and calculating the convective heat transfer coefficient, the water film thickness is also calculated along with the ice shape. Furthermore, the model is validated by comparing the local collection efficiency distribution and ice shape with experimental data, and the results show that the implemented approach provides acceptable predictions of ice accretion profiles and rates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A control-oriented combustion model framework for compression ignition engines operating on low-reactivity fuel

This work focuses on zero-dimensional modeling of the heat release rate in a compression ignition engine operating on gasoline-like fuels. Due to the properties of gasoline, such as high volatility and longer ignition delay than diesel, the injection strategies can vary significantly from the operation with conventional diesel fuel. Different injection strategies are commonly used to achieve varying degrees of in-cylinder stratification in order to shape the combustion event and maximize efficiency. The proposed zero-dimensional combustion model was developed to account for the different stages in combustion caused by the fuel stratification. As the ignition delay model is an integral part of the entire combustion process and significantly affects the prediction accuracy, special attention has been paid to local phenomena influencing ignition delay. A one-dimensional spray model by Musculus and Kattke was employed in conjunction with a Lagrangian tracking approach in order to estimate the local air–fuel ratio within the spray tip, as a proxy for reactivity. The local air–fuel ratio, in-cylinder temperature and pressure were used in an integral fashion to estimate the ignition delay. Heat release rates were modeled using first-order non-linear differential equations. The proposed combustion model was validated against experimental data of a heavy-duty compression ignition engine with up to three injection events at mostly 1038 r/min and 14 bar brake mean effective pressure. Further validation of the model was carried out at other engine loads and speeds. Model prediction errors in CA50 of less than 1 °CA across all conditions were found. Modeling results of other combustion metrics such as combustion duration and indicated mean effective pressure are also highly satisfactory. In addition, the model has been shown to be capable of estimating the ringing intensity for most conditions.

Pamminger, Michael↗

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗

Near-Field Radiative Heat Transfer Eigenmodes

The near-field electromagnetic interaction between nanoscale objects produces enhanced radiative heat transfer that can greatly surpass the limits established by far-field blackbody radiation. Here, we present a theoretical framework to describe the temporal dynamics of the radiative heat transfer in ensembles of nanostructures, which is based on the use of an eigenmode expansion of the equations that govern this process. Using this formalism, we identify the fundamental principles that determine the thermalization of collections of nanostructures, revealing general but often unintuitive dynamics. Overall, our results provide an elegant and precise approach to efficiently analyze the temporal dynamics of the near-field radiative heat transfer in systems containing a large number of nanoparticles.

74 ATOMIC AND MOLECULAR PHYSICS↗

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

Melting rate correlation with batch properties and melter operating conditions during conversion of nuclear waste melter feeds to glasses

The rate of conversion of nuclear waste melter feed to glass is affected by the selection of melter feed materials and by melter design and operation. The melting rate correlation (MRC) is an equation that relates the glass production rate with two types of variables: (1) feed and melt properties: conversion heat, cold-cap bottom temperature, and glass melt viscosity; and (2) melter design and operation parameters: melter geometry, melter operating temperature, and gas bubbling rate. The MRC shows good agreement for an extended melting-rate data set of high-level waste (HLW) melter feeds and a data set generated for low-activity waste (LAW) melter feeds. Laboratory observation of heated melter feed samples is often used to assess the cold-cap bottom temperature of HLW melter feeds (moderately foaming feeds), but this technique appears inadequate for LAW melter feeds (vigorously foaming feeds). For LAW feeds, an adequate assessment of the cold-cap bottom temperature was achieved using evolved gas analysis, which allows identification of the collapse of primary foam for oxidized feeds. This assessment shows that the cold-cap bottom temperature for vigorously foaming LAW feeds is higher than that for moderately foaming HLW feeds. When the results of MRC are compared, LAW feeds are generally less sensitive to the bubbling rate and melt viscosity, and more sensitive to the cold-cap bottom temperature than HLW feeds. The MRC qualifies as a promising tool to support the selection of melter feed materials and melter operating conditions, which is determined from expensive independent scaled melter experiments, and sophisticated mathematical models.

Lee, Seung Min↗

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

36 MATERIALS SCIENCE↗

Design optimization of gap distance for the capillary limitation of a heat pipe with annular-type wick structure

In this study, an experimental investigation was conducted on the rising height and contact angle of fluid in an annular wick-type heat pipe. The annular wick-type heat pipe was characterized by a small gap between the wick structure and tube wall, which compensated for the pressure drop along the porous media and created additional capillary force. To describe and model the advantage of this gap, the rising of a wetting liquid in the gap between a vertical solid plate and a mesh (with a small angle between them) was experimentally measured and analyzed. An additional experiment was performed to investigate the effect of curvature on the capillary rise using tubes and meshes of varying radii. Resultantly, we confirmed that the linear combination of the contact angles of the solid plate and mesh could be applied to calculate the rising height from the Laplace–Young equation. Furthermore, the effect of curvature on the rising height of the liquid was negligible. These results were extended to the investigation of finding the optimal gap distance for the annular wick-type heat pipe by referring to previous studies. We observed that a gap distance of 1.27 mm provided the largest permeability (K) over the effective pore radius (r eff ) value for a heat pipe with ethanol, which in turn resulted in the highest capillary limitation. For a sodium heat pipe, a gap distance of 0.84 mm resulted in the highest capillary limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Equation of state for 304L stainless steel

Here, a new equation of state for 304L stainless steel has been constructed and entered in the Los Alamos National Laboratory SESAME database with material ID 4273. The equation of state is based on a multiphase approach. The solid and liquid phases were fitted to available experimental data for thermal expansion, heat capacity, and shock Hugoniot, as well as to data from density functional theory calculations of the cold curve and quantum molecular dynamics of the liquid regime. The electronic thermal contribution over the entire equation of state range is based on Tartarus (Green’s functions) average atom calculations. For this low-carbon stainless steel both the density functional and electronic thermal calculations are comprised of 19 wt% Cr, 10 wt% Ni, and the remainder Fe.

36 MATERIALS SCIENCE↗

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems using MOOSE

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗