ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ПРОСТРАНСТВЕННО-ВРЕМЕННЫХ ЧАСТОТНО-МОДУЛИРОВАННЫХ РАДИОВОЛН В АНИЗОТРОПНОЙ СРЕДЕ
Крюковский Андрей Сергеевич,
д.ф-м.н, профессор, декан факультета ИСиКТ,
НОУ ВПО Российский новый университет, Москва, Россия,
Лукин Дмитрий Сергеевич,
д.ф-м.н, профессор кафедры ТиИСиСС,
НОУ ВПО Российский новый университет, Москва, Россия,
Растягаев Дмитрий Владимирович,
к.ф-м.н, доцент, проректор по ИТ,
НОУ ВПО Российский новый университет, Москва, Россия, [email protected]
Скворцова Юлия Игоревна,
заместитель декана факультета ИСиКТ, НОУ ВПО Российский новый университет, Москва, Россия, [email protected]
Работа выполнена при поддержке РФФИ (гранты № 15-02-04206-а, № 13-07-00937-а, ОФИ_М № 13-02-12121)
Ключевые слова: обыкновенная и необыкновенная волна; численное моделирование; распространение частотно-модулированных сигналов; анизотропная ионосфера; бихарактеристическая система; волнообразные возмущения; пространственно-временные каустики; доплеровское изменение частоты.
В случае обыкновенной и необыкновенной волны выполнено численное моделирование двухскач-кового распространения радиосигналов в анизотропной ионосферной плазме с частотной модуляцией. Активное применение радиоволн декаметрового диапазона для обеспечения дальней радиосвязи, радиолокации, радионавигации, загоризонтного радиозондирования, а также изучения структуры верхней атмосферы Земли - ионосферы определяет актуальность исследования. Несмотря на значительный объем работ направленных на изучения распространения частотно-модулированного радиоизлучения, распространение таких сигналов в анизотропной ионосферной плазме изучено ещё недостаточно. Целью работы является восполнение данного пробела. Основой исследования является численное решение пространственно-временной бихарактеристической системы Гамильтона, неизвестными в которой являются компоненты волнового вектора, координаты луча, а также частота и время. Считается, что источник излучения точечный, находится вне магнитоактивной плазмы, частотная модуляция линейная, а амплитуда магнитного поля постоянная. Рассмотрена модель двухслойной ионосферы как при наличии, так и в случае отсутствия нестационарного волнообразного возмущения, моделирующего гравитационные волны в плазме. Известно, что внешние магнитное поле приводит к выходу лучевой траектории из плоскости распространения. Сопоставлены проекции лучевых структур обыкновенной и необыкновенной волны на вертикальную плоскость первоначально распространения, на боковую поверхность и на горизонтальную плоскость для двух скачкового распространения, то есть с учетом однократного отражения от поверхности земли. Численное моделирование показало, что каустика (огибающая лучей второго скачка) модулируется волнообразным возмущением. При этом возникают точки заострения (точки возврата траекторий), а также петли траекторий и сложные фокусировки. Эти эффекты, характерные для обыкновенной волны, еще более заметны для необыкновенной волны. Показано, что можно подобрать такую частоту, для которой сдвиг лучей, прошедших ионосферную плазму, в боковую плоскость будет любой наперед заданной величиной, однако в силу расходимости лучевого потока интенсивность поля, образованного такими лучами будет крайне мала. Следует отметить характерный "выброс вправо", описывающий выход из основной плоскости распространения лучей в результате взаимодействия со слоем Е. Проанализировано образование каустик пространственно временных геометрооптических лучей и допле-ровское изменение частоты, вызванное нестационарностью ионосферной плазмы.
Для цитирования:
Крюковский А.С., Лукин Д.С., Растягаев Д.В., Скворцова Ю.И. Численное моделирование распространения пространственно-временных частотно-модулированных радиоволн в анизотропной среде // T-Comm: Телекоммуникации и транспорт. - 2015. - Том 9. - №9. - С. 40-47.
For citation:
Kryukovsky A.S., Lukin D.S., Rastyagaev D.V., Skvortsova Yu. I. The numerical modeling of propagation of space-time frequency-modulated radiowaves in ionospheric to medium. T-Comm. 2015. Vol 9. No.9, рр. 40-47. (in Russian).
I. Введение
В работе выполнено исследование распространения частотно-модулированных (ЧМ) сигналов в магнитоактив-ной ионосферной плазме. Активное применение радиоволн декаметрового диапазона для обеспечения дальней радиосвязи, радиолокации, радионавигации, загоризонт-ного радиозондирования, а также изучения структуры верхней атмосферы Земли - ионосферы определяет актуальность исследования. Несмотря на значительный объем работ в данном направлении (см., например, [1-3]), распространение ЧМ сигналов в анизотропной ионосферной плазме изучено ещё недостаточно. Целью работы является исследование распространения обыкновенных и необыкновенных радиоволн в условиях как стационарных, так и нестационарных двухслоевых моделей ионосферного слоя,
В основе работы лежат лучевые методы, являющиеся традиционными при описании процессов распространения коротких электромагнитных волн в различных средах. При моделировании лучевых траекторий в анизотропной, неоднородной, нестационарной среде бихарактеристическая система уравнений с гамильтонианом: 2
Г =к\ + к\ + к* -^e(r,Íc,&) с~
имеет вид [4]:
dk _ ЭГ\ г,
с/т дк с!т дг ilr ёсо c/r 8t где F ~(x,y,z) ~ координаты точки наблюдения, к = (kJC,ky,kI) ~ волновой вектор, о> - круговая частота
излучения, í - групповое время, т - параметр вдоль лучевой траектории, s{r,k,m) ~ эффективная диэлектрическая проницаемость среды распространения, с = 2,997925 I08 м/с - скорость света.
Система (1-2) (см. [5]) обобщает бихарактеристиче-скую систему, использованную нами в [6] для лучевого описания распространения стационарного сигнала в ионосфере. Эффективная диэлектрическая проницаемость среды для неоднородной анизотропной ионосферы имеет вид (см., например, [4, 7]):
*±=i----, , ■ (3)
2(1 - v) — u sin2 a ± sin4 a + 4г/(1 - v)~ cos1 a
В формуле (3) знак «+» соответствует обыкновенной волне, а знак «-» - необыкновенной волне. Параметры и и V обозначают отношения квадратов гирочастоты и плазменной частоты к квадрату рабочей частоты:
(4)
Щ^К+Щк+НМ .
<5)
(I)
(2)
СО,
со I т . с~ со \ со ¡ т.в> 1-Ю
со,,
4 ne'N
В приведённых ниже расчетах амплитуда магнитного поля предполагается постоянной, а ориентация напряжённости магнитного поля относительно локальной системы координат задается двумя постоянными углами г и ц
Нцх = Ип cos/cosp. Ип - H^Qosysmcpi
H0:=H0smr- (6)
Предположим, что начальный волновой вектор к(0) параметрически зависит от углов выхода луча:
СО I— „ , (О
£Д0) = —Js^ cosacos/;» к (0) = —Je^ sin £ cost] <
с с
Aj(0) = — ^e^smij' (J)
а источник излучения точечный, находится вне магнитоак-тивной плазмы и расположен в начале координат (х=у=2=0):
?|,=0 = (8)
Величина £0 в выражениях (7) это значение эффективной диэлектрической проницаемости среды в источнике.
Предположим, что сигнал линейно частотно-модулированный (ЛЧМ), то есть [5]:
й>(0) = й>0(1+ /(0) = г0 - (8)
2. Модель плазменного ионосферного слоя
Предположим, что распространение электромагнитной волны изначально осуществляется в плоскости (х, г). Пусть модель ионосферной плазмы двухслойная, а профиль электронной концентрации имеет вид [6-8]:
ехр(- в)
N(r)= ^gjexp-
1-0—
cosx
Г г Vil
+ /?ехр -
. ' J
6=
z™i/2
(9)
где
где е = 4,8029 10 СГСЭ - заряд электрона, те = 9,108
10~2" г - масса электрона, N - величина электронной концентрации в фиксированной точке пространства, а Н, -величина напряженности магнитного поля Земли, Величина а это угол между волновым вектором и напряженностью магнитного поля Земли яо = (#0,,Я0 ,Я01)- Выражение для квадрата косинуса угла а имеет вид:
8 = 1-Гхсо^(2ф-У 0/Ля), (Ю)
г0] - высота максимума слоя £2, гт] - условная полутолщина слоя р - безразмерный коэффициент, характеризующий степень ионизации нижнего слоя по отношению к основному, г02 - высота максимума нижнего слоя Е, 2т2 ~ условная полутолщина нижнего слоя Е, величина это электронная концентрация в максимуме основного слоя Р2.
Функция g описывает влияние горизонтального градиента с глубиной у{ и волнообразного возмущения в плазме. В формуле (10) Х1 - длина гравитационной волны, а V - скорость.
При вычислении применялись следующие значения параметров: Ы0 = 2 I О6 см"5, Н0 = 0,36 Э, гт1 = 140 км,
г01 = 300 км, 2„7 = 40 км, Г(Н = ' Ю км, С, -0, у =-72°,
ф =90°, V = 230 м/с, X = 50 км, ¡3 = 0,1 I, % = 0.
43
2. Иванов В. Л., Иванов Д В., Лащевский А Р., Рябова М. И, Исследование дисперсионных искажений широкополосных элементов непрерывного ЛЧМ-сигнала при изменении их длительности сверх критической II Вестник Поволжского государственного технологического университета. Серия: Радиотехнические и инфокоммуникационные системы, 20!4. №1(20). С. 43-53,
3. Зохорое В.Е., Котоеа Д.С. Моделирование дисперсионных искажений ЛЧМ-импульсов радиоволн в ионосфере Н Труды XXIV Всероссийской научной конференции «Распространение радиоволн», (29 июня-5 июля 2014; Иркутск). Иркутск: ИСЗФ СО РАН. 2014, Т.4. С, 97-100.
4. Лукин Д.С., Патин Е.А. Численный канонический метод в задачах дифракции и распространения электромагнитных волн в неоднородных средах. М.: МФТИ, 1982. 159 с,
5. Крюковский A.C., Растягаев ДВ,, Скворцов о Ю.И. Исследование распространения частотно-модулированных пространственно-временных сигналов в неоднородной анизотропной ионосфере Н Вестник Российского нового университета. Серия «Управление, вычислительная техника и информатика». М.: Рос-НОУ, 2013. Выпуск 4. - С. 47-52.
6. Крюковский A.C., Лукин Д.С., Кирьянова К.С. Метод расширенной бихарактеристической системы при моделировании распространения радиоволн в ионосферной плазме. II Радиотехника и электроника. М.: Наука. 2012. Т.57. №9. С. 1028-1034,
7. Крюковский A.C., Лукин Д.С., Растягаев ДВ. Исследование особенностей распространения коротких радиоволн в неоднородной анизотропной ионосфере II Электромагнитные волны и электронные системы. 2009. Т.14. Ыэ 8, С. 17-26.
8. Крюковский A.C., Лукин Д.С, Растягаев ДВ. Исследование влияния локальных неоднородностей ионосферной плазмы на распространение коротких радиоволн // Вестник Российского нового университета. Серия «Управление, вычислительная техника и информатика». М.; РосНОУ, 2010. Выпуск 3. С. 17-25.
9. Крюковский A.C., Растягаев ДВ., Скворцова Ю.И. Распространение частотно-модулированных пространственно-временных радиоволн в анизотропной ионосфере И Труды XXIV Всероссийской научной конференции «Распространение радиоволн», (29 июня - 5 июля 2014; Иркутск). Иркутск: ИСЗФ СО РАН, 2014. Т.4. С. 126-129.
10. Крюковский A.C., Скворцова Ю.И. Применение теории катастроф для описания пространственно-временной структуры частотно-модулированного сигнала е плазме II Электромагнитные волны и электронные системы. 20i3. Т.18. №8. С.18-23.
11. Крюковский АС., Лукин Д.С., Палкин Е.А Применение теории краевых катастроф для построения равномерных асимптотик быстроосцилпирующих интегралов II Дифракция и распространение волн, Междувед. Сборник. МФТИ. И., 1985. С.4-21.
12. Крюковский A.C. Равномерная асимптотическая теория краевых и угловых волновых катастроф. Монографий. М.: РосНОУ. 2013.- 368 с.
13. Крюковский АС., Скворцова Ю.И. О распространении частотно-модулированного излучения в ионосферной плазме II Труды XV международной научной конференции «Цивилизация знаний: российские реалии», Москва, 25-26 апреля 2014 г. / М.: РосНОУ, 2014. Часть 2.
14. И патов ЕЁ., Крюковский A.C., Лукин Д. С., Палкин Е.А., Растягаев ДВ. Методы моделирования распространения электромагнитных волн в ионосфере с учетом распределений электронной концентрации и магнитного поля Земли II Радиотехника и электроника, 2014, Т. 59. №12. С. I 180-1187,
15. Крюковский A.C. Необходимые и достаточные условия образования основных волновых катастроф с корангом, равным двум II Распространение и дифракция электромагнитных волн, Междуведсб. М.: МФТИ. 1993. С. 4-19.
16. Крюковский A.C., Кирьянова К.С. Динамическое моделирование распространения радиоволн в окрестности экваториальной аномалии на основе метода бихарактеристик И Электромагнитные волны и электронные системы. 201 I. Т. 16. №8. С. 21-25.
17. Лукин Д.С., Крюковский A.C., Черняк Я.М. Анализ влияния моделей магнитного поля при численном моделировании распространения коротких волн в ионосфере Земли II T-Comm: Телекоммуникации и транспорт. 2014. Т.8. №12. С.55-58.
18. Крюковский A.C., Лукин Д.С., Растягаев Д Б. Моделирование лучевой и каустической структуры электромагнитных полей по данным радиотомографии ионосферы в окрестности экваториальной аномалии. II Электромагнитные волны и электронные системы. 2010. Т.! 5. №8. С. 5-1 I.
19. Кирьянова К.С., Крюковский A.C. Особенности лучевого распространения радиоволн в окрестности пространственно-модулированного ионосферного слоя II Труды XIII международной научной конференции «Цивилизация знаний: проблемы и перспективы социальных коммуникаций», Москва, 20-21 апреля 2012 г. - Часть 2. - М.: РосНОУ, 2012. - С. 30-34.
THE NUMERICAL MODELING OF PROPAGATION OF SPACE-TIME FREQUENCY-MODULATED
RADIOWAVES IN IONOSPHERIC TO MEDIUM
Andrey Sergeevich Kryukovsky, Moscow NOU VPO Russian new University, doctor of physico-mathematical Sciences, Professor,
Dean of the faculty IICT, [email protected] Dmitry Sergeevich Lukin, Moscow NOU VPO Russian new University Professor of Tiiis, D. F.-M. N. [email protected] Dmitry Vladimirovich. Rastegaev, Moscow NOU VPO Russian new University, candidate of physico-mathematical Sciences,
associate Professor, Vice-rector for it [email protected] Yulia Igorevna Skvortsova, Moscow NOU VPO Russian new University, Deputy Dean of the faculty IICT [email protected]
Abstract. In the case of the ordinary and extraordinary waves the numerical modeling two hop radio propagation in anisotropic ionospheric plasma frequency modulation. Active use of decameter radio waves to provide long-distance radio communications, radar, navigation, horizon radiosonde and study the structure of Earth's upper atmosphere - the ionosphere determines the relevance of the study. Despite the considerable amount of work aimed at studying the spread frequency modulated radio waves, the spread of such signals in an anisotropic ionospheric plasma is still insufficiently studied. The aim of the work is to fill this gap. The basis of the study is a numerical solution of the space-time bicharacteristic Hamiltonian system, in which the unknowns are the components of the wave vector, the coordinates of the beam, as well as the frequency and time. It is believed that the source of the radiation spot is outside the magnetized plasma, linear frequency modulation and amplitude of the magnetic field is constant. A model of a two-layer of the ionosphere both in the presence and in the absence of non-stationary wave-like disturbances modeling gravitational waves in the plasma. It is known that the external magnetic field leads to the exit path of the beam propagation plane. The paper compares the projection beam structures ordinary and extraordinary waves on the vertical plane of the original distribution, on the side surface and the horizontal plane for two intermittent distribution, i.e. taking into account the single reflection from the surface. Numerical simulations show that the caustic (envelope of the second rays jump) modulated wave-like disturbance. Herewith the cusps (the cusps of trajectories), and a loop trajectory and complex focusing. These effects are specific to the ordinary wave even more noticeable for the extraordinary wave. It has been shown that it is possible to select a frequency for which the shift rays passing ionospheric plasma in the side plane is any predetermined value, however, due to the divergence of radiation flux intensity of the field formed by these rays is very small. Note the characteristic "emission right", describing a way out of the main plane of propagation of rays resulting from the interaction with the layer E. This work analyzes the formation of the geometrical space-time caustics rays and Doppler frequency change caused by the nonstationary ionospheric plasma.
Keywords: an ordinary and unusual wave; numerical modeling; the propagation of frequency-modulated signals; anisotropic ionosphere; bi-charac-teristic system; wavy perturbations; space- time caustics; Doppler change of frequency.
References
1. Ivanov D.V. Methods and mathematical model study of propagation in the ionosphere complex decameter signals and correction of dispersion distortions. Monograph. Yoshkar-Ola: Mari state technical University. 266 p. (2006). (in Russian)
2. Ivanov V.A., Ivanov D.V., Ladewski A.R., Ryabov M.I. Study of the dispersion distortion broadband continuous LFM signal when changing their duration beyond critical / Bulletin of the Volga region state technological University. Series: Electronic and information and communication systems. No.1(20). Pp. 43-53 (2014).
3. Zakharov V.E., Kotov, D.S. Modeling of dispersion distortions LFM pulses of radio waves in the ionosphere. / Proceedings of the XXIV all-Russian scientific conference "radio wave Propagation", (29 June-5 July 2014; Irkutsk / Irkutsk: ISSF SB RAS. Vol. 4. Pp. 97-100. (2014). (in Russian)
4. Lukin, D.S., Palkin E.A. Canonical numerical method in problems of diffraction and propagation of electromagnetic waves in inhomogeneous media. Moscow, MIPT, 159 p. (2014). (in Russian)
5. Kryukovsky A.S., Rastyagaev D.V., Skvortsova Y.I. The study of the propagation of frequency modulated space-time signals in an inhomogeneous anisotropic ionosphere / Journal of the Russian new university. Series "Management, computing and Informatics), Moscow, RosNOU, Issue 4. Pp. 47-52 (2013).
6. Kryukovsky A.S., Lukin D.S., Kiryanova K.S. Extended bicharacteristics system for radio wave propagation modeling in the ionospheric plasma. // Radio engineering and electronics, Moscow: Nauka. Vol. 57. No. 9. Pp. 1028-1034 (2012). (in Russian)
7. Kryukovsky A.S., Lukin D.S., Rastyagaev D.V. The study of propagation of short radio waves in an inhomogeneous anisotropic ionosphere. / Electromagnetic waves and electronic systems. Vol. 14. No. 8. Pp. 17-26 (2009). (in Russian)
8. Kryukovsky A.S., Lukin D.S., Rastyagaev D.V. Research of influence of local inhomogeneities of the ionospheric plasma on the propagation of short radio waves. //Bulletin of the Russian new University. Series "Management, computing and Informatics), Moscow, RosNOU, Issue 3. Pр. 17-25 (2010). (in Russian)
9. Kryukovsky A.S., Rastyagaev D.V., Skvortsova Yu.I. The propagation of frequency modulated space-time of radio waves in an anisotropic ionosphere / Proceedings of the XXIV all-Russian scientific conference "radio wave Propagation", (29 June-5 July 2014; Irkutsk / Irkutsk: ISSF SB RAS, Vol. 4. Pр. 126-129 (2014).
10. Kryukovsky A.S., Skvortsova Y.I. Application of catastrophe theory to description of the space-time structure of the frequency-modulated signal in plasma / Electromagnetic waves and electronic systems. 2013. Vol. 18. No. 8. Рр. 18-23. (in Russian)
11. Kryukovsky A.S., Lukin D.S., Palkin E.A. The application of the theory of boundary value catastrophes for the construction of uniform asymptotics of fast oscillating integrals / Diffraction and wave propagation. Midwived. collection / MIPT. M., 1985. Pp. 4-21. (in Russian)
12. Kryukovsky A.S. Uniform asymptotic theory of edge and corner wave catastrophes. Moscow. RosNOU, 2013. 368 p. (in Russian)
13. Kryukovsky A.S., Skvortsova, Y.I. About the spread frequency modulated radiation in the ionospheric plasma / Proceedings of XV international scientific conference "Civilization of knowledge: the realities of Russia", Moscow, 25-26 April 2014. Moscow, RosNOU, 2014. Part 2. (in Russian)
14. Ipatov E.B., Kryukovsky A.S., Lukin D.S., Palkin EA., Rastyagaev D.V. modeling and simulation of electromagnetic wave propagation in the ionosphere taking into account distributions of the electron concentration and the magnetic field of the Earth / radio engineering and electronics. 2014. Vol. 59. No. 12. Pp. 1180-1187.
15. Kryukovsky A.S. Necessary and sufficient conditions for the formation of the main wave of disasters by koranga equal to two / Propagation and diffraction of electromagnetic waves. Midwived.sat./Moscow, MIPT. 1993. Pp. 4-19. (in Russian)
16. Kryukovsky A.S., Kiryanova K.S. Dynamic modeling of radio wave propagation in the vicinity of the Equatorial anomaly on the basis of the method of bicharacteristic / Electromagnetic waves and electronic systems. 2011. Vol. 16. No. 8. Pp. 21-25. (in Russian)
17. Lukin, D.S., Kryukovsky A.S., Chernyak, Ya.M. analysis of the influence of magnetic field models for numerical simulation of the propagation of short waves in the Earth's ionosphere / T-Comm. 2014. Vol. 8. No. 12. Pp. 55-58. (in Russian)
18. Kryukovsky A.S., Lukin D.S., Rastyagaev D.V. The modeling of radiation and caustic structure of electromagnetic fields according to radiotelegraphy of the ionosphere in the vicinity of the Equatorial anomaly / Electromagnetic waves and electronic systems. 2010. Vol. 15. No. 8. Pp. 5-11. (in Russian)
19. Kiryanova K.S., Kryukovsky A.S. Singularities of the radiation propagation in the vicinity of spatially modulated ionospheric layer / Proceedings of XIII international scientific conference "Civilization of knowledge: problems and prospects of social communication", Moscow, 20-21 April 2012- Part 2. Moscow. RosNOU, 2012. Pp. 30-34. (in Russian)