Научная статья на тему 'Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (I. ТЕ-волны)'

Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (I. ТЕ-волны) Текст научной статьи по специальности «Математика»

CC BY
212
37
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УРАВНЕНИЯ МАКСВЕЛЛА / ЗАДАЧА ДИФРАКЦИИ / НЕЛИНЕЙНАЯ КРАЕВАЯ ЗАДАЧА НА СОБСТВЕННЫЕ ЗНАЧЕНИЯ / MAXWELL EQUATION / DIFFRACTION PROBLEM / NONLINEAR BOUNDARY EIGENVALUE PROBLEM

Аннотация научной статьи по математике, автор научной работы — Валовик Дмитрий Викторович

Рассматривается краевая задача для системы уравнений Максвелла, описывающая распространение электромагнитных ТЕ-волн в нелинейном диэлектрическом слое с произвольной нелинейностью. Проблема приводит к нелинейной краевой задаче на собственные значения для обыкновенного нелинейного дифференциального уравнения второго порядка. Получено дисперсионное уравнение для собственных значений задачи (постоянных распространения).

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Валовик Дмитрий Викторович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (I. ТЕ-волны)»

УДК 517.958; 517.927.4

Д. В. Валовик

ЗАДАЧА О РАСПРОСТРАНЕНИИ ЭЛЕКТРОМАГНИТНЫХ ВОЛН В СЛОЕ С ПРОИЗВОЛЬНОЙ НЕЛИНЕЙНОСТЬЮ (I. ТЕ-ВОЛНЫ)1

Аннотация. Рассматривается краевая задача для системы уравнений Максвелла, описывающая распространение электромагнитных ТЕ-волн в нелинейном диэлектрическом слое с произвольной нелинейностью. Проблема приводит к нелинейной краевой задаче на собственные значения для обыкновенного нелинейного дифференциального уравнения второго порядка. Получено дисперсионное уравнение для собственных значений задачи (постоянных распространения).

Ключевые слова: уравнения Максвелла, задача дифракции, нелинейная краевая задача на собственные значения.

Abstract. Electromagnetic propagation problem of TE-waves in a layer with arbitrary nonlinearity is considered. It was shown that the boundary value problem for Maxwell equations is reduced to a nonlinear boundary eigenvalue problem for an ordinary nonlinear differential equation of the 2nd order. Dispersion equation for the eigenvalues of the problem (propagation’s constants) is obtained.

Keywords: Maxwell equation, diffraction problem, nonlinear boundary eigenvalue problem.

Введение

Одна из первых известных автору публикаций по распространению электромагнитных волн в нелинейных волноведущих структурах, в которой задачи распространения волн рассматриваются в строгой электродинамической постановке, - это работа [1], вышедшая в 1972 г. По крайней мере, с этого времени (см., например [2, 3]) проблемы нелинейной оптики и, в частности, проблемы распространения волн в нелинейных волноведущих структурах вызывают неослабевающий интерес. Однако если задачи распространения электромагнитных волн в нелинейных волноведущих структурах рассматривать в строгой электродинамической постановке, то они приводят к нелинейным краевым задачам на собственные значения, для которых на настоящий момент не существует общих методов решения. И даже в случае конкретных проблем решению они поддаются с большим трудом. Например, в работе [4] получено дисперсионное уравнение для собственных значений в случае распространения ТЕ-волн в волноводе с керровской нелинейностью при достаточно малом коэффициенте нелинейности. Однако дисперсионного уравнения для случая произвольного коэффициента до сих пор не получено, несмотря на то, что в течение всего этого периода появлялись публикации по этой проблематике (см. [4]). Эта задача была поставлена еще в работе [1], а ведь с момента ее публикации прошло более 30-ти лет. В то же время это одна из простейших задач для нелинейного волновода. Задача еще более усложняется, если рассматривается распространение ТМ-волн. В этом случае

1 Работа выполнена при поддержке гранта Минобрнауки РФ по ФЦП «Развитие потенциала высшей школы» № 2.1.1/1647.

также удалось получить дисперсионное уравнение для достаточно малых значений коэффициента нелинейности. Автору неизвестны какие-либо строгие математические результаты относительно разрешимости краевых задач для электромагнитных волн, распространяющихся в цилиндрических волноводах с более сложными нелинейностями. К этим задачам тесно примыкают задачи распространения электромагнитных волн в нелинейном слое. Фактически это более простые задачи. Но даже такая задача в случае керровской нелинейности и ТМ-волн долгое время не поддавалась решению и была решена лишь недавно (см. [5, 6]). Для решения задачи использовался подход, разработанный в [5, 6]. Сейчас оказалось, что, используя этот метод, можно рассматривать большинство задач распространения электромагнитных волн в нелинейных слоях с единой точки зрения. Применению этого метода для случая распространения ТЕ-волн в слое с произвольной нелинейностью и посвящена настоящая работа. Задача для ТЕ-волн в слое с керровской нелинейностью была решена в работе [7], однако там использовалась существенно иная техника. Фактически там было построено в явном виде решение получающегося уравнения (в эллиптических функциях), что вообще удается сделать редко.

Исследование таких математических моделей важно, поскольку уже довольно давно [2, 3] известны вещества, в которых наблюдаются нелинейные эффекты третьего и четвертого порядков. Ясно, что с возрастанием мощности излучения (например, лазеров) становятся доступными для экспериментального изучения эффекты все более высоких порядков. С некоторыми новыми результатами и достижениями в различных отраслях нелинейной оптики можно ознакомиться по работам [8-10].

1. Задача дифракции и уравнения Максвелла

Рассмотрим электромагнитные волны, проходящие через однородный, изотропный, немагнитный диэлектрический слой с нелинейностью произвольного типа, расположенный между двумя полупространствами х < 0 и х > h в декартовой системе координат Oxyz . Полупространства заполнены изотропной немагнитной средой без источников и имеют постоянную диэлектрическую проницаемость £1 > £о и £3 > £о соответственно, где £о - диэлектрическая проницаемость вакуума. Считаем, что всюду Ц = |М) , где Мю -магнитная проницаемость вакуума.

Электрическое поле гармонически зависит от времени t и удовлетворяет уравнениям Максвелла

rotH = -7'юеЕ; rotE = 7'юцН, (1)

где Е (х, y, z) и Н (х, y, z) - комплексные амплитуды [1].

Диэлектрическая проницаемость внутри слоя описывается следующим законом:

£ = £2 + £0 f (|E\2),

где f е C[0,h] - произвольная функция (некоторые возможные ограничения на функцию f будут указаны далее); £2 > max(£1, £3)- положительные константы. Будем искать решение уравнений Максвелла во всем пространстве.

Электромагнитное поле Е , H удовлетворяет уравнениям Максвелла (1), условию непрерывности касательных составляющих компонент поля на границе раздела сред х = 0 , х = h и условию излучения на бесконечности: электромагнитное поле экспоненциально затухает при |х| в областях

х < 0 и х > h .

/ \T T

Рассмотрим ТЕ-поляризованные волны Е = (0, Ey ,0) , Н = (х ,0, Hz) ,

где ( • ) - операция транспонирования. Предполагая, что компоненты поля гармонически зависят от z, Hy = Hy (х)е1уг, Ех = Ех (х)е1уг , Ez = Ez (х)el^z, из (1) получаем систему уравнений [1, 5, 6]

1уНх (х)-H'z (х) = -lro£Ey (х),

<-lyEy (х) = 1юцНх (х), (2)

E'y (х) = lropHz (х),

где у - неизвестный спектральный параметр - постоянная распространения электромагнитной волны.

После простейших преобразований из системы (2) получаем

у2Ey (х) - Ey (х) = ro2|M£Ey (х).

2 2

Введем обозначения к =ю Ц£0, М = М0 , и выполним нормировку в соответствии с формулами х = кх, d = kd, у = ~, £ j =— (j = 1, 2, 3). Полу-

dr dS к £0

чаем из последнего уравнения

У2Ey (х)- ^(х) = £Ey (х).

Обозначая Ey (х) = Y(х) и опуская значок тильды, получаем

Y "(х) = y2Y (х )-£Y (х). (3)

Будем искать действительные решения Y(х) для уравнения (3), пола-

2

гая у действительным (так что |Е| не зависит от z) и считая

^1, х < 0,

£2 +f (Y2), 0 < х <h, (4)

£3, х > h.

Также будем полагать, что функция Y (х) дифференцируема в слое так, что Y(х)е C(-«>; + го)п C2 (-^; 0)п C2 (0; h)п C2 (h; + ^).

Будем искать такие у, что max (£1, £3) < у2 < £2 .

£ =

2. Решение системы дифференциальных уравнений

Для £ = £1 в полупространстве х < 0 из (3) и (4) получаем уравнение

Y” = (у2 -£1 )y , его общее решение Y(х) = A1e ^у -£1 х + Ае^у -£1 х , в силу условия на бесконечности получаем

Y (х) = A exp ^х^у2 -£1 j. (5)

Для £ = £3 в полупространстве х > h из (3) и (4) получаем уравнение

Y" = (у2 -£3)y , его общее решение Y(х) = B^e^ h) -£3 + Be (х h)у -£3 , в силу условия на бесконечности получаем

Y(х) = Bexp^-(х-h)^/у2 -£3 j . (6)

Постоянные A и B в (5) и (6) будут определяться граничными условиями. Внутри слоя 0 < х < h система (3) принимает вид

d 2Y d? = 'Г-Е 2

(ї2-E2 - f (y2 ))y . (7)

Будем рассматривать последнее уравнение как систему, пусть У ' (х) = 2 (х), получаем

2'(х ) = (у2-82 - / (2 ))У (х ), (8)

У' (х ) = 2 (х ).

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

Последнее

й2 у2 -е2-/(У2) / 2 I 2\\

— =-^—'-У, далее, 2й2 + (е2 -у2 + /(У2 )УйУ = 0.

уравнение, очевидно, является уравнением в полных дифференциалах. Пусть ¥ (У, 2) - его решение, тогда ¥’г (У, 2 ) = 2 , далее получаем

2 2 ф( 2) 1 , V й ф(у 2) й ф(у 2)

¥(У,2)=т+-Ц Отсюда 2Ф'(У2)=-ф)У и ф) = е2-у2 +

+/(у2), тогда

Y 2

9(Y2)= i (Є2-ї2 + f (ud) •

Yo

Первый интеграл ¥ (У, 2) принимает вид

22 +ф(У2) = С ,

(9)

где С - постоянная интегрирования. Введем новые переменные:

У (х),

:(х) = е2 + У (х), -(х) = -Щx(x),

(10)

из них имеем

2

У2 = (х-е2)/- и 22 =(х-е2)т2/- . (11)

Перейдем в системе (8) к новым переменным, последовательно получа-

ем

(мы ввели обозначение То = е^у2 )

(22) = 2(у2-е2 -/(У2)у2,

(е 2 + У2 ) = 2У2;

2 (т-е2 )тТ-2 (т-е2 )т2 Л = 2у 2 (1 -Т - / ( Т - е2 )

2 Т + 2 2 - -

Т' = 2 (Т-е2 )-. -

- = 2у

-

Окончательно получаем

т' = 2 (т-82 )

-

- = у

Тп -

1+/Ут)'

(12)

- + 3т- 2е 2.

Первый интеграл (9) в новых переменных примет вид

-+ф(т-е2) = С .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

(13)

Уравнение (13), вообще говоря, есть трансцендентное уравнение относительно т. Его решение т = т(-) легко может быть выписано явно лишь

в исключительных случаях.

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

2

-

степеням |Е|; считая, что нелинейность выражается в виде многочлена, мы

просто обрываем соответствующий ряд. Однако можно подбирать и другие функции нелинейности, но так, чтобы выполнялись некоторые условия (они будут приведены далее).

3. Граничные условия и дисперсионное уравнение

Для того чтобы выписать дисперсионное уравнение для постоянных распространения электромагнитных волн, необходимо найти значения — (0),

— (И).

Из непрерывности касательных составляющих полей Е и Н получаем Еу (И + 0 ) = Еу (И - 0) = 4И)= У (И); Еу (0 - 0 ) = Еу (0 + 0 ) = 40)= У (0);

И (И + 0) = Иг (И - 0); Иг (0 - 0 ) = Иг (0 + 0);

У'(И ) = (И + 0) = и(И) = 2 (И); У '(0) = т\И2 (0 - 0) = И(0) = 2 (0).

где постоянная Е^И) считается известной.

Отсюда мы получаем, что В = Е2), А = Е(0), У '(И) = -/у2 - 83В =

= И2), У'(0) = д/у2 -81А = И20. И окончательно, И?И) = -у2 -83Е2) и

И20) = д/у2 -81Е20).

Из всего вышесказанного получаем искомые выражения в виде

Т(0) = Є2 + (40)) ’ Х(Н) = (Еу )) ^2;

е(0))2 + е2 ( е(а))2 + £2

, -у ^ ь2 \^у т ь2

— (0) / 2 > 0, —(И )= Г~1^~< 0. (14)

л/у -81 л/у -8з

Из первого интеграла в форме (13) находим значение постоянной С при х = И :

\2 , „ V Л /, ч\2^

С = (Е^) (у2 -Е3 ) + ф[(Е())

(15)

\ /

Теперь из первого интеграла (13), используя (15), мы можем найти уравнение для т(0):

^(20)-^ ) = ( Е?-))2 (у2-83 ) (Е<И))2 ),

поскольку значение —(0) мы можем выразить через Т(0) (см. формулы (14)), то из последнего уравнения получаем окончательно

(х(°) - £2 Ну2 - £1) + ф(х(°) - £2) = (Е^) (у2 - £3) + Ф

(16)

( 2 ^

Формально из (16) получаем т(°) = О £1,£3,£2,у2,(Е(ИН) .

V )

Ясно, что должно быть т( °)> £2 , поскольку х( ° ) = £2 + У 2( ° )> °. При некоторых условиях на функцию ф корень т( ° )>£2 существует. Например,

если ф(х) = а2и+2х2п+2 + а2„х2п +... + а2х2 + а° - многочлен с неотрицательными коэффициентами. Конечно, можно и по-другому выбирать функцию / так, чтобы функция ф обладала нужным свойством. Однако, по-видимому, многочлен в качестве / с положительными коэффициентами является одной из наиболее общих нелинейностей для обычных материалов.

Кроме того, мы считаем функцию / таковой, что величина

Х° - 1 +

f( х)

> ° . Это заведомо справедливо, если в качестве f выступает

многочлен с положительными коэффициентами.

Рассмотрим систему (12), правая часть второго уравнения при сделанных предположениях (относительно функции f ) неотрицательна, это значит,

что функция ^ - неубывающая, но из формул (14) видно, что ^( ° )> °, а

^( И )< ° . Из этого следует, что функция ^ необходимо имеет точку разрыва.

2 ( т-£2 )

Из первого интеграла (13) имеем ^ =—---------—

. .•. Точками разрыва,

С-ф( Х-£2 ) ^ ^

очевидно, являются нули знаменателя последнего выражения. Причем в этих

* / * \ * I Т1 V.»

точках х = XI х ): ^ . Естественно среди всех корней знаменателя рас-

сматриваемого выражения выбирать только те х* , которые > £2 .

Предположим, что имеется ( N + 1) точка разрыва х°,...,хN на промежутке х е( °, И). Как видно, мы не связываем себя какими-то конкретными

ограничениями на количество нулей знаменателя, таким образом, высказанное предположение является достаточно общим. Позднее будет видно, что N -конечное число.

Из свойств функции Л = л( х) следует, что

л(х° - °) =... = л(XN - °) = +“, л(х° + °) = ... = л(XN + °) = -~ . (17)

-1

Введем обозначение ю = ю( ^) =

У

х° -1 +

f( х)

Л 2

+3х-2£2

X

здесь

мы подразумеваем, что х = х( ^), которое можно получить из первого интеграла (13). Будем искать решение на каждом отрезке [°,х°], [х°,х1 ], ..., [XN,И]:

л (*0)

— | тё л = х + С0, 0 < х < х0;

л(х)

л( х)

| тё л = х + Сі, хі + < х <

л (хг)

л( х)

| тё л = х + сN, XN < х < И.

„л (хм)

(18)

Из уравнений (18), учитывая (17), подставляя х = 0, х = ху+1, х = хN, найдем необходимые константы С1, с2,...,CN+1:

с0 =- | тёл;

л(0)

сі+і = | тёл - хг+і, і = 0, N -1; —^ л( И)

(19)

с^+1 =

С учетом (19) уравнения (18) примут вид

л( хо) +^

| тёл = -х + | тёл, 0<х<хо; л(х) л(0)

л (х) +^

| тё л = х + | тё л- х1+1, х1 < х < х1+1, у = 0, N -1; л( х) -с~

л(х) л(И)

| тё л = х + | тё л- И, XN < х < И.

_л( ^) -~

(20)

Введем обозначение | тёл = Т. Из формулы (20) следует, что

—^

хг+1 - ху = Т > 0 . Отсюда получается, что число точек, в которых функция л обращается в бесконечность, конечно на интервале ( 0, И). То что Т > 0 , можно усмотреть из следующего рассуждения. Поскольку интеграл берется от ограниченной неотрицательной функции и не обращающейся тождественно в нуль (это следует из того, что ни т,^ 0, ни X ^ 0 при х е (0, И)), он не мо-

жет равняться нулю. Теперь, полагая в уравнениях (20) x таковым, чтобы все

интегралы слева обратились в нуль (т.е. x = xo, x = Xj, x = xn ), сложим все

уравнения (20), получим

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

+o л(h)

0 = _x0 + ! rod л + xo + T — x^ +... + x^ _i + T — xN + x^ + ! rod л — h .

^( 0) -oo

И окончательно получаем

л(h) +o л(h)

! rodл+ ! rodл + NT = h или ! rodл + ( N + i)T = h . (21)

-o л( 0) л( 0)

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

Список литературы

1. Eleonskii, P. N. Cylindrical Nonlinear Waveguides / P. N. Eleonskii, L. G. Oga-nes’yants, V. P. Silin // Soviet Physics Jetp. - 1972. - V. 35. - № i. - P. 44-47.

2. Цернике, Ф. Прикладная нелинейная оптика / Ф. Цернике, Дж. Мидвинтер. -М. : Мир, 1976.

3. Шуберт, М. Введение в нелинейную оптику : в 2-х т. / М. Шуберт, Б. Виль-гельми. - М. : Мир, 1973.

4. Смирнов, Ю. Г. Распространение электромагнитных волн в цилиндрических волноводах, заполненных нелинейной средой / Ю. Г. Смирнов, С. Н. Куприянова // Журнал вычислительной математики и математической физики. - 2004. -Т. 44. - № 10. - С. 1850-1860.

5. Валовик, Д. В. Расчет постоянных распространения ТМ-поляризованных электромагнитных волн в нелинейном слое / Д. В. Валовик, Ю. Г. Смирнов // Радиотехника и электроника. - 2008. - T. 53. - № 8. - С. 934-940.

6. Валовик, Д. В. О распространении ТМ-поляризованных электромагнитных волн в нелинейном слое с нелинейностью, выраженной законом Керра / Д. В. Валовик, Ю. Г. Смирнов // Журнал вычислительной математики и математической физики. - 2008. - T. 48. - № 12. - С. 2186-2194.

7. Schurmann, H. W. Reflection and transmission of a plane TE-wave at a lossless nonlinear dielectric film / H. W. Schurmann, V. S. Serov, Yu. V. Shestopalov // Physica D. - 2001. - № 158 (2001). - P. 197-215.

8. Шен, И. Р. Принципы нелинейной оптики / И. Р. Шен. - М. : Наука, 1989.

9. Агравал, Г. Нелинейная волоконная оптика / Г. Агравал. - М. : Мир, 1996.

10. Boyd, R. W. Nonlinear optics / R. W. Boyd. - USA. : Academic Press, 2003.

Валовик Дмитрий Викторович

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

Valovik Dmitry Viktorovich Candidate of physico-mathematical sciences, associate professor, sub-department of mathematics and supercomputer modeling, Penza State University

E-mail: dvalovik@mail.ru

УДК 517.958; 517.927.4 Валовик, Д. В.

Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (I. ТЕ-волны) / Д. В. Валовик // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. -2010. - № 1 (13). - С. 18-27.

i Надоели баннеры? Вы всегда можете отключить рекламу.