【要約】【Python】ガウスの発散定理を3D可視化してみた [Zenn_Python] | Summary by TechDistill
> Source: Zenn_Python
Execute Primary Source
// Problem
著者は、ガウスの発散定理が数式のみでは物理的なイメージを掴みにくいという課題に着目した。数式としての理解はできても、実際の空間での挙動を直感的に捉えることは困難である。
- ・抽象的な数式により、空間内の流れと境界の挙動の関係が理解しにくい。
- ・電磁気学や流体力学の実務において、理論の物理的意味を誤解するリスクがある。
- ・連続的な現象をコンピュータで扱うための、離散化プロセスがブラックボックス化しやすい。
// Approach
著者は、連続的な数式を離散的なメッシュに変換し、コンピュータによる数値計算を行う手法を採用した。理論を理解する最良の方法として、数式を離散化し、力づくで足し合わせるプロセスを実装している。
- ・Matplotlibを用い、ベクトル場と閉曲面を3D空間に描画する。
- ・空間を立方体メッシュに分割し、領域内の発散の総和（体積分）を算出する。
- ・球面を緯度・経度方向に分割し、法線ベクトルとの内積を用いた面積分を算出する。
- ・「湧き出し」と「渦巻き」の2パターンで、理論値との一致を検証する。
// Result
著者は、Pythonを用いた数値計算により、理論値と計算値の整合性を確認した。これにより、数式が物理的な現象として正しく記述されていることを数値的に証明した。
- ・湧き出しの場において、体積分と面積分の差は0.00055という極小値に留まった。
- ・渦巻きの場では、発散が0であることと、法線ベクトルとの直交性が数値的に示された。
- ・分割数を増やすことで、計算結果が理論値へ収束するプロセスを実証した。
Senior Engineer Insight
> 本記事は、物理現象を計算機で扱うための「離散化」の本質を突いている。シミュレーション実務において、理論値と数値解の乖離（誤差）を制御する重要性を再認識させる。FEMやSPH法などの高度な数値解析手法の基礎となる思考プロセスが示されている。また、可視化による直感的な理解は、複雑な物理モデルのデバッグにおいて極めて有効である。