Teoremas integrales de campo

Teoremas de Gauss-Divergencia, Stokes y Green que conectan integrales diferenciales y de contorno.

Que es

Los teoremas integrales del cálculo vectorial establecen puentes fundamentales entre el comportamiento diferencial local de un campo en el interior de una región geométrica y la acumulación integral en su frontera orientada:

  1. Teorema de Green en el plano: ∮∂D(P dx+Q dy)=∬D(∂Q∂x−∂P∂y)dx dy\oint_{\partial D} (P \, dx + Q \, dy) = \iint_D \left( \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} \right) dx \, dy
  2. Teorema de Stokes (Teorema del Rotacional): ∮∂SF⃗⋅dr⃗=∬S(∇⃗×F⃗)⋅dS⃗\oint_{\partial S} \vec{F} \cdot d\vec{r} = \iint_S (\vec{\nabla} \times \vec{F}) \cdot d\vec{S} La circulación de un campo vectorial a lo largo de una curva cerrada orientada ∂S\partial S es igual al flujo de su rotacional a través de cualquier superficie orientable SS que tenga a dicha curva como borde.
  3. Teorema de Gauss (Teorema de la Divergencia): ∮∂VF⃗⋅dS⃗=∭V(∇⃗⋅F⃗) dV\oint_{\partial V} \vec{F} \cdot d\vec{S} = \iiint_V (\vec{\nabla} \cdot \vec{F}) \, dV El flujo neto saliente a través de una superficie cerrada ∂V\partial V es igual a la integral de volumen de su divergencia en el interior de VV.

Todos ellos son manifestaciones particulares del Teorema General de Stokes para formas diferenciales: ∫∂Ωω=∫Ωdω\int_{\partial \Omega} \omega = \int_\Omega d\omega

Por qué aparece en física

Permite transformar las leyes físicas entre su forma diferencial local (válida en cada punto del espacio) y su forma integral macroscópica (medible en el laboratorio):

  • Ecuaciones de Maxwell: El teorema de Gauss transforma ∇⃗⋅E⃗=ρ/ε0\vec{\nabla}\cdot\vec{E} = \rho/\varepsilon_0 en ∮E⃗⋅dS⃗=Qenc/ε0\oint \vec{E}\cdot d\vec{S} = Q_{\text{enc}}/\varepsilon_0; el teorema de Stokes transforma ∇⃗×E⃗=−∂tB⃗\vec{\nabla}\times\vec{E} = -\partial_t \vec{B} en ∮E⃗⋅dr⃗=−dΦBdt\oint \vec{E}\cdot d\vec{r} = -\frac{d\Phi_B}{dt}.
  • Conservación de masa en fluidos: Transforma la ecuación diferencial de continuidad en el balance macroscópico de caudales.
  • Teorema de Helmholtz: Cualquier campo vectorial físico suave que decae en el infinito se descompone de forma única en un componente irrotacional (gradiente de escalar) y un componente solenoidal (rotacional de vector): F⃗=−∇⃗ϕ+∇⃗×A⃗\vec{F} = -\vec{\nabla}\phi + \vec{\nabla}\times\vec{A}.

Idea visual o geométrica

  • En el Teorema de Gauss: Imagina el volumen como un panal con millones de celdas microscópicas. En cada celda, la divergencia mide cuánto fluido se crea. En las paredes internas compartidas, el flujo que sale de una celda entra en la vecina y se anula mutuamente. Solo sobreviven las paredes exteriores de la superficie frontera exterior ∂V\partial V.
  • En el Teorema de Stokes: Imagina una red de pequeños lazos giratorios contiguos en una superficie. Los bordes internos circulan en sentidos opuestos y se cancelan mutuamente; únicamente sobrevive la circulación en el borde perimetral exterior ∂S\partial S.

Lenguaje simbólico mínimo

  • Teorema de Green: ∮∂D(Pdx+Qdy)=∬D(∂xQ−∂yP)dA\oint_{\partial D} (P dx + Q dy) = \iint_D (\partial_x Q - \partial_y P) dA.
  • Teorema de Stokes: ∮CF⃗⋅dr⃗=∬S(∇⃗×F⃗)⋅dS⃗\oint_C \vec{F}\cdot d\vec{r} = \iint_S (\vec{\nabla}\times\vec{F})\cdot d\vec{S}.
  • Teorema de Gauss: ∮SF⃗⋅dS⃗=∭V(∇⃗⋅F⃗)dV\oint_S \vec{F}\cdot d\vec{S} = \iiint_V (\vec{\nabla}\cdot\vec{F}) dV.
  • Identidades de Green:
    1. ∫V(ϕ∇2ψ+∇⃗ϕ⋅∇⃗ψ)dV=∮S(ϕ∇⃗ψ)⋅dS⃗\int_V (\phi \nabla^2 \psi + \vec{\nabla}\phi \cdot \vec{\nabla}\psi) dV = \oint_S (\phi \vec{\nabla}\psi) \cdot d\vec{S}.
    2. ∫V(ϕ∇2ψ−ψ∇2ϕ)dV=∮S(ϕ∇⃗ψ−ψ∇⃗ϕ)⋅dS⃗\int_V (\phi \nabla^2 \psi - \psi \nabla^2 \phi) dV = \oint_S (\phi \vec{\nabla}\psi - \psi \vec{\nabla}\phi) \cdot d\vec{S}.

Que significa el resultado

Permiten sustituir integrales sobre geometrías complicadas por integrales sobre geometrías triviales con la misma frontera, simplificando drásticamente los cálculos analíticos y numéricos.

Cuando no debes usarlo

  • No uses el Teorema de Stokes si la curva no es cerrada o si la superficie contiene singularidades (como cortes o agujeros donde el campo no es diferenciable).
  • No uses el Teorema de Gauss en una superficie abierta (como un cono sin base plana).
  • En el Teorema de Green, la curva frontera ∂D\partial D debe recorrerse en sentido antihorario (dejando la región a la izquierda) para que los signos coincidan.

Resumen operativo

Para resolver una integral de línea difícil en una curva cerrada plana: aplica el Teorema de Green. Para una integral de línea cerrada en el espacio 3D: aplica Stokes integrando el rotacional en la superficie más sencilla que tenga a la curva como borde. Para una integral de flujo a través de una superficie cerrada: aplica el Teorema de Gauss integrando la divergencia en el volumen interior.

Procedimiento

Antes de calcular

Identifica si el problema involucra una curva cerrada ∮C\oint_C o una superficie cerrada ∮S\oint_S. Verifica si el campo y sus derivadas son continuos en todo el dominio interior (sin singularidades como 1/r1/r en el origen).

Paso 1: identifica el objeto matemático

Clasifica la dimensión geométrica:

  • Si es una integral de línea cerrada en 2D: Teorema de Green.
  • Si es una integral de línea cerrada en 3D: Teorema de Stokes.
  • Si es una integral de flujo sobre una superficie cerrada en 3D: Teorema de Gauss.

Paso 2: elige la operación

Calcula la derivada diferencial del integrando:

  • Para Green: calcula (∂xQ−∂yP)(\partial_x Q - \partial_y P).
  • Para Stokes: calcula ∇⃗×F⃗\vec{\nabla}\times\vec{F}.
  • Para Gauss: calcula ∇⃗⋅F⃗\vec{\nabla}\cdot\vec{F}.

Paso 3: conserva unidades y signos

Determina los límites de la región interior:

  • En Green/Stokes: confirma que la normal a la superficie cumpla la regla de la mano derecha con la circulación.
  • En Gauss: confirma que la normal apunte hacia el exterior del volumen.

Paso 4: calcula

Integra la función diferencial sobre la región interior (integral doble para Green/Stokes, integral triple para Gauss). Aprovecha simetrías polares, cilíndricas o esféricas para simplificar el cálculo.

Paso 5: interpreta el resultado

Verifica la equivalencia física: el resultado numérico obtenido debe coincidir exactamente con el que se habría obtenido integrando laboriosamente cada tramo de la frontera.

Checklist de autocontrol

  • ¿La curva o superficie es verdaderamente cerrada?
  • ¿Se comprobó la ausencia de singularidades en el interior del recinto?
  • ¿La orientación de la frontera respeta la convención estándar (mano derecha / normal saliente)?

Ejemplo guiado

Situación física

Se desea calcular el flujo saliente total del campo de fuerzas: F⃗(x,y,z)=(x3+sin⁡(yz))i^+(y3+exz)j^+(z3+ln⁡(1+x2y2))k^\vec{F}(x,y,z) = (x^3 + \sin(yz)) \hat{i} + (y^3 + e^{xz}) \hat{j} + (z^3 + \ln(1+x^2 y^2)) \hat{k} a través de la superficie esférica cerrada SS de radio R=2.0 mR = 2.0 \text{ m} centrada en el origen (x2+y2+z2=4x^2 + y^2 + z^2 = 4). Calcular esta integral de superficie directamente sería intratable debido a los términos trascendentes sin⁡(yz),exz,ln⁡(1+x2y2)\sin(yz), e^{xz}, \ln(1+x^2 y^2).

Traducción matemática

Aplicamos el Teorema de la Divergencia de Gauss: ∮SF⃗⋅dS⃗=∭V(∇⃗⋅F⃗) dV\oint_S \vec{F} \cdot d\vec{S} = \iiint_V (\vec{\nabla} \cdot \vec{F}) \, dV donde VV es la bola sólida x2+y2+z2≤4x^2 + y^2 + z^2 \le 4.

Calculo paso a paso

  1. Cálculo de la divergencia de F⃗\vec{F}: ∂Fx∂x=∂∂x(x3+sin⁡(yz))=3x2\frac{\partial F_x}{\partial x} = \frac{\partial}{\partial x}(x^3 + \sin(yz)) = 3x^2 ∂Fy∂y=∂∂y(y3+exz)=3y2\frac{\partial F_y}{\partial y} = \frac{\partial}{\partial y}(y^3 + e^{xz}) = 3y^2 ∂Fz∂z=∂∂z(z3+ln⁡(1+x2y2))=3z2\frac{\partial F_z}{\partial z} = \frac{\partial}{\partial z}(z^3 + \ln(1+x^2 y^2)) = 3z^2 ∇⃗⋅F⃗=3x2+3y2+3z2=3(x2+y2+z2)\vec{\nabla} \cdot \vec{F} = 3x^2 + 3y^2 + 3z^2 = 3(x^2 + y^2 + z^2) ¡Todos los términos trascendentes complejos se han anulado idénticamente en las derivadas parciales!

  2. Planteamiento de la integral de volumen en coordenadas esféricas: x2+y2+z2=r2,dV=r2sin⁡θ dr dθ dφx^2 + y^2 + z^2 = r^2, \quad dV = r^2 \sin\theta \, dr \, d\theta \, d\varphi Límites: r∈[0,2]r \in [0, 2], θ∈[0,π]\theta \in [0, \pi], φ∈[0,2π]\varphi \in [0, 2\pi].

  3. Evaluación de la integral triple: ∭V3r2(r2sin⁡θ dr dθ dφ)=3∫02πdφ∫0πsin⁡θ dθ∫02r4 dr\iiint_V 3r^2 (r^2 \sin\theta \, dr \, d\theta \, d\varphi) = 3 \int_0^{2\pi} d\varphi \int_0^\pi \sin\theta \, d\theta \int_0^2 r^4 \, dr

    • Integral angular: ∫02πdφ∫0πsin⁡θ dθ=(2π)(2)=4π\int_0^{2\pi} d\varphi \int_0^\pi \sin\theta \, d\theta = (2\pi)(2) = 4\pi.
    • Integral radial: ∫02r4 dr=[r55]02=325=6.4\int_0^2 r^4 \, dr = \left[ \frac{r^5}{5} \right]_0^2 = \frac{32}{5} = 6.4.
  4. Flujo total: Φ=3×(4π)×325=384π5=76.8π≈241.274[unidades de flujo]\Phi = 3 \times (4\pi) \times \frac{32}{5} = \frac{384\pi}{5} = 76.8\pi \approx 241.274 \quad [\text{unidades de flujo}]

Lectura del resultado

El teorema de Gauss transformó una integral de superficie no resoluble analíticamente en una integral de volumen polinómica elemental. El flujo saliente positivo confirma que el campo se expande desde el interior del volumen esférico.

Error común evitado

Intentar parametrizar la esfera en ángulos (θ,φ)(\theta, \varphi) e integrar directamente las funciones trigonométricas complejas. Siempre debe calcularse primero la divergencia para comprobar si el integrando volumétrico colapsa a una expresión sencilla.


Ejemplo inverso

Resultado dado

La circulación de un campo de velocidades alrededor de una circunferencia unitaria en el plano xyxy es ∮Cv⃗⋅dr⃗=8π\oint_C \vec{v} \cdot d\vec{r} = 8\pi.

Que pregunta física permite responder

¿Cuál es la vorticidad media ⟨ω⃗⟩\langle \vec{\omega} \rangle en el disco encerrado? Por el teorema de Stokes: ∮Cv⃗⋅dr⃗=∬S(∇⃗×v⃗)⋅k^ dS=⟨ωz⟩×Aˊrea\oint_C \vec{v} \cdot d\vec{r} = \iint_S (\vec{\nabla} \times \vec{v}) \cdot \hat{k} \, dS = \langle \omega_z \rangle \times \text{Área} El área del disco unitario es π(1)2=π\pi (1)^2 = \pi. ⟨ωz⟩=8ππ=8 rad/s\langle \omega_z \rangle = \frac{8\pi}{\pi} = 8 \text{ rad/s} El fluido en el disco gira con una vorticidad media neta de 8 s−18 \text{ s}^{-1}.

Errores frecuentes

Error 1: Aplicar el Teorema de Gauss a una superficie abierta

Intentar calcular el flujo a través de un tazón hemisférico abierto usando la integral de volumen de la divergencia sin cerrar la superficie con la tapa circular.

Error 2: Incompatibilidad de orientación en Stokes

No seguir la regla de la mano derecha entre el sentido de recorrido de la curva frontera ∂S\partial S y el vector normal n^\hat{n} de la superficie. Produce un error de signo global.

Error 3: Ignorar singularidades en el interior del volumen

Aplicar el Teorema de Gauss al campo de Coulomb E⃗=q4πε0r2r^\vec{E} = \frac{q}{4\pi\varepsilon_0 r^2}\hat{r} asumiendo que como ∇⃗⋅E⃗=0\vec{\nabla}\cdot\vec{E} = 0 para r>0r>0, el flujo total es cero. El origen contiene una singularidad tipo delta de Dirac que encierra toda la carga.

Error 4: Confundir el Teorema de Green con el Teorema de Stokes

Olvidar que Green es simplemente el caso particular bidimensional plano del Teorema de Stokes general.

Error 5: Integrar en sentido horario en el Teorema de Green

Recorrer la curva en sentido horario, lo que invierte el signo del resultado: ∮Chorario=−∬D(∂xQ−∂yP)dA\oint_{C_{\text{horario}}} = -\iint_D (\partial_x Q - \partial_y P) dA.

Regla rápida de detección

Si al aplicar Gauss o Stokes sobre un campo libre de singularidades el resultado difiere en un signo menos del cálculo por definición, la orientación de la normal o de la frontera fue invertida.

Aplicaciones físicas

1. Unificación de las Ecuaciones de Maxwell

Los teoremas de Gauss y Stokes permiten deducir la forma integral de las cuatro leyes de Maxwell a partir de su formulación diferencial local, permitiendo diseñar antenas, radares y cavidades resonantes.

2. Conservación global de energía y fluidos

En mecánica de fluidos computacional (CFD) y meteorología, los teoremas integrales transforman balances microscópicos de conservación de masa y momento en balances de flujos en las fronteras computacionales.

3. Descomposición de Helmholtz en campos físicos

Todo campo físico se descompone en un modo puramente longitudinal sin rotacional (−∇⃗ϕ-\vec{\nabla}\phi, como la fuerza electrostática) y un modo puramente transversal solenoidal (∇⃗×A⃗\vec{\nabla}\times\vec{A}, como el campo magnético).

4. Astrofísica estelar y teorema virial

El teorema de Gauss permite calcular la energía de auto-gravitación de estrellas y galaxias integrando la divergencia del tensor de tensiones gravitatorio sobre volúmenes cósmicos.

5. Dinámica de vórtices y tornados

El teorema de Stokes establece que la circulación de la velocidad alrededor de un vórtice atmosférico o remolino hidrodinámico es igual al flujo total de vorticidad que atraviesa su núcleo.

Practica

Problema 1: Teorema de Stokes en un lazo triangular

Sea el campo F⃗(x,y,z)=(z)i^+(x)j^+(y)k^\vec{F}(x,y,z) = (z) \hat{i} + (x) \hat{j} + (y) \hat{k}. Calcule la circulación ∮CF⃗⋅dr⃗\oint_C \vec{F} \cdot d\vec{r} a lo largo del triángulo con vértices en (1,0,0)(1,0,0), (0,1,0)(0,1,0) y (0,0,1)(0,0,1) recorrido en ese orden, aplicando el Teorema de Stokes.

Solución:

  1. Rotacional de F⃗\vec{F}: ∇⃗×F⃗=∣i^j^k^∂x∂y∂zzxy∣=(1−0)i^+(1−0)j^+(1−0)k^=i^+j^+k^\vec{\nabla}\times\vec{F} = \begin{vmatrix} \hat{i} & \hat{j} & \hat{k} \\ \partial_x & \partial_y & \partial_z \\ z & x & y \end{vmatrix} = (1 - 0)\hat{i} + (1 - 0)\hat{j} + (1 - 0)\hat{k} = \hat{i} + \hat{j} + \hat{k}
  2. Superficie plana triangular SS: La ecuación del plano que pasa por los tres puntos es x+y+z=1x + y + z = 1. El vector normal unitario hacia el primer octante es n^=i^+j^+k^3\hat{n} = \frac{\hat{i} + \hat{j} + \hat{k}}{\sqrt{3}}. Vector dS⃗=n^dS=(i^+j^+k^)dxdyd\vec{S} = \hat{n} dS = (\hat{i} + \hat{j} + \hat{k}) dx dy proyectado sobre el plano xyxy.
  3. Flujo del rotacional: (∇⃗×F⃗)⋅n^=(i^+j^+k^)⋅i^+j^+k^3=33=3(\vec{\nabla}\times\vec{F}) \cdot \hat{n} = (\hat{i}+\hat{j}+\hat{k}) \cdot \frac{\hat{i}+\hat{j}+\hat{k}}{\sqrt{3}} = \frac{3}{\sqrt{3}} = \sqrt{3} Área del triángulo equilátero con lados de longitud 2\sqrt{2}: Aˊrea=34(2)2=32\text{Área} = \frac{\sqrt{3}}{4}(\sqrt{2})^2 = \frac{\sqrt{3}}{2} ∮CF⃗⋅dr⃗=∬S(∇⃗×F⃗)⋅dS⃗=3(32)=32=1.50\oint_C \vec{F}\cdot d\vec{r} = \iint_S (\vec{\nabla}\times\vec{F})\cdot d\vec{S} = \sqrt{3} \left( \frac{\sqrt{3}}{2} \right) = \frac{3}{2} = 1.50

Problema 2: Teorema de Green para el cálculo de áreas

El Teorema de Green permite calcular el área de cualquier región plana mediante la integral de línea Aˊrea(D)=12∮∂D(xdy−ydx)\text{Área}(D) = \frac{1}{2} \oint_{\partial D} (x dy - y dx). Calcule el área encerrada por la elipse x2a2+y2b2=1\frac{x^2}{a^2} + \frac{y^2}{b^2} = 1 parametrizada como x(t)=acos⁡(t)x(t) = a \cos(t), y(t)=bsin⁡(t)y(t) = b \sin(t) para t∈[0,2π]t \in [0, 2\pi].

Solución:

  1. Diferenciales: dx=−asin⁡(t)dt,dy=bcos⁡(t)dtdx = -a \sin(t) dt, \quad dy = b \cos(t) dt
  2. Forma diferencial de Green: xdy−ydx=[acos⁡(t)](bcos⁡(t)dt)−[bsin⁡(t)](−asin⁡(t)dt)=ab(cos⁡2t+sin⁡2t)dt=ab dtx dy - y dx = [a \cos(t)](b \cos(t) dt) - [b \sin(t)](-a \sin(t) dt) = a b (\cos^2 t + \sin^2 t) dt = a b \, dt
  3. Integración: Aˊrea=12∮∂D(xdy−ydx)=12∫02πab dt=12ab(2π)=πab\text{Área} = \frac{1}{2} \oint_{\partial D} (x dy - y dx) = \frac{1}{2} \int_0^{2\pi} a b \, dt = \frac{1}{2} a b (2\pi) = \pi a b

Problema 3: Teorema de Gauss en un cilindro

Calcule el flujo del campo F⃗=(xy2)i^+(yx2)j^+(z)k^\vec{F} = (x y^2) \hat{i} + (y x^2) \hat{j} + (z) \hat{k} a través de la superficie cilíndrica cerrada x2+y2≤4x^2 + y^2 \le 4, 0≤z≤30 \le z \le 3.

Solución:

  1. Divergencia: ∇⃗⋅F⃗=∂(xy2)∂x+∂(yx2)∂y+∂z∂z=y2+x2+1\vec{\nabla}\cdot\vec{F} = \frac{\partial(x y^2)}{\partial x} + \frac{\partial(y x^2)}{\partial y} + \frac{\partial z}{\partial z} = y^2 + x^2 + 1
  2. Por Gauss, en coordenadas cilíndricas (x2+y2=r2,dV=rdrdθdzx^2+y^2=r^2, dV = r dr d\theta dz): Φ=∫03dz∫02πdθ∫02(r2+1)r dr=(3)(2π)[r44+r22]02=6π(164+42)=6π(4+2)=36π≈113.10\Phi = \int_0^3 dz \int_0^{2\pi} d\theta \int_0^2 (r^2 + 1) r \, dr = (3)(2\pi) \left[ \frac{r^4}{4} + \frac{r^2}{2} \right]_0^2 = 6\pi \left( \frac{16}{4} + \frac{4}{2} \right) = 6\pi (4 + 2) = 36\pi \approx 113.10