ポアソン方程式（静電場）
ポアソン方程式（静電場）の理論基礎
ポアソン方程式
先生、静電場のポアソン方程式って？
ガウスの法則＋電界と電位の関係を組み合わせた結果:
一般的な誘電体では:
静電場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$ | 表面電荷密度 |
まとめ
リチウムイオン電池の設計でポアソン方程式が活躍する理由
リチウムイオン電池の充放電シミュレーションでは、電解質中のリチウムイオン濃度分布と電位分布を同時に解く必要があります。その電位分布を支配するのがまさにポアソン方程式 $\nabla^2 \phi = -\rho/\varepsilon$ です。電池内部では正極・負極・電解質で誘電率が大きく異なり、境界条件の設定が精度を左右します。実際にパナソニックなどの電池メーカーがCAEを使って「電極の多孔質構造の最適化」をする際、ポアソン方程式を解いた電位分布から局所電流密度を推定し、劣化の起きやすい場所を特定しています。地味な二階微分方程式が、最新の電池技術を陰で支えているんです。
ポアソン方程式（静電場）の数値計算手法
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でポアソン方程式を実装するエンジニアは多く、行列の「疎性」を活かした数値技法の入り口として最高の題材です。
ポアソン方程式（静電場）の実務適用
実務
高電圧ブッシング、ケーブルジョイント、半導体のドーピング分布が典型的なポアソン方程式の問題。
チェックリスト
EV充電システムの絶縁設計とポアソン方程式
急速充電器（CHAdeMO、CCS）で使われる高電圧コネクタの絶縁設計にはポアソン方程式を解いた電位分布解析が欠かせません。充電中は数百アンペアの電流が流れるので、コネクタ内部の電位勾配が急峻になる箇所があり、そこが絶縁材の最弱点になります。設計者はポアソン方程式の解から「電界が集中する角部や段差」を特定し、面取り加工や絶縁材の厚みで対策します。実際の設計フローでは、まずポアソン方程式で電位分布を出し、次にその結果を熱解析の入力に使う連成解析になることも多い。1本の方程式が複数の設計フェーズにわたって活躍するんですね。
ポアソン方程式（静電場）のソフトウェア比較
ツール
全静電場ソルバーがポアソン方程式を解く。特別なツール選定は不要。COMSOL、Maxwell、FEMMのいずれでも。
ポアソン方程式ソルバーの「マルチグリッド法」が速い本当の理由
商用ツールのポアソン方程式ソルバーには「マルチグリッド法」が採用されているものが多くあります。なぜ速いかというと、通常の反復解法はメッシュの細かさに応じた「波長の短い誤差」は素早く消えるけど、「波長の長い誤差」は消えるのがとても遅いからです。マルチグリッド法は粗いメッシュと細かいメッシュを交互に使うことで、どんな波長の誤差も効率よく消せる。計算量が問題サイズNに対して $O(N)$ になるという驚異的な性能が実現できます。COMSOLやANSYSで「大規模静電解析が意外と速く終わった」というのは、こういうソルバーの工夫のおかげです。
ポアソン方程式（静電場）の先端研究
先端
ニューラルネットワークでポアソン方程式を解く——Physics-Informed NNの衝撃
2019年頃から「Physics-Informed Neural Network（PINN）」という手法が注目されています。ニューラルネットワークに支配方程式（ポアソン方程式を含む）の残差も損失関数に組み込んで学習させることで、境界値問題を解くというアプローチです。メッシュを作らなくていい、任意の点で予測できる、というのが大きなメリット。一方で収束の遅さや精度保証の難しさが課題で「FEMを完全に置き換えるのはまだ先」というのが現在の実情です。でも将来的に電池や半導体のデバイスシミュレーションで「大量の設計パラメータを高速にスキャン」する用途では、PINNが主流になる日が来るかもしれません。
ポアソン方程式（静電場）のトラブル対応
トラブル
「ソルバーが収束しない」——ポアソン方程式でよくある罠
ポアソン方程式のFEM解析で「反復解法が収束しない」というトラブルはよく起きます。原因の鉄板は「境界条件の設定ミス」。ディリクレ境界（電位固定）が一箇所もないと連立方程式の係数行列が特異になって解が定まりません。電位の絶対値が決まらず、電界だけが意味を持つ場合は「参照電位点を一点設けてアース」するのがお作法です。もう一つよくある罠が誘電率の単位ミス——相対誘電率 $\varepsilon_r$ を入れるべき場所に真空誘電率 $\varepsilon_0 = 8.85\times10^{-12}$ F/mをそのまま入れてしまって電位が12桁ずれる、という定番のミスがあります。
関連トピック
なった
詳しく
報告