Càlcul de volums a partir del teorema de Gauss

De la mateixa manera que es fa amb el teorema de Green, utilitzarem una potent eina del càlcul integral per calcular volums, el teorema de la divergència o teorema de Gauss.

Aquest teorema ens diu que si tenim un sòlid D limitat per una superfície tancada S i F(x,y,z) un camp vectorial, llavors:

SF dS=Vdiv(F) dx dy dz

on, si F(x,y,z)=(P(x,y,z),Q(x,y,z),R(x,y,z)), aleshores div(F)=Px+Qy+Rz és la divergència del camp.

Recordem que per calcular el primer terme, si r(u,v) és una parametrització de la superfície,

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

on ru×rv denota el producte vectorial de les derivades de la parametrització segons u i v.

Per tant, si trobem un camp que tingui divergència idènticament igual a 1, en el membre de la dreta tindrem el volum del sòlid i per tant tenim un mètode per calcular volums.

De camps que compleixin el que necessitem n'hi ha molts, però els més típics són:

  • 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)

Exemple

Per exemple, anem a calcular el volum delimitat per mitja esfera i el pla equatorial, és a dir:

imagen

Pel teorema de Gauss tenim que podem integrar un dels camps donats al llarg de la superfície que tanca el volum.

Parametritzem primer el casquet amb coordenades esfèriques:

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

Calculem els vectors derivades i el seu producte 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θ)

Prenent el camp vectorial F(x,y,z)=(0,0,z), tenim:

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π

Ara hem d'integrar el mateix camp sobre la tapa inferior ja que el teorema és vàlid per superfícies tancades. Una parametrització del cercle és:

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

Calculem també els vectors derivades i el seu producte 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)

Ara, tenim:

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

Per tant,

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π