Frost heave and phase front instability in freezing soils

The main equations and conditions at the phase transition front are presented for a generalized model of secondary frost heave in freezing fine-grained soils. The analytical criterion for the stability/instability of the freezing phase front in porous media is derived. This criterion is obtained for the occurrence of the frost heave process by using the perturbation method in a two-dimensional, coupled heat and mass transfer model. This model assumes that the non-instantaneous crystallization process takes place in the kinetic zone, and that the rate of crystallization is a function of supercooling. This corresponds to the Arrhenius form equation and agrees with experimental investigations. The perturbation analysis of the freezing front shows that the stability criterion depends upon 1) the Stefan and Peclet numbers, 2) a parameter describing the phase transition kinetics and also 3) dimensionless parameters which characterize the frost heave process. Employing Fourier synthesis, actual front shape evolution is calculated. It is seen that the front displays a periodic morphology whose scale is essentially unrelated to that of the initial (starting) perturbation. The effect of the non-instantaneous kinetics on the front shape evolution is described. As is shown in results, the kinetics has a stabilizing effect and, in this case, the perturbations grow more slowly. The theoretical stability/instability conditions as predicted from the derived criterion were found to be in agreement with experimental investigations of the formation of soil cryogenic structure in the freezing process. On the basis of the asymptotic solution the engineering approach for the calculation of the heave rate and maximal frost penetration depth values — main characteristics for design and construction in cold regions, is presented. The good agreement between calculated values and experimental data is observed.