Численное решение задачи переноса вещества в двухзонной среде с нелинейной кинетикой | Статья в журнале «Молодой ученый»

Отправьте статью сегодня! Журнал выйдет 2 ноября, печатный экземпляр отправим 6 ноября.

Опубликовать статью в журнале

Авторы: ,

Рубрика: Математика

Опубликовано в Молодой учёный №18 (204) май 2018 г.

Дата публикации: 03.05.2018

Статья просмотрена: 27 раз

Библиографическое описание:

Джиянов Т. О., Хотамов О. Ю. Численное решение задачи переноса вещества в двухзонной среде с нелинейной кинетикой // Молодой ученый. — 2018. — №18. — С. 4-7. — URL https://moluch.ru/archive/204/49790/ (дата обращения: 22.10.2019).



В [1] рассмотрен перенос вещества в среде с двойной пористостью, учитывающий обратимое и необратимое осаждение вещества в обеих зонах и равновесный обмен первого порядка между зонами. В каждой зоне, т. е. в трещинах и пористых блоках, происходит обратимое и необратимое осаждение (отложение) вещества с различными характеристиками, описываемые линейными уравнениями. Получено аналитическое решение задачи, которое использовано для описания результатов ранее проведенных экспериментов [2].

Используем аналогичную с [3] схему среды с двойной пористостью (Рис.1). В биконтинуальных средах, таких как трещиновато-пористые среды, необходимо использовать двухзонный подход, где в обеих зонах жидкость считается подвижной. В таких средах с двойной пористостью или двойной проницаемостью перенос вещества, как и движение жидкости происходит с различной интенсивностью, порою контрастной. Заметим что такой подход используется и для макроскопически неоднородных сред, где в обеих зонах может происходить конвективный перенос вещества.

Первая зона с индексом 1 в обозначениях имеет высокую проницаемость, а вторая зона — низкую. В каждой зоне есть по два участка, в каждой из которых происходит осаждение вещества с необратимой неравновесной кинетикой.

Рис.1. Схема переноса вещества в двухзонной среде

Здесь в отличие от [4] рассмотрим нелинейные уравнения кинетики.

Уравнения переноса вещества в одномерном случае записываем в виде [4]

,(1)

,

где — время, с, x-расстояние, м,- коэффициент продольный дисперсии, — скорости движения жидкости, м/с,,- объемная концентрация вещества в жидкости,и -концентрации отложенного вещества,, - пористости зон, , - плотность среды, ,- коэффициент массообмена между зонами, .

Осаждение вещества в каждом из участков зон происходит обратимо в соответствии с кинетическими уравнениями

,(2)

,(3)

где ,- коэффициенты отложения вещества от жидкой фазы lна твердую фазу, , - коэффициенты отрыва вещества из твердой фазы и перехода в жидкость, ,

Пусть в первоначально насыщенную чистой (без вещества) жидкостью среду с начального момента времени закачивается жидкость с постоянной концентрацией вещества . Рассмотрим такие периоды времени, где концентрационное поле не достигает правой границы среды, . При отмеченных допущениях начальные и граничные условия для задачи имеют вид

(4)

(5)

(6)

Задача (1) — (6) хотя и является линейной, получение аналитического решения является сложным, т. к. необходимо найти одновременно три поля в каждой из зон. Поэтому для решения задачи применяем метод конечных разностей. В рассматриваемой области введена равномерная по направлениям сетка

,

где I — достаточно большое целое число, выбираемое так, чтобы отрезок [],, перекрывал область расчетного изменения полей Ci, Sai и Ssi, h — шаг сетки по направлению х.

В открытой сеточной области

уравнения (1), (2), (3) аппроксимировались следующим образом

(5)

,

,(6)

,(7)

где ,, — сеточные значения функций , , ,в точке .

Из явных сеточных уравнений (6), (7) определяем ,

,(8)

,(9)

где

, ,

, .

Сеточные уравнения (5) приводятся к виду

,(10)

где , , ,

.

Устанавливается следующий порядок расчета решений. По (8), (9) определяются ,, затем решая систему линейных уравнений (10) методом прогонки — Поскольку, схемы (8), (9) устойчивы, а для (10) условия устойчивости метода прогонки выполняются.

В расчетах использованы следующие значения исходных параметров: , ,,,

дисперсность .

Некоторые характерные результаты показаны на рис. 2–4. Из рис.3 видно, что уменьшение показателя от единицы приводит к замедлению развития профилей концентрации (при неизменных остальных параметрах). Концентрации осажденного вещества при этом имеют опережающее развитие (рис.3,4). Другим словами, уменьшение показателя при прочих неизменных значениях остальных параметров приводит к интенсификации осаждения вещества в обоих участках зон. Как следствие этого в распределении концентраций вещества в подвижной жидкости обеих зон происходит отставание.

Рис.2. Профили концентраций Сl приt=3600 c,

Рис. 3. Профили концентраций Sal при t=3600 c

Рис. 4. Профили концентраций Ssl приt=3600 c

Литература:

  1. Leij F. J., Bradford S. A. Combined physical and chemical nonequilibrium transport model: analytical solution, moments, and application to colloids // Journal of Contaminant Hydrology. 110. 2009. Pp. 87–99.
  2. Bradford S. A., Simunek J., Bettahar M., van Genuchten M. T., Yates S. R. Modeling colloid attachment, straining, and exclusion in saturated porous media // Environmental Science & Technology. 37. 2003. Pp. 2242–2250.
  3. Leij F. L., Bradford S. A. Colloid transport in dual-permeability media // Journal of Contaminant Hydrology. 150. 2013. Pp. 65–76.
  4. Feike J. L., Bradford S. A. Colloid transport in dual-permeability media // Journal of Contaminant Hydrology. 150. 2013.Pp.65−76.
Основные термины (генерируются автоматически): зона, двойная пористость, профиль концентраций, твердая фаза, осаждение вещества, перенос вещества, уменьшение показателя, жидкость, уравнение, среда.


Похожие статьи

Перенос вещества в неоднородной пористой среде с учетом...

В [1] рассмотрен перенос вещества в среде с двойной пористостью, учитывающий обратимое и необратимое осаждение вещества в обеих зонах и равновесный обмен первого порядка между зонами.

Моделирование процессов фильтрации суспензии в пористой...

Перенос вещества в неоднородной пористой среде с учетом...

время релаксации, скорость фильтрации, кубическое уравнение, неравновесная жидкость, пористая среда, вязкость жидкости, характер сопротивления, вид, жидкость.

Модель с распределенными параметрами для описания динамики...

Начнем со статической задачи — движение отсутствует, рассматриваемое веществораствор соли или двухфазная смесь (раствор

1. Исходное утверждение — состояние среды двухфазное, т. е. смесь жидкого раствора при предельной концентрации и твердой фазы...

Закономерности формирования структуры пеноминеральных...

определенной зависимости профиля каналов пены и градиента давления дисперсной фазы.

Влияние биологически активных веществ на физико-химические... Образованная система будет состоять из воды (дисперсионная среда) и твердых частиц глины (дисперсная фаза).

К вопросу определения скорости фильтрации и времени...

Перенос вещества в неоднородной пористой среде с учетом...

Жидкости и пористая среда считаются несжимаемыми. где — пористость среды, — водонасыщенность, — суммарная скорость фильтрации воды и нефти, — общая подвижность, — давление в пласте, — сумма...

Диффузия органических растворителей в полимеры

Соотношение между потоком диффундирующего вещества J и градиентом концентрации в изотропной среде описывается дифференциальным уравнением Фика

Диффузия в металлах и сплавах в твердой фазе. Н. Бромберг. Журнал коллоидной химии.

Определение активности компонентов в биметаллическом расплаве

уравнение давления насыщенного пара: (1). , ; уравнение плотности кипящей жидкости. Единственно то, что ее применение ограничивает

При ректификации нефтегазоконденсатных смесей скорость переноса вещества из одной фазы в другую dM пропорциональна движущей...

Исследование размерных характеристик нанопорошков Со и Fe

В настоящее время активно исследуются методы синтеза веществ в наноразмерном

Все эти факторы приводят к уменьшению удельной поверхности образцов и росту значения Dагр.

При заданной температуре концентрацию можно заменить парциальным давлением .

Похожие статьи

Перенос вещества в неоднородной пористой среде с учетом...

В [1] рассмотрен перенос вещества в среде с двойной пористостью, учитывающий обратимое и необратимое осаждение вещества в обеих зонах и равновесный обмен первого порядка между зонами.

Моделирование процессов фильтрации суспензии в пористой...

Перенос вещества в неоднородной пористой среде с учетом...

время релаксации, скорость фильтрации, кубическое уравнение, неравновесная жидкость, пористая среда, вязкость жидкости, характер сопротивления, вид, жидкость.

Модель с распределенными параметрами для описания динамики...

Начнем со статической задачи — движение отсутствует, рассматриваемое веществораствор соли или двухфазная смесь (раствор

1. Исходное утверждение — состояние среды двухфазное, т. е. смесь жидкого раствора при предельной концентрации и твердой фазы...

Закономерности формирования структуры пеноминеральных...

определенной зависимости профиля каналов пены и градиента давления дисперсной фазы.

Влияние биологически активных веществ на физико-химические... Образованная система будет состоять из воды (дисперсионная среда) и твердых частиц глины (дисперсная фаза).

К вопросу определения скорости фильтрации и времени...

Перенос вещества в неоднородной пористой среде с учетом...

Жидкости и пористая среда считаются несжимаемыми. где — пористость среды, — водонасыщенность, — суммарная скорость фильтрации воды и нефти, — общая подвижность, — давление в пласте, — сумма...

Диффузия органических растворителей в полимеры

Соотношение между потоком диффундирующего вещества J и градиентом концентрации в изотропной среде описывается дифференциальным уравнением Фика

Диффузия в металлах и сплавах в твердой фазе. Н. Бромберг. Журнал коллоидной химии.

Определение активности компонентов в биметаллическом расплаве

уравнение давления насыщенного пара: (1). , ; уравнение плотности кипящей жидкости. Единственно то, что ее применение ограничивает

При ректификации нефтегазоконденсатных смесей скорость переноса вещества из одной фазы в другую dM пропорциональна движущей...

Исследование размерных характеристик нанопорошков Со и Fe

В настоящее время активно исследуются методы синтеза веществ в наноразмерном

Все эти факторы приводят к уменьшению удельной поверхности образцов и росту значения Dагр.

При заданной температуре концентрацию можно заменить парциальным давлением .

Задать вопрос