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

Lections_V

.pdf
Скачиваний:
11
Добавлен:
20.03.2015
Размер:
1.21 Mб
Скачать

1

 

N

2

 

+

1

 

 

 

N 2 +

1

 

 

 

 

N

2

= 0

 

 

 

 

 

n2

 

 

 

 

 

 

 

 

n2

 

 

 

n2

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

z

.

(12.19)

1− nx2

 

 

 

 

 

 

 

 

 

 

1− ny2

 

 

 

 

 

 

1− nz2

 

 

 

 

 

 

 

 

Следовательно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nx2

 

 

 

 

 

+

 

 

 

Ny2

 

 

 

+

 

 

 

Nz2

 

 

 

 

 

 

= 0

 

 

 

1

 

1

 

 

 

 

1

 

 

1

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(12.20)

 

 

 

n2

nx2

 

 

 

 

 

 

 

 

 

n2

 

ny2

 

 

 

n2

nz2

 

 

 

 

 

Разделим (12.20) на с и учтём vi(P)

=

 

 

c

, следовательно:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

Nx2

 

+

 

 

 

 

 

Ny2

+

 

Nz2

= 0 .

 

 

(12.21)

 

2

 

 

2

 

 

2 2

 

 

 

2

 

 

2

 

 

 

 

 

vP vx

 

 

 

 

 

vP vy

 

 

 

vP vz

 

 

 

 

 

 

 

 

 

 

 

Уравнения (12.20) и (12.21) эквивалентны между собой и называются уравнения-

ми волновых нормалей Френеля, где nx, ny, nz константы данного кристалла. Покажем, что одному и тому же направлению волновой нормали соответствуют

два значения фазовой скорости волны. Запишем уравнение (12.21) в виде:

Nx2 (v2y vP2 )(vz2 vP2 ) + Ny2(vx2 vP2 )(vz2 vP2 ) + Nz2 (vx2 v2p )(vy2 v2p ) = 0(12.

22)

Пусть фронт волны распространяется в x направлении, тогда: Nx = 1, Ny = Nz = 0. Следовательно:

(vy2 vP2 )(vz2 vP2 ) = 0 . Или: vP(1) = vy , v(2)p = vz ; (12.23)

Таким образом, одному направлению волновой нормали соответствуют две ско-

рости v(Py) ¹ v(Pz) , направленные вдоль компонент вектора смещения Dy || vP( y) ,

Dz || v(Pz) .

Точно также как для единичного вектора волновой нормали было получено урав- нение (12.21), можно получить аналогичное уравнение для единичного вектора направ- ления вектора Пойтинга P :

 

 

 

P2

 

 

 

 

Py2

 

 

 

P2

 

 

 

 

 

 

x

 

+

 

 

 

 

 

+

 

 

 

z

 

= 0

 

 

1

1

 

1

1

 

1

1

 

(12.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

v2

v2

 

 

 

v2

v2

 

 

 

 

v2

v2

 

 

 

 

 

Λ

 

x

 

 

 

Λ

 

y

 

 

Λ

 

z

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

. Наша задача связать компоненты p и N .

где vЛ скорость вдоль луча p

Поскольку каждому направлению волновой нормали соответствует две фазовые

скорости v(p1,2) , то следовательно для каждого направления волновой нормали имеется

два направления луча. Каждому из этих лучей соответствует своё направление линей- ной поляризации.

Рассогласование направления колебаний напряжённости и индукции электрического поля в оптических кристаллах

Скорости вдоль волновой нормали могут быть определены вспомогательным методом поверхности волновых нормалей (оптической индикатрисы). Найдём эту поверхность.

Запишем выражения для плотности энергии в кристалле:

 

1

r r

1

æ

D 2

 

D y2

 

D 2

ö

 

w =

 

( D E ) =

 

ç

x

+

 

+

z

÷

(12.25)

 

 

 

 

 

 

ç

ε x

 

ε y

 

ε z

÷ ,

 

 

è

 

 

ø

 

При заданной интенсивности света ω величина постоянная:

D2

Dy2

D2

 

 

 

 

x +

 

+

z = const

12.26)

 

 

 

 

 

 

 

εx

εy

εz

 

 

r

 

r

 

Заменим Dx2

/ const x2 и т.д.,

Z

N

D

 

получим уравнение оптической

 

 

 

индикатрисы (эллипсоид волновых

 

 

 

нормалей):

 

 

 

 

 

 

x2

 

y2

 

z2

 

 

 

 

 

X

2

+

2 +

2 = 1. (12.27)

 

 

 

 

nx

 

ny

nz

 

 

 

 

Чтобы определить показатель

Y

преломления света, распростра-

 

 

няющегося в заданном направлении,

 

 

необходимо рассмотреть сечение

 

 

оптической индикатрисы,

Рис. 12.4. Оптическая индикатриса (эллипсоид

перпендикулярное этому

 

волновых нормалей).

направлению и проходящему через

 

 

начало координат (рис. 12.4). Сечение это представляет собой эллипс, направление главных осей которого указывает направление колебаний вектора D. Показатель пре- ломления равен длине полуоси этого сечения, измеренного в направлении вектора ин-

дукции.

 

Аналогично получаем поверхность для лучей:

 

ε x x2 + ε y y2 + ε z z2 = 1.

(12.28)

Вопрос №38

Одноосные и двухосные кристаллы

Пусть из некоторой точки внутри анизотропного кристалла распространяется свет по разным направлениям. Отложим из этой точки в произвольном направлении отрезки,

Y

Z

v x

v y

v z

X

X

о п т и ч е с к и е о с и

о п т и ч е с к а я о с ь

а )

б )

Р и с . 1 3 .1 . Д в у хо с н ы й ( а ) и о д н о о с н ы й ( б ) к р и с т а л л ы .

равные vst, где vs лучевые скорости в данном направлении, а t – время распростране- ния света внутри кристалла. Геометрические места концов отрезков образуют лучевую поверхность. Двухполостная лучевая поверхность обладает четырьмя точками встречи внешней и внутренней полости. Две линии, соединяющие эти четыре точки, обладают особым свойством, вдоль них свет распространяется с единственной лучевой скоро- стью. Аналогичным образом строится поверхность нормалей, представляющих гео- метрическое место концов отрезков, равных в данном направлении vNt, где vN скоро- сти по нормали. Если скорости вдоль осей x, y, z различны (например, vx>vy>vz), то се- чение лучевой поверхности плоскостью xy будет иметь вид см. рис. 13.1,а. Эллипс

расположен внутри окружности и имеется две оптические оси. Такой кристалл является двухосным.

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

Интерференция поляризованного света в одноосных кристаллах

Запишем уравнение Френеля:

N2x (v2p − v2y )(v2p − v2z ) + N2y (v2p − v2x )(v2p − v2z ) + N2z (v2p − v2x )(v2p − v2y ) = 0

(13.1)

Оптическая ось ориентирована вдоль оси z (рис. 13.1,а) и vx = vy = vo , где vo

ordinary velocity (обыкновенная скорость). Тогда из (13.1) следует:

[N 2x (v 2p - v 2z ) + N 2y (v 2p - v 2z ) + N 2z (v 2p - v 02 )] × (v 2p - v 20 ) = 0 ,

гдеv2z = ve2 ; е – extraordinary velocity (необыкновенная скорость).

(v 2p − vо2 )[(N 2x + N 2y )(v 2p − v e2 ) + N 2z (v 2p − v 02 )] = 0

(13.2)

Обозначим через υ угол между оптической осью z и нормалью N , тогда:

N z2

= cos 2 х; N x2 + N 2y + N z2

= 1;

 

N

N x2

+ N 2y = 1 − cos 2 х; N2x + N 2y = sin 2 х.

 

υ

z

Отсюда:

− vо2 )[(vp2 − ve2 )sin 2 х + (vp2

− vо2 )cos2 х] = 0

 

 

 

 

 

 

(vp2

 

(13.3)

или

о

е

е

о

а)

б)

Ри с . 13 .2 .

О б р азо ван и е э л ли п ти ч ес ки п о ляр и зо ван н о го с вета п р и с ло ж е н и и

 

о б ы кн о вен н о го и н ео б ы кн о вен н о го лу ч ей в кр и ста л ле (а ). П о с тр о е н и е

 

во лн о вы х ф р о н то в дл я о б ы кн о вен н о го и н ео б ы кн о вен н о го лу ч ей (б ) .

vp'2 − vo2 = 0

(13.4) v'p'2

= vo2 cos2 х + ve2 sin 2 х

(13.5)

 

 

Таблицы кристаллов.

 

 

 

Изотропные кристаллы

 

 

 

 

 

вещество

n

 

 

 

 

 

 

CdTe

2,69

 

 

 

 

 

 

NaCl

1,544

 

 

 

 

 

 

Алмаз

 

2,417

 

 

 

 

 

CaAs

3,40

 

 

 

 

 

 

Положительные одноосные кристаллы

 

x

 

вещество

 

no

ne

 

z

Лёд

 

1,309

1,310

 

 

 

 

 

 

 

 

 

Кварц

 

1,544

1,553

 

 

v0 > ve

Рутил

 

2,616

2,903

 

 

 

 

 

 

 

 

Отрицательные одноосные кристаллы

вещество

 

n0

ne

LiNbO3

 

2,300

2,208

BaTiO3

 

2,416

2,364

(NH4 )H2 PO4

(ADP)

1,552

1,478

KH2PO4

(KDP)

1,507

1,467

Двухосные кристаллы.

 

 

вещество

 

n x

n y

Топаз

 

1,619

1,620

Слюда

 

1,552

1,582

4Al2O3 : Cr3+

1,923

1,938

Полевой шпат

1,522

1,526

Е

Р

r

 

E p

 

 

 

к р и с т а л л

 

 

 

x

z

v0 < ve

n z

1,651

1,588

1,947

1,530

А

r

 

E A

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

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

ции.

Y

α

Z

X

h

Р и с . 13 .4 С ло ж ен и е о б ы кн о в ен н о го и н ео б ы кн о в ен н о го луч ей в н утр и кр и с та л ла.

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

лим амплитуды и разность фаз между интерферирующими волнами при произвольных ориентациях плоскостей поляризатора и анализатора и двух взаимно ортогональных направлений Еy и Еx колебаний в кристалле (рис. 13.5). Пройдя кристалл, компоненты поля будут:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

 

 

 

 

cos φe

iδ x

 

 

 

 

D||

 

 

 

 

 

 

 

 

 

 

 

Е0

Р

 

ï E x

= E 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

 

= E

 

sin φe iδ y

 

,

(13.6)

Еy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ï E

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А

î

y

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где разность фаз между

 

 

 

n``

 

 

 

 

 

 

 

 

 

 

 

 

обыкновенной и необыкновенной

 

 

 

 

 

 

 

 

 

 

 

χ

 

 

 

 

δ =

2π

 

¢¢

 

 

¢

)

× h

 

 

 

 

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ (n

 

- n

 

 

 

 

 

 

 

 

 

 

 

 

n`

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проецируя волну на анализатор А,

 

 

 

 

 

 

 

 

 

 

 

 

Еx

D

|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ì

( Ax )

= E

0 e

iδ x

cos φ cos( φ - χ )

 

Рис. 13.6

Проекции обыкновенной и

ïE x

 

iδ y

,

í

( Ay )

= E

 

e

sin φ sin( φ - χ )

 

 

 

 

 

 

необыкновенной составляющей на

ïE

y

 

0

 

 

 

 

 

 

 

 

оси поляризатора и анализатора.

î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

но E( A) = Ex( Ax ) + Ey( Ay ) .

На экране после анализатора наблюдается интерференция лучей:

лучами:

получаем:

(13.7)

I = I1 + I2 + 2

I1I2

cos(δ x −δ y ),

гдеδx −δ y = δ.

 

 

Подставляя в уравнение значения из (13.7) получим:

 

I = E2{cos2

φcos2 (φ - χ) +sin2 φsin2 (φ - χ) + 2cosφsinφcos(φ - χ)sin(φ - χ)cosδ} =

0

 

 

 

 

 

 

 

 

I = E

2 ì

 

2

χ - sin 2φ sin 2(φ - χ)sin

2 δ

ü

(13.8)

ícos

 

2

ý

 

î

 

 

 

 

þ

 

Очевидно, что при δ = 0 для x = 0 и χ = π/2 картина сдвигается на половину по-

лосы.

 

I(δ =0) = E02 cos2 χ закон Малюса.

 

Рассмотрим частные случаи:

 

 

 

 

1) анализатор и поляризатор параллельны: ч = 0.

 

 

I|| = E2

æ

 

 

 

 

 

 

2φ sin2

δ

ö

(13.9)

 

ç1- sin2

2

÷

 

 

р

è

 

 

3

 

 

 

 

ø

 

Если

ц = 0,

, р,

 

р... – направления анализатора совпадают с одной из главных осей

 

2

 

2

 

 

 

 

 

 

 

 

 

кристалла, то результат интерференции:

 

 

I|| = Imax

= E2 .

 

 

 

 

Если ц =

р

 

;

3

р;

5

р... - направления анализатора с кристаллографическими осями

 

 

 

 

4

 

 

 

4

4

 

 

 

 

 

составляет 450, то:

I|| = E

2

æ

- sin

2

δ ö

= E

2

cos

2 δ

 

(13.10)

 

ç1

 

÷

 

2

 

 

 

 

è

 

 

2 ø

 

 

 

 

 

 

 

2) анализатор и поляризатор взаимно перпендикулярны

ч =

 

р

.

2

I = E2 sin 2 2φ sin 2 δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(13.11)

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Интерференция света после одноосного кристалла. Каноническая картина.

Пусть конус сходящихся лучей па- дает на одноосный кристалл, вы- резанный перпендикулярно опти- ческой оси (рис. 13.6). Ось конуса

лучей совпадает с оптической осью кристалла. В этом случае разность фаз между интерфери-

рующими лучами будет зависеть от величины угла падения.

Из уравнений Френеля (13.4) - (13.5) следует:

ϕ

 

χ

 

υ

Z

 

 

 

кристалл

 

Р А

Рис. 13.6. Интерференция света в сходящихся пучках.

(v'p )2 - (v'p' )2 = (v02 - ve2 )sin2 υ

(13.12)

 

 

 

 

 

 

 

c

 

 

 

1

 

 

 

1

 

 

 

æ

1

 

 

1

ö

2

υ

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

ç

 

 

÷

 

, но vp

 

 

 

. Откуда:

 

 

 

-

 

 

 

 

= ç

 

 

-

 

 

÷sin

 

(13.13)

n

 

(n¢)

2

(n¢¢)

2

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

n0

 

 

ne

ø

 

 

 

Поскольку n0

ne , то:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

n2 - n2

(n - n )(n + n )

 

2n

(n - n )

 

 

 

 

 

-

 

 

=

 

 

e

0

=

 

e

o

 

 

e

o

 

»

 

 

 

 

 

 

n2

n2

 

 

 

 

 

 

 

 

n4

 

 

 

 

 

 

 

 

 

n2n2

 

 

n2n2

 

 

 

 

 

 

 

e

o

 

 

 

0

 

e

 

 

 

 

 

o

e

 

 

 

e o

 

 

 

 

 

r

 

 

 

 

 

 

 

и n¢¢ - n¢ = (ne

- n0 ) sin2 υ,

 

 

 

 

 

 

 

 

 

Тогда разность фаз:

 

где δ = (k × r) .

δ = (n

e

 

- n ) sin2

υρ

 

 

 

 

 

(13.14)

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

с sin 2 х = const = c

 

Поверхность равной фазы д = const примет вид:

(13.15)

Что это за поверхность? Исследуем вид такой поверхности. Ниже представлен пример

наиболее распространённого метода анализа интерференционных картин.

 

В декартовых координатах:

 

 

 

 

 

 

 

 

 

ì

 

2

= x

2

+ y

2

+ z

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

í

 

2

 

 

 

 

2

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

(13.16)

 

ï

 

sin

х

= x

+ y

 

 

 

 

 

 

 

 

 

 

 

 

îс

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условия (13.15) с учетом (13.16) следует:

 

 

 

 

 

(x2 + y2 )2

 

= c2 (x2 + y2 + z2 )

 

 

 

 

 

 

 

(13.17)

 

Поверхность (13.17) можно получить при вращении кривых.

 

Пусть у=0, тогда из (13.17):

x 4

= c2 (x 2

+ z2 )

 

 

(13.18)

 

Вращая кривую (13.18) вокруг оси z получим:

 

 

 

 

для z2

>> x2 - уравнение параболы: x2 ≈ cz

 

 

(13.19)

 

для z вблизи оси x получим уравнение гиперболы:

 

 

 

 

2

 

 

 

 

 

 

x4

 

 

 

 

 

 

x2

 

 

2

æ

 

z

2 ö

 

2

 

2

 

 

c =

 

 

2

 

 

 

 

2 =

 

 

 

2

» x

 

ç

-

 

÷

» x - z

 

 

 

x

+ z

 

 

z

 

ç1

x

2 ÷

 

(13.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

è

 

ø

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ x2

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, теперь можно просто изобразить все линии равной разности фаз,

взяв сечения полученной поверхности плоскостями, находящимися на разных рас-

стояниях от начала координат (рис.13.7). Кривые, вдоль которых разность фаз посто-

янная, называется изохроматами. Для нашего случая: изохроматы имеют вид окружно-

стей. Изохромата разрезана по двум взаимно перпендикулярным направлениям.

Кривые, для которых sin 2ц = 0 главные изотермы. Картина интерференции,

представленная на экране имеет название мальтийский крест.

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п а р а б о л а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ги п е р б о л а

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

а )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б )

 

 

Р и с . 1 3 .7 . П о л у ч е н и е к р и в ы х п о в е р х н о с те й р а в н ы х ф а з (а ).

 

 

 

 

 

 

 

 

 

 

М а л ь ти й с к и й кр е с т (б ) .

 

 

 

 

 

Вопрос № 39

Дифракция света.

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

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

волны. Принцип Гюйгенса - Френеля. Интеграл Киргофа.

Световое возмущение в точке Р возникает вследствие суперпозиции вторичных волн, испускаемых поверхностью волнового фронта. Будем рассматривать монохроматичное

скалярное поле: V ( x , y , z , t ) = U ( x , y , z )e iω t

(14.1)

Тогда волновое уравнение принимает вид: (Ñ2 + k2 )U = 0 , где k = 2π / λ .

(14.2),

Уравнение (14.2) – приведенное волновое уравнение или уравнение Гельмгольца. Предположим, что U′ также удовлетворяет волновому уравнению: (Ñ2 + k 2 )U ¢ = 0 (14.3) Умножим (14.2) на U′ , а (14.3) на U и вычтем уравнения друг из друга:

 

 

¢

2

- UÑ

2

U

¢

= 0

 

 

 

 

 

 

 

(14.4)

 

 

 

 

 

 

 

 

 

U

Ñ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окружим точку наблюдения Р (рис. 14.1) гладкой поверхностью S, охватывающий объ-

ём и проинтегрируем U по всему объёму :

 

 

 

 

 

 

 

 

 

ò(U¢Ñ 2 U - UÑ 2 U¢)dV = const

(14.5)

 

 

 

 

 

 

 

 

 

A

n

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (14.3) симметрично относительно коор-

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

динат x,y,z и ему всегда удовлетворяет сферическая

 

 

 

 

r`

 

 

V

 

волна:

 

 

U

 

=

 

 

e ik r

 

 

(14.6)

 

S

 

 

 

P

S`

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Но эта функция в точке Р обращается в беско-

 

 

 

 

 

 

 

 

 

нечность. Чтобы при решении задачи избежать не-

 

Поверхность S окружает току

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 14.1

 

определённостей, окружим точку Р, сферой S дос-

 

наблюденияР. Поверхность S`

 

таточно малого радиуса r′ . А теперь воспользуемся

 

 

для левой части уравнения (14.5) теоремой Грина:

 

принадлежитповерхности S.

 

 

 

 

 

 

 

 

 

 

 

ò(UÑ

2

¢

 

¢

 

2

 

 

 

 

 

æ

U¢

 

-U

¢ U ö

 

 

 

 

 

 

 

 

 

 

U

-U Ñ

U )dV = - ò

çU

n

 

÷dS

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

S +Sè

 

 

n ø

 

 

 

 

 

 

 

 

 

 

 

(14.7), где n внутренняя нормаль по отношению к пространству V. Этот вектор

направлен наружу от объема. Тогда из (14.5) следует: -

æ

U

¢

-

U ¢

U ö

=

ç U

n

 

n

÷dS = 0

 

 

 

 

 

(14.8)

 

 

 

 

 

 

 

 

 

 

S+òS′ è

 

 

 

 

ø

 

 

 

 

¢

 

 

 

 

 

 

U¢

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

 

U

- U¢

U

ö

 

 

 

òS

æ

-U

¢

U

ö

(14.9)

 

 

 

 

 

 

 

çU

n

 

n

÷dS¢ = -

çU

n

n

÷dS

 

 

 

 

 

 

 

Sòè

 

 

 

 

ø

 

 

 

è

 

 

ø

 

 

 

 

 

 

 

 

 

Перейдём в (14.9) к пределу, устремляя радиус r′ → 0 . Правая часть (14.9) не из- меняется. В этом случае:

lim

e

ik r

= 1 , следовательно: U'~

1

é

¶ æ

1 ö

 

1 U ù

 

и: ò êU

 

ç

 

÷

-

 

údS¢ (14.10)

 

 

 

r ¢

 

r'

r ′ → 0

 

 

 

S′ ë

r ¢ è

r ¢ ø

 

r ¢ û

Величина U и

∂U

в окрестности точки Р конечна, и

∂U

медленно меняется в ок-

∂r′

∂r′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рестности Р, то можно приблизительно считать:

 

 

 

 

 

 

 

1

 

U

dS¢ »

 

U

ò

4р

dS = 4рr

U

® 0

при r¢ ® 0

(14.11) Первое сла-

 

 

 

 

 

r¢

 

r¢

Sòr¢ ¶r¢

 

 

 

r¢

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

æ

1

ö

æ

1

ö

 

 

 

1

 

гаемое в (14.10): òU

 

 

 

 

ç

 

÷dS ¢ = òU ç -

 

 

÷dS¢ = -òU

 

 

 

dS ¢

 

 

 

 

r¢

r¢

2

r¢

2

 

 

 

 

S

 

 

 

r¢ è

ø

è

 

ø

 

 

 

 

Поскольку интегрирование проводится по площади сферы радиуса r′ , то r′ не изменя- ется и его можно вынести за знак интеграла. Получаем:

1

 

òdS = -4рr¢

2

 

1

 

 

 

 

 

 

 

 

 

- U

 

 

 

 

 

 

 

U = -4рU

 

 

 

(14.12)

r¢2

 

 

 

r¢2

 

 

 

 

Таким образом (14.9) перепишется в виде:

 

 

 

 

1

æ

æ

 

e

ikr

ö

 

e ikr

 

U

ö

 

U P =

 

 

ç

 

 

ç

 

 

 

 

÷

-

 

 

 

÷

(14.13)Это и есть интеграл Кирхго-

 

 

 

 

 

 

 

 

 

 

 

ò ç U

n

ç

 

 

r

÷

r

 

n

÷ dS

 

 

S è

è

 

 

ø

 

 

ø

 

фа.

Если нам известно распределение поля на поверхности S, то с помощью (14.13) мы всегда можем найти величину поля в точке Р внутри этой поверхности.

Вопрос № 40

Теория дифракции Кирхгофа.

Рассмотрим монохроматическую волну, распространяющуюся от то- чечного источника P 0 и, проходя-

щую через отверстие в экране рис. 14.2. Будем находить поле в токе Р за экраном. Предположим, что размеры

отверстия велики по сравнению с длиной волны.

Рассмотрим интеграл Кирхгофа по замкнутой поверхности S, образован- ной: 1) отверстием А; 2) участком В неосвещённого экрана; 3) частью

большой сферы о с центром в точке Р и радиусом R.

По теореме Кирхгофа имеем:

 

 

 

ξ

 

 

B

 

r

Q

s

R

 

A

 

P

P0

 

B

n

 

 

 

Рис. 14.2. Дифракция сферической волны на

бесконечном непрозрачном экране

 

1

é

ò+ò+ò

ù

ì

¶ æ exp( ikS ) öü

 

exp(ikS ) U ù

 

 

 

U (P) =

 

ê

ú

íU

 

ç

 

÷ý

-

 

údS

(14.14) , где

 

s

s

 

ê

 

ξ

ú

î

n è

øþ

 

n û

 

∂ n

 

 

 

ë A B

û

 

 

 

 

 

 

 

 

 

 

 

производная берётся вдоль внутренней нормали к поверхности. Дело в том, что U и

на А, В, и о никогда точно не известны. Поэтому сделаем ряд приближений:

∂U

∂n

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]