ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЗАДАЧИ АТМОСФЕРНОГО ЭЛЕКТРИЧЕСТВА С НЕИЗВЕСТНЫМ ИОНОСФЕРНЫМ ПОТЕНЦИАЛОМ (2023)
В работе рассматриваются подходы к численному решению задачи о распределении электрического потенциала в рамках двумерной модели атмосферного участка глобальной электрической цепи. Для этой модели формулируется нестандартная стационарная эллиптическая краевая задача с неклассическим граничным условием. Для численного решения этой задачи, с целью изучения возможности и эффективности распараллеливания вычислений, используются два численных алгоритма на основе метода конечных элементов. Приводятся результаты расчетов для модельной задачи, в которой не учитываются особенности рельефа земной поверхности, используется простая модель проводимости и токов.
Идентификаторы и классификаторы
- eLIBRARY ID
- 54630017
В настоящее время основной моделью атмосферного электричества является концепция глобальной электрической цепи [1, 2]. Основные физические аспекты для задачи, рассматриваемой в данной работе, приведены в [3, 4].
Рассматривается двумерная модель, в которой эквипотенциальными границами являются поверхность Земли и условная граница, разделяющая атмосферу и ионосферу. Атмосфера является проводящей средой с неоднородной проводимостью. Существует устойчивая разность потенциалов между поверхностью Земли и условной границей, разделяющей атмосферу и ионосферу. Также внутри области протекают токи проводимости. Кроме того, в состав модели входят облака, которые являются независимыми источниками токов. Таким образом, потенциалы поверхности Земли и условной границы, разделяющей атмосферу и ионосферу, токи проводимости в атмосфере и источники тока в облаках образуют распределенную электрическую цепь.
Математическое моделирование для этой задачи требует суперкомпьютерных вычислений и эффективного распараллеливания, так как площадь рассматриваемой области в двумерном случае более 2 миллионов квадратных километров. В нашей работе мы рассматриваем два базовых подхода к численному моделированию скалярного электрического потенциала в атмосфере. Токи и проводимость считаются заданными. Оба подхода основаны на методе конечных элементов, но один из них основан на регулярных сетках и многопоточных вычислениях для итерационного решения сеточных уравнений, а другой, реализованный в пакете FreeFem++, использует нерегулярные адаптивные сетки и может использовать методы разделения областей, MPI. Целью нашей работы является исследование двух численных алгоритмов на основе метода конечных элементов для выявления их преимуществ и недостатков.
Список литературы
- F. Leblanc, K. L. Aplin, Y. Yair, et al. (Eds.), Planetary Atmospheric Electricity (Springer, New York, 2008). DOI: 10.1007/978-0-387-87664-1
- C. T. R. Wilson, “Investigations on Lighting Discharges and on the Electric Field of Thunderstorms”, Philos. Trans. Roy. Soc. Lon. A. 221, 73-115 (1921). DOI: 10.1098/rsta.1921.0003
- S. V. Anisimov, S. S. Bakastov, and E. A. Mareev, “Spatiotemporal Structures of Electric Field and Space Charge in the Surface Atmospheric Layer”, J. Geophys. Res. Atmos. 99 (D5), 10603-10610 (1994). DOI: 10.1029/93JD03519
- A. V. Kalinin and N. N. Slyunyaev, “Initial-Boundary Value Problems for the Equations of the Global Atmospheric Electric Circuit”, J. Math. Anal. App. 450 (1), 112-136 (2017). DOI: 10.1016/j.jmaa.2017.01.025 EDN: YVAIYH
- N. N. Slyunyaev, E. A. Mareev, A. V. Kalinin, and A. A. Zhidkov, “Influence of Large-Scale Conductivity Inhomogeneities in the Atmosphere on the Global Electric Circuit”, J. Atmos. Sci. 71 (11), 4382-4396 (2014). DOI: 10.1175/JAS-D-14-0001.1 EDN: UFHKYD
- A. A. Zhidkov and A. V. Kalinin, “Correctness of One Mathematical Problem of Atmospheric Electricity”, Vestn. Lobachevskii Univ. Nizhni Novgorod, No. 4, 123-129 (2009).
- J.-M. Jin, The Finite Element Method in Electromagnetics (Wiley, Hoboken, 2014).
- Eigenhttps://eigen.tuxfamily.org/index.php Cited August 2, 2023.
Выпуск
Методы и алгоритмы вычислительной математики и их приложения.
Параллельные программные средства и технологии.
Другие статьи выпуска
В настоящее время обнаружение аномалий в длинных временных рядах возникает в широком спектре предметных областей: цифровая индустрия, здравоохранение, моделирование климата, финансовая аналитика и др. Диссонанс формализует понятие аномалии и определяется как подпоследовательность ряда, которая имеет расстояние до своего ближайшего соседа, не превышающее наперед заданного аналитиком порога. Ближайшим соседом подпоследовательности является та подпоследовательность ряда, которая не пересекается с данной и имеет минимальное расстояние до нее. В статье представлен новый алгоритм поиска диссонансов временн´ого ряда на вычислительном кластере, каждый узел которого оснащен графическим процессором. Алгоритм применяет параллелизм по данным: временн´ой ряд разбивается на непересекающиеся фрагменты, обрабатываемые графическими процессорами узлов вычислительного кластера. С помощью ранее разработанного авторами параллельного алгоритма на каждом узле выполняется отбор локальных кандидатов в диссонансы. Далее с помощью обменов на каждом узле формируется множество глобальных кандидатов как объединение всех локальных кандидатов. Затем каждый узел выполняет глобальную очистку, удаляя из множества глобальных кандидатов ложноположительные диссонансы. Глобальная очистка распараллеливается на основе блочного умножения матрицы кандидатов и матрицы подпоследовательностей фрагмента. Результирующее множество диссонансов формируется как пересечение множеств, полученных узлами по итогу глобальной очистки. Вычислительные эксперименты с синтетическими и реальными временными рядами, проведенные на платформе суперкомпьютеров Ломоносов-2 и Лобачевский, оснащенных 48-64 графическими процессорами, показывают высокую масштабируемость разработанного алгоритма.
Статья посвящена разработке эффективных численных методов решения прямых задач распространения волн в твердых телах в векторных математических моделях. Итерационные методы решения обратных задач волновой томографии используют на каждой итерации решение прямой задачи распространения волн как в прямом, так и в обратном времени для вычисления градиента функционала невязки. Поэтому решение прямой задачи распространения волн в упругих средах является неотъемлемой частью решения обратных задач волновой томографии. Целью статьи также является определение с помощью методов математического моделирования характеристик волн Лэмба для ультразвуковой диагностики дефектов в тонких пластинах, определение диапазонов значений характерных параметров эксперимента по томографической диагностике в тонких пластинах на волнах Лэмба. Инструментом для проведения математического моделирования являются разрабатываемые численные методы и программы решения прямых задач. Конечной целью исследований является разработка методов решения обратных задач томографического неразрушающего ультразвукового контроля как на волнах Лэмба, так и на объемных волнах.
Исследуется задача выразимости всех функций x1(t), x2(t), … , xn(t), входящих в заданную однородную систему линейных дифференциальных уравнений с постоянными коэффициентами x′(t) = A·x(t), в виде линейных комбинаций производных только одной неизвестной функции xк(t), входящей в эту систему. Найден простой критерий выразимости всех функций системы x′(t) = A·x(t) в виде линейных комбинаций производных xк(t) и доказана его корректность. На основе доказанного критерия разработан соответствующий алгоритм и обоснована его корректность.
В настоящее время обработка данных временных рядов осуществляется в широком спектре научных и практических приложений, в которых актуальной является задача восстановления единичных точек или блоков значений временного ряда, пропущенных из-за аппаратных или программных сбоев либо ввиду человеческого фактора. В статье представлен метод SANNI (Snippet and Artificial Neural Network-based Imputation) для восстановления пропущенных значений временного ряда, обрабатываемого в режиме офлайн. SANNI включает в себя две нейросетевые модели: Распознаватель и Реконструктор. Распознаватель определяет сниппет (типичную подпоследовательность) ряда, на который наиболее похожа данная подпоследовательность с пропущенной точкой, и состоит из следующих трех групп слоев: сверточные, рекуррентный и полносвязные. Реконструктор, используя выход Распознавателя и входную подпоследовательность c пропуском, восстанавливает пропущенную точку. Реконструктор состоит из трех групп слоев: сверточные, рекуррентные и полносвязные. Топологии слоев Распознавателя и Реконструктора параметризуются относительно соответственно количества сниппетов и длины сниппета. Представлены методы подготовки обучающих выборок указанных нейросетевых моделей. Проведены вычислительные эксперименты, показавшие, что среди передовых аналитических и нейросетевых методов SANNI входит в тройку лучших.
Издательство
- Издательство
- МГУ
- Регион
- Россия, Москва
- Почтовый адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- Юр. адрес
- оссийская Федерация, 119991, Москва, Ленинские горы, д. 1
- ФИО
- Садовничий Виктор Антонович (РЕКТОР)
- E-mail адрес
- info@rector.msu.ru
- Контактный телефон
- +7 (495) 9391000
- Сайт
- https://msu.ru/