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