FENOMENOS DE TRANSPORTE

Luis Carrasco Venegas

Unidad de Posgrado de la Facultad de Ingeniería Química, Universidad Nacional del Callao
Unidad de Investigación de la Facultad de Ingeniería Química, Facultad de Ingeniería Química, Universidad Nacional del Callao
Correo: lacarrascov@[Link]
ORCID: [Link]
Web: [Link]
Docente Renacyt

CAPITULO VII: ANEXOS

7.1.- INTRODUCCION AL ALGEBRA VECTORIAL

Las magnitudes físicas
Las magnitudes físicas conocidas son las escalares, las vectoriales y las tensoriales, cuya
notación es la siguiente:
Escalar: s

Ejemplos de magnitudes escalares son el tiempo (t ) , la masa ( m) , el volumen (v) , la energía

( E ) , la presión ( P) , etc. A estas magnitudes se les conoce también como tensores de orden cero.

Vectorial: V (negrita; en este caso se utilizará V en color azul en lugar de negrita); también se
utiliza V (con línea inferior en la parte baja de la letra que identifica al vector) o V (con línea
superior sobre la letra que identifica al vector).

Ejemplos de magnitudes vectoriales son la velocidad (V ) , el peso (W ) , la fuerza ( F ) , el campo

eléctrico ( E ) , el campo magnético ( B) , etc. A estas magnitudes se les conoce también como
tensores de orden uno.

Tensorial:  (con dos líneas en la parte baja de la letra que identifica al tensor)

Son ejemplos de magnitudes tensoriales el producto simple de vectores o producto diádico (VW )
, el tensor de esfuerzos ( ) , la gradiente de velocidad (V ) , etc. Estas magnitudes son tensores
de orden dos o simplemente tensores.

Producto vectorial y magnitud física del resultado

En la primera columna de la tabla 1 se observa los productos que se aplican a los vectores y
tensores; al lado derecho se observa el resultado de dichas operaciones, las mismas que pueden
resultar ser magnitudes escalares, vectoriales o tensoriales. Este criterio es importante para darse
cuenta rápidamente de la magnitud física obtenida del desarrollo de dichos productos.
Tabla 1: Producto a aplicar y orden del resultado

PRODUCTO RESULTADO
•  −2
  −1
 :
 −4
Ninguno


Nota: El producto doble punto solo es aplicable a los tensores.

A continuación, se ilustra mediante ejemplos la aplicación del algebra vectorial en coordenadas
rectangulares, así como en coordenadas cilíndricas y esféricas.
Ejemplo 1
Utilizando la tabla 1, indique la magnitud física resultante al aplicar el producto escalar
(producto punto) a los vectores V y W .

Solución:
El producto punto de dos vectores esta dado por:

(V • W )

Orden 1 Orden 1

El resultado según la tabla 1 es  −2 , es decir (1+1)-2=0. Esto significa que es un tensor de

orden cero, es decir, un escalar.
Ejemplo 2
Haciendo uso de la tabla 1, indique la magnitud física resultante al aplicar el producto vectorial
de vectores (producto aspa) V y W .

Solución:
El producto aspa esta dado por:

(V W )

Orden 1 Orden 1
Con la ayuda de la tabla 1, para el producto aspa se tiene la suma de ordenes menos 1, es decir:
 −1 ; en consecuencia (1+1)-1=1. Este resultado nos produce un tensor de orden 1 que equivale
a un vector.
Ejemplo 3
Utilizando la tabla 1, indique la magnitud física resultante al aplicar el producto escalar
(producto doble punto) a los tensores  y 

Solución:
El producto doble punto está dado por:

 :

Orden 2 Orden 2
En la Tabla 1 se observa que el resultado de multiplicar escalarmente dos tensores, el orden del
resultado es  −4 ; es decir (2+2)-4=0; esto implica que es un tensor de orden cero; es decir, un
escalar
Ejemplo 4
Utilizando la tabla 1, indique la magnitud física resultante al aplicar el producto diádico (producto
simple) de dos vectores V y W

Solución:
El producto diádico de dos vectores esta dado por:

(VW )
Orden 1 Orden 1

En la tabla 1 se observa que el resultado para el producto diádico es  es decir (1+1) =2; este
resultado implica que el producto es un tensor de orden 2 o simplemente un tensor.

Representación de las magnitudes físicas desde el punto de vista del algebra vectorial
Desde el punto de vista del algebra vectorial, las magnitudes físicas son representadas
convenientemente usando la terminología de sumatorias en el caso vectorial y tensorial, para
facilitar las operaciones entre sus vectores unitarios.
Escalar: s , x , t , v , etc.
3 3 3

Vector: V =   iVi ; F = F =   k Fk ;  =  j
i =1 k =1 j =1 x j
3 3 3 3

Tensor:  =  i j ij
i =1 j =1
; V =   i
i =1 j =1 xi
 jV j

Operaciones con vectores unitarios

Producto escalar de dos vectores unitarios

1 si i = j 
 i . j =  ij  
0 si i  j 

 ij es el operador delta Kronecker y es un escalar que solo toma dos valores posibles.

Producto vectorial de dos vectores

 i   j =   ijk  k
k
Figura 1: Esquema para la determinación del signo del tensor de alternancia

+1 si ijk = 123, 231,312 

 
 ijk −1 si ijk = 132,321, 213 
0 si dos indices cualquiera son iguales 
 

 ijk es el tensor de alternancia o tensor de Levi-Civita; este también puede ser obtenido
mediante:  ijk = (1/ 2)(i − j )( j − k )(k − i)

Para realizar algunas demostraciones, las siguientes igualdades resultan ser útiles
3 3 3

  ijk hjk = 2ih y

j =1 k =1

k =1

ijk mnk =  im jn −  in jm

Estas ecuaciones también pueden ser aplicadas al producto punto y producto aspa a vectores
unitarios en los diferentes sistemas de coordenadas; a manera de ejemplo, se muestran algunos
productos.

 i . i = 1  i . j = 0 i  i = 0 i   j =  k  j   i = − k

En el producto aspa se considera positivo si va en dirección horaria y negativo en dirección

antihoraria.

Resumen de las operaciones con vectores unitarios

Al aplicar los diferentes productos a los vectores y tensores, estas operaciones se aplican
realmente a sus vectores unitarios; por tanto, es necesario conocer el resultado de dichas
operaciones, cuyo resumen se muestra. Con este conjunto de operaciones sobre los vectores
unitarios en coordenadas rectangulares y sus equivalente en otros sistemas de coordenadas, se
pueden realizar todas las operaciones requeridas utilizadas en el curso de fenómenos de
transporte.

1 si i = j 
 i . j =  ij  
0 si i  j 

 i   j =   ijk  k
k

i j :  kl = ( j . k )(i . l ) =  jkil

 i j . k  =  i ( j . k ) =  i jk

 i . j k  = ( i . j ) k =  ij k

Nota: Para operaciones, cuya forma primaria predice un resultado vectorial, se utiliza los
corchetes; para escalares no se usa ningún símbolo de colección o paréntesis; para el caso de
tensores se utiliza la llave.

  .   =  ( .
i j k l i j k ) l  =  jk  i l 

 
 
i j   l

  k  =  i  j   k  =   jkl i l 


 
     =   
i j k i j 
  k =   ijl l k 
 l 
Ejemplo 5
Realice el producto escalar de dos vectores, utilizando los conceptos del algebra vectorial.
Solución:

Escribimos los vectores V y W luego aplicamos el producto escalar correspondiente, haciendo

uso del producto escalar de vectores unitarios, cuyo resultado se expresa en términos del operador
delta Kronecker.
3 3 3 3 3 3
V .W =   iVi .  jW j = VW
i j i . j = VW
i j ij
i =1 j =1 i =1 j =1 i =1 j =1

Ahora aplicando el operador delta Kronecker; i = j se tiene

3
V .W =  VW
i i = V1W1 + V2W2 + V3W3
i =1

El resultado es el típico valor escalar obtenido, expresado en la suma de tres términos, cuyo
resultado se obtiene también fácilmente por los métodos clásicos aplicado a este ejemplo.
Ejemplo 6
Supongamos que se desea calcular el trabajo requerido para mover un cuerpo dese el punto
(1,0,0) al punto (2,1,1). Si la fuerza F está dada por:

32 32
F= 2 − 3
2 2
Encuentre dicho trabajo.
Solución:

Por definición se tiene: dW = F .ds . Ahora integramos miembro a miembro esta ecuación con
lo que se obtiene:
 Q

W =  F .ds
P

Obtenemos el vector de desplazamiento

P = (1,0,0) Q = (2,1,1) , luego se tiene:

ds = (2 − 1)dx x + (1 − 0)dy y + (1 − 0)dz z =  x dx +  y dy +  z dz

Ahora realizamos el producto escalar

32 32 32 32
F .ds = ( y −  z ).( x dx +  y dy +  z dz ) = ( dy − dz)
2 2 2 2
Luego:
Q y =1 z =1
32 32 32 32 32 32
W = ( dy − dz ) =  dy −  dz = − =0
P 2 2 y =0 2 z =0 2 2 2

Ejemplo 7
Realice el producto vectorial de dos vectores, utilizando los conceptos del algebra vectorial.
Solución:

Escribimos nuevamente los vectores V y W y aplicamos el producto vectorial usando la regla

correspondiente descrita anteriormente para el producto vectorial de vectores, del siguiente modo:
3 3 3 3 3
V W =   iVi    jW j = ViW j  ijk  k
i =1 j =1 i =1 j =1 k =1

El resultado final puede final puede quedar como tal o expandirlo para cada i, j y k, utilizando la
regla de la multiplicación expuesta, de tal manera que se obtiene los mismos resultados reportados
en la literatura para el producto vectorial de vectores dado por:

i j k
V  W =  V1 V2 V3  = (V2W3 − W2V3 ) i + (VW

1 3 ) j + (VW
1 3 − WV 1 2 ) k
1 2 − WV

W1 W2 W3 

Que es equivalente a:

V  W = (V2W3 − W2V3 )1 + (VW 1 3 ) 2 + (VW

1 3 − WV 1 2 ) 3
1 2 − WV

Ejemplo 8

Realice el producto escalar de dos tensores  y V

Solución:

Escribimos el producto doble punto entre los tensores  y V , luego aplicamos las reglas de
multiplicación antes expuestas del siguiente modo:
  
 : V  =  ij i j :   kVk   l =   i j :  k  l ijVk
i j k l xl i j k l xl

 1 si i = l y j = k 
 V  il jk  
xl 0 si i  l y j  k 
ij k
ijkl

Luego, usando esta igualdad se tiene


 : V  =  ijV j
ij xi

Esta ecuación puede ser expandida para cada i y j, según los requerimientos al aplicar a algún
problema específico.
Ejemplo 9

Realice el producto diádico de dos vectores V y W

Solución:
Dado que se trata de un producto simple de vectores (no hay ninguna operación vectorial entre
ambos vectores), el producto solo queda indicado, de tal forma que también puede ser expandido
para cada i y j, lo cual permite obtener nueve componentes.
3 3 3 3
VW =   iVi   jW j = VW
i j  i j
i =1 j =1 i =1 j =1

Ejemplo 10
Usando los criterios del algebra vectorial, demuestre la siguiente identidad:

U  V  W  = V (U .W ) − W (U .V )

Solución:
El termino i del miembro de la izquierda se puede expresar mediante:

 
U  V  W  i =   ijkU j V  W k  i =   ijkU j i   klmVW
l m
j k j k  l m 

 
j k l m
ijk klm l m =   ijk  lmkU jVW
 U jVW
j k l m
l m

U  V  W  i =  ( il jm −  im jl )U jVlWm = Vi   jmU jWm − Wi   jlU jVl
j l m j m j l

Vi U jW j − Wi U jV j = Vi (U .W ) − Wi (U .V )
j j

Es la i-estima componente del miembro de la derecha de la expresión; que, sin embargo, se

puede sumar todos los términos y obtener el valor solicitado.
7.2. LONGITUDES, AREAS Y VOLUMENES EN LOS DIFERENTES SISTEMAS
DE COORDENADAS
Sistema de coordenadas rectangulares
La longitud diferencial en las direcciones x, y, z desde el origen al punto ( x, y, z ) , está dada
por:

ds1 = dx, ds2 = dy, ds3 = dz (7.1)

El área diferencial de cada uno de los lados del cubo unitario se obtiene mediante

dAx = dydz, dAy = dxdz, dAz = dxdy (7.2)

El volumen diferencial del cubo unitario está dado por:

dv = dxdydz (7.3)

Ejemplo 11
Calcule el área y el volumen de un paralelepípedo limitado por los siguientes dominios:

0  x  1; 0  y  2; 0  z  4
Solución:
Las areas del paralelepípedo, limitado por los dominios dados en la propuesta del problema,
pueden ser calculadas mediante:
4 2 4 1
Ax =   dydz = y 0 z 0 = 2  4 = 8u 2 Ay =   dxdz = x 0 z 0 = 4 1 = 4u 2
2 4 1 4

0 0 0 0

2 1
Az =   dxdy = x 0 y 0 = 1 2 = 2u 2
1 2

0 0

Del mismo modo, el volumen se obtiene mediante:

4 2 1
v =    dxdydz = x 0 y 0 z 0 = 1 2  4 = 8u 3
1 2 4

0 0 0

Generalizando
z2 y2

Sobre la superficie x = x0 A( x0 , y, z ) =   f (x0 , y, z )dydz

z1 y1

z2 x2

Sobre la superficie y = y0 A( x, y0 , z ) =   f (x, y0 , z )dxdz

z1 x1

y2 x2

Sobre la superficie z = z0 A( x, y, z0 ) =   f (x, y, z )dxdy

y1 x1
0

Sistema de coordenadas polares

Las coordenadas polares de un punto P del plano son la distancia r al origen y el ángulo  desde
el semieje x positivo a la dirección OP en sentido antihorario, tal como se muestra en la figura
2.

Figura 2: Sistema de coordenadas polares en el plano

Si el punto P tiene coordenadas ( x, y ) , con x  0 , entonces, sus coordenadas polares en (r ,  )

se obtiene como

r = x2 + y 2 (7.4)

 = arctan( y / x) = tan −1 ( y / x) (7.5)

La transformación inversa de coordenadas polares a rectangulares es:

x = r cos  (7.6)

y = r sin  (7.7)

Mientras que x e y puede tomar cualquier valor real  tiene limitaciones; para determinar la
coordenada angular  , se deben distinguir dos casos:

Para r = 0 , el ángulo  puede tomar cualquier valor real

Para r  0 , para obtener un único valor de  debe limitarse a un intervalo de tamaño 2 .Por
convención, los intervalos utilizados son  0, 2  y  − ,  

Para obtener  en el intervalo  0, 2  , se deben usar las siguientes formulas

 3
arctan( y / x) si x  0 y y  0 ; si x = 0 y y  0 ; arctan( y / x) +  si x  0 ; si
2 2
x = 0 y y  0 ; arctan( y / x) + 2 si x  0 y y  0

  
Para obtener  en el intervalo  − ,   , se considera que arctan( y / x)   − ,  es una
 2 2
función creciente en su dominio

 
arctan( y / x) −  sí x  0 y y  0 ; − si x = 0 y y  0 ; arctan( y / x) si x  0 ; si
2 2
x = 0 y y  0 ; arctan( y / x) +  si x  0 y y  0
Sistema de coordenadas cilíndricas
Para el desarrollo de problemas en coordenadas cilíndricas, es importante tener en cuenta los ejes
de referencia tanto el eje cartesiano como el eje cilíndrico, esto tiene importancia para la ubicación
de los planos respectivos cuyo valor es constante. En la figura 3 se muestra la correspondencia
respectiva para ser usada en las transformaciones de sistemas de coordenadas.

Figura 3: Relacion entre el sistema de coordenadas rectangulares y cilindricas

Figura 4: Transformación de un sólido en coordenadas rectangulares a coordenadas cilíndricas.

Equivalencia entre coordenadas cilíndricas y rectangulares

El eje de coordenadas está dado por:

x = r cos  y = rsen z=z (7.8)

Con r  0 y 0    2 , −  z  

r = + x2 + y 2 (7.9)

 = arctan( y / x) (7.10)
z=z (7.11)
Las diferenciales de longitud, en este caso están dadas por:

ds1 = dr ds2 = rd ds3 = dz (7.12)

La diferencial del área lateral está dada por:

dAL = (rd )(dz ) = rd dz (7.13)

La diferencial del área de la base viene dada por:

dAB = (rdr )(d ) = rdrd (7.14)

La diferencial de un plano vertical paralelo a la dirección longitudinal inscrito dentro del

cilindro es:

dAV = drdz (7.15)

La diferencial de volumen

dv = (dr )(rd )(dz ) = rdrd dz (7.16)

Ejemplo 12
Un cilindro tiene un radio R y una longitud L. Considerando las ecuaciones anteriores,
encuentre:
a) El área lateral
b) El área de la base
c) El volumen
Solución:
a) Para encontrar el área lateral integramos la ecuación (7.13)
2 L 2 L
AL = 
0 0
 dAL =   rd dz = 2 RL
0 0

b) Para encontrar el área de la base integramos la ecuación (7.14)

2 R 2 R
R2
AB = 
0 0
 dAB =   rdrd = 2
0 0
2
=  R2

En forma general
z2  2

Sobre la superficie r = r0 A(r0 , , z ) =   f (r0 ,  , z )r0 d dz

z1 1

z2 r2

Sobre la superficie  =  0 A(r ,  0 , z ) =   f (r ,  0 , z )drdz

z1 r1

 2 r2
Sobre la superficie z = z0 A(r ,  , z0 ) =   f (r, , z )rdrd
1 r1
0
c) Para encontrar el volumen integramos la ecuación (7.15)
2 R L 2 R L
R2
v=    dv =
0 0 0
   rdrd dz = 2
0 0 0
2
L =  R2 L

z2  2 r2

Generalizando: v =    f (r , , z )rdrd dz
z1 1 r1

Sistema de coordenadas esféricas

Al igual que en el caso de las coordenadas cilíndricas, se debe tener bastante cuidado en el uso de
la nomenclatura adecuada; cada autor utiliza distinta nomenclatura referida principalmente a la
ubicación de los ángulos  y  .En este trabajo se usa  en lugar de  , siguiendo la
nomenclatura desarrollada por Bird et al. (1996). En la figura 4 se muestra el esquema para la
transformación de sistema de coordenadas cartesianas a esféricas y viceversa.

Figura 5: Relación entre los sistemas de coordenadas rectangulares y esféricas

 Ojo: En esta figura 4 cambiar tetha por phi

Figura 6: Transformación de un sólido en coordenadas rectangulares a coordenadas esféricas.

Las equivalencia entre coordenadas esféricas y rectangulares se muestran en las ecuaciones

subsiguientes.

x = rsen cos  y = rsen sen z = r cos (7.17)

Con r  0 , 0     , 0    2

r = + x2 + y 2 + z 2 (7.18)

x2 + y2
 = arctan (7.19)
z
y
 = arctan (7.20)
x
Las diferenciales de longitud, en este caso están dadas por:

ds1 = dr ds2 = rd ds3 = rsen d (7.21)

La diferencial del área lateral está dada por:

dA = (rd )(rsen d ) = r 2 sen d d (7.22)

La diferencial de volumen

dv = (dr )(rd )(rsen d ) = r 2 sen drd d (7.23)

Ejemplo 13
Obtenga el área y el volumen de una esfera de radio R:
Solución:
Para encontrar el área de la esfera integramos la ecuación (7.22)
 2  2  
A=  
 =0  =0
dA =  
 =0  =0
r 2 sen d d = 2 R 2  sen d = 4 R 2
0

Generalizando
2  2
Sobre la superficie r = r0 A(r0 , ,  ) =   f (r , , )r sen d d
2
0 0
 1 1

2 r2
Sobre la superficie  =  0 A(r , 0 ,  ) =   f (r, , )sen rdrd
1 r1
0 0

2 r2
Sobre la superficie  = 0 A(r , , 0 ) =   f (r, ,  )rdrd
0
1 r1

Para encontrar el volumen de la esfera integramos la ecuación (1.23)

2  R 2  R
/ 3 cos  0  0 = −( R3 / 3)(cos(180) − cos(0))  2
 2
   dv =    r sen drd d = − r
R
v= 2 3
0
0 0 0 0 0 0

R3 4
v=− (−1 − 1)  2 =  R3
3 3
2  2 r2
Generalizando: v =    f (r , ,  )r sen drd d
2

1 1 r1

Sistema de coordenadas ortogonales generales

A continuación, se presenta el esquema de un sistema de coordenadas generalizado que permite
realizar operaciones vectoriales, usando las equivalencias respectivas; para lo cual se hace uso del
esquema de la figura 7.

Figura 7: Esquema del sistema de coordenadas ortogonales generalizado

Introduzcamos los coeficientes de transformación o factores escalares h1 , h2 y h3 , con las
siguientes igualdades

d = h1ds1 , d  = h2 ds2 y d  = h3ds3

donde  ,  y  son las componentes del sistema de coordenadas ortogonales generalizado y

s1 , s2 y s3 son las longitudes respectivas.

En el sistema de coordenadas rectangulares, los coeficientes de transformación están dados por:

h1 = 1, h2 = 1 y h3 = 1

En el sistema de coordenadas cilíndricas

1
h1 = 1, h2 = y h3 = 1
r
En el sistema de coordenadas esféricas

1 1
h1 = 1, h2 = y h3 =
r rsen
En la tabla 2 se muestran las equivalencias en los diferentes sistemas de coordenadas para su
transformación y aplicación.
Tabla 2: Equivalencias de los diferentes parámetros en el sistema de coordenadas ortogonales
generalizado y los sistemas de coordenadas clásicas.

Geometria Rectangular Cilíndrica Esférica

 x r r
 y  
 z z 
d dx dr dr
d dy d d
d dz dz d
ds1 dx dr dr
ds2 dy rd rd
ds3 dz dz rsen d
h1 1 1 1
h2 1 1/ r 1/ r
h3 1 1 1/ rsen

Aplicaciones del sistema de coordenadas ortogonales generalizado

Gradiente de un escalar

s s s
s = 1h1 +  2 h2 +  3h3 (7.24)
  
Ejemplo 14
Haciendo uso de la ecuación (7.24) y las equivalencias de la tabla 1 encuentre s en los tres
sistemas de coordenadas.
Solución:
En la ecuación (1.24) se reemplaza los valores respectivos de sus parámetros, haciendo uso de la
tabla 2
a) Coordenadas rectangulares
s s s
s =  x +  y + z
x y z
b) Coordenadas cilíndricas
s 1 s s
s =  r +  + z
r r  z
c) Coordenadas esféricas
s 1 s 1 s
s =  r +  + 
r r  rsen 
Divergencia de un vector

   V    V2    V3 
divV = .V = h1h2 h3   1  +  +   (7.25)

  2 3
h h   1 3  
h h  h1h2 
Ejemplo 16

Haciendo uso de la ecuación (7.25) y sus equivalencias dadas en la tabla 2 encuentre .V en los
tres sistemas de coordenadas.
Solución:
En la ecuación (7.25) se reemplaza los valores respectivos de sus parámetros, haciendo uso de la
tabla 2
a) Coordenadas rectangulares
Vx Vy Vz
.V = + +
x y z
Es la ecuación de continuidad para un fluido compresible

b) Coordenadas cilíndricas
1 (rVr ) 1 V Vz
.V = + +
r r r  z
c) Coordenadas esféricas
1 (r 2Vr ) 1 (V sen ) 1 V
.V = + +
r 2
r rsen  rsen 
Rotacional de un vector

   V    V     V    V3     V2    V1   (7.26)
 V = 1h3h2   3  −  2  +  2 h1h3   1  −   +  3h2 h1    −  
   h3    h2     h1    h3      h2    h1  

Ejemplo 16
Haciendo uso de la ecuación (7.26) y sus equivalencias dadas en la tabla 2 encuentre  V en
los tres sistemas de coordenadas
Solución:
a) Coordenadas rectangulares
Vz Vy V V V V
V =  x ( − ) +  y ( x − z ) + z ( y − x )
y z z y x y
b) Coordenadas cilíndricas
1 Vz Vz V V 1 (rV ) 1 Vr
 V =  r ( − ) +  ( r − z ) +  z ( − )
r  z z r r r r 
c) Coordenadas esféricas
1 (V sen ) 1 V 1 Vr 1 (rV )
 V =  r ( − ) +  ( − )+
rsen  rsen  rsen  r r
1 (rV ) 1 Vr
 ( − )
r r r 
Laplaciano de un escalar

   h s    h2 s    h3 s 
2 s = h1h2 h3   1 +  +   (7.27)

  2 3
h h     1 3
h h     h1h2  
Ejemplo 17

Haciendo uso de la ecuación (7.27) y las equivalencias de la tabla 1 encuentre  s para los tres
2

sistemas de coordenadas
Solución:
a) Coordenadas rectangulares

2s 2s 2s

2 s = + +
x 2 y 2 z 2
b) Coordenadas cilíndricas
1  s 1 2s 2s
2 s = (r ) + 2 +
r r r r  2 z 2
c) Coordenadas esféricas
1  2 s 1  s 1 2s
2 s = ( r ) + ( sen ) +
r 2 r r r 2 sen   r 2 sen  2
Nota: Para hallar el Laplaciano de un vector y otras operaciones como la divergencia de un tensor,
por lo general se obtiene en coordenadas rectangulares y luego se aplica las transformaciones
respectivas.
Elemento de línea
2 2 2
 d   d    d  
ds = 
2
 +  +  (7.28)
 1   2   h3 
h h
Ejemplo 18
Haciendo uso de la ecuación (7.28) y las equivalencias de la tabla 2 obtenga los elementos de
línea en los tres sistemas de coordenadas.
Solución:
Se sustituye los parámetros respectivos de la tabla 2 en la ecuación (7.28) y se obtiene los
siguientes resultados.
a) Sistema de coordenadas rectangulares
dx 2 dy 2 dz 2
ds 2 = ( ) + ( ) + ( ) = (dx)2 + (dy)2 + (dz) 2
1 1 1
b) Coordenadas cilíndricas
dr 2 d 2 dz 2
ds 2 = ( ) +( ) + ( ) = (dr )2 + (rd )2 + (dz)2
1 1/ r 1
c) Coordenadas esféricas
dr 2 d 2 d
ds 2 = ( ) +( ) +( )2 = (dr )2 + (rd )2 + (rsen d )2
1 1/ r 1/ rsen
Elemento de área diferencial

d  d
dA = (7.29)
h2 h3

Ejemplo 19
Obtenga el área diferencial del plano xy en el sistema de coordenadas rectangulares

Solución:
Se hace uso de la ecuación (7.29) y los parámetros de la tabla 2 y se obtiene:

dy dx
dA = = dydx
1 1
Ejemplo 20
Obtenga el diferencial del área lateral de un cilindro
Solución:
Se hace uso de la ecuación (7.29) y los parámetros de la tabla 2 y se obtiene:

d dz
dAL = = rd dz
1/ r 1
En el caso de las coordenadas cilíndricas se puede modificar la ecuación anterior dependiendo
que es lo que se desea calcular. Para e área lateral se usa la formula tal cual; sin embargo, para el
área de la base se tendría que emplear la siguiente expresión.

d d 
dA =
h1 h2

Ejemplo 21
Obtenga el área diferencial de la base de un cilindro
Solución:
El área de la base de un cilindro esta dado por

dr d
dAB = = rdrd
1 1/ r
Ejemplo 22
Obtenga el área diferencial de una esfera de radio r
Solución:
Usamos la ecuación (7.29) y los parámetros de la tabla 2 y se obtiene lo siguiente:

d d
dA = = r 2 sen d d
1/ r 1/ rsen
Elemento de volumen diferencial

d d  d 
dv = (7.30)
h1 h2 h3

Ejemplo 23
Encuentre el volumen diferencial de un cilindro de radio r y longitud z
Solución:
A partir de la ecuación (7.30) y haciendo uso de los parámetros de la tabla 2 se tiene:

dr d dz
dv = = rdrd dz
1 1/ r 1
Ejemplo 24
Encuentre la diferencial de volumen de una esfera de radio r
Solución:
A partir de la ecuación (7.30) y los parámetros de la Tabla 30 se tiene:

dr d d
dv = = r 2 sen drd d
1 1/ r 1/ rsen
7.3 ALGEBRA DE VECTORES Y TENSORES EN COORDENADAS
CURVILINEAS
Previamente se debe conocer la aplicación de la regla de la cadena para hacer las transformaciones
de coordenadas rectangulares a coordenadas curvilíneas.

 F   r            z    
 x  =  x   r  +  x     +  x   z  (7.31)
y,z y,z  ,z y,z r,z y,z r ,
Transformación de coordenadas rectangulares a cilíndricas
x = r cos  y = rsen z=z (7.32)

r = + x2 + y 2  = arctan( y / x) z=z (7.33)

Al aplicar la regla de la cadena, se obtiene

  sen  
= (cos  ) + (− ) + (0) (7.34)
x r r  z
  cos   
= ( sen ) + ( ) + (0) (7.35)
y r r  z

   
= (0) + (0) + (1) (7.36)
z r  z

Figura 7: Relación de vectores unitarios en coordenadas rectangulares y cilíndricas. En este

esquema no se presenta el vector unitario en la dirección z, para tener un efecto visual más
adecuado.

Relación entre vectores unitarios en coordenadas cilíndricas

De coordenadas rectangulares a coordenadas cilíndricas

 r = (cos ) x + (sen ) y + (0) z (7.37)

 = (−sen ) x + (cos ) y + (0) z (7.38)
 z = (0) x + (0) y + (1) z (7.39)

De coordenadas cilíndricas a coordenadas rectangulares

 x = (cos  ) r + (− sen ) + (0) z (7.40)

 y = (sen ) r + (cos ) + (0) z (7.41)
 z = (0) r + (0) + (1) z (7.42)
Ejemplo 25
Obtener  en coordenadas cilíndricas

Solución:
La gradiente en coordenadas rectangulares está dada por:

  
 = x +  y + z (a)
x y z
Las derivadas parciales se expresan mediante:

  sen  
= (cos  ) + (− ) + (0) y (b)
x r r  z
  cos   
= ( sen ) + ( ) + (0) (c)
y r r  z
   
= (0) + (0) + (1) (d)
z r  z
Por otro lado, la relación entre vectores unitarios esta dada por:

 x = (cos  ) r + (− sen ) + (0) z (e)

 y = (sen ) r + (cos ) + (0) z (f)
 z = (0) r + (0) + (1) z (g)

Reemplazando (b), (c), (d) y (f) en (a) se tiene la siguiente expresión:

 sen 
 = ( cos  r − sen ) (cos  − ) + ( sen r + cos  )
r r 
 cos   
(sen + ) +z
r r  z
Desarrollando la multiplicación termino a término y simplificando se tiene

 1  
 = r +  + z
r r  z
Transformación de coordenadas rectangulares a esféricas
x = rsen cos  y = rsen sen z = r cos  (7.43)

x2 + y2 y
r = + x2 + y 2 + z 2  = arctan  = arctan (7.44)
z x
A partir de la regla de la cadena se obtiene:

  cos  cos   sen 

= ( sen cos  ) + ( ) + (− ) (7.45)
x r r  rsen 
  cos  sen  cos  
= ( sen sen ) + ( ) +( ) (7.46)
y r r  rsen 

  sen  
= (cos  ) + (− ) + (0) (7.47)
z r r  

Relación entre vectores unitarios en coordenadas esféricas

De coordenadas rectangulares a coordenadas esféricas

 r = (sen cos  ) x + (sen sen ) y + (cos ) z (7.48)

 = (cos cos  ) x + (cos sen ) y + (−sen ) z (7.49)

 = (−sen ) x + (cos  ) y + (0) z (7.50)

De coordenadas esféricas a coordenadas rectangulares

 x = (sen cos  ) r + (cos cos  ) + (−sen ) (7.51)

 y = (sen sen ) r + (cos sen ) + (cos  ) (7.52)

 z = (cos ) r + (−sen ) + (0) (7.53)

Siguiendo el mismo procedimiento anterior, para coordenadas esféricas se obtiene:

 1  1 
 = r +  + 
r r  rsen 

7.4. OPERACIONES DIFERENCIALES EN COORDENADAS CURVILINEAS

Derivada de vectores unitarios en coordenadas cilíndricas
 r   z
=0 =0 =0 (7.54)
r r r
 r   z
=  = − r =0 (7.55)
  
 r   z
=0 =0 =0 (7.56)
z z z
4.2 Derivada de vectores unitarios en coordenadas esféricas

 r  
=0 =0 =0 (7.57)
r r r

 r  
=  = − r =0 (7.58)
  
 r  
= sen = cos  = − sen r − cos  (7.59)
  

Ejemplo 26

Encontrar .V en coordenadas cilíndricas

Solución:
Utilizamos las equivalencias del vector gradiente y del vector velocidad

  1  
(.V ) =  r +  +  z  . rVr +  V +  zVz 
 r r  z 

Multiplicando miembro a miembro el producto escalar

 rVr  V  V
(.V ) =  r . + r .   + r . z z +
r r r
1  rVr 1 V 1  zVz  V  V  V
+ . +  . +  . +z. r r +z.   +z. z z
r  r  r  z z z
Luego se lleva a cabo cada una de las derivadas de los productos, es decir:

 rVr  V  V V V
r . =  r .(Vr r +  r r ) =  r .Vr r +  r . r r =  r .Vr (0) + (1) r = r
r r r r r r r
V  V  V V
r . =  r .(V  +   ) =  r .V  +  r .  =  r .V (0) + (0)  = 0
r r r r r r
 zVz  V  V V
r . =  r .(Vz z +  z z ) =  r .Vz z +  r . z z =  r .Vz (0) + (0) z = 0
r r r r r r
1  rVr 1  V 1  1 Vr 1 1 Vr 1 1 Vr 1
 . =  .( (Vr r +  r r )) =  .( Vr r ) +  . r =  .( Vr ( )) + (0) =  . ( Vr + (0) = Vr
r  r   r  r  r r  r r  r

1 V 1  V 1  1 V 1 1 V 1 V

 . =  .( (V  +  r )) =  .( Vr  ) +  . =  .(− r )( Vr ) + (1) =
r  r   r  r  r r  r 

1  zVz 1  V 1  1 V 1 1 V

 . =  .( (Vz z +  z )) =  .( Vz z ) +  . z =  .( Vz (0)) + (0) =0
r  r   r  r  r r 
 rVz  V V V
z. =  z .(Vr r +  r r ) =  z .(Vr (0)) +  z . r r =  z .(Vr (0)) + (0) r = 0
z z z z z
V  V V V
z. =  z .(V  +   ) =  z .(V (0)) +  z . 0 =  z .(V (0)) + (0)  = 0
z z z z z
 zVz  V V V V
z. =  z .(Vz z +  z z ) =  z .(V (0)) +  z . z z =  z .(Vz (0)) + (1) z = z
z z z z z z
Sumando todos los términos se obtiene
 Vr 1 V Vr Vz
.V = + + +
r r  r z
Siguiendo este procedimiento se puede encontrar el valor de .V e coordenadas esféricas, cuyo
resultado es el siguiente:

1 (r 2Vr ) 1  1 V

.V = + (V sen ) +
r 2
r rsen  rsen 

Ejemplo 27

Encontrar V en coordenadas cilíndricas

Solución:
La gradiente en coordenadas cilíndricas está dada por:

 1  
 = r +  + z
r r  z
Y el vector velocidad

V =  rVr + V +  zVz

Luego aplicamos el producto simple:

 1  
V = ( r +  +  z )( rVr + V +  zVz )
r r  z
Multiplicando cada uno de los términos se tiene:

   1  1 
V =  r  rVr +  r V +  r  zVz +   rVr +  V +
r r r r  r 
1    
  zVz +  z  rV +  z V +  z  zVz
r  z z z
Resolviendo los productos y las derivadas correspondientes

  V V V
r  rVr =  r (Vr  r +  r r ) =  r (Vr (0) +  r r ) =  r r r
r r r r r
  V V V
r V =  r (V  +   ) =  r (V (0) +   ) =  r 
r r r r r
  V V V
r  zVz =  r (Vz  z +  z z ) =  r (Vz (0) +  z z ) =  r z z
r r r r r
1  1  rVr 1  V 1 V V 1 Vr
  rVr =  +  (Vr r +  r r ) =  (Vr +  r r ) =   r +   r
r  r  r   r  r r 
1  1 V 1  V 1 V V 1 V
 V =  +  (V  +   ) =  (V (− r ) +   ) = −  r  +  
r  r  r   r  r r 
 1  1  zVz 1  V 1 V 1 Vz
  zVz =  +  (Vz z +  z z ) =  (Vz (0) +  z z ) =   z
r  r  r   r  r 
  V V V
z  rVr =  z (Vr r +  r r ) =  z (Vr (0) +  r r ) =  z r r
z z z z z
  V  V V
z V =  z (V  +   ) =  z (Vz  +   ) =  z 
z z z z z z
  V V V
z  zVz =  z (Vz  +  z z ) =  z (Vz (0) +  z z ) =  z z z
z z z z z
Sumando el resultado obtenido:
Vr V V 1 Vr V 1 V Vr 1 Vz V V V
V =  r r +  r   +  r  z z +    r ( − ) +   ( + ) +   z +  z r r +  z  +  z z z
r r r r  r r  r r  z z z

Siguiendo el mismo procedimiento, se obtiene V en coordenadas esféricas

Ejemplo 28

Encuentre .  en coordenadas esféricas

r

Solución:
La gradiente en coordenadas esféricas está dada por:

 1  1 
 = r +  + 
r r  rsen 
El Tensor de esfuerzos esta dado por

 rr r r  r  r  r r 

 
 =   r  r        
          
  r  r   
El producto total se obtendría multiplicando escalarmente el vector gradiente por el tensor de
esfuerzos.

 rr r r  r  r  r r 

  1  1    
.  =  r +  +  .   r  r        
 r r  rsen   
  r r       

De los 27 términos que se obtiene al multiplicar la gradiente por el tensor, se debe considerar solo
aquellos términos que involucre la dirección radial.
El resultado del primer producto permite obtener lo siguiente

 (  )    
r . r r rr =  r .( r r  rr +  r r rr ) =  r .(( r  r +  r r ) rr +  r r rr ))
r r r r r r
  rr  
=  r .(((0) r +  r (0)) rr +  r r )) =  r . r r rr =  r rr
r r r
Cabe mencionar en este tipo de demostraciones que se debe aplicar estrictamente la regla de la
derivación de un producto, es decir: d (uv) = uv + uv . Desde el punto de vista matemático no
tendría importancia este detalle, es decir la derivada también se puede escribir como:
d (uv) = vu + uv ; sin embargo, al momento de operar tensores es importante el orden de
aparición de los vectores unitarios y como consecuencia de ello también el resultado será otro
valor no esperado.

1  1    1   
 .  r r =  .(  r  r +   r  r ) =  .((   r +  r )  r +   r  r )
r  r   r   
1  1 
=  .((− r ) r +   ) r +   r  r ) = ((− . r r +  .  ) r +  .  r  r )
r  r 
1  1  r 1
= (0 +  r +  r  r ) =  r +    r
r  r  r
De acuerdo con Bird et al (2006), dado que solo estamos considerando la componente r el termino
de color verde corresponde a “otros términos”; en este caso se utiliza para la componente  .

Continuando con el desarrollo

1  1  1   
 .  r  r =  .  r  r =  .(  r   r +  r  r )
rsen  rsen  rsen  

1    1 
=  .((   r +  r )  r +  r  r )) =  .((−sen r r − cos   r +  sen )  r +  r  r ))
rsen    rsen 

1 
 .((−sen . r r − cos  .  r +  . sen )  r +  . r  r )) =
rsen 

1  1   r 1
= ((0 − 0 + (1) sen )  r + (1) r  r )) =  r +    r
rsen  rsen  r
1          
 . r r rr = rr  .( r r ) = rr  .( r  r +  r r ) = rr  .(  r +  r ) = rr ( r + 0) =  r rr
r  r  r   r r r
1        
 .  r r rr = rr  . r r = rr  .( r  r +  r r ) = rr  .( sen r +  r sen )
rsen  rsen  rsen   rsen

 rr  rr  rr
= (sen . r +  . r sen ) = (sen r + (0)sen ) =  r
rsen rsen r
1  1  1    1   
 .  =  . .  =  .(     +    ) =  .((   +   )  +    )
r  r  r   r   

1  1   1 
=  .((− r ) −   r ) +    ) = ((0 −  r ) +   )) = − r  + 
r  r  r r 
 1      
 .   r  r =  r  .   r =  r  .(   r +  r )
rsen  rsen  rsen  

 r  r   r cos
=  .(cos r +   sen ) = (cos r + (0)sen ) =  r
rsen rsen r sen
En este último término tampoco aparece el vector unitario en la dirección r ; sin embargo, las
operaciones entre los vectores unitarios permiten la aparición de las componentes asociadas al
radio.

1  1  1   
 .   =  .(   ) =  .(     +   )
rsen  rsen  rsen  
1    1 
=  .((   +   )  +   ) =  .((− sen r − cos   +  (−sen r − cos  )  +    )
rsen    rsen 

1 
= ((−sen . r − cos  .  +  . (−sen r − cos  )  +  .  )
rsen 

1  1 1 cos  1  
= ((0 − 0 + (− sen r − cos  )  +   ) = − r   −    + 
rsen  r r sen rsen 

Sumando los términos respectivos se tiene:

1   1   r 1   r   +  
.  = 2 (r 2 rr ) +  r cot  + + −
r r r r r  rsen  r
Otras ecuaciones de interés
Cuando se produce disipación viscosa (generación de calor por rozamiento de las moléculas
expuestas a altas velocidades, aparece este término de rozamiento que contribuye a la tasa de
generación de calor. Este producto doble punto de los tensores  y V es una magnitud escalar
y se representa en los tres sistemas de coordenadas.
Coordenadas rectangulares

Vx V V V V V
( : V ) =  xx ( ) +  xy ( x ) +  xz ( x ) +  yx ( y ) +  yy ( y ) +  yz ( y )
x y z x y z
V V V
 zx ( z ) +  zy ( z ) +  zz ( z )
x y z
Coordenadas cilíndricas

Vr 1 Vr V V V 1 V Vr V

( : V ) =  rr ( ) +  r ( − ) +  rz ( r ) +   r (  ) +   ( + ) +  z (  ) +
r r  r z r r  r z

Vz 1 Vz V
 zr ( ) +  z ( ) +  zz ( z )
r r  z
Coordenadas esféricas
 Vr 1 Vr V 1 Vr V
( : V ) =  rr ( ) +  r ( − ) +  r ( − )+
r r  r rsen  r

V 1 V Vr 1 V V

 r ( ) +   ( + ) +   ( − cot  ) +
r r  r rsen  r

V 1 V 1 V Vr V

 r ( ) +   ( ) +   ( + + cot  )
r r  rsen  r r

7.5. TEOREMAS INTEGRALES PARA VECTORES Y TENSORES

Teorema de Green

Si F es un campo vectorial continuo y diferenciable en R 2 , la integral de circulación de F en

sentido contrario a las agujas del reloj a lo largo de una curva simple cerrada C en el plano, es
igual a la integral de superficie de la componente k del rotacional del campo vectorial sobre la
región S encerrada en la curva C

Figura 8.- Esquema del campo vectorial para la aplicación del Teorema de Green

N M
 C
F .dr =   (
S
−
x y
)dA (7.60)

La integral de la izquierda representa la circulación del campo vectorial a lo largo de la curva C,

mientras que el término de la derecha representa la superficie S de la componente k del rotacional;
2
dado que estamos analizando en R , solo depende de x e y ; es decir que F = f ( M , N ) . El
rotacional esta dado por:

 i j k 
 
    N M
 F =  =( − ) k
 x y z  x y
 
 M ( x , y ) N ( x, y ) 0 
Ejemplo 29
Obtener la integral de circulación del campo vectorial definido por: F = xy i + y 2 j en la curva
C definida por la frontera de la región acotada por las gráficas y = x 2 y y = x en sentido contrario
a las agujas del reloj.
Solución:
La figura correspondiente es:

Figura 9: Superficie de aplicación de Teorema de Green

Aplicamos la ecuación correspondiente al teorema de Green

N M
 C
F .dr =   (
S
−
x y
)dA

M ( x, y ) = xy
F = f (M , N )
N ( x, y ) = y 2

N M y 2 xy
 F = ( − ) k = ( − ) k = − x k
x y x y
Dado que el vector apunta en la dirección k, en las operaciones posteriores omitiremos la escritura
del vector unitario en dicha dirección.

La intersección entre x e y se obtiene igualando las ecuaciones: y = x y y = x , es decir; x = 0

2

y x = 1 , por lo tanto:
1 x 1 1

 F .dr =   − xdA = .  − xdydx =  − x  y x2 dx =  − x( x − x 2 )dx

x
C S
0 x2 0 0

1
1
 x3 x 4  1 1 1
=  (− x + x )dx = − +  = − + = −
2 3

0  3 4 0 3 4 12

El significado del de los signos en el resultado es el siguiente:

Positivo: Nos movemos por la trayectoria en el campo mayormente a favor; es como andar con el
viento a favor.
Negativo: Nos movemos por la trayectoria en el campo mayormente en contra; es como andar con
el viento en contra.

Ejemplo 30

Obtener la integral de circulación del campo vectorial F definida por la siguiente expresión:
F = ( x + e x seny ) i + ( x + e x cos y ) j en la curva C definida por el lazo derecho de la lemniscata
 2 = cos 2 en el sentido contrario a las agujas del reloj.
Solución
La lemniscata se muestra en la figura de la izquierda y en la figura de la derecha se muestra el la
hoja derecha de la lemniscata-

 = cos 2
1

2

Figura 10: Superficie de aplicación del Teorema de Green

Se trata de obtener la integral del rotacional en el interior de la curva (color verde). De esta
manera, el rotacional esta dado por:

N M ( x + e x cos y) ( x + e x seny)
 F = − = − = (1 + e x cos y) − e x seny = 1
x y x y
El resultado es un vector en la dirección k; sin embargo, dado que la función esta en el plano el
vector direccional apuntara únicamente en la dirección k por lo que se omite esta dirección.
Previamente, hallamos los puntos de intersección para obtener los límites de integración

 
 2 = cos 2 si  = 0 se tiene cos 2 = 0 2 =  luego  = 
2 4

En consecuencia, se debe calcular la siguiente igualdad

 /4 cos 2  /4 cos 2  /4
2
 C
F .dr =   1dA =
S  
− /4 0
 d  d = 
− /4
2
d = 
− /4
cos 2 d
0

1  /4 1 1
= sen2 − /4 = (1 + 1) =
4 4 2
Nota: En la integral se ha agregado  debido a la transformación de coordenadas

Este teorema puede ser aplicado para calcular el flujo de fluidos a través de campos vectoriales.
Teorema de Stokes
Sea S una superficie orientada, simple y suave por tramos. Sea C su curva frontera, regular a
trozos, cerrada y simple, con orientación positiva. Si F es un campo vectorial cuyas componentes
tienen derivadas parciales continuas en una región R 3 que contiene a S. Entonces:

C
F .dr =  rotF .dS
S
(7.61)

Para determinar la orientación positiva de la curva C frontera de S, convenimos en que, al recorrer

C en sentido positivo con la cabeza apuntando al vector normal n que indica la orientación
positiva de S, la superficie queda a la izquierda.
El teorema de Stokes proporciona otra extensión del teorema fundamental de la integral al
relacionar una integral de superficie con la integral de línea sobre la curva frontera a dicha
superficie. Es decir, es una generalización del Teorema de Green

Ejemplo 31
Verificar que se que se cumple el Teorema de Stokes en el campo de flujo dado por la siguiente
expresión: F ( x, y, z) = − y i + x j − 2 k y la superficie S en el cono z = x + y , 0  z  4
2 2 2

con orientación hacia abajo.

Solución
Trataremos de probar el término de la izquierda del Teorema de Stokes

Para probar el primer término del Teorema de Stokes consideremos a F un campo vectorial
continuo definida sobre una curva suave C dada por una función vectorial r (t ) , cuyo intervalo
de variación de su parámetro es: a  t  b , entonces la integral de línea de F a lo largo de C es
b

 C
F .dr =  F (r (t ))r (t )
a

r (t ) y r (t ) son funciones vectoriales

Figura 11: Esquema para la aplicación del Teorema de Stokes

Cuando se hace la integral de línea lo que interesa es la curva que se genera en z = 4 (borde de
color negro) y no importa la superficie. Como esta orientado hacia abajo (en contra de las
manecillas del reloj, el flujo mayoritariamente sale del cono)

Figura 12: Superficie de circulación para la aplicación del Teorema de Stokes

Si consideramos en contra de las agujas del reloj luego de transformar a coordenadas polares se
tiene:

x = 4cos t y = 4sent z=4

Si consideramos en dirección de las agujas del reloj luego de transformar a coordenadas polares
se tiene:

x = 4cos t y = −4sent z = 4 con 0  t  2

r(t ) = −4senti − 4cos t j + 0 k

F (r (t )) = −(−4sent )i + (4cos t ) j − 2 k

2

 C
F .dr =  ((4sent )
0
i + (4cos t ) j − 2 k ).(−4sent i − 4cos t j + 0 k )dt

2 2 2

 F .dr =  (−16sen t − 16cos t )dt = −16  (sen t + cos t )dt = −16  tdt = −32
2 2 2 2
C
0 0 0

Esto significa que el flujo es a favor de las manecillas del reloj; significa que mayoritariamente el
flujo sale del cono.
Ejemplo 32
Verificar que se que se cumple el Teorema de Stokes en el campo de flujo dado por la siguiente
expresión: F ( x, y, z) = − y i + x j − 2 k y la superficie S en el cono z = x + y , 0  z  4
2 2 2

con orientación hacia abajo.

Solución
Trataremos de probar el término de la derecha del Teorema de Stokes, es decir utilizando la
integral de superficie.
La notación es la siguiente:
 rotF . dS =  (rotF ). ndS
S S

El cono del ejemplo 31 visto desde la parte superior se puede graficar del siguiente modo:

Figura 13: Región R para la aplicación del Teorema de Stokes

El rotacional de F ( x, y, z) = M  i + N j + P k es: rotF ( x, y, z ) =   F ( x, y, z)

 P N   P M   N M 
 F =  −  i +  −  j +  − k
 y z   x z   x y 
Si la superficie fuese plana, entonces tendríamos solo el componente de color rosado; en estas
condiciones representa el Teorema de Green; es decir este último teorema es un caso particular
del Teorema de Stokes.
De acuerdo con los datos suministrados se tiene:

M = −y N=x P = −2

Luego, aplicando las derivadas parciales

P N P M N M
=0 =0 =0 =0 =1 = −1
y z x z x y
Luego:

rotF ( x, y, z) = (0)i + (0) j − (−1) k = 2 k

Para calcular el vector normal n , sea S una superficie orientada por z = g ( x, y ) y sea R su
proyección sobre el plano entonces:

 F . nds = F. − g ( x, y)

S S
x i − g y ( x, y) j +  k  ds cuando está orientada hacia arriba

 F. nds = F.  g ( x, y)

S S
x i + g y ( x, y) j −  k  ds cuando está orientada hacia abajo

En esta caso, está orientada hacia abajo

g ( x, y ) = z = x 2 + y 2

Derivando parcialmente:
 2x 2y
g x ( x, y ) = g y ( x, y ) =
2 x2 + y 2 2 x2 + y 2

x i y j
 rotF . nds =  ((0)
S R
i + (0) j +(2 k )).(
x2 + y 2
+
x2 + y 2
−  k )ds

2 4

 rotF . nds =  (0 + 0 − 2)ds = −2   rdrd = −2 (4) = −32

2

S R 0 0

Dado que los dos métodos desarrollados en los ejemplos 30 y 31 nos han reportado el mismo
valor, queda demostrado el Teorema de Stokes.

Teorema de Gauss
Este teorema define una forma de calcular una integral de un campo vectorial sobre una
superficie, a través de una integral de volumen.
Supongamos que:

F ( x, y, z ) es un campo vectorial tridimensional

v es un volumen del elemento de fluido
S es la superficie de v
n es una función cuya salida son vectores unitarios normales en la superficie
Entonces el teorema de la divergencia de Gauss afirma lo siguiente:

 divFdv = 
v s
F .nds (7.62)

La idea intuitiva aquí es que ambas integrales miden la tasa a la que un fluido que fluye a lo largo
del campo vectorial F sale de la región v (o entra en v si los valores de ambas integrales son
negativos). La integración triple de la divergencia implica contar todos los pedacitos de flujo hacia
afuera de v mientras que evaluar la integral de flujo mide esto al verificar cuánto fluido sale o
entra a lo largo de la frontera de v

Expresada en forma matemática: dado un volumen v de fluido en una región cerrada, acotada
por una superficie A , entonces se tiene:

 (.V )dv = 

v S
(nV
. )dA (7.63)

En la práctica, es mucho mas complicado utilizar la integral cerrada y en su lugar se calcula el

flujo volumétrico utilizando la integral volumétrica, es decir:

 S
. )dA =  (.V )dv
(nV
v
(7.64)

Ejemplo 29

La velocidad de un fluido que atraviesa un cubo unitario limitado por 0  x  1 , 0  y  1

0  z  1 está dada por: V = 10 xy1 − 5xz 23 . Encuentre el flujo volumétrico utilizando la
integral cerrada.
Solución:
El caudal o flujo volumétrico se calcula mediante:

Q=  S
. )ds =
(nV  S
V .d A

Aquí va la figura

dAxz

Figura 14: Elemento de volumen para la aplicación del Teorema de Gauss

Superficie (1): Plano xz a y = 0 dA = − 2 dxdz (Es negativo por que el fluido entra al volumen
de control).
x =1 z =1
Qx ,0, z =   (10 xy
x =0 z =0
1 − 5 xz 2 3 ).(−1dxdz ) = 0

Superficie (2): Plano xz a y = 1 dA =  2 dxdz (Es positivo porque el fluido sale del volumen de
control).
x =1 z =1
Qx ,0, z =   (10 xy
x =0 z =0
1 − 5 xz 2 3 ).(1dxdz ) = 0

Superficie (3): Plano xy a z = 0 dA = − 3 dxdz (Es negativo porque el fluido entra al volumen
de control).
x =1 y =1
Qx , y ,0 =   (10 xy
x =0 y =0
1 − 5 xz 2 3 ).(− 3dxdz ) = 0

Superficie (4): Plano xy a z = 1 dA =  3dxdz (Es positivo porque el fluido sale volumen de
control).
x =1 y =1 x =1 y =1

    ( − 5x)dx = −2.5x
2 1
Qx , y ,0 = (10 xy1 − 5 xz 2 3 ).( 3dxdz ) =  = −2.5 m3/s
0
x =0 y =0 x =0 y =0

Superficie (5): Plano yz a x = 0 dA = −1dxdz (Es negativo porque el fluido entra al volumen
de control).
y =1 z =1

Q0, y , z =   (10 xy

y =0 z =0
1 − 5 xz 2 3 ).(−1dxdz ) = 0
Superficie (6): Plano yz a x = 1 dA = 1dxdz (Es positivo porque el fluido sale volumen de
control).
y =1 z =1 y =1 z =1 1
1
10 y 2 
Q1, y , z =   (10 xy1 − 5 xz  3 ).(1dxdz ) =   (10 y)dydz =  10 ydy = 
2
 = 5 m /s
3

y =0 z =0 y =0 z =0 0  2 0
6
Q =  Qi = −2.5 + 5 = 2.5 m3/s
1

Ejemplo 30

La velocidad de un fluido que atraviesa un cubo unitario limitado por 0  x  1 , 0  y  1

0  z  1 está dada por: V = 10 xy1 − 5xz 23 . Encuentre el flujo volumétrico utilizando la
integral cerrada.
Solución:
El caudal o flujo volumétrico se calcula mediante la integral volumétrica como consecuencia de
la aplicación del Teorema de Gauss, es decir:

Q =  (.V )dv

v

Ahora podemos calcular la divergencia de la velocidad ( .V )

   (10 xy) (0) (−5 xz 2 )

.V = (1 +  2 +  3 ).(10 xy1 − 5 xz  3 ) =
2
+ +
x y z x y z

.V = 10 y + 0 − 10 xz
Luego:
1 1 1 1 1
Q =  (.V )dv =    (10 y − 10 xz )dxdydz =   (10 y − 10 xz )dxdy
v
0 0 0 0 0

1
 5 
1
Q =  (5 − 5 x)dx = 5 x − x 2  = 5 − 2.5 = 2.5 m3/s
0  2 0

Hay dos teoremas estrechamente relacionados derivados del Teorema de Gauss-Ostrogradsky y

es el que se aplica los escalares y tensores.

 (s)dv =  nsds

v s
(7.65)

 (. )dv = 

v s
 n.  ds (7.66)

Nuevamente estas ecuaciones tienen gran aplicabilidad al resolver la integral volumétrica, en

lugar de la integral de superficie, debido a su simplicidad.
3
El teorema indica lo siguiente: si F es un campo vectorial continuo y diferenciable en R , el
flujo del campo F a través de una superficie S en la dirección del vector unitario normal exterior
a la superficie, es igual a la integral de la divergencia del campo sobre la región v encerrada en la
superficie.

 F . nds = .Fdv

S v
(7.67)

El primer término representa el flujo del campo a través de la superficie S; el segundo término es
la integral de la divergencia del campo en el volumen v encerrado por la superficie. Hay que
recordar que la divergencia indica cuanto sale un campo de un punto es decir una fuente o cuanto
se dirige hacia el siendo un sumidero.
Si el flujo es positivo entonces el campo sale de la superficie
Si el flujo es negativo entonces el campo entra en la superficie.

ds

v
ds

Figura 15: Campo de flujo y dirección de desplazamiento del fluido

En base a la figura 16 el volumen total se divide en infinitos diferenciales de volumen que son tan
pequeños que se reducen a puntos; en alguno de ellos el campo sale, mientras que, en otros, el
campo es entrante. Si sumamos estas infinitas entradas y salidas de campo de diferenciales de
volumen encerrado por S, tendremos el flujo total del campo vectorial.
 ds
ds

Figura 16: Representación de los infinitos diferenciales de volumen dentro del campo de flujo
vectorial

Ejemplo 31

Obtener el flujo del campo F = y i + xy j − z k a través de la superficie cerrada formada por
el cilindro x 2 + y 2  4 , el plano z = 0 y el paraboloide z = x 2 + y 2 .

Solución
El esquema es el siguiente:

Luego aplicamos el Teorema de Gauss; calcular la integral del primer miembro es bastante
compleja pues tendriamos que calcular la integral del cilindro, del paraboloide y para el plano
base. Lo mas sencllo es obtener la segunda integral de la ecuacion (a)

 F . nds = .Fdv (a)

S v

Luego la divergencia de F esta dada por:

Fx Fy Fz y xy (− z )

.F = + + = + + = x −1
x y z x y z
 4− x2 x + y
2 2
2

 .Fdv =  ( x − 1)dxdydz =  

v v −2 − 4 − x 2

0
( x − 1)dzdydx

Debido a la dificultad de integrar la ecuación lo transformaremos a coordenadas cilíndricas, para

lo cual hacemos los siguientes cambios:

x =  cos  0    2
y =  cos  0    2
z=z 0  z  2

Hay que recordar que

z = x2 + y 2 z = (  cos  )2 + (  cos  )2 =  2 por otor lado del resultado del integrando se

tiene: x − 1 =  cos  − 1

Luego, la integral se transforma en:

4− x2 x + y 2 2 
2 2 2
2

 .Fdv =  
v −2 − 4 − x 2

0
( x − 1)dzdydx =    (  cos − 1)  dzd d 
0 0 0

Se introduce  para hacer el cambio de coordenadas rectangulares a cilíndricas

2 2 2 2 2 2
2
=   (  cos  − 1)  z 0 d d  =   (  cos  −  )d d  =  (  cos  −   )
4 3 4 3

0 0 0 0 0 0

2
2
 2 4 
=  −2  d  =  −
3
= −8
0  4  0

Dado el signo negativo, el campo de flujo entra a través de la superficie

7.6. VALORES PROMEDIO DE FUNCIONES

Valor medio respecto al tiempo
Dada una función del tipo F ( x, y, z, t ) , el valor medio denotado por:  F  en un intervalo de
tiempo, viene dado por:
t1 t1

Fdt   Fdt
Fdt
 F =
 =
t t0
=
t0
(7.68)
 dt  dt (t1 − t0 )
t1
t

t0

Ejemplo 32
La velocidad de una partícula en un sistema de flujo turbulento durante el intervalo de tiempo
0  t   / 2 está dado por:
V = 4 − 1.2sent + 0.8cos t
Determine la velocidad media de la partícula durante ese intervalo de tiempo
Solución
Aplicamos la ecuación (1)

t  /2

 Vdt  (4 − 1.2sent + 0.8cos t )dt (4t − 1.2cos t − 0.8sent ) 0

 /2

 V = 0
= 0
 /2
=
t
 /2
 dt0
 dt
0

2 − 1.2(0 − 1) − 0.8(1 − 0)
= = 4 + 0.8 /  = 4.2546 m/s
 /2
Valor promedio de una función F a lo largo de una trayectoria
El valor promedio de una función sobre una trayectoria s viene en el intervalo 0  s  L dada
por:
L

 Fds = 
Fds
 F = s 0
(7.69)
 ds  ds
L
s
0

6.3 Valor promedio de una función F sobre una superficie A

  FdA   FdA
 F = A
= A
(7.70)
  dA A
A

Ejemplo 33
La velocidad media de un fluido en un cilindro horizontal de radio r viene dado por:

  r 2 
Vz = 100 1 −   
  R  

Encuentre la velocidad media

Solución
Calculamos la velocidad media usando la ecuación anterior, para lo cual se deberá integrar en los
limites correspondientes del tubo.
2 R
  r 2  25r 4
R

  V (r )dA
z 0 0 1 −  R   rdrd
100 (50r − 2 )
2

2 (50 − 25) R 2
 V = A
=   = 2
R 0
= = 50 m/s
2 R
  dA  R2  R2
A   rdrd
0 0
Valor promedio de una función F dentro de un volumen v
El valor medio de una función dentro de contenido en un volumen se calcula mediante la siguiente
relación:

 F =
 Fdv =  Fdv
v v
(7.71)
 dv
v
v

Ejemplo 34

El perfil de temperatura de un sólido esférico de radio R esta dado por T = 50r . Encuentre la
2

temperatura media del solido esférico.

Solución
La temperatura media se obtiene haciendo uso de la ecuación anterior y con la diferencial de
volumen obtenido en el ejemplo 23.
2  R R
50r 5  2

T (r )dv    
(50r )r sen drd d
2 2
− cos  0  0
 T =
 v
= =0 =0 r =0
=
5 0
 
 dv
2 R R
r3
   r sen drd d
 2
cos  0  0
2
v −
  =0 =0 r =0 3 0

−10 R5 (cos  − cos 0)  2

= 3
= 30 R 2
R
− (cos  − cos 0)  2
3
7.7. TEOREMA DE TRANSPORTE DE REYNOLDS
El Teorema de Reynolds es muy importante pues nos permite relacionar los enfoques
lagrangianos y eulerianos del flujo de fluidos. El enfoque lagrangiano sigue el movimiento de las
partícula; el enfoque euleriano permite fijar un volumen de control referido al fluido; como
resultado podemos encontrar la conservación de masa total (ecuación de continuidad),
conservación de cantidad de movimiento, conservación de energía y conservación de momento
lineal.
En esta primera parte obtendremos la ecuación de conservacion macroscópica de materia
expresada en términos de ecuaciones diferenciales y luego bajo ciertas condiciones obtener la
ecuación de continuidad en su forma macroscópica conocida en régimen estacionario.

Para el enfoque lagrangiano la figura delimitada con color rosado se denomina sistema y lo que
está contenido en esa área es una propiedad extensiva B (masa, cantidad de movimiento, energía
o momento angular); el enfoque euleriano toma un volumen de control fijo desde el punto de vista
del observador.
Consideremos que esto ocurre en un instante t; en decir:
Figura 18: Representación del sistema y volumen de control
En el instante t se tiene:
t
Bsist = BVct (7.72)

Ahora cuando el tiempo es t + t

Figura 19: Desplazamiento del sistema debido al flujo

El sistema se desplaza (color rosado), mientras el volumen de control permanece fijo (color azul)
Ahora estamos interesados en estudiar las propiedades del fluido contenido en el sistema marcado
con color rosado, como consecuencia de ello aparecen dos regiones (I) y (II)

Figura 20: Esquema que relaciona el sistema y el volumen de control

Por tanto, en base a la figura se tiene:
t +t t +t
Bsist = BVc − BIt +t + BIIt +t (7.73)

Restando la ecuación (7.2) de la ecuación (7.1)

 t +t
Bsist − Bsist
t
BVct +t − BVct BIt +t BIIt +t
lim = lim − + (7.74)
x →0 t x →0 t t t
Aplicando la derivada correspondiente se tiene

dBsist dBVc BIt +t BIIt +t dBVc • •

= − + = − BI + BII (7.75)
dt dt t t dt
• •
BI y BII representan flujos.

El termino de flujo se puede transformar considerando m = B la propiedad extensiva del

siguiente modo:

B m  v Qt
= = = = VA (7.76)
t t t t
De esta manera, la ecuación (3) escrita en términos de la masa, permite obtener

dmsist dmVc • •
= − mI + mII (7.77)
dt dt
De acuerdo con la ecuación (II), la masa del sistema no ha reaccionado, ni se ha agregado ni
quitado nada, por lo tanto, no hay cambios; en consecuencia, desde el punto de vista lagrangiano,
dmsist
se tiene: = 0 ; es decir, la masa del sistema se conserva.
dt
Desde el punto de vita euleriano se tiene:

dmVc • •
= − mI + mII (7.78)
dt
Esta ecuación se puede hacer extensiva para varias entradas y salidas considerando el termino (I)
correspondiente a los flujos de entrada y el termino (II) flujos de salida.

dmVc • •
= − me +  ms (7.79)
dt
dmVc
Considerando que se opera en régimen estacionario, entonces = 0 , y además, un solo flujo
dt
de entrada y un solo flujo de salida; en consecuencia, la ecuación (7.8) se transforma en:
1  V1  A1 = 2  V2  A2 ; para el caso de un fluido incompresible 1 = 2 ; por tanto
finalmente se tiene:  V1  A1 = V2  A2 y de esta ecuación se deduce que la relación de
velocidades medias en un tubo son inversamente proporcionales al cuadrado de sus cocientes:
2
 V1   d 2 
=  (7.80)
 V2   d1 

En esta segunda parte encontraremos una ecuación mas generalizada del Teorema de Reynolds,
las mismas que pueden ser aplicadas a los tres fenómenos de transporte, referido a una propiedad
extensiva.
Desde la visión euleriana el movimiento del fluido es analizado desde el punto de vista de un
sistema inercial de referencia; es decir observamos el movimiento de fluido desde un punto fijo
en movimiento rectilíneo uniforme tal como se observa en la figura.

Figura 21.- Puntos de vista euleriano y lagrangiano del flujo de una porción de fluido
Desde la visión lagrangiana no se considera un sistema inercial de referencia; es decir se considera
un sistema no inercial de referencia; es decir el observador se desplaza a la misma aceleración
que la porción de fluido. Un ejemplo típico considerado es aquel sistema donde se encuentra una
canoa con el observador en el cause de un rio y a partir del movimiento de la canoa se analiza el
movimiento del fluido. Esta descripción lagrangiana se aplica principalmente a la mecánica de
los sólidos.
También es importante recordar el concepto de caudal o flujo volumétrico

Q=  V .dA =  V . ndA =  V  A
sc sc
(7.81)

Cabe recordar que el vector velocidad y el vector unitario normal al área de entrada o salida del
fluido, deben ser paralelos. En el caso de la entrada son paralelos pero opuestos; por tanto, el
signo correspondiente al caudal será negativo. A la salida son paralelos y tienen la misma
dirección, por tanto, el signo del caudal será positivo.
Del mismo modo se puede obtener el flujo de masa.
•
m =   dQ =  V . ndA =  Q (7.82)
sc

El Teorema de Reynolds relaciona el sistema y el volumen de control.

Volumen de control: Porción de fluido, limitado por los extremos, inmóvil
Sistema cerrado: Porción de fluido que se mueve con el fluido como consecuencia de ello no hay
un cambio neto.
Consideremos a B como una propiedad extensiva, la cual puede se la masa ( m ), la cantidad de
movimiento ( p = mV ), la energía total ( E = U + Ec + E p ), momento angular (
L = r  p = r  mV )
 B
Sea  una propiedad extensiva especifica (propiedad por unidad de masa), es decir,  =
m
entonces aplicando este criterio para cada una de las propiedades extensivas se tiene:

B m
Para la masa:  = = =1
m m
B mV
Para la cantidad de movimiento:  = = =V
m m

B (U + Ec + E p ) E
Para la energía:  = = = =e
m m m
B L r  p r  mV
Para el momento angular:  = = = = = r V
m m m m
Ahora consideremos el siguiente esquema

Figura 22.- Volumen de control, superficie de control y sistema

En el esquema, el volumen de control es fijo o estático (color azul); (aunque hay casos donde el
volumen de control se mueve) el borde se denomina superficie de control (azul) y el sistema es
móvil, se desplaza del punto 1 al punto 2 (color rosado achurado); sin embargo, se mantiene
constante.
Se trata de relacionar el sistema con el volumen de control con el sistema. El Teorema de
Reynolds, está dado por:

 dB  d
 dt 
sist
=
dt Vc
 dv +   (V . n)dA
Sc
(7.83)

Si aplicamos esta ecuación a la conservacion de materia  = 1 , luego la ecuación (7.12) se

transforma en:

 dm  d
dt Vc
 dt  =  dv +   (V . n)dA (7.84)
Sc
sist

Se observa que, si bien el sistema se mueve, ingresa la misma masa que la masa que sale, por
tanto (dm dt ) sist = 0 . Se hace esta consideración por dos causas únicamente: Cuando la masa es
constante o cuando la masa no es función del tiempo.
Luego la ecuación se reduce a:
d
dt Vc
 dv +   (V . n)dA = 0 (7.85)
Sc

Esta ecuación también se denomina ecuación de continuidad

d d
En estado estacionario se tiene: 
dt Vc
 dv =  dv =0 , puesto que el volumen de control
dt Vc
es fijo. Con esta condición se obtiene:

 Sc
 (V . n)dA = 0 (7.86)

Esta ecuación se reduce a: 2V2 A2 − 1V1 A1 = 0 , que es la ecuación de continuidad macroscópica

para referido a la conservacion de masa total en régimen estacionario.
Desde el punto de vista del enfoque euleriano consideramos para nuestro estudio el volumen de
control debido a que las moléculas se mueven desordenadamente, no se puede utilizar un sistema
de masa fija; por tanto, debemos deducir las ecuaciones de dinámica de fluidos a partir de los
cambios de masa, momentun, energía, etc. A medida que el fluido pasa por el volumen de control.
La frontera del volumen se denomina superficie de control.

Figura 23: La evolución del sistema y el volumen de control

Quizá, entender este fenómeno es relativamente complejo; sin embargo, su aplicación es mucho
más sencillo ya sea en régimen estacionario o en régimen no estacionario.

Vemos que el sistema se mueve con el fluido mientras que el volumen de control permanece en
su posición original. Sea B una propiedad extensiva cualquiera (masa, energía, momentun, etc.)
dentro del sistema en un tiempo t ; y  la misma propiedad por unidad de masa. La tasa temporal
de incremento de B para el sistema se formula en base a términos de volumen de control. En el
tiempo t + t el sistema ocupa el volumen II y III, mientras que en t ocupaba solo II. El
incremento de la propiedad B del sistema en el tiempo t está dado por:

t +t
Bsist − Bsist
t (   dv +   dv)t +t − (   dv)t
lim = II III II
(7.87)
t →0 t t

 dB  (   dv)t +t − (   dv)t   dv)t +t

 dt  = II II
+ III (7.87)
sist  t t
En la ecuación (7.87) por aplicación del teorema de Gauss, el ultimo término de la derecha, una
integral volumétrica, se convierte en una integral de superficie.

 dB  d
 dt  = dt Vc  dv +  Sc  (V . n)dA (7.88)
sist

Si consideramos que B sea la masa y  = B / m , entonces la ecuación (7.89) se transforma en:

 dm 

d
 dt  = dt Vc  dv + Sc
 (V . n)dA (7.89)
sist

La ecuación (7.89) es el Teorema de Transporte de Reynolds aplicado al transporte total de masa.


d
dt Vc
 dv +  (V . n)dA = 0 (7.90)
Sc

La ecuación (7.90) representa la ecuación de conservacion de masa total.

Podemos también relacionar el enfoque de sistemas y el enfoque de volúmenes de control del
siguiente modo:
Al fijar la atención sobre una región finita del espacio, como en el punto de vista euleriano, estaría
utilizándose el método del volumen de control. Ahora, podremos relacionar el enfoque del sistema
y el enfoque del volumen de control para ciertos fluidos y ciertas propiedades del flujo que se
describirán a continuación.

Consideremos un campo de flujo arbitrario V ( x, y, z, t ) visto desde una referencia xyz , en el

cual se observa un sistema fluido de masa finita en los tiempos t y t + t como se muestra en una
forma ideal en la figura 24 mediante la línea continua y la línea punteada, respectivamente. Las
líneas de corriente corresponden al tiempo t . Además de este sistema, se considerará que el
volumen del espacio ocupado por el sistema en el tiempo t es un volumen de control fijo en xyz
. Luego, en el tiempo t el sistema es idéntico al fluido contenido en el volumen de control que se
muestra con la línea continua. Considérese alguna propiedad extensiva arbitraria B del fluido
con el propósito de relacionar la tasa de cambio de esta propiedad para el sistema, con las
variaciones de esa
misma propiedad asociadas con el volumen de control. La distribución de B por unidad de masa
estará dada como  de manera que B =    dv donde dv representa un elemento de volumen
 Figura 24 Vista simplificada de un sistema en movimiento

 dB 
=
DB
=
III  dv + II  dv)t +t − (I  dv + II  dv)t )
  lim( (7.91)
 dt  sist Dt t →0 t

(   dv −   dv )   dv   dv

DB II II III I
= lim t +t t
+ lim t +t
− lim t
(7.92)
Dt t →0 t t →0 t t →0 t
Luego estas integrales deben ser tratadas por partes:

(   dv −   dv )
d
dt 
 dv
II II
lim t +t t
= (7.93)
t →0 t vc

El segundo miembro de la ecuación (7.93) representa la tasa de cambio de B dentro del volumen
de control que tiene la forma del sistema en el tiempo t

El segundo límite de la ecuación (7.92) puede considerarse que la integral aproxima la cantidad
de la propiedad B que atraviesa parte de la superficie de control mostrada esquemáticamente
como ARB en la figura 24 durante el tiempo t , de manera que  dv / t aproxima la tasa 
III t +t

promedio de salida de B a través de ARB durante el intervalo t . En el límite cuando t → 0

, esta relación se convierte en la tasa exacta de flujo de salida de B a través de la superficie de
control. De modo similar el ultimo limite de (7.92) puede considerarse que para flujos con
propiedades y características de flujo continuo, la integral   dv
I t
se aproxima a la cantidad

de B que ha ingresado al volumen de control en el intervalo t a través de la porción restante de

la superficie de control, mostrada en la figura 24 como ALB. En el límite   dv / t se
I
convierte en la tasa exacta de flujo de entrada de B hacia el volumen de control en el tiempo t
indicado. Luego las dos últimas integrales de la ecuación (7.93) dan la tasa neta de salida de B a
través del volumen de control en el tiempo t como:

(   dv ) (   dv )
III I
lim t +t
− lim t
(7.94)
t →0 t t →0 t
La ecuación (7.94) es la tasa neta de flujo de salida de  a través de la superficie de control y
puede ser expresada de una forma más compacta y útil. En la figura 25 se tiene un campo de
velocidad de flujo permanente y una porción de superficie de control, en el que se muestra un dA
es esta superficie. Esta área también es la interfaz de fluido que esta apenas en contacto con la
superficie de control en el instante t , tal como se muestra en la figura 25 a). En la figura 25 b)
también se muestra la interfaz de fluido en el tiempo t + t . Se observa que la interfaz se ha
movido una distancia Vdt a lo largo de la distancia tangente a la línea de corriente en el punto. El
volumen de fluido que ocupa la región cubierta por dA en el tiempo dt , forma un “tubo” de
corriente cuyo volumen es:

dv = (Vdt )(dA cos  ) = V . dAdt (7.95)

 (a) (b)

Figura 25. a) Interfaz dA en la superficie de control en el tiempo t , b) Interfaz dA en la superficie

de control en el tiempo t + t

El dv es el volumen defluido que ha cruzado dA en la superficie de control en el tiempo dt .

Multiplicando por  y dividiendo por dt se obtiene la tasa instantánea de flujo de masa de fluido,
V . dA que sale del volumen de control a través del área dA indicada.

Figura 26: Superficie de control que muestra el flujo de entrada de masa

La tasa neta de flujo de salida de B a través de la superficie de control elegida puede darse
aproximadamente mediante:

Tasa de flujo de salida de la SC =   ( V . dA) , que es el flujo de

ARB
B por unidad de tiempo a

tardes de dA .

En base a la figura 26 se observa que para el fluido que entra al volumen de control, expresión
V . dA es negativa debido al producto punto. Luego la expresión de tasa de flujo de entrada B a
través de la superficie de control requiere un signo negativo para que el valor de dicho resultado
sea positivo, es decir:

Tasa de flujo de entrada a través de la SC = −   ( V . dA) , que es el flujo de B por unidad de

ALB

tiempo a tardes de dA .

Luego la tasa neta de flujo de salida de B esta dada aproximadamente por la relación:
Tasa neta de flujo de salida= tasa de flujo de salida en ARB  Tasa de flujo de entrada en ALB

  ( V . dA) = −(−   ( V . dA))

ARB ALB
(7.96)

En el límite cuando, las aproximaciones se vuelven exactas, de manera que el lado derecho de la
ecuación (7.96) puede ser expresada como:   ( V . dA) donde la integral se hace sobre una
SC

superficie cerrada en toda la superficie de control. Luego la ecuación (7.94) puede expresarse
mediante:

Tasa neta de flujo de salida de B desde VC =   ( V . dA)

SC
(7.97)

Ahora utilizando las ecuaciones (7.93) y (7.97) se obtiene:

DB 
Dt
=   ( V . dA) + t   dv
SC vc
(7.98)

La ecuación (7.98) se conoce como el Teorema de Reynolds. Para el caso de transporte de masa
se tiene  = 1 .

7.8 TABLAS DE LAS ECUACIONES DE TRANSPORTE DE CANTIDAD DE

MOVIMIENTO EN LOS DIFERENTES SISTEMAS DE COORDENADAS.
FORMA MICROSCOPICA A TRAVES DE ECUACIONES DIFERENCIALES

7.8.1 ECUACION DE CONTINUIDAD

Una de las formas de la ecuación de continuidad, las cuales pueden ser escritas en los tres sistemas
de coordenadas.


= − . V (7.99)
t
Ecuación de continuidad en coordenadas rectangulares

 ( Vx ) ( Vy ) ( Vz )

+ + + =0 (7.100)
t x y z
Ecuación de continuidad en coordenadas cilíndricas

 1 (  rVr ) 1 ( V ) ( Vz )

+ + + =0 (7.101)
t r r r  z
Ecuación de continuidad en coordenadas esféricas

 1 (  r 2Vr ) 1 ( V sen ) 1 ( V )

+ 2 + + =0 (7.102)
t r r rsen  rsen 
7.8.2 ECUACION DE CONSERVACION DE CANTIDAD DE MOVIMIENTO
La ecuación de movimiento que es equivalente a la segunda ley de Newton aplicada a la dinámica
de fluidos puede ser aplicado tanto a fluidos newtonianos como a fluidos no newtonianos; de
forma general puede ser escrita del siguiente modo:

DV
 = −p −   +  g (7.103)
Dt
Coordenadas rectangulares: Direcciones x , y y z respectivamente.

Vx V V V p   
( + Vx x + Vy x + Vz x ) = − − ( xx + yx + zx ) +  g x (7.104)
t x y z x x y z

Vy Vy Vy Vy p  xy  yy  zy

( + Vx + Vy + Vz )=− −( + + ) + gy (7.105)
t x y z y x y z

Vz V V V p   
( + Vx z + Vy z + Vz z ) = − − ( xz + yz + zz ) +  g z (7.106)
t x y z z x y z
Coordenadas cilíndricas: Direcciones r ,  y z respectivamente.

Vr V V V V 2 V p 1 (r rr ) 1  r    rz
( + Vr r +  r −  + Vz r ) = − − ( + − + ) +  gr
t r r  r z r r r r  r z
(7.107)

V V V V V V V 1 p 1 (r 2 r ) 1     z
( + Vr  +   + r  + Vz  ) = − −( 2 + + ) +  g
t r r  r z r  r r r  z
(7.108)
 Vz V V Vz V p 1 (r rz ) 1  z  zz
( + Vr z +  + Vz z ) = − − ( + + ) +  gz (7.109)
t r r  z z r r r  z

Coordenadas esféricas: Direcciones r ,  y  , respectivamente

Vr V V Vr V Vr V2 + V2 p 1  (r 2 rr ) 1  ( r sen )

( + Vr r +  +  − ) = − −( 2 +
t r r  rsen  r r r r rsen 

1  r   +  
+ − ) +  gr (7.110)
rsen  r

V V V V V V VrV V2 cot  1 p 1  (r 2 r )

( + Vr + + + − )=− −( 2
t r r  rsen  r r r  r r

1 (  sen ) 1    r cot 

+ + + −   ) +  g (7.111)
rsen  rsen  r r

V V
V V V V VVr VV 1 p 1  (r  r )
2

( + Vr + + + + cot  ) = − −(
t r r  rsen  r r rsen  r 2 r

1   1    r 2cot 
+ + + +   ) +  g (7.112)
r  rsen  r r
7.8.3 ECUACIÓN DE NAVIER-STOKES.
Es una ecuación modificada para ser aplicada al caso de flujo de fluidos newtonianos de
viscosidad y densidad constantes; respecto a la ecuación anterior, se ha modificado el termino
correspondiente a las fuerzas viscosas, empezando por el esfuerzo de corte generalizado dada por
la ecuación siguiente

V j Vi 2 V V V
 i j = − ( + ) + (  − k )( x + y + z ) ij (7.113)
xi x j 3 x y z

La ecuación (7.103) para cada uno de los sistemas de coordenadas, omitiendo el termino que
contiene el producto de k por el tensor V está dada por:

Coordenadas rectangulares: Direcciones x , y y z respectivamente.

Vx 2 Vy 2
 xx = −2 +  ( V )  yy = −2 +  ( V )
x 3 y 3
Vz 2  Vx Vy 
 zz = −2 +  ( V )  xy =  yx = −   + 
z 3  y x 

 Vy Vz   Vz Vx 

 yz =  zy = −   +   zx =  xz = −   + 
 z y   x z 
 Vx Vy Vz
V = + +
x y z

Coordenadas cilíndricas: Direcciones r ,  y z respectivamente.

 Vr 2 (   V 1 Vz 

 rr = −   2 −   V )   z =  z = −   + 
 r 3   z r  

 Vz 2 (    V 1 Vr 
 zz = −   2 −   V )  r =   r = −   r ( )+
 z 3   r r r  
  1 V Vr  2   V V 
  = − 2  +  −  (  V )  zr =  rz = −   z + r 
  r  r  3   r z 

1  1  V
V = (rVr ) + (V ) + z
r r r  z
Coordenadas cilíndricas: Direcciones r ,  y  respectivamente.

 Vr 2 (    V 1 Vr 
 rr = −   2 −   V )  r =   r = −   r ( )+
 r 3   r r r  

  1 V Vr  2 (   1 Vr  V 
  = −   2  +  −   V )   r =  r = −   + r (  )
  r  r  3   rsen  r r 

  1 V Vr V cot   2   sen  V 1 V 

  = − 2  + +  −  (  V )   =   = −   )+
rsen  
(
  rsen  r r  3   r  sen

1  2 1  1 V
V = (r Vr ) + (V sen ) +
r r
2
rsen  rsen 

Para un fluido newtoniano incompresible, el segundo miembro de la ecuación (1.113) es nula, por
lo cual, podemos escribir el modelo de viscosidad como:

 = −V (7.114)

Luego tenemos el producto de la gradiente por el tensor, es decir

. = −.V = −  2V (7.115)

Lo que al reemplazar en la ecuación (7.103) conduce a la obtención de la ecuación de Navier-

Stokes dada por (7.113)

DV
 = −p + 2V +  g (7.116)
Dt

Vx V V V p  2V  2V  2V
( + Vx x + Vy x + Vz x ) = − +  ( 2x + 2x + 2x ) +  g x (7.117)
t x y z x x y z
 Vy Vy Vy
p Vy
 2Vy  2Vy  2Vy
( + Vx + Vy + Vz ) = − −  ( 2 + 2 + 2 ) +  g y (7.118)
t x y z y x y z

Vz Vz Vz Vz p  2Vz  2Vz  2Vz

( + Vx + Vy + Vz ) = − −  ( 2 + 2 + 2 ) +  g z (7.119)
t x y z z x y z

Coordenadas cilíndricas: Direcciones r ,  y z respectivamente.

Vr V V Vr V2 V p    1  (rVr )  1  2Vr 2 V  2Vr 

( + Vr r +  − + Vz r ) = − +    + 2 − 2 + 2  +  gr
t r r  r z r  r  r r  r 
2
r  z 

(7.120)

V V V V V V V 1 p    1  (rV )  1  2V 2 Vr  2V 

( + Vr  +   + r  + Vz  ) = − −  + 2 + 2 + 2  +  g
t r r  r z r   r  r r  r 
2
r  z 

(7.121)

Vz V V Vz V p    1  (rVz )  1  2Vz  2V 

( + Vr z +  + Vz z ) = − −    + 2 + 2  +  g z (7.122)
t r r  z z  r  r r  r  2
z 

Nota: El termino V2 / r es la fuerza centrífuga. Corresponde a la fuerza efectiva en la dirección

r que resulta del movimiento en la dirección  . El termino VrV / r es la fuerza de Coriolis; es
una fuerza efectiva en la dirección  cuando existe flujo en ambas direcciones r y  e interviene
en problemas de flujo en las inmediaciones de un disco que gira.

Coordenadas esféricas: Direcciones r ,  y  , respectivamente

Vr Vr V Vr V Vr V2 + V2 p 2 2 V

( + Vr + + − ) = − +  ( 2Vr − 2 Vr − 2 
t r r  rsen  r r r r 

2 2 V
− V cot  − 2 2
2 
) +  gr (7.123)
r r sen  

V V V V V V VrV V2 cot  1 p 2 V

( + Vr + + + − )=− +  ( 2V + 2 r
t r r  rsen  r r r  r 

V 2cos  V
− − 2 2 ) +  g (7.124)
r sen  r sen  
2 2

V V V V V V VVr VV 1 p

( + Vr + +  + + cot  ) = − +  ( 2V
t r r  rsen  r r rsen 

V 2 Vr 2cos  V

− + + ) +  g (7.125)
r sen 
2 2
r sen   r 2 sen 2 
2 2
 1  2  1   1 2
2 = ( r ) + ( sen ) + (7.126)
r 2 r r r 2 sen   r 2 sen2  2
7.8.4 TRANSPORTE DE ENERGIA EN FUNCION A LA DENSIDAD DE FLUJO DE CALOR Y DENSIDAD
DE FLUJO VISCOSO DE CANTIDAD DE MOVIMIENTO

Coordenadas rectangulares

T T T T q q q     Vx Vy Vz 

 CV ( + Vx + Vy + Vz ) = −( x + y + z ) − T   + + 
t x y z x y z  T   x y z 


 V V V     V V   V V  
 Vy Vz   
−  xx x +  yy y +  zz z   −  xy  x + y  +  xz  x + z  +  yz  +  

 x y z    y x   z x   z y   


Coordenadas cilíndricas

T T V T T 1 (rqr ) 1 q qz    1  1 V Vz

 CV ( + Vr + + Vz ) = −( + + ) −T  ( (rVr ) + + )
t r r  z r r r  z  T  r r r  z
 V 1 V V     V 1 Vr   Vz Vr   1 Vz V 
−  rr r +   ( + Vr ) +  zz z  −  r  r r ( r ) + r   +  rz  r + z  +  z  + 
 r r  z       r  z 

Coordenadas esféricas

T T V T V T 1 (r 2 qr ) 1 ( sen q ) 1 q

CV ( + Vr + +  ) = −( 2 + + )−
t r r  rsen  r r rsen  rsen 

   1  2 1  ( senV ) 1 V  V 1 V Vr 1 V Vr V cot  

T  ( 2 (r Vr ) + + ) −  rr +   ( + ) +   ( + + )
 T  p r r rsen  rsen   r r  r rsen  r r 


  V 1 Vr V   V 1 Vr   1 V 1 V cot    
−  r +
 r r  − +   +  +    + − V  
r 

 r rsen   r  rsen 
r

   r 
7.8.5 COMPONENTES DE LA DENSIDAD DE FLUJO DE CALOR

FORMA MACROSCOPICA A TRAVES DE ECUACIONES INTEGRALES


Vdv = −  V (V .n)dA +  F
t 
(7.127)
v A

En esta ecuación se observa que la acumulación de cantidad de movimiento viene a ser la entrada
de cantidad de movimiento menos la salida de cantidad de movimiento mas la suma de fuerzas.
Las componentes de las fuerzas son:

Fuerza corporal o gravitacional o volumétrica: Fg

Fuerzas superficiales o de presión: Fp

Fuerza de fricción: Fs

Fuerza ejercida por la superficie solida sobre el fluido: B

7.9 TABLAS DE LAS ECUACIONES DE TRANSPORTE DE ENERGIA EN LOS

DIFERENTES SISTEMAS DE COORDENADAS.

7.10 TABLAS DE LAS ECUACIONES DE TRANSPORTE DE MATERIA EN

LOS DIFERENTES SISTEMAS DE COORDENADAS.
Analogías de los fenómenos de transporte
[Link]
perfiles de velocidad
[Link]