Численное моделирование эффектов дифракции Текст научной статьи по специальности «Математика»

(N

тН

О

(N

УДК 535.42; 517.9 ББК 22.311

А.В. Боровский А.Н. Бородин А.Л. Галкин

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЭФФЕКТОВ ДИФРАКЦИИ

Проведено численное моделирование эффектов дифракции в строгой математической постановке вне рамок теории дифракции Френеля-Кирхгофа. Рассмотрение основывается на решении краевой задачи дифракции для уравнения Гельмгольца, восходящей к трудам А. Зом-мерфельда, Ф. Франка и М. Борна. Представлены расчеты дифракции плоской электромагнитной TE-волны, падающей под произвольным углом на металлическую идеально-проводящую полуплоскость. Определены угловые и радиальные распределения действительной и мнимой частей строгого решения уравнения Гельмгольца. Вычислен вектор Пойнтинга для электромагнитной волны. Определены угловые и радиальные распределения для интенсивности излучения.

Ключевые слова: уравнение Гельмгольца, дифракционный интеграл, краевые дифракционные задачи, численное моделирование дифракции, дифракция на полуплоскости.

A.V. Borovsky A.N. Borodin A.L. Galkin

NUMERICAL SIMULATION OF DIFFRACTION EFFECTS

ssi

2sê Es =

Ép

НА

Ой"

^ 1

5s Ü

«Л*

« р I 8sî ЕнО а

£м§|

МО I

аи-|

нее

«

н

н

о

н

PQ

СО

Н

The article presents numerical simulation of diffraction effects in strict mathematical statement beyond the frameworks of the Fresnel-Kirchhoff diffraction theory. The research is based on the solution of the boundary-value problem of diffraction for the Helmholtz equation, which traces its roots to the works of A. Sommerfeld, F. Frank and M. Born. The article contains calculations of diffraction of flat electromagnetic TE-wave that falls on any angle on a metal idealy conducting half-plane. The authors determine angular and radial distributions of the real and imaginary parts of the rigorous solution of the Helmholtz equation, calculate the Poynting vector for electromagnetic wave, and determine angular and radial distributions for radiation intensity.

Keywords: Helmholtz equation, diffraction integral, boundary-value diffraction problems, numerical simulation of diffraction, diffraction on half-plane.

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

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

© А.В. Боровский, А.Н. Бородин, А.Л. Галкин, 2012

<М*

тН

О

<М

«К1

2нё Яч >

ИдЗ!

нйе

«31

Ен^ I

Ой"

^ 1

5к I

^ Л £ ё§|

« ¡г I

8Н! й§ о

ул* ЕнО а

£м§ц

МО I £И4

НОв

И

н

н

о

н

м

со

Н

наблюдении мало контрастных объектов, где наиболее информативные участки находятся в поле зрения рядом с источниками, имеющими высокую яркость. В связи с этим, приобретают интерес задачи строгой (математической) дифракции на металлических поверхностях. При этом необходимо выйти за рамки приближенной теории Френеля-Кирхгофа, которая неудовлетворительно описывает оптическое поле вблизи металлических поверхностей, а также отражение и рассеяние электромагнитного излучения на них.

Наиболее строгий подход дает применение стационарного волнового уравнения (уравнения Гельмгольца). Для большинства применений достаточным является пространственно двумерное приближение, в рамках которого возможно применение элегантных методов теории функций комплексной переменной для решения уравнения Гельмгольца. Исследование такой задачи носит фундаментальный характер в сфере теоретической и математической физики и восходит к трудам А. Зоммерфель-да [6]. Впоследствии задачами строгой теории дифракции занимались Ф. Франк и Р. Мизес [9], М. Борн и Э. Вольф [1], значительный вклад внес Г. Гринберг [5], раздел, посвященный дифракции на клине имеется в книге Л. Ландау и М. Лифшица [7]. Теорией уравнения Гельмгольца занимались Дж. Грин [10], Дж. Кирхгоф [11], А. Тихонов, А. Самарский [8] и В. Владимиров [4].

Задача дифракции представляет интерес в разных геометриях. В данной статье рассмотрим результаты численного моделирования задачи Зоммерфельда — дифракции на металлической идеально-проводящей полуплоскости.

Физическая постановка задачи Зоммерфельда иллюстрируется (рис. 1). Из бесконечности на металлическую идеально-проводящую полуплоскость под углом у к ней падает плоская электромагнитная волна. Задача изучается в полярных координатах. Радиус-вектор точки наблюдения имеет полярные координаты г, ф. Волновой вектор падающей волны обладает полярными координатами (^ а). Угол падения плоской волны у, который называют в литературе углом Зоммерфельда, равен у = а - п. В задаче для углов 0 < ф < п - у имеется падающая и отраженная волны геометрической оптики, для углов п - у < ф < п + у имеется только падающая (она же прошедшая волна), и для углов п + у < ф < 2п имеет место геометрическая тень. Кроме того, на волны геометрической оптики во всех диапазонах углов накладывается дифрагировавшее (или рассеянное) излучение. Известно [1], что причиной, вызывающей появление отраженной от металлической поверхности волны геометрической оптики, а также рассеянных волн, являются поверхностные токи, индуцированные в пограничном слое металла падающей волной.

Рис. 1. Дифракция на полуплоскости под углом Зоммерфельда у

В работах классиков современной физики [1; 6; 9] получены аналитические решения указанной задачи. Они являются строгими решениями краевой задачи для уравнения Гельмгольца. С точки зрения математики

%

(N

тН

О

(N

SS1

2нё Яч > g||' «31

оК

^ 1

5s s

^ Л £ ё§!

« Р I

8s! й§ о

ЕнО а

£м§|

МО I £И4

нее

К

н

н

о

н

PQ

со

Н

эти решения носят интегральный характер, и их аналитическое исследование крайне затруднено. В частности, в курсе теоретической физики Л.Д. Ландау и М.Е. Лифшица приведена постановка данной задачи, однако, ее решение отсутствует по причине значительной сложности и незавершенности, как отмечают авторы [7]. Применение современных компьютерных и программных средств позволяет провести эффективное численное моделирование решений Зоммерфельда, Франка и Борна.

Начнем изложение с обсуждения уравнений модели. Монохроматические действительные скалярные поля в теории Максвелла в общем случае записываются в виде (вектора отмечаем жирными буквами):

Е(г, ^ = _р(г)соз(о^) + д^вт^),

Е(г, t) = Кеи^ехр^о^)),

U(r) = р(г) + iq(r),

КeU(r) = р(г),

ImU(r) = д^).

Для стационарных задач, в которых функции р и д не зависят от времени, комплексная функция Щг) удовлетворяет уравнению Гельмгольца:

1 _д_ Р др

Р-

dU

+ -

1 э2и

Эр) р2 Эф2

+ U = 0.

Здесь

, 2п

р = дг = — г, А

где X — длина волны. На верхней и нижней поверхностях полуплоскости ставятся условия первого типа

и и = 0, и и = 0;

dU / Эф ф=2п = 0.

или условия второго типа:

ди / Эф |ф=0 = 0,

Остальные условия краевой задачи для уравнения Гельмгольца (условия излучения на бесконечности) учтены в используемом аналитическом решении.

Как показали М. Борн и Э. Вольф [1, с. 517], двумерная задача, когда решение зависит от двух координат р и ф, распадается на линейную суперпозицию двух скалярных задач для волн TE (transverse E) и TH (transverse H) поляризаций. Имеется в виду поперечность электрического или магнитного поля к плоскости падения волны, параллельной волновому вектору и нормали к полуплоскости. Волна TE-поляризации содержит три компоненты поля:

f(0, 0, Ег)

f(Hx Hy0

Из них компонента электрического поля Ez параллельна краю полуплоскости. Компоненты магнитного поля Hx, Hy перпендикулярны краю полуплоскости, причем компонента Hx параллельна полуплоскости, а компонента Hy перпендикулярна ей.

Волна TH-поляризации также содержит три компоненты поля, но других. Причем направления магнитного и электрического векторов меняются местами в сравнении с волной TE-поляризации, описанной выше:

(N

тН

О

(N

SS1

2нё ЯЧ г

MhSI

HSS

«Зі

Ен^ І

Ой"

°i<s ^ 1

5s І

^ Л £ g§! h gl S I

8s!

M g! о

ул*

СнО!

£я§|

E40 f ЙИ-І SfflE

К

н

н

о

н

PQ

со

Н

|(Е ,Е ,0)

тн* у .

1 (0,0, Иг)

Для каждой поляризации возникает скалярная задача для комплексных функций:

Е2(г) = иЕ(т), И2(т) = ин(т).

Граничными условиями для ТЕ волны являются

иЕ ІГ = 0.

Граничными условиями для ТН волны являются

dU

dn

= 0.

В данной статье ограничимся рассмотрением случая ТЕ-волны. Строгое решение задачи дифракции плоской ТЕ-волны на идеально-проводящей полуплоскости имеет вид [1; 6; 9]:

U(p,Ф,у) = U2C(р, ф-у)-U2C(р, ф+у).

Здесь,

exp| і

U2(р, 0) = exp(ip cos0)-

л/П

I ехр(-ід2)^д —

- cos^Vzp

дифракционный интеграл Франка-Зоммерфельда, а у-угол Зоммерфель-да. Эти формулы являются частным случаем более сложной теории, которая рассмотрена авторами в работе [2].

Преобразуем дифракционный интеграл к виду удобному для проведения численных расчетов. Для этого воспользуемся равенством ^ 0 ^

/ ЄХр(-ІД2)^Д = | ЄХр(-ІД2)^Д +1 ехр(-ід2)^д.

0 пт~ 0 ПТ- 0

- сое—у2р -сое—у2р

В первом интеграле в правой части сделаем замену переменных д ^ —д. Второй интеграл вычислим. В результате получим

cos0,/2p

J ехр(-гд2)^д= J ехр(Чд2)ф,+ expf -iП].

Sir- 0 2 V 4 J

-cos^2p v '

Подставим выведенную формулу в исходную формулу для дифракционного интеграла:

exp| іП | cos0,/2p

и£(р, 0) = — ехр(фcos0) + ехр(фcos0)-----| ехр(-гд2)^д. (1)

2 л/п о

Формула (1) удобна для проведения численного моделирования, так как не содержит несобственного интеграла.

Функция (1) допускает выделение действительной и мнимой частей:

ReUC (р,0) = —cos(р cos0) +

2

-у/П

cos] р cos0+nj I cos д2^д + sin ^ р cos0 + nj I sin д2^д

Г

4

1

(?4

тН

О

<М

«К1

2нё ЯЧ г

В«!'

НД

ЬА

ОЙ"

^ 1

5к І

^ Л £ ё§! ь Е

« ^ I

8Н! й§ о

ЕнО а

£м§|

МО I £И4

нее

И

н

н

о

н

м

со

Н

ІтиС(р, 0) = ^іп(рcos0) +

2

л/П

СО80^л/2р

sin| рcos0 + I cosд2^д-сов^рcos0 + | sin

СОЗ^ч/ЇЇр" 2

которые необходимы для вычисления интенсивности излучения.

Вычислим интенсивность излучения. Будем исходить из известного выражения для усредненной по времени плотности потока энергии (вектора Пойнтинга):

5 =— Re(£ • Н). 8п

(2)

Ниже коэффициент — опустим.

8п

Поля Е и Н следующим образом выражаются через решение краевой задачи для уравнения Гельмгольца:

Ег = и,

Нх = совфН,, - этфНф,

Ну = sinфHr + совфНф.

1 дЕ 1 ди

И =

ір Эф ір Эф ’ дЕ , ди

Н = і—^ = і .

ф др др

Вычисление векторного произведения (2) приводит к выражению:

5 = Иє(Е • Н ) = -е* Re(EгИ;) + еу Иє(ЕХ).

Найдем ЕгИ* и ЕН* , выражая поля через решение уравнения Гельмгольца и

ЕгИ* = и^іпфИг + cosфИ )* = sinфU| --1 ди

.віпф ди .

= і-----и—--------1 cosфU

ір дф ди*

+ cosфU| -і

ди*

др

дф

др

ЕгИх =U(cosфИг - sinфИф) = cosфU| -^

1 ди1

ір дф

совф^ди* . ТТди*

= і----и—------+1 віпфи-

- віпфи| -і

.ди*

др

р дф др Вектор Пойнтинга принимает следующий вид:

5 = е Иє

-і

віпфгтди* ттди

■и—--------совфи

Р

дф

др

(

+ е Иє і

у

'ч

со8фгтди* ттди

■V—— + віпфи

Р

дф

др

Взятие реальной части приводит к следующему выражению:

S = е*

+ е.

БіПф

„ ттдІти Т ттдИєУ ИєУ--------Іти-

совф

дф

„ ттдІти Т Г7

ИєУ-----------Іт V

дф

дф

д ИєУ дф

+ совф

„ ттдІти Т ттдИєУ ИєУ--------Іти-

+ БіПф

др

и ттдІти тт

ИєУ-----------Іти

др

др

д ИєУ др

1

*

■к

+

Р

Р

(N

тН

О

(N

SS1

2нё ЯЧ г

MhSI

HSS

«Зі

Ен^ І

Ой"

°i<s ^ 1

5s І

«Л*

g§!

« ¡Ґ I

8s!

M g! о

ул*

СнО!

£м§|

ЙО | £йі

Н®Є

К

н

н

о

н

PQ

со

Н

Интенсивность пропорциональна модулю вектора 5. Для нее получим:

I = 51 = ,

1

р

Reu dlmU - Imu дReU

+

Reu dlmu - Imu d ReU

др

др

. (3)

дф дф

Для расчета интенсивности необходимо определить производные. После достаточно громоздких аналитических вычислений находим:

д

— ReU2 (р, 0) = -cos0 • ImU2 (р, 0) + cos| р - — | cos| —

1

—Ітис(р, 0) = cos0^ ReUf^, 0) - sin | р-П | cos (0| Д—; др V 4 J ^ 2^2пр

д0ReU(р, 0) = р sin0 • ImU2C (р, 0) -^2П cos(р- 4j sin|;

tC,„ o\_ „„.¡„о T3„rrC,„ , І р • ( П 1 • 0

—1тУ„ (p, 0) = -psinS • Rey„ (p, 0) + .— sini p — I sin—.

Э0 2 2 M/2n V 4 J 2

Для определения интенсивности необходимо вычислить Re^(p, ф, у) и 1тЩр, ф, у), а также их производные для функции:

и(р, ф, у) = U(р, ф-у) -U(р, ф +Y),

(4)

1

U2 (р, 0) = — єхр(ір cos0) + єхр(ір cos0)

exp| l^- j cos^

I exp(-i^2)d^. (5)

2 * 4 ' ' " 4 ' ' -у/л о

Опуская подробности аналитических вычислений, выпишем резуль тирующие формулы:

Reи(р, ф, у) = ReUC(р, ф -у) -ReU2C(р, ф +у),

(6)

ImU^, ф, у) = ImU2C(р, ф -у) -Im UC(р, ф +у), д

—ReU^, ф, у) = - cos^ -у )ImU2C (р, ф-у) +

др 2

+ cos^ +у )ImU2 (р, ф +у) + — cos| р— I sin—sin—

4

22

д

—ImU^, ф, у) = cos^ -у )ReU2C (р, ф-у) -др 2 -cos^ +у)ReU2C(р,ф +у)- — sin| р-П|sinфsin—

4

22

—ReU(р, ф, у) = р sin^-у )ImU2C (р, ф-у) -дф 2 -р sin^ +у )ImU|C (р, ф+у) + .)j|— cos( р - n"J cos ф|віп |, д

—ImU(р, ф, у) = -р sin^ -у )ReU2C (р, ф-у) + дф 2 +р sin(ф + у ^eU^^, ф + у) - J— sin | р-П | cos ф sin у.

V п V 4 J 2 2

с4

тН

О

(М

«К1

2Н^ Ма г

н!®1

яЗ|

Ен< |

Ъ<1

и- 1 ё§| «¡л

Он®

15881

^м|| ио г ли» К®«

и н

н

о

н

м

со

н

В данной работе проведено численное моделирование дифракции ТЕ-волны на металлической идеально-проводящей полуплоскости. В основу расчетов положены формулы (3-7). Для проведения расчетов привлекались программные пакеты объектного программирования «Математика» и «MathCad».

В качестве первой задачи был рассчитан дифракционный интеграл Франка-Зоммерфельда. На рис. 2 приведены угловые распределения дифракционных интегралов ReUD(р, 0), ReU2c(р, 0) в диапазоне углов -8п < 0 < 8п для значения р — 1. Оба дифракционных интеграла для задачи дифракции на полуплоскости являются четными функциями, обладающими периодом 4п. Это подтверждают графики этих функций. На рис. 2 так же представлен график дифракционного интеграла Ке и^ (р, 0) для значения р — 10. С увеличением радиуса растет количество максимумов и минимумов в структуре дифракционного интеграла.

В качестве второй задачи было рассмотрено нормальное падение плоской волны. Здесь угол Зоммерфельда равен у — п / 2.

Рис. 2. Дифракционный интеграл Франка-Зоммерфельда — четная функция с периодом 4п: а — Кеи%(р, 0) при р — 1; б — Иеи^р, 0) при р — 1;

в — КеИ^(р, 0) при р — 10

с4

тН

О

(М

«К1

2нё Ма >

М^^!3|

«ЗІ

ёК

н- 1

ё§! ь Е

« V I

§н! Й§ о

и^«

£н0 а

>*М| МО I РнИ-1 К®«

и н

н

о

н

м

со

н

На рис. 3 представлены угловые распределения в полярных координатах действительной части решения уравнения Гельмгольца в диапазоне углов 0 < ф < 2п при различных фиксированных значениях радиуса р = 1, 10, 50.

Рис. 3. Угловое распределение Кеи(ф) при нормальном (у — п / 2) падении плоской волны на полуплоскость в полярных координатах: а — р — 1; б — р — 10; в — р — 50

На рис. 4 представлены угловые распределения действительной части решения уравнения Гельмгольца в диапазоне углов 0 < ф < 2п при различных фиксированных значения радиуса р = 1, 10, 50.

Рис. 4. Угловое распределение Кеи(ф) при нормальном (у — п / 2) падении плоской волны на полуплоскость: а — р — 1; б — р — 10; в — р — 50

с4

тН

О

(М

На рис. 5 представлены угловые распределения модуля решения.

Ма а

Мнз

«51

И- 1

«¡Л Он®

«и

ио г ли» К®«

и

н

н

о

н

м

со

н

1 2 3 4 5 6 ^

в

Рис. 5. Угловое распределение |и(ф)| при нормальном (у — п / 2) падении плоской волны на полуплоскость: а — р — 1; б — р — 10; в — р — 50

На рис. 6 представлены радиальные распределения действительной части решения при фиксированных значениях угла ф — п / 4, п / 2, п.

Яеи

0,4 Л Л Ш Ш Ш Л Ш 1

0,2 ^ 11У11У11 /ШIIУ

20 40 60 80 ,0°

-0,2

-0,4 V V У НННнуНн

Яеи

Рис. 6. Радиальное распределение Кеи(р) при нормальном (у — п / 2) падении плоской волны на полуплоскость: а — ф — п / 4; б — ф — п / 2; в — ф — п

б

с4

тН

О

(М

На рис. 7 представлены угловые распределения интенсивности излучения в диапазоне углов 0 < ф < 2п при различных фиксированных значения радиуса р — 1, 10, 50.

I (Р, Ф, У)

Рис. 7. Угловые распределения интенсивности 1(ф) при нормальном (у — п / 2) падении плоской волны на полуплоскость: а — р — 1; б — р — 10; в — р — 50

«К1 Ез =

Ен< |

Ъ<1

н- 1

ь Е

«¡г I

о

ьм!|

М0|

РнИ-1

К®«

и

н

н

о

н

м

со

н

На рис. 8 представлены радиальные распределения интенсивности излучения в диапазоне радиусов 0 < р < 50 при различных фиксированных значениях угла наблюдения ф — п / 2, п, 3п / 2.

I(р,ф, У)

Рис. 8. Радиальные распределения интенсивности 1(р) при нормальном (у — п / 2) падении плоской волны на полуплоскость: а — ф — п / 2; б — ф — п; в — ф — 3п / 2

В качестве третьей задачи рассмотрено наклонное падение плоской волны на полуплоскость под углом у — п / 3. Результаты численного моделирования представлены на рис. 9.

а

с4

тН

О

(М

«К1 Ма г

МнЗ|

и«!

И- 1

*о!

«¡Л Он®

РбМ|| МО Г ли» К®«

И

н

н

о

н

м

со

н

Рис. 9. Угловое распределение ReU(ф) при наклонном (у — п / 3) падении плоской волны на полуплоскость а — р — 1; б — р — 10; в — р — 50

Результаты численного моделирования эффектов дифракции представляют значительный интерес.

В процессе расчетов определяется комплексное решение уравнения Гельмгольца — функция и(р, ф, у). Физический смысл имеет, прежде всего, действительное значение электрического поля:

Е2(г, £) — Re(U(r)exp(-iюí)) — ReU(r)cos(юí) + 1ти(г^т(ю£).

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

В процессе численного моделирования мы определяли угловые и радиальные распределения для решения уравнения Гельмгольца ReU(r), 1ти(г). Рассмотрены два случая падения плоской волны на полуплоскость — под прямым и наклонным углом к полуплоскости.

Весьма характерным является поведение ReU и |и| как функции угла при больших значениях радиусов, рис. 9в, рис. 5в. На этих графиках отчетливо просматриваются три области:

Для углов 0 < ф < п / 2 имеет место интерференция падающей и отраженной волн. Амплитуда решения изменяется от -2 до 2 для ReU и от 0 до 2 для |и|.

Для углов п / 2 < ф < 3п / 2 имеет место дифракция для падающей и прошедшей волны.

Для углов 3п / 2 < ф < 2п возникает геометрическая тень. В зону геометрической тени происходит диффузия электромагнитного поля из освещенной области.

Второй величиной, представляющей интерес, является интенсивность излучения, для стационарной задачи пропорциональная усредненной по времени объемной плотности энергии электромагнитного поля. Интенсивность излучения в электродинамике Максвелла описывается вектором Пойнтинга (2). Аналитика приводит к выражению (3), в которое также на равных правах входит ReU, 1ти. Не представляет труда показать, что интенсивность невозмущенной плоской волны рав-

<М*

тН

О

<М

«К1

2нё Яч >

ИдЗ!

нйе

«31

Ен^ I

Ой"

^ 1

5к I

^ Л £ ё§|

« ¡г I

8Н! й§ о

ул* ЕнО а

£м§ц

МО I £И4

НОв

И

н

н

о

н

м

со

Н

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

Авторы предлагают ввести характеристику:

8 = ^1 [!(ф) -1]2 dф ,

представляющую собой среднеквадратичное отклонение от единицы углового распределения интенсивности, искаженного дифракцией. Вычисления величины s представлены в таблице.

Среднеквадратичное отклонение дифрагировавшей волны от невозмущенной плоской волны

р У в

2п ~3~ 0,331

3п Т 0,336

30 п 2 0,330

п 3 0,336

п 4 0,341

п 6 0,351

2п ~э~ 0,295

3п Т 0,295

50 п 2 0,290

п 3 0,294

п 4 0,298

п 6 0,303

В данной работе проведено численное моделирование дифракции ТЕ-волны, падающей из бесконечности на идеально-проводящую полуплоскость при нормальном и наклонном падении. В основу математической модели положено строгое решение краевой задачи для уравнения Гельмгольца с условиями Зоммерфельда на бесконечности. Решение краевой задачи для ТЕ-волны выражается в виде разности

(N

тН

О

(N

двух дифракционных контурных интегралов со сдвинутыми углами, которые в случае полуплоскости были сведены А. Зоммерфельдом, Ф. Франком и М. Борном к интегралам Френеля с переменным верхним пределом. В данной работе рассчитаны угловые распределения для дифракционных интегралов. Получены угловые и радиальные распределения для действительной части и модуля решения краевой задачи для уравнения Гельмгольца. Получено аналитическое выражение для интенсивности излучения (пропорциональной усредненной по времени объемной плотности энергии электромагнитного поля). Вычислены угловые и радиальные распределения интенсивности в ближней и дальней зонах дифракции. Рассмотрены случаи нормального и наклонного падения плоской волны на полуплоскость. Введена в рассмотрение и вычислена «дисперсия» интенсивности для прошедшей плоской волны, искаженной эффектом дифракции. Показано, что величина дисперсии возрастает с уменьшением угла Зоммерфельда.

SS1

Яз =

Ин5|

n3i

Ен А

оК

^ 1

Р« S

«Л*

^1

I

8s!

о

ЕнО а

£м§|

МО I £И4

нее

н

н

о

н

PQ

со

н

Список использованной литературы

1. Борн M. Основы оптики / M. Борн, Э. Вольф. — М.: Наука, 1973. — 720 с.

2. Боровский А.В. Дифракционный интеграл Франка-Зоммерфельда / А.В. Боровский, А.Н. Бородин. — Иркутск: Изд-во БГУЭП, 2012. — 64 с.

3. Бородин А.Н. Измерение светорассеяния в зависимости от радиуса кривизны края апертурной диафрагмы / А.Н. Бородин, А.Н. Малов, С.А. Чупра-ков // Компьютерная оптика. — 2005. — Вып. 28. — С. 87-88.

4. Владимиров В.С. Уравнения математической физики / В.С. Владимиров. — М.: Наука, 1982. — 512 с.

5. Гринберг Г.А. К вопросу о дифракции электромагнитных волн на бесконечно-тонких идеально-проводящих экранах / Г.А. Гринберг, Ю.В. Пименов // Журнал теоретической физики. — 1957. — Т. 27. — Вып. 10. — С. 2826-2339.

6. Зоммерфельд А. Оптика / А. Зоммерфельд. — М.: ИЛ, 1953. — 486 с.

7. Ландау Л.Д. Электродинамика сплошных сред: 8 т. / Л.Д. Ландау, Е.М. Лифшиц. — М.: Наука, 1982. — Т. 8. — 621 с.

8. Тихонов А.Н. Уравнения математической физики / А.Н. Тихонов, А.А. Самарский. — М.: Наука, 1972. — 735 с.

9. Франк Ф. Дифференциальные и интегральные уравнения математической физики / Ф. Франк, Р. Мизес. — М.; Л.: ОНТИ, 1937. — 998 с.

10. Green G. Mathematical papers of the late George Green. London: Ferrers’, 1871.

11. Kirchhoff G. Ann. d. physik. 1883. 18. p. 663.

Referenses

1. Born M. Osnovy optiki / M. Born, E. Vol’f. — M.: Nauka, 1973. — 720 s.

2. Borovskii A.V. Difraktsionnyi integral Franka-Zommerfel’da / A.V. Bo-rovskii, A.N. Borodin. — Irkutsk: Izd-vo BGUEP, 2012. — 64 s.

3. Borodin A.N. Izmerenie svetorasseyaniya v zavisimosti ot radiusa krivizny kraya aperturnoi diafragmy / A.N. Borodin, A.N. Malov, S.A. Chuprakov // Komp’yuternaya optika. — 2005. — Vyp. 28. — S. 87-88.

4. Vladimirov V.S. Uravneniya matematicheskoi fiziki / V.S. Vladimirov. — M.: Nauka, 1982. — 512 s.

5. Grinberg G.A. K voprosu o difraktsii elektromagnitnykh voln na beskonechno-tonkikh ideal’no-provodyashchikh ekranakh / G.A. Grinberg, Yu.V. Pimenov // Zhurnal teoreticheskoi fiziki. — 1957. — T. 27. — Vyp. 10. — S. 2826-2339.

6. Zommerfel’d A. Optika / A. Zommerfel’d. — M.: IL, 1953. — 486 s.

7. Landau L.D. Elektrodinamika sploshnykh sred: 8 t. / L.D. Landau, E.M. Lif-shits. — M.: Nauka, 1982. — T. 8. — 621 s.

8. Tikhonov A.N. Uravneniya matematicheskoi fiziki / A.N. Tikhonov, A.A. Samarskii. — M.: Nauka, 1972. — 735 s.

(N

тН

О

(N

9. Frank F. Differentsial’nye i integral’nye uravneniya matematicheskoi fizi-ki / F. Frank, R. Mizes. — M.; L.: ONTI, 1937. — 998 s.

10. Green G. Mathematical papers of the late George Green. London: Ferrers’, 1871.

11. Kirchhoff G. Ann. d. physik. 1883. 18. p. 663.

Информация об авторах

Боровский Андрей Викторович — доктор физико-математических наук, профессор, кафедра информатики и кибернетики, Байкальский государственный университет экономики и права, г. Иркутск, e-mail: [email protected].

Бородин Артур Николаевич — кандидат физико-математических наук, доцент, кафедра информатики и кибернетики, Байкальский государственный университет экономики и права, г. Иркутск, e-mail: [email protected].

Галкин Андрей Леонидович — доктор физико-математических наук, заведующий сектором, Институт общей физики РАН им. А.М. Прохорова, г. Москва, e-mail: [email protected].

Authors

Borovsky Andrey Victorovich — Doctor of Physical and Mathematical Sciences, Professor, Chair of Computer Science and Cybernetics, Baikal State University of Economics and Law, Irkutsk, e-mail: [email protected].

Borodin Artur Nikolaevich — PhD in Physical and Mathematical Sciences, Associate Professor, Chair of Computer Science and Cybernetics, Baikal State University of Economics and Law, Irkutsk, e-mail: [email protected].

Galkin Andrey Leonidovich — Doctor of Physical and Mathematical Sciences, Chair of the Sector, General Physics Institute of RAS named after A.M. Prochorov, Moskow, e-mail: [email protected].

SS?

2s!

Ой"

°¡<s ^ I

5s ¡

«Л i£

«F i

8st

M gl is

EhOH

kB!

^0¡

H

H

o

H P3

со

H