A solution to the problem of freezing of a poroelastic material is derived and analysed in the case of one-dimensional deformation. The solution is sought within the framework of thermo-poroelasticity, with specific account of the behaviour of freezing materials. The governing equations of the problem can be combined into a pair of coupled partial differential equations for the temperature and the fluid pressure, with particular forms in the freezing and the unfrozen regions. In the freezing region, the equations are highly non-linear, partly due to the dependence of thermal and hydraulic properties on water saturation, which varies with temperature. Consequently, the solution is obtained through numerical methods, with special attention to the propagation of the freezing front boundary. The response to one-dimensional freezing is illustrated for the case of cement paste. Finally, the influence on the solution of varying selected parameters is analysed, such as the temperature boundary conditions, the parameters characterizing the geometry of the porous system, the ratio of fluid and thermal diffusivities, and the rate of cooling applied at the freezing end. Copyright (C) 2008 John Wiley & Sons, Ltd.
|Journal||International Journal for Numerical and Analytical Methods in Geomechanics|
|Publication status||Published - 2008|