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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

How does humidity data impact land surface modeling of hydrothermal regimes at a permafrost site in Utqia$\mathrm{\dot{g}}$vik, Alaska?

Humidity is a basic and crucial meteorological indicator commonly measured in several forms, including specific humidity, relative humidity, and absolute humidity. These different forms can be inter-derived based on the saturation vapor pressure (SVP). In past decades, dozens of formulae have been developed to calculate the SVP with respect to, and in equilibrium with, liquid water and solid ice surfaces, but many prior studies use a single function for all temperature ranges, without considering the distinction between over the liquid water and ice surfaces. These different approaches can result in humidity estimates that may impact our understanding of surface-subsurface thermal-hydrological dynamics in cold regions. In this study, we compared the relative humidity (RH) downloaded and calculated from four data sources in Alaska based on five commonly used SVP formulas. These RHs, along with other meteorological indicators, were then used to drive physics-rich land surface models at a permafrost-affected site. Here, we found that higher values of RH (up to 40%) were obtained if the SVP was calculated with the over-ice formulation when air temperatures were below freezing, which could lead to a 30% maximum difference in snow depths. The choice of whether to separately calculate the SVP over an ice surface in winter also produced a significant range (up to 0.2 m) in simulated annual maximum thaw depths. The sensitivity of seasonal thaw depth to the formulation of SVP increases with the rainfall rate and the height of above-ground ponded water, while it diminishes with warmer air temperatures. These results show that RH variations based on the calculation of SVP with or without over-ice calculation meaningfully impact physically-based predictions of snow depth, sublimation, soil temperature, and active layer thickness. Under particular conditions, when severe flooding (inundation) and cool air temperatures are present, care should be taken to evaluate how humidity data is estimated for land surface and earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Ion–ion interactions enhance aluminum solubility in alkaline suspensions of nano-gibbsite (α-Al(OH) 3 ) with sodium nitrite/nitrate

Despite widespread industrial importance, predicting metal solubilities in highly concentrated, multicomponent aqueous solutions is difficult due to poorly understood ion–ion and ion–solvent interactions. Aluminum hydroxide solid phase solubility in concentrated sodium hydroxide (NaOH) solutions is one such case, with major implications for ore refining, as well as processing of radioactive waste stored at U.S. Department of Energy legacy sites, such as the Hanford Site, Washington State. The solubility of gibbsite (α-Al(OH) 3 ) is often not well predicted because other ions affect the activity of hydroxide (OH – ) and aluminate (Al(OH) 4 – ) anions. In the study reported herein, we systematically examined the influence of key anions, nitrite (NO 2 – ) and nitrate (NO 3 – ), as sodium salts on the solubility of α-Al(OH) 3 in NaOH solutions taking care to establish equilibrium from both under- and oversaturation. Rapid equilibration was enabled by use of a highly pure and crystalline synthetic nano-gibbsite of well-defined particle size and shape. Measured dissolved aluminum concentrations were compared with those predicted by an α-Al(OH) 3 solubility model derived for simple Al(OH) 4 – /OH – systems. Specific anion effects were expressed as an enhancement factor (Al enhc ) conveying the excess of dissolved aluminum. At 45 °C, NaNO 2 and NaNO 3 -containing systems exhibited Al enhc values of 2.70 and 1.88, respectively, indicating significant enhancement. The solutions were examined by Raman and high-field 27 Al NMR spectroscopy, indicating specific interactions including Al(OH) 4 – –Na + contact ion pairing and Al(OH) 4 – –NO 2 – /NO 3 – ion–ion interactions. Dynamic evolution of the α-Al(OH) 3 particles including growth and agglomeration was observed revealing the importance of dissolution/reprecipitation in establishing equilibrium. These studies indicate that incomplete ion hydration, as a result of the low water activity in these concentrated electrolytes, results in: (i) enhanced reactivity of the hydroxide ion with respect to α-Al(OH) 3 ; (ii) increased concentrations of Al(OH) 4 – in solution; and (iii) stronger ion–ion interactions that act to stabilize the supersaturated solutions. This information on the mechanisms by which α-Al(OH) 3 becomes supersaturated is essential for more energy-efficient aluminum processing technologies, including the treatment of millions of gallons of Al(OH) 4 – -rich high-level radioactive waste.

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

Positional ordering induced by dynamic steric interactions in superparamagnetic rods

The dynamic motion produced by precessing magnetic fields can drive matter into far-from-equilibrium states. We predict 1D periodic ordering in systems of precessing rods when magnetic interactions between rods remain insignificant. The precession angle of the rods is completely determined by the field's precession angle and the ratio of the field's precession frequency and the characteristic response frequency of the rods. Here, we develop a molecular dynamics model that explicitly calculates magnetic interactions between particles, and we also simulate rods in the limit of a strong and fast precessing magnetic field where inter-rod magnetic interactions are negligible, using a purely steric model. Our simulations show how steric interactions drive the rods from a positionally disordered phase (nematic) to a layered (smectic) phase. As the rod precession angle increases, the nematic–smectic transition density significantly decreases. The minimization of unfavorable steric interactions also induces phase separation in binary mixtures of rods of different lengths. This effect is general to any force that produces precession in elongated particles. This work will advance the understanding and control of out-of-equilibrium soft matter systems.

42 ENGINEERING↗
DOE OSTI2022

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

DNA Reaction–Diffusion Attractor Patterns

Living systems can form and recover complex chemical patterns with precisely sized features in the ranges of tens or hundreds of microns. We show how designed reaction-diffusion processes can likewise produce precise patterns, termed attractor patterns, that reform their precise shape after being perturbed. We use oligonucleotide reaction networks, photolithography and microfluidic delivery to form precisely controlled attractor patterns and study the responses of these patterns to different localized perturbations. Linear and ‘hill’-shaped patterns formed and stabilized into shapes and at time scales consistent with reaction-diffusion models. When patterns were perturbed in particular locations with UV light, they reliably reformed their steady state profiles. Recovery also occurred after repeated perturbations. As a result, by designing the far-from-equilibrium dynamics of a chemical system, this study shows how it is possible to design spatial patterns of molecules that are sustained and regenerated by continually evolving towards a specific steady state configuration

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

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effects of neoclassical tearing modes and toroidal field ripple on lost alpha power in the SPARC tokamak

Here, using the SPIRAL Monte Carlo, full particle-orbit simulation code, we investigate the effects of neoclassical tearing modes (NTMs) and toroidal field (TF) ripple on alpha power losses during steady-state operation of the SPARC primary reference discharge. Model perturbations for TF ripple and the $m/n = 2/1$ and $3/2$ NTMs with exaggerated widths selected based on an H-mode plasma approaching thermal quench are added to a simulated SPARC magnetic equilibrium through which marker particles are tracked. The 3/2 and 2/1 NTMs are located at $\rho_{\mathrm{pol}} \sim 0.76$ and $\rho_\mathrm{pol} \sim 0.86$ respectively, well positioned to increase alpha particle transport into and within an outer lossy region of the plasma beyond $\rho_\mathrm{pol} \sim 0.8$ where over 95% of lost alpha particles are born. Total alpha power losses are shown to increase modestly from 1.73% lost at a minimum to 2.34% lost at a maximum, and alpha particle surface power densities form localized hotspots on the first-wall near the low-field side midplane due to NTMs and TF ripple. We establish a conservative upper limit for first-wall alpha surface power densities on a toroidally symmetric wall for typical, flattop operation and motivate the consideration of NTMs in the design of three dimensional limiter surfaces for SPARC.

alpha power losses↗
DOE OSTI2023

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Kinetic Monte Carlo Tool for Kinetic Modeling of Linear Step-Growth Polymerization: Insight into Recycling of Polyurethanes

A kinetic Monte Carlo model of polyurethane polymerization which explicitly tracks the polymer sequences is developed and shared. Here this model is benchmarked against theoretical and experimental polyurethane data and used to investigate the effect on oligomer distributions of unequal reactivity of the first and second isocyanate to react. The reverse reactions using thermodynamic consistency are then added to the framework, and analogous to the addition polymerization concept of ceiling temperature, equilibrium chain length distributions at various temperatures are calculated. For a mixture of three monomers AA, BB, and CC, where BB and CC do not react with one another, are present in stoichiometric proportions, and have different enthalpies of reaction with AA, an odd-even effect emerges. Odd length chains are more likely than even length chains for temperatures at which BB and CC have significantly different equilibrium conversions. The concept of ceiling temperature that is typically cited for addition polymers is extended here to provide a measure of conditions under which depolymerization for recycling is favored.

condensation polymerization↗
DOE OSTI2025

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗
DOE OSTI2025

Alkali Metal Control of Triplet-Mediated C–H Activation in Iron-Mediated Coupling of Dinitrogen and Benzene

Direct coupling of N 2 with abundant feedstocks like benzene to form N-containing organic compounds is a promising strategy for N 2 fixation pathways. The challenge of coupling N 2 activation and C–H bond oxidative addition was recently solved by introducing a reversible benzene C–H bond activation process mediated by a low-valent Fe(0) complex, which gave an organometallic product that could couple with partially reduced N 2 . Interestingly, the energetics of the C–H oxidative addition/reductive elimination step depends on the choice of alkali metal. However, the reason why the alkali metal influences the C–H bond activation remained elusive. In this work, we present a comprehensive study on this Fe(0)-mediated reversible C–H activation. Through density functional theory combined with high-level coupled cluster calculations, we discovered that the intermediate-spin triplet (S = 1) controls the energy of the transition state for C–H cleavage, while the high-spin quintet (S = 2) controls the position of the equilibrium. Na + drives the equilibrium toward oxidative addition due to an electrostatic effect, while K + and Rb + are dominated by a steric effect that favors the iron(0) species. The key role played by nonbonding interactions in the Fe-mediated C–H activation provides a conceptual model for alkali control over organometallic transformations.

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

Lightwave terahertz quantum manipulation of nonequilibrium superconductor phases and their collective modes

In this work, we present a gauge-invariant density matrix description of nonequilibrium superconductor (SC) states with spatial and temporal correlations driven by intense terahertz (THz) lightwaves. We derive superconductor Bloch-Maxwell equations of motion that extend Anderson pseudospin models to include the Cooper pair center-of-mass motion and electromagnetic propagation effects. We thus describe quantum control of dynamical phases, collective modes, quasiparticle coherence, and high nonlinearities during cycles of carrier wave oscillations, which relates to our recent experiments. Coherent photogeneration of a nonlinear supercurrent with a dc component, achieved via condensate acceleration by an effective lightwave field, dynamically breaks the equilibrium inversion symmetry. Experimental signatures include high harmonic light emission at equilibrium-symmetry-forbidden frequencies, Rabi-Higgs collective modes and quasiparticle coherence, and nonequilibrium moving condensate states tuned by few-cycle THz fields. We use such lightwaves as an oscillating accelerating force that drives strong nonlinearities and anisotropic quasiparticle populations to control and amplify different classes of collective modes, e.g., damped oscillations, persistent oscillations, and overdamped dynamics via Rabi flopping. Recent phase-coherent nonlinear spectroscopy experiments can be modeled by solving the full nonlinear quantum dynamics including self-consistent light-matter coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

The Community Earth System Model version 2 (CESM2)

An overview of the Community Earth System Model version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long pre-industrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low-top” (40 km, with limited chemistry) and “high-top” (140 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low latitude precipitation and short-wave cloud forcing biases; better representation of the Madden-Julian Oscillation; better El Niño – Southern Oscillation-related teleconnections; and a global land carbon accumulation trend that agrees well with observationally-based estimates. Most tropospheric and surface features of the low- and high-top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1-5.3°C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2’s transient climate response of 1.9-2.0°C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and theyrepresent CESM2’s contributions to the Coupled Model Intercomparison Project phase 6 (CMIP6).

Danabasoglu, Gokhan↗
DOE OSTI2025

On High-Order/Low-Order and Micro-Macro Methods for Implicit Time-Stepping of the BGK Model

In this paper, a high-order/low-order (HOLO) method is combined with a micro-macro (MM) decomposition to accelerate iterative solvers in fully implicit time-stepping of the Bhatnagar–Gross–Krook (BGK) equation for gas dynamics. The MM formulation represents a kinetic distribution as the sum of a local Maxwellian and a perturbation. In highly collisional regimes, the perturbation away from initial and boundary layers is small and can be compressed to reduce the overall storage cost of the distribution. The convergence behavior of the MM methods, the usual HOLO method, and the standard source iteration method is analyzed on a linear BGK model. Both the HOLO and MM methods are implemented using a discontinuous Galerkin (DG) discretization in phase space, which naturally preserves the consistency between high- and low-order models required by the HOLO approach. Furthermore, the accuracy and performance of these methods are compared on the Sod shock tube problem and a sudden wall heating boundary layer problem. Overall, the results demonstrate the robustness of the MM and HOLO approaches and illustrate the compression benefits enabled by the MM formulation when the kinetic distribution is near equilibrium.

BGK model↗
DOE OSTI2023

Atomistic modeling of metal–nonmetal interphase boundary diffusion

Atomistic computer simulations are applied to investigate the atomic structure, thermal stability, and diffusion processes in Al–Si interphase boundaries as a prototype of metal–ceramic interfaces in composite materials. Some of the most stable orientation relationships between the phases found in this work were previously observed in epitaxy experiments. Here, a non-equilibrium interface can transform to a more stable state by a mechanism that we call interface-induced recrystallization. Diffusion of both Al and Si atoms in stable Al–Si interfaces is surprisingly slow compared with diffusion of both elements in Al grain boundaries but can be accelerated in the presence of interface disconnections. A qualitative explanation of the sluggish interphase boundary diffusion is proposed. Atomic mechanisms of interphase boundary diffusion are similar to those in metallic grain boundaries and are dominated by correlated atomic rearrangements in the form of strings and rings of collectively moving atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of the Anions on the Interaction Energy between Water and Ionic Liquids

It is essential to understand the interaction energy between water and ionic liquids (ILs) in search of ILs for dehumidification applications. In this paper, the intermolecular interaction energy between water and 1-ethyl-3-methylimidazolium ([EMIM])-based ILs was investigated. The vapor-liquid equilibrium data were first measured for solutions of [EMIM][BF 4 ], [EMIM][N(CN) 2 ], [EMIM][CF 3 SO 3 ], and [EMIM][CH 3 SO 3 ] with water molar fractions of 8–95 % at temperatures of 303–343 K. The IL-IL and water-IL intermolecular interaction energies were determined by using the correlated non-random two-liquid model. The results revealed that fluorine atoms in anions reduce the attraction energy between IL and water, and thus, decrease the water affinity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗