4. Integral de flujo.

La discusión anterior sugiere que la suma ∑ ⁡ k=1mF ⋅NΔS k puede ser una aproximación aceptable de la masa total de fluido que atraviesa Sk en la unidad de tiempo.

Propósito: Visualizar como la suma ∑ ⁡ k=1mF ⋅NΔS k aproxima la la masa total de fluido que atraviesa Sk

Si F es la densidad de flujo de una corriente de fluido y N es el vector unitario normal a S definido por

N = ∂r ∂u ×∂r ∂v | |∂r ∂u ×∂r ∂v||,

entonces la masa total de fluido (fluido neto) que pasa por S por unidad de tiempo “en la dirección” de N es

∬SF ⋅NdS = ∬DF(r(u,v)) ⋅∂r ∂u ×∂r ∂v | |∂r ∂u ×∂r ∂v|| | |∂r ∂u ×∂r ∂v||dA = ∬DF(r(u,v)) ⋅∂r ∂u ×∂r ∂vdA

Orientación. Nuestra expresión para el flujo total lleva implícita la escongencia de uno de los dos vectores normales unitarios. Escoger un vector unitario para la región S es equivalente a "orientar" la región (como veremos má adelante). Esta escogencia de N decide el signo de F ⋅N. En lo que sigue, siempre vamos a escoger

N = ∂r ∂u ×∂r ∂v | |∂r ∂u ×∂r ∂v||

PIC

Fig. 8.22: S se orienta con el N escogido

En las aplicaciones es frecuente convenir en cuál N se escoge, para poder interpretar el resultado de ∬SF ⋅NdS como “flujo neto” en la dirección escogida.

La integral ∬SF ⋅NdS calcula el flujo neto: El flujo que pasa en la “dirección” de el vector N escogido, “suma” y el flujo que pasa en la “dirección contraria”, “resta“. La suma de todo esto es el flujo neto.

Caso z = f(x,y) Como consecuencia tenemos que si S : z = f(x,y) con f de clase C1 sobre D¯, se puede parametrizar S con r(x,y) = xi + yj + f(x,y)k y entonces

∬SF ⋅NdS = ∬DF(r(x,y)) ⋅∂r ∂x ×∂r ∂ydA = ∬DxyF(x,y,z) ⋅ (−fx, − fy,1)dA

Integral de Flujo−Proyectando sobre varios varios planos

a
Proyectando sobre XY : Si S : z = z(x,y) o S : G(x,y,z) = 0, con (x,y) ∈Dxy

∬SF ⋅NdS =∬ DxyF(x,y,z(x,y)) ⋅ (−zx,−zy,1)dA

o, en "versión implícita",

∬SF ⋅NdS =∬ DxyF(x,y,z(x,y)) ⋅ (Gx,Gy,Gz) 1 GzdA con Gz≠0 en Dxy

b
Proyectando sobre XZ: Si S : y = y(x,z) o S : G(x,y,z) = 0, con (x,z) ∈Dxz

∬SF ⋅NdS =∬ DxzFdA(x,y(x,z),z) ⋅ (−yx,1,−yz)dA

o, en "versión implícita",

∬SF ⋅NdS =∬ DxzF(x,y(x,z),z) ⋅ (Gx,Gy,Gz) 1 GydA con Gy≠0 en Dxz

c
Proyectando sobre Y Z: Si S : x = x(y,z) o S : G(x,y,z) = 0, con (y,z) ∈Dyz

∬SF ⋅NdS =∬ DyzF(x(y,z),y,z) ⋅ (1,−xy,−xz)dA

o, en "versión implícita",

∬SF ⋅NdS =∬ DyzF(x(y,z),y,z) ⋅ (Gx,Gy,Gz) 1 GxdA con Gx≠0 en Dyz

La integral ∬SF ⋅NdS puede diferir en el signo al cambiar el plano de proyección. Esto es así porque los vectores (−zx,−zy,1), (−yx,1,−yz) y (1,−xy,−xz) pueden invertir la dirección. En todo caso, ∬SF ⋅NdS se interpreta como el flujo neto en la dirección del vector N escogido.

Ejemplo259Flujo
Calcular ∬SF ⋅NdS si
a
F(x,y,z) = (z + 1)k
b
S es la superficie z = 2 + y con dominio Dxy : x2 + y2 = 1.

PIC

Fig. 8.23: Superficie S y su proyección


Solución. La superficie S tiene ecuación z = 2 + y. Además, proyectando sobre le plano XY , Dxy es el círculo de radio 1. Entonces,

∬SF ⋅NdS = ∬Dxy(0,0,z + 1) ⋅ (−zx, − zy,1)dA = ∬Dxy(0,0,2 + y + 1) ⋅ (0,−1,1)dA = ∬Dxyy + 3dA = ∫ 02π ∫ 01(3 + rsen ⁡ 𝜃)rdrd𝜃 = ∫ 02π9 + 2sen ⁡ (𝜃) 6 d𝜃 = 3π

Ejemplo260Flujo
Sea
a
F(x,y,z) = (0,y,0)
b
S el cilindro z = 2 −x2 2 , desde y = 0 hasta y = 2, como se ve en la Figura 8.24.

Calcular ∬SF ⋅NdS

PIC

Fig. 8.24: Superficie S.


Solución.

Intuitivamente, el flujo no pasa a través de la superficie S, así que la integral de flujo debería ser 0.

En este caso solo se puede proyectar sobre Y Z o XY . La proyección sobre Y Z es un rectángulo.

Sea S : G(x,y,z) = 0 con G(x,y,z) = z − 2 + x2 2 . Entonces, 

∬SF ⋅NdS = ∬DyzFdA ⋅∇ ⁡G Gx dA = ∫ 02 ∫ 03∕2(0,y,0) ⋅ (x,0,1)1 xdzdy = ∫ 02 ∫ 03∕20dzdy = 0

Ejemplo261Flujo
Calcular ∬SF ⋅NdS si
a
F(x,y,z) = (x,yz,xy)
b
S el cilindro de ecuación y = 4 − x2 limitado por el plano x + z = 2, tal y como se muestra en la Figura 8.25.

PIC

Fig. 8.25: Superficie S


Solución.

La superficie S tiene ecuación G(x,y,z) = 0 con G(x,y,z) = y + x2 − 4. La proyección de S sobre el plano XZ es un triángulo.

∬SF ⋅NdS = ∬Dxz(x,yz,xy) ⋅∇ ⁡G Gy dA = ∬Dxz(x,yz,xy) ⋅ (2x,1,0)dA = ∫ 02 ∫ 02−x (2x2 + yz)dzdx = ∫ 02 ∫ 02−x (2x2 + (4 − x2)z)dzdx = 112 15

PIC

Fig. 8.26: Superficie S

Si calculamos proyectando sobre Y Z tenemos S : x = 4 − y.

∬SF ⋅NdS = ∫ 02 ∫ 4−(z−2)24 (4 − y + yz 24 − y)dydz

Ejemplo262flujo
Calcular ∬SF ⋅NdS si
a
F(x,y,z) = 4xzi + yz3j + z2k
b
S es la superficie z2 + y2 + x2 = 4 entre z = 1 y z = 2.

PIC

Fig. 8.27: Superficie S


Solución.

La superficie S tiene ecuación G(x,y,z) = 0 con G(x,y,z) = x2 + y2 + z2 − 4. La proyección Dxy es el círculo x2 + y2 = 3. Entonces,

∬SF ⋅NdS = ∬Dxy(4xz,yz3,z2) ⋅∇ ⁡G Gz dA = ∬Dxy(4xz,yz3,z2) ⋅ (x z,y z,1)dA = ∬Dxy (4x2 + y2z2 + z2) dA = ∬Dxy (4x2 + y2(4 − x2 − y2) + 4 − x2 − y2) dA = 1 2∫ 02π ∫ 0 3(8 + 6r2 − r4 + r4 cos ⁡ 2𝜃)rdrd𝜃 = 21π.

5.
Ejercicios

5.1
Calcule I = ∬SF ⋅NdS donde F es el campo vectorial dado por F(x,y,z) = yi −xj + 8zk y S la parte de la esfera de ecuación x2 + y2 + z2 = 9 que se encuentra dentro del cilindro x2 + y2 = 4, como se observa en la figura.

PIC

Fig. 8.28: Superficie S
Vamos a proyectar sobre el plano XY . La proyección es el círculo x2 + y2 ≤ 2.

Como z = 9 − x2 − y2 entonces podemos poner N = (−x∕z, − y∕z,1) ||(−x∕z, − y∕z,1)||. Luego,

∬SF ⋅NdS = ∬S(y,−x,8z) ⋅ (−x∕z, − y∕z,1) ||(−x∕z, − y∕z,1)||||(−x∕z, − y∕z,1)||dA = ∬D8zdA = ∫ 02π ∫ 0289 − r2rdrd𝜃 = 16π(53∕2 − 93∕2) −3

5.2
Sea F(x,y,z) = xy2i + x2yj + yk. Sea S es la frontera del sólido E limitado por
S1 : x2 + y2 = 1,S2 : z = 1 y S3 : z = −1

como se ve en la figura. Calcule ∬SF ⋅NdS donde N es el vector normal unitario exterior al sólido E.

PIC

Fig. 8.29: Curva C
Observe que debe calcular con N siempre exterior al sólido E.

∬SF ⋅NdS = π

5.3
Sea F(x,y,z) = xi + yj + zk y S es la superficie de ecuación z = 1 + x2 + y2, con 1 ≤ z ≤ 3, tal y como se muestra en la figura. Calcular ∬SF ⋅NdS.

PIC

Fig. 8.30: Superficie E.
∬SF ⋅NdS = ∬Rxy(x,y,1 + x2 + y2) ⋅ (−2x,−2y,1)dA = ∫ 02π ∫ 02(1 − r2)rdrd𝜃

5.4
Sea F(x,y,z) = (0,x + y,z). Calcule la integral de superficie ∬SF ⋅NdS, donde S es el trozo de cilindro de ecuación z = 1 − (x − 2)3 que está limitado por los planos y = 0, y = 3, z = 1 y z = 2 tal y como se muestra en la figura.

PIC

∬SF ⋅NdS = ∫ 13 ∫ 03(0,x + y,1 − (x − 2)3) ⋅ (−3(x − 2)2,0,1)dydx = ∫ 13 ∫ 03(1 − (x − 2)3)dydx = ∫ 13 ∫ 03(−x3 + 6x2 − 12x + 9)dydx = 3∫ 13(−x3 + 6x2 − 12x + 9)dx = 3 (−x4 4 + 2x3 − 6x2 + 9x)| 03 = 6

5.5
Sea S la superficie de ecuación z = x2 + y2 que se encuentra limitada por los planos z = 1, z = 3, y = x y el plano x = 0, tal y como se muestra en la figura. Calcule ∬SF ⋅NdS si F(x,y,z) = (x,y,z2)

PIC

Fig. 8.31: Superficie S
Proyectamos sobre XY

AS = ∬SF ⋅NdS = ∫ π∕4π∕2 ∫ 1 3 1 + 4r2rdrd𝜃

5.6
Sea F(x,y,z) = (xy,x,z + 1) y sea S la porción de superficie de ecuación z = 4 − y2 limitada por las superficies z = 3, x = 4, z = 0 y x = y, tal y como se muestra en la Figura 8.32 . Calcular ∬SF ⋅NdS.

PIC

Fig. 8.32: Superficie S
Proyectamos sobre XY .
∬SF ⋅NdS = ∫ 12 ∫ y4(xy,x,z + 1) ⋅ (0,2y,1)dxdy = ∫ 12 ∫ y42xy + z + 1dxdy = ∫ 12 ∫ y42xy + 4 − y2 + 1dxdy

PIC

Fig. 8.33: Superficie S