Numero de Brinkman

Universidad de Guanajuato
División de Ciencias Naturales y Exactas
Ingeniería Química
Transferencia de Calor
Mtro.: Jesús Isaac Minchaca Mojica
Tarea
Alumno: Juan Carlos Yebra Pérez
Fecha: 11 de octubre de 2016
Caso de estudio: “Análisis por colocación ortogonal de los

efectos de conductividad axial y disipación viscosa en la
transferencia de calor de fluidos no newtonianos”
Planteamiento
El caso en el que se enfoca en el desarrollo de un modelo teórico de transferencia
de calor análisis de fluidos no-Newtonianos, en el cual se hace inclusión de varias
contribuciones tales como la convección axial, conducción axial, conducción radial
y disipación viscosa, dando que el modelo hace ilusión a un cilindro, por lo cual las
ecuaciones que se van a ir presentando van a tener en cuenta el radio y el valor
de la altura.
El fluido en el cual se está llevando el estudio es un sistema polímero-agua, en el
cual el polímero a usar fue:
 wsr n-3000, a concentraciones de 2000 y 4000 ppm
Además se usaron para el estudio:
 Bentonita, a composiciones 4,6 y 8% P/V

 Mezclas de Crudo-Gas Oil en relaciones volumétricas del 1:1, 2:1, 3:1 y 3:2
 Aceite-Diesel Lubricante.
Consideraciones
Las consideraciones en las cuales están trabajando son en las cuales se está un
tubo de radio R y de longitud L:
 Que el perfil de velocidades que está completamente desarrollado y en
régimen laminar.
 Que la generación de calor interna no existe.
 Las propiedades de densidad, capacidad calorífica y conductividad térmica
se consideran valores promedios.
Desarrollo
A partir de las consideraciones que se plantearon se llegó a la ecuación diferencial
de temperatura
∂T ∂2T k ∂ r ∂ T 1 ∂V z
ρC V z
∂z
=k 2 +
∂z r ∂ z ∂r( )+
778
τ rz
∂r
Donde T (r, z) es la temperatura del fluido a cualquier posición (r, z), Vz es el

componente axial de la velocidad del fluido, τ rz es el componente tangencial del
esfuerzo viscoso.
Se sustituye el término de disipación viscosa en función de la temperatura media
por la aproximación de Cochrane, y combinando la definición de viscosidad
aparente en el modelo de potencia:
τ rz gc
μ=
∂Vz
∂r( ) , τ rz gc=m
∂Vz
∂r ( )
η
A
Aproximado μ=μ ∞ e T m resulta la expresión para la velocidad de deformación:
1 A
∂V z μ ∞
∂r
=
m ( ) n−1
e T m (n−1)
Y el componente de disipación viscosa:

n+1
∂V z m ∂Vz
τ rz
∂r
=
gc ( )( ) ∂r
n+1 A (n +1)
m μ∞
¿ ( )( )
gc m
n−1
e T m(n−1)
Donde Tm representa la temperatura medida del fluido, dada por:

R
∫ V z τ rz ( r , z ) rdr
0
T m ( z )= R
∫ V z rdr
0
El coeficiente convectivo de transferencia de calor, al combinar la ley de

enfriamiento de Newton con un balance diferencial de calor:
d T m −T ( r , 0 )
h=
(
ρC V z D dz T w −T ( r , 0 ) )
4 T m−T ( r , 0 )
−1
T w −T ( r , 0 )
Las condiciones de transferencia:

1) Perfil de velocidad aproximado a un modelo de potencia
Vz r 1+
1
V z (0)
=1−
z () n
2) Temperatura de entrada constante, T(r,0); Temperatura del fluido a la pared

constante en la región de transferencia T(R,Z); flujo axial de calor al
∂T
extremo de la tubería, q z =−k =h(T m −T w )
∂z
Vz z r T ( r , z )−T ( r , 0 )
Adimensionando con u=
V z (0)
, ź= ,
L
x=
R
, θ= (
T ( R , z )−T ( r , 0 ) ) e
introduciendo los números adimensionales

V ( 0) L
 Peclet ( Pe )= z
α
μ∞
 Prandtl ( Pr❑ )=
ρα
V z ( 0 ) Dρ
 Reynolds ( ℜ )=
μ∞
1 mV z ( 0 )
 Brinkman ( Br )=
778 gk Ln+ 1 ΔT
L

Parametro Ω=
V z ( 0) (
n+1 μ
m
∞ n +1
n−1
( ) )
Se puede escribir la ecuación diferencial de temperatura como:
2 N +1 A ( n+1 )
dθ 1 d θj 2L 1 BrΩ
U j j=
dz Pe d z 2
+
R ( )( RePr ) ∑ B ji θ ij + Pe e (n −1 )Tj ; j=1 … N
j=1
Donde
N +1
W i U i θij
T j−T ( r ,0 ) ∑ j=1
= N+1
T w −T ( r , 0 )
∑ W iUi
j=1
M +1 N +1
RePr
∑ W j U j ∑ Aui θ ji
j=1 i=1
Nu , k= M +1 M +1
4L
D ∑ W j U j θ jk − ∑ W j U j
j=1 j=1
1
1+
De igual manera, las condicione de transferencia: u =1−x n , donde xj es el punto
j j
de colocación ortogonal
θ M+1 ,i =1 a x M+1=1
θ j ,0=0 a z=0
M +1
N+1
∑ A N +1 , L θ j , L=Nu
L=1
( ∑ W j U j θ j , N +1
j=1
M +1
∑ WjU j
j=1
)
−1 ; j=1 … . M
Mediante la utilización de un programa de computo se logró resolver el sistema de

ecuaciones diferenciales obtenido al aplicar el método de colocación ortogonal
(asumiendo una geometría plana asimétrica y geometría cilíndrica simétrica) a la
ecuación diferencial en derivadas parciales provenientes del balance diferencial de
temperatura para un fluido viscoso a través de una tubería re radio R y longitud L.
El programa calcula el perfil de temperatura tanto en la dirección radial como axial

en los diferentes puntos de colocación ortogonal, así como también la variación
del número de Nusselt en la posición axial.
Resultados
El análisis es realizado tomando N= 5 y N= 6 puntos de colocación ortogonal, para
geometría plan r asimétrica y geometría cilíndrica simétrica respectivamente;
considerando un factor pesante W= l, lo que permite un mejor enfoque del
proceso, ya que se estudian siete (7) posiciones radiales en la tubería; implicando
esto una mayor exactitud de los cálculos.
Se estudiaron soluciones que fueron previamente descritas, para lo cual se
promediaron valores de los números de Reynolds, Prandtl y Brinkman.
La disminución del no. de Prandtl hace que el valor del no. De Nusselt sea menor
para una misma posición axial. Como consecuencia el perfil de temperatura tanto
en la dirección axial como radial es creciente, esto determina una disminución en
la conductividad térmica, un incremento en el de no. de Prandtl y a su vez en el
no. de Nusselt.
Los efectos asociadles de los no. de Reynolds y Prandlt sobre el número de

Nusselt se correlacionan así:
 Para la mezcla Polímero-Agua:
Nu= 1.803 Re= 0.2693 Pr= 0.2789;
(L/D) ≤ 30 y Br ≤ 0.00348
 Para la mezcla Bentonita-Agua:

Nu= 0.492 Re= 0.3402 Pr= 0.3408;
(L/D) ≤ 30 y Br ≤ 0.000228
Se muestran los efectos del no. de Brinkman sobre el no. de Nusselt a través de la
longitud de transferencia. Así mismo también, a medida que se aumenta la
posición axial, el valor del no. de Nusselt se hace cada vez menor. El incremento
del no. de Brinkman implica una mayor conversón de los efectos friccionales del
fluido en energía térmica y en consecuencia un mayor gradiente de temperatura a
la pared del fluido, y disminuye en la medida que resulte completamente
desarrollado el perfil de temperatura o resulten menos apreciables los efectos
viscosos frente a la contribución convectiva de calor.
Conclusiones
La utilización del método de colocación ortogonal en combinación con un método
numérico, permite resolver y analizar en forma rápida el modelo diferencial que
rige el proceso de transferencia de calor.
El comportamiento de fluidos no-newtonianos, de naturaleza orgánica e
inorgánica, con propiedades dependientes de temperatura fueron analizadas
incluyendo convección axial, conducción radial y axial y disipación viscosa para la
componente tangencial del esfuerzo viscoso.
Bibliografía
 Corone, J., Zerpa, R., Arrieta, I., & García, C. (1987). Análisis por
colocación ortogonal de los efectos de conductividad axial y disipación
viscosa en la transferencia de calor de fluidos no-newtonianos. Revista
Técnica de la Facultad de Ingeniería, 10, pp.119-127.

Numero de Brinkman

Cargado por

Copyright:

Formatos disponibles

Numero de Brinkman

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Numero de Brinkman

Cargado por

Copyright:

Formatos disponibles

Universidad de Guanajuato

División de Ciencias Naturales y Exactas

Mtro.: Jesús Isaac Minchaca Mojica

Alumno: Juan Carlos Yebra Pérez

Fecha: 11 de octubre de 2016

Caso de estudio: “Análisis por colocación ortogonal de los

 Bentonita, a composiciones 4,6 y 8% P/V

Donde T (r, z) es la temperatura del fluido a cualquier posición (r, z), Vz es el

Y el componente de disipación viscosa:

Donde Tm representa la temperatura medida del fluido, dada por:

El coeficiente convectivo de transferencia de calor, al combinar la ley de

Las condiciones de transferencia:

2) Temperatura de entrada constante, T(r,0); Temperatura del fluido a la pared

introduciendo los números adimensionales

Mediante la utilización de un programa de computo se logró resolver el sistema de

El programa calcula el perfil de temperatura tanto en la dirección radial como axial

Los efectos asociadles de los no. de Reynolds y Prandlt sobre el número de

 Para la mezcla Bentonita-Agua:

También podría gustarte