Рассмотрите случай, когда тонкий проводник с током помещают в переменное магнитное поле: как по закону Фарадея и уравнениям Максвелла предсказать распределение наведенных токов и поле, и какие эффекты возникают при высокочастотном возбуждении (скин-эффект, излучение)
Кратко — как предсказать распределение наведённых токов и поля и какие эффекты при высоких частотах. 1) Базовые уравнения (в дифференциальной форме) и закон Фарадея: ∇×E=−∂tB(Фарадей),
\nabla\times\mathbf{E}=-\partial_t\mathbf{B}\quad(\text{Фарадей}), ∇×E=−∂tB(Фарадей),∇×H=J+∂tD(Ампер–Максвелл),
\nabla\times\mathbf{H}=\mathbf{J}+\partial_t\mathbf{D}\quad(\text{Ампер–Максвелл}), ∇×H=J+∂tD(Ампер–Максвелл),
в проводнике используем закон Ома в местной форме J=σE.
\mathbf{J}=\sigma\mathbf{E}. J=σE.
Для замкнутой тонкой проволочной петли индуцированная ЭДС равна E=−dΦdt,Φ=∫SB⋅dS,
\mathcal{E}=-\frac{d\Phi}{dt},\qquad \Phi=\int_{S}\mathbf{B}\cdot d\mathbf{S}, E=−dtdΦ,Φ=∫SB⋅dS,
а в простом электрическом описании цепи LdIdt+RI=E,
L\frac{dI}{dt}+RI=\mathcal{E}, LdtdI+RI=E,
(при учёте вынужденного излучения добавляется радиационное сопротивление RradR_{rad}Rrad: LI˙+(R+Rrad)I=EL\dot I+(R+R_{rad})I=\mathcal{E}LI˙+(R+Rrad)I=E). 2) Поле и распределение наведённых токов в проводнике: - Если можно пренебречь перемещением зарядов (∂tD\partial_t\mathbf{D}∂tD мал) и волновыми эффектами (габариты << длины волн), то комбинируя ∇×E=−∂tB\nabla\times\mathbf{E}=-\partial_t\mathbf{B}∇×E=−∂tB, J=σE\mathbf{J}=\sigma\mathbf{E}J=σE и B=μH\mathbf{B}=\mu\mathbf{H}B=μH получаем диффузионное уравнение для поля внутри проводника: ∇2H=μσ ∂tH.
\nabla^2\mathbf{H}=\mu\sigma\,\partial_t\mathbf{H}. ∇2H=μσ∂tH.
При гармоническом возбуждении ∼e−iωt\sim e^{-i\omega t}∼e−iωt это даёт уравнение Гельмгольца с комплексным параметром: ∇2H+iωμσ H=0.
\nabla^2\mathbf{H}+i\omega\mu\sigma\,\mathbf{H}=0. ∇2H+iωμσH=0. - Для плоской границы (полупространство проводника) решение даёт экспоненциальное затухание и фазовый сдвиг. Вводят глубину проникновения (skin depth) δ=2ωμσ.
\delta=\sqrt{\frac{2}{\omega\mu\sigma}}. δ=ωμσ2.
Поле и плотность тока при гармонике распределяются по нормали как J(x)=J0 e−(1+i)x/δ(x≥0),
\mathbf{J}(x)=\mathbf{J}_0\,e^{-(1+i)x/\delta}\quad(x\ge 0), J(x)=J0e−(1+i)x/δ(x≥0),
т.е. амплитуда убывает как e−x/δe^{-x/\delta}e−x/δ и имеет фазовый сдвиг −x/δ-x/\delta−x/δ. - Для тонкого проводника: если радиус a≪δa\ll\deltaa≪δ ток по сечению почти равномерный; если a≳δa\gtrsim\deltaa≳δ — ток концентрируется у поверхности (скин-эффект). В многопроводных системах добавляется эффект близости (proximity effect) — перераспределение тока под воздействием полей соседних проводников. 3) Влияние на сопротивление и импеданс: - Поверхностный импеданс проводника Zs=1+iσδ=(1+i)ωμ2σ.
Z_s=\frac{1+i}{\sigma\delta}=(1+i)\sqrt{\frac{\omega\mu}{2\sigma}}. Zs=σδ1+i=(1+i)2σωμ.
- Для тонкого цилиндрического провода при a≫δa\gg\deltaa≫δ эффективное сопротивление на единицу длины примерно пропорционально Zs/(2πa)Z_s/(2\pi a)Zs/(2πa), т. е. сопротивление растёт с частотой пропорционально ω\sqrt{\omega}ω. 4) Переход к излучению (когда важны волновые эффекты): - Критерий: если габариты проводника (длина петли, отрезок) сравнимы с длиной волны λ=2πc/ω\lambda=2\pi c/\omegaλ=2πc/ω — необходимо решать полные уравнения Максвелла с учётом ретардации. Наведённые токи становятся распределёнными вдоль провода с учётом запаздывания и формируют излучающее поле — антенная теория (решают интегральные уравнения типа уравнения Галёркина / Hallén для тонкого провода). - Последствия: возникает радиационное сопротивление RradR_{rad}Rrad и излучаемая мощность; для короткого диполя (l≪λl\ll\lambdal≪λ) радиационное сопротивление малое и порядка ∼(kl)2\sim (kl)^2∼(kl)2 (где k=2π/λk=2\pi/\lambdak=2π/λ), для резонансных размеров излучение доминирует. - Общая балансная цепочка: наведённая ЭДС питает ток, который теряет энергию на омическую рассеивающую (повышенное из-за скин- и proximity-эффектов) и на излучение: E=LdIdt+(Rohm+Rrad)I.
\mathcal{E}=L\frac{dI}{dt}+(R_{ohm}+R_{rad})I. E=LdtdI+(Rohm+Rrad)I. 5) Практическое предсказание и методы расчёта: - Для низких частот и тонких проводов (габариты << λ\lambdaλ) — квазистатический подход: решить уравнение диффузии (аналитически или численно) для распределения тока; можно использовать эквивалентную цепь с частотно-зависимым R(ω)R(\omega)R(ω) и L(ω)L(\omega)L(ω). - Для высоких частот/резонансных размеров — полный волновой расчёт: метод моментов (MoM) для тонких проводов, численные методы (FDTD, FEM) для сложной геометрии; получают распределение тока, поля и поток рассеиваемой/излучаемой энергии. - Учитывать граничные условия на поверхности проводника: n^×(E2−E1)=0,n^×(H2−H1)=Ks,
\hat{n}\times(\mathbf{E}_2-\mathbf{E}_1)=0,\qquad \hat{n}\times(\mathbf{H}_2-\mathbf{H}_1)=\mathbf{K}_s, n^×(E2−E1)=0,n^×(H2−H1)=Ks,
где Ks\mathbf{K}_sKs — поверхностная плотность тока (при сильном скин-эффекте ток удобно моделировать как поверхностный). Итого: предсказание начинается с Фарадея и уравнений Максвелла + J=σE\mathbf{J}=\sigma\mathbf{E}J=σE. В квази-статике получаем диффузионное уравнение и скин-эффект с глубиной δ=2/(ωμσ)\delta=\sqrt{2/(\omega\mu\sigma)}δ=2/(ωμσ). При росте частоты сопротивление растёт (поверхностный импеданс Zs=(1+i)/(σδ)Z_s=(1+i)/(\sigma\delta)Zs=(1+i)/(σδ)), а при размерах, сравнимых с λ\lambdaλ, на первый план выходит излучение — требуется полный волновой расчёт и ввод радиационного сопротивления.
1) Базовые уравнения (в дифференциальной форме) и закон Фарадея:
∇×E=−∂tB(Фарадей), \nabla\times\mathbf{E}=-\partial_t\mathbf{B}\quad(\text{Фарадей}),
∇×E=−∂t B(Фарадей), ∇×H=J+∂tD(Ампер–Максвелл), \nabla\times\mathbf{H}=\mathbf{J}+\partial_t\mathbf{D}\quad(\text{Ампер–Максвелл}),
∇×H=J+∂t D(Ампер–Максвелл), в проводнике используем закон Ома в местной форме
J=σE. \mathbf{J}=\sigma\mathbf{E}.
J=σE. Для замкнутой тонкой проволочной петли индуцированная ЭДС равна
E=−dΦdt,Φ=∫SB⋅dS, \mathcal{E}=-\frac{d\Phi}{dt},\qquad \Phi=\int_{S}\mathbf{B}\cdot d\mathbf{S},
E=−dtdΦ ,Φ=∫S B⋅dS, а в простом электрическом описании цепи
LdIdt+RI=E, L\frac{dI}{dt}+RI=\mathcal{E},
LdtdI +RI=E, (при учёте вынужденного излучения добавляется радиационное сопротивление RradR_{rad}Rrad : LI˙+(R+Rrad)I=EL\dot I+(R+R_{rad})I=\mathcal{E}LI˙+(R+Rrad )I=E).
2) Поле и распределение наведённых токов в проводнике:
- Если можно пренебречь перемещением зарядов (∂tD\partial_t\mathbf{D}∂t D мал) и волновыми эффектами (габариты << длины волн), то комбинируя ∇×E=−∂tB\nabla\times\mathbf{E}=-\partial_t\mathbf{B}∇×E=−∂t B, J=σE\mathbf{J}=\sigma\mathbf{E}J=σE и B=μH\mathbf{B}=\mu\mathbf{H}B=μH получаем диффузионное уравнение для поля внутри проводника:
∇2H=μσ ∂tH. \nabla^2\mathbf{H}=\mu\sigma\,\partial_t\mathbf{H}.
∇2H=μσ∂t H. При гармоническом возбуждении ∼e−iωt\sim e^{-i\omega t}∼e−iωt это даёт уравнение Гельмгольца с комплексным параметром:
∇2H+iωμσ H=0. \nabla^2\mathbf{H}+i\omega\mu\sigma\,\mathbf{H}=0.
∇2H+iωμσH=0.
- Для плоской границы (полупространство проводника) решение даёт экспоненциальное затухание и фазовый сдвиг. Вводят глубину проникновения (skin depth)
δ=2ωμσ. \delta=\sqrt{\frac{2}{\omega\mu\sigma}}.
δ=ωμσ2 . Поле и плотность тока при гармонике распределяются по нормали как
J(x)=J0 e−(1+i)x/δ(x≥0), \mathbf{J}(x)=\mathbf{J}_0\,e^{-(1+i)x/\delta}\quad(x\ge 0),
J(x)=J0 e−(1+i)x/δ(x≥0), т.е. амплитуда убывает как e−x/δe^{-x/\delta}e−x/δ и имеет фазовый сдвиг −x/δ-x/\delta−x/δ.
- Для тонкого проводника: если радиус a≪δa\ll\deltaa≪δ ток по сечению почти равномерный; если a≳δa\gtrsim\deltaa≳δ — ток концентрируется у поверхности (скин-эффект). В многопроводных системах добавляется эффект близости (proximity effect) — перераспределение тока под воздействием полей соседних проводников.
3) Влияние на сопротивление и импеданс:
- Поверхностный импеданс проводника
Zs=1+iσδ=(1+i)ωμ2σ. Z_s=\frac{1+i}{\sigma\delta}=(1+i)\sqrt{\frac{\omega\mu}{2\sigma}}.
Zs =σδ1+i =(1+i)2σωμ . - Для тонкого цилиндрического провода при a≫δa\gg\deltaa≫δ эффективное сопротивление на единицу длины примерно пропорционально Zs/(2πa)Z_s/(2\pi a)Zs /(2πa), т. е. сопротивление растёт с частотой пропорционально ω\sqrt{\omega}ω .
4) Переход к излучению (когда важны волновые эффекты):
- Критерий: если габариты проводника (длина петли, отрезок) сравнимы с длиной волны λ=2πc/ω\lambda=2\pi c/\omegaλ=2πc/ω — необходимо решать полные уравнения Максвелла с учётом ретардации. Наведённые токи становятся распределёнными вдоль провода с учётом запаздывания и формируют излучающее поле — антенная теория (решают интегральные уравнения типа уравнения Галёркина / Hallén для тонкого провода).
- Последствия: возникает радиационное сопротивление RradR_{rad}Rrad и излучаемая мощность; для короткого диполя (l≪λl\ll\lambdal≪λ) радиационное сопротивление малое и порядка ∼(kl)2\sim (kl)^2∼(kl)2 (где k=2π/λk=2\pi/\lambdak=2π/λ), для резонансных размеров излучение доминирует.
- Общая балансная цепочка: наведённая ЭДС питает ток, который теряет энергию на омическую рассеивающую (повышенное из-за скин- и proximity-эффектов) и на излучение:
E=LdIdt+(Rohm+Rrad)I. \mathcal{E}=L\frac{dI}{dt}+(R_{ohm}+R_{rad})I.
E=LdtdI +(Rohm +Rrad )I.
5) Практическое предсказание и методы расчёта:
- Для низких частот и тонких проводов (габариты << λ\lambdaλ) — квазистатический подход: решить уравнение диффузии (аналитически или численно) для распределения тока; можно использовать эквивалентную цепь с частотно-зависимым R(ω)R(\omega)R(ω) и L(ω)L(\omega)L(ω).
- Для высоких частот/резонансных размеров — полный волновой расчёт: метод моментов (MoM) для тонких проводов, численные методы (FDTD, FEM) для сложной геометрии; получают распределение тока, поля и поток рассеиваемой/излучаемой энергии.
- Учитывать граничные условия на поверхности проводника:
n^×(E2−E1)=0,n^×(H2−H1)=Ks, \hat{n}\times(\mathbf{E}_2-\mathbf{E}_1)=0,\qquad \hat{n}\times(\mathbf{H}_2-\mathbf{H}_1)=\mathbf{K}_s,
n^×(E2 −E1 )=0,n^×(H2 −H1 )=Ks , где Ks\mathbf{K}_sKs — поверхностная плотность тока (при сильном скин-эффекте ток удобно моделировать как поверхностный).
Итого: предсказание начинается с Фарадея и уравнений Максвелла + J=σE\mathbf{J}=\sigma\mathbf{E}J=σE. В квази-статике получаем диффузионное уравнение и скин-эффект с глубиной δ=2/(ωμσ)\delta=\sqrt{2/(\omega\mu\sigma)}δ=2/(ωμσ) . При росте частоты сопротивление растёт (поверхностный импеданс Zs=(1+i)/(σδ)Z_s=(1+i)/(\sigma\delta)Zs =(1+i)/(σδ)), а при размерах, сравнимых с λ\lambdaλ, на первый план выходит излучение — требуется полный волновой расчёт и ввод радиационного сопротивления.