Математические модели диффузии в пороупругих средах Текст научной статьи по специальности «Математика»

УДК 517.958:531.72, 517.958:539.3(4)

МАТЕМАТИЧЕСКИЕ МОДЕЛИ ДИФФУЗИИ В ПОРОУПРУГИХ СРЕДАХ и)

А.М. Мейрманов, Р.Н. Зимин, О.В. Гальцева, О.А. Гальцев

Белгородский государственный университет,

ул.Студенческая, 14, Белгород, 308007, Россия, e-mail: reshat85@mail.ru

Аннотация. Настоящая работа посвящена получению новых корректных математических моделей, моделирующих процесс диффузии в пороупругой среде посредством усреднения и метода двухмасштабной сходимости Нгуетсенга.

Ключевые слова: усреднение, метод двухмасштабной сходимости.

1. Постановка задачи

Пусть П Е Я3 ограниченная связная область с липшицевой границей Б, полученная периодическим повторением элементарной ячейки еУ, где е > 0 малый параметр,

У = У/ и У и 7 и дУ , У = (0,1) х (0,1) х (0,1), еУ = (0,е) х (0,е) х (0,е),

и 7 = дУ/ П дУ8 - липшицева граница между множествами У/ и У,. Область У/ будем считать симметричной относительно поворотов на п/2 (рис 1.).

Через П уг обозначим периодическое повторение элементарной ячейки еУ/, а через - периодическое повторение еУ 3. Тогда

П = П/ и Пе8 и г,

где Г = дП / П дП £3 периодическое повторение границы е7.

ь, ь ь

Рис. 1. Геометрия элементарной ячейки.

11 Работа выполнена при поддержке Министерства образования Российской Федерации, соглашение M4.A18.21.0357

Мы рассмотрим две модели диффузии в пороупругой среде, состоящие из системы уравнений Стокса, которые описывают движение несжимаемой вязкой жидкости в по-ровом пространстве и системы уравнений Ламе, описывающих колебания несжимаемого твердого скелета. Рассматриваемая система дополняется конвективным уравнением диффузии для примеси в жидкости. Задача замыкается начальными и краевыми условиями. При этом считается, что плотность жидкости зависит от концентрации примеси.

Модель РЕЬ

V- (хЗД^я,^^ +(1-хе)АоВ(х^£)-р£^)+р(с£)Г = 0, жеП, ¿>0;(1)

V ■ w£ = 0 , х е П, г > 0; (2)

дс£ дw

+ • Ус£ = А)Лс£, х е щ, г>о. (з)

дг дг

Система уравнений дополняется следующими условиями:

w £(х,г) = 0 при х е 5 = дп, (4)

х £w £(х, 0) = 0 при х е П , (5)

/ р£ йх = 0 , (6)

дс £

— = о при X еГ£ и Б, г > о, (7)

дп

с£(х, 0) = с0(х), при х е п ^ , (8)

где w£(х,г) = (/ш£(х,г),'ш£(х,г),'ш£(х,г)) - перемещение сплошной среды, р£(х,г) -давление в сплошной среде, с£(х,г) - концентрация примеси, 0(ж, V) - симметрическая часть градиента вектора V (тензор напряжений), I - единичная матрица, х £(х) -характеристическая функция порового пространства,

р(с£) = х£(х)8с£(х,г),

^о - безразмерная вязкость жидкости, Ао - безразмерная постоянная Ламэ, 8 - положительная постоянная, п - единичный вектор внешней нормали к Г £, О0 - коэффициент диффузии.

Уравнение (3) понимается в смысле следующего интегрального равенства

[ I (V ^ - (Ус£) • ^ - Д^с£ • щ) сЬсЫ =- [ со(х)ф(х,0 )с1х ./о * п с V дг дг ) * п е

для произвольной гладкой функции ф(х,г), равной нулю при г = Т.

Модель PE2

= жеО, t> 0; (9) V - w£ = 0 , x G Q , t > 0; (10)

dc£ , dw

V- (D0Vc£- xeïl£f,t> 0, (11)

ôt V dt

где функция <^(c£) G C2(-o, oo) такая, что

f -1/2 , c£ < -1/2; ^(c£) = < c£ , 0 < c£ < 1; [ 3/2 , c£ > 3/2 .

Задача замыкается следующими условиями:

w £(x, t) = 0 при x G S = dQ , t> 0; (12)

f p£dx = 0 , t > 0; (13)

Jn

dc £ dw £

L>0—- ^(c£) • n = 0 при ж G Г£, t>0, (14)

dn dt

X £w £(x, 0) = 0 при x G Q , (15)

c£(x, 0) = co(x) при x G Q f . (16)

Уравнение (11) понимается в смысле теории распределения, то есть в смысле следующего интегрального равенства

T r / dw£ \ f

с.£ + </?(с£) • VV' - D0Vc.£ ■ Щ\dxdt = - со(х)ф(х, 0)dx ,

'0 и О. £ V д£ д£ / и О £

для произвольной гладкой функции ^(ж,£), равной нулю при £ = Т.

2. Основной результат

Определение 1. Тройка {ш £(х,£), р£(х,£), с£(ж,£)} называется обобщенным решением модели РЕ! в области Пт = П х (0, Т), если

1) р£ Е Ь2(Пт), ш£ Е Ж21Д(П / х (0,Т)) П ^21,0(П 8 х (0,Т)), с£ Е Ьте(П/ х (0,Т)) П ^21>0(П/ х (0,Т));

2) почти всюду в области Пт выполнено уравнение (2) и условие (6);

3) функции ш £, р£ и с £ удовлетворяют интегральным равенствам

У (^хЧФ : + (! - ги£) : 0(;г, ф) - р£У • </>) (М

= р(с£)^ ■ ф вжИ (1)

£

для произвольной гладкой вектор-функции ф(х,г), равной нулю на границе Б и при г = Т, а также

[ [ (с/ ^Г - ' Чф)(1х(И= - [ Со(х)ф(х,0)(1х (2)

]о Уп ^ V дг дг ) Jn е

для произвольной гладкой функции ф(х,г), равной нулю при г = Т.

Здесь используется обозначение: А : В = Ьт(ЛБт), где А и В - квадратные матрицы.

Определение 2. Тройка ^ £(х,г),р £(х,г), с£(х,г)} называется обобщенным решением модели РЕ2 в области Пт = П х (0, Т), если

1) р£ е ь2(Пт), w£ е ж21,1(п£ х (0,Т)) П ^21,о(п8 х (0,Т)), с£ е Ь^(П £ х (0, Т)) П ^(П £ х (0, Т));

2) почти всюду в области Пт выполнено уравнение (10) и условие (13);

3) функции w £, р £ и с £ удовлетворяют интегральным уравнениям

I пт

X £^оЩ х

дwí

: Р(х, ф) + (1 - X £)АоО(х, w £) : Р(х, ф) - р£''V ■ ф) йхйг

р(<р(с£))* ■ фйхйг (3)

пт

для произвольной гладкой вектор-функции ф(х,г), равной нулю на границе Б и при г=Ти

гТ г / дф дw£ \ Г

с.£ ^ + -Чф^ЛхЛ^ - I со(х)ф(х,0)(1х (4)

/о Jпе

дг

дг

'п е

для произвольной гладкой функции ф(х,г), равной нулю при г = Т.

Разрешимость модели РЕ1. Теорема 1. Пусть

0 ^ со(х) ^ 1,

пт

(V |2 +

д¥

дг

)йхйг ^ г2, IV*1 (х,г)| ^ г.

Тогда модель РЕ1 имеет хотя бы одно обобщенное решение и для него справедливы оценки:

тах Ао (1 - х £) |0(х, w £)|2 йх+

о<ь<т

^х

пт

Ш х

дw{

+ |р£|2 \йхйг ^ МГ2 , (5)

2

п

2

0 ^ с£(х,£) ^ 1

т

|У с£^ М^2 ,

/0 JО £

тах / х

0<*<т] О

Ш х

д'

¿х + / (1 - х£)

иот

Ш х

д'

¿х^ ^ М^2

|Уи £(х,^) - Vv £(ж,^2)|2^х ^ С|£1 - ¿211/2

где функция V £ = ЕО £ (дш £/д£) описана ниже.

Доказательство разрешимости приведено в ( [7]) Разрешимость модели РЕ2. Теорема 1. Пусть

0 ^ с0(х) ^ 1,

/от

1 д£ 1 От 1 1 л

(6)

(7)

(8) (9)

Тогда модель РЕ2 имеет хотя бы одно обобщенное решение и для него справедливы оценки:

тах А0 / (1 - х£)|0(х, ш£)|2^х+

0<ь<т

^0Х

От

Щ х,

дш

+ |р£|2 М; ^ М^2 , (10)

тах / |с£| ¿х + Д

0<ь<т

т

/0 ./о ^

|У с£|2^х^ ^ М^2

(11)

Легко заметить, что вектор-функция ш£ имеет различную гладкость в областях П/ и П8. Поэтому мы воспользуемся известными теоремами о продолжении (см. [1]- [3]): Существует линейный оператор

Ео.в : ВДП;) ^ ^(П), v£ = Ео

дш

' \~дГ

такой, что

дш

~дГ

(х,£) = V £(х,£), х Е П/ , £ Е (0,Т)

(12)

IV £(х,£)|2 ¿х ^ С0

/О/

дш

д£

(х,£)

¿х,

|Р(х, V £(х,£)) |2^х ^ С0

дш£

(ж,*))

¿х, £ Е (0,Т) , (13)

2

2

О

О

2

О

£

/

2

О

2

О

где Со не зависит от £ и г е (0,Т).

Для усреднения необходимо рассматривать функции в фиксированной области, но концентрация примеси с£ определена только в области П£, которая зависит от малого параметра £. Поэтому концентрацию примесей с£ также продолжаем из П£ в П (см. [1]):

С£ = Е п= (с £), (14)

такой, что

££

с£ = С£, х е П£, г е (0,Т), (15)

и

|с£(х,г)|2^х ^ Со ^£(x,г)|2dx,

|2 Л™ ^ П I 1У7 „£(„. 4-М2,

Vc£(x,г)|2dx ^ Со c£(x,г)|2dx, г е (0,Т), (16)

где Со не зависит £ и г е (0, Т).

Под усреднением мы будем понимать предельный переход в интегральных тождествах (1), (2), (3) и (4) при £ ^ 0. Для этого воспользуемся понятием двухмасштабной сходимости ( [5]). Напомним, что последовательность {и£к} ограниченная в Ь2(П) сходится двухмасштабно к 1-периодической по переменной у функции и(х, у):

У и£к(х)(р ^ж,—^ в,х —> И и(х,у)(р(х,у)с1ус1х

п У

для любой 1-периодической по переменной у функции р(х, у). В частности, это верно для <р(х, у) = ^о(у) ■ к(х), где ^о е Ь2(У) и к е Ьте(П).

Усреднение модели РЕ1. Верна следующая

Теорема 3. Пусть тройка ^£, р£, с£} является обобщенным решением модели РЕ1. Тогда

1) последовательности ^£}, {Vw£}, {V£}, {Vv£}, {р£}, (¡£ и Vсс сходятся слабо в Ь2(Пт) к функциям w, Vw, V = Еп^ (дw/дг), Vv = V(Eп£ (дw/дг)), р, с и Vc соответственно;

2) Предельные функции являются решением усредненной системы уравнений в области пт, состоящей из уравнения неразрывности

V■ w = 0 , (17)

усредненных уравнений баланса

п

V ■ Р + шбсЕ = 0

(18)

где

дш Г*

Р = -р! + ^ : Щх, —) + : Щх, го) + / т3{1 - т) : Щх, ю{х, т))йт ,

и усредненного конвективного уравнения диффузии

дс

т— + Ви • Ус = V • (Д)ВУс) , (19)

д£

дополненных граничными

ш(х,£) = 0 , х Е Б, £ Е (0,Т), дс

— = 0, же ге (о,т),

дп

и начальными условиями

ш(х, 0) = 0 , х Е П , с(х, 0) = тс0(х) , х Е П;

3) тензор 4-го ранга N1 - симметричный и положительно определенный. Матрица В - симметричная и положительно определенная.

Усреднение модели РЕ2.

Теорема 4. Пусть тройка {ш £, р£, с£} является обобщенным решением модели РЕ2. Тогда

1) последовательности {ш£}, {V}, }, {Vv£}, {р£}, с£ и Ус£ сходятся слабо в Ь2(Пт) к функциям ш, Vw, V = ЕОв (дш/д£), Vv = V(EО£(дш/д£)), р, с и Vc соответственно;

2) Предельные функции являются решением усредненной системы уравнений в области Пт, состоящей из уравнения неразрывности

V ■ ш = 0 , (20)

усредненных уравнений баланса

V- Р + = 0 , (21)

где

/ дш \ Г*

Р = -р! + : О ( х, + : Щх, го) + J - г) : го (ж, г))йг ,

и усредненного конвективного уравнения диффузии

дс

т,— = V • (Д)ВУс) - V • (р(с)В«) (22)

дополненных граничными

ю(х,г) = 0, х е Б, г е (0,т)

дс

— = 0, же ¿е(о,т), (23)

дп

и начальными условиями

ю(х, 0) = 0 , х е П

с(х, 0) = тсо(х) , х е П; (24)

3) тензор 4-го ранга N - симметричный и положительно определенный. Матрица В - симметричная и положительно определенная.

Так же верна следующая

Лемма 1. Матрица В путем ортогонального преобразования координат приводима к виду В = А1.

Тогда, с учетом данной леммы усредненное конвективное уравнение диффузии (37) примет вид

дс

т— + АУ (<р(с)) ■ V = АУ • (Д)Ус), (25)

дополненное граничными и начальными условиями (23) и (24).

Для указанной задачи, согласно ( [4]), справедлив принцип максимума. Учитывая все вышесказанное система (37)-(33) примет вид

V ■ Р + тбсГ = 0 , (26)

V ■ ю = 0 , (27)

дс

т,— + АУс • V = АУ • (Д)Ус). (28)

3. Доказательство Теоремы 3

Сначала перепишем интегральные уравнения (1) и (2) для продолженных функций дф дм1

/7

/о ио

( дф дм \

Xе -с£^~+Ф ■ V с£ + И0 V с£ • V ф)(1х(И =

\ дг дг /

[ х£ с о (х) ф(х, 0) дх, (1)

/7

/о Зо

'п

(хе№ »(х, V е) + (1 - х£)А0 »(х, ю е) - р£1) : »(х, ф)дхдг

ГТ р

/ / х£$сеГ ■ ф дхдг , (2) 'о ио

где v£ = Еп, (дш£/д£) и с£ = Еп, (с£).

Также запишем уравнение неразрывности (2) в виде интегрального уравнения

т

/ / ш£ -V ^ж^ = 0 , (3)

для произвольной гладкой функции Согласно оценкам (5)-(9), имеем

дш

С, Р, V С, V, V V Е ^((0,Т); ¿2(П х У)) . Поэтому, с точностью до подпоследовательностей,

с£(х, ¿) ^ с(х, £) слабо в Ь2((0,Т); И2(П)) ,

с£(х,£) — с(х,£) двухмасштабно в Ь2 ((0,Т); Ь2(П)) , V с£ — V с + V С двухмасштабно в Ь2 ((0, Т); Ь2(П)),

р£(х,£) ^р(х,£) слабо в Ь2((0,Т); Ь2(П)) , р£(х,£) — Р(х, у,£) двухмасштабно в Ь2 ((0,Т); Ь2(П)),

ш£(х,£) ^ ш(х,£) слабо в Ь2 ((0,Т); И (П)),

Р(ж, ш£) — Р(ж, ш) + Р(у, двухмасштабно вЬ2((0,Т); Ь2(П)),

О

v£(х,£) ^ v(x,í) слабо в Ь2 ((0,Т); И (П)), Р(ж, v£) — Р(ж, V) + Р(у, V) двухмасштабно в Ь2((0,Т); Ь2(П)),

ч д^.

У(х,у,г) = —(ж,у,*), v£(х,£) — V(х, £) сильно в Ь2 ((0,Т); Ь2(П))

при е — 0.

Лемма 2. Предельные функции удовлетворяют в области П при £ > 0 системе макроскопических уравнений

V ■ ш = 0 , (4)

V ■ Р + т£с ^ = 0 , (5)

Р = ц0(тВ(х,ь) + ) + Ао((1 -т)В(х,Уз) + (Р(у,1¥))Гз)

дс ( ) ( )

•т — + V • (т V с + (Уу С)У/) = Р0 V • (т V с + (Уу С)У/). (6)

Для доказательства леммы перейдем к пределу при е — 0 в интегральном тождестве (1) с пробной функцией ф = ф(х,г), в (2) с пробной функцией ф = ф(х,г), и в (3) с пробной функцией £ = £(х,г).

Лемма 3. Предельные функции удовлетворяют, в области У для всех х е П и г > 0, системе микроскопических уравнений

Vy ■ Ш = 0 , (7)

Vy ■ Р = 0 , (8)

Р = »оХ(У) (у, ^^ ) + А0 (1 - Х(У)) + В(у, Ш)) - РI,

Vy ■ (х(y)(V с + Vy С)) =0 . (9)

Для доказательства леммы перейдем к пределу при е — 0 в интегральном равенстве (1) с пробной функцией ф = еф0(х,г)ф1(х/е), в (2) - с пробной функцией ф = еф0(х,г)ф1(х/е), и в (4) - с пробной функцией £ = е£0(х,г)£1(х/е).

Лемма 4. 1-периодическиерешения {Ш(х, у, г), Р(х, у, г)} периодической начально-краевой задачи (7), (8), дополненной начальным условием

х(у)Ш(х, у, 0) = 0 , (10)

имеет вид

3 г г

Ш = ^ / Ш(гз)(у,г - т)Zij(х,т)дт,

1,3 = 1

3 ,1

Р = х(у)£ Ро(гз)(у^(х,г) + £ / Р(гз)(у,г - т)%(х,т)дт

г,з=1

г,з=1'

где

(х,Ь) = цоЩх, - )ф(х,го) =

г,з=1

и {Ш(гз)(у,г), Р(гз)(у,г)}, {Ш(г3)(у), Р(гз)(у)}, г,з = 1, 2, 3 являются решениями периодических задач

Vy ■ х^о» У,

д\¥{гз) дЬ

+ Ао (1 - х)»(у, Ш(гз)) - Р(гз) I) = 0,

Vy ■ Ш(гз) = 0, х(у)Ш(гз)(у, 0) = Ш(гз)(у); J

в области У. Здесь

V ■ (х (^0 0(У, + I7 - Р. I)) = 0, V ■ ^= 0, ^ х(у)^0гЛ(у)^ = 0,

Г-7 = ^(е» О е^ + е,- О е»),

(12)

где (е1, е2, е3) вектора декартова базиса. Простые вычисления показывают, что

Р(у, ^) = А0(у) : Р(ж, ш)+ [ А1(у,£ - т): Р(ж, ш(х,т))^т.

где

и

Ао(у) = »(У, ^0

¿,.7 = 1

(13)

(14)

= Е (/'«В у-—(!!-<) - А0О(у, И^(У^))) о (15)

¿,7=1

Из равенств (13) - (15) следует, что ( / дш \

-^-(у, I — т) : т))с?т .

* дА

Тогда

N = т ^^

¿7) 0 I7

+ ^о (Ас) у.

¿,7=1

N2 = Ао (1 - т) ^ I7 0I7 + Ао (Ао)у3 + ^о (^(у, 0))У/,

¿,7=1

д£

(16)

Лемма 5. 1-периодическое решение С(х, у,£) периодической начально-краевой задачи (9) имеет вид

¿=1

где С^(у), г = 1, 2, 3 является решением краевой задачи

V ■ (х(у)(е + V C(¿))) =0 (17)

0

88 НАУЧНЫЕ ВЕДОМОСТИ ЩД Серия: Математика. Физика. 2012. №17(136). Вып. 28 в области У.

Задача (17) является уравнением Лапласа в области У^, дополненное краевым условием

(ег + VyС(г)) ■ п = 0 на границе 7, где п вектор нормали к 7. Тогда

т V с + V С = В(с^ с = |т1 + ^ С(г)(у))У/ 0 ег V с,

и

( 3 )

(с) = т1 + ^уС(г) (у))у} 0 е%) = т1 + Во . (18)

г=1

4. Доказательства утверждений

Доказательство Теоремы 4 аналогично доказательству Теоремы 3.

Доказательство Леммы 1. Матрица В - симметричная и положительно определенная. Действительно, умножим уравнение (17) на С(з)(у) и проинтегрируем по области Уf

(VС(г) ■ VC(з) + ег ■ VC(j))dy = 0 . (1)

/у

f

Пусть п е Я3. Рассмотрим следующее выражение

3

П = Дот| ■ п + А,^/ VC(%)£l ■ пду = Ао / х(£ ■ П + ■ п)ду, (2)

г=1 и ^ ^

где

33

^С(")£1, = £С(з)пз.

П =

г=1 3=1

Умножим (1) на Б0пз£г и просуммируем по индексам г и ]. В итоге, получим

I х^щ ^фг, + I )ду = 0. (3)

Прибавим теперь к выражению (2) равенство (3). Имеем

П = Ао I х(й ■ П + VЩ ■ П + VVt ■ VФп + 1 ■ VФп)дУ

или

П = Ао уу х(£ + Vфí:) ■ (п + Vфп)ду .

Очевидно, что (BI) • п = (®п) • I, что, в свою очередь, и означает симметричность матрицы B, и

(BI) • | = Do J x(I + V^)2dy ^ 0 .

Более того

I (I) = Do J x(I + )2dy ^ ao > 0

для III = 1.

В самом деле, пусть последовательность {|n} такая, что I(|) М 0 при n М- то. Можно считать, что М |0 при n М то. Тогда, в силу непрерывности I(|), следует, что

I (I0) = 0,

то есть

I0 + = 0

или

^0 = -|° • у + const. (4)

Функция ^0 = C+ C(2)^0 + C(3)является 1-периодической по y, так как функции C(i)(y) - 1-периодические по у.

Так как поровое пространство является связным, то (4) имеет смысл только для I = 0, но II | = 1. Получили противоречие.

Пусть Т2 - поворот на угол п/2 вокруг оси y3, направленной по вектору e3, такой что Т2• e2 = e1. Здесь {e1, e2, e3} декартов базис. Если z = T2• у, у G Y, то T2 : Yf М Yf и х(у) = x(z), так как поровое пространство симметрично относительно поворотов на п/2 относительно любой оси координат.

Функция (7(2)(z) = C(2)(у) является решением следующей периодической начально-краевой задачи

Vz • (x(z)(T2 • e2 + T2 • T2 • VzC(2))) = Vz • (x(z)(ei + VzC(2))) = 0 (5)

в единичном кубе Z = Y.

Задача (5) совпадает с задачей (17) и, в силу единственности решения задачи (17), имеем

C(2) (у) = <5(2) (z) = C(1) (z) = C(1) (T2 • у). (6)

Аналогично получаем

C(3) (у) = С(3) (z) = C(1) (z) = C(1) (Тз • у), (7)

где Т3 поворот вокруг оси y2, такой что Т3 • e3 = e1. Для нахождения матрицы B0 положим, что

a = (VyC(1)(у))у,,

и вычислим

а2 = (VyC(\y))Yf = T2(VZC(1\z))Yf = T2 ■ oi,

и

Оз = (VyC(3)(y))Yf = T£<V*C(1)(z))Yf = T3 ■ oi.

Пусть o1 = a1e1 + a2e2 + a3e3. Тогда

02 = a{T*2 ■ e1 + «2^ ■ e2 + азТ2 ■ ез = a,^ — a2el + азез ,

03 = a^ ■ e1 + a2T3 ■ e2 + азТ2 ■ eз = аleз + a2e2 — aзel. Здесь мы воспользовались очевидными свойствами преобразований Т2 и Тз:

Т ■ e1 = e2, Т ■ e2 = — e1, T2 ■ eз = eз,

и

Т ■ e1 = eз, Тз ■ e2 = e2, Тз ■ eзel. Учитывая введенные выше обозначения, матрица B0 имеет следующий вид

Bo = 01 ® e1 + 02 ® e2 + Оз ® eз =

a1e1 ® e1 + a1e2 ® e2 + aleз ® eз + a2e2 ® e1 — a2e1 ® e2+ aзeз ® e1 — aзel ® eз + aзeз ® e2 + a2e2 ® eз =

a1(e1 ® el + e2 ® e2 + eз ® eз),

где a2 = аз = 0 силу симметричности матрицы B0.

Литература

1. Acerbi E., Chiado Piat V., Dal Maso G. and Percivale D. An extension theorem from connected sets and homogenization in general periodic domains // Nonlinear Anal. - 1992. - 18. - P.481-496.

2. Nguetseng G. Asymptotic analysis for a stiff variational problem arising in mechanics // SIAM J. Math. Anal. - 1990. - 21.- P.1394-1414.

3. Conca C. On the application of the homogenization theory to a class of problems arising in fluid mechanics // Math. Pures et Appl. - 1985. - 64. - P.31-75.

4. Ладыженская О.А., Солонников В.А., Уральцева Н.Н. Линейные и квазилинейные уравнения параболического вида / М.: Мир, 1967. - 736 с.

5. Nguetseng G. A general convergence result for a functional related to the theory of homogenization // SIAM J. Math. Anal. - 1989. - 20. - P.608-623.

6. Мейрманов А.М., Зимин Р.Н., Гальцева О.В., Гальцев О.А. Корректная разрешимость задачи о нелинейной диффузии в несжимаемой пороупругой среде на микроскопическом уровне / Научные ведомости БелГУ. Математика. Физика. - 2012. - 5№124;26.

MATHEMATICAL MODELS OF DIFFUSION IN POROUS ELASTIC MEDIA A.M. Meirmanov R.N. Zimin, O.V. Galtseva, O.A. Galtsev

Belgorod State University, Studencheskaja St., 14, Belgorod, 308007, Russia, e-mail: reshat85@mail.ru

Abstract. The paper is devoted to obtaining of some new correct mathematical models of diffusion process in porous media. They are derived on the basis of the averaging and Nguetseng's method of two-scale convergence.

Key words: homogenization, two-scale convergence.