Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Маркузе, Д. Оптические волноводы

.pdf
Скачиваний:
45
Добавлен:
22.10.2023
Размер:
17.73 Mб
Скачать

3 7 0 Гаана 8

вины, например 0,5 мкм. Однако внутренняя сордцевниа может быть довольно большой и тем не менее будет воз­ можен одномодовый режим работы, если отношение пл1пг окажется достаточно близким к единице. Можно допустить размер 2at в несколько микронов и при этом все же получить одномодовый волновод.

Оптические волокна имеют некоторые общие черты с полыми металлическими волноводами. И те и другие могут поддерживать ограниченное число направляемых мод на любой заданной частоте. В обеих структурах, воз­ можно преобразование мод, если волноводы отклоняются от идеальной прямолинейной геометрии [82]. Однако, если в металлических волноводах распространяются толь­ ко направляемые моды и преобразование мод связано с обменом мощности между счетным числом этих мод, спектр волн диэлектрического волновода и оптических волокон наряду с конечным числом направляемых мод с дискретным спектром имеет континуум иенаправляемых мод излучения. Моды излучения также являются решения­ ми уравнений Максвелла и удовлетворяют граничным условиям, налагаемым на поле при наличии диэлектри­ ческих поверхностей раздела. Одиако если в металличе­ ских волноводах направляемые моды имеют только дискретный набор постоянных распространения, то в диэлектрическом волокне, кроме них, существует беско­ нечное множество волн излучения с постоянными рас­ пространения, образующими непрерывный спектр. При отклонении от идеальной геометрии происходит не только перераспределение мощности между направляемыми мода­ ми диэлектрического волновода, но и часть мощности рассеивается в непрерывный спектр мод излучения. Рас­ сеяние мощности в непрерывный спектр представляет собой излучение поля во внешнюю среду.

Рассмотрим моды оптического диэлектрического волно­ вода на основе уравнений Максвелла. Поляризационные эффекты оказываются важными для этих структур и в дальнейшем учитываются. Было бы неправильно изучать распространение мод в таком слоистом диэлектрическом волноводе, используя лишь скалярное волновое уравне­ ние. Наряду с изучением свойств регулярных структур мы коснемся в гл. 9 вопросов преобразования мод и эф-

Оптические шит);HU

371

фонтов излучения, вызванных нерегулярностью1) границ волновода.

Математическое рассмотрение потерь на излучение в оптических волокнах с цилиндрической геометрией оказывается довольно сложным. Существуют родствен­ ные структуры, геометрия н математическое рассмотре­ ние которых намного проще, чем цилиндрических опти­ ческих волокон. Такой структурой является плоский

Ф и г . 8.1.2. Продольное сечение плоского волновода.

волновод, пли слой. Схематически плоский волновод показан на фиг. 8.1.2. Предполагается, что волны рас­ пространяются в направлении оси z. Слой неограниченно простирается в направлении осей z и у (перпендикулярно к плоскости фигуры). Как в плоском, так и в круглом оптических волокнах могут распространяться направ­ ляемые моды, но, кроме них, существуют также моды излучения. Явления излучения поля и преобразования мод легче изучать на модели более простого плоского волновода. Результаты, полученные при рассмотрении плоского волновода, непосредственно применимы к круг­ лым оптическим волокнам. Величины потерь на излуче­ ние, вызванные шероховатостью поверхности, почти оди­ наковы для плоского волновода и цилиндрического во­ локна. По этой причине рассмотрим распространение мод в плоском волноводе достаточно подробно, чтобы на его примере изучить свойства диэлектрического оптического волокна. Такой подход позволяет обойтись без сложных математических выкладок и продемонстрировать основные

*) Под нерегулярностью в теории волноводов понимается изменение параметров волновода вдоль ехю оси.— Прим. ред.

24*

'№

I'jiatm 8

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

S.2. НАПРАВЛЯЕМЫЕ ВОЛНЫ (МОДЫ) КРУГЛЫХ ОПТИЧЕСКИХ ВОЛОКОН

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

Ф и г . 8.2.1. Декартова и по­ лярная системы координат.

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

Уравнения поля для оптических систем цилиндриче­ ской геометрии выведены в разд. 1.4. Запишем эти урав­ нения в цилиндрических координатах т, ф, z. Для этого воспользуемся преобразованием координат. Координата z в направлении оптической оси системы является общей для декартовой н цилиндрической систем координат и нет необходимости преобразовывать ее из одной системы в другую. Соотношение между двумя системами коорди­ нат в поперечной плоскости показано на фиг. 8.2.1. Из фигуры видно, что между координатами имеют место сле­

дующие соотношения:

г cos ф,

(8.2.1)

х =

у =

г sin ф.

(8.2.2)

Преобразование компонент Fx, Fy вектора F из декарто­ вой системы координат в компоненты Fr, /Д цилиндриче-

Оптические волокна

373

скоы системы координат показано на фиг. 8.2.2. Мате­ матически это преобразование имеет вид

Fr = Fx c,os Ф+

Fj/Sin ф,

(8.2.3)

Гф= — Fx sin

cos ф.

(8.2,4)

Производные компонент поля по г и ф получаются сле­ дующим образом:

df

дф

d f _ _

'

дх

,

df

ду

дг

дх

дг

ду

дг

д/ дх , 3/ ду дх дф "Т" ду дф

df

, , df . ,

(8.2.5)

 

ду

 

 

df

ду

( 8.2.6)

 

 

Функция / обозначает либо E z, либо H z. С помощью пре­ образований (8.2.3) — (8.2.6) можно в выражениях

Фп г. 8.2.2. Вектор F п его компоненты, параллельные декартовым

иполярным координатам.

(1.4.16) (1.4.19) перейти к цилиндрическим коордипатам:

 

-сор,

1

дН,

)•

(8.2.7)

Е' = - Ы * дг

 

г

дф

 

1

дЕ^

-сор

дПг

)•

(8.2.8)

г

дф

дг

дН,

-сое-

1

дЕг

)■

(8.2.9)

1

дг

 

г

дф

(8. 2. 10)

днг

-сое

дЕ

 

 

 

 

 

 

 

 

дф

 

 

дг )•

 

374

Глава

S

 

Здесь

 

к2 -

р2

(8.2.11)

Xs =

и

=

ю2ец.

(8.2.12)

к2

Параметр (5 является постоянной распространения (волно­ вым числом) в направлении оси z. В итоге мы должны перейти в волновых уравнениях (1.4.22) и (1.4.23) к ци­ линдрическим координатам. Осуществим обратное преоб­ разование в (8.2.1) н (8.2.2)

r = v *2+*/2,

(8.2.13)

ф = arctg-|-

(8.2.14)

и с помощью полученных выражений ВЫЧИСЛИМ первые производные от Ь г

 

dEz

 

dEz

дг .

dEz дф

х

dEz

 

У dEz

 

дх

 

дг

дх

1

дф

дх

г

дг

 

г2 дф

дЕг __

дЕг

дг

.

dEz

дф

__ у

дЕг

,

X dEz

 

 

ду

 

дг

ду

*

дф

ду

 

г

дг

'

Г2 дф

а также

вторые

производные

 

 

 

 

 

 

a2E z

/ 1

г

х 2 \ d E z

1 х (' X д2Е г

а а2я г

 

дх'2-

\

г2 /

дг

1

г [ г

дг2-

 

/*2

аг аФ

 

 

 

 

1 2ху дЕ г

 

У ( х д -Е г

а д^Е

 

 

 

 

1 г4 дф

 

Г2 \ г дг дф

 

 

J - ) ,

 

 

 

 

г2 (Эф2

 

d2E z

/

1

а2 \

dEz

1 у (

У д2Е г .

х

a2£z

 

ду'2-

\ г

г3 / дг

г дг2 ' г2 дг дф ) -

 

 

 

2ху

dE z

1

1

/ у

д2Е г ,

х

a2z?z

 

 

 

 

г‘

дф 1 Г2 \ г дф дг I г2 дф2 )•

 

Подстановка этих выражений в (1.4. 22) приво

 

вому

уравнению в цилиндрической системе

 

 

 

 

a2E z

.

1

дЕ г

,

1

д2Е г> l_v 27?

--Л

 

 

 

 

дг"2-

'

Г

дг

Г2

дф 2

 

 

 

 

Аналогичное уравнение

для

H z следует

из

 

 

 

 

дЧ 1г 1 1 д Н г , 1 д°-Нг

 

 

 

 

 

 

 

дг2 • ^ г дг

Г2 5ф2 4-v2 М —I

 

(8.2.15)

(8.2.16)

(8.2.17)

(8.2.18)

(8.2.19)

Оптические колокиа

375

Волновые уравнения для E Z\\IIZявляются здесь строгими, так как они применяются в областях с однородным пока­ зателем преломления. Итак, получены все уравнения, необходимые для решения задачи о распространении мод в оптическом волокне с оболочкой. Для z-й компоненты электрического поля решение будем искать в виде

Ez= AF(r) eiv*.

( 8

. 2. 21)

Множитель, зависящий от координаты z и времени,

gi(a>t—pz)

(8

.2.22)

умножается па все компоненты поля и поэтому опущен здесь п далее во всех выражениях. Постоянная v в формуле (8.2.21) может быть как положительным, так п отрицатель­ ным числом, но обязательно целым, чтобы имела место периодичность по ф с периодом 2л.

Подстановка выражения (8.2.21) в (8.2.19) приводит к дифференциальному уравнению для F (г)

(8.2.23)

которое является хорошо изученным уравнением Бесселя [11]. Поскольку это дифференциальное уравнение второго порядка, то должны иметь место два линейно независимых решения. Существует несколько способов выбора двух независимых решений (8.2.23). Можно использовать функ­ цию Бесселя J v (кг) п функцию Неймана N v (кг). Эти функции для очень больших действительных значений их аргумента ведут себя подобно косинусам и синусам. Функция Бесселя J v остается конечной в начале коорди­

нат,

тогда

как функция Неймана N v имеет особенность

при /• =

0.

Другой набор линейно независимых решений

состоит

из

функций Ханкеля первого и второго рода

/Д 1’

(к?-) и Щ ’ (кг). Эти функции имеют особенность при

г =

0. Для больших действительных значений аргумента

#vn представляет бегущую волну (экспоненту), которую при выборе времеппой зависимости (8.2.22) рассматри­ вают как волну, распространяющуюся в направлении уменьшения г. Функция Ханкеля второго рода /Д2>

вэтом случае представляет волну, распространяющуюся

внаправлении увеличения г, т. е. от оси структуры.

3 7 li Глава S

Представляют интерес также функции Ханкеля для мни­

мых значений

х

 

(8.2.24)

 

 

х

— iy.

В этом

случае Н(Р (iyr) становится

пропорциональной

e~vr для

больших значений

аргумента,

а IIрр (iyr) пропор­

циональна е^т.

Очевидно,

что только

функция II()1(iyr)

с экспоненциальным спаданием является с физической точки зрения подходящей для описания направляемых мод вне сердцевины волокна. Функция //(?’ (iyr) должна быть отброшена, поскольку она экспоненциально растет с ростом г и поэтому не может правильно описывать рас­ пределение поля волны, которое концентрируется у сердце­ вины волокна.

Все четыре рассмотренные функции н любая их линей­ ная комбинация известны как цилиндрические функций. Если специально не выделять какую-либо из них, то реше­ ние (8.2.23) можно записать как Z (хг). Обозначения, используемые здесь для функций Неймана, такие же, как у Яике и Эмде [11] и Градштейиа и Рыжика 101]. Другие авторы иногда используют для функции Неймана символ Y v. Функцию Ханкеля с мнимым аргументом часто обозначают новым символом. Она пропорциональна так называемой модифицированной функции Ханкеля K v. Различные обозначения цилиндрических функций могут внести путаницу, поэтому более разумно использовать одинаковые символы для одинаковых функций, ясно указывая нх аргументы, вместо того чтобы менять симво­ лы, используемые для самих функций. Каждая из четы­ рех функций может быть выражена как линейная комби­ нация двух других функций. Свойства этих функций изложены в ряде книг, например [8, 11, 61].

Для описания поведения поля во внутренних и внеш­ них областях волокна требуются различные функции. Существование физической направляемой волны накла­ дывает на эти функции, являющиеся решениями (8.2.23), ограничения, которые состоят в том, что внутри сердцеви­ ны при г = 0 они должны оставаться конечными, а вне сердцевины при г оо описывать спадающее поле.

Имея теперь набор всех необходимых величин, можно записать формулы для электромагнитного поля оптического

Оптические волокна

377

волокна. Компоненты, соответствующие

координатам г

и <{>, сразу получаются из 2-й компоненты с помощью

выражений

(8.2.7) — (8.2.10).

 

г < а

Обозначим радиус сердцевины через а. Для

получим

Ez= A JV (кг) еи’ф,

 

(8.2.25)

 

 

 

 

 

Hz — BJV(от) е^ф,

 

(8.2.26)

Ег= -^-fpx/l/;(xr) +

ico|.i0^ - S /v(w )]eiv‘f>I

(8.2.27)

 

X" l

'

J

 

Еф= — ^ [ i P y i / v(w)-xcoj.i0S/v (xr)] е™ф, (8.2.28)

Я г =

— ~ j]— icoe, у A J V

 

(xr)] е”’ф,

(8.2.29)

Яф=

^coeH./v(w)-l-ij3-^-S/v(w)J е;',ф.

(8.2.30)

Штрих означает дифференцирование относительно аргу­ мента кг фупкцин Бесселя (а не ?■). Соотношение между к, Р и ку имеет внд

х2 = А-'(-Р2,

(8.2.31а)

где

(8.2.316)

1с\ = агвщо-

Диэлектрическая постоянная сердцевины связана с пока­ зателем преломления соотношением

 

 

nj = -g-.

(8.2.32)

Поле вне сердцевины, г > а1),

определяется следующим

образом:

Ez= CH?(iyr)e*v+,

'(8-2-33)

 

 

 

 

Hz = DH? (iyr) е«ф,

(8.2.34)

E r = - j r

[$ уС Н Г (1у г )+ щ ю у О Н У (iyr)] e**,

(8.2.35)

£ ф=

_ ^ [ р ^ В Д Д ^ О - у с о ц о Я Я ^ Д п я - ) ]^ ,

(8.2.36)

Я г=

[ - с о е ^ С Я ; 1’ (гуг)-(-урЯЯ^ ’ (iyr] «*♦,

(8.2.37)

Я Ф=

— ^

[уш^СЯ'” ' (iyr) + p

ЯЯ(,1) (iyr)] е**.

(8.2.38)

*) В уравнениях везде опущен множитель (8,2,22),

378 Глава 8

Штрих опять означает дифференцирование относительно аргумента, которым в этом случае является iyr. Величины

у, р и /г2 связаны соотношением

 

Vs = Р2- ^ ,

(8.2.39а)

где

(8.2.396)

A-, =(o2e2(.io-

Постоянные А, В, С п D не определяются уравнениями Максвелла. Выражения для компонент поля, разумеется, являются решениями уравнений Максвелла, но, чтобы они правильно описывали моды оптического волокна [80], они должны еще удовлетворять граничным условиям (1.5.3) и (1.5.4). Так как существуют две тангенциальные компоненты для электрического поля п две для магнитного поля, применение граничных условий дает четыре урав­ нения. Поскольку в выражениях для компонент поля четыре неопределенные константы, то число уравнений совпадает с числом неизвестных. Однако наряду с ампли­ тудными коэффициентами существует постоянная распро­ странения р, которую также необходимо определить. Ее нахождение не представляет трудности. Граничные условия приводят к четырем однородным уравнениям. Однородные уравнения имеют решения только в том случае, если определитель системы уравнений равен нулю. Это дополнительное условие является достаточным для определения постоянной распространения р. Условие равенства определителя пулю называется уравнением собственных значений, так как из него определяются собственные значения р оператора волнового уравнения1). Из уравнения собственных значений определяются соб­ ственные значения постоянных распространения направ­ ляемых мод, пли, что то же самое, волновых чисел этих мод.

Граничные условия позволяют приравнять выражения для компонент поля Ez, Е ф, Н г и Н ф на границе при

J) Правильное сказать: «собственные значения оператора, включающего дифференциальный оператор волнового уравнения и граппчпые условия». Термин «собственное значение» заимствован из математики, точнее, из се раздела, называемого функциональ­ ным анализом, где решается задача о нахождении собственных значений, а также собственных функций линейных операторов. Здесь рассматривается частный случай такой задачи.— Прим, ред.

Оптические волокна

379

г = а. Получается следующая система уравнений:

 

AJ v (na)

 

 

- С И ™ (iya)

= 0,

 

 

 

 

 

 

(8.2.40)

“ 7 -2- ,U V (ха) + i

BJ'V ( я а ) + ^ - СИ*,1' ( i v « ) - ^ DU™' (iTn)=0,

 

 

 

 

 

 

(8.2.41)

 

 

 

 

 

-HHv ' (iva)=0,

 

 

 

 

 

 

(8.2.42)

- i

,u ( , (x a )+ A -£■ HJV ( x a ) + 2 p С Я ' (iv « )+ ^ -

(8.2.43)

 

 

 

 

 

 

Уравнения

(8.2.40) и (8.2.42) связывают коэффициенты

А,

С и В,

D\

J v (xa.)

,

 

 

 

 

 

(8.2.44)

 

 

 

 

U™ (iya)

и

 

 

 

 

 

 

 

J v (ха)

 

(8.2.45)

 

 

 

 

В.

 

 

 

 

II™ (iY“)

 

 

Соотношение

между коэффициентами И и 5

получается

подстановкой

(8.2.44) и

(8.2.45)

в (8.2.43):

 

В

i аху [еп>/;(ха) Н™ (iya)+ ie2x/v {на) Н™ (iya)]

П. (8.2.46)

v”-

 

со (в!—е2)-роРЛ> (на) II™ (1Уа)

 

 

 

Равенство

 

и“-|-у2 = к\

= со" (8i — eg) Цо

(8.2.47)

 

 

 

было использовано для упрощения этого выражения. Уравнение (8.2.41) можно было бы использовать вместо (8.2.43). Полученное при этом выражение для В/А может быть преобразовано к виду (8.2.46) с помощью уравнения собственных значений [см. (8.2.53)]. Уравнение собствен­ ных значений получается из требования, чтобы определи­ тель системы уравнений (8.2.40) —(8.2.43) равнялся нулю:

Jv (на) 0

.соро

^ J r -МжО д а /(,(**)

И

0 Jv (на)

. COEj

v Р

■I----

 

Н

 

-II™

 

 

i ^ tI I ^ (iya)

II™'(iya)

 

У2

 

= 0.

0

—Н™ (iya)

 

v р

(8.2.48)

Соседние файлы в папке книги из ГПНТБ