面積分と流束
mathvector-calculussurface-integrallecture
1導入
この講義では、曲面上のスカラー量の総和と、曲面を貫くベクトル場の法線成分の総和を区別して定義する。
2用語と定義
流束 は、ベクトル場の法線成分を積分する面積分であり、曲面をどれだけ通過するかを表す。
面積分 は、曲面上のスカラー量またはベクトル場の法線成分を積分する操作である。
3方針
曲面を \boldsymbol{r}(u,v) で表示し、接ベクトル \boldsymbol{r}_u,\boldsymbol{r}_v から法線方向を構成する。流束では、ベクトル場のうち法線方向の成分だけを数える。
4向き付き曲面
流束では、曲面の表裏を区別する必要がある。法線方向を反転すると、流束の符号も反転する。閉曲面では、通常は外向き法線を選択する。
曲面全体で連続な単位法線場 \boldsymbol n を選べるとき、その曲面を向き付け可能という。連結な曲面では 2 つの向きは \boldsymbol n と -\boldsymbol n である。向き付け不可能な曲面では、曲面全体に一貫した流束を定義できない。領域の境界として与えられた閉曲面では、外向きが標準の向きである。
5厳密な説明
\Omega\subset\mathbb R^2 を境界が区分的 C^1 の有界 Jordan 領域、D=\overline\Omega とし、\boldsymbol r を D の開近傍で C^1 とする。\Omega 内で \boldsymbol r_u\times\boldsymbol r_v\ne\boldsymbol0 かつ単射であるとき、S=\boldsymbol r(D) を正則曲面パッチという。境界上だけの重複は面積 0 なので積分値へ影響しない。内部で多重被覆するパラメータ表示の右辺は、曲面上の通常の積分ではなく、重複度を含む積分を表す。
S の近傍で定義された連続スカラー場 g の面積分を
\iint_S g\,dS
=\iint_D g(\boldsymbol r(u,v))\|\boldsymbol r_u\times\boldsymbol r_v\|\,du\,dv
と定義する。これは向きに依存せず、g=1 なら曲面の面積になる。
S の近傍で定義された連続ベクトル場 \boldsymbol F に対し、\boldsymbol r_u\times\boldsymbol r_v が選んだ向きと一致するとき、流束を
\iint_S \boldsymbol{F}\cdot \boldsymbol{n}\,dS
=\iint_D \boldsymbol{F}(\boldsymbol{r}(u,v))\cdot(\boldsymbol{r}_u\times\boldsymbol{r}_v)\,du\,dv
と定義する。パッチの内部では、単位法線は \boldsymbol n=(\boldsymbol r_u\times\boldsymbol r_v)/\|\boldsymbol r_u\times\boldsymbol r_v\| である。向きを反転すると流束の符号も反転する。
\boldsymbol{r}_u\times\boldsymbol{r}_v は、接平面に垂直な方向と微小面積を同時に表す。\|\boldsymbol{r}_u\times\boldsymbol{r}_v\|\,du\,dv は面積要素であり、向き付き流束では絶対値を取らずに \boldsymbol{r}_u\times\boldsymbol{r}_v\,du\,dv を使用する。
本講義で扱うコンパクトな区分的正則曲面を 1 個のパッチで正則かつ 1 対 1 に表示できない場合は、有限個の正則パッチに分割し、共通部分を面積 0 の境界に限定して、各パッチの積分を加算する。流束では、すべてのパッチの法線を選択した曲面の向きと整合させる。この構成により、後述する球面および閉曲面の積分も定義される。
\Omega' を境界が区分的 C^1 の別の有界 Jordan 領域、E=\overline{\Omega'} とする。\boldsymbol\varphi:E\to D を、E と D の開近傍間の C^1 微分同相の制限で \boldsymbol\varphi(E)=D を満たすものとし、\boldsymbol\rho=\boldsymbol r\circ\boldsymbol\varphi とする。スカラー面積分は Jacobian の絶対値により再パラメータ化で不変である。流束は \det D\boldsymbol\varphi>0 なら向きを保って不変である。\det D\boldsymbol\varphi<0 のとき、外積をそのまま使用すると反対向きの曲面を表すため、流束の符号が反転する。選択済みの幾何学的な向きを維持する場合は、外積の符号を補正する。
6具体計算
上向きの平面 z=0、0\le x\le 1,\ 0\le y\le 1 を考える。\boldsymbol{F}=(0,0,z+1) なら、法線は (0,0,1) であり、流束は
\iint_S \boldsymbol{F}\cdot\boldsymbol{n}\,dS=\int_0^1\int_0^1 1\,dx\,dy=1
である。下向きの法線を選択すると値は -1 になる。
7球面の例
\boldsymbol{F}=(x,y,z) とし、原点を中心とする半径 R>0 の球面を外向きに向き付ける。球面では \boldsymbol{F}\cdot\boldsymbol{n}=R であり、面積が 4\pi R^2 なので、流束は 4\pi R^3 になる。後続講義では Gauss 定理を使い、発散が 3、体積が 4\pi R^3/3 であることから同じ値を確認する。
8円柱の例
半径 R、高さ H の円柱側面を
\boldsymbol{r}(\theta,z)=(R\cos\theta,R\sin\theta,z)
で表示する。これは上下面を含まない側面だけの表示で、0\le\theta\le 2\pi、0\le z\le H である。\theta=0 と \theta=2\pi は同一の継ぎ目を表すが、パラメータ領域の境界での重複なので面積分の値に影響しない。外向きの向き付き面要素は
\boldsymbol{r}_\theta\times\boldsymbol{r}_z=(R\cos\theta,R\sin\theta,0)
である。\boldsymbol{F}=(x,y,0) なら、流束は
\int_0^H\int_0^{2\pi}R^2\,d\theta\,dz=2\pi R^2H
である。
9比較例: 面に平行な場
z=0 の上向き平面片に対し、\boldsymbol{F}=(1,0,0) は面に平行である。したがって \boldsymbol{F}\cdot\boldsymbol{n}=0 であり、流束は 0 になる。矢印が曲面上に存在しても、法線成分がなければ貫通量は発生しない。
10Gauss 定理への橋渡し
有界領域 V が区分的 C^1 境界を持ち、\boldsymbol F が \overline V の近傍で C^1 であるとする。外向きに向き付けた閉曲面 \partial V を通過する総流束は、V 内部の発散の体積積分と一致する。この事実が Gauss の発散定理である。閉曲面 \partial V は体積領域の境界であり、次のページで境界上の流束と内部の局所量である発散を接続する。特異点を V 内に含む場には、この条件のままでは適用できない。
11よくある誤り
- 法線方向を指定せずに流束を計算する。
- 面積要素 \|\boldsymbol{r}_u\times\boldsymbol{r}_v\| と向き付き面要素 \boldsymbol{r}_u\times\boldsymbol{r}_v を混同する。
- 曲面の境界を確認せず、Stokes 定理の向きを誤る。
12関連リンク
data/lecture/math/vector-calculus/vector-calculus-portal.lecture.n.md
data/lecture/math/vector-calculus/gradient-divergence-and-curl.lecture.n.md
data/lecture/math/multivariable-calculus/parametrized-curves-and-surfaces.lecture.n.md
data/lecture/math/analysis/introduction-to-line-and-surface-integrals.lecture.n.md
data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md
Surface Integrals and Flux
1Introduction
This lecture distinguishes the surface integral of a scalar quantity from the flux of a vector field through a surface. The first accumulates a quantity distributed over area, whereas the second accumulates the component normal to an oriented surface.
2Terminology and Definitions
The flux of a vector field is the surface integral of its normal component and measures the net passage of the field through a surface.
A surface integral accumulates either a scalar quantity over a surface or the normal component of a vector field over that surface.
3Strategy
Represent the surface by \boldsymbol r(u,v) and construct a normal direction from the tangent vectors \boldsymbol r_u and \boldsymbol r_v. In a flux integral, only the component of the vector field in this normal direction contributes.
4Oriented Surfaces
A flux computation must distinguish the two sides of a surface. Reversing the normal direction reverses the sign of the flux. For a closed surface, the outward normal is the standard choice.
A surface is orientable if it admits a continuous unit normal field \boldsymbol n over the entire surface. A connected orientable surface has the two orientations \boldsymbol n and -\boldsymbol n. A globally consistent flux cannot be defined on a nonorientable surface. For a closed surface given as the boundary of a region, the outward orientation is standard.
5Precise Formulation
Let \Omega\subset\mathbb R^2 be a bounded Jordan domain with piecewise C^1 boundary, set D=\overline\Omega, and let \boldsymbol r be C^1 on an open neighborhood of D. If \boldsymbol r_u\times\boldsymbol r_v\ne\boldsymbol0 throughout \Omega and \boldsymbol r is injective on \Omega, then S=\boldsymbol r(D) is a regular surface patch. Overlap restricted to the boundary has area zero and does not affect the integral. If a parameterization covers interior points more than once, the formula over its parameter domain represents an integral with multiplicity rather than the ordinary integral over the geometric surface.
For a continuous scalar field g defined on a neighborhood of S, define
\iint_S g\,dS
=\iint_D g(\boldsymbol r(u,v))\|\boldsymbol r_u\times\boldsymbol r_v\|\,du\,dv.
This integral is independent of orientation. When g=1, it equals the area of the surface.
For a continuous vector field \boldsymbol F defined on a neighborhood of S, suppose that \boldsymbol r_u\times\boldsymbol r_v agrees with the selected orientation. Define the flux by
\iint_S \boldsymbol F\cdot\boldsymbol n\,dS
=\iint_D \boldsymbol F(\boldsymbol r(u,v))\cdot(\boldsymbol r_u\times\boldsymbol r_v)\,du\,dv.
At interior points of the patch, the corresponding unit normal is
\boldsymbol n
=\frac{\boldsymbol r_u\times\boldsymbol r_v}{\|\boldsymbol r_u\times\boldsymbol r_v\|}.
Reversing the orientation reverses the sign of the flux. The vector \boldsymbol r_u\times\boldsymbol r_v simultaneously specifies a direction perpendicular to the tangent plane and the infinitesimal area scale. Thus \|\boldsymbol r_u\times\boldsymbol r_v\|\,du\,dv is the unoriented area element, while the oriented flux uses \boldsymbol r_u\times\boldsymbol r_v\,du\,dv without taking its norm.
If a compact piecewise regular surface of the type considered in this lecture cannot be represented regularly and one-to-one by a single patch, it is divided into finitely many regular patches whose overlaps are confined to boundaries of area zero. Its integral is the sum of the integrals over these patches. For flux, the normal on every patch must agree with the selected orientation of the surface. This construction defines the integrals over the spheres and closed surfaces used below.
Let \Omega' be another bounded Jordan domain with piecewise C^1 boundary and set E=\overline{\Omega'}. Let \boldsymbol\varphi:E\to D be the restriction of a C^1 diffeomorphism between open neighborhoods of E and D, with \boldsymbol\varphi(E)=D, and set \boldsymbol\rho=\boldsymbol r\circ\boldsymbol\varphi. The scalar surface integral is invariant under reparameterization because the change-of-variables formula uses the absolute value of the Jacobian. For flux, \det D\boldsymbol\varphi>0 preserves the orientation and the value. If \det D\boldsymbol\varphi<0 and the resulting cross product is used without correction, the parameterization represents the oppositely oriented surface and the flux changes sign. To preserve the previously selected geometric orientation, one must correct the sign of the cross product.
6Example 1: A Planar Patch
Consider the upward-oriented portion of the plane z=0 with 0\le x\le1 and 0\le y\le1. For \boldsymbol F=(0,0,z+1), the unit normal is (0,0,1), and
\iint_S\boldsymbol F\cdot\boldsymbol n\,dS
=\int_0^1\int_0^1 1\,dx\,dy=1.
Choosing the downward normal changes the value to -1.
7Example 2: A Sphere
Let \boldsymbol F=(x,y,z) and orient the sphere of radius R>0 centered at the origin outward. On the sphere, \boldsymbol F\cdot\boldsymbol n=R. Since its area is 4\pi R^2, the flux is 4\pi R^3. The subsequent lecture verifies the same value by the divergence theorem: the divergence is 3, and the volume is 4\pi R^3/3.
8Example 3: A Cylindrical Side
Parameterize the lateral surface of a cylinder of radius R and height H by
\boldsymbol r(\theta,z)=(R\cos\theta,R\sin\theta,z),
where 0\le\theta\le2\pi and 0\le z\le H. This parameterization excludes the top and bottom disks. The values \theta=0 and \theta=2\pi describe the same seam, but this overlap lies on the boundary of the parameter domain and therefore does not affect the surface integral. The outward-oriented area vector is
\boldsymbol r_\theta\times\boldsymbol r_z
=(R\cos\theta,R\sin\theta,0).
For \boldsymbol F=(x,y,0), the flux is
\int_0^H\int_0^{2\pi}R^2\,d\theta\,dz
=2\pi R^2H.
9Comparison: A Field Tangent to the Surface
On an upward-oriented patch of the plane z=0, the field \boldsymbol F=(1,0,0) is tangent to the plane. Consequently, \boldsymbol F\cdot\boldsymbol n=0, and the flux is zero. A vector field may be nonzero along a surface while contributing no flux if its normal component vanishes.
10Connection to the Divergence Theorem
Let a bounded region V have a piecewise C^1 boundary, and suppose that \boldsymbol F is C^1 on a neighborhood of \overline V. The total outward flux through the closed surface \partial V equals the volume integral of the divergence over V. This is the divergence theorem. The next lecture relates flux through the boundary to the local interior quantity \nabla\cdot\boldsymbol F. The theorem does not apply under these hypotheses if \boldsymbol F has a singularity inside V.
11Common Errors
- Computing flux without specifying a normal orientation.
- Confusing the scalar area factor \|\boldsymbol r_u\times\boldsymbol r_v\| with the oriented area vector \boldsymbol r_u\times\boldsymbol r_v.
- Applying Stokes' theorem without checking the induced orientation of the boundary curve.
12Related Lectures
data/lecture/math/vector-calculus/vector-calculus-portal.lecture.n.md
data/lecture/math/vector-calculus/gradient-divergence-and-curl.lecture.n.md
data/lecture/math/multivariable-calculus/parametrized-curves-and-surfaces.lecture.n.md
data/lecture/math/analysis/introduction-to-line-and-surface-integrals.lecture.n.md
data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md