Статья посвящена развитию гибридного метода крупных частиц применительно к двумерным течениям с развитием физической неустойчивости на поверхностях раздела неоднородных газовых смесей. Высокая разрешающая способность метода продемонстрирована при решении задач взаимодействия ударной волны с цилиндрическим пузырем легкого или тяжелого газов в сравнении с экспериментом и расчетами по другим схемам повышенного порядка аппроксимации.
Идентификаторы и классификаторы
- eLIBRARY ID
- 43988501
Задачи с развитием неустойчивости (Рихтмайера–Мешкова, Кельвина–Гельмгольца и др.) на контактных поверхностях относятся к традиционно сложным проблемам вычислительной гидродинамики. Моделирование течений многокомпонентных неоднородных газовых смесей с возможностью многомасштабного разрешения вихрей предъявляют к численным методам ряд требований: малая численная диссипация и монотонность схемы, высокая разрешающая способность на гладких решениях, а также консервативность дискретных аналогов законов сохранения. Современные численные методы основаны на римановских солверах, схемах TVD (Total Variation Diminishing), переменных шаблонах (ENO/WENO — Weighted Essentially Non-Oscillatory), алгоритме КАБАРЕ и др. [1–9]. Подробный анализ подходов при моделировании сжимаемых многокомпонентных течений приводится в [3, 5, 8].
В настоящей статье развивается вычислительная технология, основанная на общих положениях метода крупных частиц с расщеплением алгоритма на лагранжев, эйлеров и заключительные этапы [10]. Метод модифицирован до второго порядка аппроксимации по пространству и времени. Регуляризация численного решения обеспечивается двумя способами. На лагранжевом этапе давление на гранях ячейки корректируется нелинейной искусственной вязкостью таким образом, чтобы ее величина стремилась к нулю на гладких решениях независимо от разрешения сетки. На эйлеровом и заключительном этапах вычисляются примитивные переменные (плотность, скорость и полная энергия) на гранях путем взвешенной ограничителем квазилинейной комбинации центральной и противопоточной аппроксимаций, из которых формируются дивергентные потоки. При этом выполняются дискретные аналоги законов сохранения. Повышение порядка аппроксимации по времени осуществляется методом предиктор–корректор.
К достоинствам гибридного метода крупных частиц относятся простота и однородность вычислительного алгоритма, высокая разрешающая способность и экономичность. Схема не “испытывает” проблем при прохождении звуковых точек, т.е. не требуется специальных процедур для их разрешения. Метод отличается универсальностью и возможностью решения задач негиперболического типа с учетом неоднородных законов сохранения с разномасштабными во времени компонентами решения (жесткостью) [11–14].
Список литературы
- Marquina A., Mulet P. A flux-split algorithm applied to conservative models for multicomponent compressible flows // Journal of Computational Physics. 2003. 185, N 1. 120-138. EDN: MTSQDL
- Niederhaus J.H.J. A computation parameter study for three-dimensional shock-bubble interactions. Ph.D. thesis. Madison: University of Wisconsin, 2007.
- Shankar S.K., Kawai S., Lele S. Numerical simulation of multicomponent shock accelerated flows and mixing using localized artificial diffusivity method // Proc. 48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. Reston: American Institute of Aeronautics and Astronautics, 2010. DOI: 10.2514/6.2010-352 EDN: OECRYX
- Иванов И.Э., Крюков И.А. Численный алгоритм моделирования двухфазных течений, содержащих границы раздела фаз // Физико-химическая кинетика в газовой динамике. 2012. 13, вып. 4. http://chemphys.edu.ru/issues/2012-13-4/articles/369/. EDN: SCZTDR
- Coralic V., Colonius T. Finite-volume WENO scheme for viscous compressible multicomponent flows // Journal of Computational Physics. 2014. 274. 95-121.
- Shyue K.-M., Xiao F. An Eulerian interface sharpening algorithm for compressible two-phase flow: the algebraic THINC approach // Journal of Computational Physics. 2014. 268. 326-354.
- Данилин А.В., Соловьев А.В., Зайцев А.М. Модификация схемы “КАБАРЕ” для численного моделирования течений многокомпонентных газовых смесей в двумерных областях // Вычислительные методы и программирование. 2015. 16. 436-445. EDN: YTTYQH
- Wong M.L., Lele S.K. High-order localized dissipation weighted compact nonlinear scheme for shock- and interfacecapturing in compressible flows // Journal of Computational Physics. 2017. 339. 179-209.
- Wang B., Xiang G., Hu X.Y. An incremental-stencil WENO reconstruction for simulation of compressible two-phase flows // International Journal of Multiphase Flow. 2018. 104. 20-31.
-
Белоцерковский О.М., Давыдов Ю.М. Нестационарный метод "крупных частиц" для газодинамических расчетов // Журнал вычислительной математики и математической физики. 1971. 11, № 1. 182-207.
-
Садин Д.В. TVD-схема для жестких задач волновой динамики гетерогенных сред негиперболического неконсервативного типа // Журнал вычислительной математики и математической физики. 2016. 56, № 12. 2098-2109. EDN: XGWCPL
-
Садин Д.В. Схемы с настраиваемыми диссипативными свойствами для численного моделирования течений газа и газовзвесей // Математическое моделирование. 2017. 29, № 12. 89-104. EDN: ZVGPBV
-
Садин Д.В. Применение схемы с настраиваемыми диссипативными свойствами к расчету течений газа с развитием неустойчивости на контактной границе // Научно-технический вестник информационных технологий, механики и оптики. 2018. 18, № 1. 153-157. EDN: YQLEEA
-
Садин Д.В. Использование гибридного метода крупных частиц для расчета течений многокомпонентных газовых смесей // Вычислительные методы и программирование. 2019. 20. 489-497. EDN: GTQPEA
-
Shi J., Zhang Y.-T., Shu C.-W. Resolution of high order WENO schemes for complicated flow structures // Journal of Computational Physics. 2003. 186, N 2. 690-696. EDN: MTSSSJ
-
Abgrall R. How to prevent pressure oscillations in multicomponent flow calculations: a quasi conservative approach // Journal of Computational Physics. 1996. 125, N 1. 150-160.
-
Abgrall R., Karni S. Computations of compressible multifluids // Journal of Computational Physics. 2001. 169, N 2. 594-623. EDN: MTRNIP
-
Садин Д.В., Беляев Б.В., Давидчук В.А. Сравнение модифицированного метода крупных частиц с некоторыми схемами высокой разрешающей способности. Двумерные тесты // Вычислительные методы и программирование. 2019. 20. 337-345. EDN: YEFCNT
-
Годунов С.К. Разностный метод численного расчета разрывных решений уравнений гидродинамики // Мат. сборник. 1959. 47, № 3. 271-306. EDN: SIYXCJ
-
Садин Д.В. О сходимости одного класса разностных схем для уравнений нестационарного движения газа в дисперсной среде // Журнал вычислительной математики и математической физики. 1998. 38, № 9. 1572-1577.
-
Садин Д.В. Модификация метода крупных частиц до схемы второго порядка точности по пространству и времени для ударно-волновых течений газовзвеси // Вестник Южно-Уральского государственного университета. Серия: Математическое моделирование и программирование. 2019. 12, № 2. 112-122. EDN: DHMEXP
-
Haas J.-F., Sturtevant B. Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities // J. Fluid Mech. 1987. 181. 41-76.
-
Quirk J.J., Karni S. On the dynamics of a shock-bubble interaction // J. Fluid Mech. 1996. 318. 129-163.
Выпуск
Другие статьи выпуска
В статье рассматриваются результаты экспериментальной оценки производительности и энерго-эффективности многоядерных процессоров 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/