Труды КНЦ вып.5. ГЕЛИОГЕОФИЗИКА. 8/2019(10)

На подшаге затухания и вращения сигнала магнитное поле сигнала не изменяется, и учитываются только затухание электрического поля и его вращение за счет внешнего геомагнитного поля. При этом в каждой точке расчетной сетки аналитически интегрируется система уравнений: Эта система разделяется на две независимых системы для продольных Системы (4) и (5) являются автономными линейными системами ОДУ с постоянными коэффициентами. Решения задачи Коши для этих систем выражаются достаточно громоздкими аналитическими формулами. В той части ионосферы, где выполнено неравенство Ѵе Т > 20 , системы (4) и (5) заменяются более простые на уравнение и систему которые также имеют аналитические решения задачи Коши. Изложенная схема позволяет учесть частотную дисперсию и моделировать распространение широкополосных сигналов произвольной формы. Среда распространения, модель, источник В качестве среды распространения для численного эксперимента авторы использовали участок волновода Земля-ионосфера с параметрами концентрации электронов и частотой столкновения электронов, полученных на основе данных модели IRI2016 за 26 декабря 2014 в 00:00 UT. Центр участка распологался в городе в точке 67°34'03" с. ш. 33°23'36" в.д., а его горизонтальный размер составлял 1400x384 км. По высоте в атмосфере и ионосфере - 250 км, в глубину в литосфере - 50 км. Шаг сетки по горизонтали составлял 2 км. Шаг сетки по вертикали над поверхностью Земли - 1км, по вертикали в литосфере 500 м. Центр источника сигнала располагался на расстоянии 192 км от трех боковых граней полученного параллелепипеда. Для всех внешних граней кроме нижней использовалось условие свободного ухода волны совместно с адаптированными поглощающими слоями и профилем потерь, как в работе [6]. Профиль проводимости литосферы был аппроксимирован на основе ( 3 ) электрического поля Е , = ( Ь , Е ) и тока 7ц = ( b j ) (4) и поперечных электрического поля Е ± = Е - Ъ Е ^ и тока j ± = j - b j \ \ 195