CC BY
22
www.volsu.ru
001: https://doi.Org/10.15688/jvolsu1.2017.2.5
УДК 532.516:532.546 ББК 22.253.3
ПРОНИЦАЕМОСТЬ ДВУМЕРНОЙ ПОРИСТОЙ СРЕДЫ ИЗ ВОЛОКОН КВАДРАТНОГО СЕЧЕНИЯ (ЯЧЕЕЧНАЯ МОДЕЛЬ)
Екатерина Владимировна Мосина
Кандидат физико-математических наук, старший преподаватель кафедры фундаментальной информатики и оптимального управления, Волгоградский государственный университет katefrkate@mail.ru, fiou@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Игорь Викторович Чернышев
Кандидат физико-математических наук, доцент кафедры фундаментальной информатики и оптимального управления, Волгоградский государственный университет igor_chernyshev@mail.ru, fiou@volsu.ru
просп. Университетский, 100, 400062 г. Волгоград, Российская Федерация
Аннотация. В рамках ячеечной модели численно решена задача о поперечном обтекании системы волокон (стержней) квадратного сечения. Решение для отдельной ячейки получено в приближении малых чисел Рейнольдса с краевыми условиями Хаппеля (равенство нулю нормальной составляющей скорости и касательного напряжения) на границе ячейки. Вычислена гидродинамическая проницаемость такой модельной волокнистой пористой среды. 3 Проведено сопоставление с результатами других авторов для систем стержней
., квадратного и круглого сечения в широком диапазоне изменения концентра-
ции включений. Предложены аналитические зависимости, наиболее адекватно ш описывающие проницаемость системы квадратных волокон.
§ Ключевые слова: модельная пористая среда, волокна квадратного сече-
& ния, приближение Стокса, ячеечная модель Хаппеля, проницаемость.
™ Введение
га к
§ При математическом моделировании фильтрационных задач движения жидкости
сквозь пористые среды зачастую используют уравнение для усредненной (фильтраци-@ онной) скорости, полученное на основе эмпирического закона Дарси. Более точный
микроскопический анализ течения сквозь пористую среду — задача сложная и требует дополнительных упрощений: линеаризации гидродинамических уравнений, предположения об одинаковости размеров частиц, пренебрежения броуновским движением и др. [3;7]. Важным макроскопическим параметром, определяющим характер течения и результирующее гидравлическое сопротивление, является коэффициент проницаемости, зависящий от пористости среды, структуры порового пространства, физико-химических свойств материала среды.
В настоящей работе решена задача о течении вязкой несжимаемой жидкости сквозь пористую среду, моделируемую совокупностью одинаковых стержней квадратного поперечного сечения. Такая задача может быть сведена к задаче обтекания одиночного стержня в рамках ячеечной модели [7], то есть решается краевая задача о течении в жидкой оболочке, окружающей одиночный стержень. При этом пористость для отдельной ячейки совпадает со средней пористостью всей моделируемой проницаемой среды.
Этот метод может быть эффективно применен к сравнительно концентрированным системам (с объемной долей твердых включений ф ~ 0,5 ^ 0,6). Ячеечные модели дают удовлетворительное приближение к усредненной картине течения вблизи твердых частиц в реальной физической системе, но, тем не менее, нельзя рассчитывать на то, что они хорошо описывают гидродинамические поля вблизи воображаемых границ ячеек. Вопрос о выборе формы и размеров жидкой ячейки, а также условий на ее границе, остается открытым и до сих пор привлекает исследователей (например, [15]).
Проведено сопоставление результатов, полученных в данной работ, с аналитическими [9; 11], полуаналитическими [10; 13; 14], расчетными [8; 12] и экспериментальными [13] данными других авторов, определены наиболее оптимальные аналитические зависимости для коэффициента проницаемости системы квадратных стержней в различных диапазонах объемной концентрации частиц.
1. Поперечное обтекание квадратного стержня (ячеечная модель)
Рассматривается стационарное двумерное однородное течение несжимаемой ньютоновской вязкой жидкости со средней скоростью и через ячейку квадратного сечения размера 2ЬН с расположенным в ее центре квадратным включение (стержнем) размера 2аН (см. рис. 1). Свяжем с жидкой ячейкой декартову систему координат: ось Ох в направлении скорости и, ось Оу перпендикулярно ей, ось Ог поперек потоку и вдоль образующих стержня. После обезразмеривания пространственных переменных х и у на Н, компонент скорости и, V на и и давления р на ци/Н (где ц — коэффициент динамической вязкости жидкости) система уравнений движения в приближении Стокса примет вид
Дv = Чр, V- V = 0. (1)
На поверхности включения задаются граничные условия прилипания у = 0; на границах жидкой ячейки — условия Хаппеля: непроницаемости (равенства нулю нормальной компоненты скорости ьп = 0) и равенства нулю касательного напряжения т =0.
Исключая из системы (1) давление р и учитывая двумерность постановки, задача сводится к бигармоническому уравнению на функцию тока
ДДф = 0. (2)
Линейность уравнения движения, а также симметрия постановки задачи позволяют ограничить рассмотрение четвертью ячейки (х,у) € [0,6] х [0,6] (см. рис. 2). Краевые
условия для функции тока на границах расчетной области будут следующими. Условия прилипания на твердой поверхности включения:
0 < У < а (с = Ъ — а) (3)
с < х < Ь. (4)
Мс,У) = 0, 4х{с,у) = 0 ф(ж, а) = 0, -фуа) = 0,
Симметрия относительно у = 0 и х = Ь приводит к условиям:
ф(х, 0) = 0, фуу(х, 0) = 0, 0 < ж < с, фл(6, у) = 0, фххх(Ь, у) = 0, а < у < Ь.
Граничные условия на поверхности жидкой ячейки:
Ф(0, у) = у/Ь, фхх(0,у) = 0, ф(х,Ь) = 1, фуу (х,Ь) = 0,
Для наглядности условия (3)-(8) вынесены на рисунок 2.
0 < у < Ь, 0 < х<Ь.
(5)
(6)
(7)
(8)
Рис. 1. Жидкая ячейка с твердым включением
Рис. 2. Расчетная область
Рис. 3. Линии тока при ф = 0,1
2. Функция тока и сила сопротивления
Краевая задача для функции тока решена численно с помощью симметричной конечно-разностной аппроксимации бигармонического уравнения (2) и граничных условий (3)-(8) второго порядка точности [1]. Для достижения абсолютной ошибки Ю-3 гидродинамических полей разбиение расчетной области обычно проводилось на 300 отрезков по каждому направлению. Система разностных уравнений решена методом последовательной верхней релаксации с параметром релаксации 1, 5 <ш < 2, 0. Расчеты представлены для широкого диапазона концентрации твердых включений ф = а2/Ь2 = = 0, 001 ^ 0, 8. Характерная картина линий тока изображена на рисунке 3.
По функции тока находим распределение давления в рассматриваемой области. Поскольку рх = (Дф)у, ру = — (Дф)ж, то, полагая р(0,у) = 0, получаем
х
Р(Х,у)= фХу (х,у)+У фууу (1,у) (Ц.
0
Сила сопротивления (обезразмеренная на циН) для четверти стержня равна сумме, проинтегрированных по соответствующим площадкам, нормального напряжения на
левой границе (х = с, 0 < у < а) и касательного напряжения на верхней границе
(■у = а, с < х < Ь)
4Fx = - [-р(с,у) + 24ху(с, у)] dy + [4уу(х,а) - ^хх(х,а)} dx
4хУ(с,у) - ^ууу(х,у) dx
dy + [4уу(х,а) - 4хх(х,а)] dx.
Из условий прилипания на твердых границах (3), (4) следует фху(с,у) = 0 при 0 < < у < а, фхх(х,а) = 0 при с < х < Ь, а из условия симметрии (5) фуу (х, 0) = 0 при 0 < х < с, тогда выражение для силы сопротивления будет следующим
Fx = 4 4уу (х,а) dx.
(9)
То есть расчет силы сопротивления твердого включения сводится к интегрированию касательного напряжения на уровне у = а.
ь
а
b
а
с
b
3. Проницаемость системы квадратных стержней
В рамках ячеечной модели сила сопротивления, отнесенная к площади жидкой ячейки, численно равна падению давления на длине L, обусловленному протеканием жидкости сквозь систему включений
Fх = Ар Ь2 L .
Полагая объемную долю ф твердого включения в жидкой ячейке равной объемной доли твердой фазы во всей совокупности включений и объединяя полученное выражение для силы сопротивления с законом Дарси, получим выражение для коэффициента проницаемости системы стержней (обезразмерен на Н2)
-2 = А. (10)
а2 Ьхф
На рисунке 4 изображена зависимость проницаемости квадратных стержней ks/a2, полученной численно по формуле (10) в рамках ячеечной модели «квадрат в квадрате», от проницаемости круговых стержней kc/R2 (R — радиус поперечного сечения круговых стержней), задаваемой аналитически [7], в рамках ячеечной модели «круг в круге». Сопоставление проведено в широком диапазоне концентрации ф и замечено, что при 0,1 < ф < 0, 6 эта зависимость линейная [2]
ь = ± £ = 1 ± (U - 1-ф), (11)
а2 а2 К2 а2 8ф V Ф 1 + Ф2/
где коэффициент пропорциональности а2 ~ п/4 получен методом наименьших квадратов и равен отношению площадей кругового и квадратного сечения одинакового размера.
В работах [4-6] проведено прямое численное моделирование течения вязкой несжимаемой жидкости в плоском канале, частично заполненном регулярной решеткой квадратных стержней. В качестве вспомогательного результата рассчитаны значения проницаемости такой системы стержней путем усреднения микроскопических гидродинамических полей. В ходе сравнения полученных данных для квадратных стержней с имеющимися в литературе данными для круговых стержней обнаружено, что для разреженной системы введением коэффициента формы а2 можно использовать аналитическую зависимость
к 11/1 \
^ = — 1,476+ 2ф , ф < 0, 2, (12)
а2 а2 8ф \ Ф )
полученную методом сингулярностей для круговых цилиндров [11]. Справедливость этой формулы для круговых цилиндров подтверждается рядом исследований и для более плотных упаковок (ф < 0, 6) [10; 12; 13].
Отметим, что для разреженной системы волокон тип их упаковки (квадратная, треугольная или гексагональная) практически не влияет на проницаемость [11; 13;14]. В случае стесненного расположения твердых включений влияние упаковки уже более существенно, что связано с различием извилистости поровых каналов.
Для умеренно разреженных решеток волокон хорошей аппроксимацией является полуэмпирическая формула Кармана — Козени к = (1 — ф)К2н/ск, где = (Ь2 — а2)/2а — гидравлический радиус жидкой ячейки, с коэффициентом Козени ск, рассчитанным на основе ячеечной модели Хаппеля двух концентрических круговых цилиндров [7]. Некоторые исследователи [8; 13] отмечают, что ячеечная модель дает завышенные значения проницаемости по сравнению с их данными. Для стержней квадратного сечения наблюдается такая же картина. Использование корректирующего множителя 1/а2 [5] вносит ошибку более 25% в вычисление проницаемости. С целью уменьшения этой ошибки следует использовать меньший коэффициент 8/9 [6]
к 1 / 1 1 — ф2\
¡¡2 = 9ф (1п ф — т—У • 01 < ф < 06 (13)
Для больших значений объемной концентрации коэффициент Козени ск следует брать равным 6 из интервала 4,5 < ск < 6, соответствующего экспериментальным данным для реальных пористых сред [7]. Выражение для проницаемости тогда примет вид
к (1 — ф)3
^ = • 0 6« ф < 08 (14)
На рисунке 5 представлены значения коэффициента проницаемости к/а2, полученные в результате прямого численного моделирования. Данные практически совпадают с данными Ванга [14] для проницаемости бесконечной упорядоченной совокупности квадратных стержней, полученной с использованием полуаналитического метода коллокации для функции тока.
Таким образом, формулы (12)-(14) дают хорошее приближение для коэффициента проницаемости системы квадратных стержней, в то время как ячеечная модель с условиями Хаппеля на границе жидкой ячейки дает завышенные значения проницаемости по сравнению с прямыми численными расчетами.
Рис. 4. Проницаемости системы квадратных кь/а2 и круговых кс/В? стержней, вычисленные по ячеечной модели (•), и линейная зависимость (11) (сплошная линия)
Рис. 5. Проницаемость системы квадратных стержней. Точки — численные значения: (о) — прямые численные расчеты [4-6]; (+) — работа [14]; (о) — работа [8] для ф < 0, 85. Кривые: 1 — формула (12) для ф < 0,1; 2 — формула (11) с коэффициентом формы а2 = п/4 для 0,1 < ф < 0, 6; 3 — формула (13) для 0,1 < ф < 0, 6; 4 — формула (14) для 0, 6 < ф < 0, 8; 5 — работа [14] для ф < 0, 03; 6 — работа [14] для ф > 0, 6
Заключение
В работе предложены зависимости для коэффициента проницаемости модельной пористой среды с волокнами (стержнями) квадратного сечения в широком диапазоне объемной концентрации включений. Несмотря на значительное различие формы поперечных сечений (квадрата и круга), обнаружено, что для вычисления проницаемости системы квадратных стержней могут быть использованы те же зависимости, что и для круговых цилиндров, но требуются некоторые корректирующие коэффициенты, учитывающие различия формы поперечного сечения. Полученные аппроксимационные выражения могут быть полезны в дальнейшем при решении макроскопических фильтрационных задач обтекания волокнистых пористых сред.
СПИСОК ЛИТЕРАТУРЫ
1. Андерсон, Д. Вычислительная гидромеханика и теплообмен / Д. Андерсон, Дж. Таннехил, Р. Плетчер. — М. : Мир, 1990. — Т. 2. — 392 с.
2. Кирдан, П. С. Коэффициент проницаемости двумерной решетки твердых частиц в медленном поступательном потоке : магистерская диссертация по направлению подготовки магистров «Прикладная математика и информатика» / Кирдан Павел Сергеевич. — Волгоград, 2012. — 27 с.
3. Лейбензон, Л. С. Движение природных жидкостей и газов в пористой среде / Л. С. Лейбензон. — М. ; Л. : Гостехиздат, 1947. — 245 с.
4. Мосина, Е. В. Течение жидкости в окрестности пористой границы / Е. В. Мосина, И. В. Чернышев // Вестник Нижегородского университета им. Н.И. Лобачевского. — 2011. — Вып. 4. — № 3. — C. 999-1001.
5. Мосина, Е. В. Условие скольжения на поверхности модельной волокнистой пористой среды / Е. В. Мосина, И. В. Чернышев // Письма в ЖТФ. — 2009. — Т. 35, вып. 5. — C. 103-110.
6. Мосина, Е. В. Численное исследование течения на границе жидкость — пористая среда / Е. В. Мосина // ТОХТ. — 2010. — Т. 44, вып. 5. — C. 536-542.
7. Хаппель, Дж. Гидродинамика при малых числах Рейнольдса / Дж. Хаппель, Г. Бреннер. — М. : Мир, 1976. — 631 с.
8. Analytical determination of the permeability for slow flow past periodic arrays of cylinders with different cross sections / A. C. Bayta§, D. Erdem§, H. Acar, O. Cetiner, H. Ba§ci // J. Porous Med. — 2012. — Vol. 15. — P. 1009-1018. — DOI: http://dx.doi.org/10.1615/JPorMedia.v15.i11.20.
9. Drummond, J. E. Laminar viscous flow through regular arrays of parallel solid cylinders / J. E. Drummond, M. I. Tahir // Int. J. Multiphase Flow. — 1984. — Vol. 10. — P. 515-540. — DOI: http://dx.doi.org/10.1016/0301-9322(84)90079-X.
10. Hellou, M. Stokes flow through microstructural model of fibrous media / M. Hellou, J. Martinez, M. El Yazidi // Mech. Res. Commun. — 2004. — Vol. 31. — P. 97-103. — DOI: http://dx.doi.org/10.1016/S0093-6413(03)00081-8.
11. James, D. F. Flow at the interface of a model fibrous porous medium / D. F. James, A. M. J. Davis // J. Fluid Mech. — 2001. — Vol. 426. — P. 47-72. — DOI: http://dx.doi.org/10.1017/S0022112000002160.
12. Tahir, M. A. Influence of fiber orientation on the transverse permeability of porous media / M. A. Tahir, H. V. Tafreshi // Phys. Fluids. — 2009. — Vol. 21. — Article ID: 083604-1-083604-5. — DOI: http://dx.doi.org/10.1063/1.3211192.
13. Tamayol, A. Transverse permeability of fibrous porous media / A. Tamayol, M. Bahrami // Phys. Review. — 2011. — Vol. 83 (E). — Article ID: 046314-1-046314-9. — DOI: http://dx.doi.org/10.1016/j-.ijheatmasstransfer.2008.09.032.
14. Wang, C. Y. Stokes flow through an array of rectangular fibers / C. Y. Wang // Int. J. Multiphase Flow. — 1996. — Vol. 22. — P. 185-194. — DOI: http://dx.doi.org/10.1016/0301-9322(95)00063-1.
15. Yadav, P. K. Slow motion of a porous cylindrical shell in a concentric cylindrical cavity / P. K. Yadav // Meccanica. — 2013. — Vol. 48. — P. 1607-1622. — DOI: http://dx.doi.org/10.1007/s11012-012-9689-0.
REFERENCES
1. Anderson D., Tannehill J., Pletcher R. Vychislitelnaya gidromekhanika i teploobmen [Computational Fluid Mechanics and Heat Transfer]. Moscow, Mir Publ., 1990, vol. 2. 392 p.
2. Kirdan P.S. Koeffitsient pronitsaemosti dvumernoy reshetki tverdykh chastits v medlennom postupatelnom potoke: magisterskaya dissertatsiya po napravleniyu podgotovki magistrov «Prikladnaya matematika i informatika» [The Coefficient of Permeability of Two Dimentional Array of Rigit Particles in Slow Parallel Flow: Master Thesis for Degree on Applied Mathematics and Computer Science]. Volgograd, 2012. 27 p.
3. Leybenzon L.S. Dvizhenie prirodnykh zhidkostey i gazov v poristoy srede [The Motion of Nature Fluid and Gas Through Porous Medium]. Moscow, Leningrad, Gostekhizdat Publ., 1947. 245 p.
4. Mosina E.V., Chernyshev I.V. Techenie zhidkosti v okrestnosti poristoy granitsy [Fluid Flow Near the Porous Boundary]. Vestnik Nizhegorodskogo universiteta im. N.I. Lobachevskogo, 2011, iss. 4, no. 3, pp. 999-1001.
5. Mosina E.V., Chernyshev I.V. Uslovie skolzheniya na poverkhnosti modelnoy voloknistoy poristoy sredy [Slip Condition on the Surface of a Model Fibrous Porous Medium]. Pisma vZhTF [Technical Physics Letters], 2009, vol. 35, iss. 5, pp. 103-110.
6. Mosina E.V. Chislennoe issledovanie techeniya na granitse zhidkost — poristaya sreda [Numerical Study of Flow at a Liquid — Porous Medium Interface]. TOKhT [Theoretical Foundations of Chemical Engineering], 2010, vol. 44, iss. 5, pp. 536-542.
7. Happel J., Brenner H. Gidrodinamika pri malykh chislakh Reynoldsa [Low Reynolds Numbers Hydrodynamics]. Moscow, Mir Publ., 1976. 631 p.
8. Bayta§ A.C., Erdem§ D., Acar H., Cetiner O., Ba§ci H. Analytical Determination of the Permeability for Slow Flow Past Periodic Arrays of Cylinders with Different Cross Sections. J. Porous Med., 2012, vol. 15, pp. 1009-1018. DOI: http://dx.doi.org/10.1615/JPorMedia.v15.i11.20.
9. Drummond J.E., Tahir M.I. Laminar Viscous Flow Through Regular Arrays of Parallel Solid Cylinders. Int. J. Multiphase Flow, 1984, vol. 10, pp. 515-540. DOI: http://dx.doi.org/10.1016/0301-9322(84)90079-X.
10. Hellou M., Martinez J., El Yazidi M. Stokes Flow Through Microstructural Model of Fibrous Media. Mech. Res. Commun., 2004, vol. 31, pp. 97-103. DOI: http://dx.doi.org/10.1016/S0093-6413(03)00081-8.
11. James D.F., Davis A.M.J. Flow at the Interface of a Model Fibrous Porous Medium. J. Fluid Mech, 2001, vol. 426, pp. 47-72. DOI: http://dx.doi.org/10.1017/S0022112000002160.
12. Tahir M.A., Tafreshi H.V. Influence of Fiber Orientation on the Transverse Permeability of Porous Media. Phys. Fluids, 2009, vol. 21, article ID: 083604-1-083604-5. DOI: http://dx.doi.org/10.1063/L3211192.
13. Tamayol A., Bahrami M. Transverse Permeability of Fibrous Porous media. Phys. Review, 2011, vol. 83 (E), article ID: 046314-1-046314-9. DOI: http://dx.doi.org/10.1016/j-.ijheatmasstransfer.2008.09.032.
14. Wang C.Y. Stokes Flow Through an Array of Rectangular Fibers. Int. J. Multiphase Flow, 1996, vol. 22, pp. 185-194. DOI: http://dx.doi.org/10.1016/0301-9322(95)00063-1.
15. Yadav P.K. Slow Motion of a Porous Cylindrical Shell in a Concentric Cylindrical Cavity. Meccanica, 2013, vol. 48, pp. 1607-1622. DOI: http://dx.doi.org/10.1007/s11012-012-9689-0.
THE PERMEABILITY OF TWO-DIMENTIONAL POROUS MEDIUM OF SQUARE FIBERS (CELL MODEL)
Ekaterina Vladimirovna Mosina
Candidate of Physical and Mathematical Sciences, Senior Lecturer, Department of Fundamental Computer Science and Optimal Control, Volgograd State University katefrkate@mail.ru, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Igor Viktorovich Chernyshev
Candidate of Physical and Mathematical Sciences, Associate Professor, Department of Fundamental Computer Science and Optimal Control, Volgograd State University igor_chernyshev@mail.ru, fiou@volsu.ru
Prosp. Universitetsky, 100, 400062 Volgograd, Russian Federation
Abstract. In the present study we have solved numerically the problem of a creeping two-dimensional viscous incompressible flow through a fibrous porous medium. The medium is modeled by a system of square fibers (rods) arranged across the flow. The viscous flow around the fibers is solved under the Stokes approximation using the cell model. The Happel boundary conditions (zero normal
velocity and tangential shear stresses) have been posted on the boundary of the fluid cell. The numerical realization is based on finite differences of the second order accuracy with successive over-relaxation method for solving a system of linear equations.
For regular lattice of square cylinders the hydrodynamic microscopic velocity fields were determined also by the direct numerical simulations. The permeability coefficient of the model porous medium is evaluated by macroscopic averaging and Darcy low assumption. It was found that the cell model with the Happel boundary conditions gives higher values of permeability in comparison with direct numerical calculations.
A comparison between these results obtained by cell model and other direct numerical calculations, empirical and analytical relationships in a wide range of the solid concentration * made it possible to make the following inferences.
The permeability of a rarefied system of square cylinders can be described by an explicit expression that was obtained for circular cylinders [11] but with form factor a2 that reflects the difference between the square and circular cross sections of the rods
k 11/1 \ "2 = ln Z - L476 + 2* ' * ^ 0.2.
a2 a2 8* V * J
For moderately rarefied porous media the semi-empirical Kozeny — Carman formula is appropriate. In this case corrective factor 1/a2 gives error of the permeability more than 25%. To decrease the disagreement of the numerical values with the Kozeny — Carman formula this relation is better to use with the factor 8/9. The specified expression for permeability is
k 1 / 1 1 - *2 \
-7 = - ln---* , 0.1 < * < 0.6.
a2 9* V * 1 + *2 /
At higher fiber volume fraction the Kozeny coefficient is cu = 6.0, which corresponds to a large number of experimental data on the real porous beds consisting of particles of various forms [7]
k (1 - *)3
^ = LM*2L • 06 « * < o-8-
These formulas are good approximations for the permeability coefficient for the system of square rods. Obtained relations for the permeability can be useful in the solutions of the macroscopic filtration problems with fibrous porous media.
Key words: model porous medium, square fibers, Stokes approximation, Happel cell model, permeability.
