Integrales de superficie y flujo

Cálculo de flujo vectorial a través de superficies abiertas y cerradas en física y electromagnetismo.

Que es

Una integral de superficie extiende la integración a variedades bidimensionales orientables SS curvadas en el espacio tridimensional R3\mathbb{R}^3. Se distinguen dos formas fundamentales:

  1. Integral de superficie de un campo escalar ff: ∬Sf dS=∬Df(r⃗(u,v))∥∂r⃗∂u×∂r⃗∂v∥du dv\iint_S f \, dS = \iint_D f(\vec{r}(u,v)) \left\| \frac{\partial \vec{r}}{\partial u} \times \frac{\partial \vec{r}}{\partial v} \right\| du \, dv Permite calcular el área total de una superficie curvada (f=1f=1), la masa de una cáscara con densidad superficial σ(x,y,z)\sigma(x,y,z) o la carga electrostática neta depositada sobre un conductor.
  2. Integral de flujo de un campo vectorial F⃗\vec{F}: Φ=∬SF⃗⋅dS⃗=∬S(F⃗⋅n^) dS\Phi = \iint_S \vec{F} \cdot d\vec{S} = \iint_S (\vec{F} \cdot \hat{n}) \, dS donde n^\hat{n} es el vector unitario normal a la superficie en cada punto y dS⃗=n^dSd\vec{S} = \hat{n} dS es el elemento vectorial de área. Cuantifica el paso neto o transferencia de campo a través de la superficie por unidad de tiempo.

Por qué aparece en física

El flujo a través de superficies es el pilar cuantitativo de las leyes de campo y la termofísica:

  • Flujo eléctrico y Ley de Gauss: ΦE=∮SE⃗⋅dS⃗=Qencε0\Phi_E = \oint_S \vec{E} \cdot d\vec{S} = \frac{Q_{\text{enc}}}{\varepsilon_0}. El flujo neto saliente a través de cualquier superficie cerrada es proporcional a la carga eléctrica neta confinada en su interior.
  • Flujo magnético: ΦB=∬SB⃗⋅dS⃗\Phi_B = \iint_S \vec{B} \cdot d\vec{S}. Su variación temporal induce corriente eléctrica (Ley de Faraday: E=−dΦB/dt\mathcal{E} = -d\Phi_B/dt).
  • Caudal volumétrico de fluidos: Q=∬Sv⃗⋅dS⃗Q = \iint_S \vec{v} \cdot d\vec{S}. Mide los metros cúbicos por segundo de líquido que atraviesan una sección transversal de tubería o canal.
  • Potencia térmica radiada o conducida: P=∬Sq⃗⋅dS⃗P = \iint_S \vec{q} \cdot d\vec{S}. Energía térmica total transferida a través de la piel de un reactor o nave espacial.

Idea visual o geométrica

Imagina una red de pesca sumergida en la corriente de un río caudaloso:

  • Si colocas la red perpendicular a la corriente, el caudal de agua que la atraviesa es máximo (v⃗∥n^  ⟹  v⃗⋅n^=v\vec{v} \parallel \hat{n} \implies \vec{v}\cdot\hat{n} = v).
  • Si orientas la red paralela a la corriente, el agua fluye a los lados sin atravesar la malla (v⃗⊥n^  ⟹  v⃗⋅n^=0\vec{v} \perp \hat{n} \implies \vec{v}\cdot\hat{n} = 0).
  • Si inclinas la red un ángulo θ\theta, el flujo efectivo es proporcional a cos⁡θ\cos\theta. La integral de flujo suma las componentes normales de la velocidad del agua sobre cada parche elemental de la red.

Lenguaje simbólico mínimo

  • Parametrización de superficie: r⃗(u,v)=x(u,v)i^+y(u,v)j^+z(u,v)k^\vec{r}(u,v) = x(u,v)\hat{i} + y(u,v)\hat{j} + z(u,v)\hat{k} con (u,v)∈D(u,v) \in D.
  • Vector normal fundamental no normalizado: N⃗=∂r⃗∂u×∂r⃗∂v\vec{N} = \frac{\partial \vec{r}}{\partial u} \times \frac{\partial \vec{r}}{\partial v}.
  • Elemento de área vectorial: dS⃗=N⃗ du dv=n^ dSd\vec{S} = \vec{N} \, du \, dv = \hat{n} \, dS.
  • Para superficies dadas explícitamente como z=g(x,y)z = g(x,y): dS⃗=(−∂g∂xi^−∂g∂yj^+k^)dx dyd\vec{S} = \left( -\frac{\partial g}{\partial x} \hat{i} - \frac{\partial g}{\partial y} \hat{j} + \hat{k} \right) dx \, dy
  • Para esferas de radio RR: dS⃗=R2sin⁡θ dθ dφ e^rd\vec{S} = R^2 \sin\theta \, d\theta \, d\varphi \, \hat{e}_r.

Que significa el resultado

El escalar de flujo Φ\Phi tiene unidades de [F⃗]×[m2][\vec{F}] \times [\text{m}^2]. Si Φ>0\Phi > 0, hay un transporte o irradiación neta en la dirección del vector normal elegido; si Φ<0\Phi < 0, el flujo entra en sentido contrario. Para una superficie cerrada, ∮SF⃗⋅dS⃗=0\oint_S \vec{F}\cdot d\vec{S} = 0 indica que todo el flujo que ingresa sale sin retención.

Cuando no debes usarlo

  • No olvides definir explícitamente la orientación de la superficie. En superficies cerradas, la convención universal exige que n^\hat{n} apunte hacia el exterior. En superficies abiertas, la orientación debe ser coherente con la regla de la mano derecha respecto a su curva frontera.
  • No intentes proyectar una superficie vertical (como un cilindro con paredes paralelas al eje zz) sobre el plano xyxy usando dxdydx dy; la normal fundamental sería ortogonal al eje zz y el diferencial proyectivo se indetermina. Usa coordenadas cilíndricas.

Resumen operativo

Para calcular el flujo ∬SF⃗⋅dS⃗\iint_S \vec{F}\cdot d\vec{S}: parametriza la superficie en dos variables (u,v)(u,v), calcula el producto vectorial fundamental ∂ur⃗×∂vr⃗\partial_u \vec{r} \times \partial_v \vec{r}, evalúa F⃗\vec{F} sobre la superficie, calcula el producto escalar F⃗⋅(∂ur⃗×∂vr⃗)\vec{F} \cdot (\partial_u \vec{r} \times \partial_v \vec{r}) e integra la función resultante sobre la región paramétrica bidimensional DD.

Procedimiento

Antes de calcular

Verifica si la superficie SS es cerrada. Si es cerrada, evalúa si es más rápido aplicar 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 en lugar de integrar directamente sobre múltiples caras.

Paso 1: identifica el objeto matemático

Especifica la ecuación geométrica de la superficie SS, su orientación (normal hacia arriba, exterior, o radial) y la expresión del campo F⃗\vec{F}.

Paso 2: elige la operación

Selecciona la parametrización más conveniente:

  • Si z=g(x,y)z = g(x,y): dS⃗=(−∂xg,−∂yg,1)dxdyd\vec{S} = (-\partial_x g, -\partial_y g, 1) dx dy.
  • Si es una esfera de radio RR: dS⃗=R2sin⁡θ dθ dφ e^rd\vec{S} = R^2 \sin\theta \, d\theta \, d\varphi \, \hat{e}_r.
  • Si es un cilindro de radio RR: dS⃗=R dθ dz e^rd\vec{S} = R \, d\theta \, dz \, \hat{e}_r.

Paso 3: conserva unidades y signos

Calcula el vector normal fundamental y confirma que apunte en la dirección de la orientación requerida. Si apunta en sentido opuesto, invierte su signo multiplicando por −1-1.

Paso 4: calcula

Sustituye la parametrización en el campo vectorial F⃗\vec{F} y evalúa el producto escalar: Φ=∬D[F⃗(r⃗(u,v))⋅N⃗(u,v)] du dv\Phi = \iint_D [\vec{F}(\vec{r}(u,v)) \cdot \vec{N}(u,v)] \, du \, dv Calcula la integral doble resultante aplicando el Teorema de Fubini en límites iterados.

Paso 5: interpreta el resultado

Analiza el flujo obtenido:

  • Φ>0\Phi > 0: flujo saliente neto o emisión.
  • Φ<0\Phi < 0: flujo entrante neto o absorción.

Checklist de autocontrol

  • ¿El vector normal está correctamente orientado (hacia afuera si la superficie es cerrada)?
  • ¿Se incluyó el jacobiano métrico (R2sin⁡θR^2 \sin\theta en esferas, RR en cilindros)?
  • ¿Las unidades del flujo son [campo]×[m2][\text{campo}] \times [\text{m}^2]?

Ejemplo guiado

Situación física

Un campo de velocidades de un líquido incompresible en una tobera viene modelado por: v⃗(x,y,z)=(2x)i^+(2y)j^+(3)k^[m/s]\vec{v}(x,y,z) = (2x) \hat{i} + (2y) \hat{j} + (3) \hat{k} \quad [\text{m/s}] Se desea calcular el caudal volumétrico total Q=∬Sv⃗⋅dS⃗Q = \iint_S \vec{v} \cdot d\vec{S} que fluye a través de la superficie paraboloide abierta z=4−x2−y2z = 4 - x^2 - y^2 para z≥0z \ge 0, orientada con vector normal apuntando hacia arriba (n^⋅k^>0\hat{n}\cdot\hat{k} > 0).

Traducción matemática

La superficie está dada explícitamente por z=g(x,y)=4−x2−y2z = g(x,y) = 4 - x^2 - y^2. Su dominio de proyección DD en el plano xyxy (donde z≥0z \ge 0) es el círculo x2+y2≤4x^2 + y^2 \le 4 (radio R=2R = 2). El vector de área orientado hacia arriba es: dS⃗=(−∂g∂xi^−∂g∂yj^+k^)dx dy=(2xi^+2yj^+k^) dx dyd\vec{S} = \left( -\frac{\partial g}{\partial x} \hat{i} - \frac{\partial g}{\partial y} \hat{j} + \hat{k} \right) dx \, dy = (2x \hat{i} + 2y \hat{j} + \hat{k}) \, dx \, dy

Calculo paso a paso

  1. Producto escalar del campo con el elemento de área: v⃗⋅dS⃗=[(2x)i^+(2y)j^+3k^]⋅[(2x)i^+(2y)j^+k^] dx dy\vec{v} \cdot d\vec{S} = [(2x) \hat{i} + (2y) \hat{j} + 3 \hat{k}] \cdot [(2x) \hat{i} + (2y) \hat{j} + \hat{k}] \, dx \, dy v⃗⋅dS⃗=[4x2+4y2+3] dx dy\vec{v} \cdot d\vec{S} = [4x^2 + 4y^2 + 3] \, dx \, dy

  2. Transformación a coordenadas polares en el plano xyxy: x=rcos⁡θ,y=rsin⁡θ,x2+y2=r2x = r \cos\theta, \quad y = r \sin\theta, \quad x^2 + y^2 = r^2 dx dy=r dr dθdx \, dy = r \, dr \, d\theta Límites de integración: r∈[0,2]r \in [0, 2], θ∈[0,2π]\theta \in [0, 2\pi].

  3. Planteamiento de la integral doble: Q=∫02πdθ∫02(4r2+3)r dr=2π∫02(4r3+3r) drQ = \int_0^{2\pi} d\theta \int_0^2 (4r^2 + 3) r \, dr = 2\pi \int_0^2 (4r^3 + 3r) \, dr

  4. Evaluación de la integral radial: ∫02(4r3+3r) dr=[r4+3r22]02=24+3(22)2=16+122=16+6=22\int_0^2 (4r^3 + 3r) \, dr = \left[ r^4 + \frac{3r^2}{2} \right]_0^2 = 2^4 + \frac{3(2^2)}{2} = 16 + \frac{12}{2} = 16 + 6 = 22

  5. Caudal total final: Q=2π×22=44π≈138.23 m3/sQ = 2\pi \times 22 = 44\pi \approx 138.23 \text{ m}^3/\text{s}

Lectura del resultado

A través de la cúpula paraboloide fluyen 138.23138.23 metros cúbicos de fluido por segundo. La normal hacia arriba captura el escape vertical y radial del fluido.

Error común evitado

Olvidar el jacobiano rr al transformar a polares, lo que llevaría a integrar (4r2+3)dr(4r^2+3)dr en vez de (4r3+3r)dr(4r^3+3r)dr, falseando completamente el cálculo volumétrico.

Ejemplo inverso

Resultado dado

El flujo electrostático total a través de una superficie esférica cerrada de radio R=0.50 mR = 0.50 \text{ m} que encierra un dispositivo electrónico es ΦE=∮SE⃗⋅dS⃗=+4.50×104 N⋅m2/C\Phi_E = \oint_S \vec{E} \cdot d\vec{S} = +4.50 \times 10^4 \text{ N}\cdot\text{m}^2/\text{C}.

Que pregunta física permite responder

¿Cuál es la carga eléctrica neta confinada dentro del dispositivo? Por la ley de Gauss: ΦE=Qencε0  ⟹  Qenc=ε0ΦE\Phi_E = \frac{Q_{\text{enc}}}{\varepsilon_0} \implies Q_{\text{enc}} = \varepsilon_0 \Phi_E Sustituyendo ε0≈8.854×10−12 F/m\varepsilon_0 \approx 8.854 \times 10^{-12} \text{ F/m}: Qenc=(8.854×10−12)×(4.50×104)≈3.984×10−7 C=+0.398 μCQ_{\text{enc}} = (8.854 \times 10^{-12}) \times (4.50 \times 10^4) \approx 3.984 \times 10^{-7} \text{ C} = +0.398 \,\mu\text{C} La integral de flujo permite cuantificar la carga interna total sin necesidad de abrir el dispositivo ni conocer su geometría interna.

Errores frecuentes

Error 1: Confundir la normal hacia arriba con la normal hacia abajo

Invertir el signo del producto vectorial fundamental ∂ur⃗×∂vr⃗\partial_u \vec{r} \times \partial_v \vec{r}. Produce un cambio espurio de signo en todo el flujo.

Error 2: Olvidar cerrar la superficie al aplicar teoremas

Intentar aplicar la Ley de Gauss a una superficie abierta (como un hemisferio o un cilindro sin tapas) olvidando incluir el flujo a través de las bases planas.

Error 3: Omitir los factores métricos en coordenadas esféricas

Escribir dS=dθdφdS = d\theta d\varphi en vez de R2sin⁡θdθdφR^2 \sin\theta d\theta d\varphi. Destruye la homogeneidad dimensional y el valor del área.

Error 4: Confundir flujo escalar con flujo vectorial

Olvidar que Φ=∬F⃗⋅dS⃗\Phi = \iint \vec{F}\cdot d\vec{S} es un número escalar puro, no un vector.

Error 5: Integrar sobre un dominio de proyección incorrecto

Integrar x,yx, y en un rectángulo cuando la superficie es un cono o casquete con base circular, violando los límites geométricos.

Regla rápida de detección

Si el flujo tiene componentes vectoriales o no tiene unidades de [campo]×[m2][\text{campo}] \times [\text{m}^2], el planteamiento algebraico es erróneo.

Aplicaciones físicas

1. Electromagnetismo: Ley de Gauss y flujo de campo

La ley de Gauss ∮E⃗⋅dS⃗=Q/ε0\oint \vec{E}\cdot d\vec{S} = Q/\varepsilon_0 permite calcular campos eléctricos analíticamente en geometrías de alta simetría (esferas, cilindros coaxiales, planos cargados).

2. Inducción electromagnética y transformadores

La ley de Faraday E=−dΦB/dt\mathcal{E} = -d\Phi_B/dt calcula la tensión generada en bobinas giratorias en alternadores hidroeléctricos y eólicos.

3. Dinámica de fluidos y caudales en ingeniería

En hidráulica, el aforo de ríos y el cálculo de pérdida de carga en tuberías se fundamenta en la integral de flujo de velocidad Q=∬v⃗⋅dS⃗Q = \iint \vec{v}\cdot d\vec{S}.

4. Radiación estelar y paneles solares espaciales

La potencia solar interceptada por satélites orbitales se calcula integrando el vector de Poynting S⃗=1μ0(E⃗×B⃗)\vec{S} = \frac{1}{\mu_0}(\vec{E}\times\vec{B}) sobre la superficie de las celdas fotovoltaicas.

5. Termodinámica de edificios y aislamiento térmico

La pérdida energética invernal de una vivienda se calcula mediante la integral de flujo térmico ∬q⃗⋅dS⃗\iint \vec{q}\cdot d\vec{S} sobre muros, ventanas y techos.

Practica

Problema 1: Flujo eléctrico a través de una superficie esférica (Ley de Coulomb)

Una carga puntual q=+5.0 nCq = +5.0 \text{ nC} se sitúa en el origen. El campo eléctrico es E⃗=q4πε0r2e^r\vec{E} = \frac{q}{4\pi\varepsilon_0 r^2} \hat{e}_r. Calcule el flujo eléctrico ∮SE⃗⋅dS⃗\oint_S \vec{E} \cdot d\vec{S} a través de una esfera concéntrica de radio R=0.20 mR = 0.20 \text{ m} mediante integración directa de superficie.

Solución:

  1. En la esfera de radio RR, el campo es puramente radial y de magnitud constante: E⃗(R)=q4πε0R2e^r\vec{E}(R) = \frac{q}{4\pi\varepsilon_0 R^2} \hat{e}_r
  2. El elemento de superficie vectorial es dS⃗=R2sin⁡θ dθ dφ e^rd\vec{S} = R^2 \sin\theta \, d\theta \, d\varphi \, \hat{e}_r.
  3. Producto escalar: E⃗⋅dS⃗=(q4πε0R2)(R2sin⁡θ dθ dφ)=q4πε0sin⁡θ dθ dφ\vec{E} \cdot d\vec{S} = \left( \frac{q}{4\pi\varepsilon_0 R^2} \right) (R^2 \sin\theta \, d\theta \, d\varphi) = \frac{q}{4\pi\varepsilon_0} \sin\theta \, d\theta \, d\varphi
  4. Integrando en θ∈[0,π]\theta \in [0, \pi] y φ∈[0,2π]\varphi \in [0, 2\pi]: ΦE=q4πε0∫02πdφ∫0πsin⁡θ dθ=q4πε0(2π)(2)=qε0\Phi_E = \frac{q}{4\pi\varepsilon_0} \int_0^{2\pi} d\varphi \int_0^\pi \sin\theta \, d\theta = \frac{q}{4\pi\varepsilon_0} (2\pi)(2) = \frac{q}{\varepsilon_0} ΦE=5.0×10−9 C8.854×10−12 F/m≈564.7 N⋅m2/C\Phi_E = \frac{5.0 \times 10^{-9} \text{ C}}{8.854 \times 10^{-12} \text{ F/m}} \approx 564.7 \text{ N}\cdot\text{m}^2/\text{C}

Problema 2: Flujo magnético en un plano inclinado

Un campo magnético uniforme B⃗=0.40k^ T\vec{B} = 0.40 \hat{k} \text{ T} atraviesa una espira triangular plana con vértices en (0,0,0)(0,0,0), (3,0,0)(3,0,0) y (0,2,2)(0,2,2) (en metros). Calcule el flujo magnético ΦB=∬SB⃗⋅dS⃗\Phi_B = \iint_S \vec{B} \cdot d\vec{S}.

Solución:

  1. Vectores directores de los lados: u⃗=(3,0,0)−(0,0,0)=3i^\vec{u} = (3,0,0) - (0,0,0) = 3\hat{i} v⃗=(0,2,2)−(0,0,0)=2j^+2k^\vec{v} = (0,2,2) - (0,0,0) = 2\hat{j} + 2\hat{k}
  2. Vector normal fundamental y área vectorial: A⃗=12(u⃗×v⃗)=12∣i^j^k^300022∣=12(−6j^+6k^)=−3j^+3k^[m2]\vec{A} = \frac{1}{2}(\vec{u} \times \vec{v}) = \frac{1}{2} \begin{vmatrix} \hat{i} & \hat{j} & \hat{k} \\ 3 & 0 & 0 \\ 0 & 2 & 2 \end{vmatrix} = \frac{1}{2} (-6\hat{j} + 6\hat{k}) = -3\hat{j} + 3\hat{k} \quad [\text{m}^2]
  3. Flujo para campo uniforme: ΦB=B⃗⋅A⃗=(0.40k^)⋅(−3j^+3k^)=0.40×3=1.20 Wb(Tesla⋅m2)\Phi_B = \vec{B} \cdot \vec{A} = (0.40 \hat{k}) \cdot (-3\hat{j} + 3\hat{k}) = 0.40 \times 3 = 1.20 \text{ Wb} \quad (\text{Tesla}\cdot\text{m}^2)

Problema 3: Flujo de calor a través de la pared de un tubo cilíndrico

Un tubo de vapor de radio R=0.10 mR = 0.10 \text{ m} y longitud L=2.0 mL = 2.0 \text{ m} disipa calor con densidad de flujo radial q⃗(r)=500re^r W/m2\vec{q}(r) = \frac{500}{r} \hat{e}_r \text{ W/m}^2. Calcule la potencia térmica total disipada P=∬Sq⃗⋅dS⃗P = \iint_S \vec{q} \cdot d\vec{S} a través de la pared lateral.

Solución:

  1. En la pared lateral cilíndrica (r=Rr = R): q⃗(R)=5000.10e^r=5000e^r[W/m2]\vec{q}(R) = \frac{500}{0.10} \hat{e}_r = 5000 \hat{e}_r \quad [\text{W/m}^2]
  2. Elemento de área normal exterior: dS⃗=R dθ dz e^r=(0.10)dθ dz e^rd\vec{S} = R \, d\theta \, dz \, \hat{e}_r = (0.10) d\theta \, dz \, \hat{e}_r
  3. Integrando: P=∫0Ldz∫02πdθ(5000)(0.10)=500(2π)(L)=500×2π×2.0=2000π≈6283.18 W≈6.28 kWP = \int_0^L dz \int_0^{2\pi} d\theta (5000)(0.10) = 500 (2\pi)(L) = 500 \times 2\pi \times 2.0 = 2000\pi \approx 6283.18 \text{ W} \approx 6.28 \text{ kW}