Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “three phase 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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Phase coexistence implications of violating Newton’s third law

Newton’s third law, action = reaction, is a foundational statement of classical mechanics. However, in natural and living systems, this law appears to be routinely violated for constituents interacting in a nonequilibrium environment. Here, we use computer simulations to explore the macroscopic phase behavior implications of breaking microscopic interaction reciprocity for a simple model system. We consider a binary mixture of attractive particles and introduce a parameter that is a continuous measure of the degree to which interaction reciprocity is broken. In the reciprocal limit, the species are indistinguishable, and the system phase separates into domains with distinct densities and identical compositions. Furthermore, increasing nonreciprocity is found to drive the system to explore a rich assortment of phases, including phases with strong composition asymmetries and three-phase coexistence. Many of the states induced by these forces, including traveling crystals and liquids, have no equilibrium analogs. By mapping the complete phase diagram for this model system and characterizing these unique phases, our findings offer a concrete path forward toward understanding how nonreciprocity shapes the structures found in living systems and how this might be leveraged in the design of synthetic materials.

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

Ultrahigh-pressure crystallographic passage towards metallic hydrogen

The structural evolution of molecular hydrogen H 2 under multi-megabar compression and its relation to atomic metallic hydrogen is a key unsolved problem in condensed-matter physics. Although dozens of crystal structures have been proposed by theory, only one, the simple hexagonal-close-packed (hcp) structure of only spherical disordered H 2 , has been previously confirmed in experiments. Through advancing nano-focused synchrotron X-ray probes, here we report the observation of the transition from hcp H 2 to a post-hcp structure with a six-fold larger supercell at pressures above 212 GPa, indicating the change of spherical H 2 to various ordered configurations. Theoretical calculations based on our XRD results found a time-averaged structure model in the space group $P\bar{6}2c$ with alternating layers of spherically disordered H 2 and new graphene-like layers consisting of H 2 trimers (H 6 ) formed by the association of three H 2 molecules. Here, this supercell has not been reported by any previous theoretical study for the post-hcp phase, but is close to a number of theoretical models with mixed-layer structures. The evidence of a structural transition beyond hcp establishes the trend of H 2 molecular association towards polymerization at extreme pressures, giving clues about the nature of the molecular-to-atomic transition of metallic hydrogen. Considering the spectroscopic behaviours that show strong vibrational and bending peaks of H 2 up to 400 GPa, it would be prudent to speculate the continuation of hydrogen molecular polymerization up to its metallization.

Ji, Cheng [Center for High Pressure Science and Te↗
DOE OSTI2024

Gravitational production of heavy particles during and after inflation

Abstract We investigate the gravitational production of a scalar fieldχwith a mass exceeding the Hubble scale during inflationm χ ≳H I , employing both analytical and numerical approaches. We demonstrate that the steepest descent method effectively captures the epochs and yields of gravitational production in a compact and simple analytical framework. These analytical results align with the numerical solutions of the field equation. Our study covers three spacetime backgrounds: de Sitter, power-law inflation, and the Starobinsky inflation model. Within these models, we identify two distinct phases of particle production: during and after inflation. During inflation, we derive an accurate analytic expression for the particle production rate, accounting for a varying Hubble rate. After inflation, the additional burst of particle production depends on the inflaton mass around its minimum. When this mass is smaller than the Hubble scale during inflation,H I , there is no significant extra production. However, if the inflaton mass is larger, post-inflation production becomes the dominant contribution. Furthermore, we explore the implications of gravitationally produced heavy fields for dark matter abundance, assuming their cosmological stability.

Physics↗
DOE OSTI2025

Isoconversional Kinetic Analysis of Oxygen Desorption for Sr 0.75 Ca 0.25 FeO 3-δ Perovskite in Dry and Steam-Based Environments

The desorption kinetics of Sr 0.75 Ca 0.25 FeO 3-δ perovskite material are examined in both dry and steam-based environments using redox gaseous products of a laboratory-scale fixed bed. First, desorption kinetics associated with the dry environment are proposed based on Friedman’s isoconversional method. It is found from the reconstructed reaction model that the desorption kinetics of this perovskite are controlled by a three-step mechanism, where hypothetically, the phase-boundary reaction rapidly prevails first, followed by diffusion and/or nucleation, and finally, unimolecular decay or random nucleation-controlled reaction. Subsequent analyses of the Arrhenius parameters indicated that both the early and later desorption stages, respectively, controlled by the phase boundary and unimolecular/random nucleation reactions, are associated with low energetic demand. In contrast, a high energetic penalty is required for the diffusion/nucleation reaction. Further, a satisfactory a priori verification of the proposed kinetics suggested that the three-step reaction mechanism reasonably describes the desorption of this perovskite in a dry environment. With the presence of steam in the desorption environment, oxygen production is substantially inhibited, but the hypothetical three-step mechanism is still found to control the desorption. Meanwhile, the primary effects of the steam on these mechanisms include the widening of the conversion extent associated with the diffusion/nucleation reaction, along with a higher energetic demand compared to the dry environment.

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

Seasonal representation of extreme precipitation indices over the United States in CMIP6 present-day simulations

Realistically representing the present-day characteristics of extreme precipitation has been a challenge for global climate models, which is due in part to deficiencies in model resolution and physics, but is also due to a lack of consistency in gridded observations. In this study, we use three observation datasets, including gridded rain gauge and satellite data, to assess historical simulations from sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models. We separately evaluate summer and winter precipitation over the United States (US) with a comprehensive set of extreme precipitation indices, including an assessment of precipitation frequency, intensity and spatial structure. The observations exhibit significant differences in their estimates of area-average intensity distributions and spatial patterns of the mean and extremes of precipitation over the US. In general, the CMIP6 multi-model mean performs better than most individual models at capturing daily precipitation distributions and extreme precipitation indices, particularly in comparison to gauge-based data. Also, the representation of the extreme precipitation indices by the CMIP6 models is better in the summer than winter. Although the ‘standard’ horizontal-resolution can vary significantly across CMIP6 models, from ~0.7° to ~2.8°, we find that resolution is not a good indicator of model performance. Overall, our results highlight common biases in CMIP6 models and demonstrate that no single model is consistently the most reliable across all indices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Field-induced intermediate ordered phase and anisotropic interlayer interactions in α-RuCl 3

In α-RuCl 3 , an external magnetic field applied within the honeycomb plane can induce a transition from a magnetically ordered state to a disordered state that is potentially related to the Kitaev quantum spin liquid. In zero field, single crystals with minimal stacking faults display a low-temperature state with in-plane zigzag antiferromagnetic order and a three-layer periodicity in the direction perpendicular to the honeycomb planes. In this work, we present angle-dependent magnetization, ac susceptibility, and thermal transport data that demonstrate the presence of an additional intermediate-field ordered state at fields below the transition to the disordered phase. Neutron-diffraction results show that the magnetic structure in this phase is characterized by a six-layer periodicity in the direction perpendicular to the honeycomb planes. Theoretically, the intermediate ordered phase can be accounted for by including spin-anisotropic couplings between the layers in a three-dimensional spin model. Together, this demonstrates the importance of interlayer exchange interactions in α-RuCl 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Relationships between Immersion Freezing and Crystal Habit for Arctic Mixed-Phase Clouds—A Numerical Study

Abstract The number concentration of ice particles in Arctic mixed-phase clouds is a major controlling factor of cloud lifetime. The relationships between ice nucleation mode and ice crystal habit development are not yet constrained by observations. This study uses a habit-predicting microphysical scheme within a 3D large-eddy simulation model to evaluate the relationship between immersion freezing and ice habit in a simulated Arctic mixed-phase cloud case. Three immersion freezing parameterizations are considered: a volume-dependent freezing scheme (VF), a parameterization limited to activated droplets (C-AC), and a parameterization limited to coarse aerosol particles (C-CM). Both C-AC and C-CM are based on classical nucleation theory. The freezing rate with VF is found to be greater in downdraft regions than in updraft regions due to the downdraft having a higher number concentration of large droplets. The C-AC cases show active freezing of small droplets near cloud top, whereas in the C-CM cases, mainly the 8–32- μ m-sized droplets freeze in updraft regions near the cloud base. Because the initial crystal size is assumed to affect the axis ratio of hexagonal plates, the VF cases produce crystals with larger axis ratios, resulting in smaller mode radii than the C-AC cases. In all cases, irregular polycrystals dominate near cloud top and a band-like structure develops within the cloud, which qualitatively agrees with previous observations. In the VF and C-CM cases, unactivated large droplets arising from coarse-mode aerosol particles contributed significantly to the freezing rate, producing an important influence on crystal habit.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2021

SU(4) Chiral Spin Liquid, Exciton Supersolid, and Electric Detection in Moiré Bilayers

We propose a moiré bilayer as a platform where exotic quantum phases can be stabilized and electrically detected. Moiré bilayers consist of two separate moiré superlattice layers coupled through the interlayer Coulomb repulsion. In the small distance limit, an SU(4) spin can be formed by combining layer pseudospin and the real spin. As a concrete example, we study an SU(4) spin model on triangular lattice in the fundamental representation. By tuning a three-site ring exchange term K ~ ( t 3 /U 2 ), we find the SU(4) symmetric crystallized phase and an SU (4) 1 chiral spin liquid at the balanced filling. We also predict two different exciton supersolid phases with interlayer coherence at imbalanced filling under displacement field. Especially, the system can simulate an SU(2) Bose-Einstein condensation by injecting interlayer excitons into the magnetically ordered Mott insulator at the layer polarized limit. Smoking gun evidences of these phases can be obtained by measuring the pseudospin transport in the counterflow channel.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

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

Exotic Phases of a Higgs-Yukawa Model with Reduced Staggered Fermions

We investigate the phase structure of a four dimensional SO(4) invariant lattice Higgs-Yukawa model comprising four reduced staggered fermions interacting with a real scalar field. The fermions belong to the fundamental representation of the symmetry group while the three scalar field components transform in the self-dual representation of SO(4). We explore the phase diagram and find evidence of a continuous transition between a phase where the fermions are massless to one where the fermions acquire mass. This transition is not associated with symmetry breaking and there is no obvious local order parameter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

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

Causality-respecting adaptive refinement for PINNs: enabling precise interface evolution in phase field modeling

Physics-informed neural networks (PINNs) have emerged as a powerful tool for solving physical systems described by partial differential equations (PDEs). However, their accuracy in dynamical systems, particularly those involving sharp moving boundaries with complex initial morphologies, remains a challenge. Here, this study introduces an approach combining residual-based adaptive refinement (RBAR) with causality-informed training to enhance the performance of PINNs in solving spatio-temporal PDEs. Our method employs a three-step iterative process: initial causality-based training, RBAR-guided domain refinement, and subsequent causality training on the refined mesh. Applied to the Allen-Cahn equation, a widely-used model in phase field simulations, our approach demonstrates significant improvements in solution accuracy and computational efficiency over traditional PINNs. Notably, we observe an ‘overshoot and relocate’ phenomenon in dynamic cases with complex morphologies, showcasing the method’s adaptive error correction capabilities. This synergistic interaction between RBAR and causality training enables accurate capture of interface evolution, even in challenging scenarios where traditional PINNs fail. Our framework not only resolves the limitations of uniform refinement strategies but also provides a generalizable methodology for solving a broad range of spatio-temporal PDEs. The enhanced performance of the RBAR–causality combined framework demonstrates its strong potential for advancing PINN-based modeling of physical systems characterized by complex, evolving interfaces.

Allen-Cahn equations↗
DOE OSTI2023

Elucidating the active phases of CoOx films on Au(111) in the CO oxidation reaction

Abstract Noble metals supported on reducible oxides, like CoO x and TiO x , exhibit superior activity in many chemical reactions, but the origin of the increased activity is not well understood. To answer this question we studied thin films of CoO x supported on an Au(111) single crystal surface as a model for the CO oxidation reaction. We show that three reaction regimes exist in response to chemical and topographic restructuring of the CoO x catalyst as a function of reactant gas phase CO/O 2 stoichiometry and temperature. Under oxygen-lean conditions and moderate temperatures (≤150 °C), partially oxidized films (CoO x<1 ) containing Co 0 were found to be efficient catalysts. In contrast, stoichiometric CoO films containing only Co 2+ form carbonates in the presence of CO that poison the reaction below 300 °C. Under oxygen-rich conditions a more oxidized catalyst phase (CoO x>1 ) forms containing Co 3+ species that are effective in a wide temperature range. Resonant photoemission spectroscopy (ResPES) revealed the unique role of Co 3+ sites in catalyzing the CO oxidation. Density function theory (DFT) calculations provided deeper insights into the pathway and free energy barriers for the reactions on these oxide phases. These findings in this work highlight the versatility of catalysts and their evolution to form different active phases, both topological and chemically, in response to reaction conditions exposing a new paradigm in the catalyst structure during operation.

Chen, Hao (ORCID:0000000337598352)↗
DOE OSTI2022

Loop Thermosyphon Enhanced Solar Collector

A two-phase loop thermosyphon solar collector was developed to transfer thermal energy efficiently and passively from a concentrated solar collector to a thermal desalination process. The two-phase flow in the loop thermosyphon is driven by a difference in gravitational pressure head between the liquid return line and the two-phase riser located between the evaporator and condenser. A numerical model incorporating mass, energy, and pressure drop balances was developed to predict system performance and use as a design tool. Three iterations of loop thermosyphon solar collectors were designed, fabricated, and tested to validate the numerical analysis and demonstrate system integration. Additionally, a novel solar receiver design was evaluated based on volumetric absorption of the concentrated sunlight by an optically absorptive two-phase working fluid. Life tests were performed to evaluate the working fluid thermal stability and resistance to ultraviolet degradation. A glass receiver tube and associated glass-to-metal junctions were developed. Finally, a full-scale loop thermosyphon was integrated with a parabolic trough solar collector and tested on-sun.

14 SOLAR ENERGY↗
DOE OSTI2025

Data‐Driven Insights into Rare Earth Mineralization: Machine Learning Applications Using Functional Material Synthesis Data

Understanding rare‐earth element (REE) mineralization mechanisms is essential for developing efficient separation strategies. Although the geochemical pathways that generate REE deposits are qualitatively known, quantitative links between specific conditions and mineralization outcomes remain limited. Herein, the repurpose laboratory REE hydrothermal synthesis data—originally collected for functional‐materials fabrication—as a surrogate for studying mineralization with data‐driven methods. The compiled 1,200+ hydrothermal reaction records and trained three machine‐learning models—K‐nearest neighbors (KNN), random forest (RF), and extreme gradient boosting (XGB)—to predict product elements and phases from precursors, additives, reaction conditions, and engineered features. Validation shows XGB achieves the highest accuracy. Feature importance indicates thermodynamic properties of cations and anions dominate model decisions. Correlations reveal positive relationships among precursor concentration, reaction time, pH, and temperature, consistent with classical crystallization behavior. XGB‐based regressors are built to predict crystallization temperature and pH from precursor/product attributes. Performance is strongest when similar training examples exist, while accuracy declines for underrepresented reactions, notably REE carbonates and heavy‐REE systems. Overall, the study shows that functional‐materials datasets can illuminate REE mineralization and provide priors for exploration and processing. Expanding datasets with less‐studied chemistries and conditions will improve generality and support deposit discovery and more efficient REE recovery.

feature importance analysis↗
DOE OSTI2023

Tuning the Radius Ratio to Enhance Thermoelectric Properties in the Zintl Compounds AM 2 Sb 2 (A = Ba, Sr; M = Zn, Cd)

Five novel Zintl phase solid solutions in the Ba 1–x Sr x Zn 2–y Cd y Sb 2 (0 ≤ x ≤ 0.13(1); 0 ≤ y ≤ 0.32(2)) system were successfully synthesized by the molten Pb metal-flux method, and the powder X-ray diffraction and single-crystal X-ray diffraction analyses proved that all five title compounds adopted the BaCu 2 S 2 -type phase having the orthorhombic Pnma space group (Z = 4, Pearson code oP20) with five crystallographically independent atomic sites. The previously studied BaCu 2 S 2 -type antimonides demonstrated a limited tolerance for doping in contrast to the CaAl 2 Si 2 -type antimonides. To understand the relatively narrower phase width and limited dopability of the title BaCu 2 S 2 -type phase than the CaAl 2 Si 2 -type phase in the overall Ba 1–x Sr x Zn 2–y Cd y Sb 2 system, the radius ratio of cations and anionic elements r + /r – for two structure types were thoroughly investigated. For the first time, the r + /r – ratio was identified as a critical factor for the phase selectivity: (1) r + /r – > 1 favored the BaCu 2 S 2 -type phase, and (2) r + /r – < 1 favored the CaAl 2 Si 2 -type phase. Further, we also revealed the structural transformation mechanism from the more widely observed CaAl 2 Si 2 -type phase to the title BaCu 2 S 2 -type phase as the relatively larger cationic elements were introduced to the system. A series of DFT calculations using the three hypothetical models indicated that a resonance peak near EF in the density of states curves was descended from the relatively flat band structure at several special symmetry points rationalizing the enhanced Seebeck coefficients of Ba 0.94(1) Sr 0.06 Zn 1.86(3) Cd 0.14 Sb 2 and Ba 0.96(1) Sr 0.04 Zn 1.68(2) Cd 0.32 Sb 2 . Electron localization function analysis rationalized the correlation between the polarity change of anionic Zn/Cd–Sb bonds and the charge carrier mobility on the anionic frameworks. Temperature-dependent thermoelectric properties were studied for the four title compounds, and the results proved that the Sr and Cd doping in the title Ba 1–x Sr x Zn 2–y Cd y Sb 2 system successfully enhanced the ZT values through the increased Seebeck coefficients and the reduced total thermal conductivities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗