Cálculo de volúmenes a partir del teorema de Gauss

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 F(x,y,z) un campo vectorial, entonces:

SF dS=Vdiv(F) dx dy dz

donde, si F(x,y,z)=(P(x,y,z),Q(x,y,z),R(x,y,z)), entonces div(F)=Px+Qy+Rz es la divergencia del campo.

Recordamos que para calcular el primer término, si r(u,v) ées una parametrización de la superficie,

SF dS=F(r(u,v))(ru×rv) du dv

donde ru×rv 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:

  • F(x,y,z)=(x,0,0)
  • F(x,y,z)=(0,y,0)
  • F(x,y,z)=(0,0,z)
  • F(x,y,z)=13(x,y,z)

Ejemplo

Vamos a calcular el volumen delimitado por media esfera y el plano ecuatorial, es decir:

imagen

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:

x=sinθcosφy=sinθsinφz=cosθθ[0,π2]φ[0,2π]

Calculemos los vectores derivadas y su producto vectorial:

rθ=(cosθcosφ,sinφcosθ,sinφ)rφ=(sinθsinφ,sinθcosφ,0)}

 rθ×rφ=|ijkcosθcosφsinφcosθsinφsinθsinφsinθcosφ0|=(sin2θcosφ,sin2θsinφ,sinθcosθ)

Tomando el campo vecotorial F(x,y,z)=(0,0,z), tenemos:

S(0,0,z) dS=0π202π(0,0,cosθ)(sin2θcosφ,sin2θsinφ,sinθcosθ) dφ dθ=0π202πcos2θsinθ dφ dθ=0π2cos2θsinθ(02πdφ) dθ=2π0π2cos2θsinθ dθ=2π[cos3θ3]0π2=23π

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:

x=rcos(t)y=rsin(t)z=0t[0,2π]r[0,1]

Calculemos también los vectores derivadas y su producto vectorial:

rr=(rsin(t),rcos(t),0)rt=(cos(t),sin(t),0)}

 rt×rr=|ijkrsin(t)rcos(t)0cos(t)sin(t)0|=(0,0,r2)

Ahora, tenemos:

C(0,0,z) dS=02π01(0,0,0)(0,0,r2) dr dt=0

Por lo tanto,

Vol(E)=S1 dx dy dz=Ediv(0,0,z) dx dy dz=S(0,0,z) dS+C(0,0,z) dS=23π+0=23π