Численное решение задач фильтрации неоднородных жидкостей в пористых средах

Размещено на http://www.allbest.ru/

1

МИНИСТЕРСТВО ВЫСШЕГО И СРЕДНЕ-СПЕЦИАЛЬНОГО ОБРАЗОВАНИЯ РЕСПУБЛИКИ УЗБЕКИСТАН

САМАРКАНДСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ им. А. НАВОИ

МЕХАНИКО-МАТЕМАТИЧЕСКИЙ ФАКУЛЬТЕТ

КАФЕДРА «ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ»

УДК

ДИССЕРТАЦИЯ

на соискание степени магистра по специальности «Вычислительная математика»

ЧИСЛЕННОЕ РЕШЕНИЕ ЗАДАЧ ФИЛЬТРАЦИИ НЕОДНОРОДНЫХ ЖИДКОСТЕЙ В ПОРИСТЫХ СРЕДАХ

Специальность: 5А480101 «Вычислительная математика»

ШАИМОВ КОМИЛЖОН МИРЗАКАБУЛОВИЧ

Работа обсуждена и

представлена к защите.

Научный руководитель:

Заведующий кафедрой

д.ф.-м.н., проф. Хужаёров Б.Х.

«Вычислительные методы»

к.ф.-м.н., доц. Абдирашидов А.

Самарканд - 2012

Оглавление

Введение

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

1.1 Особенности процесса осветления суспензий путем фильтрования

1.2 Модифицированные кинетические уравнения осаждения и освобождения твердых частиц суспензии

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

2.1 Решение задачи фильтрации суспензии с заданным скоростным режимом

2.2 Решение задачи фильтрации в конечном фильтре с заданным режимом подавлению

3. Задачи фильтрации суспензии в пористой среде с учетом динамических факторов

3.1 Решение задачи фильтрации суспензии в полубесконечном пласте

3.2 Решение задачи фильтрации в конечном фильтре

Заключения

Литература

Приложение

Введение

Задачи фильтрации жидкости с твердыми взвешенными частицами (суспензии) в пористой среде встречаются при исследовании различных задач техники. Известно, что в природе многие жидкости являются неоднородными, т.е. состоять из различного количества компонентов и фаз, свойства которых могут значительно отличаться друг от друга. В частности суспензии состоят из жидкой и твердой фазы, при этом мелкие частицы твердой фазы находятся во взвешенном состоянии в жидкой фазе. При течении таких суспензий наблюдается ряд явлений, которые не наблюдаются при течении однородных жидкостей. За счет взаимодействия частиц твердой фазы с частицами жидкой фазы суспензии приобретают неньютоновские реологические свойства. При движении таких суспензий в пористой среде помимо отмеченных неньютоновских свойств наблюдается еще ряд аномальных явлений. В частности частицы твердой фазы могут осаждается в порах, что приводит к выключению этой поры из процесса фильтрации. В процессах движения суспензий в пористой среде частицы жидкой и твердой фаз находится под воздействием различных сил. Силы давления стремятся протолкнуть осажденные твердые частицы из одних пор в другие. При этом градиент давления должен быть таким, чтобы преодолеть силы сцепления частиц в порах. Таким образом при фильтрации суспензии с образованием осадка существенно влияют на емкостные и фильтрационные характеристики пористой среды.

Для качественного и количественного исследования процессов фильтрации с образованием осадка должны быть созданы эффективные математические модели. Анализ показывает, что разработанные до настоящего времени модели фильтрации суспензии с образованием осадка в пористой среде в той или иной мере описывают основные характерные особенности процесса фильтрации. Небольшой анализ некоторых моделей фильтрации суспензии с образованием осадка приведен в первом параграфе данной работы. В работе [1] рассмотрено особенности процесса осветления суспензий путем фильтрования. В [1] модель не учитывает диффузионный перенос вещества в пористой среде. Кроме того, кинетика осаждения и отрыва частиц в поровом пространстве определяется от характеристик насыщения порового пространства - концентрациями взвешенных частиц в жидкости и осадка. В [12] предлагаются модифицированные уравнения кинетики процесса осадка частиц, учитывающие динамические факторы. Однако в достаточном объеме задачи фильтрации с образованием осадка для этой модели, ещё не исследовались. Кроме того, численные методы решения этих задач также не отработаны достаточно хорошо. Поскольку модель состоит изсистему нелинейных дифференциальных уравнений, разработка эффективных численных алгоритмов решения задач является важным элементом исследования этой модели. При этом имеется в виду метод конечных разностей, как наиболее универсальный и широко используемый метод решения дифференциальных уравнений. Исходя из сказанного можно заключить, что исследование моделей фильтрации суспензии с образованием осадка с учетом предельного градиента давления и разработка эффективных численных алгоритмов реализации этих моделей является актуальной задачей.

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

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

построение конечно-разностных аппроксимаций поставленных задач и исследование устойчивости разностных задач;

определение алгоритмов решения задач, составление программ расчетов на ЭВМ и получение численных результатов;

анализ результатов и разработка выводов и заключений.

Объектом исследования являются насыщенные однородной и дисперсной жидкостью пористые среды.

Методика исследования.

Для составления моделей фильтрации суспензии с образованием осадка используются классические закон механики, в частности, закон сохранения массы взвешенных твердых частиц в элементарном объеме. При составлении систем уравнений используется также кинетическое уравнение изменения насыщенности пористой среды осевшими частицами, обобщенный закон Дарси. Для решения системы уравнений используется метод конечных разностей с проверкой устойчивости метода. При анализе результатов принимается во внимание физическая сущность изменения расчетных параметров. На полученные новые явления даются соответствующие физические интерпретации.

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

Научная новизна исследований.

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

Работа состоит из введения, три параграфов, выводов и заключений, списка литературы и приложения.

В первом параграфе дается краткий анализ моделей фильтрации суспензии с образованием осадка. Рассматриваются особенности процесса осветления суспензий путем фильтрования. Затем на основе [12] дается модифицированные кинетические уравнения, учитывающие динамические факторы.

В отдельном параграфе решена задача фильтрования суспензии в пористой среде, на основе кинетического уравнения без учета динамических факторов. Затем та же задача решена в конечном фильтре. Задачи решены численно конечно-разностным методом. Определены профили концентраций взвешенных частиц в жидкости, осадке и градиента давлении.

Решение задачи фильтрации суспензии в пористой среде на основе кинетического уравнения с учетом динамических факторов даны в параграфе 3. Определены профили концентраций взвешенных частиц в жидкости, осадке и градиента давлении для различных исходных данных.

На основе полученные результатов разработаны выводы и заключения.

В конце диссертации программы расчетов решений задач приведены в качестве приложения.

По материалам диссертации опубликована следующая работа:

Махмудов Ж.М., Пардаев А.И., Шаимов К.М.Численное решение задачи фильтрации суспензии в пористой среде // Труды научной конференции “Проблемы современной математики”, 22-23.04.2011. г.Карши. С.431-434.

алгоритм задача уравнение фильтрация суспензия

§1. Учет динамических факторов в кинетике осаждения (захвата) и освобождения частиц при фильтрации суспензий

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

1.1 Особенности процесса осветления суспензий путем фильтрования

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

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

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

В [2] уравнение кинетики представлялось в форме

, (1.1)

где с - концентрация осадка, с- концентрация суспензии в жидкой фазе, б, в -постоянные.

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

Для описания динамики осветления суспензии широко применяются модели [3-6], которых в обобщенном виде можно представить как [1]

, (1.2)

где - скорость фильтрации; - скорость прилипания частиц; - скорость отрыва частиц; - потери напора на фильтре; - начальные потери напора; - гидравлический уклон; - длина фильтра.

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

Уравнение изотермы следует из уравнения кинетики системы (1.2) при условии , т.е.

, (1.3)

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

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

В линейной модели фильтрования в (1.2) принимается

, , (1.4)

где - постоянная.

Модель фильтрации состоит из уравнения баланса и кинетики (1.2) с (1.4). В модели предполагается, что каждая частица, осевшая на поверхности пористого слоя, может оторваться независимо от концентрации осадка. Очевидно, что такое предположение справедливо лишь для осадков, образующих на поверхности загрузки монослой из частиц суспензии. Кроме того, полагается, что любая частица имеет равную вероятность отрыва независимо от места адгезии. Это предположение опровергается последними данными о структуре потока в пористом слое.

В нелинейной модели фильтрования [5] полагается, что

, ,(1.5)

где - постоянные. Эта система совпадает с уравнениями динамики сорбции при внешнедиффузионной кинетике и изотерме Ленгмюра. По сравнению с линейной моделью здесь учитывалось влияние формирующегося осадка на изменение локальной скорости, что в свою очередь влияет на скорость отрыва. В рамках нелинейной модели удается доказать формирование концентрационного фронта, перемещающегося по слою с постоянной скоростью [1].

В моделях безотрывного фильтрования полагается, что , а различными авторами аппроксимируется по-разному. Так, в моделях Айвеса [6]

, ,(1.6)

в модели Шехтмана [3]

,(1.7)

где - постоянные (). Уравнения динамики аналогичны (хотя и не тождественны) системе динамики сорбции при внешней диффузии и прямоугольной изотерме [1].

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

В моделях (1.6), (1.7) учитывается формирование в пористом слое невымываемой части осадка, но не учитывается образование вымываемой части, что несомненно доказано экспериментально [7].

Осадок в пористом пространстве глобулярной структуры фильтра имеет две формы - вымываемую и невымываемую. Соответственно зоны фильтра назовем активными и пассивными. Активные зоны, омываемые струйной составляющей потока, формируют вымываемый осадок с концентрацией , пассивные зоны, являющиеся застойными, формируют невымываемый осадок с концентрацией . Обозначим полную емкость фильтра через. Из изложенного следуют

,

где и - емкости соответственно активных и пассивных зон. Подчеркнем, что указанные емкости являются динамическими характеристиками. Они зависят не только от «качества» дисперсной фазы, но и от скорости и структуры потока, а также геометрии слоя [1].

Степень вымываемости осадка характеризует величина

.

При фильтровании суспензий уравнение кинетики необходимо записывать для обеих зон: активной и пассивной. Приведем такой подход в соответствии с [1]. Уравнение кинетики для вымываемой формы осадка



совпадает с линейной моделью фильтрации (1.5). Кинетика образования невымываемой формы осадка описывается уравнением

и аналогична моделям безотрывного фильтрования (1.6). Изменение связано с эффектом уплотнения (старения) осадка.

В [1] описан процесс старения осадка пассивной зоны. В соответствии с этим описанием пока объем пассивной зоны не заполнен, каждая частица, подошедшая к ее границе, переходит в невымываемый осадок. Если же зона заполнена, то вероятность такого перехода меньше единицы. Обозначим вакантную часть поверхности через. Тогда при частично вакантной зоне кинетический коэффициент равен , если для свободной поверхности он был равен . Уменьшение происходит непрерывно, поскольку за счет старения осадка появляются новые вакансии в пассивной зоне. Наиболее точно экспериментальные данные описывает следующая аппроксимация :

(1.8)

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

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

Уравнение баланса

(1.9)

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

, .(1.10)

Решение уравнения получено сначала для малых времен. Основной вывод для этого диапазона времени заключается в том, что на достаточно коротких слоях можно наблюдать активных зон еще до того, когда начнет «ухудшаться» кинетика пассивных зон. Решение получено также при предположении отсутствия активной составляющей скорости формирования осадка.

1.2 Модифицированные кинетические уравнения осаждения и освобождения твердых частиц суспензии

Рассмотренная в §1.1 модель не учитывает диффузионный перенос вещества в пористой среде. Кроме того, кинетика осаждения и отрыва частиц в поровом пространстве определяется от характеристик насыщения порового пространства - концентрациями и . Коэффициенты в кинетическом уравнении , , , а также сами функции , считаются не зависящими от интенсивности гидродинамического процесса фильтрации. Эти предположения, возможно, справедливы, когда скорости фильтрации и, соответственно, градиенты давления не высоки. В [8-11] предлагаются модифицированные уравнения кинетики кольматации и суффозии, учитывающие вероятностный подход. Уравнение кинетики осаждения частиц с учетом кольматации и суффозии записано в виде

,(1.11)

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

Напишем уравнение кинетики [1] с учетом динамических факторов. Сначала кинетическое уравнение представим в виде

.(1.12)

Первый член справа в (1.12) учитывает, как отмечалось в §1.1, эффект «зарядки» фильтра, характеризующий улучшение эффективности массопереноса из-за увеличения удельной поверхности зарядки при образовании осадка.

Второй член справа (1.12) соответствует нелинейной модели фильтрования [5].

Когда в уравнении кинетики (1.12) образование осадка в пассивной зоне принимаем в соответствии с моделями безотрывного фильтрования, получаем следующее

, (1.13)

,(1.14)

, .

Первый член в правой части (1.13) соответствует модели Айвеса [6], а в (1.14) - модели Шехтмана [3].

Кинетические уравнения (1.13), (1.14) необходимо рассматривать совместно с уравнением баланса, которого в одномерном случае с учетом диффузионных эффектов запишем в виде

, (1.15)

где - пористость среды, - коэффициент эффективной диффузии, - скорость фильтрации, const.

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

, , (1.16)

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

,(1.17)

получаем

,(1.18)

где .

Для конечного фильтра длиной граничные условия для (1.18) могут быть заданы в виде

,(1.19)

, - заданные давления.

В соответствии с [8-11] кинетические уравнения (1.12), (1.13), (1.14) с учетом динамических факторов в работе [12] записываются в виде

,(1.20)

, (1.21)

,(1.22)

где , - постоянные коэффициенты.

Уравнения (1.20), (1.21), (1.22) учитывают динамические факторы, что представляют модификации уравнений (1.12), (1.13), (1.14), соответственно.

§ 2. Задачи фильтрации суспензии в пористой среде

В этом параграфе решаются задачи фильтрации суспензии в пористой среде на основе предложенной в 1.2 модели.

2.1 Решение задачи фильтрации суспензии с заданным скоростным режимом

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

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

I вариант. Система уравнений состоит из (1.15) и (1.12). Тогда в одномерном случае она представляется в виде [13]

,

.(2.1)

Начальные и граничные условия имеют вид

(2.2)

Для решения задачи (2.1) - (2.2) применим метод конечных разностей [14] В области вводим сетку, где - максимальное время, в течение которого исследуется процесс Для этого интервал разбиваем с шагом , а разбиваем на J частей с шагом В результате имеем сетку

Вместо функций , будем рассматривать сеточные функции, значения которых в узлах соответственно обозначим ,

Первое уравнение системы (2.1) аппроксимируется на сетке вследующим виде

.(2.3)

Разностная схема для второго уравнения системы (2.1) будет иметь вид:

.(2.4)

Начальные и граничные условия (2.2) также представляются в сеточном виде

, , ,, , ,

, , ,(2.5)

, , ,

где - достаточно большое число, для которого приблизительно выполняется .

Разностная схема (2.3) приведена к системе линейных алгебраических уравнений (СЛАУ)

, , ,(2.6)

где

, , ,

.

Схема (2.4) представляется в виде

,, ,(2.7)

Где

, .

Систему (2.6) решаем методом прогонки при известных, определенных в соответствии (2.7)

,, ,(2.8)

где

, , , .

Из граничного условия имеем

, откуда , .(2.9)

Вычисления проводятся в следующей последовательности. Согласно (2.7) определяются значения через известные величины и нижнего слоя в соответствующих точках, из (2.8) находятся В качестве исходных параметров примем следующие числовые величины: , , м/c, м²/c. Проведем анализ численных результатов. С течением времени значения , в фиксированных точках пласта возрастают (Рис.2.1 - 2.4). При увеличении параметра , графики показывают увеличение значений концентрации осадка и уменьшение концентрации взвешенных твердых частиц в жидкости в соответствующих точках пласта (Рис.2.1, 1 и 2 случаи). Сравнивая рис.2.1 (1 случай) с рис. 2.2 (1 случай) можно сказать, что увеличение параметра также приводит к увеличению концентрации осадка и уменьшению в фиксированных точках пласта. Уменьшение значения параметра приводит к увеличению отрыва частиц, за счет этого уменьшаются абсолютные значения в одинаковых точках пласта. В случае с, в точке было (Рис.2.1б, 1 случай), а при с, имеем (Рис.2.2б, 2 случай). Сравнивая графические результаты рис.2.1 (1 случай) с рис.2.3 (2 случай) можно сказать, что при увеличении параметра уменьшаются значения , но не очень значительно. При увеличении значения в 50 раз, значения понижаются приблизительно на 20 % от предыдущего случая (при). При увеличении коэффициента эффективной диффузии в уравнении баланса профили показателей и распространяются более широким фронтом. Увеличив скорость фильтрации замечаем ускоренное продвижение фронта в глубь пласта (Рис.2.4).

II вариант. Система уравнений состоит из (1.15) и (1.13)

,

.(2.10)

В этом варианте уравнение кинетики системы (2.10) после аппроксимация принимается, как в (2.7), в виде

,, ,(2.11)

где,

.

Результаты численных экспериментов представлены на рис. 2.5-2.6. В этом варианте, как и в предыдущем, с увеличением параметра увеличивается концентрация осадки и уменьшается концентрация взвешенных частиц (Рис. 2.5). С увеличением значений параметров (Рис.2.6, 1 случай) и (Рис. 2.6, 2 случай) также происходит увеличение и уменьшение в фиксированных точках пластах.

III вариант. Система уравнений состоит из (1.15) и (1.14)

,

.(2.12)

В этом варианте кинетическое уравнение аппроксимируется аналогично (2.11), а и в (2.11) выражаются следующими соотношениями

, .(2.13)

Результаты численных экспериментов представлены на рис. 2.7 - 2.8. Сравнивая 1 и 2 случаи рис.2.8 можно сказать, что при увеличении значений параметра скорость прилипания частиц увеличивается, за счет этого уменьшается значение концентрации взвешенных частиц, а с увеличением значения параметра уменьшается , увеличивается и наблюдается продвижение фронта в глубь пласта.

2.2 Решение задачи фильтрации в конечном фильтре с заданным режимом по давлению

Здесь рассматривается задача для конечной области с граничным режимом по давлению

В однородный пласт длины l, пористостью , насыщенный чистой (без частиц) жидкостью с некоторого момента времени начинает поступать дисперсная система с концентрацией взвешенных твердых частиц Течение происходит под действием градиента давления, создаваемого постоянными граничными давлениями и Исходя из сформулированных условий задачи начальные и граничные условия записываются так

, .(2.14)

Таким образом, для определения неизвестных , , необходимо решить систему, состоящую из уравнения баланса, уравнения кинетики, законафильтрации (уравнение (1.16) в 1.2) и уравнения (1.18) в 1.2 Как в параграфе 2.1 рассмотрим следующие варианты, зависящие от вида кинетического уравнения.

I вариант. Система уравнений состоит из (1.15), (1.12), (1.16), (1.18)

,

,(2.15)

, ,

, .

Система решается при условиях (2.14) методом конечных разностей

В области вводим сетку

,

где T - максимальное время, в течение которого исследуется процесс, I, J - количество точек деления сетки.

Последнее уравнение системы (2.15), когда коэффициент фильтрации выражается законом Кармана-Козени (), представим в виде

.(2.16)

Сеточное решение, определенное в точках для давления p обозначим через В результате получается следующая разностная схема

(2.17)

где .

Условия (2.14) для давления принимают вид:

(2.18)

Схема (2.17) также приводится к СЛАУ

, , ,(2.19)

где

, , , .

Решая систему уравнений (2.19) методом прогонки определяется давление на ()-ом слое и в соответствующих точках находится градиент давления по формуле:

(2.20)

Затем по схеме (2.7) вычисляются . Скорость фильтрации определяется как

.(2.21)

Первое уравнение системы (2.15) аппроксимируется на сетке в следующем виде

.(2.22)

Схема (2.22) приводится к СЛАУ и решается методом прогонки при известных значениях и

Результаты численных расчетов представлены на рис.2.9 - 2.11. Из этих рисунков видно, что с увеличением времени в одних и тех же точках пласта значения и возрастают. Анализируя профили отметим, что они также носят монотонно убывающий характер, а затем стабилизируются достигнув определенной величины. С увеличением времени вблизи

Рис. 2.9. Профили (а), (б) и p (в) в различные моменты времени при м/с, м²/с, (1), 0,2 (2) с^-1, с, , .

Рис. 2.10. Профили (а), (б) и p (в) в различные моменты времени при м/с; м²/с; с^-1; ;

с, (1); с, (2).

Рис. 2.11. Профили концентраций (а), (б), (с) в различные моменты времени при м/с; с^-1; с; ; м²/с, (1); м²/с, (2). уменьшается значение при увеличении значения параметров , и . Когда увеличиваются значения параметров , , в рассмотренной задаче в полубесконечной области, значение увеличивается, уменьшается и в этом случае вблизи свободной границы области () наблюдается резкое возрастание абсолютных значений градиента давления (Рис.2.9в,2 случай, рис.2.10в,2 случай).

II вариант. Система уравнений состоит из (1.15), (1.13), (1.16), (1.18)

,

,(2.23)

, ,

, .

Эта система решается при начальных и граничных условиях (2.14). Для определения сеточных значений , , и использованы схемы (2.17), (2.20), (2.11), (2.21) и (2.22).

Результаты численных расчетов представлены графически на рис.2.12, 2.13. Увеличение параметров и , как и в задаче, рассмотренной в полубесконечной области, приводит к увеличению концентрации осадка и это можно рассматривать как положительный фактор при фильтровании суспензий в пористой среде.

Рис. 2.12. Профили (а), (б) и p (в) в различные моменты времени при м²/с, (1), 0,2 (2) с^-1, с, , , .

Рис. 2.13. Профили (а), (б) и p (в) в различные моменты времени при м²/с; с^-1; ; , (1); , (2).

III вариант. Система уравнений состоит из (1.15), (1.14), (1.16), (1.18)

,

,(2.24)

, ,

, .

Эта система также решается при условиях (2.14). В этом варианте , , определяются как в предыдущем варианте, только для вычисления использована схема (2.13).

Результаты численных расчетов представлены на рис.2.14, 2.15.

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

Рис. 2.14. Профили (а), (б) и p (в) в различные моменты времени при м²/с, (1), 0,2 (2) с^-1, с, , , .

Рис. 2.15. Профили (а), (б) и p (в) в различные моменты времени при м²/с; с^-1; ; , (1); , (2).

§ 3. Задачи фильтрации суспензии в пористой среде с учетом динамических факторов

В этом параграфе решаются задачи фильтрования суспензии в пористой среде на основе новых кинетических уравнений, которые учитывают динамические факторы (уравнения (1.20) - (1.22)).

3.1 Решение задачи фильтрации суспензии в полубесконечном пласте

Объектом исследования, как и в параграфе 2.1, является полубесконечный однородный пласт.

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

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

I вариант. Система уравнений состоит из (1.15), (1.20) и (1.16):

,

.(3.1)

, , .

Начальные и граничные условия имеют вид:

(3.2)

Для решения задачи (3.1) - (3.2) применим метод конечных разностей

Первое уравнение системы (3.1) аппроксимируется на сетке как в схеме (2.3), а второе уравнение - в следующем виде

.(3.3)

Эта схема может быть представлена в виде

,, .(3.4)

Здесь

, ,

где

, .

Из третьего уравнения системы (3.1) определим

.(3.5)

Аппроксимируя (3.5) получим

.(3.6)

Вычисления проводятся в следующей последовательности. Согласно (3.6) определяются значения через известные величины нижнего слоя в соответствующих точках, подставляя найденные в (3.4) находятся и согласно (2.3) определяются . По полученным решениям определены профили , и в различные моменты времени (Рис. 3.1 - 3.4). Здесь также увеличение значений параметров (Рис.3.1, 2 случай), (Рис.3.2, 2 случай), как в решенной задаче в §2.1, действует как положительный фактор к прилипанию частиц, т.е. концентрация осадки возрастает. Уменьшение значения параметра (Рис.3.2, 1 случай) и увеличение (Рис.3.3, 2 случай) приводит к возрастанию отрыва осевших частиц. Сравнивая рис.2.1 и рис.3.1 оценить влияние градиента давления на процесс кинетики захвата и отрыва частиц.

Рис. 3.1. Профили (а), (б) и p (в) в различные моменты времени при м/с, м²/с, (1), 0,2 (2) с^-1, с, , , м/МПа, м/МПа.

Рис. 3.2. Профили (а), (б) и p (в) в различные моменты времени при м/с; м²/с; с^-1; ; м/МПа; м/МПа; c, (1);c, (2).

Рис. 3.3. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; ; м²/с; м/МПа; м/МПа, (1); м/МПа, (2).

Рис. 3.4. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; ; ; м/МПа;

 м/МПа, м²/с (1); м/МПа, м²/с (2).

Если при, в точке было (Рис.2.1, 1 случай), то при , имеем (Рис.3.1, 1 случай). Таким образом, увеличение параметров (Рис.3.3, 1 случай), (Рис.3.4, 1 случай) приводит к относительному уменьшению значения и относительному увеличению .

II вариант. Система уравнений состоит из (1.15), (1.21) и (1.16), т.е.

,

,(3.7)

, , .

В этом варианте уравнение кинетики системы (3.7) после аппроксимации приобретает (как в (3.4)) вид

,, .(3.8)

Здесь

, ,

где

, .

Вычисления проводятся, аналогично предыдущему варианту. Только для определения вместо схемы (3.4) используется схема (3.8). Результаты нескольких численных расчетов представлены на рис.3.5 - 3.7. Сравнение рис.3.6 с рис.3.5 (1 случай) показывает, что увеличение значений параметров , приводит к увеличению концентрации и уменьшению . Сравнение графических результатов рис.3.7 (2 случай) с рис.3.5 (1 случай) дает, что увеличение параметров , в 10 раз приводит к резкому

Рис. 3.5. Профили (а), (б) и p (в) в различные моменты времени при м/с, м²/с, (1), 0,2 (2) с^-1, с, , , , , , м/МПа, м/МПа.

Рис. 3.6. Профили (а), (б) и p (в) в различные моменты времени при м/с; м²/с; с^-1; ; c; ; ;

м/МПа; м/МПа;

, (1);, (2).

Рис. 3.7. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; м²/с; ; , , ; ;

м/МПа, м/МПа (1); м/МПа, м/МПа (2).

относительному уменьшению значений , и увеличению в фиксированных точках пласта.

III вариант. Система уравнений состоит из (1.15), (1.22) и (1.16), т.е.

,

,(3.9)

, , .

Второе уравнение этой системы аппроксимируется аналогично (3.8), а и выражаются следующими соотношениями

, ,(3.10)

где

, .

Результаты этих расчетов представлены на рис.3.8 - 3.10. Сравнение рис.3.8 с рис.3.9 (1 случай) показывает, что увеличение значения параметра приводит к относительному увеличению концентрации и уменьшению , а увеличение (Рис.3.9, 2 случай) - к уменьшению и увеличению в соответствующих точках пласта. Увеличение параметров и (Рис.3.10), также приводит к возрастанию отрыва осевших частиц.

Рис. 3.8. Профили (а), (б) и p (в) в различные моменты времени при

м/с, м²/с, (1), 0,2 (2) с^-1, с, , ,

, м/МПа, м/МПа.

Рис. 3.9. Профили (а), (б) и p (в) в различные моменты времени при м/с; м²/с; с^-1; ; c;

м/МПа; м/МПа;

, (1);, (2).

Рис. 3.10. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; м²/с; ; ; ;

м/МПа, м/МПа (1); м/МПа, м/МПа (2).

3.2 Решение задачи фильтрации в конечном фильтре

Здесь рассматривается задача для конечной области с граничным режимом по давлению, как в 2.2 Аналогично рассмотрим следующие варианты, зависящие от вида кинетического уравнения.

I вариант. Система уравнений состоит из (1.15), (1.20), (1.16), (1.18)

, ,(3.11)

, , ,

, .

Начальные и граничные условия записываются как (2.14):

, .(3.12)

Первое уравнение системы (3.11) аппроксимируется как (2.22), второе - (3.4), 3-ое - (2.21) и 4-ое уравнение как (2.17).

Последовательность расчета как в параграфе 2.2. Сначала решая систему (2.19) определяются значения градиента давления, затем по схеме (3.4) вычисляется , а также скорость фильтрации по (2.21). Подставляя найденные и в (2.22) находятся .

II вариант. Система уравнений состоит из (1.15), (1.21), (1.16), (1.18)

,

,(3.13)

, , ,

, .

Система решается методом конечных разностей при условиях (3.12). Последовательность расчета как в первом варианте, только при определения используется схема (3.8).

III вариант. Система уравнений состоит из (1.15), (1.22), (1.16), (1.18)

,

,(3.14)

, , ,

, .

Здесь при определения используется схема (3.10).

Результаты численных расчетов представлены на рис.3.11 - 3.14 для I варианта, на рис.3.15 - 3.17 для II варианта, на рис.3.18 - 3.20 для III варианта.

В отличие от результатов, которые приведены в параграфе 2.2, здесь исследуется роль параметров и . Результаты показывают, что увеличение параметров и приводит к относительному уменьшению и увеличению .

Рис. 3.12. Профили (а), (б) и p (в) в различные моменты времени при м/с; м²/с; с^-1; ; м/МПа; м/МПа;

c, (1);c, (2).

Рис. 3.13. Профили (а), (б) и p (в) в различные моменты времени при

м/с; с^-1; с; ; м²/с; м/МПа;

м/МПа, (1); м/МПа, (2).

Рис. 3.14. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; ; ; м/МПа; м/МПа, м²/с (1); м/МПа, м²/с (2).

Рис. 3.15. Профили (а), (б) и p (в) в различные моменты времени при м/с, м²/с, (1), 0,2 (2) с^-1, с, , , , , , м/МПа, м/МПа.

Рис. 3.16. Профили (а), (б) и p (в) в различные моменты времени при

м/с; м²/с; с^-1; ; c; ; ;

м/МПа; м/МПа;

, (1);, (2).

Рис. 3.17. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; м²/с; ; ; ;

м/МПа, м/МПа (1); м/МПа, м/МПа (2).

Рис. 3.18. Профили (а), (б) и p (в) в различные моменты времени при

м/с, м²/с, (1), 0,2 (2) с^-1, с, , ,

, м/МПа, м/МПа.

Рис. 3.19. Профили (а), (б) и p (в) в различные моменты времени при

м/с; м²/с; с^-1; ; c;

м/МПа; м/МПа;

, (1);, (2).

Рис. 3.20. Профили (а), (б) и p (в) в различные моменты времени при м/с; с^-1; с; м²/с; ; ; ;

м/МПа, м/МПа (1); м/МПа, м/МПа (2).

Заключение

1. Поставлены начально-краевых задач фильтрации суспензии с нового кинетического уравнения, учитывающие динамических факторов для различных режимов течения.

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

3. Решены задачи переноса вещества в пористой среде с использованием кинетического уравнения, учитывающего динамические факторы, в полубесконечной и конечной пористой среде с заданными режимами по скорости фильтрации и давлений на концах фильтра. Установлено, что в рамках принятых феноменологических допущений учет динамических факторов приводит к усилению суффозионных и ослаблению кольматационных эффектов. Это, в свою очередь, изменяет характер распределения концентрации взвешенных и осевших твердых частиц - концентрационные профили твердых взвешенных и осевших частиц продвигаются более интенсивно в направлении переноса, значения концентрации осевших твердых частиц относительно уменьшаются.

Литература

1. Веницианов Е.В., Рубинштейн Р.Н. Динамика сорбции из жидких сред. - М.: Наука, 1983. - 237 с.

2. Радушкевич Л.В. Природа вторичных процессов при фильтрации аэрозолей. Сообщ. 1. Простейшая теория соосаждения частиц в фильтре и кинетика кольматации // Изд. АН СССР. Сер.хим., 1963, № 3. С. 407-414.

3. ШехтманЮ? М? Фильтрация малоконцентрированных суспензий? - М?: Изд-во АН СССР? 1961? - 212 с?(13)

4. Минц Д.М. Кинетика фильтрации малоконцентрированных водных суспензий на водоочистных фильтрах // ДАН СССР, 1951. Т. 78, № 2. С. 315-318.

5. Веницианов Е.В., Рубинштейн Р.Н., Сенявин М.М. О возможности распространения теории динамики сорбции на расчет осветления воды зернистыми фильтрами // ДАН СССР, 1970. Т. 195, № 3. С. 658-661.

6. IvesK.J. Simulationoffiltrationonelectronicdigitalcomputer // J. Amer. Water Works Assoc., 1960. V. 52, N. 7. Pp. 933-939.

7. Mints D.M. Modern theory of filtration. - In: Inter. water supply congr. Barselona, 1966.Sp. sub.N 10. 32 p.

8. Хужаёров Б? Х? Модель фильтрации суспензии с учетом кольматации и суффозии // ИФЖ. 1992. Т. 63, № 1? С? 72-79?

9. Зубаров А.Ю., Хужаёров Б.Х., К теории релаксационной фильтрации // ИФЖ. 1988. Т. 55. № 3. С. 442-447.

10. Хужаёров Б?Х?, Давиденко М?А? О моделях переноса твердых частиц в пористой среде // «Актуальные проблемы механики контактного взаимодействия»?Сб?тр?респ?науч?конф? Самарканд, 1997? С? 167-169?

11. Хужаёров Б?Х?, Давиденко М?А? О модификации модели кольматационно-суффозионной фильтрации // ДАН РУз, 1997. № 5? С? 25-28?

12. Хужаёров Б.Х., Махмудов Ж.М. Математическое моделирование процессов фильтрации суспензии в пористой среде с учетом образования осадок // Труды научной конференции “Проблемы современной математики”, 22-23.04.2011. г.Карши. С.543-547.

13. Махмудов Ж.М., Пардаев А.И., Шаимов К.М.Численное решение задачи фильтрации суспензии в пористой среде // Труды научной конференции “Проблемы современной математики”, 22-23.04.2011. г.Карши. С.431-434.

14. Самарский А.А. Теория разностных схем. М.: Наука, 1977. - 656 с.

ПРИЛОЖЕНИЕ

program Concentration;

{$APPTYPE CONSOLE}

uses

SysUtils;

Const

al=1.0{50.0}; Gam=5.0{2.5}; K=1.0{50.0};

betta=0.1{0.2};

N=80; m0=0.2;

D=2.5e-6{7.5e-6}; v=1e-4{2e-4}; c00=0.01;

tmax=2700;

h=0.05;

tau=1;

Type massiv=array [0..N] of real;

Var ro,ro0,c,c0:massiv;

alpha,bet:array[0..N] of real;

i,j:integer;

out:text;

Procedure Nach_Dan;{Boshlang'ichvachegaraviyshartlar}

begin

for i:=0 to N do

begin

ro0[i]:=0;

c0[i]:=0;

end;

c0[0]:=c00;

end;

Procedure Progonka; {x yo`nalishb`yicha D1 soxada}

varA,B,E,F:real;

Begin

A:=tau*D/(h*h)+tau*v/h;

B:=2*tau*D/(h*h)+tau*v/h+m0;

E:=tau*D/(h*h);

alpha[1]:=0;

bet[1]:=c00;

for i:=1 to N-1 do

begin

F:=m0*c0[i]-(ro[i]-ro0[i]);

bet[i+1]:=(F+A*bet[i])/(B-A*alpha[i]);

alpha[i+1]:=E/(B-A*alpha[i]);

end;

c[N]:=(F+A*bet[N])/(B-A*alpha[N]);

for i:=N-1 downto 0 do

c[i]:=alpha[i+1]*c[i+1]+bet[i+1];

END;

Procedure Otvet1(k:integer);

Begin

writeln(out);

write (out,' ':10,' c t= ',k*tau:4);

writeln(out);

For i:=0 to N do

writeln (out,i*h:4:2,' ',c[i]/c00:4:12,' ',1e2*ro[i]:4:12);

End;

Procedure Asosiy;

Varkk:real;

BEGIN

Nach_dan;

j:=0;

for i:=0 to N do

begin

ro[i]:=ro0[i];

c[i]:=c0[i];

end;

for j:=1 to tmax do

Begin

for i:=0 to N do

begin

kk:=Gam*(1-k*ro0[i])-betta*al*c0[i]*tau*Gam*(1-k*ro0[i])+tau;

ro[i]:=Gam*(1-k*ro0[i])/kk*ro0[i]+tau*betta*c0[i]*Gam*(1-k*ro0[i])/kk;

end;

Progonka;

for i:=0 to N do

begin

c0[i]:=c[i];

ro0[i]:=ro[i];

end;

if j mod 900=0 then

Otvet1(j);

end;

END;

begin

assign(out,'ORG.out');

Rewrite(out);

writeln (out,' Г=',Gam:4:2,' K=',k:4:2,' Betta=',Betta:4:2,' Alfa=',al:4:2);

writeln (out,' D=',D:4:6,' V=',V:4:6);

Asosiy;

Close(Out);

end.

program Concentration;

{$APPTYPE CONSOLE}

uses

SysUtils;

Const

al=1.0 {50.0}; Gam=5.0{2.5}; K=1.0{50.0};

betta=0.1{0.2};

gamma=0.1{1.0}; omega=0.1{1.0};

N=80; m0=0.2;

D=2.5e-6{7.5e-6}; v=1e-4{2e-4}; c00=0.01;

tmax=2700;

h=0.05;

tau=1;

k0=8e-3;

Type massiv=array [0..N] of real;

Var N_p,ro,ro0,c,c0:massiv;

alpha,bet:array[0..N] of real;

i,j:integer;

out:text;

Np:real;

Procedure Nach_Dan;{Boshlang'ichvachegaraviyshartlar}

begin

for i:=0 to N do

begin

ro0[i]:=0;

c0[i]:=0;

end;

c0[0]:=c00;

end;

Procedure Progonka; {x yo`nalishb`yicha D1 soxada}

varA,B,E,F:real;

Begin

A:=tau*D/(h*h)+tau*v/h;

B:=2*tau*D/(h*h)+tau*v/h+m0;

E:=tau*D/(h*h);

alpha[1]:=0;

bet[1]:=c00;

for i:=1 to N-1 do

begin

F:=m0*c0[i]-(ro[i]-ro0[i]);

bet[i+1]:=(F+A*bet[i])/(B-A*alpha[i]);

alpha[i+1]:=E/(B-A*alpha[i]);

end;

c[N]:=(F+A*bet[N])/(B-A*alpha[N]);

for i:=N-1 downto 0 do

c[i]:=alpha[i+1]*c[i+1]+bet[i+1];

END;

Procedure Otvet1(k:integer);

Begin

writeln(out);

write (out,' ':10,' c t= ',k*tau:4);

writeln(out);

For i:=0 to N do

writeln (out,i*h:4:2,' ',c[i]/c00:4:12,' ',1e2*ro[i]:4:12,' ',N_p[i]:4:12);

End;

Procedure Asosiy;

Varkk:real;

BEGIN

Nach_dan;

j:=0;

for i:=0 to N do

begin

ro[i]:=ro0[i];

c[i]:=c0[i];

end;

for j:=1 to tmax do

Begin

for i:=0 to N do

begin

Np:=v*sqr(1-m0+ro0[i])/(k0*sqr(m0-ro0[i])*(m0-ro0[i]));

kk:=(1-k*ro0[i])*(1+gamma*Np)-tau*betta*al*c0[i]*(1-k*ro0[i])+(1+omega*Np)/Gam*tau*(1+gamma*Np);

ro[i]:=(1+gamma*Np)*(1-K*ro0[i])/kk*ro0[i]+tau*betta*c0[i]*(1-k*ro0[i])/kk;

end;

Progonka;

for i:=0 to N do

begin

c0[i]:=c[i];

ro0[i]:=ro[i];

end;

if j mod 900=0 then

begin

for i:=0 to N do

N_p[i]:=v*sqr(1-m0+ro0[i])/(k0*sqr(m0-ro0[i])*(m0-ro0[i]));

Otvet1(j);

end;

end;

END;

begin

assign(out,'ORG.out');

Rewrite(out);

Asosiy;

Close(Out);

end.

program Concentration;

{$APPTYPE CONSOLE}

uses

SysUtils;

Const

K=1.0{50.0};

betta=0.1{0.2};

a5=0.5{1.0}; a6=1.0{100.0};

Gam=5.0{2.5};

N=80;

m0=0.2;

D=2.5e-6; v=1e-4;

c00=0.01;

tmax=2700;

h=0.05;

tau=1;

k0=8e-3;

gamma=0.1{1.0}; omega=0.1{1.0};

pn=4; pk=0;

Type massiv=array [0..N] of real;

Var W,N_p,p,ro,ro0,c,c0:massiv;

alpha,bet:array[0..N] of real;

i,j:integer;

out,tez:text;

Procedure Nach_Dan;{Boshlang'ichvachegaraviyshartlar}

begin

for i:=0 to N do

begin

ro0[i]:=0;

c0[i]:=0;

end;

c0[0]:=c00;

end;

Procedure Progonka1; {x yo`nalishb`yicha D1 soxada}

varA,B,E,F,q:real;

Begin

alpha[1]:=0;

bet[1]:=pn;

for i:=1 to N-1 do

begin

q:=(3-m0+ro0[i])/((m0-ro0[i])*(1-m0+ro0[i]))*(ro0[i]-ro0[i-1])/h;

A:=1+h/2*q;

B:=2;

E:=1-h/2*q;

F:=0;

bet[i+1]:=(F+A*bet[i])/(B-A*alpha[i]);

alpha[i+1]:=E/(B-A*alpha[i]);

end;

p[N]:=pk;

for i:=N-1 downto 0 do

p[i]:=alpha[i+1]*p[i+1]+bet[i+1];

END;

Procedure Progonka2; {x yo`nalishb`yicha D1 soxada}

varA,B,E,F,v:real;

Begin

alpha[1]:=0;

bet[1]:=c00;

for i:=1 to N-1 do

begin

v:=k0*(m0-ro0[i])*(m0-ro0[i])*(m0-ro0[i])/sqr(1-(m0-ro0[i]))*abs(p[i+1]-p[i])/h;

A:=tau*D/(h*h)+tau*v/h;

B:=2*tau*D/(h*h)+tau*v/h+m0;

E:=tau*D/(h*h);

F:=m0*c0[i]-(ro[i]-ro0[i]);

bet[i+1]:=(F+A*bet[i])/(B-A*alpha[i]);

alpha[i+1]:=E/(B-A*alpha[i]);

end;

c[N]:=(F+A*bet[N])/(B-A*alpha[N]);

for i:=N-1 downto 0 do

c[i]:=alpha[i+1]*c[i+1]+bet[i+1];

END;

Procedure Otvet(k:integer);

Begin

writeln (tez,k:4,' ',1e3*W[i]:4:12);

End;

Procedure Otvet1(k:integer);

Begin

writeln(out);

write (out,' ':10,' c t= ',k*tau:4);

writeln(out);

For i:=0 to N do

writeln (out,i*h:4:2,' ',c[i]/c00:4:12,' ',1e2*ro[i]:4:12,' ',N_p[i]:4:12,' ',1e3*W[i]:4:12);

End;

Procedure Asosiy;

Var kk,Gam1,be:real;

BEGIN

Nach_dan;

j:=0;

for i:=0 to N do

begin

ro[i]:=ro0[i];

c[i]:=c0[i];

end;

for j:=1 to tmax do

Begin

Progonka1;

for i:=0 to N do

begin

N_p[i]:=abs(p[i+1]-p[i])/h;

be:=betta/(1+Gamma*abs(N_p[i]));

Gam1:=Gam/(1+Omega*abs(N_p[i]));

kk:=(Gam1*(1-K*ro0[i])*(1+a6*betta*c0[i]*tau)+tau)/(Gam1*(1-K*ro0[i]));

ro[i]:=ro0[i]/kk+(tau*betta*a5*c0[i])/kk;

end;

Progonka2;

for i:=0 to N do

begin

c0[i]:=c[i];

ro0[i]:=ro[i];

end;

if j<50 then

begin

if (j mod 5=0) or (1 mod j=0) then

begin

i:=35;

W[i]:=abs(k0*(m0-ro0[i])*(m0-ro0[i])*(m0-ro0[i])/sqr(1-(m0-ro0[i]))*abs(p[i+1]-p[i])/h);

Otvet(j);

end;

end

else

if j mod 50=0 then

begin

i:=35;

W[i]:=abs(k0*(m0-ro0[i])*(m0-ro0[i])*(m0-ro0[i])/sqr(1-(m0-ro0[i]))*abs(p[i+1]-p[i])/h);

Otvet(j);

end;

if j mod 900=0 then

begin

for i:=0 to N do

w[i]:=N_p[i]*k0*(sqr(m0-ro0[i])*(m0-ro0[i]))/sqr(1-m0+ro0[i]);

Otvet1(j);

end;

end;

END;

begin

assign(out,'ORG.out');

assign(tez,'ORG.tez');

Rewrite(out);

Rewrite(tez);

writeln (out,' Г=',Gam:4:2,' K=',k:4:2,' Betta=',Betta:4:2);

writeln (out,' D=',D:4:8);

Asosiy;

Close(Out);

Close(tez);

end.

Размещено на Allbest.ru