Teorema de Green: Integrales de Línea, Rotacional en el Plano y Ejercicios Resueltos

¡Qué tal amigos! El día de hoy exploraremos una de las herramientas más elegantes, prácticas y fascinantes de todo el cálculo multivariable y la física matemática: el Teorema de Green. Nombrado en honor al brillante matemático y físico británico George Green, este teorema establece un puente mágico e indestructible entre dos mundos que parecían distintos: una integral de línea cerrada a lo largo de la frontera de una región y una integral doble evaluada sobre el área interior de dicha región. ¿Te imaginas tener que parametrizar cuatro curvas complicadas con funciones trigonométricas y sustituciones engorrosas para calcular el trabajo de un campo de fuerzas, y en su lugar poder resolver una integral doble elemental en apenas tres líneas de cálculo? Eso es exactamente lo que el Teorema de Green nos permite lograr. En esta guía completa aprenderás sus hipótesis rigurosas, orientación positiva, interpretación física mediante vórtices y rotacional microscópico, fórmulas de área para curvas cerradas y 9 ejercicios resueltos paso a paso con todo detalle.

¿Qué es el Teorema de Green y cuál es su enunciado formal?

El Teorema de Green es la versión bidimensional del Teorema de Stokes (y del Teorema Fundamental del Cálculo). Transforma la circulación total de un campo vectorial alrededor de una trayectoria cerrada simple en la suma continua de todas las tendencias microscópicas a rotar (el rotacional escalar) en cada punto infinitesimal del interior de la región encerrada.

Orientación Positiva de la Curva Frontera

Una curva cerrada simple \( C \) tiene orientación positiva si se recorre en sentido antihorario, de modo que al caminar a lo largo de la trayectoria en la dirección de la flecha directriz, el interior de la región \( D \) siempre permanece a nuestra mano izquierda.

Si la curva se recorre en sentido horario (orientación negativa \( -C \)), el valor de la integral de línea simplemente cambia de signo:

\[ \oint_{-C} (L\,dx + M\,dy) = - \oint_C (L\,dx + M\,dy) = - \iint_D \left( \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} \right) dA \]

Interpretación Física: Rotacional y Cancelación de Vórtices Microscópicos

Para comprender intuitivamente por qué el Teorema de Green es verdadero, imaginemos subdividir la región \( D \) en una cuadrícula fina de miles de rectángulos microscópicos. En cada micro-rectángulo, la circulación del campo a lo largo de sus 4 lados representa el rotacional bidimensional o vorticidad neta:

Al sumar la circulación de todos los rectángulos adyacentes, las aristas compartidas interiores se recorren exactamente dos veces en sentidos opuestos, cancelándose perfectamente entre sí. Las únicas aristas que no se anulan son las que se encuentran en el borde exterior, cuya suma da exactamente la integral de línea a lo largo de la frontera global \( C \).

Cálculo de Áreas Planas con Integrales de Línea

Una de las aplicaciones más espectaculares del Teorema de Green es el cálculo de áreas planas encerradas por curvas complejas (como elipses, astroides o lazos de hipocicloides). Si elegimos las funciones \( L \) y \( M \) de modo que \( \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = 1 \), la integral doble se reduce directamente al área de la región \( D \):

La tercera fórmula, \( \frac{1}{2} \oint_C (x\,dy - y\,dx) \), es especialmente popular y simétrica, permitiendo calcular áreas exactas mediante una parametrización simple de la curva frontera.

Interpretación Física: Rotacional y Vorticidad Microscópica

En cada punto interior, el rotacional escalar representa la velocidad angular o tendencia de giro infinitesimal del fluido:

Forma Normal: Teorema de la Divergencia Bidimensional

Si proyectamos el campo vectorial sobre el vector normal exterior \( \mathbf{n} \), obtenemos la equivalencia de flujo y divergencia:

Regiones Múltiplemente Conexas (Regiones con Agujeros)

El Teorema de Green también es válido para regiones con una o más cavidades (como un anillo o corona circular), siempre que la frontera exterior \( C_1 \) se recorra en sentido antihorario y las fronteras interiores \( C_2 \) en sentido horario, asegurando que la región siempre quede a la izquierda.

Campos Conservativos y Circulación Nula

Cuando el rotacional se anula en toda la región conexa, la integral de línea cerrada es idénticamente cero, lo que garantiza que el campo es conservativo y el trabajo es independiente de la trayectoria.

9 Ejercicios Resueltos Paso a Paso del Teorema de Green

Solución:

Datos:

• Función \( L(x, y) = 3y - e^x \)
• Función \( M(x, y) = 7x + \sqrt{y^4 + 1} \)
• Vértices del triángulo \( D \): \( (0, 0) \), \( (2, 0) \), \( (2, 4) \)
• Orientación: Positiva (antihorario)

Fórmulas:

Teorema de Green:

\[ \oint_C (L\,dx + M\,dy) = \iint_D \left( \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} \right) dA \]

Sustitución y Operaciones:

Calculamos las derivadas parciales:

\[ \frac{\partial M}{\partial x} = \frac{\partial}{\partial x}(7x + \sqrt{y^4 + 1}) = 7 \]

\[ \frac{\partial L}{\partial y} = \frac{\partial}{\partial y}(3y - e^x) = 3 \]

El integrando es una constante:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = 7 - 3 = 4 \]

Sustituimos en la integral doble:

\[ \iint_D 4 \, dA = 4 \iint_D dA = 4 \cdot \text{Área}(D) \]

El triángulo tiene base \( b = 2 \) y altura \( h = 4 \):

\[ \text{Área}(D) = \frac{1}{2}(2)(4) = 4 \]

Multiplicamos por la constante:

\[ 4 \cdot 4 = 16 \]

Resultado:

\[ \oint_C (L\,dx + M\,dy) = 16 \]

Interpretación Matemática: La integral de línea cerrada equivale a 4 veces el área geométrica encerrada por la poligonal triangular.

Solución:

Datos:

• \( L(x, y) = y^2 \)
• \( M(x, y) = 3xy \)
• Región \( D \): Disco unitario \( x^2 + y^2 \le 1 \)

Fórmulas:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = \frac{\partial}{\partial x}(3xy) - \frac{\partial}{\partial y}(y^2) = 3y - 2y = y \]

\[ \oint_C (L\,dx + M\,dy) = \iint_D y \, dA \]

Sustitución y Operaciones:

Pasamos a coordenadas polares: \( x = r\cos\theta \), \( y = r\sin\theta \), \( dA = r\,dr\,d\theta \), con \( r \in [0, 1] \) y \( \theta \in [0, 2\pi] \):

\[ \iint_D y \, dA = \int_0^{2\pi} \int_0^1 (r\sin\theta) r \, dr \, d\theta = \left( \int_0^{2\pi} \sin\theta \, d\theta \right) \left( \int_0^1 r^2 \, dr \right) \]

Evaluamos la integral trigonométrica periódica:

\[ \int_0^{2\pi} \sin\theta \, d\theta = [-\cos\theta]_0^{2\pi} = -1 - (-1) = 0 \]

Por tanto, el producto completo es idénticamente cero (por simetría respecto al eje \( x \)).

Resultado:

\[ \oint_C (y^2\,dx + 3xy\,dy) = 0 \]

Interpretación Matemática: La función integrando \( y \) es impar respecto al eje horizontal, anulando la integral doble en el disco completo simétrico.

Solución:

Datos:

• \( L(x, y) = x^2 y \implies \frac{\partial L}{\partial y} = x^2 \)
• \( M(x, y) = x y^2 \implies \frac{\partial M}{\partial x} = y^2 \)
• Frontera: \( y = x^2 \) y \( y = 4 \) (puntos de intersección en \( x = -2 \) y \( x = 2 \))

Fórmulas:

\[ W = \oint_C (L\,dx + M\,dy) = \iint_D (y^2 - x^2) \, dA \]

Sustitución y Operaciones:

Planteamos la integral iterada tipo I con \( x \in [-2, 2] \) e \( y \in [x^2, 4] \):

\[ \int_{-2}^2 \int_{x^2}^4 (y^2 - x^2) \, dy \, dx \]

Integrando respecto a \( y \):

\[ \left[ \frac{y^3}{3} - x^2 y \right]_{x^2}^4 = \left( \frac{64}{3} - 4x^2 \right) - \left( \frac{x^6}{3} - x^4 \right) = \frac{64}{3} - 4x^2 + x^4 - \frac{x^6}{3} \]

Aprovechando que la función es par en \( [-2, 2] \):

\[ 2 \int_0^2 \left( \frac{64}{3} - 4x^2 + x^4 - \frac{x^6}{3} \right) dx = 2 \left[ \frac{64x}{3} - \frac{4x^3}{3} + \frac{x^5}{5} - \frac{x^7}{21} \right]_0^2 \]

Evaluamos en \( x = 2 \):

\[ \frac{128}{3} - \frac{32}{3} + \frac{32}{5} - \frac{128}{21} = \frac{96}{3} + \frac{32}{5} - \frac{128}{21} = 32 + 6.4 - 6.0952 = 32.3048 = \frac{3392}{105} \]

Multiplicando por 2:

\[ W = 2 \times \frac{3392}{105} = \frac{6784}{105} \approx 64.61 \]

Resultado:

\[ W = \frac{6784}{105} \approx 64.61 \]

Interpretación Matemática: El trabajo neto ejercido por el campo de fuerzas a lo largo del bucle parabólico cerrado es positivo.

Solución:

Datos:

• Ecuación de la elipse con semiejes \( a > 0 \), \( b > 0 \)
• Parametrización en sentido antihorario: \( x(t) = a \cos t \), \( y(t) = b \sin t \) con \( t \in [0, 2\pi] \)
• Diferenciales: \( dx = -a \sin t \, dt \), \( dy = b \cos t \, dt \)

Fórmulas:

\[ \text{Área}(D) = \frac{1}{2} \oint_C (x \, dy - y \, dx) \]

Sustitución y Operaciones:

Sustituimos la parametrización en el integrando:

\[ x \, dy - y \, dx = (a \cos t)(b \cos t \, dt) - (b \sin t)(-a \sin t \, dt) \]

\[ x \, dy - y \, dx = ab (\cos^2 t + \sin^2 t) \, dt = ab \, dt \]

Integramos en el intervalo completo \( t \in [0, 2\pi] \):

\[ \text{Área}(D) = \frac{1}{2} \int_0^{2\pi} ab \, dt = \frac{1}{2} ab [t]_0^{2\pi} = \frac{1}{2} ab (2\pi) = \pi a b \]

Resultado:

\[ \text{Área}(\text{elipse}) = \pi a b \]

Interpretación Matemática: La integral de línea cerrada reproduce de forma natural y elegante la clásica fórmula de área elíptica.

Solución:

Datos:

• \( L(x, y) = 2xy + \sin x \implies \frac{\partial L}{\partial y} = 2x \)
• \( M(x, y) = x^2 + 5x \implies \frac{\partial M}{\partial x} = 2x + 5 \)
• Región rectangular \( D = [0, 3] \times [0, 2] \)

Fórmulas:

\[ \oint_C (L\,dx + M\,dy) = \iint_D \left( \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} \right) dA \]

Sustitución y Operaciones:

Evaluamos la resta de derivadas parciales:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = (2x + 5) - (2x) = 5 \]

Sustituimos en la integral doble:

\[ \iint_D 5 \, dA = 5 \iint_D dA = 5 \cdot \text{Área}(D) \]

El área del rectángulo es base por altura:

\[ \text{Área}(D) = 3 \times 2 = 6 \]

\[ \oint_C (L\,dx + M\,dy) = 5 \times 6 = 30 \]

Resultado:

\[ \oint_C (L\,dx + M\,dy) = 30 \]

Interpretación Matemática: El rotacional neto constante de magnitud 5 genera una circulación total de 30 unidades a lo largo de los 4 lados del rectángulo.

Solución:

Datos:

• \( L(x, y) = 3x^2 y + 2x \)
• \( M(x, y) = x^3 + 4y^3 \)
• Curva \( C \): Cualquier curva cerrada simple que encierra una región \( D \)

Fórmulas:

Criterio de campo conservativo en 2D:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = 0 \iff \mathbf{F} = \nabla f \]

Sustitución y Operaciones:

Derivamos componentes:

\[ \frac{\partial M}{\partial x} = \frac{\partial}{\partial x}(x^3 + 4y^3) = 3x^2 \]

\[ \frac{\partial L}{\partial y} = \frac{\partial}{\partial y}(3x^2 y + 2x) = 3x^2 \]

Restamos:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = 3x^2 - 3x^2 = 0 \]

Sustituimos en el Teorema de Green:

\[ \oint_C (L\,dx + M\,dy) = \iint_D 0 \, dA = 0 \]

Resultado:

\[ \oint_C \mathbf{F} \cdot d\mathbf{r} = 0 \]

Interpretación Matemática: El campo es estrictamente conservativo (su función potencial es \( f(x, y) = x^3 y + x^2 + y^4 \)), por lo que la integral de línea sobre cualquier camino cerrado es idénticamente nula.

Solución:

Datos:

• \( L(x, y) = y + \ln(1 + x^2) \implies \frac{\partial L}{\partial y} = 1 \)
• \( M(x, y) = 2x + e^{y^2} \implies \frac{\partial M}{\partial x} = 2 \)
• Curva orientada en sentido horario (orientación negativa: \( -C \))
• Región circular \( D \) con radio \( R = 2 \)

Fórmulas:

Para orientación negativa invertimos el signo:

\[ \oint_{-C} (L\,dx + M\,dy) = - \iint_D \left( \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} \right) dA \]

Sustitución y Operaciones:

\[ \frac{\partial M}{\partial x} - \frac{\partial L}{\partial y} = 2 - 1 = 1 \]

Sustituimos en la integral doble con el signo negativo:

\[ - \iint_D 1 \, dA = - \text{Área}(D) \]

El área del círculo con radio \( R = 2 \) es:

\[ \text{Área}(D) = \pi R^2 = \pi (2)^2 = 4\pi \]

Por tanto:

\[ I = -4\pi \]

Resultado:

\[ \oint_C (L\,dx + M\,dy) = -4\pi \]

Interpretación Matemática: El recorrido en sentido horario invierte la orientación canónica de la frontera, anteponiendo un signo negativo al área del disco.

Solución:

Datos:

• Parametrización: \( x = \cos^3 t \), \( y = \sin^3 t \)
• Diferenciales por regla de la cadena:
  \( dx = -3 \cos^2 t \sin t \, dt \),
  \( dy = 3 \sin^2 t \cos t \, dt \)

Fórmulas:

\[ \text{Área} = \frac{1}{2} \oint_C (x \, dy - y \, dx) \]

Sustitución y Operaciones:

Construimos el integrando:

\[ x \, dy = (\cos^3 t)(3 \sin^2 t \cos t) = 3 \cos^4 t \sin^2 t \]

\[ y \, dx = (\sin^3 t)(-3 \cos^2 t \sin t) = -3 \sin^4 t \cos^2 t \]

Restamos ambos términos:

\[ x \, dy - y \, dx = 3 \cos^2 t \sin^2 t (\cos^2 t + \sin^2 t) \, dt = 3 \cos^2 t \sin^2 t \, dt \]

Utilizamos la identidad del ángulo doble \( \sin(2t) = 2 \sin t \cos t \):

\[ \cos t \sin t = \frac{1}{2} \sin(2t) \implies \cos^2 t \sin^2 t = \frac{1}{4} \sin^2(2t) = \frac{1}{4} \left( \frac{1 - \cos(4t)}{2} \right) = \frac{1 - \cos(4t)}{8} \]

\[ x \, dy - y \, dx = \frac{3}{8}(1 - \cos(4t)) \, dt \]

Integramos de \( 0 \) a \( 2\pi \):

\[ \text{Área} = \frac{1}{2} \int_0^{2\pi} \frac{3}{8}(1 - \cos(4t)) \, dt = \frac{3}{16} [t - 0]_0^{2\pi} = \frac{3}{16}(2\pi) = \frac{3\pi}{8} \]

Resultado:

\[ \text{Área}(\text{astroide}) = \frac{3\pi}{8} \]

Interpretación Matemática: El astroide encierra exactamente tres octavos del área del círculo unitario que lo circunscribe.

Solución:

Datos:

• Campo: \( P(x, y) = 4x + \cos y \), \( Q(x, y) = 3y - \sin x \)
• Curva frontera: Elipse con semiejes \( a = 3 \), \( b = 2 \)

Fórmulas:

Forma normal del Teorema de Green:

\[ \Phi = \oint_C \mathbf{F} \cdot \mathbf{n} \, ds = \iint_D \left( \frac{\partial P}{\partial x} + \frac{\partial Q}{\partial y} \right) dA \]

Sustitución y Operaciones:

Calculamos la divergencia escalar:

\[ \frac{\partial P}{\partial x} = \frac{\partial}{\partial x}(4x + \cos y) = 4 \]

\[ \frac{\partial Q}{\partial y} = \frac{\partial}{\partial y}(3y - \sin x) = 3 \]

\[ \text{div}(\mathbf{F}) = 4 + 3 = 7 \]

Sustituimos en la integral doble:

\[ \Phi = \iint_D 7 \, dA = 7 \iint_D dA = 7 \cdot \text{Área}(\text{elipse}) \]

El área de la elipse es \( \pi a b = \pi (3)(2) = 6\pi \):

\[ \Phi = 7 \times 6\pi = 42\pi \]

Resultado:

\[ \Phi = 42\pi \approx 131.95 \]

Interpretación Matemática: El flujo divergente saliente total es equivalente al producto de la divergencia constante por el área contenida en la elipse.

Problemas Prácticos de Autoevaluación

Preguntas Frecuentes sobre el Teorema de Green

1. ¿Cuál es la principal ventaja de utilizar el Teorema de Green frente al cálculo directo?
El cálculo directo de una integral de línea cerrada exige dividir la trayectoria en varias curvas individuales, hallar una parametrización adecuada para cada una (con sus respectivos límites de tiempo), calcular derivadas temporales y resolver sumas de integrales compuestas. Con el Teorema de Green, todas esas etapas se reemplazan por una única integral doble sobre un recinto plano continuo, donde a menudo las derivadas parciales se simplifican de forma drástica (e incluso se vuelven constantes).

2. ¿Qué ocurre si la curva frontera se recorre en sentido horario?
Si la trayectoria se recorre en sentido horario, la región interior queda a la mano derecha del observador (orientación negativa). Por lo tanto, el Teorema de Green adquiere un signo negativo: \( \oint_{-C} (L\,dx + M\,dy) = - \iint_D (\partial M/\partial x - \partial L/\partial y)\,dA \). Siempre debes verificar la dirección de las flechas antes de integrar.

3. ¿Cómo se relaciona el Teorema de Green con el Teorema de Stokes y el Teorema de Gauss?
El Teorema de Green es exactamente el caso particular bidimensional (cuando la superficie es plana y está contenida en el plano \( xy \)) del Teorema de Stokes (para rotacionales en superficies curvas en 3D) y del Teorema de la Divergencia de Gauss (para flujos a través de fronteras en 3D). Dominar Green es el paso imprescindible para triunfar en cálculo vectorial avanzado y física de campos electromagnéticos.