ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ И ПРОГРАММИРОВАНИЕ
Архив статей журнала
В статье рассматривается явно-неявная балансно-характеристическая схема CABARETI-NH (CABARET Implicit Non-Hydrostatic), основанная на схеме КАБАРЕ, для решения гиперболизированной системы уравнений Навье-Стокса. Неявность вдоль одного пространственного направления позволяет значительно увеличить шаг по времени на вычислительных сетках с большим аспектным отношением ячеек. Для разрешения введенной неявности используется метод гиперболической прогонки. Это позволяет сохранить вычислительную эффективность алгоритма на уровне явных схем. Приводятся результаты валидации модели на лабораторном эксперименте трехмерного гравитационного течения стратифицированной жидкости.
Схема КАБАРЕ, являющаяся представителем семейства балансно-характеристических методов, широко используется при решении многих задач для систем дифференциальных уравнений гиперболического типа в эйлеровых переменных. Возрастающая актуальность задач взаимодействия деформируемых тел с потоками жидкости и газа требует адаптации этого метода на лагранжевы и смешанные эйлерово-лагранжевы переменные. Ранее схема КАБАРЕ была построена для одномерных уравнений газовой динамики в массовых лагранжевых переменных, а также для трехмерных уравнений динамической упругости. В первом случае построенную схему не удалось обобщить на многомерные задачи, а во втором - использовался необратимый по времени алгоритм передвижения сетки. В данной работе представлено обобщение метода КАБАРЕ на двумерные уравнения газовой динамики и динамической упругости в смешанных эйлерово-лагранжевых и лагранжевых переменных. Построенный метод является явным, легко масштабируемым и обладает свойством временной обратимости. Метод тестируется на различных одномерных и двумерных задачах для обеих систем уравнений (соударение упругих тел, поперечные колебания упругой балки, движение свободной границы идеального газа).
Предложен балансно-характеристический метод решения систем линейных дифференциальных уравнений в частных производных гиперболического типа, обладающий четвертым порядком аппроксимации на равномерных сетках и вторым порядком и улучшенными дисперсионными свойствами на неравномерных сетках. Метод основан на известной схеме КАБАРЕ, балансные фазы которой модифицированы путем добавления антидисперсионных членов особого вида. Ранее метод, обладающий схожими свойствами, предлагался только для простейшего одномерного линейного уравнения переноса. Приведенная модификация схемы позволяет улучшить дисперсионные свойства переноса сразу всех инвариантов Римана рассматриваемой системы уравнений. Схема бездиссипативна при отключенных процедурах монотонизации и устойчива при числах Куранта CFL ≤ 1. Точность метода и его порядок сходимости продемонстрированы на серии расчетов задачи о переносе волны, промодулированной гауссианом, на последовательности сгущающихся сеток. Предложенный метод планируется использовать в качестве основы для построения схемы КАБАРЕ с улучшенными дисперсионными свойствами для систем нелинейных дифференциальных уравнений.