В настоящее время основной моделью атмосферного электричества является концепция глобальной электрической цепи [1, 2]. Основные физические аспекты для задачи, рассматриваемой в данной работе, приведены в [3, 4].
Рассматривается двумерная модель, в которой эквипотенциальными границами являются поверхность Земли и условная граница, разделяющая атмосферу и ионосферу. Атмосфера является проводящей средой с неоднородной проводимостью. Существует устойчивая разность потенциалов между поверхностью Земли и условной границей, разделяющей атмосферу и ионосферу. Также внутри области протекают токи проводимости. Кроме того, в состав модели входят облака, которые являются независимыми источниками токов. Таким образом, потенциалы поверхности Земли и условной границы, разделяющей атмосферу и ионосферу, токи проводимости в атмосфере и источники тока в облаках образуют распределенную электрическую цепь.
Математическое моделирование для этой задачи требует суперкомпьютерных вычислений и эффективного распараллеливания, так как площадь рассматриваемой области в двумерном случае более 2 миллионов квадратных километров. В нашей работе мы рассматриваем два базовых подхода к численному моделированию скалярного электрического потенциала в атмосфере. Токи и проводимость считаются заданными. Оба подхода основаны на методе конечных элементов, но один из них основан на регулярных сетках и многопоточных вычислениях для итерационного решения сеточных уравнений, а другой, реализованный в пакете FreeFem++, использует нерегулярные адаптивные сетки и может использовать методы разделения областей, MPI. Целью нашей работы является исследование двух численных алгоритмов на основе метода конечных элементов для выявления их преимуществ и недостатков.