Представлен численный метод решения задач динамики сжимаемого совершенного невязкого газа в одно- и двумерной пространственных постановках. Метод использует подход разрывного метода Галеркина для аппроксимации решения по пространству и снос решения вдоль характеристик для интегрирования уравнений по времени. Такой подход позволяет избавиться от стандартной многостадийной процедуры метода Рунге - Кутты и существенно уменьшить вычислительную сложность метода. Вместо интегрирования по времени применено интегрирование потоков по пространственной области зависимости на предыдущем временном слое. Представлено обобщение метода на двумерный случай для произвольной формы ячейки сетки. Проведены тестовые расчеты на ряде одно- и двумерных задач, имеющих как гладкое, так и разрывное решение. Тесты показали, что метод обладает вторым порядком точности при использовании кусочно-линейной реконструкции решения на ячейке и третьим порядком при использовании квадратных полиномов. Метод обладает высокой разрешающей способностью и малой диссипацией, позволяет разрешать газовые неустойчивости, такие как неустойчивость Рихтмайера - Мешкова. Несмотря на высокий порядок точности, метод не требует применения специальных процедур монотонизации решения
Идентификаторы и классификаторы
- SCI
- Физика
Математическая модель, включающая в себя систему уравнений динамики сжимаемого совершенного невязкого газа [1], описывает многие задачи астрофизики [2], высокоскоростных течений в технических устройствах, обтекания быстро летящих в атмосфере тел. Важнейшим инструментом исследования таких моделей являются численные методы, которые должны учитывать сложный волновой характер решения уравнений газовой динамики, а также свойство переноса компонентов решения вдоль характеристик этой гиперболической системы. Для решения первой задачи широкое применение получили приближенные решения задачи Римана о распаде разрыва и кусочно-полиномиальные реконструкции решения на ячейках расчетной сетки. Среди таких методов важное место занимает разрывный метод Галеркина [3], позволяющий повышать порядок аппроксимации без увеличения шаблона.
Если у вас возникли вопросы или появились предложения по содержанию статьи, пожалуйста, направляйте их в рамках данной темы.
Список литературы
1. Зарубин В.С., Кувыркин Г.Н. Математические модели механики и электродинамики сплошной среды. М., Изд-во МГТУ им. Н.Э. Баумана, 2008. EDN: QJUOAF
2. Lukin V., Postnov K., Shakura N. Modeling of the L1 point jet formation in binary star Her X1. Lobachevskii J. Math., 2024, vol. 45, no. 1, pp. 85-94. DOI: 10.1134/S1995080224010360 EDN: HNJOPF
3. Galanin M.P., Lukin V.V., Chechetkin V.M. 3D hydrodynamical simulation of accretion disk in binary star system using RKDG CFD solver. J. Phys.: Conf. Ser., 2018, vol. 1103, art. 012019. DOI: 10.1088/1742-6596/1103/1/012019 EDN: FBBBIJ
4. Shu C.-W., Qiu J. A comparison of troubled-cell indicators for Runge - Kutta discontinuous Galerkin methods using weighted essentially nonoscillatory limiters. SIAM J. Sci.Comput., 2005, vol. 27, iss. 3, pp. 995-1013. DOI: 10.1137/04061372X
5. Qiu J., Shu C.-W.Runge - Kutta discontinuous Galerkin method using WENO limiters. SIAM J. Sci.Comput., 2005, vol. 26, iss. 3, pp. 907-929. DOI: 10.1137/S1064827503425298
6. Qiu J.-X., Shu C.-W. Hermite WENO schemes and their application as limiters for Runge - Kutta discontinuous Galerkin method: one-dimensional case. J.Comput. Phys., 2003, vol. 193, iss. 1, pp. 115-135. DOI: 10.1016/j.jcp.2003.07.026 EDN: MTTCXJ
7. Jiang G., Shu C.-W. Efficient implementation of weighted ENO schemes. J.Comput. Phys., 1996, vol. 126, iss. 1, pp. 202-228. DOI: 10.1006/jcph.1996.0130
8. Галепова В.Д., Лукин В.В., Марчевский И.К. и др. Сравнительное исследование лимитеров семейства WENO и Hermite WENO для расчета одномерных течений газа методом RKDG. Препринты ИПМ им. М.В. Келдыша, 2017, № 131. DOI: 10.20948/prepr-2017-131 EDN: ZWTJZT
9. Lukin V.V., Korchagova V.N., Sautkina S.M. On stable Runge - Kutta methods for solving hyperbolic equations by the discontinuous Galerkin method. Diff. Equat., 2021, vol. 57, no. 7, pp. 921-933. DOI: 10.1134/S0012266121070089 EDN: DANAXA
10. Лукин В.В. Устойчивые методы Рунге - Кутты второго и третьего порядка для решения задач газовой динамики разрывным методом Галеркина. Препринты ИПМ им. М.В. Келдыша, 2022, № 52. DOI: 10.20948/prepr-2022-52 EDN: PGZSQI
11. Gatsuk V.K., Lukin V.V. Mixed characteristic discontinuous Galerkin approach for perfect gas dynamics modeling. J. Phys.: Conf. Ser., 2021, vol. 2028, art. 012018. DOI: 10.1088/1742-6596/2028/1/012018 EDN: UJMBWJ
12. Попов Ю.П., Самарский А.А. Полностью консервативные разностные схемы для уравнений газодинамики в переменных Эйлера. Ж. вычисл. матем. и матем. физ., 1970, т. 10, № 3, с. 773-779.
13. Toro E.F. Riemann solvers and numerical methods for fluid dynamics. Berlin, Heidelberg, Springer, 2013. DOI: 10.1007/b79761
14. Куликовский А.Г., Погорелов Н.В., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. М., ФИЗМАТЛИТ, 2001. EDN: OUSHIA
15. Попов М.В., Устюгов С.Д. Кусочно-параболический метод на локальном шаблоне для задач газовой динамики. Ж. вычисл. матем. и матем. физ., 2007, т. 47, № 12, c. 2055-2075. EDN: IBMNPP
16. Collela P., Woodward P. The Piecewise Parabolic Method (PPM) for gasdynamical simulations. J.Comp. Phys., 1984, vol. 54, iss. 1, pp. 174-201. DOI: 10.1016/0021-9991(84)90143-8 EDN: XVGPZI
17. Уcтюгoв C.Д., Пoпoв М.B. Кycoчнo-пapaбoличecкий мeтoд нa лoкaльнoм шaблoнe. I. Уpaвнeниe пepeнoca и ypaвнeниe Бюpгepca. Препринты ИПМ им. М.В. Келдыша, 2006, № 65.
18. Уcтюгoв C.Д., Пoпoв М.B. Кусочно-параболический метод на локальном шаблоне. II. Уравнения газодинамики. Препринты ИПМ им. М.В. Келдыша, 2006, № 71.
19. Kurganov A., Tadmor E. Solution of two-dimensional Riemann problems for gas dynamics without Riemann problem solvers. Numer. Methods Partial Differ. Equ., 2002, vol. 18, iss. 5, pp. 584-608. DOI: 10.1002/num.10025
20. Helzel C., Kerkmann D., Scandurra L. A new ADER method inspired by the active flux method. J. Sci.Comput., 2019, vol. 80, no. 3, pp. 1463-1497. DOI: 10.1007/s10915-019-00988-1 EDN: QNLZWM
21. Годунов С.К. Разностный метод численного расчета разрывных решений уравнений гидродинамики. Матем. сб., 1959, т. 47, № 89, с. 271-306. EDN: SIYXCJ
Выпуск
Другие статьи выпуска
Приведены результаты исследований щелочной обработки сульфатсодержащего раствора, состоящего в основном из сульфатов алюминия и железа, с использованием сульфатизации каолиновых глин месторождения «Чашма-Санг» с целью получить гидроксид алюминия. В лаборатории изучен процесс щелочной обработки сульфатсодержащего раствора в зависимости от температуры, продолжительности процесса, концентрации гидроксида натрия и каустического модуля. Определено, что при температуре 85 °С, времени процесса 30 мин, концентрации гидроксида натрия 10 % и каустического модуля 3,0 степень извлечения тетрагидроксоалюмината натрия в раствор составляет более 99 %, осаждение гидроксида железа - более 98 %, а оставшееся количество примесного гидроксида алюминия в составе гидроксида железа - 0,9 %. Карбонизацией раствора тетрагидроксоалюмината натрия получен гидроксид алюминия, в котором содержание Al(OH)3 составляет более 99 %. Рентгенофазовым анализом установлено, что полученный гидроксид алюминия соответствует минералу байерит (номер по картотеке PDF 74-1119). Согласно результатам физико-химического анализа, полученный гидроксид алюминия соответствует требованиям ГОСТ 11841-76 для гидроксида алюминия технического и ТУ 1711-001-00658716-99. На основе лабораторных исследований разработана принципиальная технологическая схема переработки каолиновых глин месторождения «Чашма-Санг» в целях получения гидроксида алюминия. В соответствии с результатами проведенных исследований и разработанной технологической схемой в цехе газоочистки ОАО «ТАлКо» опробована технология сульфатизации каолиновых глин и щелочной обработки сульфатсодержащего раствора с получением опытного количества гидроксида алюминия
Приведены результаты исследования процесса взаимодействия фосфата цинка с биологически активными веществами с использованием метода квантово-химического моделирования. В рамках метода получены модели молекул биологически активных веществ (незаменимых аминокислот), а также молекулярных систем Zn3(PO4)2-незаменимая аминокислота, в которых взаимодействие фосфата цинка с аминокислотой происходило через атом кислорода, присоединенный к атому фосфора в фосфате цинка, и ионизированную аминогруппу аминокислоты. В результате квантово-химического моделирования установлено, что молекулярная система Zn3(PO4)2-незаменимая аминокислота является энергетически выгодной и химически стабильной. Оптимальной конфигурацией, обладающей наибольшей разностью полной энергии (73,646 ккал/моль) и химической жесткостью (0,144 эВ), является молекулярная система Zn3(PO4)2-лизин, в которой взаимодействие происходит через ионизированную Å-аминогруппу лизина. Результаты анализа микроструктуры наночастиц фосфата цинка, стабилизированных L- лизином, показали, что поверхность образцов состоит из крупных агрегатов диаметром 6…17 мкм, которые представляют собой мелкие частицы сферической формы размерами 100…150 нм. Результаты анализа элементного состава позволили установить наличие в структуре таких элементов, как цинк и фосфор. Изучен фазовый состав образцов. Установлено, что образец наночастиц фосфата цинка, стабилизированных L- лизином, имеет орторомбическую сингонию с пространственной группой P nma. Образцы исследовали методом ИК-спектроскопии. По результатам спектроскопии определено, что взаимодействие аминокислоты происходит при связывании фосфата через аминогруппы в молекуле аминокислоты
Представлены результаты численного 3D-моделирования свободно-конвективных процессов теплообмена в каналах одного из основных вариантов панельно-лучистой системы обогрева помещений. Выявлены особенности гидродинамики и теплообмена в вертикальных каналах прямоугольного сечения с наличием горизонтально расположенных труб и поперечно установленных на них ребер, при различных условиях теплоотдачи на ограничивающих каналы поверхностях адиабатной стенки и радиационного экрана. Определены пространственные распределения температуры и скорости воздушного потока в системе для двух возможных режимов теплоотдачи. Режимы реализуются в различных диапазонах изменения модифицированного числа Рэлея. Предложены критериальные зависимости чисел Рейнольдса и Нуссельта от числа Рэлея и основных конструктивных параметров системы. Получена формула для расчета мощности удельного (на единицу длины) теплосъема с греющих поверхностей внутренних элементов всей плинтусной системы. Доказано, что удельная мощность теплосъема принимает максимальное значение при определенном значении межреберного зазора, при котором имеет место переходный режим теплообмена. Установлены соотношения между мощностями конвективных и радиационных тепловых потоков, генерируемых в рассмотренной теплообменной системе. Сравнение расчетных и существующих экспериментальных данных по мощности удельного теплосъема указывает на их удовлетворительное совпадение
В этой статье алгоритм оптимизации российского хоккея с шайбой (ORIHO), основанный на противодействии, применяется для решения задачи уменьшения потерь энергии. Алгоритм оптимизации российского хоккея с шайбой основан на методах игры в хоккей с шайбой, в которых особое внимание уделяется быстрому забросу шайбы и определению местоположения соперника. Хоккей с шайбой - это командная игра connexion в середине зимы, в которую играют на ледяных скейтбордах, обычно на покрытой инеем скользкой арене с полосами и раскрасками, характерными для игры. Местоположение участника или игрока указывает на решение соперника. В промежуточный период местоположение хоккейной шайбы является вектором определения, который будет определять местоположение следующего участника. Кроме того, игроки определяют набор определений. Хаотические последовательности интегрированы в алгоритм ORIHO. Внедрены стандарты создания хаотических карт Tinkerbell. Алгоритм обучения ORIHO, основанный на противопоставлении, использует распределение Лапласа для улучшения исследовательского навыка. Затем, изучая перспективы расширения исследования, мы разработали новый метод.
Для описания нелинейного деформирования ортотропных материалов при силовом нагружении предложено использовать эндохронную теорию пластичности ортотропных материалов. Определяющие соотношения получены введением тензорного внутреннего параметра состояния и тензора времени релаксации. Эти соотношения представлены в интегральной и дифференциальной формах. Рассмотрены различные виды тензора времен релаксации и функции упрочнения и предложены три модели теории. Для каждой модели получены зависимости при простом нагружении, модели без функции упрочнения представлены аналитическими выражениями. Анализ этих зависимостей позволил установить связь между материальными параметрами моделей и их влиянием на нелинейное поведение материала, а также предложить методику их экспериментального определения. Проведен аналогичный анализ для многослойных композитных материалов. Сопоставление расчетно-теоретических результатов с экспериментальными данными (например, для трансверсально-изотропного материала бороалюминия) показало, что модель с одним временем релаксации не подходит для описания трехосного нагружения, а модель с несколькими временами релаксации и функцией упрочнения наиболее точно соответствует экспериментальным данным. Для многослойных материалов без явно выраженного линейного упрочнения (например, AS4/PEEK) возможно использование олее простых вариантов теории без функции упрочнения. Результаты проведенного анализ вариантов эндохронной теории пластичности ортотропных материалов позволяют лучше понять влияние параметров моделей на поведение композитных материалов и могут стать основой для дальнейших исследований.
Достоверность результатов математического моделирования термонапряженного состояния теплозащитного покрытия аэрокосмического летательного аппарата зависит от достоверности оценки тепловых и механических характеристик применяемых материалов, корректности выбранных математических моделей распределения тепловых потоков, полей температур и напряжений. На примере трансверсально-изотропного композитного теплозащитного материала, состоящего из изотропной матрицы, выполненной из стеклоуглерода, и включений из трансверсально-изотропного пирографита рассмотрены методы построения двусторонних (сверху и снизу) оценок характеристик этого материала. Для расчетов выбрана слоистая структура, формируемая чередующимися этапами высокотемпературного осаждения из газовой фазы слоев пирографита и нанесения слоев стеклоуглерода. На основе метода конечных элементов проведен численный анализ термонапряженного состояния теплозащитного покрытия спускаемого аппарата, возвращающегося на Землю по завершении лунной экспедиции. В качестве примера рассмотрен спускаемый аппарат класса «Аполлон». Для получения верхней оценки интенсивности тепловых нагрузок рассмотрен спуск аппарата с нулевым углом атаки, т. е. при нулевом аэродинамическом качестве. Показана эффективность применения слоя описанного композитного структурно чувствительного материала по сравнению с двухслойным пакетом, набранным из отдельных слоев тех же материалов.
На протяжении последних трех десятилетий значительное внимание уделяется созданию математических моделей высоких порядков и исследованию их теоретических свойств. Одним из наиболее распространенных примеров является обобщенная экспоненциальная авторегрессионная модель Озаки, которая активно используется в моделировании медицинских, экономических и механических процессов. Кроме того, указанная модель является стандартом для генерации тестовых данных для искусственных нейронных сетей, способных по ряду прошлых значений случайной величины предсказать ее новые значения. Выполнен расчет асимптотической ковариационной матрицы оценки наименьших модулей в модели Озаки путем разложения функции ее ошибок по формуле Тейлора. Проведен сравнительный анализ скорости приближения отдельных реализаций модели к ее асимптотическому поведению при различных распределениях обновляющего процесса: нормального распределения и загрязненного нормального распределения, распределения Стьюдента и распределения Лапласа, логистического распределения. Научная новизна работы заключается в определении асимптотической ковариационной матрицы оценки наименьших модулей модели Озаки, а практическая значимость - в использовании результатов сравнения для выбора этой модели в инженерных расчетах
Издательство
- Издательство
- МГТУ им. Н.Э. Баумана
- Регион
- Россия, Москва
- Почтовый адрес
- 105005, г. Москва, вн. тер. г. муниципальный округ Басманный, ул. 2-я Бауманская, д. 5, с. 1
- Юр. адрес
- 105005, г. Москва, вн. тер. г. муниципальный округ Басманный, ул. 2-я Бауманская, д. 5, с. 1
- ФИО
- Гордин Михаил Валерьевич (Ректор)
- E-mail адрес
- bauman@bmstu.ru
- Контактный телефон
- +7 (499) 2636377
- Сайт
- https://bmstu.ru/