Все выпуски

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

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

Предлагается алгоритм получения решения уравнений в частных производных с правой частью, заданной на сетке $\{ (x_{1})_{\mu}, (x_{2})_{\mu}, \ldots, (x_{n})_{\mu}\},$ $(\mu=1,2,\ldots,N)\colon f_{\mu}=f((x_{1})_{\mu}, (x_{2})_{\mu}, \ldots, (x_{n})_{\mu}).$ Здесь $n$ — число независимых переменных в исходном уравнении в частных производных, $N$ — число строк в сетке для правой части, $f=f( x_{1}, x_{2}, \ldots, x_{n})$ — правая часть исходного уравнения. Алгоритм реализует редукцию исходного уравнения к системе обыкновенных дифференциальных уравнений (системе ОДУ) с начальными условиями в каждой точке сетки и включает следующую последовательность действий. Ищется решение исходного уравнения, зависящее от одной независимой переменной. Исходному уравнению ставится в соответствие некоторая система соотношений, содержащая произвольные функции и включающая уравнение в частных производных первого порядка. Для уравнения первого порядка выписывается расширенная система уравнений характеристик. Присоединяя к ней остальные соотношения, содержащие произвольные функции, и требуя, чтобы эти соотношения были первыми интегралами расширенной системы уравнений характеристик, приходим к искомой системе ОДУ, завершая редукцию. Предлагаемый алгоритм позволяет в каждой точке сетки находить решение исходного уравнения в частных производных, удовлетворяющее заданным начальным и краевым условиям. Алгоритм применяется для получения решений уравнения Пуассона и уравнения нестационарной осесимметричной фильтрации в точках сетки, на которой заданы правые части соответствующих уравнений.

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

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

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

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

В работе проводился расчет генерации шума вентилятора турбореактивного двухконтурного авиационного двигателя (ТРДД) для различных режимов его работы с помощью собственного программного пакета GHOST CFD, реализованного для графических процессоров (ГПУ). Программный пакет основан на схемах типа DRP (Dispersion Relation Preserving), имеющих высокий порядок аппроксимации и высокую разрешающую способность. Для интегрирования по времени также использовалась оптимизированная схема типа LDDRK (Low Dispersion and Dissipation Runge-Kutta). Для моделирования турбулентности использовался неявный метод крупных вихрей с релаксационной фильтрацией (LES-RF). В качестве ротор-статор-интерфейса применялись пересекающиеся (CHIMERA) сетки. Ускорение за счет использования ГПУ, по сравнению с обычным центральным процессором, составило до порядка 12-20 раз, при этом было достигнуто приемлемое время счета. Расчеты в GHOST CFD проводились в постановке «вентилятор - спрямляющий аппарат наружного контура (СА) с полными колесами лопаток». Результаты расчетов сравнивались как с экспериментальными данными, так и с результатами аналогичных расчетов в коммерческом программном пакете ANSYS CFX. При этом в части расчетов в ANSYS CFX учитывался и направляющий аппарат внутреннего контура (НА).

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