Исследуйте распространение электромагнитной волны в прямоугольном волноводе: как выводятся моды, как зависят пороговые частоты и дисперсия от размеров и материала стенок, какие механизмы потерь ограничивают передачу энергии на миллиметровых и субмиллиметровых длинах волн
Кратко и строго. 1) Вывод мод (основные шаги) - Берём уравнения Максвелла в однородном заполнителе (параметры ϵ,μ\epsilon,\muϵ,μ), предполагаем продольную зависимость exp(−jβz)\exp(-j\beta z)exp(−jβz). При этом поперечные компоненты удовлетворяют уравнению Гельмгольца ∇t2Ψ+kc2Ψ=0,k2=ω2μϵ,kc2=k2−β2,
\nabla_t^2 \Psi + k_c^2 \Psi = 0, \qquad k^2=\omega^2\mu\epsilon,\quad k_c^2=k^2-\beta^2, ∇t2Ψ+kc2Ψ=0,k2=ω2μϵ,kc2=k2−β2,
где Ψ\PsiΨ — искомая проекция поля (EzE_zEz для TM, HzH_zHz для TE). - Для прямоугольного сечения 0<x<a, 0<y<b0<x<a,\;0<y<b0<x<a,0<y<b с идеальными проводящими стенками (PEC) собственные функции и собственные значения: Ψmn(x,y)∝{cos mπxacos nπyb,(TE для соответствующих индексов)sin mπxasin nπyb,(TM)
\Psi_{mn}(x,y)\propto \begin{cases} \cos\!\frac{m\pi x}{a}\cos\!\frac{n\pi y}{b}, & \text{(TE для соответствующих индексов)}\\ \sin\!\frac{m\pi x}{a}\sin\!\frac{n\pi y}{b}, & \text{(TM)} \end{cases} Ψmn(x,y)∝{cosamπxcosbnπy,sinamπxsinbnπy,(TE длясоответствующихиндексов)(TM)
и kc,mn2=(mπa)2+(nπb)2,m,n=0,1,2,…
k_{c,mn}^2=\Big(\frac{m\pi}{a}\Big)^2+\Big(\frac{n\pi}{b}\Big)^2,\qquad m,n=0,1,2,\dots kc,mn2=(amπ)2+(bnπ)2,m,n=0,1,2,…
(для TE допускается m=0m=0m=0 или n=0n=0n=0 с соответствующими формами). - Продольная постоянная: β=k2−kc,mn2.
\beta=\sqrt{k^2-k_{c,mn}^2}. β=k2−kc,mn2. 2) Пороговые частоты и их зависимость - Пороговое (отсечное) значение для режима (m,n)(m,n)(m,n): ωc,mn=kc,mn/μϵ,fc,mn=ωc,mn2π=c2ϵr(ma)2+(nb)2,
\omega_{c,mn}=k_{c,mn}/\sqrt{\mu\epsilon},\qquad f_{c,mn}=\frac{\omega_{c,mn}}{2\pi}=\frac{c}{2\sqrt{\epsilon_r}}\sqrt{\Big(\frac{m}{a}\Big)^2+\Big(\frac{n}{b}\Big)^2}, ωc,mn=kc,mn/μϵ,fc,mn=2πωc,mn=2ϵrc(am)2+(bn)2,
где ccc — скорость света в вакууме, ϵr=ϵ/ϵ0\epsilon_r=\epsilon/\epsilon_0ϵr=ϵ/ϵ0. Для чаще используемого режима TE10_{10}10 при a>ba>ba>b: fc,10=c2aϵr.
f_{c,10}=\frac{c}{2a\sqrt{\epsilon_r}}. fc,10=2aϵrc.
- Следствие: пороговые частоты определяются только геометрией сечения и диэлектрической проницаемостью заполнителя; обратная зависимость от размеров: увеличение aaa снижает fcf_cfc. 3) Дисперсия и скорости - Дисперсионное соотношение для рабочего режима: β=k1−(kc,mnk)2=k1−(fc,mnf)2.
\beta=k\sqrt{1-\Big(\frac{k_{c,mn}}{k}\Big)^2}=k\sqrt{1-\Big(\frac{f_{c,mn}}{f}\Big)^2}. β=k1−(kkc,mn)2=k1−(ffc,mn)2.
- Фазовая и групповая скорости (без потерь): vp=ωβ=cϵr11−(fc/f)2,vg=dωdβ=cϵr1−(fc/f)2,
v_p=\frac{\omega}{\beta}=\frac{c}{\sqrt{\epsilon_r}}\frac{1}{\sqrt{1-(f_{c}/f)^2}},\qquad v_g=\frac{d\omega}{d\beta}=\frac{c}{\sqrt{\epsilon_r}}\sqrt{1-(f_{c}/f)^2}, vp=βω=ϵrc1−(fc/f)21,vg=dβdω=ϵrc1−(fc/f)2,
и vpvg=(c/ϵr)2v_p v_g=(c/\sqrt{\epsilon_r})^2vpvg=(c/ϵr)2. При f→fc+f\to f_c^+f→fc+β→0\beta\to0β→0, vp→∞v_p\to\inftyvp→∞, vg→0v_g\to0vg→0. 4) Влияние материала стенок - Идеальные ПРОВОДНИКИ дают приведённые выше fcf_cfc. Для конечной проводимости сдвиг порога мал, но вводится затухание и комплексная β\betaβ. - Поверхностное сопротивление: Rs=ωμ2σ
R_s=\sqrt{\frac{\omega\mu}{2\sigma}} Rs=2σωμ
(классическая кожная глубина). При высоких частотах и хороших металлах δ≪\delta\llδ≪ размеры, и основной эффект — потери на стенках. 5) Затухание — общие формулы и масштабирование - Общая формула для вкладa потерь в стенках: αc=Rs2∫стенки∣Ht∣2 dlℜ {∬S(E×H∗)⋅z^ dS},
\alpha_c=\frac{R_s}{2}\frac{\displaystyle\int_{\text{стенки}}|H_t|^2\,dl}{\Re\!\left\{\iint_S (\mathbf{E}\times\mathbf{H}^*)\cdot\hat z\,dS\right\}}, αc=2Rsℜ{∬S(E×H∗)⋅z^dS}∫стенки∣Ht∣2dl,
где числитель — мощность, рассеиваемая по длине, знаменатель — передаваемая мощность. - Для диэлектрических потерь (с ϵ=ϵ′−jϵ′′\epsilon=\epsilon'-j\epsilon''ϵ=ϵ′−jϵ′′): αd=ωϵ′′2β∬S∣E∣2dSℜ{∬S(E×H∗)⋅z^ dS} ∝ tanδ.
\alpha_d=\frac{\omega\epsilon''}{2\beta}\frac{\iint_S |E|^2 dS}{\Re\{\iint_S (\mathbf{E}\times\mathbf{H}^*)\cdot\hat z\,dS\}}\;\;\propto\;\;\tan\delta. αd=2βωϵ′′ℜ{∬S(E×H∗)⋅z^dS}∬S∣E∣2dS∝tanδ.
- Масштабирование: Rs∝ωR_s\propto\sqrt{\omega}Rs∝ω ⇒ проводные потери обычно растут как αc∝ω\alpha_c\propto\sqrt{\omega}αc∝ω (плюс фактор 1/β1/\beta1/β, который растёт при приближении к отсечке). Диэлектрические потери ∝tanδ\propto\tan\delta∝tanδ и обычно сильнее выражены при наличии полимеров/стёкол. 6) Механизмы потерь на миллиметровых и субмиллиметровых длинах волн - Омические потери в металле (поверхностный ток): доминируют для обычных металлов; увеличиваются с частотой через RsR_sRs. - Поверхностная шероховатость и неоднородности площадки: усиливают затухание, особенно когда глубина шероховатости сравнима с кожной глубиной. - Аномальный кожный эффект и нечистоты: при очень коротких длинах волн (и низких температурах) средний путь электронов становится сравним с δ\deltaδ — классический RsR_sRs корректируется. - Оксидные / диэлектрические слои на стенках: добавляют диэлектрические потери и рассеяние. - Рассеяние и излучение на стыках, неровностях, углах и отверстиях; утечки в плохо сопряжённых соединениях; в сильно переломанных/многоотверстных линиях — радиационные потери. - Потери в объёме диэлектрика (если заполнение не вакуум): тангенс потерь tanδ\tan\deltatanδ, поглощение из-за водяных паров/полярных молекул в подмм диапазоне. - Модовые перекрёстные и конвертируемые потери: в широкополосных/overmoded волноводах возможна передача энергии в паразитные моды и последующее рассеяние. - При суб-миллиметровых частотах важны квантово-механические и поверхностные эффекты проводника, а также точность механической обработки (требуется наношероховатость). 7) Следствия для практики - Для минимизации потерь на мм/суб‑мм: использовать сверхпроводящую облицовку или высококачественное серебрение, уменьшать число стыков и резких переходов, применять гладкие поверхности, применять диэлектрические волноводы/оптические системы (квазиизоптическая передача) или волоконные/полые волоконные решения для THz. - Пороговые частоты управляются размерами и диэлектриком — если нужно уменьшить fcf_cfc, увеличивают aaa или вводят диэлектрик с большей ϵr\epsilon_rϵr. Краткое резюме: моды получают как собственные решения поперечного Гельмгольца; fcf_cfc определяется геометрией и ϵr\epsilon_rϵr; дисперсия даёт β=k1−(fc/f)2\beta=k\sqrt{1-(f_c/f)^2}β=k1−(fc/f)2, vpv_pvp и vgv_gvg как выше; материалы стенок в основном влияют на затухание (через RsR_sRs, шероховатость, аномальный эффект), а на миллиметровых и субмиллиметровых волнах основные потери — омические (усиленные шероховатостью и аномальным эффектом), диэлектрические и рассеяние на несовершенствах.
1) Вывод мод (основные шаги)
- Берём уравнения Максвелла в однородном заполнителе (параметры ϵ,μ\epsilon,\muϵ,μ), предполагаем продольную зависимость exp(−jβz)\exp(-j\beta z)exp(−jβz). При этом поперечные компоненты удовлетворяют уравнению Гельмгольца
∇t2Ψ+kc2Ψ=0,k2=ω2μϵ,kc2=k2−β2, \nabla_t^2 \Psi + k_c^2 \Psi = 0,
\qquad k^2=\omega^2\mu\epsilon,\quad k_c^2=k^2-\beta^2,
∇t2 Ψ+kc2 Ψ=0,k2=ω2μϵ,kc2 =k2−β2, где Ψ\PsiΨ — искомая проекция поля (EzE_zEz для TM, HzH_zHz для TE).
- Для прямоугольного сечения 0<x<a, 0<y<b0<x<a,\;0<y<b0<x<a,0<y<b с идеальными проводящими стенками (PEC) собственные функции и собственные значения:
Ψmn(x,y)∝{cos mπxacos nπyb,(TE для соответствующих индексов)sin mπxasin nπyb,(TM) \Psi_{mn}(x,y)\propto
\begin{cases}
\cos\!\frac{m\pi x}{a}\cos\!\frac{n\pi y}{b}, & \text{(TE для соответствующих индексов)}\\
\sin\!\frac{m\pi x}{a}\sin\!\frac{n\pi y}{b}, & \text{(TM)}
\end{cases}
Ψmn (x,y)∝{cosamπx cosbnπy ,sinamπx sinbnπy , (TE для соответствующих индексов)(TM) и
kc,mn2=(mπa)2+(nπb)2,m,n=0,1,2,… k_{c,mn}^2=\Big(\frac{m\pi}{a}\Big)^2+\Big(\frac{n\pi}{b}\Big)^2,\qquad m,n=0,1,2,\dots
kc,mn2 =(amπ )2+(bnπ )2,m,n=0,1,2,… (для TE допускается m=0m=0m=0 или n=0n=0n=0 с соответствующими формами).
- Продольная постоянная:
β=k2−kc,mn2. \beta=\sqrt{k^2-k_{c,mn}^2}.
β=k2−kc,mn2 .
2) Пороговые частоты и их зависимость
- Пороговое (отсечное) значение для режима (m,n)(m,n)(m,n):
ωc,mn=kc,mn/μϵ,fc,mn=ωc,mn2π=c2ϵr(ma)2+(nb)2, \omega_{c,mn}=k_{c,mn}/\sqrt{\mu\epsilon},\qquad
f_{c,mn}=\frac{\omega_{c,mn}}{2\pi}=\frac{c}{2\sqrt{\epsilon_r}}\sqrt{\Big(\frac{m}{a}\Big)^2+\Big(\frac{n}{b}\Big)^2},
ωc,mn =kc,mn /μϵ ,fc,mn =2πωc,mn =2ϵr c (am )2+(bn )2 , где ccc — скорость света в вакууме, ϵr=ϵ/ϵ0\epsilon_r=\epsilon/\epsilon_0ϵr =ϵ/ϵ0 . Для чаще используемого режима TE10_{10}10 при a>ba>ba>b:
fc,10=c2aϵr. f_{c,10}=\frac{c}{2a\sqrt{\epsilon_r}}.
fc,10 =2aϵr c . - Следствие: пороговые частоты определяются только геометрией сечения и диэлектрической проницаемостью заполнителя; обратная зависимость от размеров: увеличение aaa снижает fcf_cfc .
3) Дисперсия и скорости
- Дисперсионное соотношение для рабочего режима:
β=k1−(kc,mnk)2=k1−(fc,mnf)2. \beta=k\sqrt{1-\Big(\frac{k_{c,mn}}{k}\Big)^2}=k\sqrt{1-\Big(\frac{f_{c,mn}}{f}\Big)^2}.
β=k1−(kkc,mn )2 =k1−(ffc,mn )2 . - Фазовая и групповая скорости (без потерь):
vp=ωβ=cϵr11−(fc/f)2,vg=dωdβ=cϵr1−(fc/f)2, v_p=\frac{\omega}{\beta}=\frac{c}{\sqrt{\epsilon_r}}\frac{1}{\sqrt{1-(f_{c}/f)^2}},\qquad
v_g=\frac{d\omega}{d\beta}=\frac{c}{\sqrt{\epsilon_r}}\sqrt{1-(f_{c}/f)^2},
vp =βω =ϵr c 1−(fc /f)2 1 ,vg =dβdω =ϵr c 1−(fc /f)2 , и vpvg=(c/ϵr)2v_p v_g=(c/\sqrt{\epsilon_r})^2vp vg =(c/ϵr )2. При f→fc+f\to f_c^+f→fc+ β→0\beta\to0β→0, vp→∞v_p\to\inftyvp →∞, vg→0v_g\to0vg →0.
4) Влияние материала стенок
- Идеальные ПРОВОДНИКИ дают приведённые выше fcf_cfc . Для конечной проводимости сдвиг порога мал, но вводится затухание и комплексная β\betaβ.
- Поверхностное сопротивление:
Rs=ωμ2σ R_s=\sqrt{\frac{\omega\mu}{2\sigma}}
Rs =2σωμ (классическая кожная глубина). При высоких частотах и хороших металлах δ≪\delta\llδ≪ размеры, и основной эффект — потери на стенках.
5) Затухание — общие формулы и масштабирование
- Общая формула для вкладa потерь в стенках:
αc=Rs2∫стенки∣Ht∣2 dlℜ {∬S(E×H∗)⋅z^ dS}, \alpha_c=\frac{R_s}{2}\frac{\displaystyle\int_{\text{стенки}}|H_t|^2\,dl}{\Re\!\left\{\iint_S (\mathbf{E}\times\mathbf{H}^*)\cdot\hat z\,dS\right\}},
αc =2Rs ℜ{∬S (E×H∗)⋅z^dS}∫стенки ∣Ht ∣2dl , где числитель — мощность, рассеиваемая по длине, знаменатель — передаваемая мощность.
- Для диэлектрических потерь (с ϵ=ϵ′−jϵ′′\epsilon=\epsilon'-j\epsilon''ϵ=ϵ′−jϵ′′):
αd=ωϵ′′2β∬S∣E∣2dSℜ{∬S(E×H∗)⋅z^ dS} ∝ tanδ. \alpha_d=\frac{\omega\epsilon''}{2\beta}\frac{\iint_S |E|^2 dS}{\Re\{\iint_S (\mathbf{E}\times\mathbf{H}^*)\cdot\hat z\,dS\}}\;\;\propto\;\;\tan\delta.
αd =2βωϵ′′ ℜ{∬S (E×H∗)⋅z^dS}∬S ∣E∣2dS ∝tanδ. - Масштабирование: Rs∝ωR_s\propto\sqrt{\omega}Rs ∝ω ⇒ проводные потери обычно растут как αc∝ω\alpha_c\propto\sqrt{\omega}αc ∝ω (плюс фактор 1/β1/\beta1/β, который растёт при приближении к отсечке). Диэлектрические потери ∝tanδ\propto\tan\delta∝tanδ и обычно сильнее выражены при наличии полимеров/стёкол.
6) Механизмы потерь на миллиметровых и субмиллиметровых длинах волн
- Омические потери в металле (поверхностный ток): доминируют для обычных металлов; увеличиваются с частотой через RsR_sRs .
- Поверхностная шероховатость и неоднородности площадки: усиливают затухание, особенно когда глубина шероховатости сравнима с кожной глубиной.
- Аномальный кожный эффект и нечистоты: при очень коротких длинах волн (и низких температурах) средний путь электронов становится сравним с δ\deltaδ — классический RsR_sRs корректируется.
- Оксидные / диэлектрические слои на стенках: добавляют диэлектрические потери и рассеяние.
- Рассеяние и излучение на стыках, неровностях, углах и отверстиях; утечки в плохо сопряжённых соединениях; в сильно переломанных/многоотверстных линиях — радиационные потери.
- Потери в объёме диэлектрика (если заполнение не вакуум): тангенс потерь tanδ\tan\deltatanδ, поглощение из-за водяных паров/полярных молекул в подмм диапазоне.
- Модовые перекрёстные и конвертируемые потери: в широкополосных/overmoded волноводах возможна передача энергии в паразитные моды и последующее рассеяние.
- При суб-миллиметровых частотах важны квантово-механические и поверхностные эффекты проводника, а также точность механической обработки (требуется наношероховатость).
7) Следствия для практики
- Для минимизации потерь на мм/суб‑мм: использовать сверхпроводящую облицовку или высококачественное серебрение, уменьшать число стыков и резких переходов, применять гладкие поверхности, применять диэлектрические волноводы/оптические системы (квазиизоптическая передача) или волоконные/полые волоконные решения для THz.
- Пороговые частоты управляются размерами и диэлектриком — если нужно уменьшить fcf_cfc , увеличивают aaa или вводят диэлектрик с большей ϵr\epsilon_rϵr .
Краткое резюме: моды получают как собственные решения поперечного Гельмгольца; fcf_cfc определяется геометрией и ϵr\epsilon_rϵr ; дисперсия даёт β=k1−(fc/f)2\beta=k\sqrt{1-(f_c/f)^2}β=k1−(fc /f)2 , vpv_pvp и vgv_gvg как выше; материалы стенок в основном влияют на затухание (через RsR_sRs , шероховатость, аномальный эффект), а на миллиметровых и субмиллиметровых волнах основные потери — омические (усиленные шероховатостью и аномальным эффектом), диэлектрические и рассеяние на несовершенствах.