Numerical Simulation of Thawing Process in Frozen Soil

Yan, Zhaoyu; Pan, Wei; Fang, Junjie; Liu, Zihui

doi:https://doi.org/10.1155/2020/8822320

Geofluids

On this page

Abstract Introduction Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Thermal-Hydraulic-Mechanical (THM) Coupling Behaviour of Fractured Rock Masses

View this Special Issue

Research Article | Open Access

Volume 2020 | Article ID 8822320 | https://doi.org/10.1155/2020/8822320

Numerical Simulation of Thawing Process in Frozen Soil

Zhaoyu Yan,¹Wei Pan,¹Junjie Fang,¹and Zihui Liu¹

Academic Editor: Yanlin Zhao

Received08 Aug 2020

Revised11 Sept 2020

Accepted28 Sept 2020

Published23 Oct 2020

Abstract

Permafrost has been thawing faster due to climate change which would release greenhouse gases, change the hydrological regimes, affect buildings above, and so on. It is necessary to study the thawing process of frozen soil. A water-heat coupling model for frozen soil thawing is established on Darcy’s law and Heat Transfer in Porous Media interfaces in Comsol Multiphysics 5.5. Three curves of total liquid water volume, minimum temperature, and total heat flux in the thawing process are obtained from a numerical simulation. The distributions of liquid water, temperature, and pressure based on time are simulated too. The liquid water distribution is consistent with the total liquid water volume curve. The temperature distribution is confirmed by the minimum temperature and total heat flux curve. The pressure distribution represents ice in the frozen soil that generates negative pressure during the melting process. The numerical simulation research in this article deepens the understanding of the internal evolution in the process of frozen soil thawing and has a certain reference value for subsequent experimental research and related applications.

1. Introduction

There are large areas of permafrost in boreal and high-elevation regions such as Qinghai-Tibet Plateau [1]. However, hydrological and ecological changes are exacerbating because of climate change [2–7]. Permafrost has been thawing faster affected by the global temperature rise. At the same time, greenhouse gases such as carbon dioxide and methane released by the thawing of permafrost will in turn exacerbate global warming [8–10]. This trend will continue in the foreseeable decades [11]. The distributions of surface water and subface water pathways can be changed by permafrost thaw or changing patterns of seasonal subsurface ice [12–14], so that permafrost thaw makes a big difference to hydrological regimes. In addition, the freeze-thaw cycles of soil can also affect the buildings and facilities above [15, 16]. Thus, soil under freezing and thawing conditions is the key point for studying the impact of climate change on boreal and high-elevation regions [17].

It is a mainstream research direction to carry out experimental research for variations of soil under freezing and thawing conditions. Wang et al. [18] studied physicomechanical property changes of Qinghai-Tibet clay due to cyclic freezing and thawing. Zhang et al. [19] studied the variations of the temperatures and volumetric unfrozen water contents during a freezing-thawing process. Wang et al. [20–22] studied fracturing characteristics and mechanical behavior of rock. Darrow et al. [23] measured the zeta potential of cation-treated soils to find its implication on unfrozen water mobility. Cao et al. [24] studied the effects of cyclic freeze-thaw treatments on the fracture characteristics of sandstone. Zhou et al. [25] studied the mechanical behavior of frozen loess under freeze-thaw cycles. Han et al. [26] studied the effect of freeze-thaw cycles on the shear strength of saline soil. Zhao et al. [27–30] studied the dynamic strength of frozen sulfate saline silty clay under cyclic loading. De Guzman et al. [31] discovered that the shear strength of frozen soil samples was reduced by as much as 50% upon freezing and thawing conditions. In addition, some scholars studied the shear behavior of frozen soil-concrete interface under freezing and thawing conditions [32, 33]. Besides the experimental study, numerical simulation is also a feasible approach. Lin et al. [34–36] established the stress field of an open flaw tip under uniaxial compression and shear models. Li et al. [37] developed a THM coupling model to research freeze-thaw failure mechanisms and engineer safeguarding measures. Zhao et al. [38–41] established models about rock behavior which contain the THM model. Grenier et al. [42] studied groundwater flow and heat transport for systems undergoing freeze-thaw. He et al. [43] proposed a coupled heat-moisture-deformation model for saturated frozen soils. Considering the effects of water flow in porous resulted by Soret effect and segregation potential on seepage velocity and water pressure distribution, Tan et al. [44] established a TH coupling model for the construction of Galonglatunnel. Amiri et al. [45] extended the theory of soil freezing curve paradigm. Kurylyk and Watanabe [46] summarized and examined recent models simulating coupled heat and water transport in cold regions.

Soil under freeze-thaw conditions has been studied from a variety of perspectives. However, the variation in the frozen soil cannot be displayed intuitively. Because it is difficult to observe the internal evolution of frozen soil during the thawing process, numerical simulation is the approach which is chosen to study the variations of water, temperature, heat, and pressure in the process. The research conclusion has a certain reference value for the subsequent experimental research and related applications.

2. Model Building

2.1. Model Definition

The numerical simulation of the frozen soil thawing process is carried out by the Comsol Multiphysics 5.5 software. The interfaces selected to simulate the thawing process of frozen soil are Darcy’s law and Heat Transfer in Porous Media.

The numerical model is shown in Figure 1. A flow channel 3 meters in length and 1 meter in width which is from left to right contains a squared ice inclusion of 0.333 meter in length. The initial temperature of the ice inclusion is . The temperature in the porous channel is . The fluid flow is driven by a hydraulic head gradient of over the length of the channel.

In order to facilitate research work, the model is built on these assumptions: no thermal dispersion is included in the heat transfer equation, and water density and dynamic viscosity are considered constant with respect to the temperature (see Table 1).

2.2. Model Equations

Using the storage model node, Darcy’s law interface contains an implementation of Darcy’s law which explicitly includes an option to define the linearized storage (SI unit: 1/Pa) using the compressibility of the fluid and the porous matrix:

In this case, the storage is defined as where is the water saturation, is the porosity, and is the effective compressibility, which is a combined value of water, ice, and solid matrix compressibility. For simplicity, the gravity term in Equation (1) is neglected. The variable is a mass source which represents the additional liquid water due to the melting of the ice inclusion:

Here, and are the ice and liquid water densities. The liquid water saturation depends on the phase change as

Here, is the residual liquid water saturation, is a smooth step function defined in the phase change material node. The step function is zero for temperatures below the melting temperature , and it equals to 1 for temperatures above .

The heat transfer in the porous medium is described by the following equation:

Here, is the equivalent value of density (kg/m³) and heat capacity at constant pressure (J/kg·K), is the effective thermal conductivity (W/m·K), is the temperature (K), and is a heat source (W/m³). The variable is the effective fluid’s heat capacity at constant pressure.

2.3. Model Data

The benchmark example describes the ice-to-water phase change of an ice inclusion within a porous channel. The values for thermal properties of water, ice, and solid matrix are presented in Table 1.

A latent heat of melting, kJ/kg is used. Additional physical parameters used in the example are summarized in Table 2.

3. Results and Discussion

The simulation runs for 56 hours. After 19.2 hours, the ice completely melted, but the minimum temperature in the channel did not return to the initial 5°C until the 56th hour. With the melting of ice, the changes of various physical indicators in the channel can reflect the internal evolution of the thawing process in frozen soil, which is worthy of further discussion.

3.1. Numerical Simulation Curves

With the thawing of frozen soil, the total liquid water volume changed as shown in Figure 2. After the simulation started running, the ice in the soil gradually melted, and the total liquid water volume in the channel began to increase. At about 19.2 h, the total liquid water volume in the channel rose to 1.11 and remained constant which demonstrates that the ice has completely melted.

At the beginning of the simulation, the temperature of the soil in the channel is +5°C, and the temperature of the ice is -5°C; at this time, the minimum temperature in the channel is the temperature of the ice, -5°C. After 56 hours, the minimum temperature in the channel returned to the initial 5°C, and the change curve of the minimum temperature in the channel is shown in Figure 3. As the figure shows, the temperature rose rapidly to -1°C in the first two hours after the ice melt. From then on until about 19.2 h, the minimum temperature rose at a very slow rate. In the narrow interval from 19.2 h to 19.9 h, the minimum temperature increased rapidly and exceeded the freezing point. Then, the rise of minimum temperature has experienced three stages: slow, fast, and slow. The minimum temperature in the channel returned to the initial value at 56 h.

In the model, fluid flows from left to right in the channel, so the right side of the rectangle is the outflow boundary of the system. Total heat flux via the outflow boundary is shown as Figure 4. Because the ice is in the left half of the channel, there is a certain distance from the right outflow boundary; the heat flux of the outflow boundary on the right side of the channel is not consistent with the heat change caused by ice melting; therefore, there is no obvious correspondence between Figures 3 and 4. The corresponding relationship between heat flux and temperature change will be described in the following content. The change of heat flux at the outflow boundary is driven by ice melting, so the heat flux is negative. The absolute value of total heat flux at the outflow boundary first increased and then decreased.

3.2. Liquid Water Distribution

The variation of liquid water distribution in the channel is shown in Figure 5. In Figure 5, there is the liquid water distribution in numerical simulation for 0 h, 5 h, 10 h, 15 h, 20 h, and 25 h, respectively. It can be seen that the ice begins to melt from the four corners of the square. The left side of the square ice is facing the liquid flow direction, and the melting speed of the two adjacent corners is higher than that of the right side. In the process of melting, the shape of ice gradually changed from square to round, and the scope of ice gradually reduced until it completely melted. It can be seen from the figure at 20 h and 25 h that the ice in the channel has completely melted after 20 hours, the liquid water was full of the channel, and the total liquid water volume is constant with the curve in Figure 2.

3.3. Temperature Distribution

The temperature distribution in the channel develops with the melting of ice which is shown in Figure 6. After the ice melted and heated up, the low temperature diffused around and then moved to the right as the liquid flowed. Finally, the low temperature flowed out of the channel and the initial 5°C was restored in the channel. Obviously, the part below zero is ice; from the figures in the time range of 0 h to 20 h in Figure 6, we can find that the distribution of ice in the melting process inferred from Figure 6 is consistent with Figure 5. Compared with Figure 3, temperature figures at 0 h, 0.083 h, 0.25 h, and 1 h verify the rapid rise of the minimum temperature in the first hour of the simulation. The melting of ice caused the surrounding water temperature to cool down; low-temperature water flowed to the right with the hydraulic head which would change the heat flux of the outflow boundary. The flow of low-temperature water in Figure 6 also explains the variation curve of total heat flux in Figure 4. In Figure 6, there are blue low-temperature water passing through the outflow boundary between 20 h and 30 h, which is in line with the larger absolute value of heat flux in Figure 4.

3.4. Pressure Distribution

The pressure distribution in the channel develops with the melting of ice which is shown in Figure 7. The streamlines in the figure are Darcy’s velocity field, which represents the flow path of fluid in porous media. It can be seen that the pressure changes in the channel are mainly concentrated on ice. With the beginning of ice melting, the pressure in the square ice area began to decrease. First of all, it changed from positive pressure to negative pressure from four corners and then completely changed to negative pressure in the ice area. The ice melted from a square to a circle, and the negative pressure area also changed to a circle. Then, the area of ice became smaller and smaller, and the absolute value of negative pressure gradually increased until the ice melted completely, and the pressure in the channel remained constant. In addition, the change of streamlines indicates that the water flow path gradually occupied the ice melting area.

4. Conclusions

A 2D water-heat coupling model is established to analyse the thawing process of frozen soil. The numerical simulation in this paper reveals the internal evolution of frozen soil in the thawing process from four aspects of water, temperature, heat, and pressure. The main conclusions are as follows: (1)The variation of liquid water distribution in the thawing process of frozen soil is simulated, which is in full compliance with the total liquid water volume curve, and can reflect the area of unmelted ice(2)The variation of temperature distribution in the thawing process of frozen soil is simulated. On the one hand, the variation of temperature distribution explains the minimum temperature curve; on the other hand, the minimum temperature water passing through the outflow boundary is consistent with total heat flux(3)The change of pressure and the flow path of water in frozen soil are simulated. It can be seen that the ice in the frozen soil generates negative pressure during the melting process(4)Though this numerical simulation has a reference value, it is necessary to carry out more field observation and research on permafrost thawing so that the actual melting mechanism of permafrost can be clearly understood. One potential approach is deploying sensors on a large scale in permafrost to monitor the variation of permafrost thawing

Data Availability

All data used to support the study is included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This work was supported by the National Key R&D Program of China (Grant No. 2018YFC0808404).

References

K. Xia, Y. Luo, and W. P. Li, “Simulation of freezing and melting of soil on the northeast Tibetan Plateau,” Chinese Science Bulletin, vol. 56, no. 20, pp. 2145–2155, 2011.
View at: Publisher Site | Google Scholar
W. R. Rouse, M. S. V. Douglas, R. E. Hecky et al., “Effects of climate change on the Freshwaters of arctic and subarctic north America,” Hydrological Processes, vol. 11, no. 8, pp. 873–902, 1997.
View at: Publisher Site | Google Scholar
M. C. Serreze, J. E. Walsh, F. S. Chapin et al., “Observational evidence of recent change in the northern high-latitude environment,” Climatic Change, vol. 46, no. 1/2, pp. 159–207, 2000.
View at: Publisher Site | Google Scholar
M. T. Jorgenson, C. H. Racine, J. C. Walters, and T. E. Osterkamp, “Permafrost degradation and ecological changes associated with a warmingclimate in central Alaska[J],” Climatic Change, vol. 48, no. 4, pp. 551–579, 2001.
View at: Publisher Site | Google Scholar
L. D. Hinzman, N. D. Bettez, W. R. Bolton et al., “Evidence and implications of recent climate change in northern Alaska and other Arctic regions,” Climatic Change, vol. 72, no. 3, pp. 251–298, 2005.
View at: Publisher Site | Google Scholar
D. W. Schindler and J. P. Smol, “Cumulative effects of climate warming and other human activities on freshwaters of arctic and subarctic north America,” Ambio, vol. 35, no. 4, pp. 160–168, 2006.
View at: Publisher Site | Google Scholar
Z. Zhongqiong, W. Qingbai, J. Guanli, G. Siru, C. Ji, and L. Yongzhi, “Changes in the permafrost temperatures from 2003 to 2015 in the Qinghai-Tibet Plateau,” Cold Regions Science and Technology, vol. 169, article 102904, 2020.
View at: Google Scholar
A. O. Anisimov, “Potential feedback of thawing permafrost to the global climate system through methane emission,” Environmental Research Letters, vol. 2, no. 4, article 045016, 2017.
View at: Google Scholar
O. V. Masyagina and O. V. Menyailo, “The impact of permafrost on carbon dioxide and methane fluxes in Siberia: a meta-analysis,” Environmental Research, vol. 182, article 109096, 2020.
View at: Google Scholar
IPCC, “Climate change 2001– the scientific basis,” Ksce Journal of Civil Engineering, vol. 19, no. 2, pp. 359–365, 2001.
View at: Google Scholar
Ipcc, T. Stocker, D. Qin et al., “The physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change,” Climate Change, vol. 2013, 2013.
View at: Google Scholar
R. F. Connon, W. L. Quinton, J. R. Craig, and M. Hayashi, “Changing hydrologic connectivity due to permafrost thaw in the lower Liard River valley, NWT, Canada,” Hydrological Prochydrological Processesesses, vol. 28, no. 14, pp. 4163–4178, 2014.
View at: Publisher Site | Google Scholar
B. L. Kurylyk, K. T. B. Macquarrie, and J. M. Mckenzie, “Climate change impacts on groundwater and soil temperatures in cold and temperate regions: implications, mathematical theory, and emerging simulation tools,” Earth-Science Reviews, vol. 138, pp. 313–334, 2014.
View at: Publisher Site | Google Scholar
A. Frampton and G. Destouni, “Impact of degrading permafrost on subsurface solute transport pathways and travel times,” Water Resources Research, vol. 51, no. 9, pp. 7680–7701, 2015.
View at: Publisher Site | Google Scholar
B. Yang, Z. Qin, Q. Zhou, H. Li, L. Li, and X. Yang, “Pavement damage behaviour of urban roads in seasonally frozen saline ground regions,” Cold Regions Science and Technology, vol. 174, p. 103035, 2020.
View at: Publisher Site | Google Scholar
Z. Yang, U. Dutta, F. Xiong, N. Biswas, and H. Benz, “Seasonal frost effects on the dynamic behavior of a twenty-story office building,” Cold Regions Science and Technology, vol. 51, no. 1, pp. 76–84, 2008.
View at: Publisher Site | Google Scholar
D. M. Lawrence, A. G. Slater, and S. C. Swenson, “Simulation of present-day and future permafrost and seasonally frozen ground conditions in CCSM4,” Journal of Climate, vol. 25, no. 7, pp. 2207–2225, 2012.
View at: Publisher Site | Google Scholar
D.-Y. Wang, W. Ma, X.-X. Chang, Z.-Z. Sun, W.-J. Feng, and J.-W. Zhang, “Physico-mechanical properties changes of Qinghai–Tibet clay due to cyclic freezing and thawing,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, pp. 4313–4319, 2005.
View at: Google Scholar
M. Zhang, X. Zhang, Y. Lai, J. Lu, and C. Wang, “Variations of the temperatures and volumetric unfrozen water contents of fine-grained soils during a freezing-thawing process,” Acta Geotechnica, vol. 15, no. 3, pp. 595–601, 2020.
View at: Publisher Site | Google Scholar
Y. Wang, H. Zhang, H. Lin, Y. Zhao, X. Li, and Y. Liu, “Mechanical behavior and failure analysis of fracture-filled gneissic granite,” Theoretical and Applied Fracture Mechanics, vol. 108, p. 102674, 2020.
View at: Publisher Site | Google Scholar
Y. Wang, H. Lin, Y. Zhao, X. Li, P. Guo, and Y. Liu, “Analysis of fracturing characteristics of unconfined rock plate under edge-on impact loading,” European Journal of Environmental & Civil Engineering, pp. 1–16, 2019.
View at: Publisher Site | Google Scholar
Y. Wang, H. Zhang, H. Lin, Y. Zhao, and Y. Liu, “Fracture behaviour of central-flawed rock plate under uniaxial compression,” Theoretical and Applied Fracture Mechanics, vol. 106, article 102503, 2020.
View at: Publisher Site | Google Scholar
M. M. Darrow, R. Guo, and T. P. Trainor, “Zeta potential of cation-treated soils and its implication on unfrozen water mobility,” Cold Regions Science and Technology, vol. 173, p. 103029, 2020.
View at: Publisher Site | Google Scholar
R.-H. Cao, C. Wang, R. Yao et al., “Effects of cyclic freeze-thaw treatments on the fracture characteristics of sandstone under different fracture modes: laboratory testing,” Theoretical and Applied Fracture Mechanics, vol. 109, p. 102738, 2020.
View at: Publisher Site | Google Scholar
Z. Zhou, W. Ma, S. Zhang, Y. Mu, and G. Li, “Effect of freeze-thaw cycles in mechanical behaviors of frozen loess,” Cold Regions Science and Technology, vol. 146, pp. 9–18, 2018.
View at: Publisher Site | Google Scholar
Y. Han, Q. Wang, N. Wang et al., “Effect of freeze-thaw cycles on shear strength of saline soil,” Cold Regions Science and Technology, vol. 154, pp. 42–53, 2018.
View at: Publisher Site | Google Scholar
Y. Zhao, Y. Lai, J. Zhang, and M. Liao, “A dynamic strength criterion for frozen sulfate saline silty clay under cyclic loading,” Cold Regions Science and Technology, vol. 173, p. 103026, 2020.
View at: Publisher Site | Google Scholar
Y. Zhao, L. Zhang, W. Wang, J. Tang, H. Lin, and W. Wan, “Transient pulse test and morphological analysis of single rock fractures,” International Journal of Rock Mechanics and Mining Sciences, vol. 91, pp. 139–154, 2017.
View at: Publisher Site | Google Scholar
Y. Zhao, L. Zhang, W. Wang, Q. Liu, L. Tang, and G. Cheng, “Experimental study on shear behavior and a revised shear strength model for infilled rock joints,” International Journal of Geomechanics, vol. 20, no. 9, p. 17, 2020.
View at: Google Scholar
C. Wang, Y. Zhao, Y. Zhao, and W. Wan, “Study on the interaction of collinear cracks and wing cracks and cracking behavior of rock under uniaxial compression,” Advances in Civil Engineering, vol. 2018, Article ID 5459307, 10 pages, 2018.
View at: Publisher Site | Google Scholar
E. M. B. De Guzman, D. Stafford, M. C. Alfaro, G. Doré, and L. U. Arenson, “Large-scale direct shear testing of compacted frozen soil under freezing and thawing conditions,” Cold Regions Science and Technology, vol. 151, pp. 138–147, 2018.
View at: Publisher Site | Google Scholar
P. He, Y. Mu, Z. Yang, W. Ma, J. Dong, and Y. Huang, “Freeze-thaw cycling impact on the shear behavior of frozen soil-concrete interface,” Cold Regions Science and Technology, vol. 173, p. 103024, 2020.
View at: Publisher Site | Google Scholar
J. Liu, P. Lv, Y. Cui, and J. Liu, “Experimental study on direct shear behavior of frozen soil–concrete interface,” Cold Regions Science and Technology, vol. 104-105, pp. 1–6, 2014.
View at: Publisher Site | Google Scholar
H. Lin, H. Yang, Y. Wang, Y. Zhao, and R. Cao, “Determination of the stress field and crack initiation angle of an open flaw tip under uniaxial compression,” Theoretical and Applied Fracture Mechanics, vol. 104, p. 102358, 2019.
View at: Publisher Site | Google Scholar
S. Xie, H. Lin, Y. Wang et al., “A statistical damage constitutive model considering whole joint shear deformation,” International Journal of Damage Mechanics, vol. 29, no. 6, pp. 988–1008, 2020.
View at: Publisher Site | Google Scholar
H. Lin, X. Zhang, Y. Wang et al., “Improved nonlinear Nishihara shear creep model with variable parameters for rock-like Materials,” Advances in Civil Engineering, vol. 2020, 15 pages, 2020.
View at: Publisher Site | Google Scholar
G. Li, N. Li, Y. Bai, N. Liu, M. He, and M. Yang, “A novel simple practical thermal-hydraulic-mechanical (THM) coupling model with water-ice phase change,” Computers and Geotechnics, vol. 118, p. 103357, 2020.
View at: Publisher Site | Google Scholar
Y. L. Zhao, W. J. Wang, Y. S. Zhao, and W.-h. Gao, “3D dual medium model of thermal-hydro-mechanical coupling and its application,” Journal of China University of Mining & Technology, vol. 39, no. 5, pp. 709–715, 2010.
View at: Google Scholar
Y. Zhao, Y. Wang, W. Wang, L. Tang, Q. Liu, and G. Cheng, “Modeling of rheological fracture behavior of rock cracks subjected to hydraulic pressure and far field stresses,” Theoretical and Applied Fracture Mechanics, vol. 101, pp. 59–66, 2019.
View at: Publisher Site | Google Scholar
Y. Zhao, Y. Wang, W. Wang, W. Wan, and J. Tang, “Modeling of non-linear rheological behavior of hard rock using triaxial rheological experiment,” International Journal of Rock Mechanics and Mining Sciences, vol. 93, pp. 66–75, 2017.
View at: Publisher Site | Google Scholar
Y. Zhao, L. Zhang, W. Wang, W. Wan, and W. Ma, “Separation of elastoviscoplastic strains of rock and a nonlinear creep model,” International Journal of Geomechanics, vol. 18, no. 1, p. 18, 2018.
View at: Google Scholar
C. Grenier, H. Anbergen, V. Bense et al., “Groundwater flow and heat transport for systems undergoing freeze-thaw: intercomparison of numerical simulators for 2D test cases,” Advances in Water Resources, vol. 114, pp. 196–218, 2018.
View at: Publisher Site | Google Scholar
M. He, N. Li, and N. F. Liu, “Analysis and validation of coupled heat-moisture-deformation model for saturated frozen soils,” Chinese Journal of Geotechnical Engineering, vol. 34, no. 10, pp. 1858–1865, 2012.
View at: Google Scholar
X. Tan, W. Chen, H. Tian, and J. Cao, “Water flow and heat transport including ice/water phase change in porous media: numerical simulation and application,” Cold Regions Science and Technology, vol. 68, no. 1-2, pp. 74–84, 2011.
View at: Publisher Site | Google Scholar
E. A. Amiri, J. R. Craig, and B. L. Kurylyk, “A theoretical extension of the soil freezing curve paradigm,” Advances in Water Resources, vol. 111, pp. 319–328, 2018.
View at: Publisher Site | Google Scholar
B. L. Kurylyk and K. Watanabe, “The mathematical representation of freezing and thawing processes in variably-saturated, non-deformable soils,” Advances in Water Resources, vol. 60, pp. 160–177, 2013.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Zhaoyu Yan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

871

Downloads

1108

Citations