Open Access Open Access  Restricted Access Access granted  Restricted Access Subscription Access

Vol 64, No 2 (2024)

Cover Page

Full Issue

Open Access Open Access
Restricted Access Access granted
Restricted Access Subscription Access

General numerical methods

Оценка QTT рангов регулярных функций на равномерной квадратной сетке

Зыль А.В., Замарашкин Н.Л.

Abstract

В работе доказываются оценки e-рангов для TT-разложений тензоров, полученных путем тензоризации значений регулярной функции одной комплексной переменной на равномерной квадратной сетке на комплексной плоскости. Установлена связь точности приближения и геометрии области регулярности функции. Библ. 8. Фиг. 2.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):189-199
pages 189-199 views

Улучшенная квадратурная формула для прямого значения нормальной производной потенциала простого слоя

Крутицкий П.А., Резниченко И.О.

Abstract

Рассматривается потенциал простого слоя для уравнения Гельмгольца в трехмерном случае, а также потенциал простого слоя для уравнения Лапласа. Получена квадратурная формула для прямого значения нормальной производной потенциала простого слоя с непрерывной плотностью, заданной на замкнутой либо разомкнутой поверхности. Квадратурная формула, предложенная в работе, дает значительно более высокую точность, чем известные формулы, что подтверждается численными тестами. Полученная квадратурная формула может использоваться при численном решении краевых задач для уравнений Лапласа и Гельмгольца методом граничных интегральных уравнений. Библ. 17. Табл. 1.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):200-219
pages 200-219 views

Приближение непрерывных функций с помощью классических синков и значений операторов Cλ

Пасечник В.Н.

Abstract

Рассмотрены свойства синк-приближений. Используемые ранее классические синк-аппроксимации давали плохое приближение, а новый оператор, обобщающий синк-аппроксимации, справляется с приближением этой функции лучше. Приведен график численной реализации эксперимента. Библ. 22. Фиг 2.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):220-224
pages 220-224 views

Optimal control

Синтез оптимальной устойчивой аффинной системы

Ащепков Л.Т.

Abstract

Предложен метод построения обратной связи, обеспечивающей притяжение траекторий аффинной системы к состоянию равновесия и заданному многообразию. Искомая обратная связь находится в аналитической форме как решение вспомогательной задачи оптимального управления. Приведены достаточные условия существования оптимального управления. Показано применение метода к некоторым классам линейных и нелинейных систем. Библ. 10.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):225-231
pages 225-231 views

Ordinary differential equations

Гладкие многообразия Ляпунова для автономных систем нелинейных обыкновенных дифференциальных уравнений и их применение к решению сингулярных краевых задач

Конюхова Н.Б.

Abstract

Для автономной системы N нелинейных обыкновенных дифференциальных уравнений, рассматриваемой на полубесконечном интервале и обладающей точкой покоя (псевдо)гиперболического типа, изучается n-мерное устойчивое многообразие решений, или многообразие условной устойчивости по Ляпунову, которое для каждого достаточно большого t существует в фазовом пространстве переменных системы в окрестности ее седловой точки. Гладкая сепаратрисная поверхность седла для такой системы описывается с помощью решения сингулярной задачи типа Ляпунова для системы квазилинейных уравнений с частными производными первого порядка с вырождением по начальным данным. Дается применение результатов к правильной постановке граничных условий на бесконечности и их переносу в конечную точку для автономной системы нелинейных уравнений, в том числе с указанием использования этого подхода в некоторых прикладных задачах. Библ. 26.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):232-252
pages 232-252 views

Partial Differential Equations

Явные численно реализуемые формулы для операторов Пуанкаре–Стеклова

Демидов А.С., Самохин А.С.

Abstract

Представлены явные численно реализуемые формулы для операторов Пуанкаре–Стеклова, таких как операторы Дирихле–Неймана, Дирихле–Робена, Робена1–Робена2, Гринберга–Майергойза, относящихся к двумерному уравнению Лапласа. Эти формулы базируются на лемме об однолистном изометрическом отображении замкнутой аналитической кривой на окружность. Численные результаты для областей с весьма сложной геометрией получены для нескольких тестовых гармонических функций для операторов Дирихле–Неймана и Дирихле–Робена. Библ. 9. Фиг. 9.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):253-262
pages 253-262 views

Метод численного решения нестационарного уравнения Шрёдингера десятого порядка точности

Захаров М.А.

Abstract

Представлен метод численного решения нестационарного уравнения Шрёдингера десятого порядка точности, основанный на аппроксимации оператора эволюции формулой произведения. Обсуждается проблема уменьшения числа операторных экспонент в итоговой формуле за счет оптимизации их последовательности. На основе идеи, предложенной Йошида, построены два алгоритма десятого порядка точности для аппроксимации оператора эволюции. Численные тесты продемонстрировали устойчивость этих алгоритмов и их порядок точности. Метод, использованный в статье, позволил значительно уменьшить количество экспоненциальных множителей в схеме по сравнению с известной формулой Ли–Троттера–Сузуки. Библ. 25. Фиг. 2. Табл. 2.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):263-282
pages 263-282 views

Решение двумерного нелинейного параболического уравнения теплопроводности при краевом режиме, заданном на подвижном многообразии

Казаков А.Л., Нефедова О.А., Спевак Л.Ф.

Abstract

Статья посвящена исследованию вырождающегося параболического уравнения теплопроводности с нелинейностями общего вида при наличии источника (стока) в случае двух пространственных переменных. Рассмотрена задача об инициировании тепловой волны, распространяющейся по холодному (нулевому) фону с конечной скоростью, краевым режимом, заданным на подвижном многообразии — замкнутой линии. Для нее доказана новая теорема существования и единственности и предложен численный алгоритм построения решения, основанный на методе граничных элементов, методе коллокаций и разностной аппроксимации по времени, при этом использована специальная замена переменных типа преобразования годографа. Найдены новые точные решения рассматриваемого уравнения в случае нелинейностей степенного вида. Численный алгоритм реализован в виде программы, проведен комплексный вычислительный эксперимент. Сравнение построенных численных решений с точными (как найденными в работе, так и ранее известными) показало хорошее соответствие, установлена численная сходимость относительно шага по времени и числа точек коллокации. Библ. 28. Фиг. 3. Табл. 4.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):283-303
pages 283-303 views

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

Солонуха О.В.

Abstract

Доказаны достаточные условия существования обобщенного решения нелинейного эллиптического дифференциального уравнения с нелокальными краевыми условиями типа Бицадзе-Самарского. Используется условие сильной эллиптичности вспомогательного дифференциально–разностного оператора. При сформулированных условиях дифференциально–разностный оператор является деминепрерывным, коэрцитивным и обладает полуограниченной вариацией, что позволяет применять общую теорию операторов псевдомонотонного типа. Библ. 16.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):304-321
pages 304-321 views

Mathematical physics

Локальная вейвлетная адаптация декартовых сеток в вычислительных задачах газовой динамики

Афендиков А.Л., Никитин В.С.

Abstract

Представлен метод динамической локальной адаптации градуированных декартовых деревьев для численного решения задач газовой динамики. Локальный вейвлетный анализ газодинамического поля на базе неравномерных B-сплайнов применяется независимо к каждой ячейке расчетной сетки и позволяет выделить негладкие или существенно нелинейные участки решения (или наоборот, достаточно гладкие и линейные) и модифицировать сетку для расчета следущего шага по времени так, чтобы у разномасштабных особенностей течения было адекватное сеточное разрешение. В комбинации с другими методами вычислительной газовой динамики, такими как метод свободный границы, представленный метод позволяет эффективно решать нестационарные задачи с обтеканием движущихся тел. На ряде таких задач продемонстрирована работа предложенного варианта вейвлетной адаптации. Библ. 37. Фиг. 5.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):322-336
pages 322-336 views

Новые компьютерно-экономичные аппроксимации случайных функций для решения стохастических задач теории переноса

Михайлов Г.А., Медведев И.Н.

Abstract

Разработана новая сеточная аппроксимация однородного изотропного случайного поля с заданной средней корреляционной длиной. Эта аппроксимация строится путем разбиения координатного пространства на ансамбль кубиков, размер которых воспроизводит среднюю корреляционную длину при независимом выборе значения поля из заданного одномерного распределения в каждом элементе разбиения. Сформулирован также недавно предложенный авторами метод корреляционно-рандомизированного моделирования переноса частиц через случайную среду. Проведено сравнение точности и трудоемкости соответствующих алгоритмов метода Монте-Карло для решения задач о переносе гамма-квантов через случайную среду типа мозаики Вороного. Для проверки гипотезы о существенном влиянии одномерного распределения и корреляционного радиуса оптической плотности среды на перенос излучения были также проведены дополнительные расчеты для случайного пуассоновского “поля воздушных шаров” в воде. Дано обобщение сеточной аппроксимации на неизотропные случайные поля. Библ. 19. Фиг. 3. Табл. 4.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):337-349
pages 337-349 views

Исследование применения явно-итерационной схемы к моделированию дозвуковых реагирующих газовых потоков

Пескова Е.Е., Язовцева О.С.

Abstract

Статья посвящена исследованию возможности применения явно-итерационной схемы ЛИ–М для расчета диссипативных членов к решению задач дозвуковых реагирующих потоков с радикально-цепными реакциями, активными диффузионными процессами, значительной теплопередачей и энергопоглощением. Моделирование подобных течений характеризуется ограничением на шаг интегрирования по времени, связанным, в первую очередь, с преобладанием диффузионных процессов над конвективными и наличием быстрых химических реакций. Математическая модель описана с использованием многокомпонентных уравнений Навье–Стокса. Совокупность разномасштабных процессов в модели обусловила использование принципа расщепления по физическим процессам: химическая кинетика проинтегрирована методом Радо с адаптивным шагом по времени; конвективный поток рассчитывается с использованием потока Русанова и WENO схемы; диссипативные потоки с помощью явно-итерационной схемы ЛИ–М. В результате разработаны численный алгоритм и код для исследования дозвуковых реагирующих течений в осесимметричной геометрии и проведен ряд вычислительных экспериментов. Для тестирования реализованного алгоритма было решено одномерное нестационарное неоднородное уравнение. Показано, что применение схемы ЛИ–М к расчету диссипативной части позволяет избавиться от диффузионного ограничения на шаг интегрирования по времени. Было проведено численное моделирование процесса высокотемпературной конверсии метана в одномерной постановке. Данный процесс характеризуется быстрыми химическими реакциями, значительными локальными изменениями температуры, плотности газа и теплофизических характеристик, что накладывает серьезные ограничения на шаг интегрирования по времени. Показано, что алгоритм позволяет проводить расчеты с шагом, превышающим диффузионные ограничения на шаг по времени. Проведено сравнение расчетов с расчетами по ранее верифицированному алгоритму, показано хорошее совпадение результатов со значительным выигрышем по времени выполнения программы. Проведено численное моделирование течения газа в цилиндрической трубе, полученные результаты верифицированы путем сеточной сходимости. Библ. 29. Фиг.7.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):350-363
pages 350-363 views

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

Семисалов Б.В., Медведев С.Б., Назаренко С.В., Федорук М.П.

Abstract

Предложен метод численного решения четырехволновых кинетических уравнений, возникающих в задачах волновой (слабой) турбулентности при описании однородного изотропного взаимодействия волн. Для расчета интеграла столкновений разработаны быстросходящиеся кубатурные формулы, позволяющие адаптировать алгоритм к особенностям решений и ядер интегралов. Проведены эксперименты на сходимость в задачах интегрирования из реальных приложений. Для учета существенной разномасштабности задач турбулентности в алгоритме реализованы и протестированы дробно-рациональные приближения решений и новая схема итераций по времени. Эффективность разработанного алгоритма продемонстрирована при моделировании обратного каскада частиц бозе-газа при формировании конденсата Бозе–Эйнштейна. Библ. 51. Фиг. 10. Табл. 1.

Žurnal vyčislitelʹnoj matematiki i matematičeskoj fiziki. 2024;64(2):364-386
pages 364-386 views