CC BY
45
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
The Problem of a Longitudinal Crack with a Filler in a Strip
N. N. Antonenko
Antonenko Nina Nikolayevna, Zaporozhye National Technical University, 64, Zhukovsky st., 69063, Zaporizhzhya, Ukraine, [email protected]
The method of the solution the problem of the central longitudinal crack with a filler in a strip is proposed. It is assumed that the jumps of the components of displacement vector is proportional to the corresponding stresses at its upper edge. Fourier’s method of integral transformation is used. The problem is reduced to a system of integro-differential equations. The effects of influence of thickness, mechanical properties of a strip and a filler of the crack on Mode I and Mode II stresses intensity factors (SIFs) are examined. The following conclusions are made: increase the width of a strip and the coefficient that characterizes a filler leads to reduction of SIFs; increase of shear modulus and Poisson’s ratio of a strip leads to the increase of SIF.
Key words: strip, crack, filler, stress intensity factors, method of Fourier integral transformation, integro-differential equation.
References
1. Fiehter W. B. Stresses at the tip of a longitudinal crack in a plate strip. Washington, National Aeronautics and Space Administration, 1967, 55 p.
2. Aleksandrov V. M., Smetanin B. I. Equilibrium crack in a thin layer. J. Appl. Math. Mech. [Prikl.
Mat. Mekh.], vol. 29, iss. 4, 1965, pp. 926-929.
3. Smetanin B. I. Nekotorye zadachi o shheljah v uprugom kline i sloe [Some of the problem of cracks in an elastic wedge and layer]. Inzh. zh.
MTT, 1968, iss. 2, pp. 115—122 (in Russian).
4. Savruk M. P. Dvumernye zadachi uprugosti dlia tel s treshchinami [Solution of two-dimensional problems of the theory of elasticity for bodies with cracks]. Kiev, Nauk. dumka, 1981, 324 p. (in Russian).
5. Aleksandrov V. M., Smetanin B. I. A longitudinal crack in a prestressed thin elastic layer with free boundaries. J. Appl. Math. Mech. [Prikl. Mat.
Mekh.], 2005, vol. 69, iss. 1, pp. 141-150.
6. Alexandrov V. M. Longitudinal crack in an orthotropic elastic strip with free faces. Mech.
Solids [Izvestiya Akademii Nauk. Mekhanika Tverdogo Tela], 2006, vol. 41, iss. 1, pp. 88-94.
7. Pozharskiy D. A., Molchanov A. A. Asymptotic
УДК 532.591
НЕЛИНЕЙНЫЕ ВОЛНЫ НА ПОВЕРХНОСТИ СЛОЯ ВЯЗКОЙ ЖИДКОСТИ
К. Ю. Басинский
Басинский Константин Юрьевич, кандидат физико-математических наук, доцент кафедры математического моделирования, Тюменский государственный университет, [email protected]
Рассмотрена нелинейная задача о распространении волн по свободной поверхности слоя вязкой несжимаемой жидкости бесконечной глубины в плоском случае. С помощью метода малого параметра данная нелинейная задача раскладывается на задачи в первых двух приближениях, которые последовательно разрешаются. Получены нелинейные выражения для компонент вектора скорости, динамического давления и формы свободной поверхности. Изучается движение частиц вязкой жидкости, вызванное распространением волны по свободной поверхности. Установлено, что вязкость жидкости
solutions of mixed problems for elastic strip and wedge. Vestnik of DSTU [Vestnik of DSTU]. 2010, vol. 10, pp. 447-454 (in Russian).
8. Murakami Y. Stress intensity factors handbook. Pergamon Press, 1987, vol. 1, 1566 p.
9. Antonenko N.N. Modelirovanie mezhfaznoi treshchiny s napolnitelem na granitse uprugoi polosy i uprugoi poluploskosti [Modeling of a crack with a filler on the elastic strip-half-plane interface] Visnik Donets’kogo natsional’nogo universitetu. Ser. A : Prirodnichi nauki. 2013, no. 1, pp. 23-27 (in Russian).
10. Tkachenko I. G. Dvomirna mishana zadacha termopruzhnosti dlia bagatosharovoi osnovi [A two-dimensional mixed thermoelasticity problem for a multilayer foundation] Prikladni problemi mehaniki i matematiki [Applied problems of mechanics and mathematics]. 2005, vol. 3, pp. 7078 (in Ukrainian).
11. Alexandrov V. M., Pozharskii D. A. To the problem of a crack on the elastic strip-half-plane interface. Mech. Solids [Izvestiya Akademii Nauk. Mekhanika Tverdogo Tela], 2001, vol. 36, iss. 1, pp. 70-76.
© Басинский К. Ю., 2014
К. Ю. Басинский. Нелинейные волны на поверхности слоя вязкой жидкости
оказывает существенное влияние на форму траекторий жидких частиц, которое проявляется как в уменьшении амплитуды колебаний с течением времени, так и в отличии траекторий вблизи свободной поверхности и при заглублении. Исследован нелинейный эффект Стокса, который заключается в наличии приповерхностного течения.
Ключевые слова: волновые возмущения, вязкость, траектории частиц.
DOI: 10.18500/1816-9791 -2015-15-3-322-329
ВВЕДЕНИЕ
Решение нелинейной задачи о волнах на поверхности слоя идеальной жидкости приведено в работах [1,2]. В работе [2] также найдены нелинейные траектории жидких частиц. Для вязкой жидкости известно решение линейной задачи [3-5], а нелинейная задача решена только для приближения слабовязкой жидкости [6]. Данная работа посвящена определению нелинейных выражений для скорости, давления и формы свободной поверхности при волновом движении вязкой жидкости, а также определению траекторий жидких частиц.
1. НЕЛИНЕЙНАЯ ЗАДАЧА О ВОЛНАХ НА ПОВЕРХНОСТИ ВЯЗКОЙ ЖИДКОСТИ
Рассматривается слой вязкой несжимаемой жидкости бесконечной глубины. Свободная поверхность слоя граничит со средой пренебрежимо малой плотности, характеризующейся постоянным давлением Pa (в частности, атмосферным). Декартовая система координат задана так, что плоскость z* = 0 совпадает с невозмущённой поверхностью, а ось z* противоположно направлена к вектору силы тяжести g. Движение жидкости происходит в плоскости x*z* со скоростью u* = (u* (t*, x* ,z*), 0, v* (t*, x* ,z*)). Звездочкой, где это необходимо, обозначены физические (размерные) величины.
Пусть в положительном направлении оси распространяется волна длины А. Длина волны много больше ее высоты (А >> £**ax). В области, занятой жидкостью выполняются уравнения неразрывности и движения:
div u* = 0,
д u*
+ (u* V)u*
-1 Vp* + P Au*. P P
(1)
Здесь p* = P — Pa + pgz* — динамическое давление, p — плотность, P — давление, p — коэффициент динамической вязкости.
На свободной поверхности z * = £ * (t* ,x*) задаются кинематическое и динамическое условия:
v
*
pnn
д£*
" д*
— Pa,
+ u*
д£ *
дХ* ’
pnT -
0.
(2)
При бесконечном заглублении скорость жидкости должна затухать, т. е. выполнено условие
u* ^ 0, z* ^ —то.
(3)
Система уравнений (1) и граничных условий (2), (3) является замкнутой и составляет нелинейную краевую задачу для определения характеристик волнового движения.
Введем следующие безразмерные переменные и величины:
u* = zco u, p* = epc^ p, £ * = г£/k, vo = pk/pc0,
t = kct*, x = kx*, z = kz*, a = c/c0 = w/w0, c0 = g/k,
где c0 и ш0 — фазовая скорость и частота волны линейной задачи для идеальной жидкости соответственно, c и ш — фазовая скорость и частота волны, г = k£**ax — малый волновой параметр, k = 2п/А — волновое число.
В безразмерных переменных задача (1)-(3) принимает вид
div u = 0,
д u
a—----v0Au + Vp
дt
—e(uV)u,
(4)
д£ д£
V — аЖ = ^ z = г£’
Механика
323
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
dv (дu dv\ дf ofdf \ f dv
p - f - 2vqdz = -£V4dZ + d-x)dX +e fe) (f -p - 2vqdz
z = ef,
Vo | — + ^)|1 - e2( ^ ) |= -4vo'— z = ef,
du dv
dz + dx
dx u —— 0,
dv df
' dz dx
— -to.
В силу малости волнового параметра е граничные условия на свободной поверхности z = ef разложением в ряд Маклорена входящих в них функций, сводятся к условиям на фиксированной поверхности z = 0.
Решение задачи находим в виде рядов по параметру е:
ОО ОО ОО
u =^2 ei-1ui, un = (и*, 0, vi), p =^2 ei-1 pi, f = J^ ei-1 fi.
i = 1
= Ъ e
=1
i=1
Подставляя эти ряды в уравнения (4) и разложенные в окрестности нуля граничные условия и приравнивая коэффициенты при одинаковых степенях е, получим задачи в первых двух приближениях.
В первом приближении задача имеет вид при е0
div U1 = 0, - vqAu1 = -Vp1,
dt
(5)
v1 =a-
df1 dt ’
dv1
P1 - f1 - 2vq— = 0,
d U1 + dvx = 0 dz + dx = ’
=0,
для второго приближения: при е1
div U2 = 0,
U1 —— 0, z —— —to.
d
a-Q^ - vqAu2 = -Vp2 - (u1 V)u1,
(6)
dfo + d 2 dvo d (2 dv1 \
vo=a~m + dx(u1f1) ■ P2 - f2 - 2V0 ЦТ =f1 OX rQ-dT - P1J
dU2 + dv2 = 4dU1 df1 __ f _d (dU1 + dvA
dz + dx = 4 dx dx f 1 dz V dz + dx Г
z = 0,
z = 0,
dz dz dx u2 — 0, z — -to.
Решение линейной задачи (5) имеет вид [5]
u1 = Ae- “1 [ez cos (x - t) + (aV1 + bV2) cos (x - t) + (aV2 - bV1) sin (x - t)] , v1 = Ae- a 1 {ez sin (x - t) + V1 cos (x - t) + V2 sin (x - t)} ,
(7)
p1 = Aez a 1 (a cos x + в sin x), f1 = 2 ^ 2 e a 1 (a cos x + s sin x),
a2 + s2
V1 (z) = ebz (B1 cos az - B2 sin az), V2 (z) = ebz (B2 cos az + B1 sin az),
A
- в t
B1 = -2voaA/ (a2 + s2), B = -2vo sA/ (a2 + s2).
Здесь в — безразмерный декремент затухания (в^0 — размерный), A — амплитудный параметр, s = 2v0 - в. Параметры a и b связаны соотношениями a2 = b2 - 1 + в/v0, 2ab = a/v0.
Частота волны через декремент затухания выражается следующим образом:
a
2 = 1 + s2 - 4vn/s.
Для декремента затухания получено уравнение
s6 + s4 - 4vQJs2 - 4vQ| = 0,
аналитическое решение которого найдено, но не приведено здесь из-за своей громоздкости.
324
Научный отдел
К. Ю. Басинский. Нелинейные волны на поверхности слоя вязкой жидкости
Подставив выражения (7) в задачу (6), получим систему линейных неоднородных уравнений и граничных условия для определения неизвестных функций u2, v2, p2 и £2, решение которой имеет вид
U2 = e “ * ({e2zD2 + 1/2 [ai V3 + biV4 + aV5 + (b +1) V6]} cos (2x — 2t) +
+ {— e2zDi + 1/2 [ai V4 — b\V3 + aV6 — (b +1) V5]} sin (2x — 2t)) ,
V2 = e-~a1 [(e2zDi + V3 + V5) cos (2x — 21) + (e2zD2 + V4 + V3) sin (2x — 2t)\ ,
P2 = e-^* { [e2z (eDi + aD2) + Vr\ cos (2x — 2t) + [e2z (PD2 — aDi) + Vs\ x x sin (2x — 2t) — ez AV2 + AB2 + A2 [(2vo — @s) j{a2 + s2) — e2^] /2 + + (Bi + B|) (1 — e2^) /2} ,
e-2в t
p2 = 2 (a2 + в2) {[a (D2 + B4 + K2 — Q2) — в (D1 + B3 + Ki + Qi)] cos (2x — 2t) —
— [в (D2 + B4 + K2 — Q2) + a (Di + B3 + Ki + Qi)] sin (2x — 2t)} ,
где
V3 (z) = ebl z (B3 cos ai z — B4 sin ai z), V4 (z) = ebl z (B4 cos ai z + B3 sin ai z),
V5 (z) = e(b+i)z (Ki cos az — K2 sin az), V6 (z) = e(b+i)z (K2 cos az + Ki sin az),
V7 (z) = e(b+i)z (Gi cos az — G2 sin az), Vs (z) = e(b+i)z (G2 cos az + Gi sin az),
a2 = bi — 4 + 2в/vo, ai bi = a/vo,
Ki = {[(b + 1) Ri + aR2] Ci — [aRi — (b +1) R2] C2}/Ai, Ai = (C? + C22)/2,
K2 = {[aRi — (b +1) R2] Ci + [(b +1) Ri + aR2] C2}/Ai,
Gi = {[aa + в (b — 1)] Ki — [a (1 — b) + ea] K2 + 2R2}/4,
Gi = {[a (1 — b) + ea] Ki + [aa + в (b — 1)] K2 + 2Ri}/4,
Ci = [4av0 — a (в + 2vo)\ j(2vo), Ci = [8bv2 — a2 + в0 — 4vos\ j(4vo),
Ri = A[(2bvo — s) Bi + (a — 2avo) B2]/(2vo),
R2 = A[(a — 2avo) Bi — (2bvo — s) B2]/(2vo),
B3 = vo [si (L2 — 4Di) + a (Li — 4D2)]/A2, A2 = s2 + a2,
B4 = vo[si (Li — 4D2) — a (L2 — 4Di)]/A2, si = 4vo — в,
Li = A2 {a [4v^ (a2 — s2)/A3 — s — Ki\ + ab [1 — K2 + 8vos/A3] — 3a+
+ [(в — 6vo) K2 — aKi]/(2vo)}/A3, A3 = s2 + a2,
L2 = A2 {b [4v2 (a2 — s2) j A3 — s — Ki \ — aa [1 — K2 + 8vo s/A3 ] + 3s+
+ [(в — 6vo) Ki + aK2]/(2vo)}/A3,
Di = (Fi Ji + R2 J2) / A4, D2 = (Fi J2 — F2 Ji) / A4, A4 = Ji2 + Ji2,
Ji = [8v2 (aai + sibi) — si/A2] — si, J2 = [8v2 (abi — si ai) — a/A^ + a,
Fi = A2 [2vo (a — s)/A3 + 1] [2vo (a + s)/A3 — 1]/2 + /i Ki + GK — Gi + Hi Li +
+H2L2 + Wi — aQ2)/(2A5), A5 = (a2 + в2) ,
F2 = 2avoA2(1 — 2svo/A3)/A3 + /2Ki — 1iK2 + G2 — H2Li + HiL2 + (вQ2 + aQi)/(2A5),
Ii = 2bvo + (4vo — в^5 )/2, /2 = —2avo + a/(2A5),
Hi = [2av2 bi — 2si v2 ai + av^ s/A5]/A3,
H2 = [2av/ ai + 2si v/ bi + ^si — a2) j A5 \ j A3.
Механика
325
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
2. ВОЛНОВЫЕ ТРАЕКТОРИИ ЖИДКИХ ЧАСТИЦ
Физические координаты частицы x*(t*), z*(t*) удовлетворяют уравнениям
dx*
d*
dz
dt*
= v
Физическое время необходимо обезразмерить частотой колебания жидкой частицы а, так как для нелинейных волн даже в идеальной жидкости а не совпадает с частотой волны [2].
Учитывая, что
x* = ct* + x/k, z* = z/k, t* = t/а, u* = eco u,
уравнения, описывающие движение частицы в волне можно представить в виде
dx=<а-тУ- - И)'-
t ck
Положим
ГО
x = Гх„
i=0
dz
d = 4 a ■ 1 v'
kc
ХЧ xi, z = ^2 £i Zi, — = 1 + ^ £i Yi.
(8)
(9)
i=0
i=i
Подставляя ряды (8) и выражения для компонент скорости в уравнения движения (9) для определения первых трех коэффициентов, получим следующие уравнения:
dx о dt
= -1,
dzo
dt
= 0,
dxi = Ae a t {[ez0 + aVi (zo) + bV2 (zo)] cosxo + [aV2 (zo) - bVi (zo)] sinxo} - (tYi),
dt а
dzi A - a tj\azo
dx2 e"
dt
а
= — e a t {[ez0 + Vi (zo)]cos xo + V2 (zo) sin xo} , dt а
({e2z° D2 + 1/2 [ai V3 (zo) + biV4 (zo) + aV5 (zo) + (b + 1) Ve (zo)]} cos2xo+
+ {-e2zoDi + 1/2 [aiV4 (zo) - bi V3 (zo) + aVe (zo) - (b +1) V5 (zo)]} sin2xo) - (tY2)+
- в t
e a 1
+------({Aziezo + b (2azi - xi) Vi (zo) - [(a2 - b2) zi - axi] V2 (zo)} cos xo+
а
a
+ { -Axiezo + b (2azi - xi) V2 (zo) + [(a2 - b2) zi - axi] Vi (zo)} sin xo) ,
7 - M t
dz2 _ e a t
dt а
[e2zoDi + V3 (zo) + V5 (zo)) cos2xo + (e2zoD2 + V4 (zo) + Ve (zo)) sin2xo] +
e a
+------{[Axiezo + bzi Vi (zo) + (xi - azi) V2 (zo)] cos xo + [Aziezo + bzi V2 (zo) - (xi - azi) Vi (zo)] sin xo}
а
Разрешая последовательно данные уравнения, получим следующие выражения:
- в t
e a t
xi =
xo = xL - t, zo = zL, Yi = 0,
{[-eAezL - (ab + ea) Vi (zL) + (aa - ^b) V2 (zL)] cos (xL - t) +
A5
[-aAezL - (ab + ^a) V2 (zl) - (aa - ^b) Vi (zl)] sin (xl - t)} ,
,-в t
ea
zi = [(aAezL - eVi (zl) + aV2 (zl)) cos (xl - t) + (-,5AezL - № (zl) - aVi (zl)) sin (xl - t),
A5
x2 =
- 2в t
e a t
2A5
({-(aDi + eD2) e2zL - (0ai + a&i) V3 (zt)/2 + (aai - e>i) V4 (zt)/2+
= u*.
ГО
ГО
a
326
Научный отдел
К. Ю. Басинский. Нелинейные волны на поверхности слоя вязкой жидкости
+AeZL [(2a0a + а2b - в2b)/Дг> - 0/(2vo) + (а2 - в2)/Лг>] V (zl) +
+AeZL [(2аРЬ - а2a + p2a)/ДГ) + a/(2v„) - 2ар/Дг>] V2 (zl) -
- [а (b + 1) + ра] V, (zl)/2 + [аа - р (b + 1)] V6 (zl)/2} cos (2xl - t) +
+ {(pDi - «D 1) e2zL - (ра 1 + «bi) V4 (zl) - (oai - pbi) V3 (zl) +
+AeZL [(2аРа + а2b - p2b)/Д5 - Р/(2vo) + (а2 - p2)/Д,] V2 (zl) --AeZL [(2opb - а2a + p2a)/Д5 + «/(2vo) - 2аР/Д5] Vi (zl) -
- [а (b + 1) + pa] V6 (zl)/2 - [аа - p (b + 1)] V5 (zl)/2} sin (2xl - t)),
- t
Z2 = цД { KaD2 - pDi) e2zL + аУ4 (zl) - /Л3 (zl) + «Vs (zl) - pV (zl)] x x cos (2xl - 2t) + [- («Di + pD2) e2zL - «Vs (zl) - pVi (zl) - аК5 (zl) - pV6 (zl)] x x sin (2xl - 2t) + AeZL [aVi (zt) + (b +1) V2 (zl) + A2] + bebzL (B2 + B22)} ,
1 - e-¥t
Y2 = 777—л~Г [2«voA2e2zL + AeZL (а2 - p2 + 2аа&о) Vi (zl) + 2aAeZL x 4pvo Д51
x (vo + bvo - p) V2 + ае2^L (B^ + (2a2vo + 2vo - 3p) ,
где xl, zL — лагранжевы координаты.
Величина us = c - a/k представляет собой переносную скорость вдоль горизонтальной оси, Ее приближенное выражение имеет вид
us =c( 1 - ck) =c 0- г+k
Jc72.
(10)
Для примера построены траектории движения жидких частиц при малом (рис, 1) и Большом (рис, 2) значении коэффициента кинематической вязкости, Видно, что частицы вблизи поверхности движутся быстрее, чем частицы на глубине, что обусловлено наличием приповерхностного течения Стокса, которое характеризуется переносной скоростью (10), С увеличением вязкости амплитуда волновых возмущений и переносная скорость убывают быстрее, Кроме того, вблизи поверхности при Большой вязкости траектории наклонены в верхней точке в сторону движения, что не наблюдается на заглублении, Это можно объяснить влиянием вязких касательных напряжений на свободной поверхности,
Механика
327
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
Рис. 2. Траектории частиц жидкости при v = 2 ■ 10 3 м2/с
ЗАКЛЮЧЕНИЕ
Таким образом, получено асимптотическое решение нелинейной задачи с точностью до членов второго порядка по малому амплитудному параметру. Получены выражения для траекторий частиц жидкости, а также для скорости приповерхностного течения Стокса.
Библиографический список
1. Сретенский Л. Н. Теория волновых движений жидкости. М. : Наука, 1977. 816 с.
2. Алешков Ю. 3. Теория волн на поверхности тяжелой жидкости. Л. : Изд-во Ленингр. ун-та, 1981. 196 с.
3. Ламб Г. Гидродинамика. Л. : Гостехиздат, 1947. 928 с.
4. Ландау Л. Д., Лифшиц Е. М. Гидродинамика. М. : Наука, 1986. 735 с.
5. Баринов В. А. Распространение волн по свободной поверхности вязкой жидкости // Вестн. С.-Петерб. ун-та. 2010. Сер. 10. Вып. 2. С. 18-31.
6. Баринов В. А., Басинский К. Ю. Решение нелинейной задачи о волнах на поверхности слабовязкой жидкости // Вестн. С.-Петерб. ун-та. 2011. Сер. 10. Вып. 2. С. 9-16.
Nonlinear Wave on the Surface Layer of a Viscous Liquid
K. Yu. Basinsky
Basinsky Konstantin Yurjevich, Tyumen State University, 10, Semakova st., 625003, Tyumen, Russia, [email protected]
The nonlinear problem of the propagation of waves on the free surface of the layer of viscous incompressible fluid of infinite depth in the plane case. Using small parameter method, this nonlinear problem is decomposed into problems in the first two approximations that consistently allowed. Nonlinear expressions for the components of the velocity vector, the dynamic pressure and the shape of the free surface. The motion of particles of the viscous fluid caused by the spread of the wave of the free surface. It is found that the viscosity has a significant effect on the shape of the trajectories of liquid particles which is manifested as a decrease in the amplitude of oscillation over time, and in contrast to the paths near the free surface and penetration. The nonlinear effect of Stokes, which is the presence of near-surface currents.
Key words: wave disturbances, viscosity, particle trajectories.
328
Научный отдел
К. Ю. Басинский. Нелинейные волны на поверхности слоя вязкой жидкости
References
1. Sretenskiy L. N. Teoriya volnovykh dvizheniy zhidkosti [Theory of wave motions of fluid]. Moscow, Nauka, 1977, 816 p. (in Russian).
2. Aleshkov Yu. Z. Teoriya voln na poverkhnosti tyazheloy zhidkosti [The theory of waves on the surface of a heavy liquid]. Leningrad, Leningrfd Univ. Press, 1981, 196 p. (in Russian).
3. Lamb G. Gidrodinamika [Hydrodynamics]. Leningrad, Gostekhizdat, 1947, 928 p. (in Russian).
4. Landau L. D., Lifshits E. M. Gidrodinamika
[Hydrodynamics]. Moscow, Nauka, 1986, 735 p. (in Russian).
5. Barinov V. A. Wave propagation on the free surface of a viscous fluid Vestn. S.-Peterb. un-ta, 2010, vol. 10, no. 2, pp. 18-31 (in Russian).
6. Barinov V. A., Basinsky K. Yu. The solution of the problem of nonlinear waves on the surface of the low-viscosity fluid. Vestn. S.-Peterb. un-ta, 2011, vol. 10, no. 2, pp. 9-16 (in Russian).
Механика
329
