Модификация метода крупных частиц до схемы второго порядка точности по пространству и времени для ударно-волновых течений газовзвеси

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