ポアソン方程式（静電場） — 数値解法と実装
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]$。
まとめ
- 剛性マトリクス$[K] = \int \varepsilon [B]^T[B] d\Omega$ — 構造と同形
- 1次三角形要素で手計算可能 — FEMの教育に最適
電気自動車モータ開発と電磁界解析
テスラのModel 3のモータは、リラクタンストルクと磁石トルクの両方を使うIPMSM（埋込磁石型同期モータ）。この複雑な磁場分布を最適化するには数千回の電磁界FEA解析が必要です。1回の解析に数分としても、最適化ループ全体では数週間のCPU時間。それでも実機を何十台も試作するよりは圧倒的に速くて安い。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| 直接法 | 小〜中規模の電磁場問題に適する。LU分解が主流。周波数掃引で右辺のみ変更可能な場合に効率的。 |
| 反復法（GMRES、BiCGSTAB） | 大規模3D電磁場問題で必須。前処理にはILU(k)やAMGが有効。複素対称系にはCOCG法も選択肢。 |
| DOF別推奨 | 〜5×10⁴ DOF: 直接法、5×10⁴〜: GMRES+ILU前処理、10⁶ DOF〜: 並列GMRES+AMG |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
数値解法の直感的理解
辺要素のイメージ
電磁界解析で使われる辺要素（エッジ要素）は「道路に沿った交通量」をイメージすると分かりやすい。通常の節点要素が「交差点の値」を扱うのに対し、辺要素は「辺（道路）に沿った量」を直接扱う。電場や磁場は方向を持つベクトル量であり、辺に沿った成分として自然に表現できるため、スプリアス解（非物理的な偽の解）を自動的に排除できる。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
電磁界解析の精度と計算コストの両立は永遠の課題です。 — Project NovaSolverは、既存ワークフローの改善を目指す取り組みとして、この問題に向き合っています。
CAEの未来を、実務者と共に考える
Project NovaSolverは、ポアソン方程式（静電場）における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →