静電場の境界要素法
静電場の境界要素法の理論基礎
静電場BEM
先生、静電場にもBEMが使えるんですか？
開放空間の問題にはBEMが最適。FEMは周囲の空間をメッシュ化する必要があるが、BEMは導体/誘電体の表面だけでよい。
グリーン関数 $G = 1/(4\pi r)$（3D）。表面の電位と法線電界を未知数として解く。
表面メッシュだけで無限領域を扱える！
そう。送電線の周囲電界、雷サージのシールド効果、EMC問題の開放空間解析に威力を発揮する。
FEMとBEMの比較
まとめ
境界積分方程式の原点——19世紀の数学者グリーンの孤独な発見
静電場BEMの理論的基盤となる「グリーンの定理」を発表したジョージ・グリーン（1793〜1841）は、ノッティンガムのパン屋の息子で、独学で数学を学んだ人物です。彼が1828年に自費出版した論文は生前ほぼ無視されましたが、ケルビン卿が死後に再発見し、電磁気学の礎となりました。200年近く経ってもそのグリーン関数がBEMの核心に生き続けているのは、驚くべき話です。
静電場の境界要素法の数値計算手法
BEMの離散化
表面を三角形/四角形要素で離散化。各要素の電位$\phi$と法線電束密度$D_n$を未知数として連立方程式を構築:
影響行列$[H], [G]$はグリーン関数の面積分。密行列なので$O(N^2)$のメモリ。
FMM加速
大規模問題ではFMM（高速多重極法）で$O(N\log N)$に高速化。FastCapやCOMSOLのBEMモジュールに実装。
まとめ
BEMの「グリーン関数」は一筋縄ではいかない
境界要素法の核心はグリーン関数——「ある点に置いた点電荷が周囲の電位に与える影響」を数式で表したものです。静電場の場合は $1/r$（距離の逆数）という比較的シンプルな形ですが、異方性誘電体や多層媒質になると途端に複雑になります。半導体パッケージの多層基板では、各誘電体層ごとに異なるグリーン関数が必要で、その積分を正確に行うための数値技術が研究の最前線でもあります。
静電場の境界要素法の実務適用
実務
送電線の電界解析、EMCシールド効果、雷インパルス解析が主な適用先。
チェックリスト
高圧送電の碍子（がいし）設計でBEMが選ばれる理由
送電鉄塔の碍子（絶縁体）には、汚損（ほこりや海塩の付着）による「フラッシュオーバー」防止のため、精密な電位分布解析が不可欠です。碍子の形状は複雑な曲面で、周囲は無限に広がる大気——まさにBEMが得意とする状況です。実際、275kVや500kVの超高圧碍子の設計では、表面の電界強度ピークが20kV/mm以下に収まるよう、BEMで形状最適化が行われています。
静電場の境界要素法のソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| COMSOL (BEMモジュール) | FEM-BEM連成対応 |
| FastCap（MIT） | BEM静電容量。オープンソース |
| Ansys Q3D | 内部的にBEMを使用 |
| Integrated Engineering Software | BEM専用。Coulomb, Amperes |
「BEMに特化したツール」が意外と少ない理由
市場を見渡すと、静電場BEM専用の商用ツールは意外に少なく、多くの場合FEMに「BEM機能を追加」した形態です。理由は実装の難しさ——フルBEMは密行列が生成されるため、大規模問題ではメモリと計算コストが爆発的に増えます。2000年代にFMM（Fast Multipole Method、高速多重極法）が実用化されて初めて大規模BEMが現実的になりました。ToolベンダーがFEM寄りな理由の一つでもあります。
静電場の境界要素法の先端研究
先端
BEMが宇宙機設計の「静電気対策」で輝く理由
人工衛星や探査機は、宇宙空間でプラズマ環境に晒されるため、機体の帯電（スペースクラフト・チャージング）が深刻な問題になります。太陽光や荷電粒子が機体表面に電荷を蓄積し、数千ボルトに達すると部品が破壊されることも。この解析には広大な「無限空間」を扱う境界要素法（BEM）が最適で、FEMのように宇宙空間全体をメッシュで埋める必要がありません。NASAやJAXAも静電場BEMを先端ツールとして活用しています。
静電場の境界要素法のトラブル対応
トラブル
BEMで「行列が特異」になる——近傍要素の落とし穴
静電場BEMのトラブルシューティングで頻出なのが「行列が特異（Singular）になる」エラーです。原因の多くは、評価点（コロケーション点）が積分要素に重なる「自己積分」の扱いミス。理論上は特異積分として解析的に処理すべき箇所を、単純な数値積分でやると行列が壊れます。もう一つよくあるのは「閉じた導体面に内部点を評価しようとする」ケース——BEMは境界のみを離散化するため、内部点の評価には別途処理が必要です。
なった
詳しく
報告