表皮効果の理論 — 高周波電流が導体表面に集中する仕組みを導出する

太い銅線に直流を流すと、電流は断面全体を均一に流れます。ところが、周波数を上げていくと奇妙なことが起こります — 電流が導体の表面付近だけに集中し、内部にはほとんど流れなくなるのです。これが表皮効果（skin effect）です。

この現象は日常にも深く関わっています。電子レンジのマイクロ波（2.45 GHz）が金属の壁で反射されるのは、電磁波が金属内部にほとんど侵入できないためです。高周波回路の設計者が導線の太さを太くしても抵抗が下がらないと嘆くのも、表皮効果が原因です。送電線が中心を中空にしても性能が変わらないのも、高周波成分が表面しか流れないことを利用しています。

表皮効果を理解すると、以下のような分野でその知識が直接活きます。

• 高周波回路設計: GHz帯のマイクロストリップ線路やコプレーナ導波路で、導体損失を正確に見積もれるようになる
• 電磁シールド: 金属筐体でどの程度の遮蔽効果が得られるかを、周波数と材料から定量的に評価できる
• 導波管工学: 導波管の壁面損失を表面抵抗から計算し、伝送効率を設計段階で予測できる
• 電力工学: 送電線やバスバーにおける交流抵抗の増大を正しく見積もり、発熱や効率低下を防げる

本記事の内容

• 表皮効果の直感的な理解 — なぜ高周波ほど電流が表面に集中するのか
• マクスウェル方程式から拡散方程式を導出する過程の詳細
• 侵入深さ（スキンデプス） $\delta = \sqrt{2/(\omega\mu\sigma)}$ の導出と物理的意味
• 各金属・各周波数での侵入深さの具体的な数値
• 表面抵抗 $R_s$ の導出と高周波回路への影響
• 電磁シールドおよび導波管の壁面損失への応用
• Pythonによる侵入深さの周波数依存性と電流分布の可視化

前提知識

この記事を読む前に、以下の記事を読んでおくと理解が深まります。

• マクスウェル方程式の全体像 — 4つの方程式の導出と意味
• 電磁波の導出と伝搬特性 — 波動方程式と平面波

表皮効果の直感的な理解

電流が表面に追いやられるメカニズム

表皮効果を直感的に理解するには、導体内部で何が起きているかを段階的に考えてみましょう。

導体に交流電流を流すと、その電流は時間変化する磁場を作ります。ファラデーの法則により、時間変化する磁場は電場（起電力）を誘導します。この誘導された電場が、もとの電流の変化を打ち消す方向に作用します。つまり、レンツの法則の帰結として、導体内部ほど「電流を邪魔する誘導電場」が強くなるのです。

日常的なアナロジーで考えてみましょう。満員電車で人がドアから乗り込むとき、入口付近は人が詰まっていますが、奥はまだ空いています。ところが全員が奥まで進む前にドアが閉まり、また開いて次の波が入ってきます。ドアの開閉（周波数）が速いほど、人は入口付近に集中し、奥までたどり着けません。表皮効果もこれと似ています。電流の「波」が導体表面から内部へ浸透する速度と、交流の反転する速度の競争です。周波数が高いほど反転が速く、電流は表面にとどまります。

もう少し物理的に言えば、導体内部の電磁場は拡散方程式に従います。拡散方程式は、信号が距離とともに指数関数的に減衰することを意味しています。この減衰の特性長が「侵入深さ」$\delta$ であり、表面からの距離 $z$ での電流密度は $e^{-z/\delta}$ に比例して急激に小さくなります。

周波数が高いほど浅くなる理由

拡散のスピードは材料の導電率 $\sigma$ と透磁率 $\mu$ で決まり、有限の速さです。一方、交流の反転は周波数 $f$ とともに速くなります。拡散が内部に浸透しきる前に電流の向きが反転するので、高周波ほど電流が浸透できる深さが浅くなるのです。

定量的には、侵入深さは周波数の平方根に反比例します。つまり周波数を4倍にすると、侵入深さは半分になります。50 Hzの交流電力では銅の侵入深さは約9.3 mmもありますが、1 GHzのマイクロ波では約2.1 $\mu$m — 髪の毛の太さの30分の1にまで薄くなります。

この直感をしっかり持った上で、マクスウェル方程式から厳密に導出していきましょう。

マクスウェル方程式からの拡散方程式の導出

導体内部のマクスウェル方程式

真空中の電磁波は波動方程式に従いますが、導体内部では事情が異なります。導体には自由電子が大量に存在し、電場 $\bm{E}$ が加わると電流密度 $\bm{J}$ が流れます。この関係はオームの法則の微分形で記述されます。

$$ \bm{J} = \sigma \bm{E} $$

ここで $\sigma$ は導電率（単位: S/m）です。良導体の銅では $\sigma \approx 5.96 \times 10^7$ S/m と非常に大きな値を持ちます。

導体内部で考えるべきマクスウェル方程式は次の4つです。ここでは線形等方性媒質を仮定し、$\bm{B} = \mu \bm{H}$、$\bm{D} = \varepsilon \bm{E}$ とします。

$$ \nabla \cdot \bm{E} = \frac{\rho}{\varepsilon} $$

$$ \nabla \cdot \bm{B} = 0 $$

$$ \nabla \times \bm{E} = -\frac{\partial \bm{B}}{\partial t} $$

$$ \nabla \times \bm{H} = \bm{J} + \frac{\partial \bm{D}}{\partial t} = \sigma \bm{E} + \varepsilon \frac{\partial \bm{E}}{\partial t} $$

良導体近似 — 変位電流を無視できる条件

アンペール・マクスウェルの法則の右辺には、伝導電流 $\sigma \bm{E}$ と変位電流 $\varepsilon \partial \bm{E}/\partial t$ の2つの項があります。時間調和場 $\bm{E} \propto e^{j\omega t}$ を仮定すると、変位電流の大きさは $\omega \varepsilon |\bm{E}|$ となります。伝導電流の大きさ $\sigma |\bm{E}|$ との比は

$$ \frac{\omega \varepsilon}{\sigma} $$

です。銅の場合、$\sigma = 5.96 \times 10^7$ S/m、$\varepsilon \approx \varepsilon_0 = 8.854 \times 10^{-12}$ F/m なので、周波数 $f = 10$ GHz（$\omega = 2\pi \times 10^{10}$）でも

$$ \frac{\omega \varepsilon_0}{\sigma} = \frac{2\pi \times 10^{10} \times 8.854 \times 10^{-12}}{5.96 \times 10^7} \approx 9.3 \times 10^{-9} $$

と極めて小さな値です。したがって、通常の良導体では変位電流を無視できます。この近似を良導体近似（good conductor approximation）と呼びます。

$$ \nabla \times \bm{H} \approx \sigma \bm{E} $$

この近似のもとでは、電磁場は波動方程式ではなく拡散方程式に従います。波動的な性質（振動しながら伝搬）よりも拡散的な性質（指数減衰しながら浸透）が支配的になるのです。

電場に対する拡散方程式の導出

ゴールは、導体内部の電場 $\bm{E}$ が満たす偏微分方程式を導くことです。

ファラデーの法則の両辺に回転（$\nabla \times$）を作用させます。

$$ \nabla \times (\nabla \times \bm{E}) = -\frac{\partial}{\partial t}(\nabla \times \bm{B}) $$

右辺の $\nabla \times \bm{B}$ にアンペール・マクスウェルの法則（良導体近似）$\nabla \times \bm{B} = \mu \sigma \bm{E}$ を代入します。ここで $\bm{B} = \mu \bm{H}$ を使って $\nabla \times \bm{B} = \mu (\nabla \times \bm{H}) = \mu \sigma \bm{E}$ としました。

$$ \nabla \times (\nabla \times \bm{E}) = -\mu \sigma \frac{\partial \bm{E}}{\partial t} $$

左辺にベクトル三重積の公式 $\nabla \times (\nabla \times \bm{F}) = \nabla(\nabla \cdot \bm{F}) – \nabla^2 \bm{F}$ を適用します。

$$ \nabla(\nabla \cdot \bm{E}) – \nabla^2 \bm{E} = -\mu \sigma \frac{\partial \bm{E}}{\partial t} $$

導体内部で自由電荷密度 $\rho$ がゼロと仮定できる場合（電荷の緩和時間 $\varepsilon/\sigma$ は銅で $10^{-19}$ s程度と極めて短く、マクロな時間スケールでは電荷は瞬時に表面に逃げる）、$\nabla \cdot \bm{E} = 0$ となります。したがって

$$ \nabla^2 \bm{E} = \mu \sigma \frac{\partial \bm{E}}{\partial t} $$

これが拡散方程式（diffusion equation）です。波動方程式 $\nabla^2 \bm{E} = \mu \varepsilon \, \partial^2 \bm{E}/\partial t^2$ と比べてみてください。波動方程式は時間の2階微分を含みますが、拡散方程式は1階微分です。この違いが、電磁場の振る舞いを「伝搬」から「拡散」へと根本的に変えています。

同様にして、磁場 $\bm{H}$ についても全く同じ形の拡散方程式が成り立ちます。

$$ \nabla^2 \bm{H} = \mu \sigma \frac{\partial \bm{H}}{\partial t} $$

拡散方程式が得られたことで、この方程式を解けば導体内部の電磁場の分布がわかります。次に、具体的な平面波の解を求めて侵入深さを導出しましょう。

侵入深さ（スキンデプス）の導出

平面波解の仮定

導体表面を $z = 0$ とし、$z > 0$ の方向に導体内部が広がっているとします。表面に $x$ 方向の電場 $E_x$ が加わっている状況を考えます。電場は $z$ と $t$ のみの関数と仮定します（$y$ 方向には一様）。

時間調和場を仮定し、$E_x(z, t) = \hat{E}(z) \, e^{j\omega t}$ と書きます。拡散方程式に代入すると

$$ \frac{d^2 \hat{E}}{dz^2} = j\omega \mu \sigma \, \hat{E} $$

となります。ここで $\partial/\partial t \to j\omega$ の置き換えを行いました。

特性方程式の解

この常微分方程式は定数係数の2階線形ODEです。$\hat{E}(z) = A e^{\gamma z}$ と仮定して代入すると、特性方程式

$$ \gamma^2 = j\omega \mu \sigma $$

が得られます。$\gamma$ を求めるために $j$ の平方根を計算しましょう。$j = e^{j\pi/2}$ ですから

$$ \sqrt{j} = e^{j\pi/4} = \frac{1+j}{\sqrt{2}} $$

この結果を $\gamma^2 = j\omega\mu\sigma$ に適用し、$\sqrt{j} = (1+j)/\sqrt{2}$ を利用すると、$\gamma$ は次のように求まります。

$$ \gamma = \pm \sqrt{j\omega\mu\sigma} = \pm (1+j) \sqrt{\frac{\omega\mu\sigma}{2}} $$

ここで侵入深さ $\delta$ を次のように定義します。

$$ \delta = \sqrt{\frac{2}{\omega\mu\sigma}} $$

この $\delta$ を使うと、$\sqrt{\omega\mu\sigma/2} = 1/\delta$ の関係から $(1+j)/\delta$ という形が現れ、$\gamma$ は簡潔に

$$ \gamma = \pm \frac{1+j}{\delta} $$

と書けます。

物理的に意味のある解の選択

2つの解 $\gamma = +(1+j)/\delta$ と $\gamma = -(1+j)/\delta$ のうち、$z \to +\infty$ で電場が発散する解は物理的に不適切です。$z$ が増加すると減衰する解、すなわち

$$ \hat{E}(z) = E_0 \, e^{-(1+j)z/\delta} $$

を選びます。実部と虚部を分離すると

$$ E_x(z, t) = E_0 \, e^{-z/\delta} \cos\left(\omega t – \frac{z}{\delta}\right) $$

この式から、導体内部の電場について2つの重要な特徴が読み取れます。

1. 指数減衰: 振幅が $e^{-z/\delta}$ に比例して減衰する。$z = \delta$ で振幅は表面値の $1/e \approx 37\%$ に、$z = 3\delta$ では $e^{-3} \approx 5\%$ に、$z = 5\delta$ では $e^{-5} \approx 0.7\%$ にまで落ちる
2. 位相遅れ: 内部に進むほど位相が $z/\delta$ だけ遅れる。表面での位相に対して、$z = \pi\delta$ の深さで逆位相になる

オームの法則 $\bm{J} = \sigma \bm{E}$ より、電流密度も全く同じ分布に従います。

$$ J_x(z, t) = \sigma E_0 \, e^{-z/\delta} \cos\left(\omega t – \frac{z}{\delta}\right) $$

つまり、電流は表面近傍の厚さ $\delta$ 程度の薄い層に集中して流れます。これが表皮効果の数学的な表現です。

侵入深さの公式の物理的意味

改めて侵入深さの公式を眺めてみましょう。

$$ \begin{equation} \delta = \sqrt{\frac{2}{\omega\mu\sigma}} = \frac{1}{\sqrt{\pi f \mu \sigma}} \end{equation} $$

この式は3つのパラメータへの依存性を示しています。

• 角周波数 $\omega$ に反比例の平方根: 周波数が高いほど侵入深さは浅くなる。周波数を4倍にすると $\delta$ は半分に
• 透磁率 $\mu$ に反比例の平方根: 強磁性体（鉄など）は $\mu$ が大きいので、侵入深さが極めて浅い。これが鉄の電磁シールド効果が高い理由の1つ
• 導電率 $\sigma$ に反比例の平方根: 導電率が高いほど渦電流が大きく、内部の磁場変化を打ち消す力が強いため、侵入が浅くなる

直感的には、$\delta$ は「電磁場の拡散が交流の半周期の間にどれだけ浸透できるか」の尺度です。拡散定数 $D = 1/(\mu\sigma)$ を用いると、半周期 $T/2 = \pi/\omega$ の間に拡散が到達する距離は $\sqrt{D \cdot T/2} = \sqrt{\pi/(\omega\mu\sigma)}$ となり、$\delta$ と同じオーダーになります。

侵入深さの物理的意味が明確になったところで、次に具体的な金属と周波数で $\delta$ がどのような値を取るかを見ていきましょう。

各金属・各周波数での侵入深さの具体値

代表的な金属の電気的特性

表皮効果の大きさは材料の導電率 $\sigma$ と透磁率 $\mu$ で決まります。代表的な金属の特性値を以下にまとめます。

周波数ごとの侵入深さ

銅（$\sigma = 5.96 \times 10^7$ S/m, $\mu_r = 1$）を例に、各周波数での侵入深さを計算してみましょう。$\mu = \mu_0 = 4\pi \times 10^{-7}$ H/m を用いると

$$ \delta_{\text{Cu}} = \frac{1}{\sqrt{\pi f \mu_0 \sigma}} = \frac{1}{\sqrt{\pi f \times 4\pi \times 10^{-7} \times 5.96 \times 10^7}} $$

50 Hzでは約9 mmですから、直径数cmの送電線では断面全体にほぼ均等に電流が流れます。しかし1 MHzになると66 $\mu$m — 髪の毛1本程度の厚さにまで薄くなります。GHz帯では $\mu$m オーダーとなり、導体のメッキの厚さ程度しか電流は浸透しません。

鉄の場合 — 強磁性体の劇的な表皮効果

鉄は $\mu_r \approx 200$ と大きな比透磁率を持ちます。侵入深さの式は $\mu$ にも反比例の平方根で依存するため、鉄の侵入深さは銅に比べて $\sqrt{\mu_{r,\text{Fe}} \cdot \sigma_{\text{Cu}} / \sigma_{\text{Fe}}}$ 倍程度浅くなります。

50 Hzでの鉄の侵入深さは約0.36 mmです。銅の9.34 mmと比べて25分の1以下です。このため、鉄芯を用いた変圧器では渦電流損失を低減するために鉄芯を薄い積層板（ラミネーション）で構成します。各板の厚さを侵入深さの数倍以下に抑えることで、渦電流のループを分断し、損失を大幅に低減できるのです。

各金属の侵入深さの具体値を把握したところで、次に表皮効果が導体の抵抗にどのような影響を与えるかを定量化しましょう。そのために「表面抵抗」という概念を導入します。

表面抵抗 $R_s$ の導出

表面抵抗の定義と意味

表面から内部に向かって電流密度が指数関数的に減衰するということは、実効的に電流が流れている領域の厚さが侵入深さ $\delta$ 程度であることを意味します。もし電流が完全に厚さ $\delta$ の薄い層に均一に分布していると仮定すれば、幅 $w$、長さ $l$ の導体表面の抵抗は

$$ R = \frac{l}{\sigma \cdot w \cdot \delta} $$

と見積もれます。この見積もりは実は正確です。以下で厳密に計算してみましょう。

厳密な導出

導体表面（$z = 0$）での電流密度を $J_0 = \sigma E_0$ とします。単位幅あたりの全電流は、表面から無限遠まで電流密度を積分して得られます。

$$ J_{\text{total}} = \int_0^{\infty} J_0 \, e^{-(1+j)z/\delta} \, dz $$

被積分関数は複素指数関数なので、積分を実行すると

$$ J_{\text{total}} = J_0 \left[ \frac{-\delta}{1+j} \, e^{-(1+j)z/\delta} \right]_0^{\infty} $$

$z \to \infty$ で指数関数はゼロに収束し、$z = 0$ での値を代入すると

$$ J_{\text{total}} = J_0 \cdot \frac{\delta}{1+j} $$

ここで $1/(1+j)$ の有理化を行います。分子・分母に共役な $(1-j)$ を掛けると $1/(1+j) = (1-j)/|1+j|^2 = (1-j)/2$ となるので、これを代入すると

$$ J_{\text{total}} = J_0 \cdot \frac{\delta(1-j)}{2} $$

この結果の大きさ（絶対値）は

$$ |J_{\text{total}}| = J_0 \cdot \frac{\delta}{\sqrt{2}} = J_0 \cdot \frac{\delta}{\sqrt{2}} $$

となります。一方、もし電流が厚さ $\delta$ の層に均一に分布していたら $J_{\text{total}} = J_0 \cdot \delta$ になるはずです。実際には指数減衰しているため、$1/\sqrt{2}$ の係数がかかります。

表面インピーダンスの導出

より体系的に扱うために、表面インピーダンス $Z_s$ を定義しましょう。表面での電場 $E_0$ と、単位幅あたりの全電流 $J_{\text{total}}$ の比です。

$$ Z_s = \frac{E_0}{J_{\text{total}}} = \frac{E_0}{J_0 \cdot \delta/(1+j)} = \frac{1}{\sigma} \cdot \frac{1+j}{\delta} $$

$\delta = \sqrt{2/(\omega\mu\sigma)}$ を代入して整理します。まず $1/(\sigma\delta)$ を計算すると

$$ \frac{1}{\sigma\delta} = \frac{1}{\sigma} \sqrt{\frac{\omega\mu\sigma}{2}} = \sqrt{\frac{\omega\mu}{2\sigma}} $$

この結果に $(1+j)$ を掛けると、表面インピーダンスは次の形にまとまります。

$$ Z_s = (1+j) \sqrt{\frac{\omega\mu}{2\sigma}} $$

この表面インピーダンスの実部が表面抵抗 $R_s$ です。

$$ \begin{equation} R_s = \text{Re}(Z_s) = \sqrt{\frac{\omega\mu}{2\sigma}} = \frac{1}{\sigma\delta} \end{equation} $$

最後の等号 $R_s = 1/(\sigma\delta)$ は非常に直感的です。導電率 $\sigma$ の導体が厚さ $\delta$ だけ電流を通しているのですから、単位面積あたりの抵抗が $1/(\sigma\delta)$ になるのは自然な結果です。

表面抵抗の具体値

銅の場合、$R_s$ は

$$ R_s = \sqrt{\frac{\pi f \mu_0}{\sigma}} = \sqrt{\frac{\pi f \times 4\pi \times 10^{-7}}{5.96 \times 10^7}} $$

「$\Omega$/sq」（ohm per square）は面抵抗の単位で、正方形の薄膜の対辺間の抵抗値に相当します。GHz帯でも数十m$\Omega$/sq と小さく見えますが、導波管やマイクロストリップの伝送損失に直結するため、無視できない値です。

表面抵抗の概念を理解したところで、次にこの表面抵抗が高周波回路の設計にどのような実用的影響を及ぼすかを考察しましょう。

高周波回路への影響 — 実効抵抗の増大

直流抵抗から交流抵抗へ

半径 $a$ の円形断面を持つ導線の直流抵抗（単位長さあたり）は

$$ R_{\text{DC}} = \frac{1}{\sigma \pi a^2} $$

です。断面全体を電流が均一に流れるからです。

一方、表皮効果が顕著な高周波では、電流は表面の厚さ $\delta$ の薄い層に集中します。$\delta \ll a$ の場合、電流が流れる断面積は全断面積 $\pi a^2$ ではなく、表面の円環部分の面積 $2\pi a \delta$ 程度になります（円周 $2\pi a$ × 厚さ $\delta$）。したがって、高周波での実効抵抗は

$$ R_{\text{AC}} \approx \frac{1}{\sigma \cdot 2\pi a \delta} = \frac{R_s}{2\pi a} $$

直流抵抗と交流抵抗の比は

$$ \frac{R_{\text{AC}}}{R_{\text{DC}}} = \frac{\pi a^2}{2\pi a \delta} = \frac{a}{2\delta} $$

これは $\delta \ll a$ のとき非常に大きくなります。例えば、半径 $a = 1$ mm の銅線を 1 GHz で使う場合、$\delta \approx 2.09$ $\mu$m ですから

$$ \frac{R_{\text{AC}}}{R_{\text{DC}}} \approx \frac{10^{-3}}{2 \times 2.09 \times 10^{-6}} \approx 239 $$

交流抵抗は直流抵抗の約240倍にもなります。太い導線を使っても内部には電流が流れないため、高周波での抵抗低減には限界があるのです。

高周波回路での対策

この問題に対する工学的な対策として、以下のような手法が用いられます。

リッツ線（Litz wire）は、細い絶縁された素線を多数撚り合わせた特殊な導線です。各素線の直径を侵入深さ程度にすることで、それぞれの素線の断面全体に電流を流し、実効的な抵抗を下げることができます。数十kHzから数MHz程度の周波数で効果的です。

導体表面のメッキは、GHz帯の高周波回路で重要な技術です。電流は表面の数$\mu$m しか流れないため、基板の銅配線表面に導電率の高い銀やニッケルをメッキすることで表面抵抗を低減できます。逆に、酸化や腐食で表面の導電率が下がると、高周波特性が大きく劣化します。

導体の表面粗さも無視できない影響を持ちます。表面の凹凸がスキンデプスと同程度の大きさになると、電流経路が凹凸に沿って蛇行するため実効的な経路長が長くなり、抵抗が増大します。このため、高周波基板では導体表面の粗さを管理することが重要です。

高周波回路における抵抗増大の問題を見てきました。次に、表皮効果を逆に「利点」として活用する電磁シールドの原理を考えましょう。

電磁シールドへの応用

シールド効果の原理

電磁シールドは、外部からの電磁波が機器内部に侵入するのを防ぐ（あるいは内部からの電磁波の漏洩を防ぐ）技術です。金属板による電磁シールドの効果は、大きく分けて3つのメカニズムで構成されます。

1. 反射損失（reflection loss）: 電磁波が空気と金属の境界で反射される
2. 吸収損失（absorption loss）: 金属内部を通過する間に指数関数的に減衰する
3. 内部多重反射損失（multiple reflection loss）: 金属内部で前面と背面の間を多重反射しながら減衰する

表皮効果は主に2番目の吸収損失に直接関係しています。

吸収損失の計算

厚さ $t$ の金属板を透過する電磁波は、表皮効果により

$$ |\bm{E}(t)| = |\bm{E}(0)| \, e^{-t/\delta} $$

に減衰します。吸収損失を dB（デシベル）で表すと

$$ A = 20 \log_{10}\left(\frac{|\bm{E}(0)|}{|\bm{E}(t)|}\right) = 20 \log_{10}(e^{t/\delta}) = \frac{20 t}{\delta \ln 10} \approx 8.686 \frac{t}{\delta} \quad \text{[dB]} $$

侵入深さ1つ分（$t = \delta$）で約 8.7 dB の減衰が生じます。これは振幅にして $1/e \approx 0.37$ 倍、電力にして $1/e^2 \approx 0.14$ 倍（約86%の電力が吸収される）に相当します。

実用的には、$t = 5\delta$ とすると吸収損失は約43 dB（電力で約50,000分の1）に達します。

周波数と材料によるシールド設計

シールド設計では、対象周波数での侵入深さに基づいて金属板の厚さを決めます。いくつかの例を見てみましょう。

電子機器のEMC対策（100 MHz〜1 GHz）: 銅板の侵入深さは 2〜7 $\mu$m です。0.1 mmの銅板でも $t/\delta > 14$ となり、吸収損失だけで120 dB以上が得られます。ただし実際のシールド効果は、板の接合部やケーブル貫通部からの漏洩が支配的です。

低周波磁気シールド（50/60 Hz）: 銅板の侵入深さは約9 mmと深いため、薄い銅板では吸収損失がほとんど得られません。ここで鉄やパーマロイ（高透磁率合金、$\mu_r \sim 80{,}000$）が活躍します。パーマロイは侵入深さが極めて浅いため、薄い板でも低周波磁場に対して優れたシールド効果を発揮します。精密計測器やCRT（ブラウン管）ディスプレイの磁気シールドに利用されてきました。

宇宙機の電磁適合性: 宇宙機では軽量化が必須のため、アルミニウム合金の筐体が多用されます。GHz帯のRF干渉に対しては十分なシールド効果が得られますが、低周波の磁気干渉に対しては追加の対策（磁気シールド材の局所配置など）が必要になることがあります。

反射損失の簡易評価

完全を期すために、反射損失についても触れておきましょう。良導体の表面に電磁波が入射すると、金属の波動インピーダンスと自由空間のインピーダンスの大きな不整合により、大部分が反射されます。

金属の表面インピーダンスは $Z_s = (1+j)/(\sigma\delta)$ で、その大きさは

$$ |Z_s| = \frac{\sqrt{2}}{\sigma\delta} = \sqrt{\frac{2\omega\mu}{\sigma}} $$

自由空間のインピーダンス $\eta_0 = \sqrt{\mu_0/\varepsilon_0} \approx 377 \; \Omega$ に対して、良導体では $|Z_s| \ll \eta_0$ です。この不整合による反射損失は

$$ R_{\text{ref}} \approx 20 \log_{10}\left(\frac{\eta_0}{4|Z_s|}\right) \quad \text{[dB]} $$

と近似でき、高い周波数ほど（$|Z_s|$ が小さいほど）反射損失は大きくなります。電場成分に対しては低周波でも大きな反射損失が得られますが、磁場成分に対しては波源に近い近傍界で反射損失が減少する点に注意が必要です。

電磁シールドの原理を理解したところで、次に導波管という別の応用について考えましょう。導波管の壁面損失は表面抵抗から直接計算されます。

導波管の壁面損失

導波管における損失の起源

導波管（waveguide）は、金属の管の中に電磁波を閉じ込めて伝送する構造です。理想的な完全導体で作られた導波管には損失がありませんが、現実の導体には有限の導電率があるため、壁面を流れる電流がジュール熱を発生させます。

この壁面損失は表面抵抗 $R_s$ を使って体系的に計算できます。壁面上の表面電流密度 $\bm{J}_s$ が作る単位長さあたりの損失電力は

$$ P_{\text{loss}} = \frac{R_s}{2} \oint_C |\bm{J}_s|^2 \, dl $$

です。ここで積分は導波管の断面の周に沿って行います。係数 $1/2$ は時間平均に由来します。

矩形導波管の減衰定数

幅 $a$、高さ $b$ の矩形導波管で $\text{TE}_{10}$ モードを伝送する場合を考えましょう。伝送電力 $P_{\text{trans}}$ と壁面損失 $P_{\text{loss}}$ の比から、減衰定数 $\alpha$ が得られます。

$$ \alpha = \frac{P_{\text{loss}}}{2 P_{\text{trans}}} $$

$\text{TE}_{10}$ モードの場合、減衰定数は

$$ \alpha = \frac{R_s}{b \eta \sqrt{1 – (f_c/f)^2}} \left(1 + \frac{2b}{a} \left(\frac{f_c}{f}\right)^2 \right) $$

ここで $\eta = \sqrt{\mu/\varepsilon}$ は媒質の波動インピーダンス、$f_c = c/(2a)$ は遮断周波数です。

この式から、導波管の壁面損失を低減するためのいくつかの指針が読み取れます。

• $R_s$ を小さくする: 導電率の高い材料（銅、銀）を使う。壁面に銀メッキを施すのは高性能導波管の標準的な手法
• 導波管のサイズを大きくする: $a$, $b$ が大きいほど壁面電流密度が分散し、損失が減る
• 遮断周波数に近い周波数を避ける: $f \to f_c$ で $\alpha$ が発散する（群速度がゼロに近づき、エネルギーが壁面の往復に多くの時間を費やすため）

実際の衛星通信や地上マイクロ波中継では、導波管の壁面損失が受信感度や雑音温度に直接影響するため、表面抵抗の正確な評価が不可欠です。

ここまでで表皮効果の理論と応用を一通り見てきました。次に、Pythonを使って侵入深さの周波数依存性と導体内部の電流分布を可視化し、理論の理解を深めましょう。

Pythonでの実装と可視化

侵入深さの周波数依存性

まず、各金属の侵入深さが周波数とともにどう変化するかをプロットしてみましょう。対数スケールで見ることで、広い周波数範囲にわたる振る舞いが一目でわかります。

import numpy as np
import matplotlib.pyplot as plt

# 物理定数
mu_0 = 4 * np.pi * 1e-7  # 真空の透磁率 [H/m]

# 金属の特性 (導電率 [S/m], 比透磁率)
metals = {
    "Cu (銅)":   (5.96e7, 1.0),
    "Al (アルミ)": (3.50e7, 1.0),
    "Ag (銀)":   (6.17e7, 1.0),
    "Au (金)":   (4.10e7, 1.0),
    "Fe (鉄)":   (1.04e7, 200.0),
}

# 周波数範囲
freq = np.logspace(1, 11, 500)  # 10 Hz ~ 100 GHz

# 侵入深さの計算
fig, ax = plt.subplots(figsize=(10, 6))

for name, (sigma, mu_r) in metals.items():
    mu = mu_r * mu_0
    delta = 1.0 / np.sqrt(np.pi * freq * mu * sigma)
    ax.loglog(freq, delta * 1e3, label=name, linewidth=2)

# 参考線
ax.axhline(y=0.035, color='gray', linestyle=':', alpha=0.7)
ax.text(1e2, 0.04, "PCB銅箔厚 (35 μm)", fontsize=9, color='gray')
ax.axhline(y=1e-3, color='gray', linestyle=':', alpha=0.7)
ax.text(1e2, 1.2e-3, "1 μm", fontsize=9, color='gray')

ax.set_xlabel("Frequency [Hz]", fontsize=12)
ax.set_ylabel("Skin depth [mm]", fontsize=12)
ax.set_title("Skin Depth vs Frequency for Various Metals", fontsize=14)
ax.legend(fontsize=11)
ax.grid(True, which="both", alpha=0.3)
ax.set_xlim(1e1, 1e11)
ax.set_ylim(1e-4, 1e2)
plt.tight_layout()
plt.savefig("skin_depth_vs_freq.png", dpi=150, bbox_inches="tight")
plt.show()

このグラフから、いくつかの重要な特徴が読み取れます。

1. 全ての金属で侵入深さは周波数に対して傾き $-1/2$ の直線になる: 対数グラフでの傾きが $-1/2$ であることは、$\delta \propto f^{-1/2}$ という理論の予測と完全に一致しています。
2. 鉄は他の金属に比べて1桁以上浅い: これは比透磁率 $\mu_r = 200$ の効果です。低周波の電磁シールドに鉄やパーマロイが使われる理由がグラフから直感的にわかります。
3. 銀と銅はほぼ重なる: 導電率が近い（$6.17 \times 10^7$ vs $5.96 \times 10^7$）ため、侵入深さもほぼ同じです。高周波回路で銅の代わりに銀メッキを施しても劇的な改善は期待できず、主に酸化防止と表面粗さの改善が銀メッキの利点であることがわかります。

導体内部の電流分布の可視化

次に、導体内部での電流密度分布を可視化します。異なる周波数で電流がどの程度表面に集中するかを見てみましょう。

import numpy as np
import matplotlib.pyplot as plt

# 銅の特性
sigma = 5.96e7  # [S/m]
mu = 4 * np.pi * 1e-7  # [H/m]

# 異なる周波数でのスキンデプス
frequencies = [60, 1e3, 1e5, 1e6, 1e8]
labels = ["60 Hz", "1 kHz", "100 kHz", "1 MHz", "100 MHz"]
colors = plt.cm.viridis(np.linspace(0.1, 0.9, len(frequencies)))

fig, axes = plt.subplots(1, 2, figsize=(14, 6))

# --- 左図: 正規化された電流密度分布 ---
ax = axes[0]
z_mm = np.linspace(0, 25, 1000)  # 表面からの距離 [mm]

for f, label, color in zip(frequencies, labels, colors):
    delta = 1.0 / np.sqrt(np.pi * f * mu * sigma)
    delta_mm = delta * 1e3
    J_ratio = np.exp(-z_mm / delta_mm)
    ax.plot(z_mm, J_ratio, label=f"{label} (δ={delta_mm:.3f} mm)",
            color=color, linewidth=2)

ax.set_xlabel("Depth from surface [mm]", fontsize=12)
ax.set_ylabel(r"$|J(z)| / |J(0)|$", fontsize=12)
ax.set_title("Current Density Distribution in Cu", fontsize=13)
ax.legend(fontsize=9, loc="upper right")
ax.grid(True, alpha=0.3)
ax.set_xlim(0, 25)
ax.set_ylim(0, 1.05)

# --- 右図: 円形断面の電流密度（疑似カラー） ---
ax2 = axes[1]
a = 5.0  # 導線の半径 [mm]
N = 500
x = np.linspace(-a, a, N)
y = np.linspace(-a, a, N)
X, Y = np.meshgrid(x, y)
R = np.sqrt(X**2 + Y**2)

# 100 kHz での電流分布（外側からの距離で計算）
f_plot = 1e5
delta_plot = 1.0 / np.sqrt(np.pi * f_plot * mu * sigma) * 1e3  # [mm]
depth = a - R  # 表面からの距離
J_density = np.where(R <= a, np.exp(-depth / delta_plot), np.nan)

im = ax2.pcolormesh(X, Y, J_density, cmap="inferno", shading="auto")
circle = plt.Circle((0, 0), a, fill=False, color="white",
                     linewidth=1.5, linestyle="--")
ax2.add_patch(circle)
ax2.set_xlabel("x [mm]", fontsize=12)
ax2.set_ylabel("y [mm]", fontsize=12)
ax2.set_title(f"Current Density in Wire Cross-section\n"
              f"(f = 100 kHz, δ = {delta_plot:.3f} mm)", fontsize=12)
ax2.set_aspect("equal")
fig.colorbar(im, ax=ax2, label=r"$|J| / |J_{\mathrm{surface}}|$", shrink=0.8)

plt.tight_layout()
plt.savefig("current_distribution.png", dpi=150, bbox_inches="tight")
plt.show()

左のグラフからは、周波数ごとの電流の浸透の違いが明瞭にわかります。60 Hzでは数mm以上にわたって電流が浸透していますが、100 kHzになると表面から0.5 mm程度で急激に減衰し、100 MHzではグラフ上ではほぼ表面の一点にしか見えません。右の断面図では、100 kHzにおいて半径5 mmの銅線の内部にはほとんど電流が流れておらず、表面近傍の薄い環状領域にのみ電流が集中していることが視覚的に確認できます。

表面抵抗と交流抵抗比の計算

最後に、導線の交流抵抗が直流抵抗の何倍になるかを周波数の関数としてプロットします。

import numpy as np
import matplotlib.pyplot as plt

# 銅の特性
sigma = 5.96e7  # [S/m]
mu = 4 * np.pi * 1e-7  # [H/m]

# 導線半径のバリエーション
radii_mm = [0.1, 0.5, 1.0, 2.5]  # [mm]
freq = np.logspace(1, 9, 500)  # 10 Hz ~ 1 GHz

fig, ax = plt.subplots(figsize=(10, 6))

for a_mm in radii_mm:
    a = a_mm * 1e-3  # [m]
    delta = 1.0 / np.sqrt(np.pi * freq * mu * sigma)
    # 厳密解ではベッセル関数が必要だが、δ << a では以下の近似が良い
    # δ >> a ではR_AC/R_DC ≈ 1
    ratio = np.where(delta < a, a / (2 * delta), 1.0)
    ax.loglog(freq, ratio, label=f"a = {a_mm} mm", linewidth=2)

ax.axhline(y=1, color='gray', linestyle='--', alpha=0.5)
ax.set_xlabel("Frequency [Hz]", fontsize=12)
ax.set_ylabel(r"$R_{\mathrm{AC}} / R_{\mathrm{DC}}$", fontsize=12)
ax.set_title("AC/DC Resistance Ratio for Cu Wire", fontsize=14)
ax.legend(fontsize=11)
ax.grid(True, which="both", alpha=0.3)
ax.set_xlim(1e1, 1e9)
ax.set_ylim(0.8, 1e4)
plt.tight_layout()
plt.savefig("ac_dc_resistance_ratio.png", dpi=150, bbox_inches="tight")
plt.show()

このグラフから、2つの重要な工学的示唆が得られます。

1. 太い導線ほど低い周波数から交流抵抗比が増大し始める: 半径2.5 mmの導線では数百Hz程度から表皮効果が顕著になりますが、半径0.1 mmの細い導線では数十kHz付近まで直流抵抗とほぼ等しい値を保ちます。これはリッツ線が細い素線を束ねる理由を裏付けています。
2. 高周波領域では全ての導線で $R_{\text{AC}}/R_{\text{DC}} \propto \sqrt{f}$ となる: 対数グラフでの傾きが $+1/2$ に漸近することが確認できます。これは $R_{\text{AC}} \propto R_s \propto \sqrt{f}$ であり、$R_{\text{DC}}$ は周波数に依存しないためです。

電磁シールド効果のシミュレーション

金属板の厚さと周波数から、吸収損失がどの程度得られるかを可視化します。

import numpy as np
import matplotlib.pyplot as plt

# 物理定数
mu_0 = 4 * np.pi * 1e-7

# 金属: 銅とアルミニウム
metals_shield = {
    "Cu": (5.96e7, 1.0),
    "Al": (3.50e7, 1.0),
}

freq = np.logspace(3, 10, 500)  # 1 kHz ~ 10 GHz
thicknesses_mm = [0.01, 0.05, 0.1, 0.5, 1.0]  # [mm]

fig, axes = plt.subplots(1, 2, figsize=(14, 6))

for idx, (name, (sigma, mu_r)) in enumerate(metals_shield.items()):
    ax = axes[idx]
    mu = mu_r * mu_0

    for t_mm in thicknesses_mm:
        t = t_mm * 1e-3
        delta = 1.0 / np.sqrt(np.pi * freq * mu * sigma)
        absorption_dB = 8.686 * t / delta
        ax.semilogx(freq, absorption_dB,
                     label=f"t = {t_mm} mm", linewidth=2)

    ax.axhline(y=40, color='red', linestyle=':', alpha=0.5)
    ax.text(2e3, 42, "40 dB (一般的な要求)", fontsize=9, color='red')
    ax.set_xlabel("Frequency [Hz]", fontsize=12)
    ax.set_ylabel("Absorption Loss [dB]", fontsize=12)
    ax.set_title(f"Shielding Absorption Loss ({name})", fontsize=13)
    ax.legend(fontsize=10)
    ax.grid(True, which="both", alpha=0.3)
    ax.set_xlim(1e3, 1e10)
    ax.set_ylim(0, 200)

plt.tight_layout()
plt.savefig("shielding_absorption.png", dpi=150, bbox_inches="tight")
plt.show()

このグラフからは、シールド設計における重要な知見が読み取れます。

1. 高周波ほど薄い板で十分なシールド効果が得られる: 100 MHzでは0.01 mm（10 $\mu$m）の銅箔でも40 dB以上の吸収損失があります。アルミ箔（家庭用は約15 $\mu$m）でもGHz帯では数十dBの遮蔽が得られることがわかります。
2. 低周波のシールドは困難: 1 kHz付近では0.1 mmの銅板でも吸収損失は数dBにすぎません。低周波の磁気シールドには厚い鉄板や高透磁率材料が必要であることが定量的に裏付けられます。
3. 銅とアルミの差は比較的小さい: 導電率の比が約1.7倍なので、侵入深さの比は $\sqrt{1.7} \approx 1.3$ 倍程度です。重量あたりのシールド効果ではアルミの方が有利な場合も多く、航空宇宙分野でアルミ合金が好まれる理由の一端が見えます。

各金属の侵入深さの数値表出力

最後に、一覧表として具体的な数値を確認するコードを示します。

import numpy as np

mu_0 = 4 * np.pi * 1e-7

metals_table = {
    "Cu": (5.96e7, 1.0),
    "Al": (3.50e7, 1.0),
    "Ag": (6.17e7, 1.0),
    "Au": (4.10e7, 1.0),
    "Fe": (1.04e7, 200.0),
}

freqs = [50, 1e3, 1e5, 1e6, 1e8, 1e9, 1e10]
freq_labels = ["50 Hz", "1 kHz", "100 kHz", "1 MHz", "100 MHz", "1 GHz", "10 GHz"]

print(f"{'Metal':<6}", end="")
for fl in freq_labels:
    print(f"{fl:>12}", end="")
print()
print("-" * 90)

for name, (sigma, mu_r) in metals_table.items():
    mu = mu_r * mu_0
    print(f"{name:<6}", end="")
    for f in freqs:
        delta = 1.0 / np.sqrt(np.pi * f * mu * sigma)
        if delta >= 1e-3:
            print(f"{delta*1e3:>10.2f} mm", end="")
        else:
            print(f"{delta*1e6:>10.2f} um", end="")
    print()

出力結果から、銅と銀の侵入深さがほぼ等しいこと、鉄が全周波数帯で桁違いに浅い侵入深さを持つこと、そしてどの金属でも周波数が10倍になると侵入深さが $\sqrt{10} \approx 3.16$ 分の1になっていることが数値で確認できます。これらはすべて、$\delta = 1/\sqrt{\pi f \mu \sigma}$ という理論式と完全に整合しています。

補足: 拡散方程式と波動方程式の対比

本記事の締めくくりとして、導体内部の電磁場が従う拡散方程式と、真空中（あるいは誘電体中）の波動方程式を対比しておきましょう。両者の違いを明確に理解しておくことは、表皮効果の本質を把握する上で極めて重要です。

真空中の波動方程式は

$$ \nabla^2 \bm{E} = \mu_0 \varepsilon_0 \frac{\partial^2 \bm{E}}{\partial t^2} $$

であり、右辺は時間の2階微分です。この方程式の解は $e^{j(\omega t – kz)}$ の形の伝搬波であり、振幅は距離によらず一定です。エネルギーが空間を伝搬していきます。

一方、導体内部の拡散方程式は

$$ \nabla^2 \bm{E} = \mu \sigma \frac{\partial \bm{E}}{\partial t} $$

であり、右辺は時間の1階微分です。この方程式の解は $e^{-(1+j)z/\delta} e^{j\omega t}$ の形であり、振幅が $e^{-z/\delta}$ で指数減衰します。エネルギーは導体内部に「浸透」し、ジュール熱として消散します。

この違いは、拡散定数 $D = 1/(\mu\sigma)$ と波動の位相速度 $v = 1/\sqrt{\mu\varepsilon}$ の物理的意味の違いに対応しています。拡散では信号がゆっくりと広がりながら減衰するのに対し、波動では信号が一定速度で伝搬します。

良導体の条件 $\sigma \gg \omega\varepsilon$ は、まさに「導体内部では伝導電流が変位電流を圧倒する」ことを意味しており、この条件下で電磁場の振る舞いが波動から拡散に切り替わるのです。逆に、導電率がゼロの完全な誘電体では拡散項が消え、純粋な波動方程式に戻ります。

この対比を理解しておくと、不完全な導体（半導体や海水など、$\sigma$ と $\omega\varepsilon$ が同程度の媒質）での電磁波の振る舞い — 伝搬しながら減衰するという中間的な性質 — も自然に理解できるようになります。

まとめ

本記事では、表皮効果の物理的メカニズムからマクスウェル方程式に基づく厳密な導出、そして工学的応用までを一貫して解説しました。

• 表皮効果の本質: 導体に交流電流を流すと、誘導電場が内部の電流を打ち消す方向に働き、電流が表面近傍に集中する。導体内部の電磁場は拡散方程式に従い、波動ではなく拡散的に振る舞う
• 侵入深さの公式: $\delta = \sqrt{2/(\omega\mu\sigma)} = 1/\sqrt{\pi f \mu \sigma}$ であり、周波数・透磁率・導電率の平方根に反比例する。銅で50 Hzなら約9 mm、1 GHzなら約2 $\mu$m
• 表面抵抗: $R_s = 1/(\sigma\delta) = \sqrt{\omega\mu/(2\sigma)}$ であり、高周波ほど導体の実効抵抗が増大する
• 高周波回路への影響: 交流抵抗は直流抵抗の $a/(2\delta)$ 倍に増大する。リッツ線、表面メッキ、表面粗さの管理が対策として重要
• 電磁シールド: 金属板の吸収損失は $8.686 \, t/\delta$ [dB] で与えられ、高周波ほど薄い板でシールドできる。低周波の磁気シールドには高透磁率材料が必要
• 導波管の壁面損失: 表面抵抗 $R_s$ から壁面損失が直接計算でき、導波管の伝送効率を支配する

表皮効果は電磁気学の基礎理論と高周波工学を橋渡しする重要なトピックです。マクスウェル方程式という基本原理から出発して、実際の回路設計やシールド設計に直結する定量的な結論が得られるという点で、物理学と工学の結びつきを体感できる題材でもあります。

次のステップとして、以下の記事も参考にしてください。

• マイクロ波回路の基礎 — 表面抵抗が回路設計にどう影響するかを詳しく解説
• 伝送線路の理論 — 分布定数回路としての導線の扱い方