Архив статей журнала
В рамках статистического подхода, основанного на кинетическом уравнении для функции плотности вероятности распределения скорости и температуры частиц, построена континуальная модель, описывающая псевдотурбулентные течения дисперсной фазы. Введение функции плотности вероятности позволяет получить статистическое описание ансамбля частиц вместо динамического описания отдельных частиц на основе уравнений движения и теплопереноса типа Ланжевена. На основе уравнений для первых и вторых моментов дисперсной фазы проводится численное моделирование нестационарного течения газовзвеси, возникающего при взаимодействии ударной волны с облаком частиц. Основные уравнения имеют гиперболический тип, записываются в консервативной форме и решаются с использованием численного метода типа Годунова повышенного порядка точности. Обсуждается влияние двумерных эффектов на формирование ударно-волновой структуры течения и пространственно-временны´е зависимости концентрации частиц и других параметров потока.
Для моделирования колебаний холодной плазмы как в нерелятивистском случае, так и с учетом релятивизма предложены модификации классических разностных схем второго порядка точности: метода МакКормака и двухэтапного метода Лакса-Вендроффа. Ранее для подобных расчетов в эйлеровых переменных была известна только схема первого порядка точности. Для задачи о свободных плазменных колебаниях, инициированных коротким мощным лазерным импульсом, с целью тестирования представленных схем проведены численные эксперименты по сохранению энергии и других величин. Сделан вывод о достоверности численного анализа колебаний как на основе схемы МакКормака, так и на основе схемы Лакса-Вендроффа, однако для расчетов “долгоживущих” процессов первая схема более предпочтительна. Теоретическое исследование аппроксимации и устойчивости вместе с экспериментальным наблюдением за количественными характеристиками погрешности для наиболее чувствительных величин существенно повышает достоверность вычислений.
На основе модели взаимопроникающих континуумов проводится численное моделирование нестационарного течения газовзвеси, возникающего при взаимодействии ударной волны со слоем инертных частиц. Каждая фаза описывается набором уравнений, выражающих законы сохранения массы, импульса и энергии. Межфазное взаимодействие учитывается при помощи источниковых членов в уравнениях изменения количества движения и энергии. Основные уравнения для газовой и дисперсной фаз имеют гиперболический тип, допускают запись в консервативной форме и решаются с использованием численного метода типа Годунова повышенного порядка точности. Для дискретизации уравнений по времени применяется метод Рунге-Кутты 3-го порядка. Построенная модель позволяет рассчитывать широкий спектр режимов течения газовзвеси, возникающих при изменении объемной концентрации дисперсной фазы. Обсуждаются вопросы, связанные с замыканием математической модели, а также детали реализации численной модели. Приводятся ударно-волновая структура течения и пространственно-временные зависимости концентрации частиц и других параметров потока.