De la misma forma que se hace con el teorema de Green, utilizaremos una potente herramienta del cálculo integral para calcular volúmenes, el teorema de la divergencia o teorema de Gauss.
Este teorema nos dice que si tenemos un sólido $$D$$ limitado por una superficie cerrada $$S$$ y $$\vec{F}(x,y,z)$$ un campo vectorial, entonces:
$$$\int_S \vec{F} \ d\vec{S} = \iiint_V \text{div}(\vec{F}) \ dx \ dy \ dz $$$
donde, si $$\vec{F}(x,y,z)=(P(x,y,z),Q(x,y,z),R(x,y,z))$$, entonces $$\text{div}(\vec{F})=P_x+Q_y+R_z$$ es la divergencia del campo.
Recordamos que para calcular el primer término, si $$r(u, v)$$ ées una parametrización de la superficie,
$$$\int_S \vec{F} \ d\vec{S}=\iint \vec{F}\big( r(u,v) \big) \cdot (r_u\times r_v) \ du \ dv$$$
donde $$r_u\times r_v$$ denota el producto vectorial de las derivadas de la parametrización según $$u$$ y $$v$$.
Por lo tanto, si encontramos un campo que tenga divergencia idénticamente igual a $$1$$, en el miembro de la derecha tendremos el volumen del sólido y por lo tanto tenemos un método para calcular volúmenes.
De campos que cumplan lo que necesitamos hay muchos, pero los más típicos son:
- $$\vec{F}(x,y,z)=(x,0,0)$$
- $$\vec{F}(x,y,z)=(0,y,0)$$
- $$\vec{F}(x,y,z)=(0,0,z)$$
- $$\vec{F}(x,y,z)=\dfrac{1}{3}(x,y,z)$$
Vamos a calcular el volumen delimitado por media esfera y el plano ecuatorial, es decir:
Por el teorema de Gauss tenemos que podemos integrar uno de los campos dado a lo largo de la superficie que encierra el volumen.
Parametricemos primero el casquete con coordenadas esféricas:
$$$ \begin{array}{l} x=\sin\theta\cos\varphi \\ y=\sin\theta\sin\varphi \\ z=\cos\theta \end{array} \qquad \begin{array}{l} \theta\in \big[0,\dfrac{\pi}{2}\big] \\ \varphi\in[0,2\pi] \end{array}$$$
Calculemos los vectores derivadas y su producto vectorial:
$$$ \left. \begin{array}{l} r_\theta = (\cos\theta\cos\varphi,\sin\varphi\cos\theta,-\sin\varphi) \\ r_\varphi=(-\sin\theta\sin\varphi,\sin\theta\cos\varphi,0) \end{array} \right\} \Rightarrow$$$
$$$\begin{array}{rl} \Rightarrow \ r_\theta \times r_\varphi =& \begin{vmatrix} i & j & k \\ \cos\theta\cos\varphi & \sin\varphi\cos\theta & -\sin\varphi\\ -\sin\theta\sin\varphi & \sin\theta\cos\varphi & 0 \end{vmatrix} \\ =& \big( \sin^2\theta\cos\varphi, \sin^2\theta\sin\varphi, \sin\theta\cos\theta\big) \end{array}$$$
Tomando el campo vecotorial $$F(x,y,z)=(0,0,z)$$, tenemos:
$$$ \begin{array}{rl} \int_S (0,0,z) \ d\vec{S} =& \int_0^{\frac{\pi}{2}} \int_0^{2\pi} (0,0,\cos\theta)\cdot \big( \sin^2\theta\cos\varphi, \sin^2\theta\sin\varphi, \sin\theta\cos\theta\big) \ d\varphi \ d\theta \\ =& \int_0^{\frac{\pi}{2}} \int_0^{2\pi} \cos^2\theta \cdot \sin\theta \ d\varphi \ d\theta = \int_0^{\frac{\pi}{2}} \cos^2\theta \cdot \sin\theta \Big( \int_0^{2\pi} d\varphi \Big) \ d\theta \\ =& 2\pi \int_0^{\frac{\pi}{2}} \cos^2\theta \cdot \sin\theta \ d\theta = 2\pi \Big[ \dfrac{-\cos^3\theta}{3} \Big]_0^{\frac{\pi}{2}}= \dfrac{2}{3}\pi \end{array}$$$
Ahora debemos integrar el mismo campo sobre el la tapa inferior ya que el teorema es válido por superficies cerradas. Una parametrización del círculo es:
$$$ \begin{array}{l} x=r\cos(t) \\ y=r\sin(t) \\ z=0 \end{array} \qquad \begin{array}{l} t\in[0,2\pi] \\ r\in[0,1] \end{array}$$$
Calculemos también los vectores derivadas y su producto vectorial:
$$$ \left. \begin{array}{l} r_r = (-r\sin(t),r\cos(t),0) \\ r_t=(\cos(t),\sin(t),0) \end{array} \right\} \Rightarrow$$$
$$$\Rightarrow \ r_t \times r_r = \begin{vmatrix} i & j & k \\ -r\sin(t) & r\cos(t) & 0 \\ \cos(t) & \sin(t) & 0 \end{vmatrix} = ( 0, 0, -r^2)$$$
Ahora, tenemos:
$$$ \int_C (0,0,z) \ d\vec{S} = \int_0^{2\pi} \int_0^1 (0,0,0)\cdot (0,0,-r^2) \ dr \ dt =0 $$$
Por lo tanto,
$$$ \begin{array}{rl} \text{Vol}(E)=&\iiint_S 1 \ dx\ dy\ dz= \iiint_E \text{div}(0,0,z) \ dx\ dy\ dz \\ =& \int_S (0,0,z) \ d\vec{S}+ \int_C (0,0,z)\ d\vec{S} = \dfrac{2}{3}\pi +0= \dfrac{2}{3}\pi \end{array}$$$