電荷分布解析
電荷分布の理論基礎
電荷分布
先生、導体上の電荷分布はどう決まるんですか？
導体は等電位。電荷は表面にのみ分布し、内部の電界はゼロ。表面電荷密度$\sigma_s$は:
$E_n$: 表面の法線方向電界。曲率が大きい部分（尖端）ほど電荷密度と電界が大きくなる。
避雷針の先端に電荷が集中するのはそういう理由ですね。
まさにそう。曲率半径$r$の球導体の表面電界は$E = V/r$で、半径が小さいほど電界が強い。
空間電荷
絶縁体内部や気体中にも電荷が存在する場合がある（空間電荷）:
- 半導体のドーピング → ドナー/アクセプターによる固定電荷
- 気体放電 → イオン・電子の移動電荷
- 絶縁体のトラップ電荷 → 劣化の原因
まとめ
自由電荷と束縛電荷——ガウスの法則が「内部を見る」技術
静電解析の基礎方程式であるガウスの法則∇·D=ρfは、「閉曲面を通る電束D（電位移）の総量が内部の自由電荷ρfに等しい」というシンプルな関係式だ。誘電体内部の分極電荷（束縛電荷）はDの発散には現れず、電束Dで記述することで誘電体の影響を材料定数（比誘電率ε_r）に吸収できる。この定式化のおかげで、誘電体でも金属でも同じFEMフレームワークで静電場を解けるという計算の統一性が生まれている。
電荷分布の数値計算手法
FEMでの電荷計算
FEMで電位$\phi$を解いた後、導体表面の電荷密度は後処理で計算:
全電荷はガウスの法則で面積分。BEM（境界要素法）は表面電荷密度を直接未知数とするので、電荷分布の計算に特に適している。
まとめ
尖端効果のFEM——コロナ放電開始電圧の予測
高圧機器の電極エッジ（尖端）では電界集中が著しく、局所的にコロナ放電が発生する。このコロナ放電開始電圧をFEMで予測するには、尖端部に超細かいメッシュを配置して最大電界強度（Emax）を高精度に求める必要がある。BEM（境界要素法）は無限遠の境界条件を自動的に満たすため、送電線の電極設計では円筒FEMより高精度な電界計算がBEMで得られるケースがある。「電界集中係数Kt=Emax/E_mean」の正確な計算がコロナ対策設計の出発点だ。
電荷分布の実務適用
実務
静電気対策（ESD）、帯電防止設計、半導体のドーピング分布が主な適用。
チェックリスト
静電粉体塗装の最適化——電荷分布シミュレーションで塗りムラをなくす
自動車ボディの静電粉体塗装では、帯電した粉体粒子が電界に沿って部品表面に引き付けられる。複雑な3D形状の自動車ボディでは、凹部（内板等）への塗り込みが難しく「ファラデーケージ効果」で電界が弱い領域に粉体が届かない問題がある。静電場FEMで電界分布を計算し、ガン位置・電圧・エアー流量を最適化することで内板への塗膜均一性が20〜30%改善した事例が国内の自動車塗装工場から報告されている。
電荷分布のソフトウェア比較
ツール
FEM系（COMSOL, Maxwell）で電荷密度の後処理が可能。BEM系（FastCap）は電荷分布が直接得られる。半導体のドーピング分布はSentaurus/SilvacoのTCADツール。
電荷分布解析ツール——体積電荷に対応した静電FEM商用ソフト
電荷分布解析には体積電荷密度（ρ）の空間分布を正確に扱えるFEMツールが必要で、COMSOL MultiphysicsのElectrostatics Module、Ansys Maxwell 3Dが代表例だ。COMSOLはPDE定義の柔軟性が高く、空間電荷限界（SCLC）などの非線形電荷輸送問題も対応できる。Ansys Maxwellは高圧機器（変圧器・ケーブル終端）の絶縁設計に強みがあり、温度依存の誘電率・導電率を考慮した連成解析が得意だ。
電荷分布の先端研究
先端
宇宙機の帯電——静止軌道での電位差が引き起こす静電放電
静止軌道（高度36,000km）では太陽風プラズマの影響で宇宙機が−数kVに帯電し、機体各部の電位差から静電放電（ESD）が発生して搭載電子機器が故障するリスクがある。JAXAはETS-VI（技術試験衛星VI型）での帯電トラブル経験を踏まえ、衛星の電荷蓄積シミュレーション（SPIS等）を開発・標準化した。静電場FEM解析をプラズマ粒子軌道計算と連成させるマルチスケールアプローチが宇宙機帯電対策の最先端だ。
電荷分布のトラブル対応
トラブル
「絶縁体に電荷が残る」——界面電荷のFEM解析で誤差が出やすい理由
絶縁体の界面に電荷が蓄積する現象（界面分極）は、材料の導電率比が10⁶以上になると顕在化する。FEMでこの界面電荷を精度良く計算するには、界面近傍のメッシュを材料の導電率比の平方根以下のサイズに制御する必要があり、設定が不適切だと電位分布が実測と大きくずれる。「高圧ケーブルの絶縁体内部電界分布がシミュレーションと合わない」というトラブルの多くは、この界面電荷モデルの不備が原因だ。
なった
詳しく
報告