Computer Simulation of the Propagation of a Plane Electromagnetic Wave in a Waveguide Formed by the Earth’s Surface and the Ionosphere Under the Condition of Inhomogeneous Boundary Conductivity
Authors: Твёрдый Д.А. , Малкин Е.И. 
Published: November 2023

В статье приводятся результаты компьютерного моделирования распространения плоской электромагнитной волны. Грозовые разряды являются естественными источниками импульсного электромагнитного излучения (атмосферика). Атмосферик распространяется в волноводе Земля-Ионосфера как плоская электромагнитная волна с широким спектром частот с максимумом интенсивности в диапазоне 4-10 кГц. После землетрясений с магнитудами порядка 7 и более, увеличивается насыщенность подземных вод вторичными минералами, что приводит к локальному увеличению проводимости земли в данных областях. Что определяет электрофизические свойства земли, как параметры нижней границы волновода Земля-Ионосфера. Которые влияют на характеристики распространяющихся в волноводе электромагнитных волн. Предполагается, что, изучая параметры атмосферика можно установить наличие неоднородности проводимости стенки волновода. На основе системы уравнений Максвелла с краевыми условиями в виде Perfectly matched layer, задаётся математическая модель процесса. Краевые условия модели определяют область распространения электромагнитного широкополосного сигнала, как волновод с неоднородной проводимостью границ. Система модельных уравнений решается численным методом Finite-Difference Time-domain. Для решения поставленной задачи и проведения компьютерного моделирования разработан программный комплекс в среде MATLAB. С целью проверки предположения, был проведён ряд компьютерных симуляций. В результате показано, что имеет место обратное рассеяние электромагнитной волны на трассе волновода, возникающее в следствии отражения волны при её взаимодействии с неоднородностью проводимости нижней границы волновода. Показано, что с помощью математического моделирования процесса распространения атмосферика и его взаимодействий с неоднородностью в волноводе можно установить наличие неоднородности и её связь с характеристиками излучения.

The article presents the results of computer simulation of the propagation of a plane electromagnetic wave. Lightning discharges are natural sources of pulsed electromagnetic radiation (atmosphere). Atmospheric propagates in the Earth-Ionosphere waveguide as a plane electromagnetic wave with a wide frequency spectrum with an intensity maximum in the range of 4-10 kHz. After earthquakes with magnitudes of the order of 7 or more, the saturation of groundwater with secondary minerals increases, which leads to a local increase in the conductivity of the earth in these areas. What determines the electrophysical properties of the earth, as the parameters of the lower boundary of the Earth-Ionosphere waveguide. Which affect the characteristics of electromagnetic waves propagating in the waveguide. It is assumed that, by studying the parameters of the atmosferic, it is possible to establish the presence of an inhomogeneity in the conductivity of the waveguide wall. Based on the system of Maxwell equations with boundary conditions in the form of a Perfectly matched layer, a mathematical model of the process is set. The boundary conditions of the model determine the region of propagation of an electromagnetic broadband signal as a waveguide with inhomogeneous boundary conductivity. The system of model equations is solved by the numerical Finite-Difference Time-domain method. In order to solve the problem and conduct computer modeling, a software package was developed in MATLAB environment. In order to verify the assumption, a number of computer simulations were carried out. As a result, it is shown that there is a backscattering of electromagnetic wave on the waveguide trace, arising as a consequence of reflection of the wave in its interaction with the inhomogeneity of the conductivity of the lower boundary of the waveguide. It is shown that with the help of mathematical modeling of the process of atmospheric propagation and its interactions with inhomogeneity in the waveguide it is possible to establish the presence of inhomogeneity and its relationship to the characteristics of radiation.