Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transition Model”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Phase transitions in high-purity zirconium under dynamic compression

In this work, we present results from ramp compression experiments on high-purity Zr that show the α → ω, ω → β, as well as reverse β → ω phase transitions. Simulations with a multi-phase equation of state and phenomenological kinetic model match the experimental wave profiles well. While the dynamic α ! ! transition occurs ~ 9 GPa above the equilibrium phase boundary, the ω → β transition occurs within 0.9 GPa of equilibrium. We estimate that the dynamic compression path intersects the equilibrium ω – β line at P = 29:2 GPa, and T = 490 K. The thermodynamic path in the interior of the sample lies ~ 100 K above the isentrope at the point of the ω → β transition. Approximately half of this dissipative temperature rise is due to plastic work, and half is due to the non-equilibrium α → ω transition. The inferred rate of the α → ω transition is several orders of magnitude higher than that measured in dynamic diamond anvil cell (DDAC) experiments in an overlapping pressure range. We discuss a model for the influence of shear stress on the nucleation rate. We find that the shear stress s ji has the same effect on the nucleation rate as a pressure increase δP = cϵ ij s ji /(ΔV/V); where c is a geometric constant ~ 1 and, ϵ ij are the transformation shear strains. The small fractional volume change ΔV/V ≈ 0:1 at the α → ω transition amplifies the effect of shear stress, and we estimate that for this case δP is in the range of several GPa. Correcting our transition rate to a hydrostatic rate brings it approximately into line with the DDAC results, suggesting that shear stress plays a significant role in the transformation rate.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Theoretical antiferromagnetism of ordered face-centered cubic Cr-Ni alloys

Contrary to prior calculations, the Ni-rich ordered structures of the Cr-Ni alloy system are found to be antiferromagnetic under semilocal density-functional theory. The optimization of local magnetic moments significantly increases the driving force for the formation of CrNi 2 , the only experimentally observed intermetallic phase. This structure's ab initio magnetism appears well described by a Heisenberg Hamiltonian with longitudinal spin fluctuations; itinerant Cr moments are induced only by the strength of exchange interactions. The role of magnetism at temperature is less clear and several scenarios are considered based on a review of experimental literature, specifically a failure of the theory, the existence of an overlooked magnetic phase transition, and the coupling of antiferromagnetism to chemical ordering. In conclusion, implications for related commercial and high-entropy alloys are discussed for each case.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mapping causal pathways with structural modes fingerprint for perovskite oxides

Abstract Causality is innate to the determination of the fundamental mechanism controlling any physical phenomena. However, combining causality within the standard practices of computational modelling to understand structure-functionality connections is extremely rare. This work proposes a fingerprint based on key structural modes for ABO 3 -type perovskite oxides and its derivatives, combined with causal models, for predicting Kohn–Sham energies. Our study of causal models captures the inherent coupling between structural modes such as rotation, tilt and antiferroelectric displacements, responsible for phase transition, polarization, magnetization and metal–insulator transition, exhibited by these materials. Although developed for modelling specific functionality, this method is universally applicable to derive other functionalities and even different material classes while tracking hidden causal mechanisms via structural distortions.

42 ENGINEERING↗
DOE OSTI2022

Status of the NEAMS and ARC fast reactor tools integration to the NEAMS Workbench

The Workbench initiative was launched in FY-2017 within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to facilitate the transition from conventional tools to high-fidelity tools. The NEAMS Workbench provides a common user interface for model creation, real-time validation, execution, output processing, and visualization for integrated codes. The integration of the Argonne Reactor Computation (ARC) suite of codes into the NEAMS Workbench through the PyARC module was initiated in FY-2017. The ARC codes, which focus on fast reactor multiphysics analysis, contain both legacy codes like DIF3D and REBUS-3 that were developed with over 30 years of experience, and newer NEAMS additions like MC 2 -3, PERSENT, and PROTEUS. Recent work expended this integration to other tools to support the U.S. fast reactor community such as DASSH for sub-channel thermal-hydraulics, Griffin for high-fidelity deterministic transport calculations, and OpenMC for Monte Carlo simulations (with Shift integration planned for FY-2023). The integration of the “extended ARC” suite of codes into the NEAMS Workbench interface relies on the PyARC and PyGriffin modules to handles the pre- and post-processing of these codes input, and the runtime environment. The PyARC module together with the NEAMS Workbench interface are both released under Open Source Software licenses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Adsorption-based membranes for air separation using transition metal oxides

In this work, we use computational modeling to examine the viability of adsorption-based pore-flow membranes for separating gases when a purely size-based separation strategy is ineffective. Using molecular dynamics simulations of O 2 and N 2 , we model permeation through a nanoporous graphene membrane. Permeation is assumed to follow a five-step adsorption-based pathway, with desorption being the rate-limiting step. Using this model, we observe increased selectivity between O 2 and N 2 , resulting from increased adsorption energy differences. We explore the limits of this strategy, providing an initial set of constraints that need to be satisfied to allow for selectivity. Finally, we provide a preliminary exploration of some transition metal oxides that appear to satisfy those conditions. Using density functional theory calculations, we confirm that these oxides possess adsorption energies needed to operate as adsorption-based pore-flow membranes. These adsorption energies provide a suitable motivation to examine adsorption-based pore-flow membranes as a viable option for air separation.

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

Layer-Dependent Pressure Effect on the Electronic Structure of 2D Black Phosphorus

Through infrared spectroscopy, we systematically study the pressure effect on electronic structures of few-layer black phosphorus (BP) with layer number ranging from 2 to 13. We reveal that the pressure-induced shift of optical transitions exhibits strong layer dependence. In sharp contrast to the bulk counterpart which undergoes a semiconductor to semimetal transition under ~1.8 GPa, the band gap of 2 L increases with increasing pressure until beyond 2 GPa. Meanwhile, for a sample with a given layer number, the pressure-induced shift also differs for transitions with different indices. Through the tight-binding model in conjunction with a Morse potential for the interlayer coupling, this layer- and transition-index-dependent pressure effect can be fully accounted. Our study paves a way for versatile van der Waals engineering of two-dimensional BP.

2-dimensional systems↗
DOE OSTI2025

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Primordial black holes from bubble collisions during a first-order phase transition

We study the possibility of production of primordial black holes (PBHs) from bubble collisions during a first-order phase transition. While typical colliding bubbles are small and irrelevant for PBH production, we find that those that can produce PBHs must have a macroscopically thick fluid shell and have been born much before the typical nucleation time. Particularly large uncertainties arise from an exponential sensitivity of the nucleation rate on the required duration of bubble growth which depends on the details of the collisions and the evolution of the spacetime metric toward the end of the phase transition. We introduce a few parameters to be obtained from future numerical simulation to represent those unknowns, and estimate the PBH abundance in an Abelian Higgs benchmark model and show that it can be significant. We predict an approximately monochromatic PBH mass spectrum, and find regions in the parameter space where the PBHs can constitute entire dark matter or even overclose the Universe. Our result thus shows that models with a first-order phase transition can be constrained by overabundant PBHs or null results of other PBH searches. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Reduced-order condensed-phase kinetic models for polyethylene, polypropylene and polystyrene thermochemical recycling

Thermochemical recycling of plastic waste (PW) into chemicals and energy vectors requires coupling particle and reactor-scale simulations to accurate condensed phase pyrolysis mechanisms for each constituent. This work proposes a methodology to derive reduced-order condensed-phase kinetic models from validated semi-detailed kinetic mechanisms. Two types of kinetic models are obtained for polyethylene (PE), polypropylene (PP) and polystyrene (PS): reduced semi-detailed models and multi-step fully lumped ones. These families offer different compromises between accuracy and computational cost. The former employ 50–100 gas + liquid species and describe both the radical degradation and the detailed carbon distribution of the products. Conversely, the latter involves 5–10 species per polymer tracking only the main petroleum cuts. The kinetic mechanisms are complemented by the definition of thermochemical properties of gas, liquid, and solid-phase species, accounting for phase-transitions through pseudo-chemical reactions. Model validations are performed by comparison with experimental data and the original semi-detailed mechanisms in terms of mass loss, heat fluxes and product distribution profiles. The resulting CHEMKIN-like condensed-phase models are attached as Supplementary Material and as a GitHub repository. Extending the proposed approach to other polymers and coupling it with existing subsets in the CRECK kinetic framework (e.g., biomass, PVC, PET) offers a powerful tool to model thermochemical recycling of PW and biomass/PW mixtures.

kinetics↗
DOE OSTI2021

On the Order Parameter of the Continuous Phase Transition in the Classical and Quantum Mechanical Limits

The mean field theory is revisited in the classical and quantum mechanical limits. Taking into account the boundary conditions at the phase transition and the third law of the thermodynamics the physical properties of the ordered and disordered phases were reported. The equation for the order parameter predicts the occurrence of a saturation of Ψ 2 = 1 near Θ S , the temperature below the quantum mechanical ground state is reached. The theoretical predictions are also compared with high resolution thermal expansion data of SrTiO 3 monocrystalline samples and other some previous results. An excellent agreement has been found suggesting a universal behavior of the theoretical model to describe continuous structural phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Turbulence characteristics and flow dynamics impacts on the H-mode transition in favourable magnetic geometry with lower power threshold

We report the L-H transition power threshold (P LH ) in favourable magnetic geometry (ion ∇B drift pointing towards X-point) is much lower than in the unfavourable magnetic geometry (ion ∇B drift pointing away from X-point) on multiple tokamaks. In a systematic experiment on DIII-D, the ion ∇B drift direction was changed continuously from the unfavourable to favourable configuration during plasma discharges. During such process, the input neutral beam power was kept constant at a value that was above P LH for favourable configuration, but lower than P LH for unfavourable configuration. Toroidal field and plasma current were also kept constant and there was little change in the edge electron density ne and electron temperature T e profiles. The density fluctuation amplitude was reduced approaching the transition, while a large increase of turbulence Reynolds stress and flow shear were simultaneously observed. The turbulence decorrelation rate was found to increase as the ion ∇B drift direction was moving towards the favourable configuration, but the flow shear also increased and exceeded the turbulence decorrelation rate. These measurements demonstrate an important correlation between turbulence and turbulence-driven flow and a lowering of PLH, provide insights into the underlying physics behind the hidden parameters and inform a more complete physics-based model of the L-H transition power threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

3D Characterization of the Columnar-to-Equiaxed Transition in Additively Manufactured Inconel 718

Additive manufacturing (AM) provides enormous processing flexibility, enabling novel part geometries and optimized designs. Access to a local heat source further permits the potential for local microstructure control on the scale of individual melt pools, which can enable local control of part properties. In order to design tailored processing strategies for target microstructures, models predicting the columnar-to-equiaxed transition must be extended to the high solidification velocities and complex thermal histories present in AM. Here, we combine 3D characterization with advanced modeling techniques to develop a more complete understanding of the solidification process and evolution of microstructure during electron beam melting (EBM) of Inconel 718. Full calibration of existing microstructure prediction models demonstrates the differences between AM processes and more conventional welding techniques, underlying the need for accurate determination of key parameters that can only be measured directly in 3D. The ability to combine multisensor data in a consistent 3D framework via data fusion algorithms is essential to fully leverage these advanced characterization approaches. Thermal modeling provides insight on microstructure development within isolated solidification events and demonstrates the role of Marangoni effects on controlling solidification behavior.

Polonsky, Andrew↗
DOE OSTI2025

Precision study of the massive Schwinger model near quantum criticality

We perform a numerical analysis of the massive Schwinger model in the presence of a background electric field. Using the Density Matrix Renormalization Group approach, we efficiently compute the spectrum of the Schwinger model on a staggered lattice with up to 3000 qubits. As a result, we achieve a precise computation of the critical mass of the massive Schwinger model to five digits using four different “criticality criteria,” observing perfect agreement among them Additionally, we discuss the effect of a four-fermion operator deformation of the Schwinger model and compute the critical mass for various values of the deformation parameter.

Critical phenomena↗
DOE OSTI2025

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Non-analyticity in holographic complexity near critical points

Abstract The region near a critical point is studied using holographic models of second-order phase transitions. In a previous paper, we argued that the quantum circuit complexity of the vacuum ( C 0 ) is the largest at the critical point. When deforming away from the critical point by a term the complexity C ( τ ) has a piece non-analytic in τ , namely C 0 − C ( τ ) ∼ | τ − τ c | ν ( d − 1 ) + a n a l y t i c . Here, as usual, ν = 1 d − Δ and ξ is the correlation length ξ ∼ | τ − τ c | − ν and there are possible logarithmic corrections to this expression. That was derived using numerical results for the Bose–Hubbard model and general scaling considerations. In this paper, we show that the same is valid in the case of holographic complexity providing evidence that the results are universal, and at the same time providing evidence for holographic computations of complexity.

Physics↗