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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Study of Temperature Distribution and Thermal Stresses in a Hot Rock Due to Rapid Cooling

Thermal stresses may be induced in a hot dry rock when a cold fluid is injected in the well. To study this problem, we look at the thermoelastic response of a hot rock that is suddenly cooled. The cooling is assumed to be either at a constant temperature or at a constant heat flux per unit depth. Our approach is to nondimensionalize the equations and perform a parametric study and look at the temperature distribution and the induced-thermal stresses. The results indicate that depending on the extent of cooling and the cooling time, thermal stresses can be induced. Numerical simulations on sandstone, with an initial uniform temperature of 473 K, are also carried out. The results show that if the cooling is due to the surface temperature maintained at 463 K (10 °C lower than the initial temperature of the hot rock), thermal stresses that are larger than the rock tensile strength could be induced. When the cooling is due to a constant surface heat flux, this temperature can be reached after about 777 days of cooling with a minimum value of a heat flux of -20 W/m.

Tran, Phuoc↗
DOE OSTI2026

Compression behavior of lithium fluoride up to 80 GPa and 2300 K

Here, the high-temperature compression behavior of lithium fluoride (LiF) has been determined to ∼80 GPa and 2300 K by means of high-resolution synchrotron-based x-ray diffraction in a laser-heated diamond-anvil cell. The room-temperature Vinet equation of state (EOS) for LiF yields an isothermal room-temperature, ambient-pressure lattice parameter a 0 = 4.028 (±0.003) Å, volume V 0 = 65.35 (±0.03) Å 3 , bulk modulus K 0 = 67.57 (±0.34) GPa, and its pressure derivative K 0 ’ = 4.64 (±0.03). This expanded experimental dataset is in good agreement with the LiF thermal equation of state by Myint et al. (2019) up to ∼2300 K and 25 GPa. At higher pressures and temperatures our data shows reduced thermal expansion compared to the prediction and may be particularly important for experiments at more extreme conditions.

High-pressure behavior↗
DOE OSTI2022

Theory of the magnetothermal instability in coronal plasma flows

The theory of the magnetothermal instability (MTI) is revisited through the lens of the stability of uniform systems. The linear stability analysis includes flow advection and Nernst transport. The instability criteria derived distinguish between the convective and the absolute nature of the perturbation growth. It is proven that, in the region where the Nernst and plasma blowoff velocities cancel, the MTI can be absolute and wave-packet perturbations grow in situ. This instability is mediated by the internal feedback between the Biermann battery and Righi–Leduc terms. Here, the analysis is extended to derive the dispersion relation for short-wavelength perturbations developing in nonuniform profiles with the application to coronal plasmas. It is found that the condition for MTI requires the net B-field convection velocity to be small at the isothermal sonic section, and the plasma conditions in this section govern the dynamics of the instability. Analysis of hydro-equivalent implosions suggests that unstable perturbations undergo more e-foldings of growth in larger-size targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Time-Resolved Stochastic Dynamics of Quantum Thermal Machines

Steady-state quantum thermal machines are typically characterized by a continuous flow of heat between different reservoirs. However, at the level of discrete stochastic realizations, heat flow is unraveled as a series of abrupt quantum jumps, each representing an exchange of finite quanta with the environment. Here, in this work, we present a framework that resolves the dynamics of quantum thermal machines into cycles classified as enginelike, coolinglike, or idle. We analyze the statistics of individual cycle types and their durations, enabling us to determine both the fraction of cycles useful for thermodynamic tasks and the average waiting time between cycles of a given type. Central to our analysis is the notion of intermittency, which captures the operational consistency of the machine by assessing the frequency and distribution of idle cycles. Our framework offers a novel approach to characterizing thermal machines, with significant relevance to experiments involving mesoscopic transport through quantum dots.

full counting statistics↗
DOE OSTI2022

Modelling of edge plasma dynamics with active wall boundary conditions

A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Heat rectification through single and coupled quantum dots

We study heat rectification through quantum dots in the Coulomb blockade regime using a master equation approach. We consider both cases of two-terminal and four-terminal devices. In the two-terminal configuration, we analyze the case of a single quantum dot with either a doubly-degenerate level or two non-degenerate levels. In the sequential tunneling regime we analyze the behaviour of heat currents and rectification as functions of the position of the energy levels and of the temperature bias. In particular, we derive an upper bound for rectification in the closed-circuit setup with the doubly-degenerate level. We also prove the absence of a bound for the case of two non-degenerate levels and identify the ideal system parameters to achieve nearly perfect rectification. The second part of the paper deals with the effect of second-order cotunneling contributions, including both elastic and inelastic processes. In all cases we find that there exists ranges of values of parameters (such as the levels' position) where rectification is enhanced by cotunneling. In particular, in the doubly-degenerate level case we find that cotunneling corrections can enhance rectification when they reduce the magnitude of the heat currents. For the four-terminal configuration, we analyze the non-local situation of two Coulomb-coupled quantum dots, each connected to two terminals: the temperature bias is applied to the two terminals connected to one quantum dot, while the heat currents of interest are the ones flowing in the other quantum dot. Remarkably, in this situation we find that non-local rectification can be perfect as a consequence of the fact that the heat currents vanish for properly tuned parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

P-V-T equation of state of boron carbide

We report the P-V-T equation of state measurements of B 4 C to 50 GPa and approximately 2500 K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K 0 of 221(2) GPa and derivative, (dK/dP) 0 of 3.3(1). These were used in fits with both a Mie–Grüneisen–Debye model and a temperature-dependent, Birch–Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK 0 /dT). The ambient pressure thermal expansion coefficient (α 0 + α 1 T), Grüneisen γ(V) = γ 0 (V/V 0 ) q and volume-dependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Sources of error in detonation calorimeters and error analysis for neat 2,4,6-triamino-1,3,5- trinitrobenzene (TATB)

Here, a calorimeter for measuring heats of detonation at Lawrence Livermore National Laboratory is described. A calibration precision of 0.2 % at the 95 % confidence interval (CI) is reported. Sources of uncertainty are discussed, including nonequivalent sources, which are those arising from differences between calibration and experimental tests. The systematic error due to nonequivalent sources is bounded to 0.19–0.22 % of the measured heat for a standard detonation test where the confinement material is gold. The recommendation is to correct for the systematic error by adding 0.19 % to the reported value and adding 0.03 % to the uncertainty. It is demonstrated that the precision of a detonation test is variable with testing duration because a source of uncertainty resides in the thermodynamic correction factor k6, which accounts for the contribution to heat by stirring; the generated power is additive and therefore highly impacted by total test time. The relative proportion also varies with the magnitude of heat release and sample size, adding variance to the weight of the error arising from the correction factor. A full error analysis based on the described sources of uncertainty is developed. The methodology is applied to a test series on neat 2,4,6-triamino-1,3,5- trinitrobenzene (TATB), demonstrating an ultimate precision of 0.7 % (single test) for materials testing and a relative standard deviation of 1.8 %.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2020

Vibrational hierarchy leads to dual-phonon transport in low thermal conductivity crystals

Many low-thermal-conductivity ( κ L ) crystals show intriguing temperature ( T ) dependence of κ L : κ L ∝ T -1 (crystal-like) at intermediate temperatures whereas weak T -dependence (glass-like) at high temperatures. It has been in debate whether thermal transport can still be described by phonons at the Ioffe-Regel limit. In this work, we propose that most phonons are still well defined for thermal transport, whereas they carry heat via dual channels: normal phonons described by the Boltzmann transport equation theory, and diffuson-like phonons described by the diffusion theory. Three physics-based criteria are incorporated into first-principles calculations to judge mode-by-mode between the two phonon channels. Case studies on La 2 Zr 2 O 7 and Tl 3 VSe 4 show that normal phonons dominate low temperatures while diffuson-like phonons dominate high temperatures. Our present dual-phonon theory enlightens the physics of hierarchical phonon transport as approaching the Ioffe-Regel limit and provides a numerical method that should be practically applicable to many materials with vibrational hierarchy.

42 ENGINEERING↗
DOE OSTI2022

Experimental Study of Eutectic Molten Salts NaCl/KCl/ZnCl2 Heat Transfer Inside a Smooth Tube for High-Temperature Application

Abstract Eutectic salts NaCl-KCl-ZnCl2 and NaCl-KCl-MgCl2 are two of the chloride salt systems that are promising for being used as high-temperature heat transfer fluid (HTF) and thermal energy storage (TES) materials in a wide range of temperatures from 230 °C to 750 °C in concentrated solar thermal power systems. To conduct measurement of the heat transfer coefficient of the molten salt at high temperatures, a circulation system and instrumentation of flow and heat transfer was constructed. Experimental measurement of the convective heat transfer coefficients of NaCl-KCl-ZnCl2 (molar fraction: 13.8–41.9–44.3%) inside tubes has been accomplished to find the most suitable heat transfer correlations. Experience of salt processing and operation of the high-temperature heat transfer test system was obtained. Two correlations, Dittus-Boelter equation and Gnielinski’s correlation for Nusselt number against Reynolds number and Prandtl number, are evaluated using the test results, and the latter correlation is recommended due to its better agreement of prediction against tested data.

Energy & Fuels↗
DOE OSTI2020

New Yield Estimates for Nuclear Detonations Over Water

One of the primary methods used by the United States to estimate the yields of atmospheric nuclear tests performed in the 1940’s, 50’s, and 60’s was to film each detonation with a high-speed camera, measure the blast radius as a function of time, and then use Taylor’s spherical blast wave equation to calculate the source energy (i.e., the yield). For surface detonations, the yield in Taylor’s equation was assumed to be twice as large as the actual yield to account for the energy reflected from the surface, assuming the ground was a perfect reflector. Recent simulations of surface detonations over water (referred to as barge shots) have shown that Taylor’s equation must be modified further to account for material entrainment and heat losses. We have computed this lost energy using a high-order radiation-hydrodynamics code to determine the reduction in blast radius vs. time. By analyzing simulation data in the same manner as film data, we have determined that a yield correction factor of 1.27± 0.04 should be applied to the barge shots.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Description of Sensor Assignment Optimization Method as Deployed on a Multi-Node Cluster

Data analytic methods are being developed to address the problem of how to assign a sensor set in a nuclear facility such that a requisite level of process monitoring capability is realized and that the sensor set is sufficiently rich to determine the status of the individual sensors with respect to need for calibration. There is an awareness in the nuclear industry that data analytics combined with rich sensor sets represent a means to improve operations and reduce costs. In the industry the calibration problem has been previously approached as an empirical data-driven problem with several methods having been developed. However, the experience of the utilities over the past ten years with these methods indicates that the absence of physics-based information renders the data-driven approach less reliable. Complicating factors such as the inherent variability of operation (both equipment alignment and operating condition) can confound a pure data-driven approach while there are no rigorous guidelines for determining what constitutes an adequate sensor set. The solution under development to overcome these shortcomings supplements the data analytic method with process information in a so-called process-constrained data-analytic approach. Simple balance equations are written for generic components (e.g., mechanical pump, valve, and heat exchanger). These do not require a priori knowledge of process parameters, such as heat transfer coefficients or friction factors. All that is needed on the part of the utility user is to identify the components and how they are connected. This report describes the development of a parallel computing capability for determining the optimal sensor set. The optimal sensor set problem suffers from the curse of dimensionality. Computation time increases exponentially as the size of the system grows. To overcome this difficulty a pre-conditioner algorithm is developed to find an approximate solution close the actual solution. This serves as a seed for the full-blown algorithm and acts to constrain the space that must searched. The optimization algorithms are described and the implementation on a parallel computing platform is described. The application of the method to a use case we are solving in collaboration with our utility partner served to illustrate how the default sensor set in a nuclear plant may not provide sufficient coverage to infer sensor calibration status.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Modeling the cooldown of cryocooler conduction-cooled devices

Cryocooler conduction cooled devices can experience significant cooldown time due to lower available cooling capacity compares to convection cooled devices. Therefore, the cooldown time is an important design parameter for conduction cooled devices. This article introduces a framework developed in Python for calculating the cooldown profiles and cooldown time of cryocooler conduction-cooled devices such as superconducting magnets and accelerator cavities. The cooldown time estimation problem is essentially a system of ordinary first-order differential equations comprising the material properties (temperature dependent thermal conductivity and specific heat capacity) of the components intertwined with the prevailing heat transfer channels (conduction, radiation, and heat flow across pressed contacts) and the cryocooler capacity. The formulation of this ODE system is first presented. This ODE system is then solved using the in-built Python library odeint. A case study is presented comprising a small cryocooler conduction-cooled copper stabilized niobium-titanium magnet. The case study is supplemented with the Python script enabling the reader to simply tweak the device design parameters and optimize the design from the point of view of slow/fast cooldown.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Undersizing of aged African biomass burning aerosol by an ultra-high-sensitivity aerosol spectrometer

The ultra-high-sensitivity aerosol spectrometer (UHSAS) differs from most other optical particle spectrometers by using a high-power infrared (IR) laser to detect small particles and reduce the sizing ambiguity due to the non-monotonicity of scattering with particle size. During the NASA ORACLES project (ObseRvations of Aerosols above CLouds and their intEractionS) over the southeast Atlantic Ocean, the UHSAS clearly undersized particles in the biomass burning plume extending from southern Africa. Since the horizontal and vertical extent of the plume was vast, the NASA P-3B research aircraft often flew through a fairly uniform biomass burning plume for periods exceeding 30 min, sufficient time to explore the details of the UHSAS response by selecting single particle sizes with a differential mobility analyzer (DMA) and passing them to the UHSAS. This was essentially an in-flight calibration of the UHSAS using the particles of interest. Two modes of responses appeared. Most particles were undersized by moderate amounts, ranging from not at all for 70 nm aerosols to 15 % for 280 nm particles. Mie scattering calculations show that composition-dependent refractive index of the particles cannot explain the pattern. Heating of brown carbon or tarballs in the beam causing evaporation and shrinking of the particles is the most plausible explanation, though mis-sizing due to non-sphericity cannot be ruled out. A small fraction (10 %–30 %) of the particles were undersized by 25 % to 35 %. Those were apparently the particles containing refractory black carbon. Laboratory calibrations confirm that black carbon is drastically undersized by the UHSAS, because particles heat to their vaporization point and shrink. A simple empirical correction equation was implemented that dramatically improves agreement with DMA distributions between 100 and 500 nm. It raised the median particle diameter by 18 nm, from 163 to 181 nm, during the August 2017 deployment and by smaller amounts during deployments with less intense pollution. Calculated scattering from UHSAS size distributions increased by about 130 %, dramatically improving agreement with scattering measured by nephelometers. The correction is only valid in polluted instances; clean marine boundary layer and free troposphere aerosols behaved more like the calibration spheres. We were unable to directly test the correction between 500 and 1000 nm, though aerodynamic particle sizer (APS) data appear to show that the correction is poor at the largest diameters, which is no surprise as the composition of those particles is likely to be quite different than that of the accumulation mode. This adds to the evidence that UHSAS data must be treated cautiously whenever the aerosol may absorb infrared light. Similar corrections may be required whenever brown carbon aerosol is present.

54 ENVIRONMENTAL SCIENCES↗