Всероссийская открытая научная конференция «Современные проблемы дистанционного зондирования, радиолокации, распространения и дифракции волн» - Муром 2022
УДК 537.87 DOI: 10.24412/2304-0297-2022-1-36-48
Восстановление эффективной частоты соударений электронов в магнитоактивной ионосферной плазме на основе амплитудных данных
А.С. Крюковский1, Д.С. Лукин1, Е.В. Михалёва1, Д.В. Растягаев1
1 Автономная некоммерческая организация высшего образования «Российский Новый Университет»
105005, г. Москва, ул. Радио, 22.
E-mail: kryukovsky56@yandex. ru
Разработан метод реконструкции эффективной частоты соударений электронов в ионосфере Земли на основе амплитудных данных об ослаблении частотно-модулированного радиосигнала при вертикальном зондировании ионосферы. Для моделирования применена расширенная бихарактеристическая система Лукина, а также рекуррентная система решения интегрального уравнения. Сопоставлены результаты моделирования лучевых траекторий и амплитудных характеристик для волн о- и х-поляризаций. Проведено восстановление эффективной частоты соударений как по данным обыкновенной, так и необыкновенной волны. Показано, что для обыкновенной волны получено хорошее совпадение исходных значений эффективной частоты соударений и расчётных значений. Ключевые слова: ионосфера, диагностика, эффективная частота соударений, зондирование, бихарактеристики, радиосигнал, каустика, амплитуда, расходимость, обыкновенная и необыкновенная волна.
Reconstruction of the Effective Frequency of Electron Collisions in a Magnetically Active Ionospheric Plasma Based on Amplitude Data
A.S. Kryukovsky1, D.S. Lukin1, E.V. Mikhaleva1, D.V. Rastyagaev1
1 Russian New University
A method has been developed for reconstructing the effective frequency of electron collisions in the Earth's ionosphere based on amplitude data on the attenuation of a frequency-modulated radio signal during vertical sounding of the ionosphere. The extended Lukin's bi-characteristic system, as well as the recurrent system for solving the integral equation, were used for modeling. The results of modeling of ray trajectories and amplitude characteristics for waves of o- and x-polarizations are compared. The effective collision frequency has been reconstructed from both ordinary and extraordinary wave data. It is shown that for an ordinary wave a good agreement is obtained between the initial values of the effective collision frequency and the calculated values.
Keywords: ionosphere, diagnostics, effective collision frequency, sounding, bicharacteristics, radio signal, caustic, amplitude, divergence, ordinary and extraordinary wave.
Введение
Необходимость диагностики и контроля параметров ионосферной плазмы определяется влиянием состояния ионосферы Земли на функционирование систем радиолокации, радиосвязи и радионавигации (позиционирования). Также актуальность работы связана с осуществлением мониторинга экстремальных явлений в атмосфере Земли [1].
В данной работе выполнено моделирование особенностей распространения радиоволн в ионосфере Земли по оперативным данным вертикального зондирования
(ВЗ). Хорошо известно, что использование ионозондов ВЗ это одно из наиболее эффективных и апробированных и средств диагностики ионосферы [3].
Рассмотрен метод восстановления по данным ВЗ эффективной частоты соударений электронов в ионосферной плазме [3]. В отличие от работ [4, 5], в которых рассматривалась изотропная модель ионосферы, в настоящей работе учтено внешнее магнитное поле Земли и исследованы две круговые поляризации: обыкновенная волна (о-волна) и необыкновенная волна (х-волна). Как и в предыдущих работах [4, 5] предполагалось, что из данных ионозонда вертикального зондирования [6-8] известны задержки сигналов (ío) и амплитудные характеристики.
Вычисление лучевых характеристик и расходимости
Известно, что эффективная частота соударений Ve определяет поглощение щ которое в свою очередь определяет амплитуду A зондирующего сигнала, пришедшего в приемник после отражения от ионосферных слоев. Амплитуда связана с поглощением формулой [4, 9]:
A = E0D ехр[-щ] , (1)
в которой D - это расходимость лучевого потока, а Eo - начальная амплитуда поля вне ионосферной плазмы на расстоянии ro от источника.
Предположим, что источник излучения, который для простоты будем считать изотропным, расположен в начале координат. Тогда [4, 5]:
Eo (В/м). (2)
Г0
В (2) W - мощность излучателя. Далее для определенности будем считать, что W= 1 кВт.
Поскольку из эксперимента может быть известна зависимость амплитуды A принятого сигнала от частоты f а также задержка отразившегося сигнала to как функция частоты, то можно найти поглощение щпо формуле
щ = - In— (3)
Щ DEo ()
Поскольку A и Eo известны, для того чтобы найти поглощение щ необходимо вычислить расходимость D. Для этого необходимо решить систему бихарактеристических уравнений Гамильтона-Лукина [9-13]:
df 2кс2 - co2ds / дк dk a>2ds / дг ^
dt d{sco2)ld(o dt d(so)2)l дсо
в которой г = (x.y.z) - декартовы координаты луча, к (!) - волновой вектор, с -скорость света, а> = 2лf - круговая рабочая частота, t - групповое время, a sir.к.со) -эффективная диэлектрическая проницаемость среды распространения. В данной работе будем предполагать, что среда плоскослоитая, то есть s(f ,к,со) = s(z,k,c¿>). Систему (4) необходимо дополнить начальными условиями [11,13]:
кх(0) = -&**С , к v (0) = —Js^ cos С, cos r¡ , kz(0) = —Js^ cos £ sint] , г (0) = 0 , (5) с с с
где е0 -значение эффективной диэлектрической проницаемости среды в источнике
излучения [11], а величины £ r¡ - это начальные углы выхода луча. Для
рассматриваемой модели ионосферы е0& 1 .
Для вертикального луча С = 0, r = ^/2. Поэтому формулы (5) принимают вид:
*х(0) = 0, к (0) = 0 , kz(0) = -& , г(0) = 0 , (6)
7 с
что позволяет упростить бихарактеристическую систему (4), которая принимает вид:
2кгс2 -а2де / дкг а2де / дг
&Х 8 (еа2) / да &Х 8 (ею2)/ да
К (Х) = 0
ёх
а2де / дк„
к, (Х) = 0 , = -
д (ею2)/ да а2де / дк
(8)
(9)
д (ею )/да
Таким образом видим, что система (4) разбивается на три подсистемы. Сначала интегрируется система (7) с учетом равенства нулю кх и к . Затем с учетом решения
системы (7) интегрируются выражения (8) и (9). Необходимо отметить, что выражения содержат сингулярности. Результатом решения систем (7-9) являются зависимости г(х)
и к(Х) для каждой выбранной частоты. Однако этого недостаточно для вычисления расходимости. Для определения расходимость И необходимо знать производные г(х) и к(X) по углам ¡^ и //. Эти величины можно найти, решая расширенную бихарактеристическую систему Лукина [12,14]:
8
дС
(
'С
2кс2 - а2де / дк
&Х дц
(
д (£со ) / да 2кс2 -а2де / дк
\
д (еа )/ да
д
а2де / дг
да уд (еа ) / да д ( а2де / дг
&Х дц! д (еа )/ да
(10) (11)
с начальными условиями:
а
к'хс (0) = , к'хц (0) = 0, к'у( (0) = --^8111 Сес8ц, С (0) = 0 , Ц(0) = 0:
к'уц (0) = -—л/е) со$С этц, к'с (0) = С ъ'тц, к;, (0) = соъц .
(12)
Для вертикального луча формулы (12) упрощаются [4, 5, 15]:
Кс (0) = , к'хц (0) = кУс (0) = к';с (0) = кц (0) = 0,
к'уц(0) = -а4Г0 , С(0) = 0 :
Ц(0) = 0.
(13)
Анализируя (10-11) с учетом (13) можно найти, что:
к'сс(Х) = аТё к'с(Х) = 0, кХц(Х) = 0, кУц(Х) = (14)
Таким образом, в системе (10-11) из 12 остаётся 8 уравнений, которые решаются численно с учетом решения системы (7-9).
Теперь можно перейти к определению расходимости, которая вычисляется как корень из отношения якобиана Jo к J по формулам [12,14]:
Б =
J =
Т Т
У
Ус
у;
Т
г'
• = 4
где
т=-
а2де / дк„
д (еа )/ да
а2де / дк
Т =__у_
' у = д(еа2)/да
, Т =
2к с2 - а2де / дк
г_
д (еа2) / да
(15)
(16)
кх =0
0
х
х
г
ку =0
В отличие от [4,5] в настоящей работе эффективная диэлектрическая проницаемость среды распространения задавалась в виде формулы Эпплтона-Хартри-Лассена:
2 X (1 - X)
е= 1 -
2(1 - X) - и бш2 а±^и2 Бт4 а + 4и(1 - V)2 соб2 а
(17)
в которой
X =
г л 2 9 'а„) 4жв2 N (г)
а
V
и =
т а
Vю J
в2 И,2
2 2 2 т. с а
(18)
угол а это угол между волновым вектором к и вектором напряженности магнитного поля Земли Н0 =(Н0х,Н0у,Н0г), - электронная концентрация, сор - плазменная
круговая частота, —д электрона.
-круговая гирочастота, e
заряд электрона, т
масса
Модель ионосферы
В работе рассмотрена модель электронной концентрации N(z), показанная на рис. 1.
Рис. 1. Зависимость электронной концентрации от высоты
На рисунке зависимость электронной концентрации от высоты показана дважды: как в виде жирной черной линии, так и в виде жёлтого фона. Модель содержит ионосферные слои Б.
N(г) » 70254.26ехр
+ 561828.04 ехр
1 (1 - ехр0 + 0) 1
в = ^ (196 - 2), в2 = ^ (263 - 7).
- (1 - ехр 02 +02)
(19)
(20)
Максимум слоя Б1 расположен на высоте ~196 км, а максимум слоя Б2 - на высоте ~263 км.
Зависимость логарифма эффективной частоты соударений от высоты приведена на рис. 2 (белая линия) [4, 5].
Из рисунка видно, что до высоты приблизительно 250 км, эффективная частота соударений спадает примерно по гиперболическому закону (синяя линия):
2
« 0.623965 +
415.116
(21)
и на больших высотах не очень сильно влияет на ослабление сигнала. Высота 2 выбрана
в километрах.
Рис. 2. Зависимость эффективной частоты соударений от высоты
Метод восстановления эффективной частоты соударений
Зная амплитуду А и начальную амплитуду Ео можно, после вычисления расходимости Б, перейти к определению эффективной частоты соударений \\ то есть восстановить кривую, показанную на рис. 2. Считается, что поглощение волны вдоль лучевой траектории пропорционально интегралу от мнимой части показателя преломления [2]. При расчете амплитуды сигнала по формуле (1) нами было использовано выражение для поглощения радиосигнала в виде интеграла вдоль лучевой траектории:
¥ =
10
в котором ^ - время задержки отражённого от ионосферы сигнала,
,, ч о И,
¥(*) =--Хи^ ,
с И
Ио=~
1 д(о!)
с2 да
И = 2,
К + к2у +
Гк --—д! 2 2 с2 дг
Х =
))2 + (1т!,))2 - Яе(!) .
(22)
(23)
(24)
(25)
В формуле (25) х - это мнимая часть показателя преломления, а ег - это выражение для эффективной диэлектрической проницаемости с учетом частоты соударений [2]:
! = 1--Х , (26)
1 -12 —
и бш а
и2 Бт4 а
2(1 - X -12) у 4(1 - X -12)2
- + и соб а
г
0
2
где
у
г = у. (27)
а
Если предположить, что мнимая часть эффективной диэлектрической проницаемости много меньше вещественной и пренебречь влиянием внешнего магнитного поля, то из формул (23-27) можно получить упрощенное выражение:
у/'(Х) « у>'(?) = ~е2, в2 = 1тег. (28)
При расчете амплитуды А использовались формулы (23-27) и, в частности, формулу (26) для вычисления мнимой части эффективной диэлектрической проницаемости. Однако при восстановлении эффективной частоты соударений для уменьшения времени вычислений вместо (3) использовалась формула (28)
ц/ — ——— (29)
в которой [9,16]:
0
ф = - — ^ё2сИ. (30)
2 0
В формуле (30) выражение в2 (см. (26)) приближённо заменено на ё2 [4,5,15]:
ёг=-Х2. (31)
Объединяя выражения (29)и (30) получаем интегральное уравнение Вольтерра второго рода:
V 'г А 4же2
Г уемш = - ь-А., V, = ^ (32)
- •> /) А т
а 0 ПЕ0
Введем обозначение:
1 л
ь =--1п-, . (33)
_ К
Выражение Ь является функцией рабочей частоты / и вычисляется в точке приёма сигнала. В уравнении (32) - это время, за которое сигнал проходит расстояние от источника до точки отражения от ионосферы (¿о=2 >т). Как и Ь величина 1т также является функцией частоты. Электронная концентрация N и эффективна частота соударений Уе являются функциями высоты г, которая в свою очередь зависит от группового времени так как вычисляется вдоль траектории. Поэтому функция г(.) также зависит от частоты.
Для решения уравнения (32) можно применить метод итераций [4, 5, 15]. Введем обозначение О = у Ыс. Так как существует нижняя граница ионосферы г(.то) ниже
которой электронная концентрация равна нулю, то до некоторой частоты /о, для которой время .т равно .то, интеграл в левой части (32) равен нулю:
1т0
Г ОЛ = 0. (34)
о
Рассмотрим частотный интервал (/о, /м) разделим его на п частей (/о, /1, , //, , /п= /м). Здесь /м - частота, выше которой лучи проходят ионосферные слои и не возвращаются на землю. Будем считать функцию О на каждом таком интервале постоянной и равной О/. Определим величины .т/ для каждой/ В результате получим:
9 2 •_]
о2 а2 _ -^ ^ . а ^ С,--К
Г^ _ т 1 г* — Т 2 т1 т0 Г^ Т
^ = ^ -, ^ = --^ -, ^ = ^
1 г - г г - г г - $ 3 3
I___1 I___А I
3 ^ ^ т. г т.г-1
$ . $ . 1 ,—1 $ ■ $ • 1
т. 1 т, 1—1 1 т. з т. з — 1
22
т1 1т0 1т2 1т1 1т2 1т1 1т. 3 1т. 3-
(35)
После вычислений Су эффективная частота соударений находится делением Су на значение электронной концентрации М^) в соответствующей точке.
Результаты численного моделирования
Рассмотрим численную реализацию изложенных выше методов. На рис. 3а (а также на рисунках 4-8 и 10) показана лучевая структура обыкновенной волны, а на рис. 3б -лучевая структура необыкновенной волны.
На рис. 3-6 приведены результаты решения системы бихарактеристических уравнений (6)-(9). Проекции лучевых траекторий в ионосфере на координатную плоскость (х, 2) показаны на рис.3. Видно, что за счет влияния магнитного поля волна при распространении в ионосфере отклоняется от вертикали, причем для обыкновенной и необыкновенной волны это отклонение противоположно [2] и значительно меньше.
а) б)
Рис. 3. Проекции лучевых траекторий на координатную плоскость (х, I): а) -обыкновенная волна, б) - необыкновенная волна
После отражения от ионосферы луч возвращается в источник излучения по той же траектории. Цветом на рисунках отмечена рабочая частота: от фиолетового (минимальная частота) до красного (максимальная частота). Плотность электронной концентрации ионосфера показана фоном. В легенде 1 соответствует 106 см-3.
а) б)
Рис. 4. Лучевые траектории в координатах (£, I): а) - обыкновенная волна, б)
необыкновенная волна
На рис. 4 показаны лучевые траектории в координатах (¿, £) для обыкновенной и необыкновенной волны. При этом частота / меняется в соответствии с цветами спектра от 1.6 МГц (фиолетовая линия) до 6.993 МГц (красная линия) для обыкновенной волны и от 1.6 МГц (фиолетовая линия) до 7.65 МГц (красная линия) для необыкновенной волны. Время распространения сигнала возрастает с ростом частоты.
На рис. 5 показаны зависимости вертикальной компоненты волнового вектора к2 от времени X для волн о- и х-поляризаций.
а)
б)
(а) (б)
Рис. 5. Зависимость вертикальной компоненты волнового вектора от времени: а) -обыкновенная волна, б) - необыкновенная волна
По рисунку можно проследить, как меняется вертикальная компонента волнового вектора вдоль лучевой траектории во времени.
На рис. 6 показаны зависимости вертикальной компоненты волнового вектора к от высоты г.
а) б)
Рис. 6. Зависимость вертикальной компоненты волнового вектора от высоты: а) обыкновенная волна, б) - необыкновенная волна
Видно, что кривые симметричные и основное изменение волнового вектора происходит вблизи точки отражения от ионосферы. Смысл цветовой раскраски рис. 5 и 6 тот же, что и рис. 3 и 4.
Перейдем теперь к амплитудным характеристикам. Сначала оценим расходимость Б (15). Решая расширенную бихарактеристическую систему (10)-(14), найдем ослабление А сигнала, обусловленного расходимостью, по формуле
Ат = -2018
J
(36)
На рис. 7 приведены зависимости ослабления Ат сигнала как функции времени вдоль лучевых траекторий, на рис. 8 зависимости ослабления А от высоты.
а) б)
Рис. 7. Зависимость амплитуды сигнала от времени: а) - обыкновенная волна, б)
необыкновенная волна
а) б)
Рис. 8. Зависимость амплитуды сигнала от высоты: а) - обыкновенная волна, б) -
необыкновенная волна
Как и на предыдущих рисунках, разными цветами отмечены лучи, принадлежащие разным частотам. Зависимости для обыкновенной и необыкновенной волн различаются незначительно.
Из рис. 7 и 8 видно, что при движении вдоль луча к точке отражения от источника ослабление сначала резко увеличивается (более 100 дБ), в окрестности точки отражения (которая соответствует каустике [1 6]) ослабление А уменьшается, в точке отражения А , а далее ослабление опять увеличивается.
Теперь рассмотрим результаты математического моделирования в точке приема отраженного сигнала.
На рис. 9 приведены зависимости высоты отражения сигнала Ъш от частоты, вычисленные для обыкновенной и необыкновенной волны (рис. 9 а) и зависимости от частоты времени 1ш прихода луча в точку отражения (рис. 9 б). Величины Ъш и 1ш растут с ростом частоты/. Как и следовало ожидать, кривые, соответствующие обыкновенной волне, расположены выше кривых, соответствующих необыкновенной волне. При приближении к максимуму ионосферного слоя наблюдается резкий рост кривых.
а)
МЕСС б)
Рис. 9. Зависимость от частоты: а) - высоты отражения, б) - времени прихода луча в точку отражения; синяя линия - обыкновенная волна, красная линия - необыкновенная
волна
Из рис. 9а и 9б легко получаются графики (рис. 10) зависимости высоты отражения сигнала от времени отражения Ш Обе кривые гладкие и взаимно-однозначные, однако они имеют точку пересечения.
Рис. 10. Зависимость высоты отражения от задержки
На рис. 11а показана зависимость поглощения ху от частоты / полученная по формулам (22)-(25), на рис. рис. 11б амплитуды сигналов для обыкновенной (синяя линия) и необыкновенной (красная линия) волны.
а)
МНг б)
Рис. 11. Зависимость от частоты/:(а) -поглощения амплитуды сигнала А
Видно, что поглощение необыкновенной волны особенно велико на более низких частотах. Амплитуда сигнала, отраженного от ионосферы с увеличением частоты, сначала растёт, а потом убывает (рис. 11 б). В правой части рис. 11б видим, что амплитуда необыкновенной волны превышает амплитуду обыкновенной волны. Это объясняется тем, что при одной и той же частоте необыкновенная волна проходит меньшее расстояние, чем обыкновенная волна.
Перейдем теперь к результатам восстановления эффективной частоты соударений.
На рис. 12а и б показаны результаты моделирования зависимость эффективной частоты соударений Уе от высоты г для обыкновенной (а) и необыкновенной волны (б). Белой линией показана заданная (см. рис. 2) зависимость эффективной частоты соударений от высоты, а черной линией - результаты восстановления по предложенному выше алгоритму.
Видно, что, если для обыкновенной волны существует значительная область, где заданные и расчётные значения совпадают с приемлемой точностью, то для необыкновенной волны такая область очень невелика. Здесь следует иметь в виду, что для быстроты вычислений в применяемом алгоритме расчета использовались существенные приближения, которые обсуждались выше.
а) б)
Рис. 12. Сопоставление модельной и расчётной частоты соударений: а) обыкновенная
волна, б) - необыкновенная волна
Тем не менее полученные результаты для обыкновенной волны позволяют методом наименьших квадратов получить формулу для интерполяционной кривой:
450 897
1ёУе « 0.419059 + 450 97 (37)
2
которая хорошо согласуется с формулой (21), построенной по исходным данным.
На рис. 13 показаны три кривые: исходные данные (белая линия), интерполяционная кривая, построенная по исходным данным (21) (красная линия), и результат восстановления по формуле (37) (синяя линия). Видно, что совпадение красной и синей линий удовлетворительное.
501г............................-
2 3 4 5 6 7
lgvt (v„s_l)
Рис. 13. Аппроксимация эффективной частоты соударений
Заключение
Таким образом, разработан метод реконструкции эффективной частоты соударений электронов в ионосфере Земли на основе амплитудных данных об ослаблении частотно-модулированного (ЧМ) радиосигнала при вертикальном зондировании ионосферы. Для численного моделирования была применена расширенная бихарактеристическая система Лукина, а также рекуррентная система решения интегрального уравнения. Сопоставлены результаты моделирования лучевых траекторий и амплитудных характеристик для волн о- и х-поляризаций.
Исследованы параметры радиосигналов в точке приема и проведено восстановление эффективной частоты соударений как по данным обыкновенной, так и необыкновенной волны. Показано, что в рамках тех приближений, которые были использованы в работе для обыкновенной волны получено хорошее совпадение исходных значений эффективной частоты соударений и расчётных значений.
Исследование выполнялось в рамках гранта Российского научного фонда (проект № 20-12-00299).
Литература
1. Ryabova M.I., Ivanov D.V., Chernov A.A., Ovchinnikov V.V. (2020), Experimental Study of Phase-Frequency Responses of HF Radio Channels with the Use of SDR Ionosonde, 2020 Systems of Signal Synchronization, Generating and Processing in Telecommunications (SYNCHROINFO), Svetlogorsk, Russia, P. 1-4, DOI: 10.1109/SYNCHROINFO49631.2020.9166094.
2. Дэвис К. Радиоволны в ионосфере - М.: Мир, 1973. 502 с.
3. Danilkin N.P., Denisenko P.F., Barabashov B.G., Vertogradov G.G. (2005), Electron collision frequency and HF-waves attenuation in the ionosphere, International Journal of Geomagnetism and Аегопоту. V. 5. № 3. P. GI3009. DOI:10.1029/2004GI000081
4. Крюковский А.С., Лукин Д.С., Михалёва Е.В., Растягаев Д.В. Развитие метода восстановления эффективной частоты соударений электронов в ионосферной плазме // Радиотехника и электроника. 2022. Т. 67. № 2. С. 117-129. DOI: 10.31857/S0033849422020085
5. Kryukovsky A.S., Lukin D.S., Mikhaleva E.V., Rastyagaev D.V. (2021), Application of the bi-characteristic method for reconstructing the effective frequency of electron collisions in the ionosphere, Journal of Physics: Conference Series. Сер. "Russian Open Scientific Conference "Modern Problems of Remote Sensing, Radar, Wave Propagation and Diffraction", MPRSRWPD 2021", P. 012001. DOI: 10.1088/1742-6596/1991/1/012001
6. Cedrik M., Podlesnyi A., Kurkin V. (2020), The Different Scale TIDs Diagnostics Based on Chirp-signals Amplitude Measurements Data in Vertical and Near-Vertical Ionosphere Sounding, 2020 7th All-Russian Microwave Conference (RMC), Moscow, Russia, P. 260-262, DOI: 10.1109/RMC50626.2020.9312341.
7. Podlesnyi A., Kurkin V., Cedrik M. (2020), Ionosond-MS ionosonde receiving antenna system, 2020 7th All-Russian Microwave Conference (RMC), Moscow, Russia, P. 263-265, DOI: 10.1109/RMC50626.2020.9312232.
8. Крюковский А.С., Куркин В.И., Ларюнин О.А., Лукин Д.С., Подлесный А.В., Растягаев Д.В., Черняк Я.М. Численное моделирование амплитудных карт для скорректированной модели IRI-2012 с плавными возмущениями ионосферы // Радиотехника и электроника. 2016. Т. 61. № 8. С. 794-799.
9. Бова Ю.И., Крюковский А.С., Лукин Д.С. Распространение частотно-модулированного излучения электромагнитных волн в ионосфере Земли с учетом поглощения и внешнего магнитного поля // Радиотехника и электроника. 2019. Т. 64. № 1. С. 3-14.
10. Казанцев А.Н., Лукин Д.С., Спиридонов Ю.Г. Метод исследования распространения радиоволн в неоднородной магнитоактивной ионосфере. // Космические исследования, 1967. Т. 5. Вып. 4. С. 593-600.
11. Крюковский А.С., Лукин Д.С., Растягаев Д.В., Скворцова Ю. И. Математическое моделирование распространения частотно-модулированных радиоволн в ионосферной плазме // Радиотехника и электроника, 2015, Т. 60, № 10. С. 1001-1009.
12. Крюковский А.С., Лукин Д.С., Кирьянова К.С. Метод расширенной бихарактеристической системы при моделировании распространения радиоволн в ионосферной плазме. // Радиотехника и электроника, М.: Наука. 2012. Т.57. № 9. С. 1028-1034.
13. Крюковский А.С., Растягаев Д.В., Скворцова Ю.И. Исследование распространения частотно-модулированных пространственно-временных сигналов в неоднородной анизотропной ионосфере //Вестник Российского нового университета. Серия «Управление, вычислительная техника и информатика». 2013. Выпуск 4. С. 47-52.
14. Ипатов Е.Б., Крюковский А.С., Лукин Д.С., Палкин Е.А., Растягаев Д.В. Методы моделирования распространения электромагнитных волн в ионосфере с учетом распределений электронной концентрации и магнитного поля Земли // Радиотехника и электроника. 2014. Т. 59. № 12. С. 1180-1187.
15. Kryukovsky A.S., Lukin D.S., Mikhaleva E.V., Popchenko O.V., Rastyagaev D.V. (2020), Diagnostics of the Effective Frequency of Electronic Collisions in the Ionosphere Based on Analysis of the Amplitude Characteristics of Continuous Linear Frequency Modulated Radio Signals, 2020 7th All-Russian Microwave Conference (RMC), IEEE, Moscow, Russia, P. 211-215. DOI: 10.1109/RMC50626.2020.9312352.
16. Крюковский А.С., Лукин Д.С., Бова Ю.И. Моделирование поля в окрестности каустик обыкновенной и необыкновенной волн при ионосферном распространении // Радиотехника и электроника, 2020. Т. 65. № 12. C. 1160-1169. DOI: 10.31857/S0033849420120128