ポアソン方程式（静電場）
理論と物理
ポアソン方程式
先生、静電場のポアソン方程式って？
ガウスの法則＋電界と電位の関係を組み合わせた結果:
一般的な誘電体では:
静電場FEMの支配方程式はこれ。電荷分布$\rho_v$が与えられたとき、電位$\phi$を求める。
構造の平衡方程式$\nabla \cdot \sigma + f = 0$に似てますね。
数学的に同じ楕円型偏微分方程式。構造のヤング率→誘電率$\varepsilon$、外力→電荷密度$\rho_v$、変位→電位$\phi$と読み替えるだけ。
境界条件
種類 数学形式 物理的意味 Dirichlet $\phi = \phi_0$ 電極（電位固定） Neumann $\partial\phi/\partial n = 0$ 対称面、絶縁面 混合 $\varepsilon \partial\phi/\partial n = \sigma_s$ 表面電荷密度
まとめ
| 種類 | 数学形式 | 物理的意味 |
|---|---|---|
| Dirichlet | $\phi = \phi_0$ | 電極（電位固定） |
| Neumann | $\partial\phi/\partial n = 0$ | 対称面、絶縁面 |
| 混合 | $\varepsilon \partial\phi/\partial n = \sigma_s$ | 表面電荷密度 |
リチウムイオン電池の設計でポアソン方程式が活躍する理由
リチウムイオン電池の充放電シミュレーションでは、電解質中のリチウムイオン濃度分布と電位分布を同時に解く必要があります。その電位分布を支配するのがまさにポアソン方程式 $\nabla^2 \phi = -\rho/\varepsilon$ です。電池内部では正極・負極・電解質で誘電率が大きく異なり、境界条件の設定が精度を左右します。実際にパナソニックなどの電池メーカーがCAEを使って「電極の多孔質構造の最適化」をする際、ポアソン方程式を解いた電位分布から局所電流密度を推定し、劣化の起きやすい場所を特定しています。地味な二階微分方程式が、最新の電池技術を陰で支えているんです。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
FEM離散化
Galerkin法で弱形式化し、要素ごとに離散化:
$[B]$は形状関数の勾配行列。構造解析の$[B]^T [D] [B]$から$[D]$を$\varepsilon [I]$に置き換えた形。
2次元の三角形要素なら$B$は定数（1次要素）だから、手計算できますね。
その通り。CST三角形要素の静電場版は教科書の演習問題に最適。要素剛性マトリクスは$\varepsilon \cdot A_e / (4A_e) \cdot [b_i b_j + c_i c_j]$。
まとめ
FDMでポアソン方程式を解く——5点スター差分の直感
ポアソン方程式を有限差分法（FDM）で離散化すると、2Dでは「注目格子点の電位＝周囲4点の電位の平均＋電荷密度項」という形になります。これが「5点スター」差分式で、直感的には「電荷がなければ、その点の電位は隣の平均になる」という意味です。この性質はSOR（逐次過緩和法）という反復解法の速度に直結していて、過緩和係数 $\omega$ の選び方一つで収束が2倍速くなったり逆に発散したりします。FEMが主流の今でも、教育用や概念確認用にFDMでポアソン方程式を実装するエンジニアは多く、行列の「疎性」を活かした数値技法の入り口として最高の題材です。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務
高電圧ブッシング、ケーブルジョイント、半導体のドーピング分布が典型的なポアソン方程式の問題。
チェックリスト
EV充電システムの絶縁設計とポアソン方程式
急速充電器（CHAdeMO、CCS）で使われる高電圧コネクタの絶縁設計にはポアソン方程式を解いた電位分布解析が欠かせません。充電中は数百アンペアの電流が流れるので、コネクタ内部の電位勾配が急峻になる箇所があり、そこが絶縁材の最弱点になります。設計者はポアソン方程式の解から「電界が集中する角部や段差」を特定し、面取り加工や絶縁材の厚みで対策します。実際の設計フローでは、まずポアソン方程式で電位分布を出し、次にその結果を熱解析の入力に使う連成解析になることも多い。1本の方程式が複数の設計フェーズにわたって活躍するんですね。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
全静電場ソルバーがポアソン方程式を解く。特別なツール選定は不要。COMSOL、Maxwell、FEMMのいずれでも。
ポアソン方程式ソルバーの「マルチグリッド法」が速い本当の理由
商用ツールのポアソン方程式ソルバーには「マルチグリッド法」が採用されているものが多くあります。なぜ速いかというと、通常の反復解法はメッシュの細かさに応じた「波長の短い誤差」は素早く消えるけど、「波長の長い誤差」は消えるのがとても遅いからです。マルチグリッド法は粗いメッシュと細かいメッシュを交互に使うことで、どんな波長の誤差も効率よく消せる。計算量が問題サイズNに対して $O(N)$ になるという驚異的な性能が実現できます。COMSOLやANSYSで「大規模静電解析が意外と速く終わった」というのは、こういうソルバーの工夫のおかげです。
選定で最も重要な3つの問い
- 「何を解くか」：ポアソン方程式（静電場）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端
ニューラルネットワークでポアソン方程式を解く——Physics-Informed NNの衝撃
2019年頃から「Physics-Informed Neural Network（PINN）」という手法が注目されています。ニューラルネットワークに支配方程式（ポアソン方程式を含む）の残差も損失関数に組み込んで学習させることで、境界値問題を解くというアプローチです。メッシュを作らなくていい、任意の点で予測できる、というのが大きなメリット。一方で収束の遅さや精度保証の難しさが課題で「FEMを完全に置き換えるのはまだ先」というのが現在の実情です。でも将来的に電池や半導体のデバイスシミュレーションで「大量の設計パラメータを高速にスキャン」する用途では、PINNが主流になる日が来るかもしれません。
トラブルシューティング
トラブル
「ソルバーが収束しない」——ポアソン方程式でよくある罠
ポアソン方程式のFEM解析で「反復解法が収束しない」というトラブルはよく起きます。原因の鉄板は「境界条件の設定ミス」。ディリクレ境界（電位固定）が一箇所もないと連立方程式の係数行列が特異になって解が定まりません。電位の絶対値が決まらず、電界だけが意味を持つ場合は「参照電位点を一点設けてアース」するのがお作法です。もう一つよくある罠が誘電率の単位ミス——相対誘電率 $\varepsilon_r$ を入れるべき場所に真空誘電率 $\varepsilon_0 = 8.85\times10^{-12}$ F/mをそのまま入れてしまって電位が12桁ずれる、という定番のミスがあります。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ポアソン方程式（静電場）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告