Рубрика: "1.2.2. Математическое моделирование, численные методы и комплексы программ (технические науки)"
519.17 Задача о преследовании в 3D-пространстве с произвольными начальными углами прицеливания
Бодряков В. Ю. (Уральский государственный педагогический университет)
doi: 10.18698/2309-3684-2024-2-6884
В статье впервые получено аналитическое решение задачи о преследовании в системе «хищник–жертва» в евклидовом 3D-пространстве для произвольных начальных углов прицеливания. В процессе преследования жертва движется равномерно и прямолинейно, постоянный по модулю вектор скорости хищника нацелен на жертву. Точное решение задачи получено в форме параметрически заданной пространственной кривой преследования. Получены точные выражения для других существенных характеристик процесса преследования (времени преследования, координат жертвы, длины кривой преследования, и др.). Проведено реалистичное компьютерное моделирование взаимного движения хищника и жертвы в пространстве и во времени; определены характерные параметры процесса преследования. Отмечен значительный дидактический потенциал решенной задачи о 3D-преследовании для подготовки будущих специалистов в области механики и управления; задача может служить содержательной основой для выполнения обучающимися исследовательских проектов, курсовых и выпускных квалификационных работ.
Бодряков В.Ю. Задача о преследовании в 3D-пространстве с произвольными начальными углами прицеливания. Математическое моделирование и численные методы, 2024, № 2, с. 68-84.
Пархоменко В. П. (МГТУ им.Н.Э.Баумана/ФИЦ ИУ РАН)
doi: 10.18698/2309-3684-2021-2-3853
Целью данной работы является построение глобальной модели цикла углерода. Модель описывает продукционный процесс лесных экосистем с учетом сезонного хода климатических факторов. Она предназначена для моделирования длительного периода времени в составе глобальной климатической модели промежуточной сложности. Установлено, что глобальные характеристики климатической системы выходят на установившейся режим за время около 2000 лет и модель устойчиво работает. Приведены временные и пространственные распределения полученных климатических характеристик и биогеохимического углеродного цикла наземной растительности.
Пархоменко В.П. Глобальная климатическая модель с учетом биогеохимического углеродного цикла растительности суши. Математическое моделирование и численные методы, 2021, № 2, с. 38–53.
519.8 Диффузионная модель эволюции кластера в металлическом расплаве жаропрочного никелевого сплава
Тягунов А. Г. (Уральский Федеральный Университет), Зейде К. М. (Universidad Politècnica Salesiana/University of Genoa), Мильдер О. Б. (Уральский Федеральный Университет), Тарасов Д. А. (Уральский Федеральный Университет)
doi: 10.18698/2309-3684-2023-2-332
В работе производится построение математической модели термо-временной эволюции кластера в расплаве жаропрочного никелевого сплава ЖС6У. Формулируется начально-краевая задача с движущейся границей, для решения которой применяется численное моделирование методом трассировки траектории частиц, а для описания эволюционных процессов используется ряд классических физических теорий. Для проверки точности модели привлекается физический эксперимент построения политерм и изотерм электросопротивления рассматриваемого сплава. Подтверждено, что модель броуновской диффузии и теория проводимости Друде применимы для описания, как временной, так и температурной эволюции кластера. Так же оправдал себя подход к моделированию на основе «твердых шаров». По результатам моделирования, во временном диапазоне от 1690 до 1752 К количество частиц в составе кластера меняется от 5000 до 2000, средняя динамическая вязкость кластера изменяется от 3 до 2 *1010 Па*с, однако предполагается, что центральная часть существенно плотнее периферии, радиус кластера изменяется от 24 до 18, радиус свободной зоны вокруг кластера – от 56 до 43. Определены направления дальнейшего развития модели.
Тягунов А.Г., Зейде К.М., Мильдер О.Б., Тарасов Д.А. Диффузионная модель эволюции кластера в металлическом расплаве жаропрочного никелевого сплава. Математическое моделирование и численные методы, 2023, № 2, с. 3–32.
Иванов М. Ю., Бушуев А. Ю. (МГТУ им.Н.Э.Баумана), Щербаков Н. С. (МГТУ им.Н.Э.Баумана), Реш Г. Ф.
doi: 10.18698/2309-3684-2024-3-100119
В различных технических системах для обеспечения синхронного перемещения исполнительных органов широко применяются гидравлические устройства - нерегулируемые дроссели, делители потока, регуляторы и/или стабилизаторы расхода. Последние характеризуются тем, что их функционирование происходит в диапазоне перепадов давлений жидкости, составляющем несколько сотен атмосфер. Рассмотрены вопросы, связанные с численным моделированием нестационарных физических процессов в стабилизаторе расхода, конструкция которого защищена патентом Российской Федерации на изобретение. Представлены результаты компьютерного моделирования на основе теоретической модели с сосредоточенными параметрами, использования конечно-разностного неявного метода Гира для решения системы жёстких дифференциальных уравнений. Сформулирована и решена задача оптимального усовершенствования конструкции такого стабилизатора расхода в соответствии с выбранным критерием. Этим критерием оптимизации является обеспечение условия минимально возможного положительного статизма расходно-перепадной (статической) характеристики в условиях широкого изменения перепада давления на устройстве и воздействия осевой составляющей гидродинамической силы. Задача оптимального усовершенствования конструкции решалась с применением одного из широко используемых эволюционных алгоритмов оптимизации генетического алгоритма с вещественным кодированием. Результаты вычислительных экспериментов при моделировании физических процессов задачи анализа соответствуют имеющимся экспериментальным данным, которые ранее получены авторами работы. Показано, что усовершенствование существующей конструкции стабилизатора расхода возможно угол наклона расходно-перепадной характеристики к горизонтальной оси уменьшился практически в два раза. При этом удалось получить более высокую точность поддержания объёмного расхода жидкости. Эта точность составляет порядка ±7,5 % от номинального (настроечного) значения стабилизатора расхода. Для сравнения, точность поддержания объёмного расхода жидкости до выполнения процедуры оптимизации составляла порядка ±10 %.
Иванов М.Ю., Бушуев А.Ю., Щербаков Н.С., Реш Г.Ф. Компьютерное моделирование динамических процессов в гидравлическом стабилизаторе расхода и его оптимизация на основе эволюционного алгоритма. Математическое моделирование и численные методы, 2024, № 3, с. 100-119.
Крееренко С. С. (ПАО «ТАНТК им. Г.М. Бериева»), Крееренко О. Д. (ПАО «ТАНТК им. Г.М. Бериева»)
doi: 10.18698/2309-3684-2024-3-8199
Рассматривается задача моделирования продольного движения самолета транспортной категории и параметрическая идентификация аэродинамических характеристик продольного движения: составляющих безразмерных коэффициентов аэродинамической подъемной силы и момента тангажа. Задача решается в классе модульных полуэмпирических динамических моделей, созданных объединением теоретического и нейросетевого моделирования. Работоспособность и практическая значимость моделей подтверждается результатами вычислительных экспериментов. Разработка нейросетевой модели продольного движения самолета выполнена на языке Python с использованием открытой программной библиотеки Tensorflow для машинного обучения и высокоуровневого API Keras в составе Tensorflow.
Крееренко С.С., Крееренко О.Д. Моделирование и параметрическая идентификация аэродинамических характеристик самолета транспортной категории с использованием нейросетей в среде Тensorflow. Математическое моделирование и численные методы, 2024, № 3, с. 81–99.
Донской И. Г. (Институт систем энергетики им. Л.А. Мелентьева СО РАН)
doi: 10.18698/2309-3684-2021-3-2441
В статье рассматривается численная модель течения газа в пористой среде, содержащей частицы реакционноспособного компонента (полимера). При нагреве эти частицы расширяются, деформируются и заполняют порозное пространства, в результате чего проницаемость существенно снижается. Связь между пористостью и проницаемостью описывается формулой Козени-Кармана. Тогда вблизи нижней (входной) границы образуется область с низкой проницаемостью (агломерат), рост которой определяется условиями на боковой и входной границе. В результате расчетов получены характерные сценарии блокировки пористой среды при разных температурах нагрева. Показано, что при нагреве через стенку полимер разлагается, и пористая среда частично восстанавливает проницаемостью При нагреве поступающим газом агломерат намного более устойчив, поскольку он блокирует источник нагрева.
Донской И.Г. Численное моделирование процессов образования, роста и разложения агломератов в пористой среде при разных режимах нагрева. Математическое моделирование и численные методы, 2021, № 3, с. 24–41.
519.63 Молекулярно-динамическое моделирование модификации алюминия лазерной ударной волной
Перов Е. А. (Объединенный институт высоких температур РАН), Жаховский В. В. (Всероссийский научно-исследовательский институт автоматики им. Н.Л. Духова/Объединенный институт высоких температур РАН), Иногамов Н. А. (Институт теоретической физики им. Л.Д. Ландау РАН), Шепелев В. В. (Институт автоматизации проектирования РАН), Фортова С. В. (Институт автоматизации проектирования РАН), Долуденко А. Н. (Объединенный институт высоких температур РАН)
doi: 10.18698/2309-3684-2023-4-7492
Пластические деформации лежат в основе такой промышленной технологии, как лазерное термоупрочнение или лазерный пиннинг (LSP, laser shock peening). В данной работе методом классической молекулярной динамики исследована возможность упрочнения поверхностного слоя алюминиевого образца, облученного единичным фемтосекундным лазерным импульсом. Рассмотрены три ориентации кристаллической решетки — [1, 0, 0] (первая ориентация кристаллической решетки), [1, 1, 0] (вторая ориентация кристаллической решетки), [1, 1, 1] (третья ориентация кристаллической решетки). Проведено численное исследование влияния различных величин вложенной энергии в диапазоне от 120,98 Дж/м2 до 2540,01 Дж/м2 лазерного импульса на глубину залегания пластических деформаций, влияющих на упрочнение исследуемого материала. Построена зависимость максимальной глубины залегания пластических деформаций от вложенной энергии. Значения энергий подобранны таким образом, что пластический фронт УВ (ударной волны) останавливался до того, как достигнет правой границы моделируемого образца. Необходимость соблюдения этого условия обусловлена тем фактом, что отразившаяся от правой границы образца волна растяжения может тормозить пластический ударный фронт, выступая в роли волны разгрузки. С помощью построенной в работе зависимости максимальной глубины залегания пластических деформаций от вложенной энергии определено пороговое значение вложенной энергии, при превышении которого алюминий начинает пластически деформироваться.
Перов E.А., Жаховский В.В., Иногамов Н.А., Шепелев В.В., Фортова С.В., Долуденко А.Н.. Молекулярно-динамическое моделирование модификации алюминия лазерной ударной волной. Математическое моделирование и численные методы, 2023, № 4, с. 74-92
Чуев В. Ю. (МГТУ им.Н.Э.Баумана), Дубограй И. В. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2022-2-102113
С использованием метода динамики средних разработана «смешанная» модель противоборства многочисленных группировок при линейных зависимостях от времени эффективных скоростей нанесения воздействий единицами сторон. Построен алгоритм, позволяющий исследовать ход протекания процесса и вычислить его основные показатели. Установлено, что использование моделей с постоянными эффективными скоростями нанесения воздействий во многих случаях приводит к значительным ошибкам при вычислении основных показателей процесса. Исследовано влияние упреждающего воздействия одной из противоборствующих сторон на ход его протекания и окончательный итог.
Чуев В.Ю., Дубограй И.В. «Смешанная» модель противоборства многочисленных группировок при линейных зависимостях от времени эффективных скоростей воздействий единицами сторон. Математическое моделирование и численные методы, 2022, № 2, с. 104–115
Мозжорина Т. Ю. (МГТУ им.Н.Э.Баумана), Рахманкулов Д. А. (МГТУ им.Н.Э.Баумана)
doi: 10.18698/2309-3684-2021-3-7487
В данной работе рассматривается оптимизация перелета спутника малой массы с орбиты Земли на орбиту Марса под солнечным парусом. Оптимизация управления углом установки солнечного паруса проводится с использованием принципа максимума Понтрягина при минимизации времени перелета. В отличие от предшествующих работ на эту тему решение краевой задачи, к решению которой сводится принцип максимума, получено методом пристрелки. Программа расчета написана на языке программирования С++. Несмотря на вычислительные сложности, возникающие при использовании метода пристрелки, удалось добиться хорошей сходимости метода Ньютона, лежащего в основе алгоритма. Проведен анализ точности полученных результатов и показана возможность применения метода пристрелки при решении подобных задач. Проведено сравнение с данными ранее опубликованных работ. Несмотря на некоторые допущения, использованные при разработке алгоритма расчета, работа имеет свою ценность в плане оценки возможности использования метода пристрелки, дающего наиболее точные численные результаты оптимизации.
Мозжорина Т.Ю., Рахманкулов Д.А. Моделирование и оптимизация управлением спутника малой массы при перелете с орбиты Земли на орбиту Марса под солнечным парусом. Математическое моделирование и численные методы, 2021, № 3, с. 74–87.