DOI: 10.17516/1997-1397-2021-14-1-47-56 УДК 517.95
Mathematical Model of Fluids Motion in Poroelastic Snow-ice Cover
Margarita A. Tokareva* Alexander A. Papin^
Altai State University Barnaul, Russian Federation
Received 19.06.2020, received in revised form 16.10.2020, accepted 02.12.2020 Abstract. The dynamics of a snow-ice cover is considered within the theory of poroelasticity. The snow-ice cover is modeled by a three-phase medium consisting of water, air and ice. The governing equations are the equations of mass conservation for each phase with phase transitions, the equations of conservation of phase momentum in the form of Darcy's law, the equation of conservation of momentum of the whole system, the rheological equation for porosity and the equation of heat balance of snow. In the full formulation the liquid and air pressures are functions of the temperature and the corresponding densities, and the viscosity and compressibility coefficients of ice are functions of the temperature. The problem of two-dimensional nonstationary filtration of water in a thin poroelastic ice plate is considered in the model case. The solution is obtained in quadratures.
Keywords: poroelasticity, porosity, snow-ice cover, thin layer.
Citation: M.A. Tokareva, A.A. Papin, Mathematical Model of Fluids Motion in Poroelastic Snow-ice Cover, J. Sib. Fed. Univ. Math. Phys., 2021, 14(1), 47-56. DOI: 10.17516/1997-1397-2021-14-1-47-56.
A large number of studies are devoted to balance models and observations of the behavior of a snow cover [1-3]. These works primarily study snow stocks, flood risks and river flow formation. Models of a melting snow-ice cover based on the approaches of the mechanics of multiphase medium are intended to obtain a more detailed picture of the movement of water and air in melting snow, i.e. determination of the velocity field, concentration, phase pressures, etc. Only a few articles use the "full" models. The problems describing phase transitions and movement of the ice skeleton remain open.
The foundations of the theory of the movement of water and air in the melting snow within the mechanics of multiphase medium are formulated by S.C.Colbeck [4] and his followers [5,6]. In these works, snow is considered as a multiphase medium, but the variable porosity of ice, its deformation and phase transitions are not taken into account. In this paper, a mathematical model of the movement of water and air in a deformable poroelastic medium is proposed. The model includes the laws of conservation of mass, momentum in the form of Darcy's law for fluid components, the rheological relation for the porosity and the equation for temperature [7]. The obtained system of the equations is a generalization of the well-known system of Musket-Leverett equations in the case of variable porosity [8]. We use the concept of effective stress introduced by Terzagi for the porous deformation model when modeling a poroelastic medium. In this approach, the effective stress is defined as the difference between the total stress and the
* [email protected] [email protected] © Siberian Federal University. All rights reserved
pressure of the liquid phase [9]. This proposition means that the fluid carries part of the load. Fundamental in this approach is the relationship between the deformation of the skeleton of a solid matrix and the processes of fluid flow. Filtration in a thin poroelastic layer is considered in the case of one liquid. The model is based on the theory of multiphase flows in poroelastic medium. Structurally similar mathematical models are used to describe the movement of magma in the earth's crust, the flow of water and oil in a porous medium, the growth of tumors and other filtration processes in poroelastic medium [10-13].
1. Governing equations
A system of differential equations describing the movement of air and water in the snow-ice cover taking into account the phase transitions is considered in the domain Q = (x, z) = = [0, L] x [0, H]. The equations of continuity with the variable porosity $ are
-—--+ div((1 - $)pi Vi) = Iwi + Iai, -iTT--+ div(pw $Sw Vw ) = Uw + law ,
dt dt
d(PaSa$) , , , ^ S T , T , -dt--+ div(pw $SaVw) = Uiw + Uwa, Sa + Sw = 1.
Here pa, pw, pi, va, vi, vw are the true density and phase velocity, respectively (a — air, w — water, i — ice), sa,sw are the saturations of air and water, respectively, A is the porosity, Ilm is the intensity of mass transfer from /-phase to the m-phase per unit volume at time instant,
Ilm Iml.
Instead of the equation of dynamics in the theory of two-phase filtration, we use the general Darcy law for liquid phases, taking into account the movement of ice
k(A)
A(Vw - Vi) =--J— (Vpw - Pwg),
Pw (9)
A(( V V) k($)
A(Va - vi) =--TTTV iypa - Pag) ,
Pa(V)
where 9 is the temperature of the medium (9i = 9w = 9), pw ,pa are the water and air pressures, respectively, k(A) is the permeability, pw(9), pa(9) are the dynamic viscosities of water and air, g is the gravity acceleration vector. The phase pressures pa and pw differ by the magnitude capillary jump pc: pa - pw = pc(x, Sw)■
The Maxwell-type relation between the porosity and the effective pressure pe is [17]
divvi = -*(a(9)pe + ¡3(9)dpf) , d = di +(Vi • V),
where a(9), ¡(9) are the specified environment parameters. The effective pressure pe and the pressures of the liquids pw, pa and solid pi phases are related by the equations
ptot = Apf + (1 - A)pi, pe = (1 - A)(pi - pf ), pf = Sw pw + Sapa■
The balance equation for the forces of the whole system has the form Vptot = Ptotg+div(^(1 - A)v(9)(^dXi + (jx^ ) ), Ptot = APf + (1-A)Pi, Pf = Sw Pw +SaPa,
where ptot is the total pressure, ptot is the total density, n(0) is the viscosity of the porous skeleton, * is the symbol of the transposition operation.
The energy conservation equation is taken as [1]
d0
(PaCaSa 0 + Pw Cw Sw 0 + PiCi(l - 0)) — + (paCaSa&Va + Pw CWSW 0Vw + Pi Ci(l - 0)Vi)V0 =
= div (Atot(0)V0) - Lilai - Vlwi.
Here ci = const > cw = const > 0, ca = const > 0 are the heat capacity of ice, water and air at constant volume, respectively; Li = const > 0 is the specific heat of ice sublimation, v = const > 0 is the specific heat of ice melting; Xtot is the thermal conductivity of the medium as a whole (Atot = atot + btotPlot, atot = const > 0, btot = const > 0). To close the equations, one can take into account the dependence of the densities of liquid phases on pressures and temperatures, or set the densities equal to constant.
Assuming that there is no air in the pores, we shell study the system of differential equations describing the filtration of liquid in the melting ice (the phase densities are constant):
d(l - 0)Pi . ^ T d(Pw sw0) , , , ^ s T (1)
-—--+ div((l - 0)PiVi) = lwi, ---+ div(Pw0swVw) = hw, (1)
dt dt
k(0)
0(Vw - Vi) =--T^T (VPw - Pwg) , (2)
Pw (0)
divVi = -t(a(0)pe + ¡(0) d^j , = dt + ( Vi •V), (3)
Ptot = 0Pw + (l - 0)Pi, Pe = (l - 0)(Pi - Pw), (4)
Vptot = Ptot9 + div ^(l - 0)n(0) ^dx IX^^) ), Ptot = &Pw + (l - &)Pi, (5) d0
(Pw cw 0 + Pici(l - 0)) dt + (Pw cw&Vw + Pici(l - 4>)Vi)V0 = div (Atot (0)V0) - vlwi. (6)
A structurally similar model is used to describe the nonisothermal motion of magma in a poroelastic medium. The system (1)-(6), in the one-dimensional case with constant temperature, the dependence of the density of the liquid phase on the pressure, and in the absence of the phase transitions is studies and the local solvability is proved in [18]. The global solvability is proved in [19] for the same problem when the densities are constant. The numerical analysis of the initial-boundary value problem for system (1)-(6) in the isothermal case is carried out in [20].
2. Introduction of small parameter
We shell make equations (1)-(6) dimensionless. Let x, Z, t be dimensionless variables defined by the equalities [21]
x _ z i H
x =—, Z = —, t= £ Tot, £ = — C l,
L H L
where [L] = [H] = [ m], To] = [l/ s]; l is a fixed parameter.
The required functions in the dimensionless variables are
Pw (t,x,z) = PPw (t,x,z) = PPw(^£lTot, Hj,
pi(t, x, z) = ppi(t, x, z) = ppi(el Tot, -Hj,
vj(t,x,z) = v°vi(i,x,z) = v°vj ^'''o^ ^ H), j = 2 vw(t,x,z) = vivw(z,x,z) = vivw(£lтot,Jl,~H), j = 1,2
ptot(t,x,z)= pptot(t,x,z), pe(t,x,z)= ppe(t,x,z), 9 = 59, Iwi = YIwi, Iiw = ihw ■
Here [vi] = [m/ s], [p] = [Pa], [5] = [K], [a] = [1 / c], the quantities p, vj,5, y are the positive constants. Further, we use the assumption that k(A) = kAn, k = const, n(9) = n = const > 0 [14,17,22].
Equations (1)-(6) in scalar form in the dimensionless variables are
l d{1-A)^ V1 dv1(1 - A) , V2 dv2(1 - A) z + L dx + H dz =YIw
l . v1 dv w$ , vv d^w$ z £T0 dt + T^ + H= YIiw,
1/ 1-1 1-1n k , n p dpw A(v vw - v V1) =--A Y^—,
p L Ox
1/2-2 2-2\ k i n ( p dpw \
A(v vw - v v i ) = -- $ [nit + Pw
p \H dt
dvti ( _ afl dpe 1 pv1 dpe 2 pv2 dpe
2T2 dx ((1 - A) §) + H2 | {(1 - A)§) +
L2 dx
+—— ((1- A)dttl\ = p dptot + HLdt V A) dx ) Ln dx
(Pw A + Pj (1 - A))g + v2 d { dt2\ +2 v2 d { dv2\ +
-n-Ptt + L(x <v(1 - A)-ax) +2n dt - A)dt) +
+iLd f(i A)dvT\ = p dptot
+ HL dx V( A) dt J Hn dt ,
l a 1 11 ,wd9 v15 , d9
£ To5 (PwCw A + PiCi(1 - A)) dt + L P CwAv w + PiCi(1 - AM) dx +
v25 ( t 2 2) d9 5 d (. ,, M ,d(t)
+ H \Pw Cw Avw + PiCi(1 - A)vi) dt = J2 dx I (a + b(APw + (1 - A)Pi) dx) +
5 d / 00\
+ H2 dt ( (a + b(APw + (1 - A)Pi)2) dt ) - VYIwi■
We now choose
v1 = £lT0L, v2 = £l T0H■ Then equations (7)-(12) can be reduced to the form
d(1- AA + div((1 - A)v i) = -^Li,
dt £l T0
^ + div(Av w) = ^r—Iiw,
dtt £lT0
(7)
(8)
Lix + Yiy = -A{«ppe - 3 {p^n+vi^d!+vino)), (9)
(10)
(11)
(12)
kp
£l+2T0^L2
£lT0pL2 - dPw
0(vw - vi ) = -0 -dX'
kp
JJ-2 -2\ in dpw
0(vw - vi ) = -0 ~dzr -
Pw gL£
£1tq . dv— ap dx
dv[ dz
" 'u,^ ^i . "^i , _ " 'UK i " re . -1 dpe -2 dpe
TP^^ + d") = ( dV + vidX + VidZ
£lTuft ( dpe
a \ dt
dz
'
(14)
(15)
^ d (" - S)+£-2 K^ - « §)+
£' ff((i - *)^
d z \ d z
__p dptot
dx J nTu dx '
(16)
£,+2 d (a - *) i)+2£i did - *) §)+
dx
gLpw VTu
■0 +
dz
gL£pi VTu
(1 - 0) +
p dptot
nTu dz
(17)
p
l( (Pw Cw x , n dd (PwCw,-1 . ,, de , (PwCw, -2 . c d0\
0 + (1 - 0))( H0Vw +(1- 0'v) dx + (p)0Vw + (1 - 0)v') dt)
=d(t a + ((-1) 0+1)2) S)+
dx \\L2TuPiCi L2TuCi \\ p.
1 d (( a b(pw - Pi) b \ dd\ V
+ dV{{lt^ + -LTUPC-0 + LTUCJ eV) - vlIwi' (18) 3. The passage to the limit and solution by quadratures
For small £, £ ^ 0, the form of equations (13)-(18) depends on the parameter l. We consider the most physical case (for l = -2) describing a slow flow of the fluid in the poroelastic medium (below, the dashes are omitted):
d I_
dz
d(1- 0) + div((1 - 0)vi)=O, d0 — + div(0vw) = 0' dt (19)
vi1 = vw1 ' (20)
TU^L\f 2 ^ Xn dpw kp 0(vw vi ) = 0 dz' (21)
dvl + dv1 = ftpd(ptot - pw) dx dz dt ' (22)
i (d- 0)£) =«■ (23)
(1 - 0)§) + S ((1 - 0)£) =0' (24)
1 "'a+b^ - O+fl dt)-
2N -N (25)
L2ToCi dz \\ Pi \ VP
We outline the scheme of the method of the solution for system (19)-(24). Some of the equations are integrated and as the result of some transformations we arrive at an equation for the function A.
At the first stage, from equation (19) we find divvi and substitute it into equation (22), then
dG(A) , ^ T7niX\ d(ptot - pw) , ^ vy/ N
+ vi • VG(A) = -k---+ vi • V(ptot - pw),
where
dt 1 dt
dG(A) 1
dA 3pa(1 - A)'
From the last equation we get the equation for pi (if A and pw are found):
dU ! dU 2 dU
H + vi dx + vi dz = °,
where U = G(A) + (1 - A)(pw - pi).
This equation is solved by the method of characteristics, which are defined as the solution of the Cauchy problem
v! = dt, v2 = U \t=o= G(A0) + (1 - A0)(p0w -pO)
At the second stage, we find a representation for the velocity components. To do this, we twice integrate equation (23) over z and, taking into account (20), we obtain
rz 1
v\ = vw = A(x,t) --- dT + B(x,t), (26)
Jo 1 - A
where the functions A(x,t), B(x,t) are arbitrary functions of (x,t). Equation (24) can be written, using (26), as
2^ ((1 - A)di) + ^^
Integrating the last equation over z twice, we obtain a representation for vi2 in terms of A and arbitrary functions C(x,t), D(x,t):
2 1 dA(x,t) fz t , ^, , [z 1
v2 — v ' ' I ^ '
z T z 1
---dT + D(x,t) --- dT + C (x,t).
Jo 1 - A Jo 1 - A
2 dx J0 1 - A ' J0 1 - A
At the third stage, we find representations for the liquid phase. To do this, we summarize equations (19):
div(A(vw - vi)+ vi) =0^
Using equations (20) and (21), we have
^ kp d I ^pw
divvi = -A -TT-
T0pL2 dz
Substituting into the last equation the representation for the velocity components of the solid medium, we get the equation for pw:
U™ £ ih *) + ^ - '/2
_i_nf ^ 1 kpd, upw
+D(x,t)--- = —$
1 — $ ropL2 dz
Integrating this equation twice over z, we obtain a representation for pw in terms of $ and arbitrary functions E(x,t),F(x,t)
r0L r
Pw = —,- $
kP Jo
[ $-n[ jT ¿(^ i—-$ *)«-
1 dA(x,t) rc T , dB(x,t) , , fc 1 , ^ ^ .
i-jdT +—C + D(x,t) ---dT +E(x,t)] dZ + F(x,t).
Jo 1 — $ dx Jo 1 — $
2 dx Jo 1 — $ dx 'Jo 1 — $
Using a representation for pw, from (21) we derive the equation for the second component of the velocity v2w in terms of $
2 1
Vw = — $
I ¿HZ ^ dT^— 2 r—-$dT+
dB(x,t) . fz 1 , ^ .
+ —z + D(x,t) --- dT + E (x,t)
dx Jo 1 — $
dB(x,t) t , fz 1
~ ' ^ t /o 1—$
r 1
+D(x,t) --- dT + C(x, t).
1 dA(x,t) fz t
2 ~jo 1—$ dT+
Jo 1 — $
Thus, the velocity components vi,vw and the pressure pw are obtained as the functions of $.
At the last stage, the continuity equation for the solid phase is used to determine the porosity $:
' " ' ' s-ldT + sB(x, t) I +
ds d ( ¡'z \
- + gx[SA(x,t) Jo ^ + SB(x't))
+ d (—s1 dA(x,t) f TdT + sD(x,t) i s-id,T + sC(x,t))=0, s = 1 — $. dz V 2 dx Jo s Jo J
Note that the representation for pi is found after obtaining the representation for the components of the velocity of the solid phase and the pressure of the liquid.
Now we consider the following initial-boundary value problem for equations (19)-(24). At z = 0, the rigid porous skeleton satisfies a plane stress state with the stress tensor aw = a22 = ptot, al2 = a2l = 0, and the fluid pressure is also set. The velocity of the rigid skeleton is specified at z = H, and the phase velocities are equal. The porosity and phase pressures are set at the initial moment of time. Thus, the following initial boundary conditions are considered:
dv2
\z=o= 0, V2 \z=H = C = const,
QV1
V1 \z=H = B = const, \z=o = 0, (27)
$ \t=o= $o(x,t), \z=H =0,
pw \t=o= p°w(x,z), pi \t=o= po(x,z), pw \z=o= po(x,t),
where the constants B and C are given. This allow us to uniquely define the following functions:
A(x, t) = 0, B(x, t) = B = const, C(x,t) = C = const, D(x,t) = 0, E (x,t)=0, F (x,t) = p0(x,t),
v1 = vw = B, vi = vw = C\ pw = Po(x, t),
and the equation for S can be rewritten as
— + B— + C— =a
dt dx dz
The solution of this equation is
A = A0(x - Bt,z - Ct).
The solution for pi is
pi = po(x,t) - p0w (x - Bt,z - Ct)+ pO(x - Bt,z - Ct). The equation for 9 takes the form
A0(x - Bt,z - Ct) + (1 - A0x - Bt,z - Ct») (d9 + Bdx + C=
= ^ + bP' ^ -Bt■ z -Ct) Ppw+1 - A°'x - Bt-z -a?)) S)
and is supplemented by the following initial-boundary conditions
d9 d9
9 \t=o= 90(x,z), — \z=o,z=H =0, dx \x=o,x=L=0. (28)
The equation for 9 is linear uniformly parabolic. The maximum principle for the function 9 is valid, and theorem of existences is also established for this problem (q.v. [23], pp.29, 364). Thus, the solution of problem (19)-(25), (27), (28) is found.
Conclusion
A new model of the dynamics of the snow-ice cover, taking into account the deformation of the porous iced skeleton, is considered in this work. The model is a generalization of the classical Musket-Leverett two-phase filtration model and it is applicable for a wide class of problems of multiphase flows in a poroelastic media. A feature of the model is the variable porosity. The problem of the justification of this model is open. There are also no results on the construction of simple solutions. On the other hand, the study of such a general problem is important for the analysis of complicated processes of salt-mass transfer in a melting snow. For example, the speed of pollution in spring during melting of the snow cannot be described without knowing the field of velocities and water saturation. Using a sophisticated model allows to get a complete picture of the dynamics of the process. The important point is the correct choice of initial data, which guarantees the existence of a solution. Only examples of single-phase numerical solutions are known without phase transitions and temperature effects [14,22]. In this article, we consider a simplified model based on the selection of a small parameter, which is a common technique in
the study of complex problems. Solutions in quadratures are constructed. In the future, we plan to study the justification for the formulated problem and its applications in polar mechanics.
The work was carried out under the project MK-204-2020.1 "Initial-boundary value problems for the equations of fluid motion in poroelastic media and their applications in the dynamics of snow and ice cover" with the support of a grant from the President of the Russian Federation.
References
[1] L.S.Kuchment, P.Romanov, A.N.Gelfan, V.N.Demidov, Use of satellite-derived data for characterization of snow cover and simulation of snowmelt runoff through a distributed physically based model of runoff generation, Hydrol. Earth Syst. Sci., 14(2010), 339-350.
DOI: 10.5194/hess-14-339-2010
[2] A.N.Gelfan, J.W.Pomeroy, L.S. Kuchment, Modeling Forest Cover Influences on Snow Accumulation, Sublimation, and Melt, Journal of Hydrometeorology, 5(2004), 785-803.
[3] M.Sturm, J.Holmgren, G.Liston, A seasonal snow cover classification-system for local to global applications, Journal of Climate, 8(1995), no. 5, 1261-1283.
DOI: 10.1175/1520-0442(1995)008<1261:ASSCCS>2.0.C0;2
[4] S.C.Colbeck, A theory of water percolation in snow, Journal Glaciol, 11(63)(1972), 369-385.
[5] J.M.N.T.Gray, Water movement in wet snow, Philosophical Transactions: Mathematical, Physical and Engineering Sciences, 354 (1707)(1996), 465-500.
[6] S.Sellers, Theory of water transport in melting snow with a moving surface, Cold Regions Science and Technology, 31(2000), no. 1, 47-57. DOI: 10.1016/S0165-232X(00)00006-9
[7] R.I.Nigmatullin, Dynamics of multiphase media, Vol. 1, Moscow, Nauka, 1987 (in Russian).
[8] S.N.Antontsev, A.V.Kazhikhov, V.N.Monakhov, Boundary-Value Problems of the Mechanics of Inhomogeneous Fluids, Novosibirsk, Nauka, 1983 (in Russian).
[9] K.Terzaghi, Theoretical Soil Mechanics, New York, Jhon Wiley, 1943.
[10] V.V.Vedernikov, V.N.Nikolaevsii, Mechanics equations for porous medium saturated by a two-phase liquid, Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, 5(1978), 769-773.
[11] Y.Geng, L.Zhang, Bifurcations of traveling wave solutions for the magma equation, Applied Mathematics and computation, 217(2010), 1741-1748. DOI: 10.1016/j.amc.2009.11.035
[12] G.Simpson, M.Spiegelman, M.I.Weinstein, Degenerate dispersive equations arising in the study of magma dynamics, Nonlinearity, 20(2007), 21-49.
[13] A.Friedman, Cancer as Multifaceted Disease, Math. Model. Nat. Phenom., 7(2012).
[14] J.A.D.Connolly, Y.Y.Podladchikov, Temperature-dependent viscoelastic compaction and compartmentalization in sedimentary basins, Tectonophysics, 324(2000), 131-135.
[15] J.Bear, Dynamics of Fluids in Porous Media, Elsevier, New York, 1972.
[16] A.A.Papin, Y.Y.Podladchikov, Isothermal motion of two immiscible fluids in poroelastic medium, Izvestiya Altai State University, 1/2(85)(2015), 131-135 (in Russian).
DOI: 10.14258/izvasu(2015)1.2-24
[17] A.Fowler, Mathematical Geoscience, Springer-Verlag, London Limited, 2011.
[18] A.A.Papin, M.A.Tokareva, On Local Solvability of the System of the Equations of One Dimensional Motion of Magma, Journal of Siberian Federal University. Mathematics and Physics, 10(2017), no. 3, 362-371. DOI: 10.17516/1997-1397-2017-10-3-362-371
[19] M.A.Tokareva, A.A.Papin, Global solvability of a system of equations of one-dimensional motion of a viscous fluid in a deformable viscous porous medium, Journal of Applied and Industrial Mathematics, 13(2019), no. 2, 350-362. DOI: 10.33048/sibjim.2019.22.208
[20] M.N.Koleva, L.G.Vulkov, Numerical analysis of one dimensional motion of magma without mass forces, Journal of Computational and Applied Mathematics, 366(2020), 112338. DOI: 10.1016/j.cam.2019.07.003
[21] J.Escher, M.Hillairet, P.Laurencot, C.Walker, Thin film equations with soluble surfactant and gravity: modeling and stability of steady states, Mathematische Nachrichten, 285(2012), 210-222.
[22] C.Morency et al., A numerical model for coupled fluid flow and matrix deformation with applications to disequilibrium compaction and delta stability, Journal of Geophysical Research: Solid Earth, 112(2007), no. B10.
[23] O.A.Ladyzhenskaya, V.A.Solonnikov, N.N.Ural'tseva, Linear and Quasilinear Equations of Parabolic Type, Moscow, Nauka, 1967 (in Russian).
Математическая модель движения жидкостей в пороупругом снежно-леДовом покрове
Маргарита А. Токарева Александр А. Папин
Алтайский государственный университет Барнаул, Российская Федерация
Аннотация. Задача динамики снежно-ледового покрова рассматривается в рамках теории по-роупругости. Снежно-ледовый покров рассматривается как трехфазная среда, состоящая из воды, воздуха и льда. В основу математической модели положены уравнения сохранения массы для каждой из фаз с учетом фазовых переходов, уравнения сохранения импульсов фаз в форме законов Дарси, уравнение сохранения импульса системы в целом, реологического уравнения для пористости и уравнение теплового баланса снега. В полной постановке учтена зависимость давлений жидкости и воздуха от температуры и соответствующих плотностей, а также зависимость коэффициентов вязкости и сжимаемости льда от температуры. В модельном случае рассматривается фильтрация воды в тонкой пороупругой ледовой пластине. Получены решения в квадратурах.
Ключевые слова: пороупругость, пористость, снежно-ледовый покров, тонкий слой.