Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equilibrium modeling”

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

Low-Temperature Plasma Assisted Kinetics Study of Ethanol

Non-equilibrium plasma assisted kinetics of ethanol at 0.5 atm over a temperature range of 523-1203 K were studied at stoichiometric conditions by using a plasma flow reactor (PFR). The effects of plasma chemistry on fuel and its influence on neutral chemistry is decoupled with both experimental and modelling efforts. Enhanced reactivity was observed experimentally in both plasma-assisted oxidation and pyrolysis cases with fuel consumption beginning at temperatures as low as 523 K. Absolute fuel consumption was observed at temperatures nearly 200 K earlier than pure thermal reactions highlighting the effectiveness of non-equilibrium plasma to promote reactivity. Furthermore, modelling results illustrated the effect of electron impact reactions towards accelerating fuel consumption and promoting low temperature chemistry by generating excited state N2 and reactive O and OH radicals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A dissipation-induced superradiant transition in a strontium cavity-QED system

Driven-dissipative many-body systems are ubiquitous in nature and a fundamental resource for quantum technologies. However, they are also complex and hard to model because they cannot be described by the standard tools in equilibrium statistical mechanics. Probing nonequilibrium critical phenomena in pristine setups can illuminate fresh perspectives on these systems. Here, we use an ensemble of cold 88 Sr atoms coupled to a driven high-finesse cavity to study the cooperative resonance fluorescence (CRF) model, a classic driven-dissipative model describing coherently driven dipoles superradiantly emitting light. We observe its nonequilibrium phase diagram characterized by a second-order phase transition. Below a critical drive strength, the atoms quickly reach the so-called superradiant steady state featuring a macroscopic dipole moment; above the critical point, the atoms undergo persistent Rabi-like oscillations. At longer times, spontaneous emission transforms the second-order transition into a discontinuous first-order transition. Our observations pave the way for harnessing robust entangled states and exploring boundary time crystals in driven-dissipative systems.

Science & Technology - Other Topics↗
DOE OSTI2022

On the effect of mixing-driven vaporization in a homogeneous relaxation modeling framework

The homogeneous relaxation model (HRM) is one of the most widely used models to describe the liquid–gas phase transition in multiphase flows due to the occurrence of cavitation. However, in its original formulation, the HRM does not account for the presence of ambient gas species, which generally limits its applicability to the injector's internal flow where ambient gases are negligible. In this work, a mixing-driven vaporization (MDV) model was developed to extend the capability of the HRM in handling the mixing effect in the regions external to the nozzle, where vapor–liquid equilibrium for multi-species mixtures of fuel and ambient gas is considered. Herein, to assess the model performance, simulations of the Engine Combustion Network's Spray G injector were performed with the HRM and the MDV model under both flash-boiling and evaporating conditions. It was found that the MDV model led to a better match against x-ray measurements of fuel density in the near-nozzle region. In contrast to the HRM, the MDV model was able to reproduce the vaporization process in the mixing zone at the edge of the fuel jet, which aligns with the expected physics. This resulted in substantial differences in the prediction of other flow characteristics such as mixture temperature and pressure. Furthermore, this work demonstrates that evaporation timescales have a considerable effect on the MDV model's predictions, as shown by a parametric study in which a time factor was introduced to mimic the effect of different timescales due to different phase change mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Non-Equilibrium Actinide Radiation Chemistry and the Nuclear Fuel Cycle

Actinides are inherently unstable elements that frequently coexist with other radioisotopes, generating intense ionizing radiation fields that drive the formation of non equilibrium oxidation states. These transient species exert a profound mechanistic influence on the radiation response of actinide containing systems due to their unique redox chemistry. Despite their importance, they remain poorly understood, yet such insight is essential for advancing actinide science and accurately predicting radiation driven behavior. Actinide separations—critical for nuclear energy technologies, strategic deterrence, space exploration, and nuclear medicine—depend on precise control of actinide oxidation states to recover targeted elements from complex matrices such as used nuclear fuel. However, during these processes, actinides, their coordination complexes, and the separation media are all exposed to intense, multicomponent (alpha, beta, gamma, etc.) radiation fields that can alter process efficiency, selectivity, and chemical stability. Understanding, controlling, and mitigating radiation induced reactions is therefore key to innovating and optimizing next generation separation technologies. This seminar will provide an overview of the nuclear fuel cycle and non equilibrium actinide radiation chemistry in the context of recovering actinides from used nuclear fuel, with a particular emphasis on direct dissolution–based reprocessing strategies. We will explore time resolved electron pulse radiolysis and alpha and gamma dose accumulation studies, integrated with multiscale computational modeling, to elucidate the molecular level roles of radiation driven, non equilibrium actinide species in process performance and in the radiolytic stability of organic ligands used for actinide recovery. These insights offer new pathways for designing advanced separation methods and next generation solvent systems, with broad implications for the future of the nuclear fuel cycle.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2020

Thermodynamic Origin of Reaction Non-Uniformity in Battery Porous Electrodes and Its Mitigation

The development of non-uniform reaction current distribution within porous electrodes is a ubiquitous phenomenon during battery charging/discharging and frequently controls the rate performance of battery cells. Reaction inhomogeneity in porous electrodes is usually attributed to the kinetic limitation of mass transport within the electrolyte and/or solid electrode phase. In this work, however, we reveal that it is also strongly influenced by the intrinsic thermodynamic behavior of electrode materials, specifically the dependence of the equilibrium potential on the state of charge: the electrode reaction becomes increasingly non-uniform when the slope of the equilibrium potential curve is reduced. We employ numerical simulations and equivalent circuit model to elucidate such a correlation and show that the degree of reaction inhomogeneity and the resultant discharge capacity can be predicted by a dimensionless reaction uniformity number. For electrode materials that have equilibrium potentials insensitive to the state of charge and exhibit significant reaction non-uniformity, we demonstrate several approaches to spatially homogenizing the reaction current inside porous electrodes, including matching the electronic and ionic resistances, introducing graded electronic conductivity and reducing the surface reaction kinetics.

25 ENERGY STORAGE↗
DOE OSTI2021

Low-temperature domain-wall freezing and nonequilibrium dynamics in the transverse-field Ising model material CoNb 2 O 6

CoNb 2 O 6 is a rare realization of the transverse-field Ising model, making it a useful tool for studying both equilibrium and nonequilibrium many-body quantum physics. Despite a large body of work dedicated to characterizing this material, details of the ordered states in the presence of relatively weak transverse fields have not been discussed in detail. Here, we present a detailed study of CoNb 2 O 6 via ac susceptibility measurements in order to further characterize its low-temperature behavior in the presence of a transverse field. Specifically, we call attention to an unconventional freezing transition in zero field below T F = 1.2 K, occurring within the well-known commensurate antiferromagnetic (CAFM) state that onsets at T N2 = 1.9 K. We performed a series of transverse-field quenches into this frozen state, which resulted in a slowly relaxing susceptibility, χ' (t), that followed a logarithmic decay within the time range measured. Additionally, we discuss the frozen state in the context of the freezing of previously discussed “free” chains arising from domain walls between the four degenerate sublattices of the CAFM state. We also attempted to observe Kibble-Zurek scaling by quenching the transverse field into the frozen state at different rates. This produced a null result; the behavior can be fully explained by coarsening of domains over the time scale of the quenches. The absence of a clear Kibble-Zurek scaling is itself surprising, given the proposed ubiquity of the phenomenon for general second-order phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model

This report details a fully coupled neutronics thermal hydraulics reference plant model for a gas-cooled high-temperature pebble-bed reactor. The multiphysics model is developed on the Nuclear Regulatory Commission's Comprehensive Reactor Analysis Bundle (BlueCRAB) available on the Idaho National Laboratory high performance computer, which natively and seamlessly couples Griffin, Pronghorn, and the BISON MOOSE-based applications. Griffin provides the reactor physics capabilities including depletion to the equilibrium core, k-eigenvalue, adjoint, and transient solutions. Pronghorn solves the porous medium equations for the fluid regions and conduction in the solid regions. BISON solves thermal conduction problems for the average pebble and TRISO in the pebble-bed core thus providing the fuel and moderator spatial fields for each pebble burnup group. The neutronics feedback relies primarily on fuel and moderator temperatures. We present results for the coupled steady state equilibrium core and a protected loss of flow event. Although this model is to be considered prototypical with regard to the capabilities in BlueCRAB, its results are consistent with published work by the Institute of Nuclear and New Energy Technology (INET) in China and other research entities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Consensus Equilibrium Approach for 3-D Land Seismic Shots Recovery

Physical and budget constraints often result in inadequate sampling for accurate subsurface imaging. Preprocessing approaches, such as missing trace interpolation, are typically employed to enhance seismic data in such cases. The compressed sensing (CS) framework has been applied for modeling missing seismic data, which is estimated by sparsity-based computational algorithms. While existing work mainly focuses on recovering missing traces resulting from receiver subsampling, source subsampling has greater economical advantages, as sources are more expensive than receivers. Moreover, stronger image models different from sparsity have not been explored for source recovery. This work presents a consensus equilibrium (CE) approach to recover missing seismic shots, which enables to incorporate various regularization operators modeling different data priors. Here, simulation results from a real 3-D land seismic dataset demonstrate that the CE approach provides more accurate estimations of the linear and hyperbolic events in the recovered shots, compared with pure sparsity-based reconstructions.

58 GEOSCIENCES↗
DOE OSTI2021

Role of Equilibrium Fluctuations in Light-Induced Order

Engineering novel states of matter with light is at the forefront of materials research. An intensely studied direction is to realize broken-symmetry phases that are “hidden” under equilibrium conditions but can be unleashed by an ultrashort laser pulse. Despite a plethora of experimental discoveries, the nature of these orders and how they transiently appear remain unclear. To this end, we investigate a nonequilibrium charge density wave (CDW) in rare-earth tritellurides, which is suppressed in equilibrium but emerges after photoexcitation. Using a pump-pump-probe protocol implemented in ultrafast electron diffraction, we demonstrate that the light-induced CDW consists solely of order parameter fluctuations, which bear striking similarities to critical fluctuations in equilibrium despite differences in the length scale. By calculating the dynamics of CDW fluctuations in a nonperturbative model, we further show that the strength of the light-induced order is governed by the amplitude of equilibrium fluctuations. Overall, these findings highlight photoinduced fluctuations as an important ingredient for the emergence of transient orders out of equilibrium. Our results further suggest that materials with strong fluctuations in equilibrium are promising platforms to host hidden orders after laser excitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Thermodynamic Modeling of Point Defects in Triple Conducting Perovskite Ba 0.95 La 0.05 FeO 3- $_δ$ with Incorporation of the Hydride Defect Formation Reaction for Solid Oxide Cells

Distinct from the proton defect, the hydride defect species may be present in certain perovskite materials in reducing environments such as in fuel electrodes of solid oxide cells (SOCs) or on the reducing side of ceramic membranes. A generalized defect thermodynamic model was developed for the triple-conducting perovskites (La,Ba)Fe 1-x M x O 3-δ (M = Y and Zr) to allow inclusion of the hydride defect formation reaction in addition to the other three main defect reactions, namely, the oxygen vacancy formation, hydration, and charge disproportionation reactions. This comprehensive defect model also allows the incorporation of polynomial functional forms of oxygen nonstoichiometry δ to describe the defect reaction energies and entropies and to enable refinements of the defect reaction equilibrium constants in the defect thermodynamic analysis. As a first step, the developed model is applied to the Ba 0.95 La 0.05 FeO 3-δ material as an illustrative system to obtain its Brouwer diagrams with both the proton and hydride defects in relevant SOC conditions, particularly for more reducing environments. In conclusion, the results provide direct guidance on the influence of electronic and ionic defect concentrations upon thermodynamic properties and ultimately on the performance of Ba 0.95 La 0.05 FeO 3-δ and potentially other (La,Ba)Fe 1-x M x O 3-δ perovskite materials involved in SOC applications.

25 ENERGY STORAGE↗
DOE OSTI2021

High‐Efficiency Ion‐Exchange Doping of Conducting Polymers

Abstract Molecular doping—the use of redox‐active small molecules as dopants for organic semiconductors—has seen a surge in research interest driven by emerging applications in sensing, bioelectronics, and thermoelectrics. However, molecular doping carries with it several intrinsic problems stemming directly from the redox‐active character of these materials. A recent breakthrough was a doping technique based on ion‐exchange, which separates the redox and charge compensation steps of the doping process. Here, the equilibrium and kinetics of ion exchange doping in a model system, poly(2,5‐bis(3‐alkylthiophen‐2‐yl)thieno(3,2‐b)thiophene) (PBTTT) doped with FeCl 3 and an ionic liquid, is studied, reaching conductivities in excess of 1000 S cm −1 and ion exchange efficiencies above 99%. Several factors that enable such high performance, including the choice of acetonitrile as the doping solvent, which largely eliminates electrolyte association effects and dramatically increases the doping strength of FeCl 3 , are demonstrated. In this high ion exchange efficiency regime, a simple connection between electrochemical doping and ion exchange is illustrated, and it is shown that the performance and stability of highly doped PBTTT is ultimately limited by intrinsically poor stability at high redox potential.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C 2 H 2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. Furthermore, the rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH–PAH stacking with size.

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

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Edge spin transport in the disordered two-dimensional topological insulator $\mathrm{WTe_2}$

The spin conductance of two-dimensional topological insulators (2D TIs) is not expected to be quantized in the presence of perturbations that break the spin-rotational symmetry. However, the deviation from the pristine-limit quantization has yet to be studied in detail. In this paper, we define the spin current operator for the helical edge modes of a 2D TI and introduce a four-terminal setup to measure spin conductances. Using the developed formalism, we consider the effects of disorder terms that break spin-rotational symmetry or give rise to edge-to-edge coupling. We identify a key role played by spin torque in an out-of-equilibrium edge. We then utilize a tight-binding model of topological monolayer $\mathrm{WTe_2}$ and scattering matrix formalism to numerically study spin transport in a four-terminal 2D TI device. In particular, we calculate the spin conductances and characteristic spin decay length in the presence of magnetic disorder. In addition, we also study the effects of interedge scattering in a quantum point contact geometry. Here we find that the spin Hall conductance is surprisingly robust to spin symmetry-breaking perturbations, as long as time-reversal symmetry is preserved and interedge scattering is weak.

36 MATERIALS SCIENCE↗