УДК 532.51
БОТ: 10.18698/2309-3684-2019-3-1938
Моделирование проницаемости неньютоновских жидкостей в трехмерных композитных структурах на основе метода асимптотической гомогенизации
© Ю.И. Димитриенко, Шугуан Ли МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
Исследована математическая модель многомасштабного процесса фильтрации неньютоновской жидкости в трехмерных периодических пористых средах методом асимптотической гомогенизации. Сформулированы так называемая локальная задача фильтрации в отдельной поре и локальное неньютоновско-вязкое определяющее соотношение. Разработан итерационный метод конечных элементов для решения локальной задачи в 1/8 ячейке периодичности, основанный на физической симметрии структуры. Рассчитаны распределение компонент скорости фильтрации, микрополей давления и неньютоновской вязкости в отдельной поре. На основе закона Дарси проанализирован нелинейный закон фильтрации, показано влияние реологических свойств жидкости на проницаемость.
Ключевые слова: метод асимптотической гомогенизации, неньютоновская жидкость, периодические пористые среды, закон фильтрации, реологические свойства
Введение. Течение неньютоновской жидкости в пористых средах играет важную роль во многих технических приложениях — получении полимерных композитов, технологии текстильных изделий и производстве бумаги. Таким образом, для оптимизации этих процессов требуется надлежащее описание течения через пористые среды [1-4].
Если жидкость можно считать ньютоновской, прямое применение известного закона Дарси [5] к пористым средам позволяет достаточно хорошо предсказать процесс течения. Основными параметрами пористой среды, определяющими процесс фильтрации, являются вязкость жидкости и компоненты тензора проницаемости. Для того чтобы избежать сложных лабораторных экспериментов, эти компоненты тензора проницаемости могут быть определены с использованием аналитических или численных методов. В литературе отображены теоретические работы по численным расчетам о течении Стокса и газа в периодических пористых средах [5-13]. Выполненные на определенной геометрии, они представляют собой интересные шаги к пониманию течения жидкостей и газов через более сложные пористые среды [14].
Поскольку неньютоновские жидкости обладают сложными реологическими свойствами, фильтрационные свойства неньютоновских жидкостей в пористых средах также стали довольно сложными. В проведенных теоретических исследованиях [4, 15-17] предложены
законы фильтрации для поперечного ползучего течения степенных жидкостей через волокнистую среду. Однако этих исследований недостаточно для формулировки закона, описывающего процесс фильтрации неньютоновской жидкости при произвольных сложных геометриях пористых сред.
В данной работе мы изучаем характеристики течения неньютоновской жидкости в трехмерной композитной структуре. Для моделирования течения в отдельных порах трехмерной пористой среды использовали метод асимптотической гомогенизации, закон фильтрации неньютоновской жидкости в пористых средах и влияние неньютоновской вязкости на проницаемость.
Цель данной работы — разработать физико-математическую модель локального движения транспорта неньютоновской несжимаемой жидкости в пористой композитной структуре, предложить численный алгоритм для локальных задач неньютоновских жидкостей в ячейке периодичности, проанализировать роль влияния реологических свойств на проницаемость.
Принятые допущения и геометрическая модель расчетной области. Принято, что жидкость представляет собой изотропную не-ньютоновско-вязкую несжимаемую среду. Процесс фильтрации считается изотермическим, плотность массовых сил полагается равной нулю.
а б в
Рис. 1. Геометрическая модель расчетной области:
а — периодическая пористая структура; б — ячейки периодичности; в — 1/8 ячейки периодичности
Рассмотрим модель пористой среды (рис. 1). Микроскопические упругие свойства этой ортогональной пористой структуры были изучены в работах [18-20]. Предполагается, что пористая структура характеризуется периодичностью, отсутствуют тупиковые поры, ячейки периодичности геометрически и физически симметричны относительно координатных плоскостей местной декартовой системы координат. Всю область рассматриваемой трехмерной ортого-
нальной композитной структуры (см. рис. 1), которая занята порами, будем обозначать V, а ее границу (поверхность раздела с твердым телом) — Е. Введем также следующие обозначения: V — ячейка периодичности; V^p — область, занятая порой; — граница поры
с твердым телом.
Математическая модель течения неньютоновских жидкостей.
Рассмотрим движение несжимаемой неньютоновской жидкости в пористой среде, описываемое системой уравнений [21, 22]:
V- v = 0; (1)
fdv ^
р| — + v-V®v l = -Vp + V-o,
U )
где р — плотность; v — скорость; p — давление; а — тензор вязких напряжений; ® — тензорное произведение; V — набла-оператор.
Дополним систему (1) определяющими соотношениями неньютоновской жидкости. Для этой цели могут быть рассмотрены разные модели, например модели AI и AV фойгтовских изотропных вязких сред [22]. Показано, что для изотропной неньютоновско-вязкой среды фойгтовского типа тензор вязких напряжений можно представить в виде квазилинейной функции от тензора скоростей деформации D [23]:
а = 2|D = |(v® v + V® vT ), (2)
где | — коэффициент неньютоновской вязкости, зависящий от I2 ( D) — второго инварианта тензора скорости деформации [24]:
I2 (D) = V2D - D. (3)
Рассмотрим модель Carreau [21, 22, 25], в которой зависимость коэффициента вязкости |( ) имеет степенной вид
n—1
_(л , 1 2Т2\ 2
10 -1
= (1 + 122 ) 2. (4)
Здесь |0 — вязкость с нулевой скоростью сдвига; — вязкость с бесконечной скоростью сдвига; X — постоянная времени, число Carreau; n — степенной индекс. Поэтому отклонение n от единицы указывает на степень отклонения от ньютоновского поведения. При n < 1 мы имеем псевдопластичные жидкости (истончение сдвига), при n > 1 — дилатантные жидкости (утолщение сдвига). Псевдопластичные жидкости характеризуются кажущейся вязкостью, кото-
эо
рая уменьшается с увеличением скорости сдвига, однако в жидкости дилатантов кажущаяся вязкость возрастает с увеличением скорости сдвига. Очевидно, что при uo = , À = 0 или n = 1 модель Carreau описывает ньютоновские жидкости.
На поверхностях раздела «твердое тело - жидкость» Z характеризует граничное условие прилипания
v| z= (5)
В начальный момент времени при t = to считается заданным давление
H=to = Po. (6)
Введем безразмерную форму:
- v - p - t - x xo - u
v = —; p = — ; t = -; x= —; to = -°; u = —, (7) vo Po to xo vo Uo
где vo — модуль вектора характерной скорости, м/с; po — характерное давление, Па; xo — характерный размер всей области среды, м; Uo — характерное значение вязкости; to — характерное время.
Систему уравнений для течения несжимаемой неньютоновской жидкости можно записать в безразмерной форме следующим образом (символ * далее опускаем):
V • v = o;
—(— + v •V® v | = -Vp + —^V^(uD) ; (8)
Eu Ut I F EuRe V ^
D = 2 (V® v + V® vT ) ;
n-1 _
2тЛ 2 • л = */?D • D • 5 =
u = 5 + (1 -5)(l + Cu2/22) 2 ; I2 = V2D••D;
v 1 Uo
P о v x
где Eu = —— число Эйлера; Re = o o o — число Рейнольдса;
0ovo uo
Cu = — число Carreau.
xo
Применение метода асимптотического осреднения (МАО) для постановки задачи на ячейке периодичности. Применение метода асимптотического осреднения для системы уравнений Навье-Стокса
описано в работах [1, 3, 12]. В рассматриваемом случае суть метода состоит в следующем. Пусть введены: ¡о — линейный размер ячейки
периодичности V$ среды, малый параметр к = — «1, два типа без
$
хо
размерных координат — глобальные х = — = (Х1, Х2, Х3 )
хо
хх
и локальные \ = — = — = ($1, $2, $3 ). Тогда все функции (/), описы-¡о к
вающие течение жидкости в порах, можно считать квазипериодическими, т. е. зависящими от времени I локальных (£) и глобальных (х) координат. Дифференцирование этих функций можно осуществлять с помощью правила
V/ х/ + -У$/. (9)
к
Условие квазипериодичности заключается в том, что функции / медленно изменяются по аргументу х1 на расстояниях Хо и являются периодическими относительно аргумента $ j, т. е.
/ ($1 - 2, £2, $3 ] = / ($1 + 2, £2, $3 ];
/ ($1, $2 -$з ] = / ($1, $2 + 2, $3 ]; (ю)
/ ($1, $2, $3 - 2) = / ($1, $2, $3 + 2].
Кроме того, вводится операция осреднения функций по области
(/> = -L | /¿V, (11)
Ф
Р V*
$р
где -р =| ¿V — пористость среды.
Тогда, выполняя осреднение р, V, получим
= V, (р) = р, (|д) = Ц, (12)
где V, р, | — средние значения функций.
Соотношения (12) следует понимать как дополнительные условия, предъявляемые к локальным параметрам течения в У^ либо как
обозначения для вычисленных осредненных функций.
Введем соотношения между безразмерными параметрами к, Ей, Яе и Си:
Еи ке = Еи0 = О (1); Яе к" = Яе0 = О (1);
Си0 кс = Си; Си0 = О (1), (13)
где е, г, с е2 — некоторые целые числа, которые соответствуют разным режимам течения неньютоновской жидкости.
В соответствии с общей концепцией МАО решение задачи (8) находим в виде асимптотических разложений по степеням малого параметра к:
V = у(0) (х, %) + ку(1) (х, %) + к2у(2) (х, %) +. -; р = р( 0) ( X, % ) + кр(1) ( X, % ) + к2 р( 2) ( X, % ) + • • •; (14)
| = |д(0) ( X, % ) + к|Д(1) ( X, % ) + к2|Д( 2) ( X, % ) + • .
Рассмотрим случай е = 2, г = 0 и с = 0, т. е. Еи » О (1), Яе = = О (1) и Си = О (1). Подставляя разложения (14) в уравнения (8),
с учетом правила дифференцирования квазипериодических функций (9) и соотношений (10)—(12) получаем локальную задачу нулевого уровня:
^(0) = 0; ^р(0) = 0; ^0) = 0; (р^,%,/)) = р. (15)
Л
Из уравнения (15) следует, что р(0)( X, %, I ) = р( 0)( X, I) не зависит от локальных координат %. Рассматривая следующий уровень, получаем локальную задачу «на ячейке периодичности» относительно
(0) (1) скорости нулевого уровня V ; и давления первого уровня р !:
^ ¿0) = 0; -V р(1)+от ^(0)»(0)И ^р(0);
Б(0) = 2 ^0)г); (16)
п-1
1(0) = 5 + (1 -5)(1 + Си02Г(0)) 2 ; 7(0) = 2В(0)..В(0);
V
(0)
= 0;
/V 0,
V
(0)"
= 0,
,(1)"
= 0.
Неизвестными в этой системе являются функции V*-0) и р(1), гра-^ (0)
диент давления V хр ! рассматривается как один из входных параметров. Символом обозначены условия периодичности, причем
- 2 < ^ < -2 — область ячейки периодичности .
Постановка локальных задач фильтрации с учетом трехмерной структуры пор. Рассмотрим далее трехмерную пористую структуру, у которой течение среды осуществляется вдоль одной из
трех осей О^. Поскольку локальная задача является нелинейной,
нельзя воспользоваться методом разделения переменных, как в работах [9, 19]. Однако, основываясь на применяемых методах [9, 12, 19], решение задач можно упростить, если воспользоваться следующей теоремой о продолжении решения.
Считается, что ячейка периодичности V трехмерной структуры
имеет зеркальную симметрию относительно координатной плоскости
Д0)
(1) (0) р ' и | ' локальных за-
О^2, О^3, Тогда решение V
дач (16) можно получить с помощью симметричного или асимметрична) ~ (а)
р ' и | ', опре-
/(а) = (»
ного продолжения функций V(а) = (, ^ , ^ деленных в 1/8 ячейки периодичности (в первом квадранте
: 0 < ^ < 2|) (рис. 2) и являющихся решениями следующих £(а) локальных задач формы компонента:
~(а)
Л. =
г[г
(а) р[// +
Еи0 Яе0
(а) 4а))
и
= Р^Ча;
(17)
4а) =2'
1 > + &)) 2\ г/-> 1/гГ
Е
5
2
и-1
Д (а) = ^ + (1 - ^ ) (1 + Си 02У (а)) ^; У (а) = 2Д(а) —
где ^ = дХ ' ^ = д&~ ' — символ Кр°некера-
© © © © © © © ©
ы © ы © ы © ©
© 0 © 0 © 0 © 0
0 © «а 0 © «а 0 © «а 0 © §1
Рис. 2. Изменение знаков функций при симметричном или асимметричном
продолжении для локальной задачи Ь
(1)
Используя принцип симметричного и антисимметричного продолжения решения системы (17), можно записать граничные условия на граничных плоскостях 1/8 ячейки периодичности У^, удовлетворяющие условиям периодичности системы (16):
*- = 0
* 1 и = 2:
+ -25ЛI 1 -(б* + -)
д)
(а)
= 0;
Р(% -«а') = 0;
(а)
а, /, - = 1,3.
(18)
Численное решение локальных задач. Решим локальную задачу Ь(а) (17) на основе метода конечных элементов [26], для чего используем 6-узловой тетраэдр с тремя степенями свободы по скоростям в каждом узле и с одной степенью свободы по давлению в вершинах. Поскольку локальная задача Ь(а) (17) — нелинейная, для ее решения применим итерационный метод как разновидность метода упругих решений [18-20]. Согласно этому методу определяющие соотношения в системе (17) для вязкости неньютоновских жидкостей линеаризуются:
= 2
с^Н = 2|(а){т-1} Б(а^т} =
и-1
02т?(а){т-1}\ 2
(1 - 5)(
5 + (1 - 5 )(1 + Си02Т
Б
(а){т}
(19)
где Т(а){т-1} = 2 ¿(а){т-1}Б(а){т-1}; ¿(а)М = I (у^П + уНМ^); ^НМ
/?(а){т}
и п ' — значения тензора неньютоновско-вязких напряжений с(а){0} и тензора скорости деформации на т-м шаге итераци-
~(а){т}
онного цикла; | п ' — вязкость неньютоновской жидкости на т-м шаге итерации.
та ~(а){т} ~(а){т} ~(а){0}
Обозначим у п ' и 1 — соответственно скорость у п '
и давление р( а){0} на т-м шаге итерационного цикла. Тогда на т-м шаге итерации вместо задачи (16) получим следующую линеаризованную задачу:
4а ){т}=0;
Р а){т} . 2 /П(а){т-1} б( а){т}) = ~(0), (20)
+ Еи0Яе0 V Б I . - Р,г дга ■ (20)
Проверяем сходимость на т-м шаге итерации, если
' X{т+1} - X{т})^/X{т}
< в, где в < 10 5 является очень малым чис-
„ _(/) ~(а) ~(а)
лом; X = Vг !, р !, | !, затем итерационный процесс завершается. Нелинейный закон фильтрации неньютоновской жидкости.
Из системы (16) видно, что решение (V'-0) и р(1)) локальной задачи
V (0)
зависит не только от входных данных V хр , но и от вязкости неньютоновской жидкости. Поэтому
V 0) = ^ VxР 0)). (21)
После решения серии задач £(а) (17) методом конечных элементов для всех а проинтегрируем скорости v( 0) по областям, занятым неньютоновской жидкостью. В результате получаем осредненную скорость в локальной задаче фильтрации (16). Тогда эффективный
закон фильтрации, связанный с этой скоростью и градиентом
V (0)
давления Vхр , можно записать в следующем виде:
да
у<°>) = ^ (V А.р(«)).
(22)
Поскольку ранее был указан алгоритм вычисления осредненных скоростей (у(0)\ по заданным значениям градиента давления Vхр(0),
то, следовательно, фактически был указан алгоритм нахождения значения функции (22). Функция (22) может быть указана, если имеется информация о типе геометрической симметрии ячейки периодичности композита и типе анизотропии волокон и матрицы. Ввиду того, что основное допущение симметрии принято, решение всех задач
£(а) (17) будет иметь указанный тип симметрии. Это относится и к функции (22). Но перечисленные выше и тождественные преобразования образуют группу ортотропии. Согласно [22], функция (22) будет тензорной функцией, индифферентной относительно группы ортотропии. Но тогда для такой функции можно использовать представление ее в тензорном базисе группы ортотропии, которое имеет вид
(у(0)) = £ф(а)('1р, 1Р, 1Р)еа или (у{°>) = £ф(а)(/1р,/2р, 1Р, (23)
а=1 а=1
где 1р — скалярные функции от инвариантов градиента давления
V хР(0),
I? =
др
(0)
дх;
(24)
»
Функции ф ' фактически представляют собой искомую эффективную проницаемость в пористых средах. Зная значения 0)^ и
Vхр, функции ф(а) можно вычислить. Для ортогональных полей запишем формулы (23) в явном виде:
ф(а)(I?, 12?, I? ) = Ф(а){
?{а)
Поэтому
V
(°)\ =
( 3
р(а)
г?(а)
) 1р=ф(а) р,{°).
V хр
\а=1
Л°)
(25)
(26)
где еа = еа ® еа — диадный базис; К = ^ ф(а)еа — тензор прони-
а=1
цаемости.
2
3
Из анализа локальных задач видно, что направление осредненной скорости противоположно направлению градиента давления, поэтому
¿0)\ = -КУ хр( 0),
(27)
где К = К
,(0)
V хр
(V хр(0)|)= К
= К) — тензорная функция, зависящая от модуля
= xР(0)V хр( 0) — градиента давления.
^(а)
Поскольку функция V ' для локальных координат имеет симметрию или асимметрию, например, при а = 1 функция для координатной плоскости (О^1^з и О^2^з) является асимметричной, поэтому К = 0, функция у3а) для этой координатной плоскости (О^1^2 и О^2^3) также асимметрична, поэтому К = 0. Следовательно, матрица К = б1а§(К/, К22, К33) будет диагональной. Стоит отметить, что микроструктура является изотропной, имеем К11 = = К22 = К33 = К.
Аналогично эффективную вязкость, связанную со осредненной вязкостью и градиентом давления ^.р(0), можно записать в форме функции
3
\| =
¿л(а)( Лр(а),/р(а), 13р(а))
!л(а)(
а=1 V
а=1
р( а)
р( а)
р( а)
) 1(а)
' а=1 ^
(28)
Для того чтобы изучить реологические свойства жидкостей, введем понятия средней вязкости и пористости:
|> = -
Ф Ф ^
Фр фр V а=1
3 ( Л
~ (а)
йУ;
(29)
Е
Ф = Г йУ = 8/ йУ = 8У Г йУ
р У К Уе
где Е — число конечных элементов.
Осредненная вязкость описывает макроскопическое свойство ре-алогии. Значение вязкости, отличающейся от единицы, характеризует
1
2
3
прочность неньютоновской жидкости. Очевидно, что осредненная вязкость ньютоновской жидкости равна единице.
Проверка адекватности математической модели на моделирование ньютоновских жидкостей. В качестве примера для численного моделирования мы используем бензол (С6Н6). Он входит в состав бензина, является исходным сырьем для производства лекарственных средств, различных пластмасс, синтетической резины, красителей.
Таблица 1
Основные параметры расчета
Параметр Значение Параметр Значение
^0, Па-с 6,52-Ш-4 У0, м/с 0,1
Си0 1 Х0, м 1
Рo, Па 105 ^ м 10-4
к 10-4 1 (ЕиЧе0 ) 6,52-10"2
Задача течения ньютоновской жидкости (п = 1) для определения
свойств композитных структур решается при
V ХР
= 1, пористос-
ти — 0,6230, коэффициенте проницаемости К = 0,0689.
Таблица 2
Результаты локальных задач !,(1) ньютоновской жидкости
Параметр Значение
минимальное максимальное
Р(1), Па -0,3529 0,3442
-(1) / , м/с 0 0,1436
-(1) / -2 , м/с -0,0774 0
-(1) / -3 , м/с 0 0,0776
д(1), Па-с 1 1
Для установления точности численного алгоритма сначала проверяют изотропность материала. Согласно закону Дарси, закон фильтрации ньютоновской жидкости удовлетворяет линейной зависимости от градиента давления V Хр(0) (рис. 3).
О 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 Рис. 3. Зависимость средней скорости от модуля
градиента макродавления
V ХР
(0)
для ньютоновской
жидкости в локальных задачах и ':
^ _ ¿0); —А— _ ^ _ ¿(3)
Данные рис. 3 согласуются с законом Дарси для изотропных пористых сред.
Результаты моделирования локальной задачи !,(1) псевдопластичной жидкости. Выполним численное решение локальной зада-
чи (17) при
V ХР
,(о)
= 1 и п = 0,25. Результаты расчетов фильтрации
неньютоновской жидкости приведены в табл. 3. По сравнению с результатами для ньютоновской жидкости (см. табл. 2) в процессе фильтрации псевдопластичной жидкости (п = 0,25) импульс давления немного увеличился, скорость течения возросла, вязкость значительно снизилась, пористость осталась неизменной.
Таблица3
Результаты локальных задач неньютоновской жидкости при п = 0,25
Параметр Значение
минимальное максимальное
Р(1), Па -0,3570 0,3450
-(1) , , м/с 0 0,2451
, м/с -0,1197 0
-(1) , -3 , м/с 0 0,1196
Д(1), Пас 0,3136 0,9999
Распределение компонентов скорости V
¡(1)
~(1)
давления р' и вяз-
кости Д 7 неньютоновской жидкости показаны на рис. 4. На этом и других аналогичных рисунках введены обозначения для локальных координат = х, = у, = z.
1
0,20 0,4
0,15 0,3
0,10 0,2 ОД
0,05 0
0 г
0,10 0,4
0,08 0,3
0,06 0,2
0,04 ОД
0,02 0
0
-0,01 -0,02 -0,03 -0,04 -0,05 -0,06 -0,07 -0,08 -0,09 -0,10 -0,11
Рис. 4. Результаты решения локальной задачи -(') б -(')
а — компоненты скорости у ; б — компоненты скорости -2 ; в — компоненты скорости
Л1).
~ (!) Л ~
г — компоненты давления р ; о — компоненты вязкости Д
(1)
3
Влияние реологических свойств на проницаемость. Исследуем
связь между средней скоростью (у(а)), градиентом давления Vхру"! и
V хР
(0)
М-
степенным индексом п, рассмотрим п е [0,25,1,75] и Зависимости между средней скоростью (у(1) ) и градиентом давления
V хР
,(0)
при разных степенных индексах п показаны на рис. 5.
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 |УХ^0)|
Рис. 5. Зависимость средней скорости ()} от градиента
давления
V хР
,(0)
и степенного индекса п для локальной
задачи Ё
(1)
■ — п = 0,25; -•--п = 0,50; -•■ — п = 0,75; -•--п = 1,00;
---п = 1,25; -•- — п = 1,50; ---п = 1,75
Данные, приведенные на рис. 5, показывают, что осредненная скорость \у(1)) уменьшается с увеличением степенного индекса п.
При увеличении модуля градиента давления
V хР
( 0)
средняя ско-
рость возрастает, эффект степенного индекса также увеличивается. Нелинейное соотношение, приведенное на рис. 5, соответствует закону фильтрации (26) неньютоновской жидкости.
На рис. 6 показана взаимосвязь между средней вязкостью (Д
и модулем
V хР
(0)
градиента давления при разных степенных индек-
сах п. Данные, приведенные на рис. 6, показывают, что средняя вязкость (д(1)) возрастает с увеличением степенного индекса. При уве-
личении модуля
V хР
,(0)
градиента давления влияние степенного
индекса на среднюю вязкость также увеличивается.
0,70
' 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 |Ух/>>|
Рис. 6. Зависимость средней вязкости ( д(1) ) от градиента
давления
V хР
(0)
и степенного индекса п для локальной
задачи П
(1)
■ — п = 0,25; -•- — п = 0,50; - — п = 0,75; -•- — п = 1,00; ---п = 1,25; -•--п = 1,50; -.--п = 1,75
Заключение. Предложенная методика позволяет вычислять распределение микрополей давления, вязкости неньютоновской жидкости и компонент скорости фильтрации в пределах отдельной поры, а также рассчитывать основные параметры пористой среды (пористости и коэффициентов проницаемости) без проведения каких-либо дополнительных эмпирических исследований. Локальная задача, полученная на основе метода асимптотического осреднения, решается методом конечных элементов. Результаты моделирования выявляют роль реологии жидкости. Показано, что средняя скорость фильтрации
уменьшается с увеличением индекса степенного закона, при увеличении модуля градиента макродавления средняя скорость фильтрации возрастает и эффект степенного индекса также увеличивается. Кроме того, средняя вязкость жидкости возрастает с увеличением степенного индекса, при увеличении модуля градиента макродавления влияние степенного индекса на среднюю вязкость и фильтрацию также увеличивается.
ЛИТЕРАТУРА
[1] Бардзокас Д.И., Зобнин А.И. Математическое моделирование физических процессов в композиционных материалах периодической структуры. Москва, Едиториал УРСС, 2003, 376 с.
[2] Бахвалов Н.С., Панасенко Г.П. Осреднение процессов в периодических средах. Математические задачи механики композиционных материалов. Москва, Наука, 1984, 352 с.
[3] Dimitrienko Yu.I. Thermomechanics of composites under high temperatures. Springer, 2015, 434 p.
[4] Idris Z., Orgeas L., Geindreau C., et al. Microstructural effects on the flow law of power-law fluids through fibrous media. Modelling and Simulation in Materials Science and Engineering, 2004, vol. 12, pp. 995-1015.
[5] Димитриенко Ю.И., Захарова Ю.В., Богданов И.О. Математическое и численное моделирование процесса фильтрации связующего в тканевом композите при RTM-методе изготовления. Университетский научный журнал, 2016, № 19, с. 33-43.
[6] Димитриенко Ю.И., Шпакова Ю.В., Богданов И.О., Сборщиков С.В. Моделирование процесса многоуровневой фильтрации жидкого связующего в тканевом композите при RTM-методе изготовления. Инженерный журнал: наука и инновации, 2015, № 12, с. 14-54.
[7] Димитриенко Ю.И., Богданов И.О. Многомасштабное моделирование процессов фильтрации жидкого связующего в композитных конструкциях при RTM-методе изготовления. Математическое моделирование и численные методы, 2017, № 2, с. 3-27.
[8] Димитриенко Ю.И., Иванов М.Ю. Моделирование нелинейных динамических процессов переноса в пористых средах. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, № 1, 2008, с. 39-56.
[9] Димитриенко Ю.И., Левина А.И., Боженик П. Конечно-элементное моделирование локальных процессов переноса в пористых средах. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, 2008, № 3, с. 90-104.
[10] Димитриенко Ю.И., Богданов И.О. Многомасштабное моделирование процессов фильтрации в пористых средах. Инженерный журнал: наука и инновации, 2018, № 3,с. 17-38.
[11] Dimitrienko Yu.I. Dynamic transport phenomena in porous polymer materials under impulse thermal effects. Transport in porous media, 1999, vol. 35, pp. 299-326.
[12] Dimitrienko Yu.I., Dimitrienko I.D. Simulation of local transfer in periodic porous media. European Journal of Mechanics/B-Fluids, 2013, № 1, pp.174-179.
[13] Wang J.G., Leung C.F., Chow Y.K. Numerical solutions for flow in porous media. International Journal for numerical and analytical methods in geomecha-nics, 2003, vol. 27, pp. 565-583.
[14] Koponen A., Kandhai D., Hellen E., Alava M., Hoekstra A., Kataja M., Niskanen N., Sloot P., Timonen J. Permeability of three-dimensional random fiber webs. Physical Review Letters, 1998, vol. 80, pp. 716-719.
[15] Auriault J.L., Royer P., Geindreau C. Filtration law for power-law fluids in anisotropic porous media. International Journal of Engineering Science, 2002, vol. 40, pp. 1151-1163.
[16] Orgeas L., Idrisl Z., Geindreau C., et al. Modelling the flow of power-lawfluids through anisotropic porous media at low-pore Reynolds number. Chemical Engineering Science, 2006, vol. 61, pp. 4490-4502.
[17] Orgeas L., Geindreau C., Auriault J.L., et al. Upscaling the flow of generalised Newtonian fluids through anisotropic porous media. Journal of Non-Newtonian Fluid Mechanics, 2007, vol. 145, pp. 15-29.
[18] Димитриенко Ю.И., Соколов А.П. Об упругих свойствах композиционных материалов. Математическое моделирование, 2009, № 4, с. 96-110.
[19] Димитриенко Ю.И., Кашкаров А.И., Макашов А.А. Конечно-элементный расчет эффективных упругопластических характеристик композитов на основе метода асимптотического осреднения. Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки, № 1, 2007, с. 26-46.
[20] Димитриенко Ю.И., Губарева Е.А., Сборщиков С.В. Многомасштабное моделирование упругопластических композитов с учетом повреждаемости. Математическое моделирование и численные методы, 2016, № 2, с. 3-23.
[21] Димитриенко Ю.И., Ли Шугуан. Конечно-элементное моделирование неизотермического стационарного течения неньютоновской жидкости в сложных областях. Математическое моделирование и численные методы, 2018, № 2, с. 70-95.
[22] Димитриенко Ю.И. Нелинейная механика сплошной среды. Москва, Физ-матлит, 2009, 624 с.
[23] Димитриенко Ю.И. Механика сплошной среды. Т. 2: Универсальные законы механики и электродинамики сплошных сред. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2011, 560 с.
[24] Dimitrienko Yu.I. Tensor analysis and nonlinear tensor functions. Springer, 2002, 662 p.
[25] Bird R.B., Armstrong R.C., Hassager O. Dynamics of polymeric liquids. Vol. 1: Fluid mechanics. John Wiley & Sons, 1987, 649 p.
[26] Zienkiewicz O.C., Taylor R.L., Zhu J.Z. The Finite Element Method: Its Basis and Fundamentals. Elsevier, 2005, 733 p.
Статья поступила в редакцию 24.08.2019
Ссылку на эту статью просим оформлять следующим образом: Димитриенко Ю.И., Шугуан Ли. Моделирование проницаемости неньютоновских жидкостей в трехмерных композитных структурах на основе метода асимптотической гомогенизации. Математическое моделирование и численные методы. 2019. № 3. c. 19-38.
Димитриенко Юрий Иванович — д-р физ.-мат. наук, профессор, заведующий кафедрой «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана, директор Научно-образовательного центра «Суперкомпьютерное инженерное моделирование и разработка программных комплексов» МГТУ им. Н.Э. Баумана. e-mail: dimit.bmtstu@gmail.com
Шугуан Ли — аспирант кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. e-mail: shuguangli2008@gmail.com
Modeling non-Newtonian liquid permeability in three-dimensional composite structures on the base of asymptotic homogenization method
© Y.I. Dimitrienko, Sh. Li Bauman Moscow State Technical University, Moscow, Russian Federation
The paper describes investigating a mathematical model of the process of a non-Newtonian liquid multiscale filtration in three-dimensional periodic porous media by asymptotic homogenization. The so-called local problem of filtration in a single pore is formulated as well as the local non-Newtonian-viscous defining relationship. An iterative finite element method is developed for solving a local problem in 1/8 periodicity cell, based on the physical symmetry of the structure. The distribution of the components of the filtration rate, pressure micro-fields and non-Newtonian viscosity in a single pore is calculated. On the basis of Darcy's law the nonlinear filtration law is analyzed, the effect of liquid rheological properties on permeability is shown.
Keywords: asymptotic homogenization method, non-newtonian liquid, periodic porous media, filtration law, rheological properties
REFERENCES
[1] Bardzokas D.I., Zobnin A.I. Matematicheskoe modelirovanie fizicheskikh protsessov v kompozitsionnykh materialakh periodicheskoy struktury [Mathematical modeling of physical processes in composite materials of periodic structure]. Moscow, Editorial URSS Publ., 2003, 376 p.
[2] Bakhvalov N.S., Panasenko G.P. Osrednenie protsessov v periodicheskikh sredakh. Matematicheskie zadachi mekhaniki kompozitsionnykh materialov [Averaging processes in periodic media. Mathematical problems of the composite material mechanics]. Moscow, Nauka Publ., 1984, 352 p.
[3] Dimitrienko Yu.I. Thermomechanics of composites under high temperatures. Springer, 2015, 434 p.
[4] Idris Z, Orgeas L, Geindreau C, Bloch, J.F., Auriault, J.L Modelling and Simulation in Materials Science and Engineering, 2004, vol. 12, pp. 995-1015.
[5] Dimitrienko Yu.I., Zakharova Yu.V., Bogdanov I.O. Humanities and Science University Journal, 2016, no. 19, pp. 33-43.
[6] Dimitrienko Yu.I., Shpakova Yu.V., Bogdanov I.O., Sborshchikov S.V. Inzhe-nernyy zhurnal: Nauka i innovatsii — Engineering Journal: Science and Innovation, 2015, no. 12, pp. 14-54.
[7] Dimitrienko Yu.I., Bogdanov I.O. Matematicheskoe modelirovanie i chislennye menody — Mathematical Modeling and Computational Methods, 2017, no. 2, pp. 3-27.
[8] Dimitrienko Yu.I., Ivanov M.Yu. Vestnik MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2008, no. 1, pp. 39-56.
[9] Dimitrienko Yu.I., Levina A.I., Bozhenik P. Vestnik MGTU im. N.E. Baumana. Seria Estestvennye nauki - Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2008, no. 3, pp. 90-104.
[10] Dimitrienko Yu.I., Bogdanov I.O. Inzhenernyy zhurnal: Nauka i innovatsii -Engineering Journal: Science and Innovation, 2018, no. 3, pp. 17-38.
[11] Dimitrienko Yu.I. Transport in porous media, 1999, vol. 35, pp. 299-326.
[12] Dimitrienko Yu.I., Dimitrienko I.D. European Journal of Mechanics/B-Fluids, 2013, no. 1, pp.174-179.
[13] Wang J.G., Leung C.F., Chow Y.K. International Journal for Numerical and Analytical Methods in Geomechanics, 2003, vol. 27, no. 7, pp. 565-583.
[14] Koponen A., Kandhai D., Hellen E., Alava M., Hoekstra A., Kataja M., Niskanen N., Sloot P., Timonen J., Physical Review Letters, 1998, vol. 80, no. 4, pp. 716-719.
[15] Auriault J.L, Royer P., Geindreau C. International Journal of Engineering Science, 2002, vol. 40, pp. 1151-1163.
[16] Orgeas L., Idrisl Z., Geindreau C., Bloch, J.F., Auriault, J.L. Chemical Engineering Science, 2006, vol. 61, pp. 4490-4502.
[17] Orgeas L., Geindreau C., Auriault J.L., Bloch, J.F., Auriault, J.L. Journal of Non-Newtonian Fluid Mechanics, 2007, vol. 145, pp. 15-29.
[18] Dimitrienko Yu.I., Sokolov A.P. Matematicheskoe modelirovanie — Mathematical Models and Computer Simulations, 2012, vol. 24, no. 5, pp. 3-20.
[19] Dimitrienko Yu.I., Kashkarov A.I., Makashov A.A. VestnikMGTU im. N.E. Baumana. Seria Estestvennye nauki — Herald of the Bauman Moscow State Technical University. Series: Natural Sciences, 2007, no. 1, pp. 26-46.
[20] Dimitrienko Yu.I., Gubareva E.A., Sborshchikov S.V. Matematicheskoe mode-lirovanie i chislennye menody — Mathematical Modeling and Computational Methods, 2016, no. 2, pp. 3-23.
[21] Dimitrienko, Yu.I. Li S. Matematicheskoe modelirovanie i chislennye menody — Mathematical Modeling and Computational Methods, 2018. № 2. pp. 70-95.
[22] Dimitrienko Yu.I. Nelineinaya mekhanika sploshnoi sredy [Nonlinear Continuum Mechanics]. Moscow, Fizmatlit Publ., 2009, 624 p.
[23] Dimitrienko Yu.I. Mekhanika sploshnoy sredy. Tom 2. Universalnye zakony mekhaniki i elektrodinamiki sploshnoy sredy [Continuum Mechanics. Vol. 2. Universal laws of continuous media mechanics and electrodynamics]. Moscow, BMSTU Publ., 2011, 560 p.
[24] Dimitrienko Yu.I. Tensor analysis and nonlinear tensor functions. Springer, 2002, 662 p.
[25] Bird R.B, Armstrong R.C., Hassager O. Dynamics of polymeric liquids. Vol. 1: Fluid mechanics. John Wiley & Sons Publ., 1987, 649 p.
[26] Zienkiewicz O.C., Taylor R.L., Zhu J.Z., The Finite Element Method: Its Basis and Fundamentals. Elsevier Publ., 2005, 733 p.
Dimitrienko Y.I., Dr. Sc. (Phys.-Math.), Professor, Head of the Department of Computational Mathematics and Mathematical Physics, Director of Research and Education Center of Supercomputer Engineering Modeling and Software Packages (Simplex), Bauman Moscow State Technical University. e-mail: dimit.bmtstu@gmail.com.
Li Sh., post-graduate Student, Department of Computational Mathematics and Mathematical Physics, Bauman Moscow State Technical University. e-mail: shuguangli2008@gmail.com.