Numerical Study of Soil-Thawing Effect of Composite Piles Using GMsFEM

: During construction works, it is advisable to prevent strong thawing and an increase in the moisture content of the foundations of engineering structures in the summer. Since the density of water and ice differ, due to the difference bulging of the foundation sections can occur when it freezes back in winter. In this work, the effect of ﬁber-reinforced piles on the thermal ﬁeld of the surrounding soil is investigated numerically; that is, the study of the inﬂuence of aggregates with high and low thermal-physical properties on the temperature of frozen soils is conducted. Basalt and steel ﬁber reinforcement are compared. The difﬁculty of this work is that the inclusions inside piles are too small compared to the pile itself. Therefore, to solve the Stefan problem, a generalized multiscale ﬁnite element method (GMsFEM) was used. In the GMsFEM, the usual conforming partition of the domain into a coarse grid was used. It allowed reducing problem size and, consequently, accelerating the calculations. Results of the multiscale solution were compared with ﬁne-scale solution, the accuracy of GMsFEM was investigated, and the optimal solution parameters were deﬁned. Therefore, GMsFEM was shown to be well suited for the designated task. Collation of basalt and steel ﬁber reinforcement showed a beneﬁcial effect of high thermal conductive material inclusion on freezing of piles in winter. draft preparation, visualization, P.V.S.;


Introduction
Thermal calculations are important in the construction of engineering geotechnical structures and buildings in the permafrost zone. The temperature regime (a set of sequential temperature fields in the soil mass corresponding to any given points in time from the beginning of the calculation) is calculated as the forecast of the thermal effects on the upper and lower boundaries of the structure foundation set for the entire calculation period [1]. The most characteristic feature of these processes is the previously unknown ("free") boundaries between various thawed and frozen states of the soil. Due to this feature, their mathematical models are non-linear and difficult to analyze. The results of calculations of climatic phenomena, which are critical for geotechnical structures on permafrost foundations (in the permafrost zone), show a continuous repetition of the same constant seasonal cycle. This determines periodic harmonic oscillations of all quantities that determine the state of the foundation of a geotechnical structure, including temperature. However, cyclic harmonic changes in the foundation state caused by constantly repeating influences also occur during other processes. The course of temperature change during one period can have a different character. For example, the temperature can change abruptly, continuously increase or decrease, etc.
During construction works, it is advisable to prevent strong thawing and an increase in the moisture content of the foundations of engineering structures in the summer. Since the density of water and ice differ, due to the difference bulging of the foundation sections can occur when it freezes back in winter. Repetition of this phenomenon for many years causes the destruction of engineering structures and buildings. In the regions of the Far North, it leads to catastrophic emergencies.
In this work, the effect of fiber-reinforced piles on the thermal field of the surrounding soil is numerically investigated [2], that is, the study of the influence of aggregates with high and low thermal-physical properties on the temperature of frozen soils is conducted [3,4].
Basalt fibers are produced from molten basalt rock. They have very low thermal conductivity, good strength properties and, at the same time, are relatively cheap. Due to these characteristics, they are used in concrete [5][6][7]. However, in addition to the chemical and mechanical properties of basalt fibers, their cost varies, depending on the type and quality of the raw material and the production process of these fibers [8]. Nevertheless, the characteristics mentioned above and an environmentally friendly manufacturing process [9] might determine their application in high performance concrete structures instead of the most commonly used steel and polypropylene fibers. The length-to-diameter aspect for basalt fibers is about 1000.
As for steel fiber hybridization, industrial or recycled steel fibers can be used [10][11][12][13]. According to results available in the scientific literature [14,15], they show similar mechanical responses, both in terms of tensile strength and matrix-to-fiber bond. Recycled steel fibers are derived from waste tires. Their geometrical characterization can be highly variable: they are generally characterized by a nominal diameter ranging between 0.1 and 2 mm with a corresponding average aspect ratio (i.e., length-to-diameter ratio) ranging between 20 and 150. These variations mainly depend on both the original source (i.e., tires typology) and recycling processes.
The volume concentration of fibers is equal to 10%. In our article, we use the lengthto-diameter ratio for both types of fibers equal to 32 for structured and 10 for random distribution of fibers for calculation convenience. Fibers' mechanical properties are presented in Table 1. One of the topical problems of mathematical physics is the Stefan problem, which is an initial-boundary value problem for a parabolic differential equation with discontinuous coefficients, and serves as a mathematical model of the change in the phase state of a substance with unknown interfaces [16,17]. The phase transition boundary is represented as a smearing zone when the phase transition occurs in a given temperature range [16,18].
The difficulty of this work lies in the pile inclusions that are too small compared to the solution area. For such geometrically complex problems, reduction techniques such as multiscale methods are preferred [19,20]. Multiscale methods are becoming very popular at this time [21][22][23]. Paper [24] develops a multiscale eigenvalue method based on multiscale substructure technology for multiscale analysis of periodic composite structures. This provides a comparative study from a user's viewpoint for multiscale methods, including the mathematical homogenization method, heterogeneous multiscale method, and multiscale finite element method. The article [25] studies the multiscale finite element method for solving elliptic problems arising from composite materials. The influence of production porosity, diameter [26], and density of structures on the thermal conductivity of a composite are considered using multiscale finite element modeling [27]. In [28,29], mathematical modeling is applied for the thermodynamic analysis and design of composite structures. This article is organized as follows. Section 2 presents the formulation of a mathematical model that describes the dynamics of the temperature distribution with the phase transitions taken into account. In the third, finite element approximation is given. In the fourth, a generalized multiscale finite element method (GMsFEM) is considered. Section 5 presents the numerical results in a two-dimensional setting.

Mathematical Model
Let us consider a mathematical model describing the dynamics of the temperature distribution, taking into account the phase transitions of pore moisture into ice and back, at a certain given temperature of the phase transition T * in area Ω = Ω − ∪ Ω + . Where Ω + (t) = {x|x ∈ Ω, T(x, t) > T} is the area occupied by the thawed soil, where the temperature exceeds the phase transition temperature and Ω − (t) = {x|x ∈ Ω, T(x, t) < T} is an area occupied by the frozen ground. A phase transition occurs at the interface between thawed and frozen soils S = S(t) (Figure 1). ment method for solving elliptic problems arising from composite materials. The influence of production porosity, diameter [26], and density of structures on the thermal conductivity of a composite are considered using multiscale finite element modeling [27]. In [28,29], mathematical modeling is applied for the thermodynamic analysis and design of composite structures.
This article is organized as follows. Section 2 presents the formulation of a mathematical model that describes the dynamics of the temperature distribution with the phase transitions taken into account. In the third, finite element approximation is given. In the fourth, a generalized multiscale finite element method (GMsFEM) is considered. Section 5 presents the numerical results in a two-dimensional setting.

Mathematical Model
Let us consider a mathematical model describing the dynamics of the temperature distribution, taking into account the phase transitions of pore moisture into ice and back, at a certain given temperature of the phase transition * in area Ω = Ω ∪ Ω . Where is the area occupied by the thawed soil, where the temperature exceeds the phase transition temperature and Ω ( ) = | ∈ Ω, ( , ) is an area occupied by the frozen ground. A phase transition occurs at the interface between thawed and frozen soils = ( ) ( Figure 1). To simulate heat transfer processes with phase transitions, the Stefan model is used. It describes thermal processes accompanying phase transformations of the medium with absorption and release of latent heat: where is the specific heat of phase transition, m is a porosity, ( − * ) is Dirac delta function.
For the coefficients of the equation, there are following relations: where , , and , , are the specific heat, density, and thermal conductivity of thawed and frozen soil, respectively.
Since the process of heat propagation is considered in a saturated porous medium, there are following thermal-physical characteristics: To simulate heat transfer processes with phase transitions, the Stefan model is used. It describes thermal processes accompanying phase transformations of the medium with absorption and release of latent heat: where L is the specific heat of phase transition, m is a porosity, δ(T − T * ) is Dirac delta function. For the coefficients of the equation, there are following relations: where c + , ρ + , λ + and c − , ρ − , λ − are the specific heat, density, and thermal conductivity of thawed and frozen soil, respectively. Since the process of heat propagation is considered in a saturated porous medium, there are following thermal-physical characteristics: where m is a porosity. The subscripts sc, w, i denote the skeleton of the porous medium, water, and ice, respectively. For the coefficients of the thermal conductivity in the thawed and frozen zones, there are similar relations: In practice, phase transformations do not take place instantly and can proceed in a small temperature range [T − ∆, T + ∆]. The discontinuous coefficients of Equation (1) are replaced by sufficiently smooth functions of temperature: Thus, the following equation for the temperature is obtained: Thus, Equation (2) is a multidimensional quasilinear parabolic equation with smooth coefficients.
For piles and fibers the following is true: (cρ) ∆ = (cρ) p, f and λ ∆ = λ p,f , where c p, f , ρ p, f , λ p,f are specific heat, density, and thermal conductivity of piles and fibers. That is, there is no phase transition in the pile and fiber areas, respectively.
Equation (2) is supplemented with the initial and boundary conditions Here η-normal vector.

Fine-Scale Approximation
The quasilinear parabolic Equation (2) with the corresponding boundary and initial conditions is approximated using the finite element method in combination with a purely implicit linearized finite difference approximation in time. This means that during discretization, the coefficients depending on the desired function are taken from the previous time layer. Let us write down the variational problem statement for each time layer: find T ∈ H 1 such that: To solve the problem numerically, it is necessary to pass from a continuous variational problem to a discrete one. For this, finite-dimensional spaces V h ∈ H 1 ,V h ∈ H 1 0 are introduced and the following problem is defined on them: find T h ∈ V h such that: τ Ω (cρ) σ T n h w h dx + Γt αT air ds.

Generalized Multiscale Finite Element Method (GMsFEM)
Construction local reduction of a model on the snapshot space is described by solving several local spectral problems using GMsFEM. In the GMsFEM the usual conforming partition of domain into finite elements is used. This partition is called the coarse grid T H . The coarse grid is partitioned into coarse-grid blocks. That is, T H is a coarse grid in domain where h is size of fine grid. Each coarse-grid block (K i ) consists of the connected union of fine-grid blocks. Furthermore, a certain domain for this is constructed and denoted by N c the coarse nodes number, by {x i } N c i=1 the vertices of the coarse mesh and define the neighborhood of the node x i : The following steps need to be implemented for GMsFEM realization: 1.
Coarse grid generation T H ; 2.
Construction of snapshot space that will be used to compute an offline space; 4.
Construction of a small dimensional offline space by performing dimensional reduction in the space of local snapshots; 5.
Solution of a coarse-grid problem for any force term and boundary condition.
In the first step of GMsFEM (Offline stage), the "snapshots" space must be constructed, a large dimensional snapshots space of local solutions. In the "snapshots" space, the following is considered: where k α are the coefficients of the thermal conductivity, δ j (x) are certain set of function defined on ∂ω i , here j = 1, J ω . The J ω is a number of fine grid edges on ω. Therefore, following is defined: This allow reducing the snapshot space to offline space via some spectral procedure. Offline space is constructed using the following local spectral problems in the snapshots space: To generate the offline space, the smallest M where M After this procedure, conforming basis functions can be obtained in the space: Further, the projection matrix R T = ϕ ω i 1 , . . . , ϕ Using constructed multiscale space, the coarse-scale system is solved: After solving the coarse-scale solution, a solution on the fine grid T ms = R T T c can be calculated.

Numerical Results
In this section, the effect of fibers in a pile made of different materials on the temperature field of the soil is calculated. Thus, the melting effect of the pile is numerically simulated.
Conditions for the impact of piles on the temperature regime of the soil: • Cement-sand mortar for filling the sinuses between the soil and the pile with a temperature of +20 • C; Soil temperature is 1.5 °C. Thermal-physical characteristics are taken from SP 25.13330.2012 and presented in Table 2. Calculations were performed for 365 days with a time step = 1 day (24 h).
The temperature of the daylight surface of the structure base was set considering the amplitude of the air temperature fluctuations taken from the data of the Yakutsk meteorological station. Figure 3 demonstrates the changes in the air temperature. It shows the distribution of the air temperature during one year.

1.
Fibers from basalt arranged in a structured manner; 2.
Fibers made of steel arranged in a structured manner; 3.
Fibers made of randomly located steel.
Soil temperature is 1.5 • C. Thermal-physical characteristics are taken from SP 25.13330.2012 and presented in Table 2. Calculations were performed for 365 days with a time step τ = 1 day (24 h). The temperature of the daylight surface of the structure base was set considering the amplitude of the air temperature fluctuations taken from the data of the Yakutsk meteorological station. Figure 3 demonstrates the changes in the air temperature. It shows the distribution of the air temperature during one year.  Table 2. Calculations were performed for 365 days with a time step = 1 day (24 h).
The temperature of the daylight surface of the structure base was set considering the amplitude of the air temperature fluctuations taken from the data of the Yakutsk meteorological station. Figure 3 demonstrates the changes in the air temperature. It shows the distribution of the air temperature during one year.  Here are the results of numerical calculations. Figures 4 and 5 show the soil temperature distributions for different points in time. These results were obtained by the GMsFEM method for eight basis functions. Here are the results of numerical calculations. Figures 4 and 5 show the soil temperature distributions for different points in time. These results were obtained by the GMsFEM method for eight basis functions.    Here are the results of numerical calculations. Figures 4 and 5 show the soil temperature distributions for different points in time. These results were obtained by the GMsFEM method for eight basis functions.   For numerical comparison of the fine-scale and multiscale solutions, relative L 2 and energy errors are used.
where T h and T ms are the fine-scale and multiscale solutions. Tables 3 and 4 show the relative errors of L 2 and energies for different number of multiscale basis functions and corresponding degrees of freedom (DOF). Degrees of freedom for fine scale mesh is denoted as DOF f . Here are the uncertainties for fibers arranged in structured order. When using basalt fibers, you can limit yourself to using 2 multiscale basis functions, and when using steel fibers, you need to use at least 4 multiscale basis functions.  Tables 5 and 6 show the relative errors of L 2 and energies for different number of multiscale basis functions for case, when fibers are located randomly. It can be seen that using four multiscale basis functions leads to good accuracy. A comparison was also made for different fiber arrangements (Figures 6 and 7). As the figures show, random arrangement of the fibers is better. However, this may be due to the fact that they are located horizontally.

Discussion
For the solution of the Stefan problem related to the soil thawing effect of composite piles, GMsFEM was used. According to results showing the accuracy of the method, the usage of four basis functions leads to the results with error in L2 norm lower than 1%. The errors in energy norms are also very respectable. Thus, usage of GMsFEM for such problems is recommended.
Results show that the process of total freezing of piles with steel fiber inclusion in winter runs faster. Analysis shows that the length between isoclines of phase change is about 47 cm for structured and 37 for random distribution of fibers. Freezing of pile with steel fibers comes on the 239th day, and a pile with basalt inclusion freezes on the 250th day in the structured case. For random distribution, similar days are 245th and 252nd. Therefore, there are 11 and 7 days' delay between freezing of piles with steel and fiber inclusions, depending on the distribution.

Conclusions
In this work, the effect of a composite pile on the temperature field of the surrounding frozen soil using GMsFEM is considered. Results of calculations using

Discussion
For the solution of the Stefan problem related to the soil thawing effect of composite piles, GMsFEM was used. According to results showing the accuracy of the method, the usage of four basis functions leads to the results with error in L2 norm lower than 1%. The errors in energy norms are also very respectable. Thus, usage of GMsFEM for such problems is recommended.
Results show that the process of total freezing of piles with steel fiber inclusion in winter runs faster. Analysis shows that the length between isoclines of phase change is about 47 cm for structured and 37 for random distribution of fibers. Freezing of pile with steel fibers comes on the 239th day, and a pile with basalt inclusion freezes on the 250th day in the structured case. For random distribution, similar days are 245th and 252nd. Therefore, there are 11 and 7 days' delay between freezing of piles with steel and fiber inclusions, depending on the distribution.

Conclusions
In this work, the effect of a composite pile on the temperature field of the surrounding frozen soil using GMsFEM is considered. Results of calculations using

Discussion
For the solution of the Stefan problem related to the soil thawing effect of composite piles, GMsFEM was used. According to results showing the accuracy of the method, the usage of four basis functions leads to the results with error in L 2 norm lower than 1%. The errors in energy norms are also very respectable. Thus, usage of GMsFEM for such problems is recommended.
Results show that the process of total freezing of piles with steel fiber inclusion in winter runs faster. Analysis shows that the length between isoclines of phase change is about 47 cm for structured and 37 for random distribution of fibers. Freezing of pile with steel fibers comes on the 239th day, and a pile with basalt inclusion freezes on the 250th day in the structured case. For random distribution, similar days are 245th and 252nd. Therefore, there are 11 and 7 days' delay between freezing of piles with steel and fiber inclusions, depending on the distribution.

Conclusions
In this work, the effect of a composite pile on the temperature field of the surrounding frozen soil using GMsFEM is considered. Results of calculations using GMsFEM show good accuracy. It can be said that when using materials with higher thermal conductivity, the thermal regime restores faster. In the future, we plan to consider the problems in the 3D statement and work with the piles with more complex compositions.