ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ И ПРОГРАММИРОВАНИЕ

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

The objective of this paper is to construct a generalized quadratic spectrum approximation based on the Kantorovich projection method which llows us to deal with the spectral pollution problem. For this purpose, we prove that the property U (see Eq. 3) holds under weaker conditions than the norm and the collectively compact convergence. Numerical results illustrate the effectiveness and the convergence of our method.

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

В работе рассматривается предобусловливатель блочного неполного обратного треугольного разложения первого порядка “по значению” BIIC-IC1 для решения систем линейных алгебраических уравнений с симметричной положительно определенной матрицей. Рассматривается способ применения MPI+OpenMP технологии для построения и обращения предобусловливателя BIIC-IC1, при этом в предобусловливателе число блоков кратно числам используемых процессоров и используемых потоков. Предлагается способ применения MPI+OpenMP технологии для построения и обращения предобусловливателя BIIC-IC1, в котором для применения OpenMP технологии используется специальное упорядочение узлов сетки внутри подобластей, соответствующих расчетам на процессорах. Проводится сравнение времени решения задач методом сопряженных градиентов с предобусловливателем BIIC-IC1 с использованием MPI и гибридной MPI+OpenMP технологии на примере модельной задачи и ряда задач из коллекции разреженных матриц SuiteSparse.

Разработан алгоритм высокоточного численного решения эллиптического уравнения второго порядка при наличии в области нескольких интерфейсов, в том числе пересекающихся и невыпуклых. Для аппроксимации задачи в окрестности интерфейсов используются нерегулярные ячейки (н-ячейки), отсекаемые ими от регулярных ячеек прямоугольной сетки, и законтурные части этих ячеек. Для построения приближенного решения предложено: 1) выписывать дополнительные условия согласования в н-ячейках на интерфейсах, увеличивая количество согласуемых ячеек вблизи интерфейсов; 2) уменьшать общую часть интерфейса, заключенную в соседних ячейках и используемую для записи условий. Для решения краевой задачи Дирихле реализован hp-вариант метода коллокации и наименьших квадратов (hp-МКНК) в сочетании с современными алгоритмами ускорения итерационного процесса: предобуславливание; распараллеливание с помощью OpenMP; ускорение, основанное на подпространствах Крылова; многосеточный алгоритм. При решении различных тестовых задач исследованы сходимость hp-МКНК и обусловленность возникающих переопределенных систем линейных алгебраических уравнений (СЛАУ). Проведено сравнение результатов, полученных МКНК, с результатами других авторов, использовавших метод MIB (англ. matched interface and boundary).