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

В работе демонстрируется, как метод апостериорной оценки порядка точности разностной схемы по Ричардсону позволяет сделать вывод о некорректности постановки (в смысле отсутствия решения) решаемой численно начально-краевой задачи для уравнения в частных производных. Это актуально в ситуации, когда аналитическое доказательство некорректности постановки ещё не получено или принципиально невозможно.

Предложен балансно-характеристический метод решения систем линейных дифференциальных уравнений в частных производных гиперболического типа, обладающий четвертым порядком аппроксимации на равномерных сетках и вторым порядком и улучшенными дисперсионными свойствами на неравномерных сетках. Метод основан на известной схеме КАБАРЕ, балансные фазы которой модифицированы путем добавления антидисперсионных членов особого вида. Ранее метод, обладающий схожими свойствами, предлагался только для простейшего одномерного линейного уравнения переноса. Приведенная модификация схемы позволяет улучшить дисперсионные свойства переноса сразу всех инвариантов Римана рассматриваемой системы уравнений. Схема бездиссипативна при отключенных процедурах монотонизации и устойчива при числах Куранта CFL ≤ 1. Точность метода и его порядок сходимости продемонстрированы на серии расчетов задачи о переносе волны, промодулированной гауссианом, на последовательности сгущающихся сеток. Предложенный метод планируется использовать в качестве основы для построения схемы КАБАРЕ с улучшенными дисперсионными свойствами для систем нелинейных дифференциальных уравнений.

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

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

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