Научная статья на тему 'Однотемпературная неизотермическая фильтрация двухфазной жидкости в трещиновато-пористых средах'

Однотемпературная неизотермическая фильтрация двухфазной жидкости в трещиновато-пористых средах Текст научной статьи по специальности «Математика»

CC BY
17
4
i Надоели баннеры? Вы всегда можете отключить рекламу.
i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Однотемпературная неизотермическая фильтрация двухфазной жидкости в трещиновато-пористых средах»

32

Методы вычислительной алгебры и решения уравнений математической физики

Список литературы

1. Il'in V. P., Zadorin A. I. Adaptive formulas of numerical differentiation of functions with large gradients // J. of Physics: Conference Series. 2019. V. 1260, 042003.

2. Задорин А. И. Анализ формул численного дифференцирования на сетке Шишкина при наличии пограничного слоя // СибЖВМ. 2018. Т. 21, № 3. С. 243-254.

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

М. И. Иванов, И. А. Кремер, Ю. М. Лаевский

Институт вычислительной математики и математической геофизики СО РАН

Email: kremer@sscc.ru

DOI: 10.24412/cl-35065-2022-1-00-24

В работе представлена однотемпературная модель неизотермической фильтрации двухфазной жидкости в трещиновато-пористых средах. Такая модель включает двойной набор уравнений в порах и трещинах для суммарных скоростей и давлений в смешанном виде, а также переноса фаз жидкостей [1]. По аналогии с [2], в систему уравнений добавлен закон сохранения энергии, который записывается в терминах температуры и теплового потока, а зависимость динамической вязкости фаз жидкости от температуры задается формулой Вальтера. Для численного решения поставленной начально-краевой задачи используется расширение IMPES-схемы, суть которой заключается в комбинировании явно-неявных методов. На каждом шаге времени явным образом пересчитываются нелинейные зависимости и конвективные переносы масс и энергии, а для описания диффузионных изменений давлений и температуры используются неявные методы решений. Отметим, что задача для скоростей и давлений содержит одномерное ядро, которое следует учитывать в алгоритме ее решения [3]. Описанная схема реализована в 3D варианте, на примерах численных решений модельных задач обсуждаются свойства алгоритмов.

Список литературы

1. Ivanov M.I., Kremer I.A., Laevsky Yu.M. A Computational Model of Fluid Filtration in Fractured Porous Media // Num. Analysis and Appl. 2021. V. 14. № 2. P. 126-144. DOI: 10.1134/S1995423921020038.

2. M. I. Ivanov, I. A. Kremer, Yu. M. Laevsky, Schemes for solving filtration problem of a heat-conducting two-phase liquid in a porous medium // International conference Marchuk scientific readings 2021: Abstracts of the Intern. conf., October 4-8, 2021, Novosibirsk, DOI 10.24412/CL-35064-2021-022.

3. M. I. Ivanov, I. A. Kremer, Yu. M. Laevsky, Algorithm for the numerical solution of the pure Neumann problem in fractured porous media // International conference Marchuk scientific readings 2021: Abstracts of the Intern. conf., October 4-8, 2021, Novosibirsk, DOI: 10.24412/CL-35064-2021-020.

О методах бидиагонализации для решения несимметричных СЛАУ

В. П. Ильин1, Д. И. Козлов1,2, А. В. Петухов1

хИнститут вычислительной математики и математической геофизики 2Новосибирский государственный университет Email: ilin@sscc.ru

DOI: 10.24412/cl-35065-2022-1-02-24

Рассматриваются итерационные алгоритмы в подпространствах Крылова для решения несимметричных систем линейных алгебраических уравнений (СЛАУ) с большими разреженными матрицами, возникающими при аппроксимациях многомерных краевых задач с помощью конечных разностей, конечных объемов, конечных элементов и разрывных методов Галёркина различных порядков точности на неструктурированных сетках.

Конструируются методы сопряжённых невязок и минимальных ошибок на основе соответствующих процессов ортогонализации направляющих векторов с применением весовой бидиагонализации

i Надоели баннеры? Вы всегда можете отключить рекламу.