Viii Anexos Ecuaciones Auxiliares para Fenomenos de Transporte
Cargado por
Victor Alca SalcedoViii Anexos Ecuaciones Auxiliares para Fenomenos de Transporte
Cargado por
Victor Alca SalcedoUNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
FENOMENOS DE TRANSPORTE
VIII ANEXOS
ECUACIONES AUXILIARES PARA FENOMENOS DE TRANSPORTE
Pag
8.1 Introducción al algebra vectorial 03
8.1.1 Las magnitudes físicas 03
8.1.2 Producto vectorial y magnitud física del resultado 03
8.1.3 Representación de las magnitudes físicas desde el punto de vista del algebra 05
vectorial
8.1.4 Operaciones con vectores unitarios 05
8.2 Longitudes, áreas y volúmenes en los diferentes sistemas de coordenadas 10
8.2.1 Sistema de coordenadas rectangulares 10
8.2.2 Sistema de coordenadas polares 11
8.2.3 Sistema de coordenadas cilíndricas 12
8.2.4 Sistema de coordenadas esféricas 14
8.2.5 Sistema de coordenadas ortogonales generales 17
8.2.6 Aplicaciones del sistema de coordenadas ortogonales generalizado 18
8.3 Algebra de vectores y tensores en coordenadas curvilíneas 23
8.3.1 Transformación de coordenadas rectangulares a cilíndricas 23
8.3.2 Relación entre vectores unitarios en coordenadas cilíndricas 24
8.3.3 Transformación de coordenadas rectangulares a esféricas 25
8.3.4 Relación entre vectores unitarios en coordenadas esféricas 25
8.4 Operaciones diferenciales en coordenadas curvilíneas 25
8.4.1 Derivada de vectores unitarios en coordenadas cilíndricas 25
8.4.2 Derivada de vectores unitarios en coordenadas esféricas 26
8.4.3 Otras ecuaciones de interés 30
8.5 Teoremas integrales para vectores y tensores 31
8.5.1 Teorema de Green 31
8.5.2 Teorema de Stokes 34
8.5.3 Teorema de Gauss 37
8.5.4 Resumen de los teoremas integrales 42
8.6 Valores promedio de funciones 43
8.6.1 Valor medio respecto al tiempo 43
8.6.2 Valor promedio de una función F a lo largo de una trayectoria 44
8.6.3 Valor promedio de una función F sobre una superficie A 44
8.6.4 Valor promedio de una función F dentro de un volumen v 44
8.7 El teorema de Reynolds 45
8.7.1 Primera forma 45
8.7.2 Segunda forma 47
8.7.3 Tercera forma 50
8.8 Tablas de las ecuaciones de transporte de cantidad de movimiento energía y 54
materia en su forma macroscópica.
8.8.1 Balance macroscópico de cantidad de movimiento 55
8.8.2 Balance macroscópico de energía 58
1
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
8.8.3 Balance macroscópico de materia 60
8.9 Tablas de las ecuaciones de transporte de cantidad de movimiento energía y 61
materia en su forma macroscópica
8.9.1 Ecuación de conservación de cantidad de movimiento 61
8.9.2 La ecuación de conservación de energía 65
8.9.3 Transporte de materia en los diferentes sistemas de coordenadas 68
8.10 Analogías de los fenómenos de transporte difusional 69
2
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
8.1.- Introducción al algebra vectorial
8.1.1 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 y de color azul)
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.
8.1.2 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
Fuente: (Bird, et al., 1993; Bird, et al., 2006)
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.
3
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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.
4
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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.
8.1.3 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
i =1
; F = F = k Fk
k =1
; = j
j =1 x j
3 3 3 3
Tensor: =
i =1 j =1
i j ij ; V = i
i =1 j =1 xi
jV j
8.1.4 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
5
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Figura 1: Esquema para la determinación del signo del tensor de alternancia
Fuente: Elaboración propia
+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
j =1 k =1
ijk hjk = 2 ih y
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.
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 equivalentes 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 (8.1)
0 si i j
δi × δ j = ε ijk δk (8.2)
k
δi δ j : δk δl = (δ j .δk )(δi .δl ) = δ jk δil (8.3)
i j . k = i ( j . k ) = i jk (8.4)
6
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
i . j k = ( i . j ) k = ij k (8.5)
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:
7
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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:
8
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
: V = ij i j : kVk l = i j : k l ijVk
i j k l xl i j k l xl
1 si l = i y k = j
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.
9
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
8.2. Longitudes, áreas y volúmenes en los diferentes sistemas de coordenadas
8.2.1 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 (8.6)
El área diferencial de cada uno de los lados del cubo unitario se obtiene mediante
dAx = dydz, dAy = dxdz, dAz = dxdy (8.7)
El volumen diferencial del cubo unitario está dado por:
dv = dxdydz (8.8)
Figura 2: Sistema de coordenadas rectangulares
Fuente: (Bird, et. al. 1993)
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 áreas 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 Ay = dxdz = x 0 z 0 = 4 1 = 4u 2
2 4 2 1 4
0 0 0 0
2 1
Az = dxdy = x 0 y 0 = 1 2 = 2u 2
1 2
0 0
10
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
8.2.2 Sistema de coordenadas polares
Las coordenadas polares o sistema de coordenadas polares son un sistema de
coordenadas bidimensional en el que cada punto del plano se determina por una distancia y
un ángulo. Este sistema es ampliamente utilizado en física y trigonometría.
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
3.
Figura 3: Sistema de coordenadas polares en el plano
Fuente: Elaboración propia
Si el punto P tiene coordenadas ( x, y ) , con x 0 , entonces, sus coordenadas polares en (r , )
se obtiene como
r= x2 + y2 (8.9)
11
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
= arctan( y / x ) = tan −1 ( y / x) (8.10)
La transformación inversa de coordenadas polares a rectangulares es:
x = r cos (8.11)
y = r sin (8.12)
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
8.2.3 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 4: Esquema del sistema de coordenadas cilíndricas: (a) Representación de los ejes
cartesianos y cilíndricos; (b) Vista isométrica de las coordenadas cilíndricas; (c) Elementos
diferenciales de volumen, de línea y de superficie.
12
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
(a) (b)
(c)
Fuente: (Bird et. al 2006; Tosun, 2006)
Equivalencia entre coordenadas cilíndricas y rectangulares
El eje de coordenadas está dado por:
x = r cos y = rsen z=z (8.13)
Con r 0 y 0 2 , − z
r = + x 2 + y 2 , = arctan( y / x ) z=z (8.14)
Las diferenciales de longitud, en este caso están dadas por:
ds1 = dr ds2 = rd ds3 = dz (8.17)
La diferencial del área lateral está dada por:
dAL = (rd )(dz ) = rd dz (8.18)
La diferencial del área de la base viene dada por:
dAB = (rdr )(d ) = rdrd (8.19)
13
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
La diferencial de un plano vertical paralelo a la dirección longitudinal inscrito dentro del cilindro
es:
dAV = drdz (8.20)
La diferencial de volumen
dv = (dr )(rd )(dz ) = rdrd dz (8.21)
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 (8.18)
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 (8.19)
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 , , z )drdz
0
z1 r1
2 r2
Sobre la superficie z = z0 A(r , , z0 ) = f (r, , z )rdrd
0
1 r1
c) Para encontrar el volumen integramos la ecuación (8.21)
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
8.2.4 Sistema de coordenadas esféricas
14
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 5 se muestra el esquema para la
transformación de sistema de coordenadas cartesianas a esféricas y viceversa.
Figura 5: Esquema del sistema de coordenadas esféricas: (a) Representación de los ejes
cartesianos y esféricos; (b) Vista isométrica de las coordenadas esféricas; (c) Elementos
diferenciales de volumen, de línea y de superficie.
(a) (b)
(c)
La equivalencia entre coordenadas esféricas y rectangulares se muestran en las ecuaciones
subsiguientes.
x = rsen cos y = rsen sen z = r cos (8.22)
Con r 0 , 0 , 0 2
r = + x2 + y2 + z2 (8.23)
x2 + y2
= arctan (8.24)
z
y
= arctan (8.25)
x
Las diferenciales de longitud, en este caso están dadas por:
15
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
ds1 = dr ds2 = rd ds3 = rsen d (8.26)
La diferencial del área lateral está dada por:
dA = (rd )(rsen d ) = r 2 sen d d (8.27)
La diferencial de volumen
dv = (dr )(rd )(rsen d ) = r 2 sen drd d (8.28)
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 (8.27)
2 2
A= dA = r sen d d = 2 R sen d = 4 R
2 2 2
=0 =0 =0 =0 0
Generalizando
2 2
Sobre la superficie r = r0 A(r0 , , ) = f (r , , )r sen d d
0 0
2
1 1
2 r2
Sobre la superficie = 0 A(r , 0 , ) = f (r, , )sen rdrd
0 0
1 r1
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 (8.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
Tabla 1:
Cartesiana Cilíndrica Esférica
Cartesiana x = (−, +) x = r cos x = cos sen
y = (−, +) y = rsen y = sen sen
z = (−, +) z=z z = cos
Cilíndrica r = (0, r = sen
r = x2 + y 2
= (0, 2 ) =
z = cos
16
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
= tan −1 ( y / x) z = −, +
z=z
Esférica = x2 + y 2 + z 2 = r2 + z2 = 0, 2 )
= = (0, 2 )
= tan −1 ( y / x)
= (0, )
= tan −1 (r / z )
= tan −1 ( x 2 + y 2 / z )
Nota: Para diferenciar los radios del cilindro y la esfera, en el primer caso se usa r para
denominar el radio y en el segundo caso para su denominación.
8.2.5 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 6.
Figura 6: Esquema del sistema de coordenadas ortogonales generalizado
Fuente: Elaboración propia
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:
17
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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.
Geometría 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
Fuente: (Theodore,1995)
8.2.6 Aplicaciones del sistema de coordenadas ortogonales generalizado
Gradiente de un escalar
s s s
s = 1h1 + 2 h2 + 3h3 (8.29)
Ejemplo 14
Haciendo uso de la ecuación (8.29) y las equivalencias de la tabla 1 encuentre s en los tres
sistemas de coordenadas.
Solución:
En la ecuación (8.29) se reemplaza los valores respectivos de sus parámetros, haciendo uso de la
tabla 2
a) Coordenadas rectangulares
18
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 + + (8.30)
h2 h3 h1h3 h1h2
Ejemplo 15
Haciendo uso de la ecuación (8.25) y sus equivalencias dadas en la tabla 2 encuentre .V en los
tres sistemas de coordenadas.
Solución:
En la ecuación (8.30) 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 (8.31)
V = 1h3h2 3 − 2 + 2 h1h3 1 − + 3h2 h1 −
h3 h2 h1 h3 h2 h1
Ejemplo 16
Haciendo uso de la ecuación (8.31) y sus equivalencias dadas en la tabla 2 encuentre V en
los tres sistemas de coordenadas.
19
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 + + (8.32)
h2 h3 h1h3 h1h2
Ejemplo 17
Haciendo uso de la ecuación (8.32) 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 2 + 2
r r r r z
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
20
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
2 2 2
d d d
ds =
2
+ + (8.33)
h1 h2 h3
Ejemplo 18
Haciendo uso de la ecuación (8.33) 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 (8.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 2
ds 2 = ( ) +( ) +( ) = (dr )2 + (rd )2 + (rsen d )2
1 1/ r 1/ rsen
Elemento de área diferencial
d d
dA = (8.34)
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 (8.34) 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 (8.29) y los parámetros de la tabla 2 y se obtiene:
21
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 (8.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 = (8.35)
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 (8.35) y haciendo uso de los parámetros de la tabla 2 se tiene:
dr d dz
dv = = rdrd dz
1 1/ r 1
22
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ejemplo 24
Encuentre la diferencial de volumen de una esfera de radio r
Solución:
A partir de la ecuación (8.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
8.3 Algebra de vectores y tensores en coordenadas curvilíneas
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 (8.36)
y,z y,z ,z y,z r ,z y,z r ,
8.3.1 Transformación de coordenadas rectangulares a cilíndricas
x = r cos y = rsen z=z (8.37)
r = + x2 + y2 = arctan( y / x) z=z
(8.38)
Al aplicar la regla de la cadena, se obtiene
sen
= (cos ) + (− ) + (0) (8.39)
x r r z
cos
= ( sen ) + ( ) + (0) (8.40)
y r r z
= (0) + (0) + (1) (8.41)
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.
23
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Fuente: (Bird et al., 1993)
8.3.2 Relación entre vectores unitarios en coordenadas cilíndricas
De coordenadas rectangulares a coordenadas cilíndricas
r = (cos ) x + ( sen ) y + (0) z (8.42)
= (− sen ) x + (cos ) y + (0) z (8.43)
z = (0) x + (0) y + (1) z (8.44)
De coordenadas cilíndricas a coordenadas rectangulares
x = (cos ) r + (− sen ) + (0) z (8.45)
y = ( sen ) r + (cos ) + (0) z (8.46)
z = (0) r + (0) + (1) z (8.47)
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
24
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
cos
= ( sen ) + ( ) + (0) (c)
y r r z
= (0) + (0) + (1) (d)
z r z
Por otro lado, la relación entre vectores unitarios está 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
8.3.3 Transformación de coordenadas rectangulares a esféricas
x = rsen cos y = rsen sen z = r cos (8.48)
x2 + y2 y
r =+ x + y +z
2 2 2
= arctan = arctan (8.49)
z x
A partir de la regla de la cadena se obtiene:
cos cos sen
= ( sen cos ) + ( ) + (− ) (8.50)
x r r rsen
cos sen cos
= ( sen sen ) + ( ) +( ) (8.51)
y r r rsen
sen
= (cos ) + (− ) + (0) (8.52)
z r r
7.3.4 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 (8.53)
= (cos cos ) x + (cos sen ) y + (−sen ) z (8.54)
= (− sen ) x + (cos ) y + (0) z (8.55)
25
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
De coordenadas esféricas a coordenadas rectangulares
x = ( sen cos ) r + (cos cos ) + (−sen ) (8.56)
y = ( sen sen ) r + (cos sen ) + (cos ) (8.57)
z = (cos ) r + (− sen ) + (0) (8.58)
Siguiendo el mismo procedimiento anterior, para coordenadas esféricas se obtiene:
1 1
= r + +
r r rsen
8.4. Operaciones diferenciales en coordenadas curvilíneas
8.4.1 Derivada de vectores unitarios en coordenadas cilíndricas
r z
=0 =0 =0 (8.59)
r r r
r z
= = − r =0 (8.60)
r z
=0 =0 =0 (8.61)
z z z
8.4.2 Derivada de vectores unitarios en coordenadas esféricas
r
=0 =0 =0 (8.62)
r r r
r
= = − r =0 (8.63)
r
= sen = cos = − sen r − cos (8.64)
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
26
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 = 2 + (V sen ) +
r r rsen rsen
Ejemplo 27
Encontrar V en coordenadas cilíndricas
Solución:
27
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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:
28
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
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 .
29
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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:
30
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
1 1 r 1 r +
. = 2 (r 2 rr ) + r cot + + −
r r r r r rsen r
8.4.3 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 ) (8.65)
x y z
Coordenadas cilíndricas
Vr 1 Vr V V V 1 V Vr
( : V ) = rr ( ) + r ( − ) + rz ( r ) + r ( ) + ( + )
r r r z r r r
V V 1 Vz V
+ z ( ) + zr ( z ) + z ( ) + zz ( z ) (8.66)
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 ) (8.67)
r r rsen r r
8.5.1 Teorema de Green
2
Si F es un campo vectorial continuo y diferenciable en R , 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
31
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Fuente: Elaboración propia
N M
C
F .dr = (
S
−
x y
)dA (8.68)
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 y y = x en sentido
2
contrario a las agujas del reloj.
Solución:
La figura correspondiente es:
Figura E_8.29: Superficie de aplicación de Teorema de Green
Fuente: Elaboración propia
Aplicamos la ecuación correspondiente al teorema de Green
32
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 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 en la
hoja derecha de la lemniscata.
Figura E_8.30: Superficie de aplicación del Teorema de Green
33
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
= cos 2
1
2
Fuente: Elaboración propia
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 cos y = 1
x y x y
El resultado es un vector en la dirección k; sin embargo, dado que la función está 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.
8.5.2 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
3
tienen derivadas parciales continuas en una región R que contiene a S. Entonces:
C
F .dr = rotF . dS
S
(8.69)
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
34
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ejemplo 31
Verificar que sí 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 9: Esquema para la aplicación del Teorema de Stokes
Fuente: Elaboración propia
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 está orientado hacia abajo (en contra de las
manecillas del reloj, el flujo mayoritariamente sale del cono)
Figura 10: Superficie de circulación para la aplicación del Teorema de Stokes
Fuente: Elaboración propia
35
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 ) = −4senti − 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 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 E_8.32: Región R para la aplicación del Teorema de Stokes
36
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 este 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
37
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
x i y j
rotF . nds = ((0) i + (0) j +(2 k )).(
S R 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.
8.5.3 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 (8.70)
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 (8.71)
En la práctica, es mucho más 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
(8.72)
Ejemplo 33
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 xy1 − 5xz 23 . Encuentre el flujo volumétrico utilizando la
integral cerrada.
38
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Solución:
El caudal o flujo volumétrico se calcula mediante:
Q= S
. )ds =
(nV V .d A
S
Figura E_8.32: Elemento de volumen para la aplicación del Teorema de Gauss
dAxz
Fuente: Elaboración propia
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 = − 3dxdz (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
1
Qx , y ,0 = (10 xy1 − 5 xz 3 ).( 3dxdz ) =
2
( − 5 x)dx = −2.5 x 2 = −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
39
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 xy1 − 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 34
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 xy1 − 5xz 23 . 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 xy1 − 5 xz 2 3 ) = + +
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
(8.73)
(. )dv =
v s
n. ds (8.74)
Nuevamente estas ecuaciones tienen gran aplicabilidad al resolver la integral volumétrica, en
lugar de la integral de superficie, debido a su simplicidad.
40
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
El teorema indica lo siguiente: si F es un campo vectorial continuo y diferenciable en R 3 , 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
(8.75)
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.
Figura 11: Campo de flujo y dirección de desplazamiento del fluido
ds
ds
v
En base a la Figura 11 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.
Figura 12: Representación de los infinitos diferenciales de volumen dentro del campo de flujo
vectorial
ds
ds
Fuente: Elaboración propia
41
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ejemplo 35
Obtener el flujo del campo F = y i + xy j − z k a través de la superficie cerrada formada por
el cilindro x + y 4 , el plano z = 0 y el paraboloide z = x + y .
2 2 2 2
Solución:
El esquema es el siguiente:
Figura E_8.35: Representación esquemática de las ecuaciones que caracterizan las curvas
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 + ( sen ) 2 = 2 por otro lado, del resultado del integrando se
tiene: x − 1 = cos − 1
Luego, la integral se transforma en:
42
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
8.5.4 Resumen de los teoremas integrales
Teorema fundamental del calculo
b
d
dx
a
f ( x )dx = f (b) − f (a ) (8.76)
Teorema para integrales de línea
b
f ( x) dr = f ( x).dr = f (r (b)) − f (r (a))
C a
(8.77)
Teorema de Green
F ( x)dA = F ( x) dr (8.78)
Teorema de Gauss
F ( x )dv = F ( x) dS (8.79)
Teorema de Stokes
( F ( x )).dS = F ( x ) dr (8.80)
r : vector de parametrización f ( x) : Función primitiva S : Vector superficie
v : Volumen de control
8.6. Valores promedio de funciones
8.6.1Valor medio respecto al tiempo
43
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
(8.81)
dt dt (t1 − t0 )
t1
t
t0
Ejemplo 36
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 /2
(4t − 1.2cos t − 0.8sent ) 0
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
8.6.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.82)
ds ds
L
s
0
8.6.3 Valor promedio de una función F sobre una superficie A
FdA FdA
F = A
= A
(7.83)
dA A
A
44
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ejemplo 37
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
8.6.4 Valor promedio de una función F dentro de un volumen v
El valor medio de una función contenido en un volumen se calcula mediante la siguiente relación:
F =
Fdv = Fdv
v v
(8.84)
dv v
v
Ejemplo 38
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
45
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
8.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ículas; 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 conservación 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.
8.7.1 Primera forma
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 13: Representación del sistema y volumen de control
En el instante t se tiene:
t
Bsist = BVct (8.85)
Ahora cuando el tiempo es t + t
Figura 14: Desplazamiento del sistema debido al flujo
El sistema se desplaza (color rosado), mientras el volumen de control permanece fijo (color azul)
46
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 15: Esquema que relaciona el sistema y el volumen de control
Por tanto, en base a la figura se tiene:
t +t
Bsist = BVct +t − BIt +t + BIIt +t (8.86)
Restando la ecuación (8.2) de la ecuación (8.1)
t +t
Bsist − Bsist
t
B t +t − BVct B t +t B t +t
lim = lim Vc − I + II (8.87)
x →0 t x →0 t t t
Aplicando la derivada correspondiente se tiene
dBsist dBVc BIt +t BIIt +t dBVc • •
= − + = − BI + BII (8.88)
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 Qt
= = = = VA (8.89)
t t t t
De esta manera, la ecuación (3) escrita en términos de la masa, permite obtener
dmsist dmVc • •
= − m I + m II (8.90)
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 • •
= − m I + m II (8.91)
dt
47
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 (8.92)
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
= (8.93)
V2 d1
8.7.2 Segunda forma
En esta segunda parte encontraremos una ecuación más 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 16: 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 cauce 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
48
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Q= V .dA = V . ndA = V A
sc sc
(8.94)
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 (8.95)
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 ser 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 17: Volumen de control, superficie de control y sistema
49
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 = dt Vc dv + Sc
(V . n)dA (8.96)
sist
Si aplicamos esta ecuación a la conservación de materia = 1 , luego la ecuación (8.12) se
transforma en:
dm d
dt = dt Vc dv + Sc
(V . n)dA (8.97)
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 (8.98)
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 (8.99)
Esta ecuación se reduce a: 2V2 A2 − 1V1 A1 = 0 , que es la ecuación de continuidad macroscópica
para referido a la conservación de masa total en régimen estacionario.
7.7.3 Tercera forma
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 18: La evolución del sistema y el volumen de control
50
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
(8.100)
t →0 t t
dB ( dv ) t +t − ( dv ) t dv ) t +t
dt = II II
+ III
(8.101)
sist t t
En la ecuación (8.87) por aplicación del teorema de Gauss, el último 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 (8.102)
sist
Si consideramos que B sea la masa y = B / m , entonces la ecuación (8.89) se transforma en:
dm
d
dt = dt Vc dv + Sc
(V . n)dA (8.103)
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 (8.104)
Sc
La ecuación (7.90) representa la ecuación de conservación 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
51
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 19: Vista simplificada de un sistema en movimiento
dB
=
DB
=
III dv + II dv )t +t − (I dv + II dv )t ) (8.105)
lim(
dt sist Dt t →0 t
( dv − dv ) dv dv
DB II II III I
= lim t + t t
+ lim t + t
− lim t
(8.106)
Dt t →0 t t → 0 t t → 0 t
Luego estas integrales deben ser tratadas por partes:
( dv − dv )
d
dt
II II
lim t + t t
= dv (8.107)
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
52
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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 límite de (8.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 (8.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
(8.108)
t →0 t t →0 t
La ecuación (8.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 (8.109)
Figura 20: 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
(a) (b)
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 21: Superficie de control que muestra el flujo de entrada de masa
53
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
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
(8.110)
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 (8.94) puede expresarse
mediante:
Tasa neta de flujo de salida de B desde VC = ( V . dA)
SC
(8.111)
Ahora utilizando las ecuaciones (8.93) y (8.97) se obtiene:
DB
Dt
= ( V . dA) + t dv
SC vc
(8.112)
La ecuación (7.98) se conoce como el Teorema de Reynolds. Para el caso de transporte de masa
se tiene = 1 .
54
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
8.8 Tablas de las ecuaciones de transporte de cantidad de movimiento energía y materia
en su forma macroscópica.
Los balances macroscópicos de cantidad de movimiento, energía y materia reciben esa
denominación debido a que no contienen el termino de transporte molecular; en el primer
caso por diferencia de velocidades, en el segundo caso por diferencia de temperaturas y
en el tercer caso por diferencia de concentraciones; solo consideran el transporte debido
al flujo.
Figura 22: Esquema para la aplicación de los balances de cantidad de movimiento, energía
y materia
A : Área del elemento de control con su respectiva área diferencial
v : Volumen de control con su respectivo volumen de control
Donde cos dA es el área diferencial proyectada en la dirección perpendicular al vector V y
es el ángulo entre el vector velocidad y el vector n en la dirección exterior, tal como se observa
en la Figura.
8.8.1 Balance macroscópico de cantidad de movimiento
La ecuación general está dada por:
Acumulación de C.M. = Entrada de C.M. - Salida de C.M. + Generación de C.M.
Vdv = − V (V . n)dA + F
t
(8.113)
vc SC
Esta ecuación puede ser escrita para cada uno de sus componentes cartesianos:
Dirección x
Vx dv = − Vx (V . n)dA + Fx = − VxV cos dA + Fx
t
(8.114)
vc SC SC
Dirección y
55
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
V y dv = − V y (V . n)dA + Fy = − V yV cos dA + Fy
t
(8.115)
vc SC SC
Dirección z
Vz dv = − Vz (V . n)dA + Fz = − VzV cos dA + Fz
t
(8.116)
vc SC SC
C. M. Cantidad de movimiento
F = F p + Fg + Fs + B (8.117)
Fp : Fuerzas superficiales o fuerzas de presión; causada por las fuerzas de presión que actúan
sobre la superficie del sistema.
Fg : Fuerzas corporales o volumétricas; causada por la fuerza del campo gravitacional
Fs : Fuerzas de fricción; debido al rozamiento entre un fluido y la superficie del solido
VNCM : Es la variación neta de cantidad de movimiento debido al flujo
F : Es la suma de fuerzas actuando sobre el fluido
B : Fuerza de la superficie solida ejercida sobre el fluido; en los casos en los cuales la superficie
de control pasa por una superficie sólida, esta presenta una fuerza B , que viene a ser la resultante
de las fuerzas que actúan sobre el volumen de control en dichos puntos.
R : Es la fuerza de reacción del fluido sobre la superficie interna del conducto
En base a la ecuación (8.116), la variación neta de cantidad de movimiento en la dirección x , en
régimen estacionario, está dada por:
− Vx (V . n)dA = VNCM x = m(V2 cos 2 − V1cos1 )
SC
Asimismo, la suma de fuerzas se obtiene mediante:
F x = p1 A1 cos 1 − p2 A2 cos 2 + Bx
Por otro lado:
VNCM x = Fx
Bx = m(V2 cos 2 − V1 cos 1 ) + p2 A2 cos 2 − p1 A1 cos 1
Luego:
Rx = − Bx
56
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ahora en la dirección y
− Vy (V . n)dA = VNCM x = m(V2 sen 2 − V1sen1 )
SC
Asimismo, la suma de fuerzas se obtiene mediante:
F y = p1 A1sen1 − p2 A2 sen 2 + Bx − mg
Por otro lado:
VNCM y = Fy
By = m(V2 sen2 − V1sen1 ) + p2 A2 sen2 − p1 A1sen1 + mg
Luego:
Ry = − By
La fuerza total resultante se obtiene mediante:
R = Rx 2 + Ry 2
La dirección de la fuerza se obtiene a partir de las fuerzas resultantes en cada una de las
direcciones de los ejes de coordenadas.
= arctan( Ry / Rx )
Nota: La fuerza de fricción por lo general no se toma en cuenta debido a su contribución pequeña;
sin embargo, debido a que estos balances macroscópicos se llevan a cabo en régimen turbulento,
solo sería necesario conocer el perfil de velocidad. Este perfil se usa para estimar el esfuerzo de
corte en la pared y con el conocimiento del área de contacto del fluido con el sólido que lo
contiene, es posible obtener la fuerza de rozamiento o fricción en la pared.
Para el caso de la fuerza de presión, es necesario relacionarlo con la velocidad; esta relación se
obtiene haciendo uso de la ecuación de Bernoulli.
La fuerza volumétrica se obtiene mediante el conocimiento del peso del fluido y su contribución
vertical. Si el flujo fuera horizontal, no se toma en cuenta este término.
La fuerza que ejerce la superficie solida se obtiene a partir del balance macroscópico de cantidad
de movimiento.
Por ejemplo, si quisiéramos observar cómo varia la velocidad V a la salida de la Figura xxx
(tobera convergente divergente horizontal) en régimen no estacionario, conociendo los otros
parámetros de la ecuación (7.102), esta ecuación se reduce a:
( Vz v )
= AV
1 1 − A2V2 + P1 A1 − P2 A2 + Bz
2 2
(8.118)
t
Se ha eliminado la fuerza debido al campo gravitacional por tratarse de una tobera horizontal; se
ha eliminado la también la fuerza debido a la fricción por ser su contribución muy pequeña; sin
57
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
embargo, se puede estimar este valor a partir de la fuerza de rozamiento en la pared y el área de
rozamiento, o a partir de la diferencia de presiones y el área de flujo.
En (7.104) hacemos que V2 sea la variable de salida y que V1 sea conocida.
Vz
v = AV
1 1 − A2Vz + P1 A1 − P2 A2 + Bz
2 2
t
Vz AV 2
AV2 PA − P A B
= 1 1 − 2 z + 1 1 2 2+ z
t v v v v
Esta ecuación tiene la forma
Vz
= −MVz2 + N
t
La misma que puede ser integrada fácilmente
8.8.2 Balance macroscópico de energía
Figura 23: Esquema de un sistema termodinámico
Consideremos el siguiente sistema termodinámico
Suponiendo que se suministra calor y se realiza trabajo sobre el sistema; por convención, tanto el
calor como el trabajo serán positivos; por tanto, una de las formas del balance de energía
mecánico, resulta:
Tasa de acumulación de energía = Tasa de entrada de energía-Tasa de salida de energía + Tasa
de generación de energía.
E
= ( Ee + Q + W ) − Es (8.119)
t
E = m E (8.120)
Las letras que llevan un punto sobre si, están expresadas en unidades de energía por unidad de
tiempo.
Las letras que llevan una línea sobre si indican unidades de energía por unidad de masa.
Las letras que no llevan ningún símbolo están dadas en unidades de energía.
Nota:
a) La ecuación de balance de energía mecánica por lo general no contiene al termino de
generación de energía
58
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
b) El calor puede ingresar al sistema, por diferentes medios (Una chaqueta externa, un
serpentín, una resistencia u otro medio)
c) La tasa de generación de calor debido a la fricción se considera despreciable.
La ecuación (8.108) se puede reescribir como
E
= Ee − Es + Q + W (8.121)
t
Por otro lado, la energía transportada por la materia se puede expresar del siguiente modo:
E = U + Ec + E p (8.122)
Es decir, es la contribución de la energía interna ( U ), la energía cinética ( Ec ) y la energía
potencial ( E p ).
Ahora expresando la ecuación (8.110) en términos de la energía interna, cinética y potencial
E
= (U + Ec + E p )m − (U + Ec + E p )m + Q + W
t
(8.123)
e s
La tasa de trabajo realizado por los alrededores sobre el sistema puede ser expresado del siguiente
modo.
W = Ws + ( P v m)e − ( P v m) s (8.124)
El primer término de la derecha representa en trabajo de eje o de flecha (trabajo producido por la
rotación de un eje) y el segundo termino se conoce como flujo de trabajo
v : Volumen especifico del fluido
P : Presión del fluido
Reemplazando (8.113) en (8.112) y haciendo H = U + P v , la ecuación (8.112) resulta:
E
= ( H + Ec + E p )m − ( H + Ec + E p )m + Q + Ws
t
(8.125)
e s
Cuando el cambio de la energía potencial y cinética son pequeños, la ecuación (8.114) se puede
simplificar para régimen estacionario como:
Hm e − Hm s + Q + Ws = 0
(8.126)
A partir de la ecuación (8.112), se deduce la primera ley de la termodinámica para un sistema
cerrado
Ahora, debemos tener en cuenta que, en el caso de ingeniería química, se debe considerar la
energía transportada por la materia, la energía generada y la energía suministrada ya sea en forma
de calor y trabajo. Dependiendo de las aplicaciones, podemos usar alguna de las formas descritas
anteriormente.
Por ejemplo, cuando se trata del transporte de fluido usando una bomba, la ecuación (8.114) se
transforma en:
59
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
( Ec + E p + pv )m − ( Ec + E p + pv )m = Ws (8.127)
s e
Para lograr obtener la ecuación (8.116), se ha omitido el termino no estacionario, además del
término de energía interna, pues en una bomba se trabaja a temperatura constante; también se ha
eliminado el flujo de calor; más bien tendría que agregarse el termino de pérdidas de energía, es
decir
( Ec + E p )m − ( Ec + E p )m + mH f = Ws (8.128)
s e
Donde H f es la perdida de energía por unidad de masa, la cual se expresa mediante
V2
H f = f ( L / d ) (8.129)
2 gc
V : Velocidad media en el conducto
f : Factor de fricción de Moody
L : Longitud de tubería recta o su equivalente debido a sus accesorios y todos de un mismo
diámetro.
Para el caso específico, por ejemplo, de reactores, tanques calentados, enfriados, es mejor
expresarlo el balance de energía en términos de ecuaciones integrales.
Edv = − E (V .n)dA + Gdv + Q + W
t
(8.130)
En esta ecuación para el caso indicado se omite la contribución de las energías potencial y
cinética, así como el trabajo del eje y el flujo de trabajo, de tal manera que la ecuación final
resultante es:
Hdv = − H (V .n)dA + Gdv + Q
t
(8.131)
Donde
G : Unidades de energía por volumen y por tiempo; y para el caso de reactores químicos esta
expresada en términos de la tasa de generación por reacción química, es decir:
G = H R rA (8.132)
Nota: Se debe tener en cuenta que en la literatura se presentan diversas formas del balance de energía; sin embargo,
cada balance es particular que depende de la aplicación especifica y las suposiciones respectivas. Por ejemplo, si en
la ecuación (16) consideramos un agitador de alta velocidad, esa energía del agitador se disipará en forma de calor
dentro del sistema y se adicionará al calor generado por otros mecanismos (calentador) y al calor transportado por
la masa (en forma de entalpia). Si al sistema no se adiciona calor, y además este no está aislado, entonces Q se
expresará en términos de pérdidas de calor (transferencia de calor por convección)
Si la ecuación (8.119) se integra en términos de E , se adiciona Ws y se omite el termino
generación, entonces quedara igual a la ecuación (8.114)
8.8.3 Balance macroscópico de materia
60
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Balance de masa total
dv = − (V .n)dA
t
(8.133)
Ecuación de balance de masa de A para un sistema no reaccionante, donde la concentración A
esta expresada en términos de masa de soluto por volumen de solución
A dv = − A (V .n)dA
t
(8.134)
Ecuación de balance de masa de A para un sistema reaccionante
C A dv = − C A (V .n)dA + rA dv
t
(7.135)
8.9 Tablas de las ecuaciones de transporte de cantidad de movimiento energía y materia
en su forma microscópica
El termino de balance microscópico de los fenómenos de transporte se denomina así,
debido a que, en la deducción de las ecuaciones, también se consideran los términos
debido al gradiente potencial de transferencia. En el caso de transferencia de cantidad de
movimiento es la diferencia de velocidades; en el caso de transferencia de calor, es debido
a la diferencia de temperaturas y en el caso de transferencia de materia debido al gradiente
de concentración (expresado en concentración molar o diferencia de presiones o en
diferencia de composiciones molares)
8.9.1 Ecuación de conservación de cantidad de movimiento
La ecuación de continuidad
Una de las formas de la ecuación de continuidad, las cuales pueden ser escritas en los tres sistemas
de coordenadas.
= − . V (8.136)
t
Ecuación de continuidad en coordenadas rectangulares
( Vx ) ( Vy ) ( Vz )
+ + + =0 (A)
t x y z
Ecuación de continuidad en coordenadas cilíndricas
1 ( rVr ) 1 ( V ) ( Vz )
+ + + =0 (B)
t r r r z
Ecuación de continuidad en coordenadas esféricas
61
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
1 ( r 2Vr ) 1 ( V sen ) 1 ( V )
+ 2 + + =0 (C)
t r r rsen rsen
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 tal como indica la ecuación (7.128)
La ecuación de conservación de cantidad de movimiento
DV
= −p − + g (8.137)
Dt
Coordenadas rectangulares: Direcciones x , y y z respectivamente.
Vx V V V p yx zx
( + Vx x + Vy x + Vz x ) = − − ( xx + + ) + gx (A)
t x y z x x y z
Vy Vy Vy Vy p xy yy zy
( + Vx + Vy + Vz )=− −( + + ) + gy (B)
t x y z y x y z
Vz V V V p yz zz
( + Vx z + Vy z + Vz z ) = − − ( xz + + ) + gz (C)
t x y z z x y z
Coordenadas cilíndricas: Direcciones r , y z respectivamente.
Vr Vr V Vr V2 V p 1 (r rr ) 1 r rz
( + Vr + − + Vz r ) = − − ( + − + ) + gr
t r r r z r r r r r z
(D)
V V V V VrV V 1 p 1 (r 2 r ) 1 z
( + Vr + + + Vz )=− −( 2 + + ) + g
t r r r z r r r r z
(E)
Vz V V Vz V p 1 (r rz ) 1 z zz
( + Vr z + + Vz z ) = − − ( + + ) + gz (F)
t r r z z r r r z
Coordenadas esféricas: Direcciones r , y , respectivamente
Vr V V Vr V Vr V2 + V2 p 1 (r 2 rr ) 1 ( r sen )
( + Vr r + + − ) = − −( 2 +
t r r rsen r r r r rsen
1 r +
+ − ) + gr (G)
rsen r
V V V V V V VrV V2 cot 1 p 1 (r 2 r )
( + Vr + + + − )=− −( 2
t r r rsen r r r r r
62
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
1 ( sen ) 1 r cot
+ + + − ) + g (H)
rsen rsen r r
V V
V V V V VVr V V 1 p 1 (r r )
2
( + Vr + + + + cot ) = − −( 2
t r r rsen r r rsen r r
1 1 r 2cot
+ + + + ) + g (I)
r rsen r r
La 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 tensorial de la ley de viscosidad de Newton:
V j Vi 2 V V y Vz
i j = − ( + ) + ( − k )( x + + ) ij (8.138)
xi x j 3 x y z
La ecuación (8.138) para cada uno de los sistemas de coordenadas, omitiendo el término que
contiene el producto de k por el tensor V ij está dada por:
Coordenadas rectangulares: Direcciones x , y y z respectivamente.
Vx 2 Vx Vy
xx = −2 + ( V ) (A) xy = yx = − + (D)
x 3 y x
Vy 2 V V
yy = −2 + ( V ) (B) yz = zy = − y + z (E)
y 3 z y
Vz 2 Vz Vx
zz = −2 + ( V ) (C) zx = xz = − + (F)
z 3 x z
Vx Vy Vz
V = + + (G)
x y z
Nota: La ecuación (8.138) es equivalente a:
2
i j = − (V + (V )t ) + ( − k )( V ) ij
3
Coordenadas cilíndricas: Direcciones r , y z respectivamente.
Vr 2 ( V 1 Vr
rr = − 2 − V ) r = r = − r ( )+
r
(K)
r 3 r r
(H)
63
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
1 V Vr 2 V 1 Vz
= − 2 + − ( V ) (I) z = z = − +
r r 3 z r
(L)
V 2 Vz Vr
zz = − 2 z − ( V ) zr = rz = − +
z
(M)
z 3 r
(J)
1 1 V
V = (rVr ) + (V ) + z (N)
r r r z
Coordenadas esféricas: Direcciones r , y respectivamente.
Vr 2 ( V 1 Vr
rr = − 2 − V ) r = r = − r ( )+
r
(R)
r 3 r r
(O)
1 V Vr 2 ( 1 Vr V
= − 2 + − V ) (P) r = r = − + r ( ) (S)
r r 3 rsen r r
1 V Vr V cot 2 (Q) sen V 1 V (T)
= − 2 + + − ( V ) = = − )+
rsen
(
rsen r r 3 r sen
1 2 1 1 V
V = (r Vr ) + (V sen ) + (U)
r r
2
rsen rsen
Para un fluido newtoniano incompresible, el segundo miembro de la ecuación (8.138) es nula, por
lo cual, podemos escribir el modelo de viscosidad como:
= −V (8.139)
Luego tenemos el producto de la gradiente por el tensor, es decir
. = −.V = − 2V (8.140)
Lo que al reemplazar en la ecuación (8.128) conduce a la obtención de la ecuación de Navier-
Stokes dada por (7.141)
DV
= −p + 2V + g (8.141)
Dt
Coordenadas esféricas: Direcciones x , y y z , respectivamente
Vx V V V p 2V 2V 2V
( + Vx x + Vy x + Vz x ) = − + ( 2x + 2x + 2x ) + g x (A)
t x y z x x y z
Vy Vy Vy p Vy
2Vy 2Vy 2Vy
( + Vx + Vy + Vz )=− + ( 2 + + 2 ) + gy (B)
t x y z y x y 2 z
64
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Vz V V V p 2V 2V 2V
( + Vx z + Vy z + Vz z ) = − + ( 2z + 2z + 2z ) + g z (C)
t x y z z x y z
Coordenadas cilíndricas: Direcciones r , y z respectivamente.
Vr V V Vr V2 V p 1 (rVr ) 1 2Vr 2 V 2Vr
( + Vr r + − + Vz r ) = − + + 2 − 2 + 2 + gr (D)
t r r r z r r r r r
2
r z
V V V V V V V 1 p 1 (rV ) 1 2V 2 V 2V (E)
( + Vr + + r + Vz ) = − + + 2 + 2 r + 2 + g
t r r r z r r r r r
2
r z
Vz V V Vz V p 1 (rVz ) 1 2Vz 2V
( + Vr z + + Vz z ) = − + + 2 + 2 + g z (F)
t r r z z r r r r
2
z
Nota: El termino V2 / 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 V2 + V2 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 (G)
r r sen
V V V V V V VrV V2 cot 1 p 2 V
( + Vr + + + − )=− + ( 2V + 2 r
t r r rsen r r r r
V 2cos V
− − ) + g (H)
r 2 sen2 r 2 sen2
V V V V V V VVr V V 1 p
( + Vr + + + + cot ) = − + ( 2V
t r r rsen r r rsen
V 2 Vr 2cos V
− + + 2 2 ) + g (I)
r sen
2 2
r sen r sen
2 2
La ecuación (J) es aplicable a las ecuaciones (G), (H) e (I)
1 2 1 1 2
2 = ( r ) + ( sen ) + (J)
r 2 r r r 2 sen r 2 sen2 2
8.9.2 La ecuación de conservación de energía
65
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Ley de Fourier
q = −k T (8.142)
Coordenadas rectangulares Coordenadas cilíndricas Coordenadas esféricas
dT dT dT
qx = −k (A) qr = −k (D) qr = −k (G)
dx dr dr
dT 1 dT 1 dT
q y = −k (B) q = −k (E) q = −k (H)
dy r d r d
dT dT 1 dT
qz = − k (C) qz = − k (F) q = −k (I)
dz dz rsen d
Ley general de conservación de energía en términos de la densidad de flujo de calor
En primer lugar, se debe tener en cuenta que hay muchas formas de presentar la ecuación de
conservación de energía, dependiendo si se trata de un sólido, o fluido. La forma general útil
desde el punto de vista de la ingeniería química es la que se muestra en la ecuación (8.143) y es
aplicable a sólidos y fluidos incompresibles.
T
C p ( + V .T ) = −.q + G (8.143)
t
El termino de generación puede ser de origen viscoso, nuclear, eléctrico o químico.
Expandiendo:
Coordenadas rectangulares
T T T T q q q
C p ( + Vx + Vy + Vz ) = −( x + y + z ) + G (A)
t x y z x y z
Coordenadas cilíndricas
T T V T T 1 (rqr ) 1 q qz
C p ( + Vr + + Vz ) = −( + + )+G (B)
t r r z r r r z
Coordenadas esféricas
T T V T V T 1 (r 2 qr ) 1 ( sen q ) 1 q
Cp( + Vr + + ) = −( 2 + + ) + G (C)
t r r rsen r r rsen rsen
Nota: Las ecuaciones (8.154) al (8.156), son aplicables para los casos de conductividad constante y conductividad
variable.
Ley general de conservación de energía en términos de la temperatura
Reemplazando (8.152) en (8.153) se tiene la siguiente forma general
T
C p + V .T = k2T + G (8.144)
t
Expandiendo
Coordenadas rectangulares
66
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
T T T T 2T 2T 2T
C p ( + Vx + Vy + Vz ) = k( 2 + 2 + 2 ) + G (A)
t x y z x y z
Coordenadas cilíndricas
T T V T T 1 T 1 2T 2T
C p ( + Vr + + Vz )=k r + 2 2 + 2 ) + G
(B)
t r r z r r r r z
Coordenadas esféricas
T T V T V T 1 r 2T 1 T 1 2T (C)
Cp( + Vr + + )=k 2 + 2 sen + 2 2 +G
t r r rsen r r r r sen r sen 2
Nota: Las ecuaciones (8.158) al (8.160), son aplicables solo a los casos de conductividad constante.
La generación de energía en sus diferentes formas
El termino de generación se pueden expresar de los modos:
Generación de origen eléctrico: Por lo general, la disipación de origen eléctrico se da a baja
temperatura, por tanto, su tasa de generación puede considerarse constante, es decir: G = G0
Generación de origen nuclear: Este tipo de generación puede ser expresado en términos de la
temperatura o en términos de la distancia desde el núcleo hacia la parte externa, es decir:
G = G0 f (T ) ; G = G0 f ( x, y, z ) (Coordenadas rectangulares); G = G0 f (r ) (Coordenadas
cilíndricas y esféricas)
Generación de origen químico: Está relacionada íntimamente con la velocidad de reacción y la
entalpia de reacción, así que se puede usar la siguiente expresión: G = H R rA .
Generación de origen viscoso: Está relacionada directamente con las altas velocidades a las que
se desplaza un fluido, y su tasa de generación según (Bird, et al. 1993; Bird, et al. 2006); esta tasa
de generación viene dada por: G = V ; donde V es la disipación viscosa, la cual viene
expresada por las siguientes ecuaciones dependiendo del sistema de coordenadas.
Coordenadas rectangulares
V 2 V y 2 V 2 V y V 2 V V y 2 V Vz
2
V = 2 x
+ +
z
+ + x
+
z
+ +
x
+ −
x y z x y y z z x
2 Vx V y Vz
2
+ + (8.145)
3 x y z
Coordenadas cilíndricas
V 2 1 V Vr 2 Vz 2 V 2
1 Vr 1 Vz V Vr Vz
2 2
V = 2 r + + + + r + r + r + z + z + r −
r r r z r r
2 1 (rVr ) 1 V Vz
2
+ + (8.146)
3 r r r z
Coordenadas esféricas
67
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
V 2 1 V V 2 1 V V + V cot 2 V 2
1 Vr
V = 2 r +
+ r + + r
+ r
+ +
r r r rsen r
r r r
2 2
sen V 1 V 1 Vr V
+ + +r
r sen rsen rsen r r
1 V
2
2 1 (r 2Vr ) 1
− 2 + (V sen ) + (8.147)
3 r r rsen rsen
Nota: En los casos de problemas de aplicación de generación de origen viscoso, necesariamente, también
se tiene que realizar el balance de cantidad de movimiento.
8.9.3 Transporte de materia en los diferentes sistemas de coordenadas.
En primer lugar, presentamos la ley de Fick aplicable al transporte molecular
Primera forma de la ley de Fick
Densidad de flujo molar de A respecto a la velocidad media molar
J A = − DC A (8.148)
Coordenadas rectangulares Coordenadas cilíndricas Coordenadas esféricas
dCA dCA dCA
J Ax = − D (A) J Ar = − D (D) J Ar = − D (G)
dx dr dr
dC 1 dCA 1 dCA
J Ay = −D A (B) J A = −D J A = − D (H)
dy r d (E) r d
dC dC 1 dCA
J Az = −D A (C) J Az = −D A (F) J A = − D (I)
dz dz rsen d
Segunda forma de la ley de Fick
Densidad de flujo molar de A respecto a un punto de referencia fijo, para un sistema binario o
multicomponente
N A = J A + y A ( N A + N B + ...) (8.149)
Ley general de conservación de materia referido al componente A
C A
+ V .C A = −.N A + rA (8.150)
t
Nota: Estas ecuaciones pueden ser usadas para los siguientes casos:
Sistemas en contradifusion equimolar y no equimolar; sistemas con difusividad constante o
variable; sistemas con cinética compleja
Coordenadas rectangulares
68
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
CA C C C N N N
+ Vx A + Vy A + Vz A = −( Ax + Ay + Az ) + rA (A)
t x y z x y z
Coordenadas cilíndricas
CA C V CA C 1 (rN Ar ) 1 N A N Az
+ Vr A + + Vz A ) = −( + + ) + rA (B)
t r r z r r r z
Coordenadas esféricas
CA C V CA V CA 1 (r 2 N Ar ) 1 ( sen N A ) 1 N A
+ Vr A + + ) = −( 2 + + ) + rA (C)
t r r rsen r r rsen rsen
Si reemplazamos la ecuación de la ley de Fick en su primera forma en la ecuación 8.150 se obtiene
C A
+ V .C A = D 2 C A + rA (8.151)
t
Estas ecuaciones son utilizadas en los casos de contradifusion equimolar y para el caso de
difusividad constante.
Coordenadas rectangulares
CA C C C 2C 2C A 2C A
+ Vx A + Vy A + Vz A = D( 2A + + 2 ) + rA (A)
t x y z x y 2 z
Coordenadas cilíndricas
C A C A V C A C A 1 C A 1 2C A 2C A
+ Vr + + Vz ) = D r + 2 + + rA (B)
t r r z r r r r
2
z 2
Coordenadas esféricas
C A C A V C A V C A 1 C 1 C A 1 2C A
+ Vr + + ) = D 2 r2 A + 2 sen + + rA
t r r rsen r r r r sen r 2 sen 2 2
(C)
8.10 Analogías de los fenómenos de transporte difusional
= −
Transporte de cantidad de Transporte de energía Transporte de materia
movimiento
: Densidad de flujo viscoso : Densidad de flujo de calor o : Densidad de flujo de
de transporte de cantidad de flux de calor (energía calorífica) materia o flux de materia, en
movimiento o flux de cantidad términos de masa o moles.
de movimiento.
: Propiedad física (viscosidad) : Propiedad física : Propiedad física (difusividad
(conductividad térmica) de materia)
: Potencial de transferencia : Potencial de transferencia : Potencial de transferencia
(gradiente de velocidad) (gradiente de temperatura) (gradiente de concentración,
fracción molar o presión)
69
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
UNIVERSIDAD NACIONAL DEL CALLAO – FACULTAD DE INGENIERIA QUÍMICA
Cantidad de movimiento Energia Moles
= = =
sup erficie tiempo sup erficie tiempo sup erficie tiempo
mV Kg (m / s) N E J W n moles
= : = 2 = : 2 = 2 = : 2
At m2 s m At m s m At m s
= −V Ley de Viscosidad q = −kT J A = − DC A
de Newton Ley de Fourier Ley de Fick
Trabajos de laboratorio de ingeniería III
[Link]
[Link]
70
Fenómenos de Transporte Cap. VIII - Luis Carrasco y col.
También podría gustarte
- Álgebra Vectorial en Fenómenos de TransporteAún no hay calificacionesÁlgebra Vectorial en Fenómenos de Transporte56 páginas
- Introducción a Magnitudes FísicasAún no hay calificacionesIntroducción a Magnitudes Físicas26 páginas
- Introducción A Los Fenómenos de Transporte para Ingenieros Químicos88% (8)Introducción A Los Fenómenos de Transporte para Ingenieros Químicos126 páginas
- Vectores en Física: Definiciones y OperacionesAún no hay calificacionesVectores en Física: Definiciones y Operaciones14 páginas
- Soluciones Ideales. Equilibrio Líquido VaporAún no hay calificacionesSoluciones Ideales. Equilibrio Líquido Vapor27 páginas
- 3.4 Apendice A. Notacion Vectorial y TensorialAún no hay calificaciones3.4 Apendice A. Notacion Vectorial y Tensorial16 páginas
- Introducción al Análisis Vectorial en FísicaAún no hay calificacionesIntroducción al Análisis Vectorial en Física20 páginas
- Análisis Vectorial en Física GeneralAún no hay calificacionesAnálisis Vectorial en Física General50 páginas
- CM1 - Magnitudes Escalares y VectorialesAún no hay calificacionesCM1 - Magnitudes Escalares y Vectoriales27 páginas
- CEM4042 Cap 1-2 - Análisis Vectorial Sadiku-2014Aún no hay calificacionesCEM4042 Cap 1-2 - Análisis Vectorial Sadiku-201419 páginas
- Álgebra Lineal: Vectores y OperacionesAún no hay calificacionesÁlgebra Lineal: Vectores y Operaciones12 páginas
- Introducción a Vectores y MagnitudesAún no hay calificacionesIntroducción a Vectores y Magnitudes12 páginas
- Vectores en Física: Definiciones y OperacionesAún no hay calificacionesVectores en Física: Definiciones y Operaciones12 páginas
- 0 - Intro - Magnitudes Vectoriales y EscalaresAún no hay calificaciones0 - Intro - Magnitudes Vectoriales y Escalares14 páginas
- Fundamentos de Estatica y DinamicaAún no hay calificacionesFundamentos de Estatica y Dinamica110 páginas
- Álgebra Vectorial: Conceptos y AplicacionesAún no hay calificacionesÁlgebra Vectorial: Conceptos y Aplicaciones40 páginas
- Introducción a Vectores en FísicaAún no hay calificacionesIntroducción a Vectores en Física5 páginas
- Introducción al Cálculo Vectorial y TensorialAún no hay calificacionesIntroducción al Cálculo Vectorial y Tensorial35 páginas
- Teoría Electromagnética y CálculoAún no hay calificacionesTeoría Electromagnética y Cálculo275 páginas
- Semana 02 - Vectores-Teoria y PracticaAún no hay calificacionesSemana 02 - Vectores-Teoria y Practica71 páginas
- Escalares y Vectorales. Investigacion.Aún no hay calificacionesEscalares y Vectorales. Investigacion.13 páginas
- Manual de Laboratorio de PolímerosAún no hay calificacionesManual de Laboratorio de Polímeros23 páginas
- 1.5.3 Balance de Cantidad de Movimiento en Régimen Estacionario en Un Sistema de Coordenadas CilíndricasAún no hay calificaciones1.5.3 Balance de Cantidad de Movimiento en Régimen Estacionario en Un Sistema de Coordenadas Cilíndricas8 páginas
- R R H H T T: Transporte de Calor en Una Pared Cilíndrica Con Transferencia de Calor Por Convección en Ambos LadosAún no hay calificacionesR R H H T T: Transporte de Calor en Una Pared Cilíndrica Con Transferencia de Calor Por Convección en Ambos Lados50 páginas
- 3.8.-Transporte de Masa de Fluidos en Solidos PorososAún no hay calificaciones3.8.-Transporte de Masa de Fluidos en Solidos Porosos46 páginas
- Modelo UNIFAC: Datos de Acetona/AguaAún no hay calificacionesModelo UNIFAC: Datos de Acetona/Agua5 páginas
- Estrategias de Discipulado en Tito 2Aún no hay calificacionesEstrategias de Discipulado en Tito 25 páginas
- Aseguramiento de Validez de ResultadosAún no hay calificacionesAseguramiento de Validez de Resultados12 páginas
- ¿Cómo Puedes Formar Una Buena Relación o Reparar Una Que Se DañóAún no hay calificaciones¿Cómo Puedes Formar Una Buena Relación o Reparar Una Que Se Dañó11 páginas
- Metodo Runge-Kutta 2do Orden Aldo 2MM4Aún no hay calificacionesMetodo Runge-Kutta 2do Orden Aldo 2MM42 páginas
- Taller Sobre Validación de Documentos JMSD SenaAún no hay calificacionesTaller Sobre Validación de Documentos JMSD Sena9 páginas
- Válvulas en Instalaciones de AguaAún no hay calificacionesVálvulas en Instalaciones de Agua23 páginas
- PLANIFICACIÓN DE CLASE - Español 1 - UNG - Fubel y Rondinone (06 - 11)Aún no hay calificacionesPLANIFICACIÓN DE CLASE - Español 1 - UNG - Fubel y Rondinone (06 - 11)8 páginas
- Viaje A Parras de La Fuente Desde GDL - Buscar Con GoogleAún no hay calificacionesViaje A Parras de La Fuente Desde GDL - Buscar Con Google1 página
- Práctica PROPIEDADES DEL AGUA PDFAún no hay calificacionesPráctica PROPIEDADES DEL AGUA PDF3 páginas
- Blue & White Company Profile PresentationAún no hay calificacionesBlue & White Company Profile Presentation8 páginas
- Licenciatura en Educacion Mencion Orientacion AcademicaAún no hay calificacionesLicenciatura en Educacion Mencion Orientacion Academica3 páginas
- I-P-01 V3 Instructivo AlistamientoAún no hay calificacionesI-P-01 V3 Instructivo Alistamiento3 páginas
