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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations

There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Computational investigation of water glasses using machine-learning potentials

The molecular origins of water’s anomalous properties have long been a subject of scientific inquiry. The liquid–liquid phase transition hypothesis, which posits the existence of distinct low-density and high-density liquid states separated by a first-order phase transition terminating at a critical point, has gained increasing experimental and computational support and offers a thermodynamically consistent framework for many of water’s anomalies. However, experimental challenges in avoiding crystallization near the postulated liquid–liquid critical point have focused attention to water’s canonical glassy states: low-density and high-density amorphous ice. Here, we use two Deep Potential machine-learning models, trained on the Strongly Constrained and Appropriately Normed density functional and the highly accurate Many-Body Polarizable potential, to conduct an investigation of water’s glassy phenomenology based on quantum mechanical calculations. Despite not being explicitly trained on amorphous ices, both models accurately capture the structure and transformation of the water glasses, including their interconversion along different thermodynamic paths. Isobaric quenching of liquid water at various pressures generates a continuum of intermediate amorphous ices and density fluctuations increase near the liquid–liquid critical pressure. The glass transition temperatures of the amorphous ices produced at different pressures exhibit two distinct branches, corresponding to low-density and high-density amorphous ice behaviors, consistent with experiment and the liquid–liquid transition hypothesis. Extrapolating transformation pressures from isothermal compressions to experimental compression rates brings our simulations into excellent agreement with data. Our findings demonstrate that machine-learning potentials trained on equilibrium phases can effectively model nonequilibrium glassy behavior and pave the way for studying long-timescale, out-of-equilibrium processes with quantum mechanical accuracy.

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

Validation of ablation model for polyethylene using pulsed x-ray and proton exposures

The surface erosion of polyethylene is interrogated using pulsed x rays at the Z Machine (Sandia National Laboratories) and with proton beams at the Gamble II generator (Naval Research Laboratory) to validate a coupled model for volumetric thermal ablation, photoionization, finite-rate decomposition, and molecular recombination of radicals. The intense radiation pulses (up to ∼1014W/m2 over tens of nanoseconds) are used to generate one-dimensional vapor flows with low ionization fractions and a simplified geometry compared to typical laser ablation, allowing for evaluation of the model under local thermal equilibrium conditions. Areal momentum carried by the ensuing uniaxial hydrodynamic shock is used to indicate the extent of ablation. The threshold fluence for ablation is found to be in close correspondence with the bulk melt transition, and reasonable agreement with the model is obtained for peak temperatures in polyethylene up to 5500 K and heating rates up to 1011K/s where thermal decomposition reactions are also active.

Moore, N. W. (ORCID:0000000272259674)↗
DOE OSTI2024

Rotating cylinder electrode in reactive CO 2 capture: Identifying active C species via transport, VLE models and kinetics

Here, this article explores technical challenges and potential methodologies for understanding electrochemical Reactive CO 2 Capture (RCC) mechanisms. RCC offers potential energy cost advantages by directly converting captured CO 2 into fuels and chemicals, unlike traditional carbon capture and utilization (CCU) processes that require sequential capture, concentration, and compression. However, direct conversion of captured CO 2 introduces complexity due to additional equilibrium buffer reactions, making it challenging to identify active species for reduction in electrochemical studies. This article discusses methods to integrate transport, thermodynamics, and kinetics concepts to identify active carbon sources in RCC. Vapor‐Liquid Equilibrium (VLE) and transport models are validated against experimental results obtained in a gastight rotating cylinder electrode reactor and are shown as useful tools for studying RCC in heterogeneous electrocatalysts across different capture agents, solvents, and temperatures. This article establishes an experimental framework for advancing research in electrochemical RCC.

Electrocatalysis↗
DOE OSTI2022

Modeling the thermodynamics of the FeTi hydrogenation under para-equilibrium: An ab-initio and experimental study

We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.

08 HYDROGEN↗
DOE OSTI2024

Rotational state-to-state transition rate coefficients for H 2 O + H 2 O collisions at nonequilibrium conditions

Aims.The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H 2 O deviates from local thermodynamic equilibrium. Methods.A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures,T kin andT rot . Results.Upon introducing our new parameterization of the collisional rates, taking into account their dependence on bothT kin andT rot , we find a change of up to 20% in the H 2 O rotational level populations for both ortho and para-H 2 O for the part of the cometary coma where the nonequilibrium regime occurs.

Astronomy & Astrophysics↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Sensitivity Analysis, Reduced-order Modeling, and Optimization of a Gas-Cooled Pebble Bed Reactor using Equilibrium-Core and DLOFC Performance

This work presents and applies a workflow for performing design optimization on gas-cooled pebble-bed reactors. Based on previous research, a representative equilibrium core of a pebble-bed reactor and a depressurized loss-of-forced-cooling model are created. These applications are built using the Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically utilizing Griffin, Pronghorn, and Bison. After defining design-related parameters and quantities of interest regarding reactor safety and efficiency, this multiphysics model is sampled using the MOOSE stochastic tools module. The result is a comprehensive dataset of configurations, enabling sensitivity analysis and the generation of reduced-order models. Subsequently, the dataset and reduced-order models are employed in an optimization study aimed at maximizing fuel utilization while adhering to safety and operational constraints. The optimization process leads to an improvement of fuel utilization by approximately 10\%, compared to engineering-judgment-based nominal conditions.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A new step forward in realistic cluster lens mass modelling: analysis of Hubble Frontier Field Cluster Abell S1063 from joint lensing, X-ray, and galaxy kinematics data

We present a new method to simultaneously and self-consistently model the mass distribution of galaxy clusters that combines constraints from strong lensing features, X-ray emission, and galaxy kinematics measurements. We are able to successfully decompose clusters into their collisionless and collisional mass components thanks to the X-ray surface brightness, as well as use the dynamics of cluster members, to obtain more accurate masses exploiting the fundamental plane of elliptical galaxies. Knowledge from all observables is included through a consistent Bayesian approach in the likelihood or in physically motivated priors. We apply this method to the galaxy cluster Abell S1063 and produce a mass model that we publicly release with this paper. The resulting mass distribution presents different ellipticities for the intra-cluster gas and the other large-scale mass components as well as deviation from elliptical symmetry in the main halo. We assess the ability of our method to recover the masses of the different elements of the cluster using a mock cluster based on a simplified version of our Abell S1063 model. Thanks to the wealth of mutliwavelength information provided by the mass model and the detected X-ray emission, we also found evidence for an ongoing merger event with gas sloshing from a smaller infalling structure into the main cluster. In agreement with previous findings, the total mass, gas profile, and gas mass fraction are all consistent with small deviations from the hydrostatic equilibrium. This new mass model for Abell S1063 is publicly available, as the lenstool extension used to construct it.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Muted extratropical low cloud seasonal cycle is closely linked to underestimated climate sensitivity in models

Abstract A large spread in model estimates of the equilibrium climate sensitivity (ECS), defined as the global mean near-surface air-temperature increase following a doubling of atmospheric CO 2 concentration, leaves us greatly disadvantaged in guiding policy-making for climate change adaptation and mitigation. In this study, we show that the projected ECS in the latest generation of climate models is highly related to seasonal variations of extratropical low-cloud fraction (LCF) in historical simulations. Marked reduction of extratropical LCF from winter to summer is found in models with ECS > 4.75 K, in accordance with the significant reduction of extratropical LCF under a warming climate in these models. In contrast, a pronounced seasonal cycle of extratropical LCF, as supported by satellite observations, is largely absent in models with ECS < 3.3 K. The distinct seasonality in extratropical LCF in climate models is ascribed to their different prevailing cloud regimes governing the extratropical LCF variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A non-cooperative meta-modeling game for automated third-party calibrating, validating and falsifying constitutive laws with parallelized adversarial attacks

The evaluation of constitutive models, especially for high-risk and high-regret engineering applications, requires efficient and rigorous third-party calibration, validation and falsification. While there are numerous efforts to develop paradigms and standard procedures to validate models, difficulties may arise due to the sequential, manual, and often biased nature of the commonly adopted calibration and validation processes, thus slowing down data collections, hampering the progress towards discovering new physics, increasing expenses and possibly leading to misinterpretations of the credibility and application ranges of proposed models. This work attempts to introduce concepts from game theory and machine learning techniques to overcome many of these existing difficulties. Here, we introduce an automated meta-modeling game where two competing AI agents systematically generate experimental data to calibrate a given constitutive model and to explore its weakness such that the experiment design and model robustness can be improved through competitions. The two agents automatically search for the Nash equilibrium of the meta-modeling game in an adversarial reinforcement learning framework without human intervention. In particular, a protagonist agent seeks to find the more effective ways to generate data for model calibrations, while an adversary agent tries to find the most devastating test scenarios that expose the weaknesses of the constitutive model calibrated by the protagonist. By capturing all possible design options of the laboratory experiments into a single decision tree, we recast the design of experiments as a game of combinatorial moves that can be resolved through deep reinforcement learning by the two competing players. Our adversarial framework emulates idealized scientific collaborations and competitions among researchers to achieve a better understanding of the application range of the learned material laws and prevent misinterpretations caused by conventional AI-based third-party validation. Numerical examples are given to demonstrate the wide applicability of the proposed meta-modeling game with adversarial attacks on both human-crafted constitutive models and machine learning models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

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

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

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

Effect of Impurity Atoms on Thermonuclear Yield

The presence of impurity atoms in an inertial confinement fusion (ICF) implosion can greatly affect the burning of the thermonuclear fuel. This is a serious concern for attaining ignition as high atomic number atoms, such as carbon, can mix into the hot spot and reach the deuterium–tritium (DT) fuel, absorb energy from the implosion, and radiate away that energy. The temperature of the DT fuel drops and reduces the reactivity of the fuel resulting in lower thermonuclear yield. Fuel impurities can also be a useful diagnostic as when higher-Z atoms, such as Ar, Kr, and Xe, are purposefully mixed into the fuel. From the spectra emitted by these atoms, the temperature and density of the fuel can be deduced. The long history of using such tracer atoms by the national ICF program shows the great advances in experimental technique, diagnostic instruments, and atomic-physics modeling over the past 50 years. A particular set of experiments, the “High-Z” campaign, is reviewed in depth as this was the most systematic study of tracer element composition and quantity ever reported. Difficulties in interpreting the results of those experiments due to the limited code capabilities at that time are discussed. A main concern for the High-Z campaign was the then current state of non-Local Thermodynamic Equilibrium (nLTE) atomic physics models. The present state of nLTE modeling and its implementation in Los Alamos’s multi-physics codes is discussed in this report. To further the understanding of the complex interaction of high-Z atoms in the gas, a survey of relevant existing experimental platforms at the National Ignition Facility (NIF) and the Omega Laser Facility are presented along with the current state of pertinent diagnostics at NIF and Omega. This review concludes with a discussion of the trade-off available between temperature and volume of the imploded ICF capsule with the goal of maximizing the total yield.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Minimizing the Electromechanical Stresses in Poloidal Field Coils by Optimizing their Numbers and Locations using FREDA Framework

Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.

Hassan, Ehab [ORNL] (ORCID:0000000181060301)↗
DOE OSTI2022

Probing physics beyond the standard model: limits from BBN and the CMB independently and combined

We present new Big Bang Nucleosynthesis (BBN) limits on the cosmic expansion rate or relativistic energy density, quantified via the number N ν of equivalent neutrino species. We use the latest light element observations, neutron mean lifetime, and update our evaluation for the nuclear rates d + d → 3 He + n and d + d → 3 H+p. Combining this result with the independent constraints from the cosmic microwave background (CMB) yields tight limits on new physics that perturbs N ν and η prior to cosmic nucleosynthesis: a joint BBN+CMB analysis gives N ν = 2.898 ± 0.141, resulting in N ν < 3.180 at 2σ. We apply these limits to a wide variety of new physics scenarios including right-handed neutrinos, dark radiation, and a stochastic gravitational wave background. The strength of the independent BBN and CMB constraints now opens a new window: we can search for limits on potential changes in N ν and/or the baryon-to-photon ratio η between the two epochs. The present data place strong constraints on the allowed changes in N ν between BBN and CMB decoupling; for example, we find -0.708 < N ν CMB - N ν BBN < 0.328 in the case where η and the primordial helium mass fraction Yp are unchanged between the two epochs; we also give limits on the allowed variations in η or in (η, N ν ) jointly. We discuss scenarios in which such changes could occur, and show that BBN+CMB results combine to place important constraints on some early dark energy models to explain the H 0 tension. Looking to the future, we forecast the tightened precision for N ν arising from both CMB Stage 4 measurements as well as improvements in astronomical 4 He measurements. Here, we find that CMB-S4 combined with present BBN and light element observation precision can give σ(N ν ) ≃ 0.03. Such future precision would reveal the expected effect of neutrino heating (N eff -3 = 0.044) of the CMB during BBN, and would be near the level to reveal any particle species ever in thermal equilibrium with the standard model. Improved Y p measurements can push this precision even further.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2024

A Vehicle-to-Grid planning framework incorporating electric vehicle user equilibrium and distribution network flexibility enhancement

The rapid surge in electric vehicle (EV) adoption, coupled with advancements in charging technologies, emphasizes the critical necessity for expanding EV recharging infrastructure. Simultaneously, the Distribution Network (DN) encounters escalating challenges in meeting charging demand during peak traffic periods. Consequently, there is a mounting demand for the deployment of innovative Vehicle-to-Grid (V2G) technologies to augment the DN’s flexibility in power dispatch and alleviate travel costs for EV users. Hence, this paper proposes an EV-user-equilibrium-(UE)-constrained V2G planning framework that enhances flexibility in the DN. The framework aims to ascertain the optimal placement and capacity of EV charging stations (EVCSs) and V2G charging piles within the Transportation Network (TN). It takes into account the equilibrium condition stemming from competitive EV charging and routing behaviors alongside the optimal expansion of DN energy resources to accommodate the electricity supplied by the V2G piles. This study commences by analyzing EV drivers’ travel decisions, considering the influence of charging and V2G pile locations and sizes. Subsequently, we tackle the Traffic Assignment Problem with User Equilibrium (TAP-UE) model to characterize the steady-state traffic flow distribution of EVs. Following this, we formulate the optimization model for the Coordinated Power and Transportation Network (CPTN), which encompasses the optimal expansion of DN facilities and traffic flow regulation under UE conditions. To mitigate the computational complexity associated with the V2G planning model, we introduce a series of linearization methods to obtain a manageable Mixed-Integer Linear Programming (MILP) solution. Finally, to validate the efficacy of our proposed planning framework, we apply it to two test systems, including a real-world case study. Through these case studies, we explore the necessity and potential benefits of V2G technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗