Online multiscale finite element simulation of thermo-mechanical model with phase change

This paper presents a thermo-mechanical model with phase transition considering changes in the mechanical properties of the medium. The proposed thermo-mechanical model is described by a system of partial differential equations for temperature and displacements. In the model, soil deformations occur due to porosity growth caused by ice and water density differences. A finite element approximation of this model on a fine grid is presented. The linearization from the previous time step is used to handle the nonlinearity of the problem. For reducing the size of the discrete problem, offline and online multiscale approaches based on the Generalized Multiscale Finite Element Method (GMsFEM) are proposed. A two-dimensional model problem simulating the heaving process of heterogeneous soil with a stiff inclusion was considered for testing the mathematical model and the multiscale approaches. Numerical solutions depict the process of soil heaving caused by changes in porosity due to the phase transition. The movement of the phase transition interface was observed. The change of medium properties, including the elastic modulus, was traced and corresponds to the phase transition interface. The proposed multiscale approaches significantly reduce the size of the discrete problem while maintaining reasonable accuracy. However, the online multiscale approach achieves better accuracy than the offline approach with fewer degrees of freedom.