Search NASA⌕ Search

SEARCH · Search NASA

Results for “thermodynamic 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An extension of the localized artificial diffusivity method for immiscible and high density ratio flows

The localized artificial diffusivity (LAD) method is widely regarded as the preferred multi-material regularization scheme for the compact finite difference method, because it is conservative, easy to implement, and generally robust for a wide range of multi-material problems. However, traditional LAD methods face significant challenges when applied to flows with large density ratios and when maintaining thermodynamic equilibrium across material interfaces. These limitations arise from the formulation of the artificial diffusivity flux and the reliance on enthalpy diffusion for interface regularization. Additionally, traditional LAD methods struggle to ensure stability under large density ratio conditions, fail to maintain a finite interface thickness, and are therefore unsuitable for modeling immiscible interfaces. Here, in this work, we discuss the origins of these issues in traditional LAD methods and propose modifications which enable the simulation of large density ratio and immiscible flows. The proposed method targets the artificial diffusion fluxes at gradients and ringing in the volume fraction, rather than the mass fraction in traditional methods, to consistently regularize large density ratio interfaces. Furthermore, the proposed method introduces an artificial bulk density diffusion term to enforce equilibrium conditions across interfaces. To address the challenge of modeling immiscible flows, a conservative diffuse interface term is incorporated into the formulation to ensure a finite interface thickness. Specific consideration is taken in the design of the method to ensure that these crucial properties are maintained for N -material flows. The effectiveness of the proposed method is demonstrated through a series of canonical test cases, and its accuracy is validated by comparison with experimental data on micro-bubble collapse in water. These results highlight the method’s robustness and its ability to overcome the limitations of traditional LAD approaches.

Artificial diffusivity↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
DOE OSTI2026

Dopant adsorption as a function of bulk concentration at a near 42º (100) twist grain boundary in SrTiO3

The enrichment of grain boundaries with dopant atoms is of critical importance for the macroscopic physical properties of materials. In thermodynamic equilibrium the Gibbs adsorption isotherm relates grain boundary excess of dopant atoms, their chemical potential in the adjacent bulk, and the respective interface energy. This study has used bicrystals with a near 42º (100) twist grain boundary in SrTiO3 to demonstrate that different kinetic pathways of Fe dopant atom additions converge towards comparable grain boundary configurations. Despite differences in bulk chemical potentials remarkably similar grain boundary excess quantities were observed. The experimental results indicate the general experimental feasibility to establish quantitative relationships between variations of grain boundary energy, interfacial excess, and overall dopant concentration. Improved experimental counting statistics are needed to distinguish solute interface excess as a function of bulk chemical potential.

Hahn, William [University of California, Davis]↗
DOE OSTI2022

Lamellae of α phase in a metastable β -Ti alloy studied by small-angle x-ray scattering

The shapes and orientations of lamellae of the hexagonal close-packed α phase in single crystals of the body-centred cubic β phase in a metastable β -Ti alloy (Timetal LCB) were investigated by synchrotron small-angle x-ray scattering. Several steps of thermal treatment with various annealing times and temperatures were employed. The measured scattering data were fitted to a model of ellipsoidal lamellae of random sizes, and the mean lamella shapes and orientations of their axes with respect to the β -Ti matrix were determined from the fit. The number and relative volume of the lamellae were estimated from the integrated scattered intensity. The elastic energy of the lamellae as well as the energy of defects at the α / β interface (misfit dislocations) were calculated from the parameters of the lamellae. In conclusion, the results provide information on the kinetics of the β $ \rightarrow$ α transformation in metastable β -Ti alloys and allow assessing the degree of thermodynamic equilibrium of the final α / β composite.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Real-time monitoring of the orthohydrogen fraction in a liquid hydrogen moderator

Here, we present the first direct measurements of the orthohydrogen fraction in an operating liquid hydrogen moderator loop at a high-powered neutron source. We have developed a Raman immersion probe capable of characterizing the orthohydrogen and parahydrogen makeup of a liquid hydrogen loop approximately every ten seconds, which can see orthohydrogen fractions lower than 0.20 without difficulty. At the Spallation Neutron Source, where the time-averaged power deposition in the studied moderator system is around 0.02 W/cm 3 , the orthohydrogen fraction increases with beam power and equilibrates at levels of 0.40 and responds to beam trips in a dynamic fashion. Additionally, we find that without deliberate catalyst systems in place, the relaxation due to the piping in the system has insignificant effect in catalyzing the relaxation of freshly condensed hydrogen toward thermodynamic equilibrium, and that the radiation levels associated with activation in a non-operating spallation neutron source target-moderator-reflector assembly likewise has insignificant effect on the liquid hydrogen kinetics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2021

Early- and late-time dynamics of laser-produced plasmas by combining emission and absorption spectroscopy

Optical emission spectroscopy and laser absorption spectroscopy are two highly complementary techniques for diagnosing plasmas, specifically in the determination of important parameters such as temperature and density. In this work, we characterized a laser-produced plasma by combining emission and tunable absorption spectroscopy. Plasmas are produced by focusing pulses from a frequency quadrupled Nd:YAG laser onto a glass target containing a minor concentration of aluminum in an ambient atmosphere of ~100 Torr N 2 . Here, emission spectroscopy is used to investigate the line behavior of Al and Ca, electron density, and excitation temperature at times = 5 μs, while absorption spectroscopy is used to investigate the line behavior of Al and Ca, species density, and kinetic temperature at times = 5-100 μs. A comparison is made between the excitation temperature measured from emission using the Boltzmann plot method and the kinetic temperature measured from absorption; the implications of this comparison on local thermodynamic equilibrium and line-of-sight averaging are discussed. By combining the two methods of optical emission and laser absorption spectroscopy, a comprehensive characterization of the plasma over the entire lifetime of its evolution is possible, leading to a greater understanding of the plasma behavior including particle condensation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Comparison of atomized mass and crater volume in laser ablation

The extent and dynamics of laser ablation are typically studied using crater imaging or by invoking the relationship between atomic emission and the mass removed. The former is a static view of a dynamic process and requires the accumulation of multiple shots in one location. The latter is complex and not absolute without calibration of the optical system with a standard of spectral radiance. We measure the mass of the atomized plume by laser atomic-absorption spectroscopy at 2 µs without external mass calibration; mass uncertainty in laser-ablation atomic absorption spectroscopy (LA-AAS) results from uncertainties in line area fitting, tabulated oscillator strengths, and partition functions as well as a <10% underrepresentation of the mass due to treating the plasma as a single thermodynamic equilibrium. The LA-AAS-measured mass pertains to the atoms in the probed charge states and does not include condensed or molecular species. The ablation efficiency at 300 mbar of helium varies significantly between the two focusing conditions used. A comparison of crater-derived masses and LA-AAS masses suggests that more defocused ablation may result in significant redeposition of material in the crater, distorting the conclusions from drilling studies and crater imaging. More focused ablation conditions result in crater-derived masses that exceed those measured with LA-AAS and suggest melt expulsion or phase explosion.

Merten, Jonathan [Arkansas State University, Jones↗
DOE OSTI2022

Leveraging computational thermodynamics to guide SiC-ZrC chemical vapor deposition process development

Here using the CALPHAD approach to understand zirconium carbide deposition, a series of phase equilibria were calculated from a custom thermodynamic database based on a literature source, and the equilibria were used to explore the potential chemical vapor deposition (CVD) processing space in the ZrCl 4 -CH 3 SiCl 3 -CH 4 -H 2 system as a function of pressure, temperature, and gas composition. Several gas ratios were considered. At a given ZrCl 4 :CH 3 SiCl 3 ratio within the range studied, the most important factor was found to be the ratios of CH 4 :ZrCl 4 , wherein the nature of the composition – carbide vs. silicide – could be controlled. A pure binary composition of ZrC and SiC is expected to form by increasing the initial amount of methane and decreasing the amount of hydrogen from values predicted purely based on thermodynamic equilibrium. Rietveld analysis of the x-ray diffractograms from corresponding experimental depositions confirmed that increasing the CH 4 :ZrCl 4 ratio increased the fraction of carbon-containing species (SiC, ZrC) and decreased the fraction of non-carbides (ZrSi, ZrSi 2 , etc.), as predicted from the CALPHAD results.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

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

Macromolecular Crowding Affects Voltage-Dependent Alamethicin Pore Formation in Lipid Bilayer Membranes

Macromolecular crowding is known to modulate chemical equilibria, reaction rates, and molecular binding events, both in aqueous solutions and at lipid bilayer membranes, natural barriers that enclose the crowded environments of cells and their subcellular compartments. Previous studies on the effects of macromolecular crowding in aqueous compartments on conduction through membranes have focused on single-channel ionic conduction through previously formed pores at thermodynamic equilibrium. Furthermore, the effects of macromolecular crowding on the mechanism of pore formation itself were studied using the droplet interface bilayer (DIB) technique with the voltage-dependent pore-forming peptide alamethicin (alm). Macromolecular crowding was varied using 8 kDa molecular weight polyethylene glycol (PEG8k) or 500 kDa dextran (DEX500k) in two aqueous droplets on both sides of the bilayer membrane. In general, voltage thresholds for pore formation in the presence of crowders in the droplets decreased compared to their values in the absence of crowders, due to excluded volume effects, water binding by PEG, and changes in the ordering of water molecules and hydrogen-bonding interactions involving the polar lipid headgroups. In addition, asymmetric crowder loading (e.g., PEG8k–DEX500k on either side of the membrane) resulted in transmembrane osmotic pressure gradients that either enhanced or degraded the ionic conduction through the pores.

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

Experiments on the Capillary Condensation/Evaporation Hysteresis of Pure Fluids and Binary Mixtures in Cylindrical Nanopores

Here, a new experimental method, i.e., isochoric cooling/heating cycle measurement using differential scanning calorimeter (DSC), is applied to investigate the hysteretic phenomenon of pure components and mixtures confined in two different types of cylindrical nanopores (SBA-15 and KIT-6). The capillary condensation/evaporation hysteresis end point (T he ) for a specific system is determined by performing a series of measurements with different densities by increasing initial temperature and pressure. We confirm that for a pure component, T he is lower when confined in smaller nanopores, whereas T he for CH 4 /CO 2 and CH 4 /C 2 H 6 mixtures at the specific composition investigated are higher than that of pure CO 2 and C 2 H 6 , respectively, when confined in the same nanopores. To the best of our knowledge, this is the first time that T he of fluid mixture is experimentally studied, which further demonstrates that binary mixtures behave similarly to pure fluids when confined in nanoporous media. At low temperatures where hysteresis is observed, the shape as well as the duration of the exothermic peak (due to capillary condensation) and endothermic peak (due to capillary evaporation) indicate that capillary condensation, rather than capillary evaporation, takes place at/near a thermodynamic equilibrium. Within the accuracy of the measurements, we find that the total heat of capillary condensation is identical to that of capillary evaporation for all systems examined, regardless of the presence of hysteresis. By comparing the heat per unit volume among different systems, we also find that the heat per unit volume of the phase transition of confined pure fluid is lower than that of bulk, which further demonstrates that the molar density of confined pure liquid is less than that of bulk pure liquid at saturation.

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

Synthesis, Magnetic Properties, and Electronic Structure of Magnetic Topological Insulator MnBi 2 Se 4

The intrinsic magnetic topological insulators MnBi 2 Te 4 and MnBi 2 Se 4 support novel topological states related to symmetry breaking by magnetic order. Unlike MnBi 2 Te 4 , the study of MnBi 2 Se 4 has been inhibited by the lack of bulk crystals, as the van der Waals (vdW) crystal is not the thermodynamic equilibrium phase. In this study, we report the layer-by-layer synthesis of vdW MnBi 2 Se 4 crystals using nonequilibrium molecular beam epitaxy. Atomic-resolution scanning transmission electron microscopy and scanning tunneling microscopy identify a well-ordered vdW crystal with septuple-layer base units. The magnetic properties agree with the predicted layered antiferromagnetic ordering but disagree with its predicted out-of-plane orientation. Instead, our samples exhibit an easy-plane anisotropy, which is explained by including dipole–dipole interactions. Angle-resolved photoemission spectroscopy reveals the gapless Dirac-like surface state, which demonstrates that MnBi 2 Se 4 is a topological insulator above the magnetic-ordering temperature. These studies show that MnBi 2 Se 4 is a promising candidate for exploring rich topological phases of layered antiferromagnetic topological insulators.

36 MATERIALS SCIENCE↗