Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Кураев А.А., Попкова Т.Л., Синицын А.К. Электродинамика и распространение радиоволн.pdf
Скачиваний:
262
Добавлен:
02.05.2014
Размер:
4.84 Mб
Скачать

Волновой КПД (ηB ) определяется через поток энергии, образуемый

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

η

 

(T )=

ν01 Re(j (E&1 (T )H& * (T )E&1 (0)H& * (0)))R0

.

 

 

 

B

 

 

G B2

(1R )

 

 

 

 

0

0

0

 

Сопоставление ηe ,ηB

в отсутствие потерь позволяет оценить точность

численных расчетов.

8.6. Уравнения возбуждения произвольно-нерегулярного коаксиального волновода

В разделе 8.4 cформулированы уравнения возбуждения произвольнонерегулярного полого волновода (односвязная область поперечного сечения). Большой интерес, однако, представляют коаксиальные волноводы, особенно в области миллиметровых и субмиллиметровых волн, где на их основе создаются приборы и устройства, обладающие уникальными характеристиками. Последнее связано с аномальной дисперсией волн Hni (n>>1), позволяющей в нерегулярных коаксиальных волноводах эффективно осуществить селекцию паразитных колебаний и волн в полосе порядка октавы, что открывает путь к созданию сверхразмерных одномодовых коаксиальных структур с рабочей

волной Hni .

Строгой теории нерегулярных коаксиальных структур, однако, не существует; оценки свойств таких структур строятся на базе теории регулярной коаксиальной линии (например, дисперсионного уравнения для такой линии). Двусвязность области поперечного сечения коаксиальной структуры (наличие двух границ в отличие от нерегулярного полого волновода) требует при использовании наиболее естественного для рассматриваемой задачи метода преобразования координат введения новой функции отображения. В данном разделе определена такая функция и на ее основе сформулирована строгая теория произвольно (по z и ϕ) нерегулярной коаксиальной структуры, включая общий случай, когда в ней действуют сторонние негармонические источники.

Рассмотрим произвольно (по z и ϕ) нерегулярный коаксиальный волновод ( r,ϕ, z - компоненты исходной цилиндрической системы координат).

Поверхности внутреннего и внешнего проводников S1,S2 задаются соответственно как b1(ϕ, z) и b2 (ϕ, z) . Задача состоит в определении поля, возбуж-

даемого в волноводе источниками, заданными плотностью стороннего элек-

r

r

δ

 

r

 

r

δ

 

(r,ϕ, z,t)

и плотностью

трического тока δ

= r

r

(r,ϕ, z,t)+ϕ δ

(r,ϕ, z,t)+ z

z

 

0

0

ϕ

0

 

 

стороннего

 

 

 

магнитного

 

 

 

тока

δrМ = rrδ М (r,ϕ, z,t )+ϕr

δ

М (r,ϕ, z,t)+ zr δ М (r,ϕ, z,t).

 

Искомое

поле должно

0 r

 

0

ϕ

0

z

 

 

 

 

99

удовлетворять граничному условию на S1, S2 (потерями в стенках пренебрегаем, σ → ∞ ; изломы S1, S2 отсутствуют):

 

r

r

 

 

S

= 0,

(8.49)

 

 

 

 

n1,2

, E

 

 

 

 

 

 

1,2

 

 

nr1,2 - внешняя нормаль к S1 или S2.

 

Для решения поставленной задачи воспользуемся методом преобразования координат, позволяющим преобразовать граничную задачу (8.49) к элементарной. Введем следующую функцию преобразования:

ρ = (r +b b1) /b ,

(8.50)

b = (b2 b1) /(α 1), α = b2 (0) / b1(0) .

При этом ρ [1,α]

и в новых переменных внутренняя граница волно-

вода регулярна: ρ

 

r =b

=1, ρ

 

r =b =α . Учитывая обратное преобразование

 

 

 

 

1

 

 

2

r = ρb b + b1 , для радиуса-вектора точки во внутренней области волновода в новой системе координат ρ,ϕ, z имеем

r (ρ,ϕ, z) = z0z + (x0 cosϕ + y0 sinϕ)(ρb b +b1) .

(8.51)

Определим основную систему векторов косоугольной системы ρ,ϕ, z :

ar1 = arρ = rr = b(xr0 cosϕ + yr0 sinϕ) = brr0 ,

ρ

ar

= ar

=

 

rr

=[(ρ 1)

b

+

b1

]rr

+ (ρb b + b )ϕr

,

 

ϕ

 

 

2

ϕ

 

 

 

 

ϕ

ϕ 0

1 0

 

ar

= ar

=

r

= zr

+[(ρ 1)

b +

b1

]rr .

 

 

z

 

 

3

z

 

 

0

 

 

 

z

z

0

 

Взаимная система векторов ar1,ar2 ,ar3 определяется следующим обра-

зом:

[ar2,ar3]

 

 

[ar3,a1]

 

 

[a1,a2 ]

 

ar1 =

,

ar2 =

,

ar3 =

,

V

V

 

 

 

 

 

V

 

V = a1[a2,a3] = a2[a3,a1] = a3[a1,a2 ] = b(ρb b + b1) .

Производя указанные действия, имеем

100

 

 

 

 

 

 

 

( ρ 1)

 

b

+ b1

 

 

 

 

 

 

 

 

 

 

 

r

 

 

1 r

 

 

 

r

 

1

 

b

 

b

r

 

 

 

 

 

 

 

ϕ ϕ

 

 

 

 

a1

=

 

 

r

 

 

 

 

 

 

ϕ

 

 

[(ρ 1)

 

+

1

]z

 

,

b

 

b(ρb b + b )

 

b

z

 

 

 

 

0

 

 

 

0

 

 

 

z

 

0

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

r

2

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

=

 

 

 

 

 

 

 

0

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρb

b + b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ar3 = zr0 .

Найдем теперь элементы метрического тензора

 

 

 

 

 

 

 

 

 

 

 

(ρ 1)

b

+

b1

 

 

 

 

 

 

 

 

 

 

 

 

 

r r

 

 

 

1

 

 

 

 

 

 

 

1

 

b

 

b

 

 

 

 

 

 

 

 

 

 

 

 

ϕ ϕ

 

 

 

 

 

 

 

 

g11

= a1a1

=

 

+[

 

 

 

 

 

 

 

 

 

]2

+

 

 

[(ρ 1)

 

+

1

]2

,

b2

 

b(ρb b + b1)

b2

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

g22 = ar2ar2 =1/(ρb b + b )2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g33 = ar3ar3 =1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ρ 1)

b

 

+ b1

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

r

 

 

 

 

1

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

12

1

2

 

 

 

 

 

 

ϕ

 

 

21

 

 

 

 

 

 

 

g

 

= a a

 

 

= −b

 

 

 

 

 

 

 

 

= g

 

 

,

 

 

 

 

 

 

 

 

 

(

ρ

b

+

b1)

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g13 = ar1ar3 = −b1[(ρ 1) bz + bz1 ] = g31 ,

g23 = ar2ar3 = g32 = 0 .

Воспользуемся теперь определением оператора rot в косоугольной системе ρ,ϕ, z :

r

 

1

 

f

3

 

f

2

r

f

1

 

f

3

r

 

f

2

 

f

1

r

 

rotF

=

 

 

{(

 

 

)a + (

 

 

)a

+ (

 

 

)a

},

V

ϕ

z

 

 

ρ

 

 

ϕ

 

 

 

 

1

z

2

 

ρ

3

 

rr

 

 

 

 

rr

 

 

 

rr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1 = (Fa1),

f2 = (Fa2 ),

f3 = (Fa3 ) - ковариантные проекции вектора F .

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

avi arj

= δij .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.52)

Применим (8.52) и найдем ковариантные компоненты уравнений Максвелла в системе ρ,ϕ, z . Однако при записи компонент введем вспомога-

101

тельные векторы Er, Hr,δr,δ M таким образом, чтобы для них компоненты rot имели формальную запись, тождественную выражению их в ортогональной цилиндрической системе координат ρ,ϕ, z . Тогда преобразованные контравариантныеrкомпонентыr уравнений Максвелла для вспомогательных векторов Er, Hr,δ ,δ M в системе переменных ρ,ϕ, z , которые рассматриваются теперь как ортогональные, имеют вид

 

 

r

 

 

 

 

r

 

 

) r

 

 

 

 

 

 

 

 

 

 

) E

 

 

 

 

 

 

 

 

 

 

rotH

=

ε0 g

 

 

 

+ gδ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t r

) r

 

 

 

 

 

(8.53)

 

 

r

 

 

) H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gδ

М

 

 

 

 

 

 

 

rotE

= −µ0 g t

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g11

/ ρ g12

g13 / ρ

 

 

 

 

 

Здесь g) =V

 

 

 

 

 

 

 

 

 

 

 

 

 

g21

 

 

ρg22

g23

,

 

 

 

 

 

 

 

 

 

 

/ ρ g32

g33 / ρ

 

 

 

 

 

 

 

 

 

 

g31

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

r

+

r

 

 

 

r

 

 

 

 

r r

 

r r

Ez

rr

остальные

E

= ρ0Eρ

ϕ0Eϕ′ + z0Ez

, Eρ′ = (E,a1), Eϕ

= (E,a2 ) / ρ,

= (Ea3 ) ,

векторы

 

 

rrrM

 

 

конструируются

аналогичным

образом.

Причем,

 

H

,δ ,δ

 

 

 

 

ρr0,ϕr0, zr0 - тройка ортогональных векторов.

 

 

 

 

Физические векторы в исходной системе r,ϕ, z

рассчитываются через

вспомогательные следующим образом (на примере E):

Er = g11 (g11Eρ + g12ρEϕ′ + g13Ez ) ,

Er = g22 (g21Eρ + g22ρEϕ′ + g23Ez ) ,

Er = g33 (g31Eρ + g32ρEϕ′ + g33Ez ) .

Здесь

g

11

= ar

ar

= b2 , g

22

= ar

ar

=[(ρ 1)

b

+

b1

]2 + (ρb b + b )2

,

 

 

 

 

1

1

 

2

2

 

ϕ

 

ϕ

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g33 = ar3ar3 =1+[(ρ 1) bz + bz1 ]2 .

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

(8.53) в ортогональной системе ρ,ϕ, z при граничных условиях

 

 

r

r

 

 

r

r

 

= 0,

(8.54)

 

 

 

 

 

 

 

 

n1

, E

ρ=1 = 0,

n2

, E

ρ=α

nr1 = rr0, nr2 = −rr0 .

 

 

 

 

 

 

 

102

Прежде чем переходить к решению (8.53), целесообразно для улучшения его сходимости выделить в (8.53) электростатическую часть поля источ-

ников, содержащую разрыв первой производной Eи магнитостатическую, содержащую разрыв H.

При этом динамическая задача имеет вид

 

 

 

 

 

 

r

) r

 

 

 

 

r

) r

 

r

 

= ε

 

) E

r

= −µ

 

) H

(8.55)

rotH

 

 

g

1

+ gδ

, rotE

 

g

t

1 gδ,

 

1

 

0

 

t

1

1

 

0

 

1

 

 

 

r

 

r

 

 

 

 

r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, E1

 

 

 

 

 

 

 

= 0 .

 

 

(8.56)

 

r0

 

 

ρ=1

= 0, r0

, E1

 

ρ=α

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

δr =δr′−ε

0

grad (∂Φ/t ) ,

E = E

+ gradΦ.

 

 

 

 

 

 

 

1

 

 

 

 

 

 

e

1

 

e

 

 

М

=

δ

М

 

 

 

 

 

 

/t),

H1 = H

 

δ1

 

 

µ0 grad(∂ΦМ

 

+ gradΦM ,

- соответственно электрический и магнитный скалярные потенциалы

Φe

,Φм

источников.

 

 

 

 

 

r r

 

 

 

 

 

Существенно, что E1, H1

- непрерывные на границе источников векто-

ры и операция почленного дифференцирования представляющих их в реше-

нии рядов (rot E, H) допустима, поскольку эти ряды сходятся равномерно. Остановимся на решении задачи (8.55),(8.56), полагая режим устано-

вившимся и периодическим. Представим:

Er1t = ReEr&tme jmωt , Er1z = ReEr&zme jmωt , m m

r

I

N

&

e

re

&

м

rм

&

T rT

Etm = ∑ ∑(Amni (z)eni + Amni (z)eni )+ Ame ,

 

i=1n=−N

 

 

 

 

 

 

 

 

r

I

N

 

 

 

 

 

 

 

 

E&zm = ∑ ∑ C&mni (z)ϕniar3.

 

 

 

i=1 n=−N

Здесь и далее индекс «t» обозначает поперечную составляющую соответствующих компонент

r

r

jmωt

r

 

 

r

 

H1t

= ReH&tme , H1z

= ReH&zme jmωt ,

 

m

 

 

 

 

m

 

 

r

I N

 

r

+ B&

 

r

 

r

H&

= ∑ ∑(B&e

 

(z)he

м

(z)h м)+ B&T

(z)hT ,

tm

mni

ni

mni

ni

m

 

 

i=1n=−N

 

 

 

 

 

 

 

103

 

 

 

 

r

 

 

 

 

I

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H&zm = ∑ ∑H&mni (z)ψni zr0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1n=−N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

 

 

 

ϕ

ni

= F e

 

(ρ)ejnϕ , ψ

ni

= F м(ρ)ejnϕ ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

erT = ρr

1/ ρ, hrT =ϕr

1/ ρ, ere

 

=[ρr

 

 

 

e

(ρ) jϕr

 

n

 

F e

(ρ)]ejnϕ ,

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

ni

 

 

0

 

 

ni

 

 

 

 

0 χe

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

erм

= −[ρr

F м(ρ) j

n

 

 

+ϕr

 

 

 

м(ρ)]ejnϕ

, hre

=[ρr

j

n

F e (ρ) +ϕr

 

 

e (ρ)]ejnϕ ,

 

 

F

F

χ м

 

 

 

ni

 

 

 

0

 

ni

 

 

 

 

 

 

0

 

 

ni

 

 

 

 

 

 

 

ni

0

 

 

χe

ni

0

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

 

 

 

 

hrм

=[ρr

 

 

 

м(ρ) ϕr

j

n

 

F м(ρ)]ejnϕ ,

F e (ρ) =

Jn (χnie ρ)

 

Nn (χnie ρ)

 

,

F

 

ni

0

 

 

 

ni

 

0

 

 

м

 

 

ni

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

χni

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Jn (χni )

 

Nn (χni )

 

 

 

e

(ρ) =

Jn(χnie ρ)

Nn(χnie ρ)

, F м(ρ)

=

Jn (χniмρ)

Nn (χniмρ)

,

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

 

 

 

Jn (χnie )

 

 

Nn (

χnie )

ni

 

 

 

 

 

Jn(χniм)

 

 

 

 

Nn(χniм)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м

(ρ) =

Jn(χniмρ)

 

Nn(χniмρ)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ni

 

 

 

Jn(χniм)

 

 

 

 

Nn(χniм)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Собственные числа χniм, χnie определяются следующими дисперсионными уравнениями:

 

Jn(χniм)

=

Jn(αχniм)

,

Jn (χnie )

=

Jn (αχnie )

.

 

 

 

 

 

 

 

 

Nn(χniм) Nn(αχniм) Nn (χnie ) Nn (αχnie )

 

Применяя к решению (8.55), (8.56) проекционную процедуру, имеем

 

2π α

r

r

 

r

r

 

 

∫ ∫{rot(Htm + Hzm )jmωε0 g)(Etm + Ezm )}ereni ρdρdϕ =

(8.57)

01

=1 2π 2π αg)δr1ereniejmωt ρdρdϕdωt ,

π0 0 1

2π α

 

 

r

r

r

r

 

{rot(Htm + Hzm )jmωε0 g)(Etm + Ezm )}erµni ρdρdϕ =

(8.58)

0

1

 

 

 

 

 

 

 

 

 

 

1

2π 2π α

 

) rrM jmωt

 

 

 

=

 

∫ ∫ ∫

gδ e

e

ρd ρdϕdωt ,

 

 

π

 

 

 

 

 

 

 

 

1 ni

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

104

2π α

r

r

r

r

 

∫ ∫{rot(Htm + Hzm )jmωε0 g)(Etm + Ezm )} ϕniar3ρdρdϕ =

(8.59)

01

=1 2π 2π αg)δr1ar3ϕniejmωt ρdρdϕdωt ,

π0 0 1

2π α

r

r

r

r

r

 

∫ ∫{rot(Etm + Ezm )+ jmωµ0 g)(Htm + Hzm )}heni ρdρdϕ =

(8.60)

01

=1 2π 2π αg)δr1мhreniejmωt ρdρdϕdωt ,

π0 0 1

2π α

r

r

r

r

r

 

∫ ∫{rot(Etm + Ezm )+ jmωµ0 g)(Htm + Hzm )}hмni ρdρdϕ =

(8.61)

01

=1 2π 2π αg)δr1мhrмniejmωt ρdρdϕdωt ,

π0 0 1

2π α

r

r

r

r

 

 

∫ ∫{rot(Etm + Ezm )+ jmωµ0 g)(Htm + Hzm )}ψni zr0

ρdρdϕ =

(8.62)

01

=π1 2π 2π αg)δr1мψni zr0ejmωt ρd ρdϕ ,

0

0

1

 

 

 

2π α

r

r

r

r

 

∫ ∫{rot(H&tm + H&zm )jmωε0 g)(E&tm + E&zm )}ρr0dρdϕ =

(8.63)

01

=1 2π 2π 1 g)δr1ρr0dρdϕdωt ,

π0 0 0

2π α

r

r

r

r

 

∫ ∫{rot(E&tm + E&zm )+ jmωµ0 g)(H&tm + H&zm )}ϕr0dρdϕ =

(8.64)

01

=1 2π 2π αg)δr1мϕr0dρdϕdωt .

π0 0 1

Уравнения (8.57)-(8.64) образуют полную систему ОДУ первого порядка, определяющую комплексные амплитуды связанных волн

A&mnie (z), A&mniм (z), B&mnie (z), B&mniм (z),C&mni (z), H&nmi (z), A&mT (z), B&mT (z) , т.е. (8.57)-(8.64)

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

105