Все выпуски

Работа посвящена использованию бикомпактных схем для численного решения эволюционных уравнений гиперболического типа. Основным преимуществом схем этого класса является сочетание двух положительных свойств: пространственной аппроксимации высокого четного порядка на шаблоне, всегда занимающем одну ячейку сетки, и спектрального разрешения, лучшего по сравнению с классическими компактными конечно-разностными схемами того же порядка пространственной аппроксимации. Рассматривается одна особенность бикомпактных схем — жесткая привязка их пространственной аппроксимации к декартовым сеткам (с ячейками-параллелепипедами в трехмерном случае). Она делает затруднительным применение бикомпактных схем к решению задач в сложных расчетных областях в рамках подхода неструктурированных сеток. Предлагается решать эту проблему путем применения известных методов аппроксимации границ сложной формы и соответствующих им краевых условий на декартовых сетках. Обобщение бикомпактных схем на задачи в геометрически сложных областях проводится на примере задач газовой динамики и уравнений Эйлера. В качестве конкретного метода, позволяющего учесть на декартовых сетках влияние твердых границ произвольной формы на течение газа, выбирается метод свободной границы. Приводится краткое описание этого метода, выписываются его уравнения. Для них строятся бикомпактные схемы четвертого порядка аппроксимации по пространству с локально-одномерным расщеплением. Компенсационный поток метода свободной границы дискретизируется со вторым порядком точности. Для интегрирования по времени в получаемых схемах применяются неявный метод Эйлера и $L$-устойчивый жестко-точный трехстадийный однократно диагонально-неявный метод Рунге–Кутты третьего порядка точности. Разработанные бикомпактные схемы тестируются на трех двумерных задачах: о стационарном сверхзвуковом обтекании с числом Маха, равным трем, одного круглого цилиндра и группы изт рех круглых цилиндров, а также о нестационарном взаимодействии плоской ударной волны и круглого цилиндра в канале с плоскопараллельными стенками. Полученные результаты хорошо согласуются с результатами других работ: твердые тела физически корректно влияют на поток газа, давление в контрольных точках на поверхностях тел рассчитывается с точностью, в целом отвечающей выбранному разрешению сетки и уровню численной диссипации.

Представлена гиперболическая модель вязкого теплопроводного газа, в которой для гиперболизации уравнений использован подход Максвелла–Каттанео, обеспечивающий распространение волн с конечными скоростями. В модифицированной модели вместо оригинальных законов Стокса и Фурье использовались их релаксационные аналоги и показано, что при стремлении времен релаксации $\tau_\sigma^{}$ и $\tau_w^{}$ к нулю гиперболизированные уравнения приводятся к классической системе Навье–Стокса негиперболического типа с бесконечными скоростями перемещения вязких и тепловых волн. Отмечено, что рассматриваемая в работе гиперболизированная система уравнений движения вязкого теплопроводного газа инвариантна не только по отношению к преобразованиям Галилея, но и к повороту, поскольку при дифференцировании по времени компонентов тензора вязких напряжений использована производная Яуманна. Для интегрирования уравнений модели применены гибридный метод Годунова (ГМГ) и многомерный узловой метод характеристик. ГМГ предназначен для интегрирования гиперболических систем, в которых имеются как уравнения, записанные в дивергентном виде, так и уравнения, не приводящиеся к таковому (оригинальный метод Годунова применяется только для систем уравнений, представленных в дивергентной форме). При вычислении потоковых переменных на гранях смежных ячеек использован линеаризованный римановский решатель. Для дивергентных уравнений применена конечно-объемная, а для недивергентных — конечноразностная аппроксимация. Для расчета ряда задач в работе также использовался неконсервативный многомерный узловой метод характеристик, который базируется на расщеплении исходной системы уравнений на ряд одномерных подсистем, для решения которых использован одномерный узловой метод характеристик. С помощью описанных численных методов решен ряд модельных одномерных задач о распаде произвольного разрыва, а также рассчитано двумерное течение вязкого газа при взаимодействии ударного скачка с прямоугольной ступенькой, непроницаемой для газа.

В статье представлены результаты численного моделирования вихревого пространственного нестационарного движения среды, возникающего около боковой и донной поверхностей десантного модуля при его спуске в атмосфере Марса. Численное исследование проведено для высокоскоростного режима обтекания при различных углах атаки. Математическое моделирование осуществлено на основе модели Навье – Стокса и модели равновесных химических реакций для газового состава марсианской атмосферы. Результаты моделирования показали, что при рассматриваемых условиях движения спускаемого аппарата около его боковой и донной поверхностей реализуется нестационарное течение, имеющее ярко выраженный вихревой характер. Численные расчеты указывают на то, что в зависимости от угла атаки нестационарность и вихревой характер потока могут проявляться как на всей боковой и донной поверхностях аппарата, так и, частично, на их подветренной стороне. Для различных углов атаки приводятся картины вихревой структуры потока около поверхности спускаемого аппарата и в его ближнем следе, а также картины полей температуры и показателя адиабаты. Нестационарный характер обтекания подтверждается представленными временными зависимостями газодинамических параметров потока в различных точках поверхности аппарата. Проведенные параметрические расчеты позволили построить зависимости аэродинамических характеристик спускаемого аппарата от угла атаки. Математическое моделирование осуществляется на основе являющегося методом конечных объемов консервативного численного метода потоков, основанного на конечно-разностной записи законов сохранения аддитивных характеристик среды с использованием upwind-аппроксимаций потоковых переменных. Для моделирования возникающей при обтекании сложной вихревой структуры потока около спускаемого аппарата используются неравномерные вычислительные сетки, включающие до 30 миллионов конечных объемов с экспоненциальным сгущением к поверхности, что позволило выявить мелкомасштабные вихревые образования. Численные исследования проведены на базе разработанного комплекса программ, основанного на параллельных алгоритмах используемого численного метода и реализованного на современных многопроцессорных вычислительных системах. Приведенные в статье результаты численного моделирования получены при использовании до двух тысяч вычислительных ядер многопроцессорного комплекса.

Работа посвящена численному моделированию двухфазных течений, а именно расчету сверхзвукового обтекания затупленного тела потоком вязкого газа с примесью относительно крупных частиц, масса которых позволяет после отражения от поверхности выйти за пределы ударного слоя, двигаясь по инерции навстречу набегающему потоку. Натурные и вычислительные эксперименты показывают, что движение высокоинерционных частиц существенным образом изменяет структуру течения газа в ударном слое, а формирующиеся при этом направленные на тело импактные струи вызывают увеличение давления газа вблизи участков поверхности и кратный рост конвективного теплового потока.

Построена математическая модель обтекания затупленного тела сверхзвуковым потоком вязкого газа с твердыми частицами. Решение системы нестационарных уравнений Навье–Стокса в консервативных переменных осуществляется бессеточным методом, в основе которого лежит аппроксимация частных пространственных производных газодинамических величин и содержащих их функций методом наименьших квадратов на множестве распределенных в области расчета узлов. Расчет невязких потоков выполняется методом HLLC в сочетании с MUSCL-реконструкцией третьего порядка, вязких потоков — схемой второго порядка. МНК-аппроксимация частных производных параметров газа по направлению также применяется для реализации краевых условий Неймана на выходной границе области расчета, а также поверхностях обтекаемых тел, которые считаются изотермическими твердыми стенками.

Каждое движущееся тело окружено облаком расчетных узлов, принадлежащих его домену и перемещающихся вместе с ним в пространстве. Реализовано два подхода к моделированию перемещения объектов с учетом обратного влияния на течение газа: метод скользящих облаков фиксированной формы и эволюции единого облака узлов, представляющего собой объединение узлов разных доменов. Проведенные численные эксперименты подтвердили применимость предложенных методов к решению целевых задач моделирования движения крупных частиц в сверхзвуковом потоке.

Выполнена программная реализация представленных алгоритмов на основе технологии параллельных гетерогенных вычислений OpenCL. Представлены результаты моделирования движения крупной частицы вдоль оси симметрии сферы навстречу набегающему потоку с числом Маха $\mathrm{M}=6$.

Рассматривается конечномерная оптимизационная задача, постановка которой, помимо искомых переменных, содержит параметры. Ее решение есть зависимость оптимальных значений переменных от параметров. В общем случае такие зависимости не являются функциями, поскольку могут быть неоднозначными, а в функциональном случае — быть недифференцируемыми. Кроме того, область их существования может оказаться уже области определения функций в условии задачи. Эти свойства затрудняют решение как исходной задачи, так и задач, в постановку которых входят данные зависимости. Для преодоления этих затруднений обычно применяются методы типа недифференцируемой оптимизации.

В статье предлагается альтернативный подход, позволяющий получать решения параметрических задач в форме, лишенной указанных свойств. Показывается, что такие представления могут исследоваться стандартными алгоритмами, основанными на формуле Тейлора. Данная форма есть функция, гладко аппроксимирующая решение исходной задачи. При этом величина погрешности аппроксимации регулируется специальным параметром. Предлагаемые аппроксимации строятся с помощью специальных функций, устанавливающих обратные связи между переменными и множителями Лагранжа. Приводится краткое описание этого метода для линейных задач с последующим обобщением на нелинейный случай.

Построение аппроксимации сводится к отысканию седловой точки модифицированной функции Лагранжа исходной задачи. Показывается, что необходимые условия существования такой седловой точки подобны условиям теоремы Каруша – Куна – Таккера, но не содержат в явном виде ограничений типа неравенств и условий дополняющей нежесткости. Эти необходимые условия аппроксимацию определяют неявным образом. Поэтому для вычисления ее дифференциальных характеристик используется теорема о неявных функциях. Эта же теорема применяется для уменьшения погрешности аппроксимации.

Особенности практической реализации метода функций обратных связей, включая оценки скорости сходимости к точному решению, демонстрируются для нескольких конкретных классов параметрических оптимизационных задач. Конкретно: рассматриваются задачи поиска глобального экстремума функций многих переменных и задачи на кратный экстремум (максимин-минимакс). Также рассмотрены оптимизационные задачи, возникающие при использовании многокритериальных математических моделей. Для каждого из этих классов приводятся демонстрационные примеры.

Численное решение системы уравнений высокотемпературной радиационной газовой динамики (ВРГД) является вычислительно трудоемкой задачей, так как взаимодействие излучения с веществом нелинейно и нелокально. Коэффициенты поглощения излучения зависят от температуры, а поле температур определяется как газодинамическими процессами, так и переносом излучения. Обычно для решения системы ВРГД используется метод расщепления по физическим процессам, выделяется блок решения уравнения переноса совместно с уравнением баланса энергии вещества при известных давлениях и температурах. Построенные ранее разностные схемы, используемые для решения этого блока, обладают порядками сходимости не выше второго. Так как даже на современном уровне развития вычислительной техники имеются ограничения по памяти, то для решения сложных технических задач приходится применять не слишком подробные сетки. Это повышает требования к порядку аппроксимации разностных схем. В данной работе впервые реализованы бикомпактные схемы высокого порядка аппроксимации для алгоритма совместного решения уравнения переноса излучения и уравнения баланса энергии. Предложенный метод может быть применен для решения широкого круга практических задач, так как обладает высокой точностью и подходит для решения задач с разрывами коэффициентов. Нелинейность задачи и использование неявной схемы приводит к итерационному процессу, который может медленно сходиться. В данной работе используется мультипликативный HOLO-алгоритм — метод квазидиффузии В.Я. Гольдина. Ключевая идея HOLO-алгоритмов состоит в совместном решении уравнений высокого порядка (high order, HO) и низкого порядка (low order, LO). Уравнением высокого порядка (HO) является уравнение переноса излучения, которое решается в многогрупповом приближении, далее уравнение осредняется по угловой переменной и получается система уравнений квазидиффузии в многогрупповом приближении (LO₁). Следующим этапом является осреднение по энергии, при этом получается эффективная одногрупповая система уравнений квазидиффузии (LO₂), которая решается совместно с уравнением энергии. Решения, получаемые на каждом этапе HOLO-алгоритма, оказываются тесно связанными, что в итоге приводит к ускорению сходимости итерационного процесса. Для каждого из этапов HOLO-алгоритма предложены разностные схемы, построенные методом прямых в рамках одной ячейки и обладающие четвертым порядком аппроксимации по пространству и третьим порядком по времени. Схемы для уравнения переноса были разработаны Б.В. Роговым и его коллегами, схемы для уравнений LO₁ и LO₂ разработаны авторами. Предложен аналитический тест, на котором демонстрируются заявленные порядки сходимости. Рассматриваются различные варианты постановки граничных условий и исследовано их влияние на порядок сходимости по времени и пространству.

Предложен численный метод решения квазилинейных уравнений параболического типа, основанный на аппроксимации потоков. Описана реализация метода на прямоугольной сетке. Приведены результаты численных расчетов. В отличие от применяемых методов для данного метода используется аппроксимация потоков на нерасширенном шаблоне. Для каждой итерации метода Ньютона возможно решение линейной задачи с помощью метода верхней релаксации (SOR). По сравнению с методами потоковой прогонки рассмотренный метод обладает большим потенциалом для использования на современных параллельных вычислительных комплексах.

Дана постановка задачи управления движения тела в вязкой жидкости. Движение тела индуцируется перемещением внутренних материальных точек. На основе численного решения уравнений движения тела и гидродинамических уравнений получены аппроксимирующие зависимости для вязких сил. С применением аппроксимаций решается задача оптимального управления движением тела по заданной траектории с применением гибридного генетического алгоритма. Установлена возможность направленного движения тела под действием возвратно-поступательного движения внутренней точки. Оптимальное управление направлением движения осуществляется движением другой внутренней точки по круговой траектории с переменной скоростью.

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

В статье сформулирован обобщенный подход к выбору значений структурных параметров искусственной нейронной сети (ИНС) и объема обучающий выборки, основанный на принципе минимизации количества элементов структуры ИНС и объема обучающей выборки при ограничении на значение показателя качества работы нейросетевой модели динамики объекта. Реализован алгоритм выбора структурных параметров ИНС и построения нейросетевой модели.
Проведена серия вычислительных экспериментов, демонстрирующая применимость алгоритма для построения моделей динамических объектов, в основе которых лежит нелинейная автокорреляционная нейронная сеть.