II-2019 SEGUNDO TRABAJO DE TRANSFERENCIA DE CALOR

I-Problema de flujo bidimensional

Problema-C: Una lámina bidimensional de cerámica transfiere calor en flujo estacionario.

Para el problema que le corresponda tomar una separación nodal de un centímetro tanto en
dirección-x como en dirección-y. responda lo siguiente:

 Deduzca la ecuación de cada tipo de nodo.

 Aplique las ecuaciones determinadas a la red de nodos mostrada en la figura según el
tipo de nodo.
 Determine la distribución de temperatura.
 Utilice las temperaturas halladas para trazar las líneas isotermas y las líneas de flujo, y
calcule el calor total transferido.

Paula Acuña, César Bracamonte, Daniela Martínez, Laura Villalobos

Nota: Resuelva el problema que le corresponda. De lo contrario tendrá nota nula

II- Problema: Optimización Intercambiador de coraza y tubos (Diseño térmico)

Para el problema 11-1 (Kern) determine el diámetro de la coraza y el espaciamiento de

deflectores que hace óptimo el coste de operación del diseño térmico del intercambiador.
Para ello plantear función objetivo respecto de costes de operación, sin considerar gastos
fijos. El coste de operación es
  Coste   Coste 
 Coste     Coste   
 Coste     de bombeo    de bombeo
    de fluido      de fluido   fluido 
 de operacion    fluido
 caliente    frio  
 
 caliente   frio 

Los costes unitarios se definen como

Ccaliente : Coste unitario del fluido caliente $/Lb

Cb.caliente : Coste unitario de bombeo fluido caliente= Coste de energía al freno, suministrada
a la bomba en kilovatio-hora $/hora.

C frio : Coste unitario del fluido frío $/Lb

Cb. frio : Coste unitario de bombeo fluido frío= coste de energía al freno, suministrada a la
bomba en kilovatio-hora $/hora.

La potencia al freno utilizando Bernoulli modificado se calcula como

 Po Wg
Wpc  Para el fluido caliente en la coraza
 o pc

 pc : Eficiencia de la bomba

La relación empírica para calcular la caída de presión en el lado de la coraza es

f o Go Do  N  1
 Po  Lb/plg2
5.22  1010 Deo so

N : Número de bafles o deflectores, N  1 : Número de cruces verticales que el fluido hace en el lado de la
coraza.

Para el fluido frío en los tubos la potencia al freno utilizando Bernoulli modificado se
calcula como

 Pi wg
Wpf 
 i pf
pf : Eficiencia de la bomba.

La caída de presión en el lado de los tubos es

 PT   PL   Pr en Lb/plg2

Las relaciones empíricas para calcular  PL y  Pr son

 f i Gi 2 L n 4n Vi 2 62.5
 PL  y  Pr  en Lb/plg2
5.22  1010 d i si i s 2 g  144

n : número de pasos en los tubos.

Relaciones de diseño: Se trabaja con las relaciones de diseño que figura en el algoritmo de
cálculo del Kern (página 275-277).

Para reducir significativamente el número de relaciones de diseño que allí figuran y

simplificar el algoritmo de cálculo con el propósito de optimizar, se realiza lo siguiente:

1- Trabajar en lugar de las temperaturas calóricas, trabajar con temperatura media en cada
fluido.

2- La ecuación (6.1 Kern) utilizada para calcular el coeficiente de T.Calor en el lado de los
tubos, la cual está en función de “hi=f(k,di,Gi,mi,miw)”. En esta ecuación se ha sustituido
Gi por w/ai y también posteriormente se sustituyo “ai” por Nt*((pi*di^2)/4)/n. para que el
coeficiente quede finalmente

2 1
2
k i 3 w 0.8 n 0.8 c 3
hi  8185
.  10 (Se recomienda hacer el desarrollo mencionado en el
 0.8 d i 1.8 N t 0.8 i 0.327 iw 0.14
numeral 2 e incluirlo en el trabajo, para verificar si esta ecuación es correcta)

N t : Número de tubos, n : Número de pasos en los tubos.

3- La ecuación (que figura en la parte superior de la página 171 Kern) utilizada para
calcular el coeficiente de T.Calor en el lado de la coraza, la cual está en función de “ho=f(k,
De ,Gs, mo, mow)”. En esta ecuación se ha sustituido Gs por “Gs=W/as” y también
posteriormente se sustituyo “as” por as=Ds*B*cc/PT, (Donde cc es el claro o ligadura) para
que el coeficiente ho quede finalmente

2 1
k o 3W 0.55 PT 0.55 C 3
ho  0.36 0.45 0.55 0.55 0.55 0.077 0.14 (Se recomienda hacer el desarrollo mencionado en el
De Ds B cc o ow
numeral 3 e incluirlo en el trabajo, para verificar si esta ecuación es correcta)

PT : Pitch, cc : Claro o ligadura

Según lo mencionado anteriormente el número de relaciones de diseño, finalmente se

redujo a 14, por lo cual para este problema de optimización trabajar con las siguientes
relaciones de diseño:

Q=w*c*(t2-t1) (1)
Q=W*C*(T1-T2)
DMLT=((T1-t2)- (T2-t1))/log((T1-t2)/ (T2-t1)) (2)
FT=FT(T1,T2,t1,t2) (3)
A=Q/Ud/DMLT/Ft (4)
% De la tabla 10 para el diámetro nominal y el calibre que da el problema leer “di” y
espesor de pared.
% A=pi*do*L
L=A/pi/do (5)
% Número de tubos
Nt=L/16 (6)

% Leer de la tabla 9 con el número de tubos que más se aproxime al anterior el diámetro de
la coraza

Ds=f(tabla 9) (7)

% Cálculo de coeficiente hi

hi = 8.185e-2*k^(2/3)/di^1.8*w^.8/Nt^.8/pi^.8*n^.8/mi^.327*c^(1/3)/miw^.14 (8)

hio=hi*di/do (9)

% Cálculo de coeficiente ho

De=4*(PT^2-(pi/4)*do^2)/pi/do (10)

ho = .36*k^(2/3)/De^.45*W^.55/Ds^.55/B^.55/cc^.55*PT^.55/mo^7.7e-
2*c^(1/3)/mow^.14 (11)

Uc=hio*ho/(hio+ho) (12)
Rd =1/Ud-1/Uc (13)
PT=cc+do (14)

% Han resultado 14 relaciones de diseño considerando la relación de viscosidades como

“1”

Variables : Q, w, W, t1, t2, T1, T2, DMLT, A, Ud, Ft, L, Nt, hi, hio, Ds, PT, ho, hio, De, B,
cc, Uc, Rd

24 variables, 8 variables especificadas, 14 Ecuaciones

Numero de variables de diseño=24-8-14=2

- Plantear el problema de optimización: Función objetivo más restricciones. Utilice como

restricciones, la restricción de la caída de presión en cada corriente:  Po  10 psi y
 Pi  10 psi, además considere U 
D inf erior  U D  U D  sup erior (este intervalo lo da la
tabla-8 apéndice Kern).
- Trabajar como variables de diseño el diámetro de la coraza Do , y B (espaciado de
deflectores). El espaciado de los deflectores según recomienda el Kern varía en el intervalo

1
D  B  Do
5 o

Para optimizar es necesario expresar la función objetivo en términos de las variables de

diseño. Para lo cual se realiza lo siguiente

Función objetivo en términos de las variables de diseño


C0 p  Ccaliente w  C frioW  Cb.calienteWpc  Cb. frioWpf  
 Po Wg  Pi wg
Wpc  y Wpf 
 o pc  i pf

f o Go Do  N  1
En la ecuación  Po  se sustituye Go por “Go=W/(Do*B*cc/PT)” y Deo
5.22  1010 Deo so
por “Deo=Do-do” y N+1 por “N+1=L/B” resultando

f oWPT L
 Po  1916  10 11
B cc  D0  d 0 so
. 2
(Efectúe las sustituciones para verificar este resultado)

Nota: Estimen un factor de fricción f 0 y trabajen con el factor de fricción constante. O si quieren busquen en
la literatura una ecuación para calcular el factor de fricción.

De los datos de la tabla de la figura-9 (apéndice Kern) para ¾ de pulgada diámetro externo
de tubo, hacer una regresión de los datos Ds (diámetro de la coraza) VS número de tubos
(descartar el primer dato de la tabla) para la columna del intercambiador 1-2. En la
regresión utilizar la ecuación “Do=q*Nt^m” hallar el parámetro “q” y el exponente “m”.
Posteriormente utilizar la relación de diseño (6) y la ecuación de regresión “Do=q*Nt^m”
para obtener L en función de Do “ verificar si al hacer lo anterior da L=16* (Ds/q)^(1/m)”
y obtenido “L ” sustituirlo en la caída de presión. Con ello  Po queda en función de Do y B
(variables de diseño).

Caída de presión en el lado de los tubos en términos de las variables de diseño:

f i Gi 2 L n 4n Vi 2 62.5
Dado  PL  y  Pr 
5.22  1010 d i si i s 2 g  144

Nota: Estimen un factor de fricción f i y trabajen con el factor de fricción constante. O si quieren busquen en
la literatura una ecuación para calcular el factor de fricción.
- En  PL Sustituya “L” obtenido del procedimiento anterior
- En  Pr exprese Vi en términos de las variables de diseño: Para lo cual utilice las
ecuaciones:
w  Vi ai ai=Nt*((pi*di^2)/4)/n Do=q*Nt^m”

De la anterior ecuación explicite “Nt” y sustitúyalo en “ai” y luego sustituya “ai” en

w  Vi ai y de la ecuación resultante explicite Vi , la cual queda en términos de una o
ambas variables de diseño. Y finalmente sustituya Vi en  Pr , quedando en términos de
variables de diseño.

Con ello finalmente la función objetivo queda en términos de las variables de diseño

Nota: inicialmente considerar las variables de diseño como continúas y continuar el

procedimiento de optimización:

- Construya la región de restricciones en el plano Do Vs B

- Determine por derivación de la función objetivo el punto crítico de la función objetivo y

determine si el punto crítico cae en la región de restricciones. Si es así, utilice ahora el
resultado de Do óptimo y lea de la tabla 9 el diámetro de la coraza comercial que más se
acerque al valor de Do calculado y lea el número de tubos y utilice las relaciones de diseño
anteriores para determinar el diseño térmico final.

- Si el punto crítico de la función objetivo no cae en la región de restricciones, Busque el

óptimo calculando el coste de operación, para ello hay que desplazarse en las fronteras de
la región de restricciones (utilice Do y B leído de las fronteras para calcular el coste).
Obtenido el óptimo utilice ahora el resultado de Do óptimo y lea de la tabla 9 el diámetro
de la coraza comercial que más se acerque al valor de Do calculado y lea el número de
tubos y utilice las relaciones de diseño anteriores para determinar el diseño térmico óptimo
final.

Con respecto al problema de optimización responda lo siguiente:

1- ¿Cuál es el desarrollo para llegar a la ecuaciones

2 1
2
k i 3 w 0.8 n 0.8 c 3
hi  8185
.  10
 0.8 d i 1.8 N t 0.8 i 0.327 iw 0.14 y
2 0.55 0.55 1
k o 3W PT C 3
ho  0.36
De 0.45 Ds 0.55 B 0.55 cc 0.55 o 0.077 ow 0.14

2- Escriba la función objetivo en términos de las variables de diseño

 Cop=f(Diámetro de la coraza, espaciamiento entre bafles)

3- Escriba la región de restricciones.

4- Represente la región de restricciones, en el espacio coordenado de variables de diseño.

5- Dibuje en el diagrama 3D: Cop Vs Diámetro de la coraza Vs B. y determine las

condiciones optimas del diseño térmico para esta operación de transferencia de calor