Search NASA⌕ Search

SEARCH · Search NASA

Results for “phase kinetics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hydro-Code Implementation and Testing of a Kinetic Phase Transition Framework

In this report we describe the Kinetic Phase Transition (KPT) framework that has been worked out over the last 10 years (from around 2014) and the implementation of it into three different codes, the one-dimensional hydro- LASLO and the three-dimensional magneto-hydro- ALEGRA, Sandia codes, via subroutines in the LAMBDA Equations of State and constitutive models package, and Flag, an arbitrary Lagrangian-Eulerian multiphysics code developed within the Lagrangian Applications project (LAP) at LANL. We discuss the introduction of phase mass (and/or volume) fractions that are needed in a code for it to be ‘phase aware’, that is, not only the thermodynamic state is known in each point but also the mixture of the materials’ phases in that point. Further we point to the need of a full Equations of State for each phase in a material to achieve phase awareness and we review the equilibrium phase model, where a phase mixture is at its lowest Gibbs free energy state, to make this point clear. Contrasting the kinetic phase transition to this equilibrium model seamlessly introduce us to the KPT framework that is subsequently thoroughly discussed. While the determination of the total state and the states and mass fractions of phases in each point is a problem that can borrow many of its numerical details from Eulerian codes and mixture of materials (not phases), the update of mass fractions with time in a KPT framework needs a new set of considerations. General for any update model is that we need to prevent mass fractions from becoming unphysical (negative or their sum to be larger than one). We have solved this problem by implementing a subdivision of the hydro time step that prevents the phase from being fully present to not present at all in one subdivided time step by limiting the size of the subdivided time step. This scheme also corrects numerical problems from abrupt changes in parameter values, the so called Gibbs phenomena, that gives rise to slushing between phases in the KPT framework. Interspersed throughout the report are discussions on different thermodynamics considerations. EOS validity windows, limitations on the EOS phase space, are needed for the KPT framework and are discussed separately and exemplified. The KPT framework described in this report has been verified by code comparison, but validation is still an active area of research. There is room for improvement in the update model, both in the model for determination of rates and in how to prevent the mass fractions from becoming unphysical. In addition, the parameters in the KPT update model and the placement of the phase boundary in the EOS phase space, and interactions with other constitutive models, are closely related and interfering with each other. One possible way forward is to simultaneously develop KPT parameters, EOS, and constitutive models for each material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Parameter study of a Kinetic Phase Transition model

We investigate the ability of the two model parameters in the Kinetic Phase Transition (KPT) model of Carl Greeff to describe phase transitions of various types when combined with the tin Equation of states for the phases, SESAME 2162. While the two parameters are correlated, one combination can be associated with the nucleation process giving rise to the hysteresis seen in phase transition paths in pressure-temperature space and another can be associated with the overall timescale of the transition. We acknowledge the desire to have a model where the parameters are more independent of each other for easier model parameter determination.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2025

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗
DOE OSTI2025

A Semi-Detailed Pyrolytic Gas-Phase Kinetic Model for the Volatiles of Polyethylene Thermal Degradation

This work presents a semi-detailed kinetic model to address the pyrolytic gas-phase reactivity of volatiles formed during thermal degradation of polyethylene (PE). The model builds on a validated multi-step condensed-phase model and employs validated lumping approaches. Short-chain compounds are modelled with high detail, while long-chain ones are described by surrogate species representative of diesel-cuts (NC16H32) and waxes (NC30H60). The reactivity of short chains is described through the comprehensive CRECK kinetic model, updated to align C5-C7 olefins based on recent literature experimental data. Due to the lack of experimental data for longer olefins, their reactivity is modeled by analogy to the shorter ones, ensuring an asymptotic behavior with increasing carbon numbers. The semi-detailed model is validated through experimental data on PE pyrolysis, assuming an instantaneous mixing of the inert inlet flow with released volatiles, followed by a segregated plug-flow behavior. Validation across different reactor setups confirms the model’s capability to predict detailed product distributions. Despite minor discrepancies, the proposed model effectively captures experimental trends. Further work will address modelling the reactivity in oxygen-containing environments.

kinetics↗
DOE OSTI2022

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phase transition kinetics of superionic H2O ice phases revealed by Megahertz X-ray free-electron laser-heating experiments

Abstract H 2 O transforms to two forms of superionic (SI) ice at high pressures and temperatures, which contain highly mobile protons within a solid oxygen sublattice. Yet the stability field of both phases remains debated. Here, we present the results of an ultrafast X-ray heating study utilizing MHz pulse trains produced by the European X-ray Free Electron Laser to create high temperature states of H 2 O, which were probed using X-ray diffraction during dynamic cooling. We confirm an isostructural transition during heating in the 26-69 GPa range, consistent with the formation of SI-bcc. In contrast to prior work, SI-fcc was observed exclusively above ~50 GPa, despite evidence of melting at lower pressures. The absence of SI-fcc in lower pressure runs is attributed to short heating timescales and the pressure-temperature path induced by the pump-probe heating scheme in which H 2 O was heated above its melting temperature before the observation of quenched crystalline states, based on the earlier theoretical prediction that SI-bcc nucleates more readily from the fluid than SI-fcc. Our results may have implications for the stability of SI phases in ice-rich planets, for example during dynamic freezing, where the preferential crystallization of SI-bcc may result in distinct physical properties across mantle ice layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Assessment of Physics Models for Phase Transition Kinetics

The time dependence of phase diagrams and how to model rate dependent transitions remains one of the key unanswered questions in physics. When a material is loaded dynamically through equilibrium phase boundaries, it is the kinetics that determines the real time expression of a phase transition. Here we report the atomic and nanosecond-scale quantification of kinetics of shock-driven phase transition in multiple materials. We uniquely make use of a both a simple shock as well as shock-and-hold loading pathways compress different crystalline solids and induce structural phase transitions below melt. Coupling shock loading with time-resolved synchrotron x-ray diffraction (DXRD), we probe the structural transformations of these solids in the short-lived high pressure and temperature states generated. The novelty and power of using DXRD for the assessment of kinetics of phase transitions lies in the ability to discover and identify new phases and to examine kinetics without prior knowledge of a material's phase diagram. Our results provide a quantified expression and a physics model of kinetics of formation of high-pressure phases under shock loading: transition incubation time, evolution, completion time and crystallization rate.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Elucidating the Microscale Behavior and Phase Separation Kinetics of Thermally Responsive Ionic Liquid–Water Mixtures

Thermally responsive ionic liquids (ILs) exhibit liquid-liquid phase separation into a water-rich (WR) and ionic-liquid-rich (ILR) phase when heated above a lower critical solution temperature (LCST). This phase behavior has been leveraged for applications ranging from forward osmosis (FO) desalination, where the IL acts as a draw solute, to refrigeration and dehumidification cycles, where the IL acts as a liquid desiccant. While significant effort has been devoted to characterizing the thermodynamic and thermophysical properties of LCST ILs, their phase separation kinetics have not been investigated. In this work, we describe the macroscale phase separation kinetics (phase separation time) by gleaning insight into the microscale colloidal behavior of aqueous mixtures of four different materials, P 4444 TFA (tetrabutylphosphonium-2,4-trifluoroacetate), P 4444 DMBS (tetrabutylphosphonium-2,4-dimethyl-benzenesulfonate), N 4444 Sal (tetrabutylammonium salicylate), and P 4444 Sal (tetrabutylphosphonium salicylate) as a function of IL concentration at a separation temperature of 70 °C. We report the discontinuous microscale size distributions for each material and correlate their theoretical settling velocities to experimental phase separation times. The results indicate that a simple Stokes' law model can predict the phase separation time within reasonable accuracy. Overall, this work lays the foundation for understanding the micro- to macroscale phase separation behavior and kinetics of LCST ILs for various water-energy applications.

LCST↗
DOE OSTI2025

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Thickness-dependent phase transition kinetics in lithium-intercalated MoS 2

The phase transitions of two-dimensional (2D) materials are key to the operation of many devices with applications including energy storage and low power electronics. Nanoscale confinement in the form of reduced thickness can modulate the phase transitions of 2D materials both in their thermodynamics and kinetics. Here, using in situ Raman spectroscopy we demonstrate that reducing the thickness of MoS 2 below five layers slows the kinetics of the phase transition from 2H- to 1T'-MoS 2 induced by the electrochemical intercalation of lithium. In this study, we observe that the growth rate of 1T' domains is suppressed in thin MoS 2 supported by SiO 2 , and attribute this growth suppression to increased interfacial effects as the thickness is reduced below 5 nm. In this study, the suppressed kinetics can be reversed by placing MoS 2 on a 2D hexagonal boron nitride ( h BN) support, which readily facilitates the release of strain induced by the phase transition. Additionally, we show that the irreversible conversion of intercalated 1T'-MoS 2 into Li2S and Mo is also thickness-dependent and the stability of 1T'-MoS 2 is significantly increased below five layers, requiring a much higher applied electrochemical potential to break down 1T'-MoS 2 into Li 2 S and Mo nanoclusters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Inference of strength and phase transition kinetics in dynamically-compressed tin

Dynamic compression experiments in condensed matter are of interest in part because they provide opportunities to examine material response under extreme conditions; however, the inference of material behavior from dynamic experiments is challenging in the presence of phase transitions exhibiting kinetic processes. Demonstrating an approach to quantitative interpretation of such dynamic experiments, here we present a Bayesian model calibration of strength and phase transformation parameters to data drawn from pulsed power and gas gun shot experiments. The posterior predictions of the Bayesian model capture the experimental measurements and account for the various uncertainties in the experimental configurations. This holistic approach to model calibration utilizing multiple types of experimental data identifies important cross correlations among kinetics, strength, and the phase boundary. Improved insight into potential sources of current model form error is provided by comparing the differences between calibrations against different subsets of the experimental data.

42 ENGINEERING↗
DOE OSTI2025

Millisecond Phase Transition Kinetics of Lyotropic Liquid Crystalline Nanoparticles Observed by Time‐Resolved Small Angle X‐ray Solution Scattering

This study investigates the dynamic behavior of lyotropic liquid crystal nanoparticles (LCNPs), which are widely recognized for their applications in drug delivery. By employing nanosecond near‐infrared laser pulse‐induced temperature jump (T‐jump) and time‐resolved X‐ray solution scattering, the structural dynamics of phase transitions in phytantriol‐based cubosomes and hexosomes are revealed. Both cubosome and hexosome LCNPs undergo phase transitions into noncrystalline phases at high temperatures. Their phase transition kinetics, occurring within milliseconds (ms) and involving one intermediate structure, are captured. Additionally, the reverse self‐assembly processes of LCNPs were observed, occurring on the timescale of a few hundred ms. To our knowledge, this is the first observation of LCNP T‐jump induced phase transitions on the ms timescale and their reverse self‐assembly. These findings provide valuable insights into the LCNP phase transition processes, with potential implications for drug delivery applications.

Chemistry↗
DOE OSTI2023

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Comparative study of crystallization kinetics and phase segregation of triple cation and methylammonium lead iodide perovskites on moisture probing using synchrotron X-ray based radiation

3D mixed perovskites have achieved substantial success in boosting solar cell efficiency, but the complicated perovskite crystal formation pathway remains mysterious. Here we present detailed crystallization kinetics of mixed perovskites FA 0.83 MA 0.17 Pb(I 0.83 Br 0.17 ) 3 , where FA is formamidinium and MA is methylammonium, with the addition of Cs + to form a triple cation perovskite (3-CAT), in a comparison with the perovskite building block MAPbI 3 (MAPI) via static grazing-incidence wide-angle X-ray scattering (GIWAXS) and micro-diffraction measurements. Spin-coated films produced α-perovskite peaks with no PbI 2 or δ-intermediate phases, which was a promising result for the 3-CAT perovskite from micro-diffraction measurements. However, the 3-CAT did not remain stable on probing with varied relative humidity (RH) conditions as segregation back to the δ-intermediate and PbI 2 phase after 10 s of exposure to an RH value of 11% was found to occur from the GIWAXS results. When RH levels were elevated to over 100%, segregation peaks of PbI 2 and δ-intermediate (2H, 4H and 6H) became conspicuous as the α-phase intensity diminished, unlike for MAPI that remains relatively stable. The possible cause of this is hydrophilic bonds that form between the 3-CAT crystals and the small annealing window of the best composition perovskite (5% Cs + ) film.

Shatsala, Miller (ORCID:0009000927230192)↗