ПРОГРАММНЫЙ КОМПЛЕКС ДЛЯ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ РАЗРУШЕНИЯ ТЕРМОПОРОУПРУГОЙ СРЕДЫ (2020)
Приведено описание программного комплекса для математического моделирования эволюции термопороупругой среды с учетом ее разрушения. Используемая математическая модель является модификацией модели Био для случая термопороупругих сред и позволяет моделировать изменение напряженно-деформированного состояния среды, фильтрацию флюида, неизотермические эффекты, а также разрушение среды. Разрушение среды описывается с использованием подхода континуальной механики разрушения путем введения дополнительной переменной, называемой параметром повреждаемости. Этот параметр характеризует степень разрушения среды, а его эволюция определяется заданным кинетическим уравнением. Вычислительный алгоритм основан на методе конечных элементов. Дискретизация уравнений по времени производится по неявной схеме для перемещений, давления и температуры и по явной для параметра повреждаемости. В качестве конечных элементов выбраны элементы Тейлора-Худа, имеющие второй порядок аппроксимации по перемещениям и первый по давлению и температуре. Система уравнений решается в рамках “монолитной” постановки без итерационного связывания между группами уравнений. Рассмотрены результаты расчетов с использованием программного модуля на примере задачи термического воздействия на нефтяной пласт.
Идентификаторы и классификаторы
- eLIBRARY ID
- 43988502
Изучение процессов разрушения материалов под действием термомеханических нагрузок является актуальным направлением исследований в различных областях физики и имеет прикладную значимость в металлургии, строительстве, добыче углеводородов и других отраслях. В частности, в нефтяном инжиниринге одной из первостепенных задач является разработка различных методов увеличения нефтеотдачи для низкопроницаемых коллекторов. Одним из таких методов является термическое воздействие, при котором в пласт закачивается теплоноситель, который увеличивает подвижность нефти. При этом теплоноситель имеет температуру, значительно превосходящую пластовую, а его закачка производится под большим давлением. В связи с этим закачка сопровождается различными процессами, происходящими в пласте, такими как деформация породы, фильтрация и изменение физико-химических свойств пластового флюида, изменение температурного поля, а также разрушение породы. При этом для моделирования процесса термического воздействия требуется корректный учет взаимного влияния указанных эффектов.
В настоящий момент существует ряд моделей, описывающих процесс эволюции пороупругой среды с учетом ее разрушения. Основными недостатками этих моделей является то, что некоторые из них в общем случае не являются термодинамически корректными [1, 2], а другие [3] имеют в представленном виде скорее теоретическую направленность и не пригодны для проведения прикладных расчетов.
В настоящей статье представлена термодинамически согласованная математическая модель для описания эволюции термопороупругой среды с учетом деформационных, фильтрационных и неизотермических процессов, а также разрушения. Для данной модели описан вычислительный алгоритм, который реализован в виде программного комплекса на языке C++, приведены примеры прикладных расчетов
Список литературы
- Krajcinovic D., Fonseka G.U. The continuous damage theory of brittle materials, part 1: general theory // Journal of applied Mechanics. 1981. 48, N 4. 809-815.
- Murakami S. Continuum damage mechanics: a continuum mechanics approach to the analysis of damage and fracture. Dordrecht: Springer, 2012.
- Кондауров В.И., Фортов В.Е. Основы термомеханики конденсированной среды. М.: МФТИ, 2002.
- Biot M.A. General theory of three dimensional consolidation // Journal of Applied Physics. 1941. 12, N 2. 155-164.
- Noll W. A mathematical theory of the mechanical behavior of continuous media // Archive for Rational Mechanics and Analysis. 1958. 2, N 1. 197-226. EDN: VTEPUK
- Griffith A.A. The phenomena of rupture and flow in solids // Philosophical Transactions of the Royal Society of London. Ser. A. 1920. Vol. 221. 163-198.
- Качанов Л.М. О времени разрушения в условиях ползучести // Изв. АН СССР. ОТН. 1958. № 8. 26-31.
- Работнов Ю.Н. Механизм длительного разрушения // Вопросы прочности материалов и конструкций. М.: АН СССР, 1959. 5-7.
- Меретин А.С., Савенков Е.Б. Математическая модель фильтрационных процессов в термопороупругой среде с учетом континуального разрушения // Препринты ИПМ им. М.В. Келдыша. 2019. № 58. 38 с. DOI: 10.20948/prepr-2019-58 EDN: PABXYT
-
Kim J., Tchelepi H.A., Juanes R. Stability, accuracy, and efficiency of sequential methods for coupled flow and geomechanics // SPE J. 2011. 16. DOI: 10.2118/119084-PA
-
Taylor C., Hood P.A. A numerical solution of the Navier-Stokes equations using the finite element technique // Computers and Fluids. 1973. 1, N 1. 73-100.
-
Brezzi F., Fortin M. Mixed and hybrid finite element methods. New York: Springer, 1991.
-
Саад Ю. Итерационные методы для разреженных линейных систем. М.: Изд-во Моск. ун-та, 2013.
-
Neuman S.P. Saturated-unsaturated seepage by finite elements // J. Hydraul. Div. 1973. 99, N 12. 2233-2250.
-
Cuthill E., McKee J. Reducing the bandwidth of sparse symmetric matrices. New York: ACM Press, 1969. 157-172.
-
C++ Template Library for Linear Algebra. URL: http://eigen.tuxfamily.org/.
-
HYPRE: Scalable Linear Solvers and Multigrid Methods. https://computation.llnl.gov/projects/hypre-scalable-linear-solvers-multigrid-methods.
-
ParaView. https://www.paraview.org/.
-
Visualization Toolkit (VTK). https://www.vtk.org.
-
Pogacnik J., O'Sullivan M., O'Sullivan J. A damage mechanics approach to modeling permeability enhancement in thermo-hydro-mechanical simulations // Proceedings of 39th Workshop on Geothermal Reservoir Engineering. Stanford: Stanford University. 2014. https://pangea.stanford.edu/ERE/pdf/IGAstandard/SGW/2014/Pogacnik.pdf.
-
Tang C.A. et al. Coupled analysis of flow, stress and damage (FSD) in rock failure // International Journal of Rock Mechanics and Mining Sciences. 2002. 39, N 4. 477-489. EDN: AYCKZN
-
Beggs H.D., Robinson J.R. Estimating the viscosity of crude oil systems // Journal of Petroleum Technology. 1975. 27, N 9. 1140-1141.
-
Sun F., Jia P., Xue S. Continuum damage modeling of hydraulic fracture from perforations in horizontal wells // Mathematical Problems in Engineering. 2019. DOI: 10.1155/2019/9304961 EDN: CURUNF
Выпуск
Другие статьи выпуска
В статье рассматриваются результаты экспериментальной оценки производительности и энерго-эффективности многоядерных процессоров MALT в задачах обработки изображений на примере фильтрации изображения с помощью оператора Собеля. Измерения осуществлялись с использованием низкоуровневого эмулятора MALTemu, прототипа процессора в ПЛИС и экспериментальной СБИС модели MALT-Cv2 Rev1. Полученные результаты сравниваются с аналогичными результатами для процессоров общего назначения (последовательная реализация) и графических процессоров с поддержкой технологии CUDA.
Предложен новый метод восстановления изображений, имеющих три неизвестные градации яркости. Для их определения используются фрагменты изображения, гистограммы которых согласуются с заданным распределением шума. Далее все пиксели распределяются по найденным уровням яркости посредством бинарной классификации. Выполнен вычислительный эксперимент, по результатам которого оказалось, что ошибка оценки исходных яркостей не превысила 3%. При относительно низком уровне шума доля неверно классифицированных пикселей от их общего числа составила менее 0.006.
Рассматривается нелинейная модификация стохастической модели галактического динамо, в рамках которой коэффициент, отвечающий за турбулентную диффузию, полагается случайным процессом с обновлением. Показано, что при малых значениях напряженности магнитного поля его статистические моменты ведут себя примерно так же, как и в линейной модели (в частности, продемонстрировано наличие эффекта перемежаемости). Получены оценки для характерных времен выхода моментов на стабилизацию, которая наступает по мере приближения поля к равновесному значению. Проведено сопоставление результатов численного эксперимента, полученных при усреднении различных объемов выборки независисмых случайных реализаций поля.
Функционалы Минковского являются важным инструментом для изучения морфологии пористых сред. Настоящая работа посвящена построению алгоритма вычисления функционалов Минковского четырехмерных цифровых изображений, возникающих, в частности, при описании динамики изменения порового пространства среды. В работе впервые программно реализован алгоритм вычисления функционалов Минковского четырехмерных цифровых изображений.
Рассматривается задача упаковки шаров двух типов в замкнутое ограниченное множество в трехмерном пространстве как с евклидовой, так и со специальной неевклидовой метрикой. Требуется максимизировать радиус шаров при известном количестве шаров каждого типа и заданном отношении между радиусами. Предложен вычислительный алгоритм, основанный на комбинации метода бильярдного моделирования и оптико-геометрического подхода, базирующегося на фундаментальных физических принципах Ферма и Гюйгенса. Приведены результаты вычислительного эксперимента.
Статья посвящена развитию гибридного метода крупных частиц применительно к двумерным течениям с развитием физической неустойчивости на поверхностях раздела неоднородных газовых смесей. Высокая разрешающая способность метода продемонстрирована при решении задач взаимодействия ударной волны с цилиндрическим пузырем легкого или тяжелого газов в сравнении с экспериментом и расчетами по другим схемам повышенного порядка аппроксимации.
Издательство
- Издательство
- МГУ
- Регион
- Россия, Москва
- Почтовый адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- Юр. адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- ФИО
- Садовничий Виктор Антонович (РЕКТОР)
- E-mail адрес
- info@rector.msu.ru
- Контактный телефон
- +7 (495) 9391000
- Сайт
- https://msu.ru/