Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal equilibrium”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiphysics analysis of fuel fragmentation, relocation, and dispersal susceptibility–Part 3: Thermal hydraulic evaluation of large break LOCA under high-burnup conditions

Increasing the peak rod average burnup of pressurized water reactor (PWR) fuel beyond 62 GWd/tU may increase fuel fragmentation, relocation, and dispersal (FFRD) susceptibility during a large break loss of coolant accident (LBLOCA). TRACE thermal hydraulic (TH) LBLOCA analyses were performed for a realistic 24-month high-burnup PWR equilibrium cycle, to inform subsequent transient BISON high-burnup FFRD susceptibility evaluations. Realistic LBLOCA systems behavior was first established by configuring to and comparing with the BEMUSE OECD LBLOCA benchmark. Fuel and operating conditions were then applied from high-burnup VERA depletion calculations. LBLOCA simulations were performed for 281 selected high-burnup rods, for which transient TH boundary conditions were collected for later use in BISON. The TRACE results indicated that rod linear heat rate (rather than burnup) is the main predictor of peak cladding temperature (PCT) during the event. PCT typically occurred at a local burnup lower than the rod-average burnup, especially for twice-burned fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

ARC disruption physics and strategy

Commonwealth Fusion Systems (CFS) plans to operate a tokamak power plant called ARC in the early 2030s. Tokamak plasmas have stability limits that, if crossed, lead to a rapid termination of the plasma, referred to as a disruption. Disruptions pose a melt risk to the first wall resulting from thermal and non-thermal particle heat fluxes, and an electromagnetic loading risk on all metal components within the equilibrium coils. A comprehensive set of models is used herein to provide an assessment of both mitigated and unmitigated ARC disruption loads. A preliminary massive gas injection system is baselined and a runaway electron mitigation coil option is proposed to close possible gaps in the baseline. It is predicted that all ARC disruption loads are within a factor of 2 of the disruption loads in SPARC, a tokamak presently under construction by CFS, and therefore SPARC provides an opportunity to calibrate models, test solutions and inform the design of ARC. The goal for ARC is disruption-free operation, however, the pragmatic design target is to withstand one mitigated disruption per day, and to restart the plasma following mitigation in tens of seconds without interrupting the power output. Unmitigated disruptions must be rare, and experience with unmitigated disruption impacts in SPARC will better define what rare means. The implications of this strategy for plasma disruptivity and disruption prediction are discussed, and operating the ARC scenario on SPARC is expected to refine the ARC final design and operational plan.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Using Artificial Microstructures to Understand Microstructure Property Relationship-Toughening Mechanisms in Metallic Glass (Final Report)

Metallic glasses are a new class of structural materials which exhibit exciting mechanical properties including high strength and elasticity. In terms of fracture toughness, the material class of metallic glasses spans a wide range; Some metallic glasses are extremely brittle and exhibit near ideal brittle behavior whereas others can be exceptional tough with values comparable to the toughest metals out there. Such large range of observed fracture toughness within the material class of metallic glasses is surprising as they have seemingly a very similar atomic structure. Therefore, we developed “artificial microstructures” which allows to decouple the various contribution of sample geometry, imperfection and structure. Specifically, we decouple variations in the alloys’ chemical composition and the atomic structure and quantified the resulting fracture toughness. Atomic structure of a metallic glass can be modified by the fictive temperature. The fictive temperature of a glass is the temperature at which the liquid metallic glass falls out of equilibrium upon colling and forms a glass. Upon further cooling the structure is maintained only thermal oscillations decrease due to a lower absolute temperature. We found that the effect of fictive temperature (same chemistry, different structural stages of the glass) is comparable to the variations of fracture toughness when the chemistry is varied. Hence, it appears that the subtle differences in the glass structure are responsible for the large range of fracture toughness’s observed. Our results reveal that fracture toughness within the material class of metallic glasses varies significantly and we found some example alloys with exceptional high resistance to fracture and others that are almost ideally brittle. Significant influences on the fracture toughness have the structure of the glass, its chemistry, and some imperfections in the structure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Smaller than Expected Bright-spot Offsets in Spitzer Phase Curves of the Hot Jupiter Qatar-1b

We present Spitzer full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature—and intermediate surface gravity and orbital period—as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse depths of 0.21 ± 0.02% at 3.6 μm and 0.30 ± 0.02% at 4.5 μm, corresponding to dayside brightness temperatures of 1542{sub −31}{sup +32} K and 1557{sub −36}{sup +35} K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are 1117{sub −71}{sup +76} and 1167{sub −74}{sup +69} K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of 0.12{sub −0.16}{sup +0.14} and a moderate day–night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with Spitzer full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: 11° ± 7° at 3.6 μm and −4° ± 7° at 4.5 μm. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Giant Outer Transiting Exoplanet Mass (GOT ‘EM) Survey. II. Discovery of a Failed Hot Jupiter on a 2.7 Yr, Highly Eccentric Orbit

Radial velocity (RV) surveys have discovered giant exoplanets on au-scale orbits with a broad distribution of eccentricities. Those with the most eccentric orbits are valuable laboratories for testing theories of high-eccentricity migration. However, few such exoplanets transit their host stars, thus removing the ability to apply constraints on formation from their bulk internal compositions. We report the discovery of Kepler-1704 b, a transiting 4.15 M {sub J} giant planet on a 988.88 day orbit with an extreme eccentricity of 0.921{sub -0.015}{sup +0.010}. Our decade-long RV baseline from the Keck I telescope allows us to measure the orbit and bulk heavy-element composition of Kepler-1704 b and place limits on the existence of undiscovered companions. A failed hot Jupiter, Kepler-1704 b was likely excited to high eccentricity by scattering events that possibly began during its gas accretion phase. Its final periastron distance was too large to allow for tidal circularization, so now it orbits its host from distances spanning 0.16–3.9 au. The maximum difference in planetary equilibrium temperature resulting from this elongated orbit is over 700 K. A simulation of the thermal phase curve of Kepler-1704 b during periastron passage demonstrates that it is a remarkable target for atmospheric characterization from the James Webb Space Telescope, which could potentially also measure the planet’s rotational period as the hot spot from periastron rotates in and out of view. Continued characterization of the Kepler-1704 system promises to refine theories explaining the formation of hot Jupiters and cool giant planets like those in the solar system.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Stable isotope equilibria in the dihydrogen-water-methane-ethane-propane system. Part 2: Experimental determination of hydrogen isotopic equilibrium for ethane-H2 from 30 to 200 °C and propane-H2 from 75 to 200 °C

The stable isotopic compositions of light n-alkanes, including methane, ethane, and propane, are often used to identify the sources and thermal maturity of natural gas samples. Though stable isotopic compositions of these molecules are commonly assumed to be controlled by kinetic isotope effects, recent studies have proposed both carbon and hydrogen isotopic equilibrium may also occur in some samples. Assessing whether samples are in isotopic equilibrium requires knowledge of light alkane equilibrium fractionation factors over geologically relevant temperatures for formation and storage (up to ∼300 °C). In this study, we report experimental results of hydrogen isotopic equilibrium between ethane and H2 from 30 to 200 °C and propane and H2 from 75 to 200 °C. We compare these results with high-level theoretical calculations and provide a preferred polynomial fit to describe equilibrium fractionation factors. Comparison of these fractionation factors with a compilation of ∼500 compiled environmental gas samples supports the proposal that many (∼50%) of these natural gas samples exhibit hydrogen isotopic compositions consistent with having formed in or attained methane-ethane-propane hydrogen isotopic equilibrium over geologically relevant temperatures for formation and storage (50–300 °C).

Turner, Andrew C↗
DOE OSTI2025

Bulk viscosity of nuclear matter with pions in the neutrino-trapped regime

Recent work has shown that a population of thermal pions could modify the equation of state and transport properties of hot and dense neutron-rich matter and introduce new reaction pathways to change the proton fraction. In this article we study their impact on the bulk viscosity of dense matter, focusing on the neutrino-trapped regime that would be realized in neutron star mergers and supernovae. We find that the presence of a thermal population of pions alters the bulk viscosity by modifying the EoS (via the susceptibilities) and by providing new reaction pathways to achieve beta equilibrium. In neutron star merger conditions, the bulk viscosity in neutrino-trapped npeμ matter (without pions) has its peak at temperatures of at most a couple MeV and is quite small at temperatures of tens of MeV. We find that thermal pions enhance the low-temperature peak of the bulk viscosity by a factor of a few and shift it to slightly lower temperatures. At higher temperatures, where the pion abundance is large but the bulk viscosity is traditionally small, pions can increase the bulk viscosity by an order of magnitude or more, although it is still orders of magnitude smaller than its peak value.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Coupling between bulk thermal defects and surface segregation dynamics

Surface segregation is a phenomenon that depends on the delicate interplay between thermodynamic driving forces and kinetic obstacles, for which elevated temperature is often needed to enhance the atom mobility and reach equilibrium. Using the classic system of Cu 3 Au ( 100 ) under nonisothermal conditions, in this study, we show an adatom process underlying transient surface segregation dynamics through the temperature-change-driven creation and annihilation of thermal vacancies in the bulk and the resulting bulk-surface mass exchanges. This is demonstrated by monitoring the surface composition evolution of Cu 3 Au ( 100 ) with temperature changes between 250 and 500 ° C , showing that the increase in temperature decreases monotonically the surface Au concentration as a result of the transfer of more Cu than Au from the bulk to the surface to form Cu-rich clusters of adatoms. Such a bulk thermal defect effect is expected to be universal in inducing the disparity in the bulk-surface mass exchanges of dissimilar atoms in multicomponent materials because of the inherent differences in the vacancy formation energies of the constituent atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗
DOE OSTI2022

Oxygen Isotope Fractionation in U 3 O 8 during Thermal Processing in Humid Atmospheres

The incorporation of oxygen isotopes from water into uranium oxides during industrial processing presents a pathway for determining a material’s geographical origin. This study is founded on the hypothesis that oxygen isotopes from atmospheric water vapor will exchange with isotopes of oxygen in solid uranium oxides during thermal processing or calcination. Using a commonly encountered oxide, U 3 O 8 , the exchange kinetics and equilibrium fractionation with water vapor (in a concentration range of 50–55% relative humidity) were investigated using processing temperatures of 400, 600, and 800 °C. In an atmosphere containing only water vapor diluted in N 2 , oxygen isotope equilibration in U 3 O 8 occurred within 12 h at 400 °C and within 2 h at 600 and 800 °C. Fractionation factors (1000lnα, U 3 O 8 –H 2 O) between the water and oxide were –12.1, –11.0, and –8.0 at 400, 600, and 800 °C, respectively. With both humidity and O 2 present in the calcining atmosphere, isotopic equilibration is attained within 2 h at and above 400 °C. In this mixed atmosphere, which was designed to emulate Earth’s troposphere, isotopes are incorporated preferentially from water vapor at 400 °C and from O 2 at 600 and 800 °C. Rapid and temperature/species-dependent isotope exchange also elucidated the impact of retrograde exchange in humid air, showing a shift from O 2 -dependent to H 2 O-dependent fractionation as U 3 O 8 cooled from 800 °C. These results confirm that uranium oxides inherit oxygen isotopes from humidity during thermal processing, illuminating an important mechanism in the formation of this forensic signature.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SAM Two-Phase Flow Model Development and Applications for Operational Transients in Advanced Reactors

As advanced nuclear technologies continue to develop, the need for the flexible operation and generation of these advanced reactors becomes necessary to maximize economic potential. As large-scale experiments are not always feasible, modeling and simulations of advanced reactors play a crucial role in design optimization and analysis. The SAM (System Analysis Module) code developed at Argonne National Laboratory is a state-of-the-art system-level thermal-hydraulic code aimed at simulating advanced reactor systems. Recent code developments have implemented two-phase flow modeling using the homogeneous equilibrium model, and a new steam generator component has been developed to utilize the two-phase flow implementation. In addition to verification tests, a load-following simulation was performed to model a realistic load-following transient in a proposed integrated system consisting of a conceptual advanced reactor known as the Advanced Burner Test Reactor (ABTR) and thermal energy storage (TES) tanks. The integrated system model uses two large TES tanks designed for sodium and a model helical coil steam generator to simulate the operational load-following transient. The flow rates of the feedwater and secondary loops are regulated to meet a prescribed steam generator load consistent with the electricity demand over a 24-h period. In conclusion, the results found the ABTR system was able to maintain stable reactor conditions and primary- and secondary-side characteristics over the course of the load-following transient.

Advanced Burner Test Reactor (ABTR)↗
DOE OSTI2021

Examining normal modes as fundamental heat carriers in amorphous solids: The case of amorphous silicon

Normal mode decomposition of atomic vibrations has been used to provide microscopic understanding of thermal transport in amorphous solids for decades. In normal mode methods, it is naturally assumed that atoms vibrate around their equilibrium positions, and that individual normal modes are the fundamental vibrational excitations transporting heat. With the abundance of predictions from normal mode methods and experimental measurements now available, we carefully analyze these calculations in amorphous silicon, a model amorphous solid. We find a number of discrepancies, suggesting that treating individual normal modes as fundamental heat carriers may not be accurate in amorphous solids. Furthermore, our classical and ab initio molecular dynamics simulations of amorphous silicon demonstrate a large degree of atomic diffusion, especially at high temperatures, leading to the conclusion that thermal transport in amorphous solids could be better described starting from the perspectives of liquid physics rather than from crystalline solids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modified Calcium Manganites for Thermochemical Energy Storage Applications

CaAl 0.2 Mn 0.8 O 3-δ (CAM28) and CaTi 0.2 Mn 0.8 O 3-δ (CTM28) are perovskite metal oxides developed for high-temperature thermochemical energy storage (TCES) applications, e.g., in support of air Brayton power generation. Previous reports for these compounds focus on the equilibrium non-stoichiometry (δ) as a function of temperature and oxygen partial pressure (pO 2 ) and the endotherm (or exotherm) accompanying changes in δ resulting from thermal reduction (or re-oxidation). Herein, we report results for elemental substitution and doping (Al, Co, Fe, La, Sr, Ti, Y, Zn, and Zr) of calcium manganites (CM) that establish the preference for CAM28 and CTM28. Techniques employed include conventional (screening and equilibrium) and ballistically heated multi-cycle thermogravimetric analysis (TGA), conventional and high temperature ( in-situ ) X-ray diffraction (XRD), and differential scanning calorimetry (DSC). Forward-looking results for A-site Y-doped materials, e.g., Ca 0.9 Y 0.1 MnO 3-δ (CYM910), establish a route to increasing the reduction enthalpy relative to CAM28 and CTM28, albeit at the expense of increased reduction temperatures and raw materials costs. A thermodynamic model presented for CAM28, but extendable to related materials, provides values for the reaction enthalpy and extent of reduction as a function of temperature and oxygen partial pressure for use in design efforts. Taken as a whole, the results support the choice of Al-doped CaMnO 3-δ as a low-cost material for TCES in a high temperature air Brayton application, but point the way to achieving higher stored energy densities that could lead to overall cost savings.

Miller, James E.↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗