kinetic-monte-carlo-method-1759596534

Kinetic Monte Carlo method

The Kinetic Monte Carlo (kMC) method is a stochastic simulation approach that explicitly models atomic or microstructural events such as nucleation, growth, and coarsening using probabilistic rules. Unlike the JMAK model, which provides an analytical expression for transformed fraction, kMC simulates individual events over discrete time steps to capture detailed microstructural evolution.

For a comprehensive overview of kMC in the context of Finite Element Analysis, see the kMC-FEA introduction. The kMC method will be compared against the JMAK model as a benchmark for solid-state phase transformation modelling.

Event-Driven Kinetics

In kMC, the evolution of the system is governed by the rates of possible events (e.g., nucleation, growth, dissolution), which are translated into probabilities used to randomly select events at each time step.

Event Rate

Each possible event $i$ has an associated rate $r_i$ determined by Arrhenius-type kinetics:

$$r_i={\nu }_i\exp \left(-\frac{E_i}{k_{B}T}\right)$$

where

• ${\nu }_i$ is the attempt frequency (pre-exponential factor),
• $E_i$ is the activation energy for event $i$,
• $k_B$ is the Boltzmann constant,
• $T$ is the absolute temperature.

Total Rate and Selection Probability

The total rate is the sum over all possible events:

$$R=\sum _i{r}_i$$

Each event’s selection probability $p_i$ is:

$$p_i=\frac{r_i}{R}$$

Time Increment

The kMC time step $\Delta t$ advances stochastically according to:

$$\Delta t=-\frac{\ln (u)}{R}$$

where $u$ is a uniform random number in $(0,1)$.

Coupling with FEA

Within each finite element, local conditions (e.g., temperature, stress) influence $E_i$ and ${\nu }_i$. The kMC algorithm at each time step:

1. Calculate event rates $r_i$ based on local conditions.
2. Compute total rate $R$ and select an event according to $p_i$.
3. Update microstructural state variables (e.g., phase fraction, grain size) accordingly.
4. Increment simulation time by $\Delta t$.
5. Pass updated microstructure information back to the FEA solver to update material properties.

This approach enables multiscale coupling, capturing microscale stochasticity within a macroscale continuum framework. The single element transient thermal model built with deal.II provides the foundation for this coupling, while the Docker build setup ensures a consistent compilation environment for the deal.II codebase. The kMC method offers finer microstructural resolution than the JMAK model but at significantly higher computational cost.

Advantages

• Captures detailed, spatially-resolved microstructural evolution.
• Can model complex kinetics and heterogeneous nucleation sites.
• Incorporates microstructural effects such as grain boundaries and anisotropy more naturally.
• Allows direct simulation of discrete events rather than relying on averaged assumptions.

Disadvantages

• Computationally intensive, especially for large-scale FEA domains.
• Requires careful parameterization and validation of event rates and probabilities.
• Less straightforward to couple with macroscopic models due to scale differences.
• Longer simulation times can limit applicability for very large problems.