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

Non-equilibrium molecular dynamics studies of thermal diffusion of hydrogen isotopes in low concentration zirconium hydrides

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermal and dynamic partition of dumbbell interstitials in complex concentrated alloys

Complex concentrated alloys (CCAs) are promising candidates for applications in extreme conditions, such as irradiation where interstitial mediated diffusion is important. In CCAs with N principal elements, $\frac{N(N+1)}{2}$ types of dumbbell interstitials exist. Currently, there is no way to predict the thermal partition (fractional concentration at equilibrium) and the dynamic partition (fractional time an interstitial spends during diffusion) of each type of dumbbell interstitial. To mitigate this issue, this work proposes a theoretical model for computing the equilibrium concentrations and thermal partition of dumbbell interstitials in CCAs and validates the model using grand canonical Monte Carlo simulations. Lattice kinetic Monte Carlo simulations show that the thermal partition is equivalent to the dynamic partition, and both are governed by composition and formation energies of dumbbells. As a result, the model proposed provides a foundation for understanding radiation enhanced diffusion and induced segregation in CCAs under irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thermal false vacuum decay around black holes

In flat space and at finite temperature, there are two regimes of false vacuum decay in quantum field theory. At low temperature, the decay proceeds through thermally-assisted tunneling described by periodic Euclidean solutions—bounces—with nontrivial time dependence. On the other hand, at high temperatures the bounces are time independent and describe thermal jumps of the field over the potential barrier. Here, we argue that only solutions of the second type are relevant for false vacuum decay catalyzed by a black hole in equilibrium with thermal bath. The argument applies to a wide class of spherical black holes, including d-dimensional anti–de Sitter/de Sitter Schwarzschild black holes and Reissner-Nordström black holes sufficiently far from criticality. It does not rely on the thin-wall approximation and applies to multifield scalar theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

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

Initial study on cross-section generation requirements for a PBR equilibrium core

A Serpent model of the equilibrium core HTR-PM small modular nuclear reactor in China, was developed for use in cross-section preparation studies in order to guide methods development for the Griffin reactor multiphysics application. The model includes detailed isotopics for 10 distinct pebble burnup groups in 126 core zones with unique fuel and moderator temperatures obtained from a coupled neutronics-thermal-fluids equilibrium core calculation using Griffin-Pronghorn. A sensitivity study of the fuel and moderator temperatures for various core regions was performed with the MOOSE stochastic tools. The results show that the uncertainties are, not unexpectedly, dominated by the value of the fluid temperature and that the power level, heat transfer coefficient and effective conduction to neighboring pebbles and fluid constitute, at best, second order effects. The temperature uncertainty range varies from 28 K to 57 K between the core entry and exit planes, respectively, but these values are probably higher. We still have to quantify the significance of these uncertainties in the preparation of cross-sections in future work. In addition, we verified that the effective pebble approximation used in the PEBBED and V.S.O.P. computer codes works well for the preparation of region averaged cross-sections. Nevertheless, there are some discrepancies in the cross-sections when compared to the multi-pebble model, which could affect the prediction of peak values and the depletion calculation. We conclude that is highly desirable for future studies with Griffin to be able to handle both the 'effective' pebble approximation and the multi-pebble approach for various pebble burnup groups. This enables Griffin users with the flexibility to perform higher-fidelity studies. Finally, we initiated the preparation of cross-sections for various core regions from the full core Serpent reference model. We quantified the differences in 26 group cross-sections from infinite domain models versus the full core approach. These reference cross-sections will serve to verify the double heterogeneity, self-shielding, and spectrum-correction methods in Griffin. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Alternating Current Photovoltaic Effect

It is well known that the photovoltaic effect produces a direct current (DC) under solar illumination owing to the directional separation of light-excited charge carriers at the p–n junction, with holes flowing to the p-side and electrons flowing to the n-side. In this study, we found that apart from the DC generated by the conventional p–n photovoltaic effect, there is another new type of photovoltaic effect that generates alternating current (AC) in the nonequilibrium states when the illumination light periodically shines at the junction/interface of materials. The peak current of AC at high switching frequency can be much higher than that from DC. The AC cannot be explained by the established mechanisms for conventional photovoltaics; instead, it is suggested to be a result of the relative shift and realignment between the quasi-Fermi levels of the semiconductors adjacent to the junction/interface under the nonequilibrium conditions, which results in electron flow in the external circuit back and forth to balance the potential difference between two electrodes. By virtue of this effect, the device can work as a high-performance broadband photodetector with extremely high sensitivity under zero bias; it can also work as a remote power source providing extra power output in addition to the conventional photovoltaic effect.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Construction of a Code Verification Matrix for Heat Conduction With Finite Element Code Applications

When establishing the pedigree of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model. During this process, numerical results on various meshes are systematically compared to a reference analytic solution. The selection of analytic solutions can be a laborious process, as it is difficult to establish adequate code confidence without performing redundant work. In this report we address this issue by applying a physics-based process that establishes a set of reference problems. In this process, code simulation options are categorized and systematically tested, which ensures that gaps in testing are easily identified and addressed. The resulting problems are primarily intended for code verification analysis but may also be useful for comparison to other simulation codes, troubleshooting activities, or training exercises. The process is used to select fifteen code verification problems relevant for the one-dimensional steady-state heat conduction equation. These problems are applicable to a wide variety of simulation tools, but, in this work, a demonstration is performed using the finite element-based nuclear fuel performance code BISON. Convergence to the analytic solution at the theoretical rate is quantified for a selection of the problems, which establishes a baseline pedigree for the code. Not only can this standard set of conduction solutions be used for verification of other codes, but also the physics-based process for selecting problems can be utilized to quantify and expand testing for any simulation tool.

42 ENGINEERING↗
DOE OSTI2024

Phonon modal analysis of thermal transport in ThO 2 with point defects using equilibrium molecular dynamics

Defects can significantly degrade the thermal conductivity of ThO 2 , an advanced nuclear fuel material as well as a surrogate for other fluorite-structured materials. Here, we investigate how point defects in ThO 2 impact phonon mode-resolved thermal transport. By incorporating phonon modes from lattice dynamics, we decompose the trajectory and heat flux to phonon normal mode space and extract key phonon properties, including phonon relaxation times and their contributions to thermal conductivity. We implement two methods. The first method is based on the Green Kubo formalism to resolve the contribution of each phonon mode to thermal conductivity. The second resolves the lifetime of individual phonon modes and the thermal conductivity is calculated using the Boltzmann transport equation within relaxation time approximation. Notably, a lower contribution of acoustic modes is revealed compared to perturbative approaches considering only three-phonon scattering processes. The effects of four types of point defects are evaluated. The strongest impact on a reduction in thermal conductivity is from Th interstitials, followed by Th vacancies. O interstitials/vacancies have a similar impact, albeit smaller than defects on the thorium sublattice. These observations are consistent with previous studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Towards chemical equilibrium in thermochemical water splitting. Part 1: Thermal reduction

The efficiency of many processes strongly depends on their thermodynamic reversibility, i.e., proximity to equilibrium throughout the process. In thermochemical cycles for water and/or carbon dioxide splitting, thermochemical air separation, and thermochemical energy storage, operating near equilibrium means that the oxygen chemical potential of the solid and gas phases must not differ significantly. Herein we show that approaching this ideal is possible in thermal reduction only if the reaction step occurs at a specific, reaction coordinate- and material-dependent temperature. The resulting thermal reduction temperature profile also depends on the ratio of gas and solid flows.

08 HYDROGEN↗
DOE OSTI2020

Thermomechanical response of thickly tamped targets and diamond anvil cells under pulsed hard x-ray irradiation

In the laboratory study of extreme conditions of temperature and density, the exposure of matter to high intensity radiation sources has been of central importance. Here, we interrogate the performance of multi-layered targets in experiments involving high intensity, hard x-ray irradiation, motivated by the advent of extremely high brightness hard x-ray sources, such as free electron lasers and 4th-generation synchrotron facilities. Intense hard x-ray beams can deliver significant energy in targets having thick x-ray transparent layers (tampers) around samples of interest for the study of novel states of matter and materials’ dynamics. Heated-state lifetimes in such targets can approach the microsecond level, regardless of radiation pulse duration, enabling the exploration of conditions of local thermal and thermodynamic equilibrium at extreme temperature in solid density matter. The thermal and mechanical responses of such thick layered targets following x-ray heating, including hydrodynamic relaxation and heat flow on picosecond to millisecond timescales, are modeled using radiation hydrocode simulation, finite element analysis, and thermodynamic calculations. Assessing the potential for target survival over one or more exposures and resistance to damage arising from heating and resulting mechanical stresses, this study doubles as an investigation into the performance of diamond anvil high pressure cells under high x-ray fluences. Long used in conjunction with synchrotron x-ray radiation and high power optical lasers, the strong confinement afforded by such cells suggests novel applications at emerging high intensity x-ray facilities and new routes to studying thermodynamic equilibrium states of warm, very dense matter.

Meza-Galvez, J. (ORCID:000000032699849X)↗
DOE OSTI2026

Elucidating key reducing species beyond ions in hydrogen plasma smelting reduction of iron ore

Hydrogen plasma smelting reduction (HPSR) of iron ore has attracted significant attention over the past decade due to its high-temperature operation, rapid plasma mediated reduction kinetics, and simpler density-based separation of molten iron product, compared to H2-based solid-state reduction. All of these attributes enable processing of low-grade ores for downstream use in electric-arc furnaces, as virgin iron with low gangue content is required for high quality steel and improved furnace operation. While positive ions exist within the plasma arc, this work demonstrates that near the anodic ore surface, hydrogen radicals and vibrationally excited hydrogen species dominate and their densities correlate well with observed reduction rates. Species concentrations in the transferred plasma arc and at the plasma-ore interface are evaluated using coupled thermal plasma and near-wall non-equilibrium plasma models. The thermal plasma model is validated against experimental voltage data and spectroscopic measurements of plasma temperature and density for varying current inputs. Modeling of the near surface thermochemical non-equilibrium and micrometer scale anode sheath layer reveals, in addition to the expected H + , significant concentrations of ArH + and H$^+_3$ ions, typically not observed in thermal plasmas under thermodynamic equilibrium. Our results show that the inverted sheath structure at the anodic ore surface strongly suppresses reactive positive ion fluxes, while non-equilibrium electron-impact processes generate abundant hydrogen radicals and vibrationally excited species. These findings highlight the critical role of non-equilibrium effects in hydrogen arc-driven iron ore reduction and advance understanding beyond prevailing hypotheses centered on hydrogen ion-driven mechanisms.

08 HYDROGEN↗
DOE OSTI2020

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

QCD thermalization: Ab initio approaches and interdisciplinary connections

Heavy-ion collisions at BNL's Relativistic Heavy Ion Collider and CERN's Large Hadron Collider provide strong evidence for the formation of a quark-gluon plasma, with temperatures extracted from relativistic viscous hydrodynamic simulations shown to be well above the transition temperature from hadron matter. Outstanding problems in QCD include how the strongly correlated quark gluon matter forms in a heavy-ion collision, its properties off-equilibrium, and the thermalization process in the plasma. In this work, we review the theoretical progress in this field in weak coupling QCD effective field theories and in strong coupling holographic approaches based on gauge-gravity duality. We outline the interdisciplinary connections of different stages of the thermalization process to non equilibrium dynamics in other systems across energy scales ranging from inflationary cosmology, to strong field QED, to ultracold atomic gases, with emphasis on the universal dynamics of non-thermal and hydrodynamic attractors. We survey measurements in heavy-ion collisions that are sensitive to the early non-equilibrium stages of the collision and discuss the potential for future measurements. We summarize the current state-of-the art in thermalization studies and identify promising avenues for further progress.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The coset construction for non-equilibrium systems

We propose a systematic coset construction of non-equilibrium effective field theories (EFTs) governing the long-distance and late-time dynamics of relativistic, finite-temperature condensed matter systems. Our non-equilibrium coset construction makes significant advances beyond more standard coset constructions in that it takes advantage of recently-developed techniques, which allow the formulation of non-equilibrium effective actions that account for quantum and thermal fluctuations as well as dissipation. Because these systems exist at finite temperature, the EFTs live on the closed-time-path of the Schwinger-Keldysh contour. Since the coset construction and the non-equilibrium effective actions may be unfamiliar to many readers, we include brief introductions to these topics in an effort to make this paper self-contained. To demonstrate the legitimacy of this coset construction, we successfully reproduce the known EFTs for fluids and superfluids at finite temperature. Then, to demonstrate its utility, we construct novel EFTs for solids, supersolids, and four phases of liquid crystals, all at finite temperature. We thereby combine the non-equilibrium effective action and the coset construction to create a powerful tool that can be used to study many-body systems out of thermal equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A coupled equilibrium boundary layer model with stable water isotopes and its application to local water recycling

The contribution of evapotranspiration (ET) to regional precipitation, known as “local water recycling”, is a key process in the water cycle that can affect water management. However, the amount of planetary boundary layer (PBL) moisture that arises from ET is highly uncertain due to complex atmosphere and land surface conditions. In this study, an idealized two-layer equilibrium boundary layer model was coupled with a stable water isotope module including HDO and H 2 18 O to constrain PBL growing season water transport processes. The model was validated using turbulent heat fluxes and isotope ratios of water vapor (δ m ) and precipitation (δ P ) measured at a cropland site and a nearby tall tower in the Upper Midwest, United States. The results show that the PBL equilibrium features of δ m and δ P are well-constrained by thermal and moisture equilibrium in the PBL. For this study region, the summer values of rain evaporation fraction (f) and local water recycling ratio (LRR) are estimated to be 0.09 and 0.29 ± 0.12, respectively. The summer LRR values for the years 2006–2010 were 0.35, 0.36, 0.17, 0.29, and 0.29, respectively. The small value of LRR in 2008 corresponded to a drought condition with the lowest precipitation and second lowest ET among the five years. The summer magnitude of the amount effect is –2.8‰ (mm day –1 ) –1 and –0.8‰ (mm day –1 ) –1 for HDO and H 2 18 O, respectively. The local water recycling is identified as a significant factor influencing the continental isotope effect. Cropland has likely changed the regional LRR by –7.6 to 19.5% under different pre-agriculture land use scenarios. In conclusion, the feedback processes revealed here indicate that local water recycling is expected to be weakened under drought conditions, but it will be enhanced if irrigation is applied more intensely with more frequent drought events as the climate continues to warm.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Zero carbon solid-state rechargeable redox fuel for long duration and seasonal storage

Here, this work presents a unique thermochemical process for charging magnesium-manganese-oxide-based solid-state rechargeable redox fuel. High-temperature heating of the processing furnace can be driven by either renewable electricity or concentrated solar power. The simple recyclable fuel charging concept is based on a tubular falling bed reactor with countercurrent oxygen-depleted gas flow for complete heat recuperation. The main focus of this work is achieving solid flowability at high temperatures (1,450°C) and extracting chemically charged solid at ambient temperature with minimal energy loss. The operation strategies described in this work have enabled consistent magnesium-manganese-oxide particle flow up to 1,450°C. The measured extent of the thermal reduction reaction after cooling the particles is more than 90% of the fully reduced state at equilibrium. The thermal-to-chemical efficiency and overall system efficiency are 96% and 35% respectively, which are the highest reported for thermochemical fuels to date.

25 ENERGY STORAGE↗
DOE OSTI2024

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Physics↗
DOE OSTI2021

Simulations of multicomponent relativistic thermalization

Multicomponent relativistic fluids have been studied for decades. However, simulating the dynamics of the particles and fluids in such a mixture has been a challenge due to the fact that such simulations are computationally expensive in three spatial dimensions. Here, we report on the development and application of a multidimensional relativistic Monte Carlo code to explore the thermalization process in a relativistic multicomponent environment in a computationally inexpensive way. As an illustration we simulate the fully relativistic three-dimensional Brownian-motion-like solution to the thermalization of a high-mass particle (proton) in a bath of relativistic low-mass particles (electrons). Here, we follow the thermalization and ultimate equilibrium distribution of the Brownian-like particle as can happen in the cosmic plasma during big-bang nucleosynthesis. We also simulate the thermalization of energetic particles injected into the plasma as can occur, for example, by the decay of massive unstable particles during the big bang.

79 ASTRONOMY AND ASTROPHYSICS↗