CC BY
47
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
(N.S.), Ser. Math. Mech. Inform., 2012, voI. 12, iss. 1, pp. 16-22 (in Russian),
9. Bejaneu A. Kahler contact distributions, J. Geom. Phys., 2010, vol, 60, no, 12, pp, 1958-1967,
10, Bukusheva A, V., Galaev S, V, Connections over a distribution and geodesic sprays. Russian Math., 2013, vol. 57, no. 4, pp. 7-13. DOI: 10.3103/S1066369X13040026.
УДК 517.958
УСРЕДНЕНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ АКУСТИКИ
А. А. Герус1, С. А. Гриценко2
1 Герус Артур Андреевич, аспирант кафедры математики, Белгородский государственный национальный исследовательский университет, [email protected]
2 Гриценко Светлана Александровна, кандидат физико-математических наук, доцент кафедры теоретической и математической физики, Белгородский государственный национальный исследовательский университет, [email protected]
В работе исследуется математическая модель акустики в гетерогенной среде с двумя компонентами, разделенными общей границей. Одна из компонент является ограниченной жидкой областью, другая — упругим телом. Упругое тело пронизано системой пор, заполненных жидкостью. Дифференциальные уравнения модели, описывающие движение жидкости и совместное движение твердого скелета и жидкости в порах, базируются на классических законах механики сплошной среды и содержат быстро осциллирующие коэффициенты, зависящие от малого параметра, равного отношению среднего размера пор к размеру рассматриваемой области. Быстро осциллирующие коэффициенты делают невозможным применение модели для численных расчетов. В работе доказывается существование обобщенного решения начальнокраевой задачи. На основе метода двухмасштабной сходимости Г. Нгуетсенга выводятся усредненные уравнения (т. е. уравнения, не содержащие быстро осциллирующих коэффициентов) для различных случаев. Полученные приближенные модели могут быть полезны для численных расчетов.
Ключевые слова: композитные среды, периодическая структура, уравнения Стокса, уравнения Ламе, уравнения акустики, пороупругость, усреднение периодических структур, двухмасштабная сходимость.
DOI: 10.18500/1816-9791 -2015-15-3-264-272
1. ПОСТАНОВКА ЗАДАЧИ
Рассматриваемая ограниченная область Q е R3 представляет собой единичный куб: Q = (0,1) х х (0,1) х (0,1), в котором пороупругая среда занимает область П = (0,1) х (0,1) х (0,а), 0 < a < 1, а область n(f), занятая жидкостью, есть открытое дополнение области П:
q = п и n(f) и s(0), s(0) = ап п ш(/).
Движение жидкости в пороупругой области П описывается системой уравнений:
С + i-£ p+у w = 0,
Cf C J
d 2 w
(Qf Xs + (1 - xlQs) -Q2 = V- P + QsF,
P = XsЙАx, J) + (1 - Xs)aAB(x, w) - pI,
(1)
(2)
(3)
где xs(x) — характеристическая функция порового пространства nf е П, cs и Cf — скорость звука в твердой и жидкой части соответственно, р — плотность среды, F — заданный вектор распределенных массовых сил, l — средний размер пор, L — характерный размер рассматриваемой области, малый параметр е = l/L. Здесь и далее в работе используются обозначения
B : C = tr(BCT),
© Герус А. А., Гриценко С. А., 2015
А. А. Герус, С. А. Гриценко. Усреднение математической модели акустики
где B, C — тензоры второго ранга, D(x,u) = -(Vu + Vu*) — симметрическая часть Vu, I —
2
единичный тензор.
Движение жидкости в области П(/) при t > 0 описывается системой уравнений Стокса:
Дг p + V • w = 0, г/
е/^ = V^ P(/> + е,F,
p(/) = aflx’jwW) - р!•
На общей границе S(0) выполняются условия непрерывности для перемещений
lim w(x, t) = lim w(x, t)
x^x
xefi(s)
x^x
xeo
и для нормальных компонент моментов
lim P(/)(x,t) • n(x0)= lim P(x,t) • n(x0)•
x^x
xen(f)
x^x
xGO
Завершают задачу однородные граничные условия:
w(x, t) = 0, (x, t) £ ST = S x (0,T),
на границе S = dQ, и однородные начальные условия:
d w
w(x, 0) = — (x, 0) = 0, x £ Q.
(4)
(5)
(6)
(7)
(8)
(9)
(10)
Пусть
/ (|F(x, t)|2 +
JQt \
d F л ~at (x’t)
dxdt = F2 < ж
и выполнены предположения о периодичности порового пространства и твердого скелета, а также о существовании пределов при г ^ 0 коэффициентов ам, ал. ., описанные в работе [1].
При выполнении упомянутых предположений о периодичности
X£(x) = С(x)X (7) ,
где Z(x) есть характеристическая функция области П.
Обозначим
е(/) = (1 - С)q, + С (q,X + (1 - X)Qs) •
Определение. Пара функций {w(, p (} таких, что
о 1,1
w( £Wr (Qt),p( £ Lr(Qt),
называется обобщенным решением задачи (1)-(10), если эти функции удовлетворяют уравнению неразрывности:
( i /х£ i - х(\ \
p ( + V • w ( = 0 (11)
(1 - С )4г + С|£ + -Т2Г
//
почти всюду в Qt , и интегральному тождеству
Lе/) (Ж'Зф + F-) dxd. = / (сp + «■ -
: D(x, ф) dxdt
(12)
о 1,0 Qm
для всех функций ф таких, что ф £W2 (QT), £ L2(Пт) и ф(x,T) = 0 для x £ Q.
2
S
Математика
265
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
Теорема 1. Для всех г > 0 на произвольном интервале времени [0, T] существует единственное обобщенное решение (w£, p£} задачи (1)-(10) и для него справедлива оценка
max / ( |p£(x,t)|2 + o<t<Tjn 1 V '
+nm« / ipe(x,t)i2 +
0<t<T JQ(f) \
+ max / n<t<TJ n
d w£
dt
(x,t)
+ (1 - X)ал|D(x, w£)| dx+
dw£ , , 2\ , , fT f £ ( d wA
dt (x’4) 1 dx + / X a» Jn 2n(f) Чж’ *)
dp £ 2 d2 w£ . . 2 - X£)ал ( dw£\
"dT (x’f) + d t2 (x-f) + (1- 4dt)
+ max
n<t<T J Q(f)
dp £ , л 2 d 2 w% .
it (x't) + ^ (x’‘>
2\ c T r
| dx +
+ x£ а
J Пт
_ , dw£
D|x^df
+
n JQ(f) 2
X£ а»
D ( x,
d 2 w£
~8F
dxdt+
dx+
2
dx dt+
D x,
d 2 w£ dt2
dxdt < Cn F2,
(13)
где постоянная Cn не зависит от малого параметра г и коэффициентов ал, аЛп) , а». Доказательство. Априорная оценка теоремы следует из энергетических тождеств:
2 dt й -
d w£
1 d_ Г
+ 2 dt JQ(s)
{
-s
dt
д w£
1
+ (1 - x£)алD(x, w£) : D(x, w£) + — |p£|2 dx+
aP
dt
+ (1 — X£)ал0)D(x, w£) : D(x, w£) + ——|p£|2
c.Sn)
\
dx+
+ I X\ а»D | x,
dwX ( dw£\\ f dw£
D x, —— dx = - £F • —— dx,
dt
1 d_
2 dt J П
1 d Г
2 dt Jn(s)
d 2 w£
dt2
dt
d w£
+ (1 — X£)алD ( x,~dtr ) : D ( x,
JQ
dw£
dt
-s
V
d2
dt2
+ (1 — X£)<э;ЛП)D I x
dt2
dw£
^^~aT
d2w£
: D x,
dt
+
1 а p dp£ dt 2 dx+
1 dp£ 2
(<»)■ dt /
dx+
+ ^ X ( а»D ( x, ddw ) : D ( x,
dt2
7 f _£ dF d2w£ 7 dx = q — • 2 dx.
Q
dt dt2
Эти тождества получаются подстановкой в уравнение (2) явного выражения для тензора P из уравне-
ее
□
ния состояния (3), умножением (2) на dwt) и интегрированием по частям по области Q. Далее существование и единственность обобщенного решения доказывается методом Галеркина.
2. УСРЕДНЕННЫЕ МОДЕЛИ
При выводе усредненных уравнений используются результаты А. М. Мейрманова [2-4] и применяется метод двухмасштабной сходимости Г. Нгуетсенга (G. Nguetseng) [5]. Главная проблема в доказательстве усредненных уравнений состоит в условиях непрерывности на общей границе S(n) между областями ) и G. Эти условия следуют из предельного интегрального тождества:
— f Р (V • ф) dxdt = J J^Q(f )(x, y) ^ F — ddt^ (x,t, y)^ • ^(x,t) dydxdt (14)
Qt
о m
для любой гладкой функции ф GW2 (Qt) и интегрального тождества:
[((,.. ^ 1 1 — шЛ dp . dw\
Ц Д — Z} у + z(У + ~б — V ф • вГ J
dx dt = 0
(15)
2
2
2
2
2
П
2
£
в
£
266
Научный отдел
А. А. Герус, С. А. Гриценко. Усреднение математической модели акустики
для любой гладкой функции ф е W21,0(QT)• Здесь W(x, t, y) — двухмасштабный предел последовательности (w£}, Y — ячейка периодичности (см. [5]) и
Q(f)(x У) = (1 - C(x)) Qf + z(x) (Qf Х(У) + (1 - Х(У)) Qs)) •
Для всех случаев (14) и (15) влекут систему уравнений акустики (17) и (18) в области QT), условия непрерывности (23) и (24) на общей границе S(0), и уравнение неразрывности
m 1 — m \ dp „ dw ^2 + —— -£ + V — = 0
cf cs2 I dt dt
(16)
в области QT•
Все различия сконцентрированы в уравнении динамики в области QT ив представлении скорости д w
смеси —— • dt
Теорема 2. Пусть [w£, p£} — обобщенное решение задачи (1)-(10) и
р1 = А1 = то.
Тогда пределы v = “w (скорость жидкости) и p (давление) последовательностей | ““L j и {р £} удовлетворяют системе уравнений акустики:
dv
Qf dt + V P = Qf F,
1 dp
^2^7 + V^ v = 0
cf dt
в области Q(f) при t > 0, и системе уравнений акустики в области Q при t > 0:
dv
Q dt + V p = QF,
(m 1 — m \ dp „
7 + -nr) at + V v = 0'
где
m = X(v)dy.
Соотношения (17)-(20) замыкаются однородным граничным условием:
v(x,t) • n(x) = 0
на границе ST, однородными начальными условиями:
p(x, 0) = 0, v(x, 0) = 0, x е Q,
и условиями непрерывности:
lim v(x, t) • n(x0) = lim v(x, t) • n(x0),
x^x
x^x
xGO
lim p(x, t) = lim p(x, t)
x^x
xen(f)
x^x
xGft
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
на общей границе S^°. Здесь Q = mQf + (1 — m) qs, n(x) есть нормальный вектор к S в точке
0\ ______с(0) „ „0
x е S, n(x0) — нормальный вектор к S(0) в точке x0 е S(0).
Математика
267
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
Доказательство. Для этого случая W(x, t, y) = w(x, t) и интегральное тождество (14) влечет
уравнение динамики: д 2 w
Q ж = V p + QF (25)
в области . □
Теорема 3. Пусть (w£, p£} — обобщенное решение задачи (1)—(10) и
0 ^ pi, Ai < то.
дw f Qw£ 1
Тогда пределы v = (скорость жидкости) и p (давление) последовательностей I г и {р£}
удовлетворяют в области qT) системе уравнений акустики (17), (18) и системе уравнений акустики в области , состоящей из уравнения баланса моментов в форме
v(x,t)= / B(a)(pi, Ai; t - T) •Vp (x,T)dr + f(x,t) J 0
(26)
и уравнения неразрывности (20).
Дифференциальные уравнения замыкаются граничным и начальным условиями (21), (22) и условиями непрерывности (23), (24).
Матрица B(a) (p1, A1; t) и функция f (x, t) задаются формулами (33), (34).
Доказательство. Уравнение неразрывности в этом случае имеет вид
1 1 \ др dw
32 + 32 тог + V • ~ттг = 0 x E q t > °-cf cs2 I dt dt
др
Используя вложение Vp E L2 (QT x Y), т. е. V p E L2 (QT), V ^ — ) E L2 (QT), выводим микро-
скопическое уравнение баланса моментов:
d2W / / д W \ \
Q(У)_дtГ' = Vy • (pix(y)^[y, -Qf) + Ai (1 - x(y))D(y, W) - Пij -
-Vp + g(y)F, y E Y, t> 0,
(27)
где
Q(y) = Qf X(y) + Qs (1 - X(y)), и микроскопическое уравнение неразрывности:
Vy • W = 0, y E Y.
Эти уравнения замыкаются однородными начальными условиями:
д W
W(x, y, 0) = (x, y, 0) = 0, y E Y.
Мы рассматриваем периодическое решение задачи как сумму
з , t
з , t
W(x, t, y) = £ [ W(i) (y, t - r) дХ (x, r) dr + [ W^ (y, t - r)F (x, r) dr,
i=iJo дх i=iJo
n(x, t, y) = £ / n(i) (y, t - r) дХ (x,r) dr + ^ f П<? (y, t - r)F (x, r) dr,
i=iJo дх i=iJo
где F(x,t) = (Fi(x,t),F2(x,t),F3(x,t)).
(28)
268
Научный отдел
А. А. Герус, С. А. Гриценко. Усреднение математической модели акустики
В свою очередь, пары {W(i), n(i)}, и (W^, П^} для i = 1, 2, 3 есть решения периодических начально-краевых задач в области Y для t > 0:
d2W(i) / /dW(i)\ ,л х.ч \ ,л
в{у)~дё^ = V' Iх(у)уЛ —^ ) + Ai (1 - x(y)) V,w« - n(i) 0 , V, ■ w(*> = о, (29)
V dt
.л dW(i)
W(i) (y, 0) = 0, g(y) dt (У, 0) = -ei, y G Y,
(30)
, s d2w(? (
Q(y) dt2F = V ■ I Mi X(y)v,
d W
(i)
+ Ai (1 - x(y)) V, Wj? - П« I | , V, ■ WF' = 0, (31)
(i)
dt
W^) (y, 0) = 0,
d W
(i)
dt
■(y, 0) = ei, y G Y
соответственно. Таким образом,
^ = U ^ (y,t - T) dp (x,T )dT + £/^ (y,t - T )Fi (x,T) dT
3 rt ялхт(д
(32)
d w
at
Л [4 dW(i) \ dp Л [4 dW(t\ , _ л,
Ы(—)Y(t-T)dX(x’T)dT + Ш dt /уiy,t-T)Fi(x,T)dT =
= / B(a)(pi ,Ai; t - t) -V p (x, t )dT + f (x,t),
Jo
где
(a) (Mi ,Ai; t) = ^dW2) (t) ® ei,
i = 1
dt /у
f (x,t) = E
3 rt / dw(i) \
,—Q^) (y,t - T )F (x,T) dT-
i=i у
Полученные соотношения дают следующее представление для скорости смеси:
д Z1 t
~дг (x,t) = J B(a) (Mi ,Ai; t - t) ■ V p (x,T)dT + f (x,t). Обозначение: матрица a ® b определяется как
(a ® b) ■ c = a(b ■ c)
для любых векторов a, b, и c
Теорема 4. Пусть (w£, p£} — обобщенное решение задачи (1)-(10)
Mi = то, 0 ^ Ai < to,
и w£ = (w£) — продолжение из в
д w
(33)
(34)
□ (35)
Тогда пределы v = — (скорость жидкости) и p (давление) последовательностей . dt
и (p£ } где
д w£
dwf dw(s) ^ ч
v = (1 - Z)v + Zm-^f + G-dp = (1 - Z)v + Zm Vf + z v(s),
(36)
и
F
И
Математика
269
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
и w(s) и wf — пределы последовательностей {(1 — х£)w£} и (w£}, удовлетворяют в области пТ ) системе уравнений акустики (17), (18) и системе уравнений акустики в области ПТ, состоящей из уравнения баланса моментов:
mpf vf + q,s v(s) +
Г
о
(—QF + V p) (x, r) dr = 0
(37)
для жидкой компоненты, уравнения баланса моментов
v(s) — (1 — m)vf
■t
B(s) (ж, Ai; t — r) •
о
V p + Qs
df — F
dr
(x, r)dr
(38)
для твердой компоненты, и уравнения неразрывности
m + !1ir Ш + V(m vf + Vs") = 0
(39)
Задача замыкается граничным и начальными условиями (21), (22) и условиями непрерывности (23), (24).
В уравнениях (37), (38) Q = mgf + (1 — m) gs и матрица B(s) (ж, Ai; t) задается формулой (48).
В теореме используется обозначение
wf = Em (w£),
где EQe : Wl(Hf) ^ Wl(П) — оператор продолжения из Of в П, такой, что w£ = w£ в Of х (0,T),
|D(x, w£)|2 dx < C0 |D(x, w£)|2 dx.
|w£ |2 dx < C0 |w£|2 dx,
'Q JQ JQf
Корректность такого продолжения обоснована в работе С. Сопеа [6]. Доказательство. Для этого случая скорость смеси задана формулой
дW, , дwf , ЛЛ dW,
(x,t,y) = (x,t) + (1 — X(y)) (x,t,y),
dt
д w dwf дw(s) t \ f
~Qt~ = m-Q-f + —дj^, w(s) (x,t) = J (1 — Х(У)) W(x,t. У) dV-
дЬ
дЬ
Интегральное тождество (14) влечет уравнение динамики (37) для жидкой компоненты.
Получим представление (38).
Если W(s) = (1 — x(y)) W, то пара {W(s), n(s)} удовлетворяет микроскопическому уравнению динамики для твердой компоненты:
Qs = A1 Ay W(s) — Vy n(s) — Vp,
дЬ2 2
микроскопическому уравнению неразрывности:
V • W(s) = 0
в области Ys, и начальным условиям
(40)
V \ QW(s)
W(s) (x, 0, y) = -— (x, 0, y) = 0, y e Ys.
(41)
В силу теоремы Нгуетсенга W(s),Q2W(s)/дt2, VyW(s) e L2(QT х Ys). Эти условия вместе с
формулой (30) обеспечивают граничное условие:
W(s) (x, t, y) = wf (x, t), (y, t) e S(0) х (0, T).
(42)
И
270
Научный отдел
А. А. Герус, С. А. Гриценко. Усреднение математической модели акустики
Решения {W(s), n(s)} периодических начально-краевых задач (40)-(42) имеют вид
3 п t \ f dp f , d2 wn .
- T)' dX (x,T) + Q^otr (x,T) * dT’
W(s) = wf (x, t) + ^J^ w(s) (y, t - T)(^dX (x? t) +
n(s) (x, t, y) = ^ ^ n(s) (y,t - t) ^dX (x,t) + Qs ^dW± (x? t^ dr,
где {W( ), П( )}, i = 1,2, 3, в свою очередь, являются решениями периодических начально-краевых задач:
Qs^dX = A2 △ W((s) - Vyn(s), (y,t) e Ys x (0,T), (43)
Vy ■ W(s) (y,t) = 0, (y,t) e Y, X (0,T), (44)
( ) dW(s) W(( ) (y, 0)=0, Qs dti (y, 0) = ei, y e Ys, (45)
W(s) (y,t) = 0, (y,t) e S(0) X (0,T). (46)
Однозначная разрешимость задач (43)-(46) следует из энергетического тождества:
, d w(s) Qs 1 dt
(y,t)|2 + 2|VWX(y,t)|2 j dy
(1 - m) Qs
Задача (43)-(46) для соленоидальных функций W(i), равных нулю на S(0) и при t = 0, понимается как интегральное тождество:
рТ
0 JYs
Q дW(s) дv Qs dt ' ~8t
AiVW(s) : Vv
dy dt =
ei ■ v(y, 0) dy
для любой соленоидальной 1-периодической гладкой функции v, равной нулю на S(0) и при t = T. По определению
dw(s), л [ dW(s),
"dT (x’t)= L~dT (x’f’y) dy =
= (1 - m)
d Wf
~дГ
-t f 3
£
i = 1 Js
d wf
JYS
Y ~dt~ (y,t - T)d^ ® e^ ■ (V p (x,T) + Qs ddT2f (x,TX dr =
= (1 - m)~d~ - J B(s) (TO’ Ai; t - t) ■ ^Vp (x, t) + QsddT2f (x,r)^ dr,
(47)
где
-».Ai,) = Ё(l, ^
■(y,t)dy <g> ei.
□
(48)
Работа выполнена при финансовой поддержке Российского научного фонда (проект № 14-17-00556 «Математическое моделирование флюидопотоков в нефтяных резервуарах с учетом разномасштабных свойств пласта-коллектора»).
Математика
271
Изв. Сарат. ун-та. Нов. сер. Сер. Математика. Механика. Информатика. 2015. Т. 15, вып. 3
Библиографический список
1. Геру с А. А., Гриценко С. А. Модель акустики в конфигурации упругое тело - пороупругая среда// Научные ведомости БелГУ. Сер. Математика. Физика. 2014. № 25 (196), вып. 37. С. 68-75.
2. Мейрманов А. М. Метод двухмасштабной сходимости Нгуетсенга в задачах фильтрации и сейсмоакустики в упругих пористых средах // Сиб. ма-тем. журн. 2007. Т. 48, № 3. С. 645-667.
3. Мейрманов А. М. Уравнения акустики в упругих пористых средах // Сиб. журн. индустр. матем. 2010. Т. XIII, № 2. С. 98-110.
4. Мейрманов А. М. Вывод уравнений неизотермической акустики в упругих пористых средах // Сиб. матем. журн. 2010. Т. 51, № 1. С. 156-174.
5. Lukkassen D, Nguetseng G, Wall P. Two-scale convergence // Intern. J. Pure and Appl. Math. 2002. Vol. 2, № 1. P. 35-86.
6. Conca C. On the application of the homogenization theory to a class of problems arising in fluid mechanics // Math. Pures et Appl. 1985. Vol. 64. P. 31-75.
Homogenization of the Acoustics Mathematical Model
A. A. Gerus, S. A. Gritsenko
1 Gerus Artur Andreevich, Belgorod State National Research University, 85, Pobedy st., 308015, Belgorod, Russia, [email protected]
2Gritsenko Svetlana Aleksandrovna, Belgorod State National Research University, 85, Pobedy st., 308015, Belgorod, Russia, [email protected]
We consider a mathematical model of acoustics in heterogeneous medium with two different components with the common boundary. One of these is a bounded liquid domain and the other is a poroelastic medium. Poroelastic medium is perforated by pores. A pore space is filled with a viscous liquid. The motion of the liquid and the joint motion of the poroelastic media with porous space are governed by the differential equations based on the continuum mechanics laws. These equations contain rapidly oscillating terms, depending on the small parameter. The small parameter is the ratio of the average pores size to the size of domain under consideration. Rapidly oscillating terms prevent from the numerical simulations. The unique existence of the generalized solution of the boundary-value problem is proved. Homogenized equations (i.e. free from rapidly oscillating terms) are based upon the Nguetseng method of the two-scale convergence. We derived approximate models useful to the numerical calculations.
Key words: composite medium, periodic structure, Stokes equations, Lame’s equations, acoustics equations, poroelastic, homogenization of periodic structures, two-scale convergence..
This work was supported by the Russian Scientific Fond (projects no. 14-17-00556).
References
1. Gerus A. A., Gritsenko S. A. Acoustics Model in the Configuration Elastic Body — elastic porous medium. Belgorod State University Scientific Bulletin. Mathematics and Physics, 2014, no. 25 (196), iss. 37, pp. 68-75.
2. Meirmanov A. Nguetseng’s two-scale convergence method for filtration and seismic acoustic problems in elastic porous media. Siberian Mathematical Journal, 2007, vol. 48, pp. 519-538. DOI: 10.1007/s11202-007-0054-9.
3. Meirmanov A. Acoustics Equations in Elastic Porous Media. Sib. Zh. Ind. Mat., 2010, vol. XIII, iss. 2, pp. 98-110 (in Russian).
4. Meirmanov A. Derivation of the equations of nonisothermal acoustics in elastic porous media. Sib. Math. J., 2010, vol. 51, iss. 1, pp. 128-143. DOI: 10.1007/s11202-010-0014-7.
5. Lukkassen D., Nguetseng G., Wall P. Two-scale convergence. Intern. J. Pure and Appl. Math., 2002, vol. 2, no. 1, pp. 35-86.
6. Conca C. On the application of the homogenization theory to a class of problems arising in fluid mechanics. Math. Pures et Appl., 1985, vol. 64, pp. 31-75.
272
Научный отдел
