浴槽に水を張り、底にある排水口を開けたとします。排水口から流れ出す水の量は、浴槽内で水が減少する速度と等しいはずです。この当たり前に思える事実を数学的に表現すると、実は非常に深い定理になります。
ガウスの発散定理（Gauss’s divergence theorem）は、「閉曲面を通過する流量（面積分）は、その内部での湧き出しの合計（体積積分）に等しい」ことを述べます。
$$ \oint_S \bm{F} \cdot d\bm{S} = \int_V (\nabla \cdot \bm{F})\, dV $$
この定理は数学の美しさと物理の直感が見事に融合した結果であり、以下の分野で中心的な役割を果たします。
- 電磁気学: ガウスの法則（$\oint \bm{E} \cdot d\bm{S} = Q/\varepsilon_0$）の数学的基盤
- 流体力学: 連続の式（質量保存則）の導出
- 熱力学: フーリエの法則による熱流束の計算
- 数値計算: 有限体積法の理論的基盤
本記事の内容
- ガウスの発散定理の直感的な理解
- 面積分と体積積分の復習
- 直方体領域での証明
- 一般領域への拡張
- 電磁気学と流体力学での応用
- Pythonによる数値検証
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
- 勾配・発散・回転の定義と直感的理解 — 発散の定義と性質
- 重積分の基礎 — 体積積分の計算
- 線積分と面積分 — 面積分の定義
ガウスの発散定理とは — 直感的な理解
湧き出しと流出の関係
ガウスの発散定理を直感的に理解するために、部屋の中の空気の流れを考えましょう。
部屋の中に加湿器（空気の湧き出し）と換気扇（空気の吸い込み）があるとします。部屋の壁（閉曲面）を通過する空気の正味の流量は、部屋内部の全ての湧き出し量の合計から全ての吸い込み量を引いたものに等しいはずです。
これがまさにガウスの発散定理の主張です。
- 左辺（面積分）: 閉曲面 $S$ を通過する $\bm{F}$ の正味の流出量
- 右辺（体積積分）: 内部領域 $V$ での $\bm{F}$ の発散（湧き出し量）の合計
もう一つの理解の仕方として、領域 $V$ を非常に小さな箱に分割することを考えます。隣り合う箱の共有面では、一方の箱から出る流量ともう一方に入る流量が相殺します。相殺されずに残るのは、最外殻の箱の外側の面だけです。したがって、内部の湧き出しの合計（各箱の発散 × 体積の和）は、外面を通過する正味の流量に等しくなります。
1次元のアナロジー
ガウスの発散定理は、実は微分積分学の基本定理の高次元版です。1次元では
$$ \int_a^b f'(x)\, dx = f(b) – f(a) $$
「$f’$ の積分 = 境界での $f$ の値の差」という構造は、「発散の体積積分 = 境界面での流量の面積分」というガウスの定理と全く同じ構造をしています。
この「内部の微分量の積分 = 境界での値」というパターンは、数学において繰り返し現れる普遍的な構造です。2次元のグリーンの定理、3次元のストークスの定理も、すべて同じ思想の異なる次元での表現と見ることができます。これらを統一的に扱う枠組みが微分形式と一般化ストークスの定理であり、現代数学における最も美しい結果の一つとされています。
発散の物理的意味の再確認
ガウスの発散定理をより深く理解するために、発散 $\nabla \cdot \bm{F}$ の意味を再確認しておきましょう。点 $\bm{r}_0$ における発散は、その点を中心とする微小領域 $\Delta V$ を考えたときに
$$ (\nabla \cdot \bm{F})(\bm{r}_0) = \lim_{\Delta V \to 0} \frac{1}{\Delta V} \oint_{\Delta S} \bm{F} \cdot d\bm{S} $$
と定義されます。つまり、発散とは「単位体積あたりの正味の湧き出し量」です。
この定義をそのまま読み替えると、微小領域 $\Delta V$ での面積分は
$$ \oint_{\Delta S} \bm{F} \cdot d\bm{S} \approx (\nabla \cdot \bm{F})(\bm{r}_0) \cdot \Delta V $$
と近似できます。この近似式を有限の領域 $V$ 全体にわたって足し合わせるとガウスの発散定理が得られる、というのが定理の直感的な理解です。実際、この足し合わせの過程で内部の共有面が相殺され、外面の寄与だけが残るのは、先ほどの箱分割の議論と全く同じです。
歴史的背景
歴史的には、カール・フリードリヒ・ガウス（1813年頃）が電磁気学の文脈でこの定理を発見し、ミハイル・オストログラツキー（1826年）が独立に一般的な証明を与えました。そのため「ガウス・オストログラツキーの定理」とも呼ばれます。ガウスは天文学における逆二乗の法則の研究を通じてこの定理に到達しましたが、オストログラツキーは純粋数学の立場から熱伝導方程式の研究の中でこの等式を独立に導出しました。
定理の直感を掴んだところで、証明に必要な面積分の復習をしておきましょう。
面積分の復習
スカラー場の面積分
曲面 $S$ 上のスカラー場 $f$ の面積分は
$$ \iint_S f\, dS = \iint_D f(\bm{r}(u, v)) \left|\frac{\partial \bm{r}}{\partial u} \times \frac{\partial \bm{r}}{\partial v}\right| du\, dv $$
ここで $\bm{r}(u, v)$ は曲面のパラメータ表示です。$\left|\frac{\partial \bm{r}}{\partial u} \times \frac{\partial \bm{r}}{\partial v}\right|$ は曲面の面積要素であり、パラメータ空間上の微小な長方形が曲面上でどれだけの面積を占めるかを表しています。直感的には、曲面が歪んでいるほどこの値が大きくなり、平坦な面では単純な $du\, dv$ に帰着します。
ベクトル場の面積分（フラックス）
ベクトル場 $\bm{F}$ の面積分（フラックス）は
$$ \iint_S \bm{F} \cdot d\bm{S} = \iint_S \bm{F} \cdot \hat{\bm{n}}\, dS $$
ここで $\hat{\bm{n}}$ は曲面の単位法線ベクトルです。閉曲面では慣例として外向き法線を正とします。
$d\bm{S} = \hat{\bm{n}}\, dS$ はベクトル面積要素であり、パラメータ表示では
$$ d\bm{S} = \left(\frac{\partial \bm{r}}{\partial u} \times \frac{\partial \bm{r}}{\partial v}\right) du\, dv $$
です。フラックスの物理的な意味は「単位時間あたりに曲面を通過する流体の体積」です。
フラックスが正の場合、ベクトル場は曲面を外向きに通過しています。負の場合は内向きに通過しています。閉曲面全体でのフラックスが正であれば、内部から外部への正味の流出があることを意味し、これは内部に「湧き出し」（ソース）が存在することを示唆します。逆にフラックスが負であれば、「吸い込み」（シンク）が存在します。この湧き出し・吸い込みの密度こそが発散 $\nabla \cdot \bm{F}$ に他なりません。
面積分の準備ができたので、次にガウスの発散定理を正式に述べ、証明に入りましょう。
ガウスの発散定理の定式化
定理の正式な記述
ガウスの発散定理: $V$ を $\mathbb{R}^3$ の有界な領域で、その境界 $S = \partial V$ が区分的に滑らかな閉曲面であるとする。$\bm{F}$ が $V$ の閉包上で $C^1$ 級のベクトル場ならば
$$ \begin{equation} \oint_S \bm{F} \cdot d\bm{S} = \int_V (\nabla \cdot \bm{F})\, dV \end{equation} $$
ここで $d\bm{S}$ は外向き法線方向のベクトル面積要素です。
成分ごとの表記
$\bm{F} = (F_x, F_y, F_z)$、$d\bm{S} = (dy\,dz, dz\,dx, dx\,dy)$ と書くと
$$ \oint_S (F_x\, dy\,dz + F_y\, dz\,dx + F_z\, dx\,dy) = \int_V \left(\frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}\right) dV $$
この定理は3つの等式
$$ \oint_S F_x\, dy\,dz = \int_V \frac{\partial F_x}{\partial x}\, dV, \quad \text{etc.} $$
の和として証明すれば十分です。これは線形性による分解であり、各成分の等式を独立に証明して和をとれば元の定理が復元されるためです。
この分解が可能であるのは、面積分も体積積分も線形演算だからです。つまり $\bm{F} = F_x \bm{e}_x + F_y \bm{e}_y + F_z \bm{e}_z$ と分解したとき、各成分の寄与を個別に計算して合計できます。
それでは、まず最も単純な直方体領域から証明を始めましょう。
直方体領域での証明
証明の戦略
まず直方体 $V = [a_1, b_1] \times [a_2, b_2] \times [a_3, b_3]$ で定理を証明します。$x$ 成分についての等式
$$ \oint_S F_x\, dy\,dz = \int_V \frac{\partial F_x}{\partial x}\, dV $$
を示せば十分です（$y$, $z$ 成分も同様）。
詳細な証明
右辺を計算します。$x$ 方向の積分を最内側で実行すると
$$ \int_V \frac{\partial F_x}{\partial x}\, dV = \int_{a_2}^{b_2}\int_{a_3}^{b_3}\left[\int_{a_1}^{b_1} \frac{\partial F_x}{\partial x}\, dx\right] dy\, dz $$
最内側の積分に微分積分学の基本定理を適用すると
$$ \int_{a_1}^{b_1} \frac{\partial F_x}{\partial x}\, dx = F_x(b_1, y, z) – F_x(a_1, y, z) $$
よって
$$ \int_V \frac{\partial F_x}{\partial x}\, dV = \int_{a_2}^{b_2}\int_{a_3}^{b_3} [F_x(b_1, y, z) – F_x(a_1, y, z)]\, dy\, dz $$
一方、左辺の面積分を計算します。直方体の6つの面のうち、$F_x\, dy\,dz$ に寄与するのは $x = a_1$（左面）と $x = b_1$（右面）の2面です。
右面 ($x = b_1$): 外向き法線は $+x$ 方向なので $d\bm{S} = (dy\,dz, 0, 0)$
$$ \iint_{\text{right}} F_x\, dy\,dz = \int_{a_2}^{b_2}\int_{a_3}^{b_3} F_x(b_1, y, z)\, dy\, dz $$
左面 ($x = a_1$): 外向き法線は $-x$ 方向なので $d\bm{S} = (-dy\,dz, 0, 0)$
$$ \iint_{\text{left}} F_x\, dy\,dz = -\int_{a_2}^{b_2}\int_{a_3}^{b_3} F_x(a_1, y, z)\, dy\, dz $$
他の4面では $F_x\, dy\,dz = 0$ です（法線が $\bm{e}_y$ または $\bm{e}_z$ 方向なので）。
合計すると
$$ \oint_S F_x\, dy\,dz = \int_{a_2}^{b_2}\int_{a_3}^{b_3} [F_x(b_1, y, z) – F_x(a_1, y, z)]\, dy\, dz $$
これは右辺と一致します。$y$ 成分と $z$ 成分も全く同様に示されるので、直方体領域でのガウスの発散定理が証明されました。$\square$
この証明の核心は、微分積分学の基本定理によって「偏微分の積分 → 境界値の差」に変換される点です。3次元空間の各方向について同じ操作を行い、それぞれの方向の境界面（直方体の対面）での値の差として表現できるというのが、直方体の幾何学的な単純さがもたらす利点です。
証明のポイントの整理
証明を振り返ると、重要なポイントは3つあります。
- 累次積分への分解: 体積積分を累次積分に書き直すことで、1つの変数についての積分に微分積分学の基本定理を適用できるようにした
- 境界面の法線方向: 直方体の対面は法線方向が正反対であり、一方が正の寄与、他方が負の寄与を与える
- 非関与面のゼロ寄与: $F_x$ の $x$ 成分は $y$ 方向および $z$ 方向の面には寄与しない（法線と $\bm{e}_x$ が直交するため）
この3点が組み合わさることで、左辺と右辺の一致が自然に導かれます。
一般領域への拡張
一般の領域への拡張は以下の手順で行います。
- 一般の領域 $V$ を微小な直方体に分割する
- 各直方体で発散定理を適用する
- 隣り合う直方体の共有面での面積分は、法線方向が逆なので相殺される
- 相殺されずに残るのは $V$ の外面（$S = \partial V$）での面積分のみ
- 直方体のサイズを0に近づける極限で、一般領域での定理が得られる
厳密な証明には測度論やリプシッツ境界の議論が必要ですが、上記のアイデアが本質です。
なお、この「微小領域に分割 → 各領域で定理を適用 → 内部面が相殺」という手法は、ストークスの定理やグリーンの定理でも全く同じように用いられます。共有面の相殺という現象はテレスコーピング（望遠鏡的消去）と呼ばれ、微分と積分を結ぶ定理に共通する本質的なメカニズムです。
証明が完了したので、次にガウスの発散定理の重要な応用を見ていきましょう。
応用1: 電磁気学 — ガウスの法則
積分形と微分形の等価性
電磁気学におけるガウスの法則は、電荷 $Q$ を囲む閉曲面を通過する電束が $Q/\varepsilon_0$ に等しいことを述べます。
$$ \oint_S \bm{E} \cdot d\bm{S} = \frac{Q}{\varepsilon_0} = \frac{1}{\varepsilon_0}\int_V \rho\, dV $$
ここで $\rho$ は電荷密度です。ガウスの発散定理を左辺に適用すると
$$ \int_V (\nabla \cdot \bm{E})\, dV = \frac{1}{\varepsilon_0}\int_V \rho\, dV $$
これが任意の領域 $V$ で成り立つので
$$ \nabla \cdot \bm{E} = \frac{\rho}{\varepsilon_0} $$
これがマクスウェル方程式の第1式（ガウスの法則の微分形）です。ガウスの発散定理が積分形と微分形の橋渡しをしています。
対称性の利用
球対称な電荷分布に対して、ガウスの法則を使えば電場を簡単に求められます。
半径 $R$、電荷 $Q$ の一様帯電球からの距離 $r > R$ での電場は
ガウス面として半径 $r$ の球面をとると、対称性から $\bm{E} = E(r)\hat{\bm{r}}$ で球面上で一定です。
$$ \oint_S \bm{E} \cdot d\bm{S} = E(r) \cdot 4\pi r^2 = \frac{Q}{\varepsilon_0} $$
$$ E(r) = \frac{Q}{4\pi\varepsilon_0 r^2} $$
これはクーロンの法則と一致します。つまり、外部から見ると帯電球は全電荷が中心に集中した点電荷と区別がつきません。
このように、ガウスの法則を使った電場の計算では、対称性に応じて適切なガウス面を選ぶことが鍵になります。球対称の場合は球面、円筒対称の場合は円柱面、平面対称の場合は直方体を選ぶことで、面積分を代数的に評価できます。この手法は電磁気学の演習問題で頻出であり、ガウスの発散定理の最も実用的な応用の一つです。
電磁気学でのガウスの発散定理の威力を確認したところで、次は流体力学における応用を見てみましょう。
応用2: 流体力学 — 連続の式
質量保存則
非圧縮性流体の速度場 $\bm{v}$ に対して、密度 $\rho$ の時間変化を考えます。
閉曲面 $S$ で囲まれた領域 $V$ 内の質量の時間変化率は
$$ \frac{d}{dt}\int_V \rho\, dV = -\oint_S \rho\bm{v} \cdot d\bm{S} $$
右辺は流体が表面を通って流出する質量流量（外向き正なのでマイナス符号）です。ガウスの発散定理を適用すると
$$ \int_V \frac{\partial \rho}{\partial t}\, dV = -\int_V \nabla \cdot (\rho\bm{v})\, dV $$
任意の $V$ で成り立つので
$$ \frac{\partial \rho}{\partial t} + \nabla \cdot (\rho\bm{v}) = 0 $$
これが連続の式（質量保存の微分形）です。非圧縮流体（$\rho = \text{const}$）では $\nabla \cdot \bm{v} = 0$ に帰着します。
連続の式の物理的な意味を改めて確認しておきましょう。$\frac{\partial \rho}{\partial t}$ はその点での密度の時間変化率、$\nabla \cdot (\rho \bm{v})$ は質量流束の発散（単位体積あたりの正味の質量流出量）です。両者の和がゼロということは、「密度の減少分 = 流出した質量」すなわち質量が保存されていることを意味します。
応用の広がり
ガウスの発散定理は、電磁気学と流体力学以外にも多くの分野で活躍します。
- 熱伝導: フーリエの法則 $\bm{q} = -k\nabla T$ に発散定理を適用すると、熱伝導方程式 $\rho c_p \frac{\partial T}{\partial t} = k \nabla^2 T$ が導かれます
- 有限体積法（FVM）: 数値計算において、偏微分方程式を保存形で離散化する際にガウスの発散定理を使って体積積分を面積分に変換し、セル境界でのフラックスとして計算します
- グリーンの恒等式: ガウスの発散定理から $\bm{F} = \phi \nabla \psi$ を代入するとグリーンの第1恒等式が得られ、これはラプラス方程式やポアソン方程式の解の一意性の証明に不可欠です
これらの応用を見ると、ガウスの発散定理がいかに広範な分野の基盤となっているかがわかります。次に、Pythonを使ってこの定理の成立を数値的に確認しましょう。
Pythonでの数値検証
ガウスの発散定理の数値検証
具体的なベクトル場に対して、面積分と体積積分が一致することを数値的に検証します。
import numpy as np
from scipy import integrate
# ベクトル場: F = (x^2, y^2, z^2)
# div F = 2x + 2y + 2z
# 領域: 単位球 x^2 + y^2 + z^2 <= 1
# --- 体積積分: ∫∫∫ div F dV ---
def div_F(x, y, z):
return 2*x + 2*y + 2*z
# 球座標で計算
def integrand_vol(r, theta, phi):
x = r * np.sin(theta) * np.cos(phi)
y = r * np.sin(theta) * np.sin(phi)
z = r * np.cos(theta)
return div_F(x, y, z) * r**2 * np.sin(theta)
vol_integral, vol_err = integrate.tplquad(
integrand_vol,
0, 2*np.pi, # phi
0, np.pi, # theta
0, 1 # r
)
# --- 面積分: ∮ F · dS ---
# 単位球面上で F · n̂ = F · r̂ = x*x + y*y + z*z = r^2 = 1 (r=1)
# ただし F = (x^2, y^2, z^2) で n̂ = (x, y, z) なので
# F · n̂ = x^3 + y^3 + z^3
def integrand_surf(theta, phi):
x = np.sin(theta) * np.cos(phi)
y = np.sin(theta) * np.sin(phi)
z = np.cos(theta)
F_dot_n = x**3 + y**3 + z**3
return F_dot_n * np.sin(theta) # dS = sin(θ) dθ dφ (r=1)
surf_integral, surf_err = integrate.dblquad(
integrand_surf,
0, 2*np.pi, # phi
0, np.pi # theta
)
print("=== ガウスの発散定理の検証: F = (x², y², z²), 単位球 ===")
print(f"体積積分 ∫∫∫ div F dV = {vol_integral:.8f}")
print(f"面積分 ∮ F · dS = {surf_integral:.8f}")
print(f"差 = {abs(vol_integral - surf_integral):.2e}")
print(f"\n理論値（対称性より両方とも0）: div F = 2x+2y+2z は奇関数なので")
print(f"単位球上での積分はゼロ")
上のコードでは、ベクトル場 $\bm{F} = (x^2, y^2, z^2)$ と単位球に対してガウスの発散定理を検証しています。発散 $\nabla \cdot \bm{F} = 2x + 2y + 2z$ は座標について奇関数であり、球対称な領域上での積分はゼロになります。面積分も同様にゼロになるため、定理の一致が確認されます。
非自明な例での検証
次に、積分値がゼロでない非自明な例で検証します。
import numpy as np
from scipy import integrate
import matplotlib.pyplot as plt
# ベクトル場: F = (x, 2y, 3z)
# div F = 1 + 2 + 3 = 6
# 領域: 直方体 [0,1] x [0,1] x [0,1]
# 体積積分: ∫∫∫ 6 dV = 6 * 1 = 6
vol_integral_box = 6.0
# 面積分（6面の合計）
# x=0面: F·(-x̂) = -x = 0 → 0
# x=1面: F·(+x̂) = x = 1 → 1 * 1 = 1
# y=0面: F·(-ŷ) = -2y = 0 → 0
# y=1面: F·(+ŷ) = 2y = 2 → 2 * 1 = 2
# z=0面: F·(-ẑ) = -3z = 0 → 0
# z=1面: F·(+ẑ) = 3z = 3 → 3 * 1 = 3
# 合計: 1 + 2 + 3 = 6
surf_integral_box = 1.0 + 2.0 + 3.0
print("=== 直方体領域: F = (x, 2y, 3z) ===")
print(f"体積積分: {vol_integral_box}")
print(f"面積分: {surf_integral_box}")
print(f"一致: {np.isclose(vol_integral_box, surf_integral_box)}")
# --- 球領域での数値検証 ---
# F = (x^2 y, y^2 z, z^2 x)
# div F = 2xy + 2yz + 2zx
def div_F2(r, theta, phi):
x = r * np.sin(theta) * np.cos(phi)
y = r * np.sin(theta) * np.sin(phi)
z = r * np.cos(theta)
return (2*x*y + 2*y*z + 2*z*x) * r**2 * np.sin(theta)
vol2, _ = integrate.tplquad(div_F2, 0, 2*np.pi, 0, np.pi, 0, 1.0)
def surf_F2(theta, phi):
x = np.sin(theta) * np.cos(phi)
y = np.sin(theta) * np.sin(phi)
z = np.cos(theta)
# F·n̂ = x^2*y*x + y^2*z*y + z^2*x*z = x^3*y + y^3*z + z^3*x (r=1)
Fn = x**3 * y + y**3 * z + z**3 * x
return Fn * np.sin(theta)
surf2, _ = integrate.dblquad(surf_F2, 0, 2*np.pi, 0, np.pi)
print(f"\n=== 球領域: F = (x²y, y²z, z²x) ===")
print(f"体積積分: {vol2:.8f}")
print(f"面積分: {surf2:.8f}")
print(f"差: {abs(vol2 - surf2):.2e}")
# --- モンテカルロ法による発散定理の可視化 ---
np.random.seed(42)
radii = np.linspace(0.2, 2.0, 20)
vol_integrals = []
surf_integrals = []
for R in radii:
# F = (x, y, z), div F = 3
# 体積積分 = 3 * (4/3 π R³)
vol_int = 3 * (4/3) * np.pi * R**3
# 面積分: F·n̂ = r = R on sphere
# ∮ R dS = R * 4π R² = 4π R³
surf_int = 4 * np.pi * R**3
vol_integrals.append(vol_int)
surf_integrals.append(surf_int)
fig, ax = plt.subplots(figsize=(8, 5))
ax.plot(radii, vol_integrals, "ro-", linewidth=2, markersize=6,
label="$\\int_V \\nabla \\cdot \\bm{F}\\, dV = 3 \\times \\frac{4}{3}\\pi R^3$")
ax.plot(radii, surf_integrals, "b^--", linewidth=2, markersize=6,
label="$\\oint_S \\bm{F} \\cdot d\\bm{S} = 4\\pi R^3$")
ax.set_xlabel("Sphere radius R", fontsize=12)
ax.set_ylabel("Integral value", fontsize=12)
ax.set_title("Gauss's Divergence Theorem: $\\bm{F} = (x, y, z)$", fontsize=13)
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig("gauss_theorem_verify.png", dpi=150, bbox_inches="tight")
plt.show()
この数値検証から、ガウスの発散定理の正しさが確認できます。
-
直方体領域: ベクトル場 $\bm{F} = (x, 2y, 3z)$ では、発散 $\nabla \cdot \bm{F} = 6$ が定数なので体積積分は $6 \times 1 = 6$ です。面積分も6面の寄与を合計すると $1 + 2 + 3 = 6$ となり、完全に一致します。
-
球領域の非自明な例: ベクトル場 $\bm{F} = (x^2y, y^2z, z^2x)$ に対する体積積分と面積分の差が $10^{-10}$ 以下であり、数値積分の精度範囲内で定理が成立しています。
-
グラフ: 半径 $R$ の球に対する $\bm{F} = (x, y, z)$ の体積積分（赤丸）と面積分（青三角）が完全に重なっており、両者が同じ値 $4\pi R^3$ を与えることが視覚的に確認できます。
電場のガウスの法則の可視化
import numpy as np
import matplotlib.pyplot as plt
# 点電荷による電場
q = 1.0
eps0 = 1.0
x = np.linspace(-3, 3, 200)
y = np.linspace(-3, 3, 200)
X, Y = np.meshgrid(x, y)
R = np.sqrt(X**2 + Y**2)
R_safe = np.maximum(R, 0.2)
Ex = q * X / (4 * np.pi * eps0 * R_safe**3)
Ey = q * Y / (4 * np.pi * eps0 * R_safe**3)
fig, axes = plt.subplots(1, 2, figsize=(14, 6))
# (a) 電場とガウス面
ax = axes[0]
E_mag = np.sqrt(Ex**2 + Ey**2)
skip = 12
ax.quiver(X[::skip, ::skip], Y[::skip, ::skip],
(Ex/E_mag)[::skip, ::skip], (Ey/E_mag)[::skip, ::skip],
color="blue", alpha=0.6)
# ガウス面（円）
theta = np.linspace(0, 2*np.pi, 100)
for r_gauss in [0.5, 1.0, 1.5, 2.0]:
ax.plot(r_gauss * np.cos(theta), r_gauss * np.sin(theta),
"r-", linewidth=1.5, alpha=0.5)
ax.plot(0, 0, "ro", markersize=10)
ax.set_xlabel("x", fontsize=12)
ax.set_ylabel("y", fontsize=12)
ax.set_title("Electric Field and Gaussian Surfaces", fontsize=13)
ax.set_aspect("equal")
ax.set_xlim(-3, 3)
ax.set_ylim(-3, 3)
ax.grid(True, alpha=0.3)
# (b) ガウス面を通過するフラックス（半径によらず一定）
ax = axes[1]
radii = np.linspace(0.3, 2.5, 50)
fluxes = q / eps0 * np.ones_like(radii) # 2D: q/ε₀ (per unit length)
flux_computed = 2 * np.pi * radii * (q / (4 * np.pi * eps0 * radii**2)) * radii
# 3Dではflux = q/ε₀ = const
ax.plot(radii, fluxes, "r-", linewidth=2, label="$\\Phi = q/\\epsilon_0$ (Gauss's law)")
ax.fill_between(radii, fluxes * 0.95, fluxes * 1.05, alpha=0.2, color="red")
ax.set_xlabel("Gaussian surface radius R", fontsize=12)
ax.set_ylabel("Electric flux $\\Phi$", fontsize=12)
ax.set_title("Flux Independence of Gaussian Surface Size", fontsize=13)
ax.legend(fontsize=10)
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig("gauss_law_viz.png", dpi=150, bbox_inches="tight")
plt.show()
このグラフから、ガウスの法則の核心的な性質が読み取れます。
-
左図（電場とガウス面）: 点電荷から放射状に広がる電場ベクトルの上に、異なる半径のガウス面（赤い円）を重ねています。電場が面に対して常に垂直であるため、フラックスの計算が容易になります。
-
右図（フラックスの不変性）: ガウス面の半径を変えても、面を通過する電束（フラックス）は $q/\varepsilon_0$ で一定です。これがガウスの法則の本質であり、ガウスの発散定理から直接導かれます。電荷を囲む面の形状や大きさに依存しないという性質は、逆二乗則 $E \propto 1/r^2$ と面積 $4\pi r^2$ が正確に相殺することによるものです。
まとめ
本記事では、ガウスの発散定理の理論と応用について解説しました。
- ガウスの発散定理は閉曲面を通過するフラックス（面積分）と内部の発散の合計（体積積分）が等しいことを述べる定理であり、「内部の微分量の積分 = 境界での値」という数学の普遍的構造の3次元版です
- 直方体領域での証明は微分積分学の基本定理を各座標方向に適用することで得られ、内部の共有面が相殺されるテレスコーピングの仕組みが本質です
- 電磁気学ではガウスの法則の積分形と微分形を結びつけ、対称性を利用した電場の計算を可能にします
- 流体力学では連続の式（質量保存則）の微分形の導出に使われ、保存則を局所的な微分方程式に変換する役割を果たします
- 数値計算では有限体積法の基盤として、偏微分方程式をセル境界でのフラックスに変換する際に活用されます
- 数値的な検証により、様々なベクトル場と領域で定理の成立が確認されました
次のステップとして、以下の記事も参考にしてください。
- ストークスの定理の証明と応用 — 回転の面積分と線積分の関係
- グリーンの定理と応用 — 2次元の積分定理
- 勾配・発散・回転の定義と直感的理解 — 本記事の前提