- Mining Institute of the Ural Branch of the Russian Academy of Sciences, Perm, Russia
This study examines unsteady heat transfer in the soil mass surrounding an inclined access mine working. Initially, the soils are considered frozen, but due to contact with warmer air entering the mine during summer, they undergo seasonal thawing. The primary challenge in the calculation of heat transfer is the loosening of soils near the upper boundary of the freezing front, caused by thaw consolidation—an effect particularly common in water-saturated, fine-grained soils. This phenomenon affects both the propagation rate of the pore moisture phase transition front and the stability of frozen soils in this zone. A bond-based peridynamic approach is employed to model heat transfer throughout the frozen and thawing soils surrounding the mine, using a case study of seasonal thawing near an access working in a mine located in the Republic of Yakutia. Numerical simulation was conducted using integral equations describing particle interactions within a finite influence region. The loosening zone is calculated using a simplified solution of the peridynamic motion equation. Model verification was performed to ensure the independence of the solution from numerical method parameters. The temporal evolution of thawing zone depth was analyzed for various thermal properties of the soil. For the first time, it has been demonstrated that thaw consolidation significantly affects the temperature field in soils surrounding a mine working.
1 Introduction
Mining in the Arctic regions presents several challenges, the most significant of which are permafrost and extreme climatic conditions (Solovyov, 2009). The continuous supply of air through access mine workings inevitably leads to the thawing of the surrounding frozen soil mass, significantly reducing its load-bearing capacity as temperatures rise above freezing point (Wu et al., 2006; Wen et al., 2011). Under certain conditions, the thawed zone can expand indefinitely, and once it reaches a critical size, it negatively affects the stability of mine workings. According to estimates (Dyadkin, 1968), when the thawing zone extends 3–4 m, the normal operation of horizontal and inclined mine workings becomes virtually impossible.
The temperature of the air supplied to mines varies daily and seasonally (Novopashin and Kurilko, 2010). Repeated freeze-thaw cycles of soils near ventilated mine workings significantly reduce their strength within the thawed zone. According to (Dyadkin, 1968), one of the key requirements for the thermal regime of mines in permafrost regions is preventing the progressive thawing of frozen soils. Additionally, minimizing the number of freeze-thaw cycles around mine workings is crucial (Novopashin and Kurilko, 2010).
The primary measures for maintaining the stability of workings in frozen soils include passive thermal insulation of their surfaces (Galkin et al., 1992) and the use of anchoring (Wu et al., 2006) or concrete lining (Dyadkin, 1968) to prevent the collapse of thawed soil.
Developing effective methods to ensure the stability of mine workings in permafrost requires predicting the temperature field of the soil mass and determining the thawing zones. Heat transfer in frozen soils with phase transitions is described by the Stefan problem (Solovyov, 2009; Lubini Tshumuka et al., 2022). Key aspects of modeling include phase transitions of moisture, heat exchange at the boundary with the working, nonlinear changes in the thermal properties of soils when crossing the zero isotherm, and moisture migration (Alzoubi et al., 2020). Due to the complexity of these processes, an analytical solution to the Stefan problem is not feasible, so heat transfer calculations in frozen soils are typically performed using numerical methods.
The most well-known numerical approaches for solving the heat transfer equation in frozen soils are:
1. Explicit tracking of the phase transition front (Carslaw and Jaeger, 1964; Meyer, 1978; Nikolaev et al., 2024).
2. The enthalpy method (Crowley, 1978; Levin et al., 2021).
3. The effective heat capacity method (Popov et al., 2018; Samarskii and Moiseenko, 1965; Bonacina et al., 1973).
4. The source term method (Voller and Swaminathan, 1991).
The main differences between these approaches lie in how the time derivative of soil internal energy is computed in the left-hand side of the heat transfer equation, while the term responsible for heat diffusion remains the same and is discretized using finite difference (Levin et al., 2021; Voller and Swaminathan, 1991), finite volume (Alzoubi et al., 2020; Hu and Wang, 2021), or finite element methods (Spiridonov et al., 2023; Zhelnin et al., 2022). All these approaches enable a quantitative description of heat transfer patterns in frozen soils but assume the continuity of the considered medium. However, in the case of thawing frozen soils around an inclined mine working, the situation differs. According to Dyadkin (1968), thawing is inevitably accompanied by soil consolidation (as some of the generated water drains into the working through cracks). Together with lining deformation under load and the sinking of its supports into the ground, this can lead to the formation of a void near the upper boundary of the thawing zone.
Even for small thawing zones, the consolidation and subsidence of thawed soil can lead to soil loosening near the upper boundary of the phase transition front. Accounting for this phenomenon within the framework of continuum mechanics is challenging: when fractures or voids form, the classical concept of spatial derivatives becomes inapplicable (Bobaru and Duangpanya, 2010).
The nonlocal heat transfer models can be used to describe such processes. One of the widely used formulations of continuum mechanics is peridynamics (Nikolaev et al., 2022). Unlike differential spatial operators, peridynamics employs integral equations that describe interactions between particles within a finite influence region (horizon).
In numerical heat transfer simulation using the peridynamic approach, the computational domain is divided into particles, each possessing a finite volume. Interactions occur not only with the nearest neighbors (as in the finite volume method, where heat exchange occurs through adjacent faces) but also with all particles within the interaction horizon, determined by the distance δ from the center of the considered particle. This distance typically exceeds the discretization step (the initial distance between particle centers) by several times (Ladányi and Gonda, 2021). Due to this nonlocal interaction mechanism, which is not tied to specific boundaries between particles, it is possible to accurately model the gradual reduction in heat transfer when material fractures develop.
Studies on heat transfer with phase transitions within the peridynamic framework have been presented in several previous works. There are two main formulations of the peridynamic model: bond-based and state-based (Javili et al., 2019; Anicode and Madenci, 2022). The bond-based approach is more commonly applied in heat transfer problems. The first formulation of an unsteady heat conduction problem using the peridynamic formalism was proposed in Bobaru and Duangpanya (2010). In Chen and Bobaru (2015), the choice of an integral kernel in the heat conduction equation was studied to ensure consistency with the classical heat conduction solution. The peridynamic approach has also been applied to the Stefan problem in studies (Nikolaev et al., 2022). A peridynamic model of the phase transition of water in saturated porous media has been developed, considering nonlocal heat transfer effects and mass redistribution during thawing and freezing. Additionally, the problem of heat transfer in nonlocal formulations has been explored for cases of rotational and spherical symmetry.
However, heat transfer processes in geomaterials with phase transitions of water, in the presence of voids and loosening zones, have not been previously studied within the peridynamic approach. This aspect is of primary importance in our study, as the thawing of frozen soils around an inclined mine working may lead to the formation of a void near the upper boundary of the thawing zone due to the consolidation of thawed soils and the deformation of the lining under load. This introduces additional thermal resistance and reduces heat flow into the frozen soils above the mine working. On the other hand, consolidated soil at the bottom of the thawing zone can facilitate more intensive heat transfer and more rapid thawing during the summer.
To identify heat transfer patterns that account for phase transitions of pore moisture and geomaterial discontinuities in the phase transition zone, this study formulates the heat transfer problem in a nonlocal framework using the bond-based peridynamic formulation. In addition, the peridynamic equation of motion for particles in a gravitational field is solved to determine the intensity of thawed soil compaction. The analysis is carried out using a case study on the design of an access mine working in Yakutia, in the permafrost zone. A numerical solution is obtained, the model is validated, and the results are compared with those from the classical thermodynamic formulation.
2 Methodology
A two-dimensional computational domain is considered in the vertical cross-section of the mine working. The soil mass is represented as a set of particles, each characterized by a position vector
Each particle has a finite volume, specific enthalpy H, and temperature T. A particle with position vector

Figure 1. Schematic representation of the computational domain and the neighborhood of particle x within its interaction horizon.
2.1 Heat transfer
The heat transfer dynamics is simulated using the peridynamic heat conduction equation, which accounts for nonlocal interactions between particles. For each particle with position
where
The change in the specific enthalpy of the soil was used to calculate its temperature change. To do this, the equation of state, which relates specific enthalpy to temperature, was applied:
where
Using condition Equation 4, we assume that the soil is initially fully saturated with ice, prior to the onset of thawing caused by heat exchange with the excavation wall. It is further assumed that the phase transition of moisture occurs over a relatively narrow temperature interval,
Dirichlet boundary conditions within the peridynamic approach are formulated as follows. Several additional layers of particles,
a fixed temperature

Figure 2. Setting Dirichlet boundary conditions at the outer boundary using ghost particles, and setting Robin boundary conditions on the mine working wall by explicitly defining the heat flux from the air to the boundary particles.
Otherwise, Robin boundary conditions are implemented on the wall of the mine working:
where
In this case, no ghost particle layers are added. Instead, a heat flux
In addition to the peridynamic heat conduction Equation 1, we will also use the classical heat conduction equation in the form (Semin et al., 2024; Levin et al., 2023):
with boundary conditions Equations 5, 6. This is necessary for the verification of the peridynamic model. Its numerical solution will be determined using the finite difference method on a regular grid with a fixed step size h. The absence of mesh refinement at the boundary with the mine working was motivated by the desire to maintain a situation similar to the peridynamic setting, where particles are typically of the same size. Furthermore, a coarse grid at the wall of the working is adequate for sufficiently accurate calculation of weak temperature gradients, with smooth changes in air temperature across seasons.
As an initial condition in both formulations, the assumption of a homogeneous temperature field will be used:
The air temperature
where
2.2 Settlement of thawed soils
The simulation of thawed soil settlement within the peridynamic approach is performed by solving the integral motion equation (Song et al., 2020; Van Le and Bobaru, 2018):
where
where
The stiffness coefficient
The settlement of thawed wet soil occurs due to the consolidation, which is characterized by the parameter
In practice, the maximum value of the parameter P is unlikely to be achieved, as this would require applying a large compressive force. Moreover, this would necessitate the complete drainage of all thawed moisture from the soil, which is also typically not the case, as a small amount of bound moisture usually remains in the soil. In the situation we are considering, Formula 14 provides an upper estimate for the possible range of values of the coefficient P.
In this study, the most pessimistic scenario is assumed, where all the pore space of the soil is initially filled with ice. The parameter P is introduced into the expression for
As a result, particles in the thawed soil zone shift under the influence of gravity (see Figure 3), and near the upper boundary of the thawed zone, a rarefaction zone forms, leading to an increase in the denominator under the integral in Equation 1. Thus, an additional thermal resistance arises in the vicinity of the freezing front.

Figure 3. Heat transfer implementation through the discontinuity of the soil mass at the upper boundary of the thawed zone.
In general, approach Equation 15 is rather non-standard, as within its framework, the consolidation deformation of a single soil particle does not depend on the forces acting on it. Most often, consolidation deformation is accounted for using the equation of water filtration in pores proposed by Biot (1973). However, given that the characteristic times for moisture migration are significantly shorter than the characteristic times for the movement of the phase transition front in the problem under consideration, we assume that all gravitational moisture flows out of the pore space of the soil instantaneously. For thawed soils we are considering, hydraulic diffusivity is an order of magnitude higher than thermal diffusivity, which supports this conclusion. However, this may not be the case for other types of soils.
This study investigates the dynamics of the thawed zone over several years, so it is important to understand how the soils will behave after a “thaw-freeze” cycle. Due to the outflow of pore water from the thawed zone in the first year, when this zone re-freezes, a nonzero additional term will remain in the right-hand side of Equation 16. Therefore, when considering the re-freezing of the soil, the new term in Equation 15 should not be excluded. However, in this case, the compaction coefficient may differ slightly from the value of
Another important aspect of the geomechanical problem is modeling the collapse of frozen soils above the upper boundary of the thawing front. According to Dyadkin (1968), the width of the thawed zone may exceed the span of the stable frozen soil exposure, leading to periodic collapses and abrupt increases in pressure on the lining of the mine working.
To account for this effect in the calculation of the force
The analysis of the strength of the mine working lining and the conditions under which critical stresses develop in the frozen soils above the thawed zone is beyond the scope of this study and will be the subject of future research. In certain situations, the nature of collapses will affect the additional thermal resistance at the upper boundary of the thawed zone. However, for simplicity, it can be assumed that the frozen soils maintain their structural integrity, and Equation 11 in this case will be more useful for calculating the consolidation and settlement of thawed soils.
Dirichlet boundary conditions at the outer boundary of the calculation domain
in the mechanical problem are formulated in the same way as in the thermodynamic problem—by assigning several layers of ghost particles and fixing their displacements. In cases where one of the displacements is not fixed (for example, displacement along the vertical axis on the side walls of the calculation domain), the ghost particles move along this axis together with the corresponding boundary particle.
Neumann boundary conditions at the wall boundary of the mine working are naturally set without adding ghost particle layers. Physically, this is equivalent to a condition of zero pressure on the mine working wall. Neumann boundary conditions at the upper boundary of the calculation domain are set by applying a specific force to each boundary particle.
3 Results and discussion
Calculations were carried out for the conditions of an access mine working passing through frozen soils from a quarry berm. The angle of the transport incline is approximately
At the entrance to the inclined mine working, the path was cut through aleurites and sandstones. According to geological surveys, the soils are unstable and prone to delamination. The thermophysical properties of the dry soil skeleton are presented in Table 1. The porosity ranged from
The thermophysical properties of the soil mass in its thawed and frozen states were recalculated using the following formulas (Semin et al., 2024):
where
It was assumed that the soil mass was fully ice-saturated at the initial time. The specific latent heat of ice melting was taken as 330 kJ/kg. The following numerical parameters were used when calculating the dependence of air temperature on time Equation 10:
The geometric dimensions of the arched mine working are shown in Figure 4. This figure also illustrates the discretization of the computational domain into particles. Near the external boundary, a different number of layers of ghost particles (from 1 to 3) was assigned; Figure 4, as an example, presents the case with two layers.
The discretized form of Equations 1, 11 is presented below:
where
Equation 19 is linked to Equations 20, 21 through the parameter
The time derivative is omitted on the left-hand side of Equation 20 because Equations 20, 21 are used to study the static equilibrium of the particle system. In this case, the iterative solution of the transient Equation 11 with the inertial term proves to be more computationally expensive than solving the system of linear Equation 21 for the unknown displacements.
For the numerical solution of Equation 19, an explicit integration scheme (first-order Euler method) was used, where the time step
where
The equations for particle displacement Equations 20, 21 were solved independently of the heat transfer Equation 19 every
The elastic modulus of frozen soils was assumed to be 50 GPa. Poisson’s ratio is set equal to 0.2. Based on these values, the peridynamic stiffness coefficient was calculated as
At the preliminary stage, verification of the heat transfer model was carried out. The verification process included checking the independence of the numerical solution from several parameters:
− The time step size,
− The size of the computational domain,
− The number of particles (spatial discretization step),
− The interaction horizon radius.
The quantitative criterion used to assess solution independence was the time dependence of the temperature of the soil mass at a distance of 2 m from the mine working wall (see Figure 4). Figures 5, 6 present the computed temperature-time dependencies at this point for various time step sizes, computational domain sizes, numbers of particles, and interaction horizon widths.

Figure 5. Temporal variations of soil temperature at a distance of 2 m from the mine working wall for different Courant numbers (a) and different computational domain sizes (b).

Figure 6. Temporal variations of soil temperature at a distance of 2 m from the mine working wall for different numbers of particles (a) and different horizon radii (b).
Based on a series of heat transfer simulations, the optimal numerical parameters were selected for the subsequent simulation of the thawed zone dynamics around the mine working. The computational domain size was set to 25 m, the Courant number
It is important to note that solving the heat transfer problem in the peridynamic formulation depends on the parameter
Despite the calibration of the model, achieving full agreement of the time dependencies for different
Overall, Figures 5, 6 show that the peridynamic approach provides additional smoothing of the time dependencies of temperature due to the “long-range” interactions between particles. In the classical enthalpy formulation of the Stefan problem Equation 8, the dependencies
In the case of
If the verification of the thermal model was conducted only for the initial calculation area with a mine working, then the verification of the mechanical model was performed separately by considering the deformation of a vertical soil column H = 10 m high under its own weight. In this case, the solution to the peridynamic problem can be compared with the exact analytical solution, which is obtained by integrating the 1D equilibrium equation. Both scenarios, with and without considering the consolidation of thawed soil, are examined. A comparative analysis of these solutions is presented in Figure 7.

Figure 7. Comparative analysis of numerical and exact analytical solutions for vertical soil column deformation under its own weight: (a) without consolidation of thawed soil and (b) with consolidation of thawed soil.
Figure 7a shows the situation without considering the consolidation of thawed soil. In this case, the vertical displacement
It was assumed that the
The peridynamic model accurately reproduces this power-law quadratic dependence for different parameters, ℎ and δ. However, the peridynamic displacement field in absolute values depends on the choice of δ, which is understandable and has been noted in the literature (Bobaru and Duangpanya, 2010). For further analysis, we selected the model with δ = 2 and calibrated the parameter c specifically to achieve the best fit between the model and the exact analytical solution for δ = 2:
The situation considering soil consolidation is shown in Figure 7b. As a result of heating the upper half of the column to a positive temperature, a fixed value of soil consolidation deformation P was applied. In this case, for the upper half of the column (
The figure shows that the peridynamic model accurately reproduces the exact analytical solution without any additional calibrations for different values of consolidation deformation. Thus, the calibrated and verified peridynamic model was subsequently used for further numerical simualtion of soil subsidence due to the thermal influence of air in a mine working.
Figure 8 presents the thawed zone distributions around the mine working for several time points and two different porosity values. For a porosity of 30%, two cases are shown: without considering the consolidation of thawed soils (Figure 8b) and with its consideration (Figure 8c). For a porosity of 3%, the thawed zone is visualized only for the case without consolidation (Figure 8a), as in this case, the factor has a negligible effect.

Figure 8. Propagation of the thawed zone at different soil porosities without considering settlement (a,b) and with considering settlement (c) of the thawed soils.
Calculations considering the consolidation of thawed soils were conducted with a reduced Courant number (CFL = 0.75) due to the local shortening of distances between particles caused by soil compaction. Figure 8c also visualizes the “voids” that appear near the upper boundary of freezing front. Overall, Figure 8 shows that soil mass porosity is a significant factor influencing the formation of the thawed soil zone. Additionally, the consolidation of thawed soils also plays a crucial role.
At the same time, the increase in thermal diffusivity within the thawed consolidated zone has a more pronounced effect on heat transfer than the increase in thermal resistance near the phase transition front, which, on the contrary, should slow down the heat exchange process.
Figure 9a presents temperature profiles along a vertical line extending from the axis of the excavation to the upper boundary of the computational domain. The curves correspond to a time of 200 days. The dashed curve represents the case with thawed soil consolidation, while the solid curve corresponds to the case without it. The figure shows that in the thawed soil zone, the dashed curve lies above the solid one, indicating a more intense heat transfer process when consolidation is considered. Furthermore, the temperature jump near the phase transition front is also greater for the dashed curve, suggesting the presence of thermal resistance in this region.

Figure 9. Temperature (a) and vertical displacement (b) profiles of the soil along the vertical path with (P = 0.0264) and without (P = 0) considering the settlement of thawed soils.
Figure 9b shows the vertical displacement profiles along the same line. In the case without considering the consolidation of thawed soils, the displacements do not exceed 5 mm in absolute value. When consolidation is taken into account, the maximum displacements are observed in the thawed zone, reaching 0.26 m (in absolute value) at the upper boundary of the thawed zone. Such large displacements at vertical coordinates y<6 m are primarily caused by the compaction of thawed soils.
The obtained maximum jump in vertical displacements at the transition from the thawed soil zone to the frozen soil zone is smaller than the horizon radius. As a result, a nonzero heat flux through the phase transition front above the thawing zone is maintained.
This raises the question of how physically accurate it is to preserve heat transfer through an air layer or another gas fraction when such wide separation (or void) zones form in the soil. It is likely that beyond a certain threshold
Overall, the behavior of thawing soil above the excavation roof represents just one aspect of a broader set of phenomena that warrant further investigation. A more detailed consideration of various coupled processes—such as support deformation, air migration through pores, and others—is required. These aspects, in turn, demand proper parameterization based on field data or scaled laboratory testing.
4 Conclusion
A peridynamic model of heat transfer in permafrost soil around a mine working, accounting for the consolidation of thawed soil, has been developed, numerically implemented, and verified. This model highlights the significant impact of thawed soil deformations on seasonal temperature variations around the access mine working, driven by changes in the temperature of the atmospheric air supplied through it.
The numerical solution is highly sensitive to parameters such as the initial soil porosity and the degree of consolidation, emphasizing the need for their accurate determination during preliminary engineering surveys. The peridynamic approach has demonstrated its effectiveness in solving the heat transfer equation with moisture phase transitions, even under strong spatiotemporal variability in soil thermophysical properties. By incorporating long-range interactions, the peridynamic approach produces smoother temperature fields, particularly in high-moisture soils where enthalpy-temperature relationships exhibit sharp discontinuities. At the same time, it yields results that are quantitatively comparable to the classical Stefan problem solved with the Laplace diffusion operator.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
MS: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review and editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. The study was carried out with the financial support of the Ministry of Science and Higher Education of the Russian Federation within the state assignment (project no. 122030100425–6).
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Alzoubi, M. A., Xu, M., Hassani, F. P., Poncet, S., and Sasmito, A. P. (2020). Artificial ground freezing: a review of thermal and hydraulic aspects. Tunn. Undergr. Space Technol. 104, 103534. doi:10.1016/j.tust.2020.103534
Anicode, S. V. K., and Madenci, E. (2022). Bond-and state-based peridynamic analysis in a commercial finite element framework with native elements. Comput. Methods Appl. Mech. Eng. 398, 115208. doi:10.1016/j.cma.2022.115208
Biot, M. A. (1973). Nonlinear and semilinear rheology of porous solids. J. Geophys. Res. 78 (23), 4924–4937. doi:10.1029/jb078i023p04924
Bobaru, F., and Duangpanya, M. (2010). The peridynamic formulation for transient heat conduction. Int. J. Heat. Mass Transf. 53 (19–20), 4047–4059. doi:10.1016/j.ijheatmasstransfer.2010.05.024
Bonacina, C., Comini, G., Fasano, A., and Primicerio, M. (1973). Numerical solution of phase-change problems. Int. J. Heat. Mass Transf. 16 (10), 1825–1832. doi:10.1016/0017-9310(73)90202-0
Chen, Z., and Bobaru, F. (2015). Selecting the kernel in a peridynamic formulation: a study for transient heat diffusion. Comput. Phys. Commun. 197, 51–60. doi:10.1016/j.cpc.2015.08.006
Crowley, A. B. (1978). Numerical solution of Stefan problems. Int. J. Heat. Mass Transf. 21 (2), 215–219. doi:10.1016/0017-9310(78)90225-9
Dyadkin, Yu. D. (1968). Fundamentals of mining thermophysics for mines of the North. Moscow: Nedra, 255.
Galkin, A. F., Kiselev, V. V., and Kurilko, A. S. (1992). Sprayed concrete thermal protection support. Yakutsk: YANC SO RAN, 164.
Hu, T., and Wang, T. (2021). A finite volume-based model for the hydrothermal behavior of soil under freeze–thaw cycles. PLoS One 16 (6), e0252680. doi:10.1371/journal.pone.0252680
Isiet, M., Mišković, I., and Mišković, S. (2021). Review of peridynamic modelling of material failure and damage due to impact. Int. J. Impact Eng. 147, 103740. doi:10.1016/j.ijimpeng.2020.103740
Javili, A., Morasata, R., Oterkus, E., and Oterkus, S. (2019). Peridynamics review. Math. Mech. Solids 24 (11), 3714–3739. doi:10.1177/1081286518803411
Ladányi, G., and Gonda, V. (2021). Review of peridynamics: theory, applications, and future perspectives. J. Mech. Eng. 67 (12), 666–681. doi:10.5545/sv-jme.2021.7289
Levin, L., Golovatyi, I., Zaitsev, A., Pugin, A., and Semin, M. (2021). Thermal monitoring of frozen wall thawing after artificial ground freezing: case study of Petrikov Potash Mine. Tunn. Undergr. Space Technol. 107, 103685. doi:10.1016/j.tust.2020.103685
Levin, L., Semin, M., and Golovatyi, I. (2023). Analysis of the structural integrity of a frozen wall during a mine shaft excavation using temperature monitoring data. Fract. Struct. Integr. 17 (63), 1–12. doi:10.3221/igf-esis.63.01
Liu, W., and Hong, J. W. (2012). Discretized peridynamics for linear elastic solids. Comput. Mech. 50, 579–590. doi:10.1007/s00466-012-0690-1
Lubini Tshumuka, A., Krimi, A., and Fuamba, M. (2022). Modeling heat transfer through permafrost soil subjected to seasonal freeze-thaw. Land 11 (10), 1770. doi:10.3390/land11101770
Meyer, G. H. (1978). The numerical solution of Stefan problems with front-tracking and smoothing methods. Appl. Math. Comput. 4 (4), 283–306. doi:10.1016/0096-3003(78)90001-2
Nikolaev, P., Jivkov, A., Rajabi, H., Yan, H., Zhu, X., and Sedighi, M. (2024). Semi-analytical predictive model for natural and artificial thawing of circular ground-ice walls. Comput. Geotech. 171, 106394. doi:10.1016/j.compgeo.2024.106394
Nikolaev, P., Sedighi, M., Jivkov, A. P., and Margetts, L. (2022). Analysis of heat transfer and water flow with phase change in saturated porous media by bond-based peridynamics. Int. J. Heat. Mass Transf. 185, 122327. doi:10.1016/j.ijheatmasstransfer.2021.122327
Novopashin, M. D., and Kurilko, A. S. (2010). Properties of geomaterials in permafrost conditions and their impact on mining technology and stability of workings. Min. Inf. Anal. Bull. 4 (12), 21–33.
Popov, V. I., and Kurilko, A. S. (2018). Approximate method to solve problems of heat and mass transfer under water freezing in permafrost rocks. Min. Inf. Anal. Bull. 12, 57–64. doi:10.25018/0236-1493-2018-12-0-57-64
Samarskii, A. A., and Moiseenko, B. D. (1965). Economical through-count scheme for the multidimensional Stefan problem. J. Comput. Math. Math. Phys. 5 (5), 816–827.
Semin, M., Golovatyi, I., Levin, L., and Pugin, A. (2024). Enhancing efficiency in the control of artificial ground freezing for shaft construction: a case study of the Darasinsky potash mine. Eng. Technol. 18, 100710. doi:10.1016/j.clet.2023.100710
Solovyov, D. E. (2009). “Methods for calculating temperature and ventilation regimes of an unsteady mine network in permafrost zones,”. Dissertion (Yakutsk: Cand. Sci. Tech.), 161.
Song, Y., Liu, R., Li, S., Kang, Z., and Zhang, F. (2020). Peridynamic modeling and simulation of coupled thermomechanical removal of ice from frozen structures. Meccanica 55, 961–976. doi:10.1007/s11012-019-01106-z
Spiridonov, D., Stepanov, S., and Vasiliy, V. (2023). An online generalized multiscale finite element method for heat and mass transfer problem with artificial ground freezing. J. Comput. Appl. Math. 417, 114561. doi:10.1016/j.cam.2022.114561
Sun, W. K., Zhang, L. W., and Liew, K. M. (2020). A smoothed particle hydrodynamics–peridynamics coupling strategy for modeling fluid–structure interaction problems. Comput. Methods Appl. Mech. Eng. 371, 113298. doi:10.1016/j.cma.2020.113298
Van Le, Q., and Bobaru, F. (2018). Objectivity of state-based peridynamic models for elasticity. J. Elast. 131 (1), 1–17. doi:10.1007/s10659-017-9641-6
Voller, V. R., and Swaminathan, C. R. (1991). ERAL source-based method for solidification phase change, Numer. Heat Transf. Part B Fundam. 19 (2), 175–189. doi:10.1080/10407799108944962
Weißenfels, C. (2022). Peridynamics, simulation of additive manufacturing using meshfree methods: with focus on requirements for an accurate solution, 125–138.
Wen, Z., Zhang, T., Sheng, Y., Ma, W., Wu, Q., Feng, W., et al. (2011). Managing ice-rich permafrost exposed during cutting excavation along Qinghai–Tibetan railway: experiences and implementation. Eng. Geol. 122 (3–4), 316–327. doi:10.1016/j.enggeo.2011.07.012
Wu, H., Bandopadhyay, S., and Izaxon, V. U. (2006). “Optimum insulation for engineering control of mine thermal regime around a mine airway in permafrost,” in 11th US/North American Mine Ventilation Symposium 2006: Proceedings, Pennsylvania, USA, 5–7 June 2006 (Boca Raton, FL: CRC Press), 271.
Keywords: numerical simulation, peridynamics, permafrost, mine working, thermomechanical problem
Citation: Semin M (2025) Peridynamic simulation of permafrost thawing around an inclined mine working. Front. Mech. Eng. 11:1590285. doi: 10.3389/fmech.2025.1590285
Received: 09 March 2025; Accepted: 03 July 2025;
Published: 16 July 2025.
Edited by:
Eswaramoorthy Muthusamy, Shri Mata Vaishno Devi University, IndiaReviewed by:
Ravi Goyal, Uttaranchal University, IndiaAmit Kumar, Shri Mata Vaishno Devi University, India
Copyright © 2025 Semin. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mikhail Semin, c2VtaW5tYUBpbmJveC5ydQ==