О ЧИСЛЕННОМ МОДЕЛИРОВАНИИ МЕДЛЕННОЙ НЕОБЫКНОВЕННОЙ ВОЛНЫ В МАГНИТОАКТИВНОЙ ПЛАЗМЕ (2020)
Численно и аналитически исследовано влияние внешнего магнитного поля на плоские нерелятивистские нелинейные плазменные колебания. Для инициализации медленной необыкновенной волны в магнитоактивной плазме предложен способ построения недостающих начальных условий на основе решения линейной задачи методом Фурье. С целью численного моделирования нелинейной волны построена схема метода конечных разностей второго порядка точности типа МакКормака на основе эйлеровых переменных. Показано, что при учете внешнего магнитного поля ленгмюровские колебания трансформируются в медленную необыкновенную волну, энергия которой вибрирует при перемещении от начала координат. При этом скорость волны увеличивается с ростом внешнего постоянного поля, что способствует выносу энергии из первоначальной области локализации колебаний.
Идентификаторы и классификаторы
Полностью ионизованная плазма является сильно нелинейной средой, в которой даже относительно небольшие начальные коллективные смещения частиц могут приводить к колебаниям и волнам большой амплитуды [1]. Численному моделированию колебаний в холодной плазме, а также кильватерных волн, возбуждаемых коротким мощным лазерным импульсом, посвящена монография [2]. Сложность и разнообразие постановок задач резко возрастает, если рассматривать динамику магнитоактивной плазмы, т.е. плазмы, помещенной во внешнее магнитное поле. Даже в случае малых возмущений, т.е. при рассмотрении только линейных волн, в классических учебниках и монографиях вводится их специальная классификация [3, 4]. В частности, выделяют волны, распространяющиеся в продольном и поперечном направлениях по отношению к вектору заданного магнитного поля. Уже в холодной магнитоактивной плазме существует пять ветвей колебаний, распространяющихся строго вдоль внешнего магнитного поля, причем четыре из них описывают поперечные, а одна — продольные волны. При этом следует отметить, что в общем случае линейные продольные и поперечные волны не являются независимыми, как в изотропной плазме. Именно по этой причине особый интерес вызывает моделирование нелинейных колебаний и волн в магнитоактивной плазме.
Список литературы
- Davidson R.C. Methods in nonlinear plasma theory. New York: Academic Press, 1972.
- Чижонков Е.В. Математические аспекты моделирования колебаний и кильватерных волн в плазме. М.: Физматлит, 2018.
- Александров А.Ф., Богданкевич Л.С., Рухадзе А.А. Основы электродинамики плазмы. М.: Высшая школа, 1988.
- Гинзбург В.Л., Рухадзе А.А. Волны в магнитоактивной плазме. М.: Наука, 1975.
- Красовицкий В.Б., Бугримов А.Л. Ускорение электронов в поле световой волны в магнитоактивной плазме // Вестник Московского государственного областного университета. Серия: Физика-математика. 2016. № 1. 31-43. EDN: VXKOEH
- Moradi A. Energy behaviour of extraordinary waves in magnetized quantum plasmas // Physics of Plasmas. 2018. 25. DOI: 10.1063/1.5031753
- Силин В.П. Введение в кинетическую теорию газов. М.: Наука, 1971.
- Силин В.П., Рухадзе А.А. Электромагнитные свойства плазмы и плазмоподобных сред. М.: ЛИБРОКОМ, 2012. EDN: QJZTUZ
- Фролов А.А., Чижонков Е.В. О применении закона сохранения энергии в модели холодной плазмы // Журнал вычислит. математики и матем. физики. 2020. 60, № 3. 503-519. EDN: LEBHFI
-
Sheppard C.J. R. Cylindrical lenses - focusing and imaging: a review [Invited] // Applied Optics. 2013. 52, N 4. 538-545.
-
Chizhonkov E.V., Frolov A.A. Influence of electron temperature on breaking of plasma oscillations // Russ. J. Numer. Anal. Math. Modelling. 2019. 34, N 2. 71-84. EDN: WJKHMZ
-
MacCormack R.W. The effect of viscosity in hypervelocity impact cratering // J. Spacecr. Rockets. 2003. 40, N 5. DOI: 10.2514/2.6901
-
Поттер Д. Вычислительные методы в физике. М.: Мир, 1975.
-
Шокин Ю.И., Яненко Н.Н. Метод дифференциального приближения. Применение к газовой динамике. Новосибирск: Наука, 1985.
-
Андерсон Д., Таннехилл Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. Т. 1. М.: Мир, 1990.
-
Федотова З.И. О применении разностной схемы Мак-Кормака для задач длинноволновой гидродинамики // Вычислительные технологии. 2006. 11, спецвыпуск. Часть II. 53-63. EDN: MWKRUN
-
Чижонков Е.В. О схемах второго порядка точности для моделирования плазменных колебаний // Вычислительные методы и программирование. 2020. 21 (1). 115-128. EDN: OZLVQA
-
Фролов А.А., Чижонков Е.В. Численное моделирование плазменных колебаний с учетом теплового движения электронов // Вычислительные методы и программирование. 2018. 19 (2). 194-206. EDN: YRJNHV
-
Maity C. Lagrangian Fluid Technique to Study Nonlinear Plasma Dynamics. PhD Thesis. Kolkata: Saha Institute of Nuclear Physics, 2013.
-
Swarztrauber P.N. Vectorizing the FFTs // Parallel Computations. New York: Academic Press. 1982. 51-83.
Выпуск
Методы и алгоритмы вычислительной математики и их приложения.
Параллельные программные средства и технологии.
Другие статьи выпуска
Статья посвящена изучению нелинейных эффектов в динамике популяции промысловой рыбы пиленгас Азовского моря при низкой и высокой ее численности с учетом эффекта Олли, конкуренции за ресурсы, таксиса, вылова, пространственного распределения биогенных веществ и детрита на основе многовидовой модели взаимодействия планктона и рыб. Дискретный аналог разработанной модельной задачи водной экологии, входящей в состав программного комплекса, получен на основе схем второго порядка точности с учетом частичной заполненности расчетных ячеек. Возникающая в процессе дискретизации система сеточных уравнений большой размерности была решена на основе модифицированного попеременно-треугольного метода, имеющего наибольшую скорость сходимости при условии асимптотической устойчивости разностных схем для параболических уравнений, эффективность которого была улучшена на основе уточненных спектральных оценок. Разработка эффективных параллельных алгоритмов численной реализации поставленной задачи биологической кинетики, ориентированных на многопроцессорную вычислительную систему (МВС) и графический ускоритель NVIDIA Tesla K80 с модификацией формата хранения данных, позволила анализировать процессы воспроизводства популяций биогидроценоза в режиме реального и ускоренного времени.
Для изучения одного из важнейших процессов нефтепереработки - каталитического риформинга, требуется детализированная кинетическая модель. При разработке кинетической модели возникает сложность в связи с большим количеством компонентов реакционной смеси и большим количеством стадий химических превращений. Альтернативой могут быть сокращенные механизмы реакций, которые применимы для решения задачи и обеспечивают реалистичное описание процесса. В данной работе для анализа кинетической модели и получения сокращенного механизма реакции используются методы анализа чувствительности математической модели. Применение указанной методики позволяет выявить стадии каталитического риформинга бензина, наименее влияющие на общую динамику изменения концентраций значимых веществ реакции. Исследовано влияние исключения данных стадий на кинетику процесса с химической точки зрения. Предложена редуцированная схема каталитического риформинга бензина с исключением данных стадий. Редуцированная схема обеспечивает вполне удовлетворительное согласие как по профилям температуры, так и по профилям концентраций значимых веществ реакции.
Рассматривается задача численного моделирования распространения электромагнитных волн в неоднородной тропосфере на основе широкоугольных обобщений метода параболического уравнения. Используется конечно-разностная аппроксимация Паде оператора распространения. Существенно, что в предлагаемом подходе указанная аппроксимация осуществляется одновременно по продольной и поперечной координатам. При этом допускается моделирование произвольного коэффициента преломления тропосферы. Метод не накладывает ограничений на максимальный угол распространения. Для различных условий распространения радиоволн проведено сравнение с методом расщепления Фурье и методом геометрической теории дифракции. Показаны преимущества предлагаемого подхода.
Одна из основных сложностей разработки параллельной программы для кластера - необходимость принятия глобальных решений по распределению данных и вычислений с учетом свойств всей программы, а затем выполнения кропотливой работы по модификации программы и ее отладки. Большой объем программного кода, а также многомодульность, многовариантность и многоязыковость, затрудняют принятие решений по согласованному распределению данных и вычислений. Опыт использования предыдущей системы САПФОР показал, что при распараллеливании на кластер больших программ и программных комплексов необходимо уметь распараллеливать их постепенно, начиная с наиболее времяемких фрагментов и постепенно добавляя новые фрагменты, пока не достигнем желаемого уровня эффективности параллельной программы. С этой целью предыдущая система была полностью переработана, и на ее основе была создана новая система SAPFOR (System FOR Automated Parallelization). В данной статье будет рассмотрен опыт применения метода частичного распараллеливания, идея которого заключается в том, что распараллеливанию подвергается не вся программа целиком, а ее части (области распараллеливания) - в них заводятся дополнительные экземпляры требуемых данных, производится распределение этих данных и соответствующих им вычислений.
Существующая технология численного анализа устойчивости течений вязкой несжимаемой жидкости в каналах постоянного сечения была ранее расширена на случай локальных пространственных аппроксимаций на неструктурированных сетках, приводящих к задачам с большими разреженными матрицами. Для пространственной аппроксимации при этом используется метод конечных элементов, а для решения частичных проблем собственных значений, возникающих при исследовании устойчивости течений, эффективный метод ньютоновского типа. В данной работе проводится подробное численное исследование предложенного подхода на примере двумерной конфигурации - течения Пуазейля в канале эллиптического сечения. Работоспособность подхода демонстрируется для широкого диапазона отношений длин полуосей сечения вплоть до отношения, при котором данное течение становится линейно неустойчивым. Показана сходимость ведущей части спектра по шагу сетки и совпадение результатов с результатами, полученными на основе аппроксимации спектральным методом коллокаций.
Рассматриваются априорные оценки неоднозначности (погрешности) приближенных решений условно-корректных нелинейных обратных задач, основанные на модуле непрерывности обратного оператора и его модификациях. Установлена связь модуля непрерывности обратного оператора с разрешающей способностью геофизического метода. Показано, что в классе кусочно-постоянных решений, определенных на заданной сетке параметризации, модуль непрерывности обратного оператора и его модификации монотонно возрастают с увеличением размерности сетки. Предложен метод построения оптимальной сетки параметризации, которая имеет максимальную размерность при условии, что модуль непрерывности обратного оператора не превышает заданной величины. Представлен численный алгоритм расчета модуля непрерывности обратного оператора и его модификаций с использованием алгоритмов Монте-Карло, исследуются вопросы сходимости алгоритма. Предлагаемый метод применим также для расчета классических апостериорных оценок погрешности. Приводятся численные примеры для нелинейных обратных задач геоэлектрики.
Для генерации сверхподробных тетраэдральных сеток объемом до 1 миллиарда ячеек используется открытое ПО Gmsh. Пакет позволяет строить автоматически указанные сетки на ПК в параллельном режиме OpenMP за время, не превышающее 1 часа при использовании ПК Intel i7-9700K. Описан опыт применения пакета для построения пространственных сеток для задач обтекания аэрокосмических объектов. Приведены примеры суперкомпьютерных расчетов обтекания аэрокосмических объектов.
Издательство
- Издательство
- МГУ
- Регион
- Россия, Москва
- Почтовый адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- Юр. адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- ФИО
- Садовничий Виктор Антонович (РЕКТОР)
- E-mail адрес
- info@rector.msu.ru
- Контактный телефон
- +7 (495) 9391000
- Сайт
- https://msu.ru/