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); Levitas, Preston and Lee (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 $α\toε$ iron phase transition to experimental data of Smith et al. (2013). 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 characterization data for the microstructure of the sample under consideration.

💡 Research Summary

The paper presents a comprehensive kinetic framework for solid‑solid phase transformations in metals, extending the Levitas‑Preston (LP) phase‑field model—originally limited to near‑zero pressure—to arbitrary pressure conditions. The authors construct a Gibbs free‑energy functional G(σ,P,T,η) that depends on the stress tensor, pressure, temperature, and a scalar order parameter η (η = 0 for the parent phase A, η = 1 for the product phase M). By expressing the transformation strain ε_t as a sum of deviatoric and volumetric components and incorporating pressure‑dependent bulk (B) and shear (μ) moduli, they derive an explicit expression for the transformation work W_λ, which combines mechanical work and elastic energy contributions. This work term, together with the pressure‑ and temperature‑dependent free‑energy difference ΔG(P,T), defines the thermodynamic driving force for the A↔M transition.

Three nucleation pathways are modeled: homogeneous nucleation, nucleation at grain boundaries, and nucleation on dislocations. For each, the critical nucleus radius and nucleation barrier are obtained from a three‑dimensional spherical nucleus model, and nucleation rates are expressed in Arrhenius form I = I₀ exp(–ΔG*/k_BT). The dislocation‑mediated nucleation rate explicitly includes the dislocation density ρ_d, highlighting the importance of microstructural characterization.

The evolution of the transformed volume fraction X(t) is treated with the Kolmogorov‑Johnson‑Mehl‑Avrami (KJMA) theory. The authors couple the nucleation rate I with an interface propagation speed v_int, which they model as v_int = M (W_λ – ΔG), where M is a mobility parameter. This yields the classic KJMA differential equation dX/dt = 4π v_int r² I (1 – X), which is integrated to obtain the full time‑dependent phase fraction. While the current formulation handles only two phases, the authors discuss straightforward extensions to multi‑phase systems.

To validate the model, the authors simulate the α (bcc) → ε (hcp) transition in iron under ramp loading, where pressure increases linearly with time (P = α t). Using an in‑house Python code, they compare simulated pressure‑versus‑fraction curves with experimental data from Smith et al. (2013). Parameter studies reveal that multiple combinations of dislocation density and interface mobility can reproduce the experimental observations, underscoring the degeneracy of model parameters and the need for accurate microstructural data. The model successfully captures the “overshoot” pressure at which the transition initiates and the total transformation time, both of which shift with loading rate.

In conclusion, the work delivers a unified, physics‑based kinetic model that integrates pressure‑dependent phase‑field thermodynamics with classical nucleation theory and KJMA kinetics. It provides a tool for predicting transformation dynamics under high‑pressure, high‑rate conditions relevant to shock compression, dynamic material synthesis, and high‑energy density physics. Future directions include refining the treatment of dislocation‑grain boundary interactions, extending to reverse (ε → α) transitions, and incorporating shock‑wave propagation effects for a fully coupled dynamic simulation framework.