ガウスの法則（静電界） — 数値解法と実装
数値解法の詳細
ガウスの法則から導かれるポアソン方程式をFEMでどう定式化するんですか？
Galerkin法による弱形式化が基本だ。試験関数 $w$ を掛けて部分積分すると
離散化すると $[K]\{\phi\} = \{f\}$。構造FEMの剛性方程式と全く同じ形だ。
誘電率がヤング率に対応するわけですね。
その通り。だから構造FEMの経験があれば静電場FEMは理解しやすい。
ガウスの法則による電荷計算
FEM結果から導体上の電荷を後処理で求めるには？
ガウスの法則を使って導体表面の電束密度を積分する。
COMSOLでは表面積分の後処理機能で直接計算可能。この値から静電容量 $C = Q/V$ を算出する。
容量行列の抽出
多導体系の容量行列はどうやって求めるんですか？
導体 $j$ を1Vにして他をすべて0Vに設定し、各導体上の誘導電荷を計算する。これを全導体について繰り返せば容量行列 $C_{ij}$ が得られる。Ansys Q3Dではこの操作が完全に自動化されていて、PCB配線の寄生容量抽出で広く使われている。
ファラデー——「数学が苦手だった」天才
電磁誘導の法則を発見したマイケル・ファラデーは、正規の教育を受けておらず、高等数学が使えませんでした。彼は「力線」という直感的なイメージで電磁気現象を理解し、実験で次々と発見をしました。後にマクスウェルがファラデーの直感を数学で定式化したのがマクスウェル方程式です。CAEの数式の裏には、常に「物理的な直感」があることを忘れずに。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
辺要素（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は、ガウスの法則（静電界）における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →