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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Experimental Liquidus Study of the Ternary CaO-ZnO-SiO{sub 2} System

Phase equilibria of the ternary CaO-ZnO-SiO{sub 2} system have been investigated at 1170 °C to 1691 °C for oxide liquid in equilibrium with air and solid oxide phases: tridymite or cristobalite SiO{sub 2} (up to two immiscible liquids), pseudowollastonite (CS) CaSiO{sub 3}, rankinite (C{sub 3}S{sub 2}) Ca{sub 3}Si{sub 2}O{sub 7}, dicalcium silicate (C{sub 2}S) (Ca, Zn){sub 2}SiO{sub 4}, tricalcium silicate (C{sub 3}S) (Ca, Zn){sub 3}SiO{sub 5}, lime (Ca, Zn)O, zincite (Zn, Ca)O, willemite Zn{sub 2}SiO{sub 4} and hardystonite (melilite) Ca{sub 2}ZnSi{sub 2}O{sub 7}, covering the ranges of concentrations not studied before. High-temperature equilibration on primary phase (silica) or inert metal (platinum) substrates followed by quenching and direct measurement of the Ca, Zn and Si concentrations in the phases with the electron probe X-ray microanalysis (EPMA) has been used to accurately characterize the system. Liquidus phase equilibrium data of the present authors for the CaO-ZnO-SiO{sub 2} system are essential to obtain a self-consistent set of parameters of thermodynamic models for all phases.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bubbles in US gasoline prices: Assessing the role of hurricanes and anti–price gouging laws

This paper examines the presence of bubbles — price changes that cannot be attributed to market fundamentals — in weekly US retail gasoline prices for ten US cities, and whether these bubbles are affected by hurricane events and anti–price gouging laws. Rolling-window, right-tailed unit root tests were used to identify and date bubbles. The analysis found 6 to 22 bubbles, lasting 4 weeks or longer, in city-level retail gasoline prices between 2000 and 2017. A linear regression model of the adjustment speed coefficients from the bubble identification step found that hurricane events increase the rate of price explosivity (or slow the rate of return to equilibrium) following a shock. Furthermore, the regression estimates also imply that state anti–price gouging laws can counteract these effects.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Synchrotron Imaging of Pore Formation in Li Metal Solid-State Batteries Aided by Machine Learning

High-rate capable, reversible lithium metal anodes are necessary for next generation energy storage systems. In situ tomography of Li/LLZO/Li cells is carried out to track morphological transformations in Li metal electrodes. Machine learning enables tracking the lithium metal morphology during galvanostatic cycling. Nonuniform lithium electrode kinetics are observed at both electrodes during cycling. Hot spots in lithium metal are correlated with microstructural anisotropy in LLZO. Mesoscale modeling reveals that regions with lower effective properties (transport and mechanical) are nuclei for failure. Advanced visualization combined with electrochemistry represents an important pathway toward resolving non-equilibrium effects that limit rate capabilities of solid-state batteries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structure–Property Correlation of Hierarchically Porous Carbons for Fluorocarbon Adsorption

Although traditional commercially available porous carbon-fluorocarbon working pairs have shown promising applicability for adsorption cooling, advancements in synthetic carbons may further improve performance. Moreover, insights into structure-property relationships that target higher sorption capacities within these synthesized carbons may guide such materials' future design. We utilized hierarchically porous synthetic carbons (HPCs) with colossal microporous and mesoporous characterized by high surface areas (up to 2689 m 2 /g) and pore volumes (up to 10.31 cm 3 /g) towards fluorocarbon R134a adsorption. This unique pore topology leads to exceptional R134a uptake, ~250 wt.%, outperforming the highest uptake carbon material to date, Maxsorb III (~220 wt.%). Additionally, material characterizations reveal that the outstanding R134a capacity may be attributed to textural properties and oxygen-terminated functional groups more than graphitization of the material. Most importantly, HPCs are efficiently utilized in a two-bed model chiller device, where the performance shows excellent working capacity (105 wt%, ~1.5 times the value of reported carbon materials/R134a). Fluorocarbons adsorption on HPCs also displays fast kinetics (equilibrium time: ~2 min) mainly driven by physical adsorption (Qst: ~27 kJ/mol), characteristic of swiftly reversible behavior adsorption-desorption behaviors. This work provides a fundamental understanding of the applicability of HPCs/R134a working pair for adsorption cooling.

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

Methanol Decomposition on Copper Surfaces under Ambient Conditions: Mechanism, Surface Kinetics, and Structure Sensitivity

Here, we study the adsorption of methanol vapor under ambient pressure and temperature conditions on low-index Cu surfaces using surface-sensitive infrared (IR) and X-ray spectroscopy techniques. The first step of methanol decomposition, i.e. , breaking of the O—H bond to form surface-bound methoxy, readily occurs under ambient conditions. Time-lapse IR spectra clearly indicate a gradually decreasing methoxy coverage, which does not obey well established kinetic models. We rationalize the initial temperature-independent, high, nonequilibrium coverage of methoxy by a H-bonded methanol assembly in the precursor state. A temperature-dependent equilibrium coverage is achieved as the excess methoxy is eliminated gradually via further dehydrogenation to CO that desorbs to the gas phase. The kinetics of this process displays a significant structure sensitivity with considerably faster kinetics on the Cu(110) surface compared to Cu(111) and Cu(100) surfaces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Biased Estimates of Equilibrium Climate Sensitivity and Transient Climate Response Derived From Historical CMIP6 Simulations

Abstract This study assesses the effective climate sensitivity (EffCS) and transient climate response (TCR) derived from global energy budget constraints within historical simulations of eight CMIP6 global climate models (GCMs). These calculations are enabled by use of the Radiative Forcing Model Intercomparison Project (RFMIP) simulations, which permit accurate quantification of the radiative forcing. Long‐term historical energy budget constraints generally underestimate EffCS from CO 2 quadrupling and TCR from CO 2 ramping, owing to changes in radiative feedbacks and changes in ocean heat uptake efficiency. Atmospheric GCMs forced by observed warming patterns produce lower values of EffCS that are more in line with those inferred from observed historical energy budget changes. The differences in the EffCS estimates from historical energy budget constraints of models and observations are traced to discrepancies between modeled and observed historical surface warming patterns.

Dong, Yue↗
DOE OSTI2024

Visualizing the SEI formation between lithium metal and solid-state electrolyte

The solid electrolyte interphase (SEI) is regarded as the most important factor affecting the durability of lithium-metal anode in all-solid-state batteries (ASSBs). Despite its significance, the nucleation and growth mechanism of SEI is not yet well understood. Here, we elucidate the thermodynamics and kinetics governing SEI formation at the Li|β-Li 3 PS 4 interface at the atomic scale via thermodynamic phase equilibrium analysis and machine-learning-potential-assisted molecular dynamics (MD) simulations. An accurate moment tensor potential using the machine-learning method is developed for a reactive model of Li|β-Li 3 PS 4 . This potential enabled us to perform large-scale MD simulations with the model size expanded to the experimental dimensions (~40 nm) while maintaining the same level of accuracy as density functional theory calculations. The results reveal a four-stage evolution process at the Li|β-Li 3 PS 4 interface, namely (i) fast ion diffusion, (ii) nucleation, (iii) Li 2 S growth, and (iv) stabilization. Notably, we demonstrate that the SEI can be categorized into crystalline and amorphous regions. The simulated SEI thickness, structure, and composition closely match experimental findings, validating the accuracy of the MD simulations. We further disclose the significant impact of ion diffusion kinetic limitations on the phase formation and crystallization of interfacial products. Furthermore, we shed light on the detailed potential energy (PE) distribution of lithium along the direction perpendicular to the interface. Finally, this information is crucial for better understanding interfacial ion mobility.

25 ENERGY STORAGE↗
DOE OSTI2020

Modeling radiative-shocks created by laser–cluster interactions

Radiative-shocks induced by laser–cluster interactions are modeled using radiation-hydrodynamic simulations. A good agreement—in both shock velocity and density profiles—is obtained between experiment and simulations, indicating that non-local thermodynamic equilibrium (NLTE) radiative effects are important in the experimental regime examined, particularly at early times (≤30 ns) due to the elevated temperatures (≥35 eV). The enhanced NLTE radiative emission causes the shock to be reduced in amplitude, increased in width, and reduced in propagation velocity, while the amplitude of the radiative precursor is increased. As the density and temperature conditions are relatively modest, this potentially has important implications for the scalings that are used in laboratory–astrophysics to transform between laboratory and astrophysical scales, which do not hold for non-LTE systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A reduced-model (nSOLT) simulation of neutral recycling effects on plasma turbulence in the divertor region of MAST-U

The 2D scrape-off-layer turbulence code (nSOLT), which includes 1D kinetic neutral–plasma interactions, is applied to study effects of neutral recycling on plasma turbulence for parameters illustrative of the MAST-U divertor region. Neutral recycling is modeled by injecting a fraction of the parallel plasma flux to the divertor back into the simulation domain as a source of Franck–Condon-distributed neutrals. Stationary sources, concentrated at the magnetic separatrix, model plasma streaming into the divertor region from the upstream scrape-off-layer and sustain plasma turbulence absent neutral recycling. Starting from one such no-neutrals equilibrium, we initiate recycling in a numerical experiment designed to diagnose and identify the effects of various neutral–plasma interactions on the divertor plasma, divertor turbulence, and plasma exhaust. The onset of recycling triggers an initial burst of enhanced cross field plasma transport that is quelled by ionization cooling and charge–exchange (CX) friction, with growing neutral pressure, leading to a quiescent, turbulence-free state. Diagnosis of this transient burst reveals that (1) the sudden increase in plasma density due to ionization dominates the onset of the burst, (2) electron cooling due to ionization increases collisionality and disconnects blob filaments from the sheath, and (3) CX friction drives tripole polarization of a blob that can dominate the curvature-driven dipole polarization, leading to the stagnation of blob propagation and reduced radial turbulent transport. It is shown that CX friction is negligible compared to sheath physics in determining equilibrium mean flow shearing rates, for parameters considered herein (specifically a short connection length to the divertor target), while it can significantly reduce interchange-instability growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Laser pulse driven control of charge and spin order in the two-dimensional Kondo lattice

Fast and dynamical control of quantum phases is highly desired in the application of quantum phenomena in devices. On the experimental side, ultrafast laser pulses provide an ideal platform to induce femtosecond dynamics in a variety of materials. Here we show that a laser pulse driven heavy fermion system can be tuned to dynamically evolve into a phase, which is not present in equilibrium. Using the state-of-the-art time-dependent variational Monte Carlo method, we perform numerical simulations of realistic laser pulses applied to the Kondo lattice model. By tracking spin and charge degrees of freedom we identify a dynamical phase transition, within the partially Kondo screened phase, from a charge ordered into a metallic state, while preserving the spin order of the system. Furthermore, we propose using high-harmonic generation and nonequilibrium optical conductivity to identify the dynamical phase transition in an ultrafast experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Probabilistic measures for biological adaptation and resilience

This paper introduces an approach to quantifying ecological resilience in biological systems, particularly focusing on noisy systems responding to episodic disturbances with sudden adaptations. Incorporating concepts from nonequilibrium statistical mechanics, we propose a measure termed “ecological resilience through adaptation,” specifically tailored to noisy, forced systems that undergo physiological adaptation in the face of stressful environmental changes. Randomness plays a key role, accounting for model uncertainty and the inherent variability in the dynamical response among components of biological systems. Our measure of resilience is rooted in the probabilistic description of states within these systems and is defined in terms of the dynamics of the ensemble average of a model-specific observable quantifying success or well-being. Our approach utilizes stochastic linear response theory to compute how the expected success of a system, originally in statistical equilibrium, dynamically changes in response to a environmental perturbation and a subsequent adaptation. Importantly, the resulting mathematical derivations allow for the estimation of resilience in terms of ensemble averages of simulated or experimental data. Finally, through a simple but clear conceptual example, we illustrate how our resilience measure can be interpreted and compared to other existing frameworks in the literature. The methodology is general but inspired by applications in plant systems, with the potential for broader application to complex biological processes.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Lower bound for the neutrino magnetic moment from kick velocities induced at the birth of neutron stars

We show that the neutrino chirality flip, which can take place in the core of a neutron star at birth, is an efficient process to allow neutrinos to anisotropically escape, thus providing a to induce the neutron star kick velocities. The process is not subject to the no-go theorem since although the flip from left- to right-handed neutrinos happens at equilibrium, the reverse process does not take place given that right-handed neutrinos do not interact with matter and therefore detailed balance is lost. For simplicity, we model the neutron star core as being made of strange quark matter. We find that the process is efficient when the neutrino magnetic moment is not smaller than [Formula: see text], where [Formula: see text] is the Bohr magneton. When this lower bound is combined with the most stringent upper bound, which uses the luminosity data obtained from the analysis of SN 1987A, our results set a range for the neutrino magnetic moment given by [Formula: see text]. The obtained kick velocities for natal conditions are consistent with the observed ones and span the correct range of radii for typical magnetic field intensities.

Physics↗
DOE OSTI2024

Quantifying the Temperature Dependence of the Multi-Species, Multi-Reaction Model. Part 1: Parameterization for a Meso-Carbon Micro-Bead Graphite

In this work, the temperature impact on the Multi-Species, Multi-Reaction (MSMR) model is studied. This is accomplished by acquiring data from slow rate lithiation and delithiation of a meso-carbon micro-bead (MCMB) graphite. The MSMR model is used to simulate linear-sweep voltammetry data of a porous electrode composed of graphite, and because the electrode is close to a state of dynamic equilibrium, the peaks in the differential voltage spectroscopy plot can be analyzed. Through this analysis, the temperature impact on the total fraction of available host sites in a particular MSMR gallery (X j ), the impact on the reference potential (U° j ), and the impact on the parameter detailing the deviation from Nernstian behavior ( ω j ) can be found. This is the first time the temperature dependence of the MSMR parameters have been experimentally analyzed. In Part 2, the impact of the temperature dependence of the MSMR parameters on the entropy coefficient of an intercalation material will be studied.

Paul, Abigail (ORCID:0000000172892069)↗
DOE OSTI2021

ORNL/Thermo4PFM

Library to compute equilibrium compositions and equal chemical potential compositions in binary and ternary alloys. This library supports thermodynamic data formulated according to CALPHAD (Computer Coupling of Phase Diagrams and Thermochemistry) models, as well as dilute binary alloys. One major aim is to solve the Kim-Kim-Suzuki (KKS) equations to obtain local compositions given a local mixture of two phases. These quantities are then used to evaluate the phase-field model driving force (right handside of time-evolution equation), as well as the energy density

Fattebert, Jean-Luc↗
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An experimental and modeling study of the application of membrane contactor reactors to methanol synthesis using pure CO 2 feeds

We describe here a post-combustion CO 2 capture and utilization (CCU) technology that converts the CO 2 into methanol (MeOH), a valuable chemical, thus providing a way to monetize the carbon captured to offset process costs. Methanol synthesis (MeS) has been discussed recently for application to CCU, but thermodynamic limitations make it difficult to convert in a single pass a large CO 2 fraction. Conventional catalysts show slow kinetics in converting CO 2 -rich syngas (or pure CO 2 ) into MeOH. Our Group developed a novel MeS process, employing a membrane contactor reactor (MCR) system that attains carbon conversions significantly higher than equilibrium. Our focus here is to process pure CO 2 streams by combining the MCR with a separate reactor, which converts the CO 2 into a syngas via the reverse water gas shift (RWGS) reaction. In this preliminary effort, the RWGS reactor (RWGSR) is assumed to reach equilibrium. Additional MeS kinetic rate data are generated validating experimentally the ability of the MeS-MCR to process as a feed the RWGSR exit stream. The performance of the combined (RWGSR/MeS-MCR) system is then simulated using a recently developed MeS-MCR model. In conclusion, the findings are encouraging, and research is currently ongoing to experimentally validate the RWGSR/MeS-MCR system performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗