Аннотация:
Для моделирования нерелятивистских колебаний холодной плазмы предложена модификация известной схемы Русанова, имеющей третий порядок точности. Ранее для подобных расчетов в эйлеровых переменных были известны только схемы первого и второго порядков точности. Для тестовой задачи с гладким решением проведено исследование погрешностей построенной схемы, а также – сравнение с погрешностями схемы Мак-Кормака. Для задачи о свободных плазменных колебаниях, инициированных коротким мощным лазерным импульсом, приведены результаты численных экспериментов по сохранению энергии и дополнительной функции для обеих схем, а также – по точности электронной плотности в центре области. Сделан вывод о теоретическом превосходстве схемы Русанова, хотя для практических вычислений более приспособлена схема Мак-Кормака. В первую очередь это касается расчетов “долгоживущих” процессов и колебаний холодной плазмы, близких к реальным. Теоретическое исследование аппроксимации и устойчивости вместе с экспериментальным наблюдением за количественными характеристиками погрешности для наиболее чувствительных величин существенно повышает достоверность вычислений.
Библ. 20. Фиг. 3. Табл. 4.
Ключевые слова:численное моделирование, плазменные колебания, схемы Русанова и Мак-Кормака, порядок точности разностной схемы, законы сохранения, гиперболические системы уравнений.
УДК:519.633.6
Поступила в редакцию: 25.05.2022 Исправленный вариант: 25.05.2022 Принята в печать: 25.05.2022