Fast algorithm for computing non-isothermal crystallization kinetics

A fast algorithm is presented for computing fractional crystallized volume of an amorphous solid under non-isothermal conditions, as a function of the number of discrete time intervals, temperature history and (temperature dependent) nucleation and growth rates. The algorithm is a modification of the discrete Yinnon–Uhlmann approach to compute the standard double integral formula under quasi-steady-state conditions. Rather than re-computing the kinetics over the entire thermal history of all previous time intervals for each new interval, the crystallized fraction for a given time interval is computed by a sum of terms involving only the current and previous intervals, such that orders of magnitude increase in computational speed is obtained. •A fast numerical method is derived for computing non-isothermal crystallization for an arbitrary thermal history.•Derived from the standard JMAK double integral equation•Does not require re-computing the entire thermal history at each step•Accuracy and computational efficiency is demonstrated using examples.•Possible extensions to include non-QSS (quasi-steady state) are discussed.