ガウスの法則（静電界） — 理論と支配方程式
概要
先生、ガウスの法則ってMaxwell方程式の一つですよね？ 静電場解析でどう使うんですか？
ガウスの法則は「閉曲面を貫く電束の総量は、その内部の全電荷に等しい」という法則だ。Maxwell方程式の第一式に対応する。
微分形では
ここで $\mathbf{D} = \varepsilon \mathbf{E}$ は電束密度、$\rho$ は体積電荷密度だ。
これがFEMの静電場解析とどう関係するんですか？
$\mathbf{E} = -\nabla\phi$ を代入するとポアソン方程式になる。
FEMソルバーはこれを離散化して解く。つまりガウスの法則はFEM静電場解析の根幹にある方程式だ。
対称性を利用した解析解
ガウスの法則から直接電界を求められる場合もありますよね？
高い対称性がある場合に限って求まる。
| 対称性 | ガウス面 | 電界 |
|---|---|---|
| 球対称（点電荷） | 同心球面 | $E = Q/(4\pi\varepsilon r^2)$ |
| 円筒対称（線電荷） | 同軸円筒面 | $E = \lambda/(2\pi\varepsilon r)$ |
| 平面対称（面電荷） | 直方体 | $E = \sigma/(2\varepsilon)$ |
これらの解析解はFEM結果の妥当性検証に不可欠だ。COMSOLでもAnsys Maxwellでも、まず対称問題で理論値との一致を確認してから複雑な問題に進むのが鉄則だよ。
電気自動車モータ開発と電磁界解析
テスラのModel 3のモータは、リラクタンストルクと磁石トルクの両方を使うIPMSM（埋込磁石型同期モータ）。この複雑な磁場分布を最適化するには数千回の電磁界FEA解析が必要です。1回の解析に数分としても、最適化ループ全体では数週間のCPU時間。それでも実機を何十台も試作するよりは圧倒的に速くて安い。
各項の物理的意味
- 電場項 $\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仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値例：表皮深さの周波数依存性（銅, σ=5.96×10⁷ S/m, μr=1）
表皮深さ δ = √(2/(ωμσ))。周波数が上がるほど電流は表面に集中します：
GHz帯では電流はほぼ表面のみ！ 高周波回路で導体の表面粗さが性能に直結するのはこのためです。メッシュもδの1/3以下の要素サイズが必要です。
電磁界解析の精度と計算コストの両立は永遠の課題です。 — Project NovaSolverは、既存ワークフローの改善を目指す取り組みとして、この問題に向き合っています。
次世代CAEプロジェクト：開発者と実務者をつなぐ
Project NovaSolverは、ガウスの法則（静電界）を含む幅広い解析分野において、実務者の知見を最大限に活かせる環境の実現を探求しています。まだ道半ばですが、共に歩んでいただける方を募集しています。
開発パートナー登録 →