Использование интегрального и дифференциального методов теории дифракции для прогноза напряженности поля над земной поверхностью Текст научной статьи по специальности «Физика»

Электронное научно-техническое издание

НАУКА и ОБРАЗОВАНИЕ

Эя №ФС 77 - 30569. Государственная регистрация №0421100025.155Н 1994-0408_

Использование интегрального и дифференциального методов теории дифракции для прогноза напряженности поля над земной поверхностью

77-30569/247813

# 11, ноябрь 2011 Ахияров В. В.

УДК. 621.371

НИИРЭТ МГТУ имени Н.Э. Баумана [email protected]

Задачи, связанные с определением характеристик электромагнитного поля, рассеянного объектами, называются дифракционными. В самой общей постановке дифракционная задача состоит в определении комплексных векторов напряженности электрического и магнитного поля, удовлетворяющих системе уравнений Максвелла, граничным условиям на поверхности объекта и условиям излучения на бесконечности.

Наиболее общим подходом при решении дифракционных задач является метод, основанный на численном решении интегральных уравнений (ИУ) Стреттона-Чу относительно неизвестной функции плотности поверхностного тока [1]. Для численного решения интегральных уравнений используется метод разложения неизвестной функции по N базисным функциям (число N равно числу элементарных площадок, на которые разбивается поверхность объекта). Далее требуется использовать N весовых функций, что приводит к системе линейных алгебраических уравнений (СЛАУ) с матрицей размерностью N х N .

Считается, что для получения приемлемой точности линейный размер площадки не должен превышать //(5...7), где X - длина волны в свободном пространстве, поэтому при

увеличении размеров анализируемого объекта или при увеличении частоты размерность матрицы СЛАУ стремительно увеличивается. В качестве примера в таблице 1 приведены требования к объему оперативной памяти, необходимой для хранения элементов матрицы СЛАУ при расчете распределения токов на поверхности параболической антенны. Видно, что при прямом численном решении интегрального уравнения требования к оперативной памяти существенно превышают возможности современного персонального компьютера.

Таблица 1 - Оценка требуемой оперативной памяти [2]

Диаметр антенны N Память

19X 100 000 150 Гб

38Х 400 000 2 400 Гб

72Х 1 500 000 33 500 Гб

Однако, существует большое число скалярных дифракционных задач, решение которых можно получить за приемлемое время без использования распределенных вычислений. Прежде всего, это относится к численным методам расчета напряженности поля над земной поверхностью. Данные методы пришли на смену известным аналитическим формулам, которые основаны на решении модельных задач дифракции при использовании допущений о том, что земная поверхность является гладкой и сферической, вертикальная стратификация атмосферы - линейная и т.д.

Численные методы свободны от подобных ограничений и основной вопрос, возникающий при их использовании, заключается в реализации эффективных вычислительных алгоритмов. В данной работе для прогноза напряженности поля с учетом рельефа земной поверхности рассматриваются методы, основанные на численном решении интегрального и параболического уравнений.

Поскольку в УКВ диапазоне практически все типы земной поверхности являются диэлектриками, можно считать, что при скользящем падении радиоволн коэффициент отражения при вертикальной и горизонтальной поляризации Яв г ~ — 1, что соответствует

краевому условию Дирихле. Для вертикальной поляризации поля источника это соответствует отражению от поверхности идеального магнитного проводника, и в этом случае имеют место граничные условия:

п х Н = 0, = 0, (1)

1

где Е и Н - напряженности электрического и магнитного поля на рассеивающей поверхности, п - вектор нормали, направленный в область вычисления рассеянного поля (рис. 1).

Далее, для получения приемлемых по точности результатов при разумных вычислительных затратах, полагаем, что поверхность является цилиндрической с профилем, не зависящим от координаты у. При вертикальной поляризации излучения (Н0 1 Х02) в пределах первой зоны Френеля это эквивалентно условию:

п ■ Н « 0

(2)

п

п

К,

X

Рис. 1. Определение поверхностной плотности тока

Для численного решения рассматриваемой задачи интегральные формулы Стреттона-Чу, с учетом (1) и (2) можно представить в виде скалярных ИУ Фредгольма 1-го и 2-го рода. Как показали вычислительные эксперименты [3], предпочтительным является ИУ Фредгольма 1-го рода:

Н0 (р )=^ | М (р')Н 02)(№2 К

4 7

(3)

где Н0 (р) - поле источника, М (р') - искомое распределение плотности поверхностного магнитного тока, Я2 - расстояние между точками интегрирования р' и наблюдения р

(рис.1), Н02)(х) - функция Ханкеля второго рода нулевого порядка, Ь - профиль трассы, £"о - диэлектрическая проницаемость вакуума, а - круговая частота, к = 2п// - волновое число (X - длина электромагнитной волны).

Для численного решения ИУ (3) следует представить в виде системы линейных алгебраических уравнений (СЛАУ):

& ]=[гпт ][мт ]. (4)

В (4) использованы следующие обозначения: ] - дискретные значения источника поля, в качестве которого рассматривается нить магнитного тока единичной амплитуды, \Мт ] - искомая плотность поверхностного тока, элементы матрицы т ] вычисляются следующим образом:

7 =

п,т

(ОБ0

Н02\кЯ2 )А/, п ф т

(0 Н(2/1 А/Л

4 Н 0

к 2

А/, п = т

где А/ - шаг по профилю трассы.

Для получения приемлемого по точности результата требуется выполнение условия А/ < X/3, что приводит к необходимости решения СЛАУ огромного размера. Однако, если не учитывать обратное рассеяние, СЛАУ (4) решать не требуется, поскольку в этом случае матрица [7п т ] является нижней треугольной. Поэтому для вычисления плотности

поверхностного тока используется простой вычислительный алгоритм [3]:

1

Мп * — 7пп

п-1

Qn - Е 7птМт

V т=0 У

Поле в точке наблюдения ч над рассеивающей поверхностью определяется выражением [3]:

Н (ч | 02 \кЭ )с1/,

4 ь

где О = |q — q'|, q и с/ - точки наблюдения и интегрирования, q' е Ь .

Полное поле в точке q определяется не только рассеянным, но и первичным полем, поэтому для получения окончательного результата следует вычислить сумму Н (ч) и

Н0 (Ч).

При решении задачи дифракции над земной поверхностью важным является вопрос о соответствии приближения идеального магнитного проводника реальным электрическим параметрам трассы распространения. Сравнение численного расчета для задачи дифракции с краевым условием Дирихле и аналитического решения для поверхности с электрическими параметрами, соответствующими сухой почве, показало, что в УКВ диапазоне приближение идеального магнитного проводника вполне приемлемо для вычисления напряженности поля над земной поверхностью [3].

Результаты расчетов множителя ослабления для реального профиля трассы, представленного на рис. 2, показаны на рис. 3. Решение было получено при следующих исходных данных: X = 3 м, высоты источника и приемника = ^ = 10 м, А/ = //3. Невязка краевого условия 8 < —40 дБ, что свидетельствует о высокой точности вычислений.

Рис. 2. Профиль трассы

Рис. 3. Дистанционная зависимость множителя ослабления

Идея метода параболического уравнения принадлежит М.А. Леонтовичу, который преобразовал исходную краевую задачу для уравнения Гельмгольца (эллиптического типа) к уравнению параболического типа. Далее М.А. Леонтович и В.А. Фок получили аналитическое решение параболического уравнения (ПУ) для частных случаев распространения радиоволн в линейной и слоисто-однородной атмосфере [4]. В настоящее время существуют эффективные алгоритмы численного решения ПУ для произвольного профиля рельефа и показателя преломления п(х, г) [5]. Поскольку рассматриваемая задача считается плоской, при горизонтальной поляризации излучения отличной от нуля является поперечная компонента электрического поля, а при вертикальной поляризации решение ищется для поперечной компоненты магнитного поля.

Параболическое уравнение для амплитуды поля и (х, г) над сферической моделью

Земли, справедливое в малоугловом (параксиальном) приближении при |тах| < 150 (см.

рис. 4) имеет вид:

— =1 + к2 (т(х, г)2 - 1)1и(х, г) = 0, (5)

дх 2 I дг2 ^ 7 Г ;

где т(х, г) = п(х, г)+— - модифицированный показатель преломления (а - радиус Зем-

г а

ли).

Рис.4. Геометрия задачи

Для численной реализации уравнение (5) необходимо дополнить граничным условием (как и при решении ИУ, хорошим приближением будет краевое условие Дирихле).

Алгоритм решения уравнения (5) заключается в следующем: на дальности х поле и(х, г) разлагается в угловой спектр плоских волн, который умножается на передаточную

функцию слоя пространства G( p )= e

jk

2

!Л -1

k2

Л

Ах

. Далее вычисляется обратное преобра-

зование Фурье, соответствующее распределению поля по высоте на дальности х + Ах, и результат умножается на фазовый множитель, учитывающий рефракцию радиоволн:

и(х + Ах, z )= e

jkí«(x,z)2 -1+2Z 1ах 2 ^ a) р -\

F [u(x, z )]e

jk

i-4 -i k2

Л

Ах

(6)

где p = ksin O , F[...] и F 1 [...] - прямое и обратное преобразование Фурье.

Отметим, что условие Дирихле соответствует синус-преобразованию Фурье, однако этот же результат можно получить с помощью преобразования Фурье, если дополнить область расчета ее зеркальным изображением. Последний вариант является более предпочтительным, поскольку быстрое преобразование Фурье реализовано практически во всех программах компьютерной математики.

Поле источника излучения моделируется линейной апертурой с известным распределением поля. Для граничных условий Дирихле преобразование Фурье начального распределение поля u(0, z) имеет вид:

U(0, p) = F(p - pA )e-pZA - F(pa - p j ,

где Za - высота подъема антенны, pa = ksinOa (Oa - угол наклона диаграммы направленности к горизонту).

Для того, чтобы учесть влияние рельефа на напряженность поля, используется метод ступенчатого моделирования, при этом реальная трасса распространения заменяется вертикальными ступеньками, как это показано на рис. 5. При распространении над горизонтальными участками поле вычисляется по алгоритму (6), а при встрече со ступеньками $2 - приравнивается к нулю (см. рис. 6).

При численной реализации (6) возникает очевидная трудность, связанная с отражением волн от искусственной границы, расположенной на высоте ± Zтах. Для ограничения расчетной области сверху и ее зеркального изображения снизу необходимо дополнительное условие, гарантирующее отсутствие источников поля на бесконечности. Данное

условие задается на дополнительной границе < |^тах|, при этом решение ищется внутри слоя — ^ 2 < [6].

Рис. 5. Моделирование геометрии рельефа Рис. 6. Алгоритм вычисления амплитуды [5] поля [5]

Результаты расчетов множителя ослабления для реального профиля рельефа протяженностью 40 км и линейной модели атмосферы при высоте антенны к = 30 м и

X = 1 м представлены на рис. 10.

Рис. 7. Дифракция на реальном профиле рельефа

Сравнивая рассмотренные методы решения скалярных задач дифракции для прогноза напряженности поля над земной поверхностью, необходимо отметить, что метод ИУ является строгим электродинамическим подходом, который основан на использовании интегральных формул Стреттона-Чу. Основной вопрос, связанный с использованием метода ИУ, заключается в снижении времени вычислений. Метод ПУ является малоугловым приближением уравнения Гельмгольца. Однако, при распространении радиоволн над земной поверхностью малоугловое приближение выполняется автоматически: максимальный поперечный размер области расчетов, как правило, не превышает нескольких километров, в то время как продольное расстояние оказывается существенно большим. По сравнению с решением интегрального уравнения, метод ПУ требует существенно меньших вычислительных затрат и позволяет получить решение задачи дифракции радиоволн над земной поверхностью с учетом диаграммы направленности источника излучения и вертикальной стратификации атмосферы.

Литература

1. Неганов В.А., Осипов О.В., Раевский С.Б., Яровой Г.П. Электродинамика и распространение радиоволн. / Учеб. Пособие для вузов. Под ред. Неганова В.А. и Раевского С Б. - М.: Радио и связь, 2005.

2. Ашихмин А.В., Пастернак Ю.Г., Попов И.В., Рембовский Ю.А. Обзор принципов построения, возможностей и эффективности программных средств численного электродинамического моделирования // Антенны. — 2007. — № 3. — С.64-80.

3. Ахияров В.В. Методы численного решения задачи дифракции радиоволн над земной поверхностью // Электромагнитные волны и электронные системы. — 2010. — № 3. - С.39-46.

4. Фок В.А. Проблемы дифракции и распространения электромагнитных волн. — М.: Сов. Радио, 1970.

5. Levy М.Б. Parabolic equation method for electromagnetic wave propagation. London. IEE. 2000. 336 p.

6. Ахияров В.В. Метод параболического уравнения в теории дифракции. - Успехи современной радиоэлектроники. 2010. №9. с.72-80.

electronic scientific and technical periodical

SCIENCE and EDUCATION

_EL № KS 77 - 3Ü56'». .V;II421100025, ISSN 1994-jMOg_

Use of integral and differential methods of the diffraction theory for predicting field strength above the Earth surface

77-30569/247813

# 11, November 2011 Ahiyarov V.V.

Bauman Moscow State Technical University [email protected]

In this article the author considers numerical methods for field prediction above the Earth surface, based on the solution of the integral and parabolic equations. The author presents diffraction problem solutions for real profiles. It is shown that the parabolic equation method has high computational efficiency and ensures accuracy of the solution which is comparable with the rigorous method of integral equations.

Publications with keywords: integral equation, parabolic equation, radiowave propagation Publications with words: integral equation, parabolic equation, radiowave propagation

Reference

1. Neganov V.A., Osipov O.V., Raevskii S.B., Iarovoi G.P., Electrodynamics and radio waves propagation, Moscow, Radio i sviaz', 2005.

2. Ashikhmin A.V., Pasternak Iu.G., Popov I.V., Rembovskii Iu.A., Antenny 3 (2007) 64-80.

3. Akhiiarov V.V., Elektromagnitnye volny i elektronnye sistemy 3 (2010) 39-46.

4. Fok V.A., The problem of diffraction and propagation of electromagnetic waves, Moscow, Sov. Radio, 1970.

5. Levy M.F., Parabolic equation method for electromagnetic wave propagation, London, IEE, 2000, 336 p.

6. Akhiiarov V.V., Uspekhi sovremennoi radioelektroniki 9 (2010) 72-80.