Diseño Pasteurizador

CÁLCULOS JUSTIFICATIVOS 10.
Diseño del pasteurizador
10.2. Cálculo y diseño del intercambiador de calor de

placas
10.2.1. Método de cálculo
El método de cálculo que se ha realizado para cada una de las secciones del
intercambiador de calor de placas es el siguiente:
A. Balance térmico. Cálculo del caudal másico de los fluidos y calor total
intercambiado.
B. Diferencia de temperatura media logarítmica.
C. Cálculo de los coeficientes individuales de transmisión de calor de ambos
fluidos.
D. Cálculo del coeficiente global de transmisión de calor sucio; se denomina
así porque se incluyen las resistencias térmicas originadas por los
contaminantes de los fluidos caliente y frío que se depositan en las dos
superficies de las placas.
E. Cálculo del área de transferencia
F. Cálculo del número de placas por sección.
G. Cálculo de las pérdidas de carga originadas por cada fluido en cada
sección.
H. Cálculo de la longitud del tubo de mantenimiento.
10.2.2. Hipótesis de cálculo
Para el diseño del intercambiador de calor para la pasteurización de Qmix= 5000

l/h, se considera flujo en contracorriente, ya que se consigue un área de
transferencia de calor mayor que en el caso de flujo en paralelo.
Las hipótesis que se establecen para el estudio serán:
1. El intercambiador de calor opera en condiciones estacionarias.

2. El coeficiente global de transferencia de calor es constante a lo largo de
todo el intercambiador.
3. El equipo de pasteurización deberá ser lo más versátil posible y operar con
todas las composiciones de mixes posibles, desde los que el porcentaje en
sólidos totales y materia grasa son elevados hasta los que presentan
valores pequeños, como es el caso de las mezclas de sorbetes. Se partirá de
una mezcla cuya composición (40 % en sólidos totales y 20 % en materia
grasa) representa las condiciones de viscosidad más altas para las que se
realizará el diseño del intercambiador de calor, es decir, se han
seleccionado para los cálculos las condiciones más desfavorables.
4. El calor específico del mix se considerará que se mantiene constante con la
temperatura, al igual que la densidad; al ser el mix un producto de
composición variable y suma de muchos otros ingredientes, resulta difícil
encontrar tabulados la variación de dichas propiedades con la temperatura
107
CÁLCULOS JUSTIFICATIVOS 10. Diseño del pasteurizador
como ocurre con otros productos de composición constante como la leche,

la mantequilla, la nata, etc.
5. Las temperaturas y velocidades de los fluidos son uniformes en toda la
superficie de la placa. Esto no es rigurosamente cierto, ya que las
propiedades de los fluidos dependen de la temperatura, pero para
simplificar el cálculo, y puesto que no se conoce la temperatura de la
mezcla en cada punto se considerará que cada que cada fluido se encuentra
a la temperatura media entre la entrada y la salida.
6. No existe conducción de calor en la dirección del flujo.
7. El caudal de fluido está repartido equitativamente entre las distintas
superficies de las placas.
8. Las pérdidas de calor al exterior son despreciables
9. Se dispone de un intercambiador con cinco secciones a cada una de las
cuales se aplicarán los balances oportunos.
108
Cálculo de la Sección III (Sección de calentamiento)
La ecuación de energía, suponiendo que no existen pérdidas térmicas es la

siguiente:
ΔH + ΔEc + ΔEp = Q + We (Balance energético)
Para que la expresión anterior queda:
ΔH = Q
según las siguientes simplificaciones:
Formas mecánicas de energía son nulas ΔEc ≅ ΔEp ≅ We = 0
Para cualquiera que sea el tipo de equipo intercambiador empleado, si sólo se

consideran las condiciones de entrada y salida de los fluidos, se puede establecer
el balance térmico global del equipo escribiendo que la cantidad de calor Q
perdida por el fluido caliente es igual a la que gana el fluido frío, si se desprecian
las pérdidas térmicas:
Q= mf (h1 – h2)= mc (H2 – H1)
y por tanto:
mfCpf (Tfe – Tfs) = mcCpc (Tcs -Tce)
En esta sección las propiedades del fluido y según la tabla anterior, la relación de
temperaturas es la siguiente:
Propiedades del fluido Unidades Caliente Frío

Mix sin
Agua pasteurizar
Densidad kg/m3 968,3 1100
Calor Específico kJ/kg ºC 4,201 3,34
Conductividad Térmica W/m ºC 0,673 0,504
Viscosidad Entrada mPa s 0,32 58
Viscosidad Salida mPa s 0,35 42
Temperatura Entrada ºC 90 67
Temperatura Salida ºC 81 84
Tfe = Temperatura del mix sin pasteurizar a la entrada: 67 ºC

Tfs = Temperatura del mix pasteurizado a la salida: 84 ºC
Tce= Temperatura del agua caliente a la entrada: 90 ºC
Tcs =Temperatura del agua caliente a la salida: 81 ºC
109
mf: Caudal másico del fluido frío(mix a pasteurizar) Kg/h como la

ρmix= 1100 kg/m3, para una qent= 5000 l/h, tenemos que mf= 5500 kg/h
Cpf: Calor específico del fluido frío (mix a pasteurizar) = 3,34 kJ/kg ºC
mc: Caudal másico del fluido caliente (agua caliente): desconocido
Cpc: Calor específico del fluido caliente (agua caliente) = Cpc= 4,201 kJ/kg ºC
de esta ecuación sólo se desconoce la variable mc:
donde mc=[ mfCpf (Tfe – Tfs)]/ Cpc (Tcs – Tce)]
mc= [5500 kg/h ×3,34 kJ/kg ºC×(67 ºC - 84 ºC)]/ [4,201 kJ/kg ºC×(81 ºC- 90 ºC)]
=
mc= 8259,673 kg/h = 2,294 kg/s como ρc= 968,3 kg/m3
qc= 8,53 m3/h
Aplicando a continuación la ecuación de Fourier a todo el equipo:
Q= UA Δtm
De la ecuación tanto Q como ΔTm son valores fácilmente calculables.
Cálculo del calor intercambiado (Q) se define como:
Q = mf Cpf (Tfs – Tfe)= 5500 kg/h ×3,34 kJ/kg ºC× (84 ºC -67 ºC)= 312290 kJ/h=
Q =86,747 kW
Diferencia media de temperatura (ΔTm):

El estudio de ΔTm se efectuará para los fluidos en contracorriente y suponiendo
que el coeficiente global de transmisión U es constante, así como los calores
específicos de los fluidos y además no existe cambio de fase por parte de ninguno
de los fluidos.
Por la teoría general de intercambiadores la diferencia de temperaturas entre
ambos fluidos a lo largo de un intercambiador no es constante. Para el cálculo de
la diferencia efectiva se hará uso de la temperatura media logarítmica, que viene
expresada en función de las temperaturas en las secciones de entrada y salida de
los fluidos.
ΔTm = ΔT1 - ΔT2 / (ln (ΔT1 / ΔT2))
Los intercambiadores de placas se diseñan para funcionar con el flujo de los

fluidos en contracorriente, de modo que la expresión utilizada es:
ΔTm = [(Tce – Tfs) – (Tcs – Tfe) / ln ((Tce - Tfs) / (Tcs – Tfe))]
110
ΔTm = [(90 – 84) - (81 – 67) / ln ((90 – 84) / (81 – 67))] = 9,44 ºC
Los canales de los extremos del intercambiador presentan la mitad de superficie

de transferencia de los canales interiores y por tanto transfieren menos calor. Esto
provoca que la efectividad del intercambiador sea algo menor que la que se
predice mediante el uso de la temperatura media logarítmica media. Es necesario
introducir un factor F, menor que la unidad, quedando que la transferencia de
calor a lo largo del intercambiador es:
Q= UAΔTm= UAFΔTlm
Para el cálculo del factor de corrección F es necesario definir la llamada relación

de capacidades (R) y la efectividad (P) que dan información de cómo son las
capacidades caloríficas de ambos fluidos y cual es la efectividad del
intercambiador de calor. Existen tablas en las que entrando con los valores de R y
P se puede calcular el valor de F.
ΔTm= F ΔTlm
ΔTm = (T1 – t2)- (T2 – t1)/ln [(T1 – t2)- (T2 – t1)]= 9,44 ºC
F: Factor de corrección (o de aproximación)
R= (m Cp)f / (m Cp)c = (Tce – Tcs) / (Tfs – Tfe)= (90 – 81) / (84 – 67) =
0,529
P= (Tfs – Tfe)/ (Tce – Tfe) = (84 – 67) / (90 – 67) = 0,739
ΔTm= F ΔTlm
El efecto de la disminución de la efectividad, debido a que la superficie de

transferencia de los canales exteriores es menor que la de los interiores, es
importante cuando se trata con intercambiadores cuyo número de placas es menor
que 50. En estos intercambiadores es necesario el cálculo del factor F para poder
utilizar la temperatura logarítmica media como modelización de la diferencia
efectiva de temperaturas. Para aquellos cuyo número de placas es mayor a 50, los
efectos de los extremos dejan de ser significativos y se pueden despreciar. En este
caso no es necesario el cálculo del factor F y se puede tomar un valor para él de F
= 1. De este modo quedaría que:
ΔTm= ΔTlm
Sin embargo si el número de placas es menor a 50, el factor F se calcula

gráficamente, en la que se representa el factor de aproximación F, en función de
los parámetros adimensionales R y P, para distintas configuraciones geométricas
del cambiador.
Se supondrá que el número de placas en esta sección es inferior a 50.
111
A partir de la gráfica para el cálculo del factor de corrección para

intercambiadores de calor tubulares de un paso lado de la carcasa y dos pasos (o
más) del lado de los tubos y con valores de R=0,529 y P=0,739, se obtiene que
F≈0,57, así:
ΔTm= F ΔTlm = 0,57×9,44 ºC= 5,38 ºC
Obtención del coeficiente global de transmisión de calor (Us sucio)
Cálculo del coeficiente de calor sucio (Us)
El coeficiente global de transferencia de calor viene dado por la expresión:
1/Us = 1/ hc + e/km + 1/ hf + Rf + Rc
o lo que es lo mismo:
Us= 1/ [1/ hc + e/km + 1/ hf + Rf + Rc]
hc = coeficiente de transferencia de calor por convección del fluido caliente.

hm = coeficiente de transferencia de calor por convección del fluido frío.
Km= conductividad térmica del material del que están hechas las placas.
e = espesor de la placa.
Rf y Rc = resistencias térmicas originadas por los contaminantes de los
fluidos que se depositan en las dos superficies de las placas.
Su peso, en porcentaje, es significativo. Estas resistencias térmicas de

ensuciamiento se generan como consecuencia de que los fluidos pueden
transportar contaminantes, y con el paso del tiempo éstos se van depositando
sobre las superficies. De este modo se va creando una capa entre el fluido y la
superficie que va creciendo en espesor generando una resistencia térmica
adicional.
Los coeficientes de convección de los fluidos son función de sus propiedades

físicas y de la geometría de las placas. Se pueden calcular a partir de correlaciones
o bien se pueden obtener de manera aproximada de gráficas como las propuestas
por Cooper y Usher (1983 exclusivas para fluidos Newtonianos y placas de
pequeño tamaño).
112
Obtención de los coeficientes individuales de convección a partir de

correlaciones
Cálculo de coeficiente individual de convección para el fluido caliente

(agua):
Los coeficientes individuales se estiman a partir de correlaciones específicas para

flujo de fluidos entre placas planas, representándose a continuación las más
utilizadas en función de los números adimensionales de Reynolds (Re) y Prandtl
(Pr). (“Ingeniería de la industria alimentaria”, 1999. Conceptos Básicos.
Volumen I. José Aguado. Ed. Síntesis):
Re>400 hc= 0,2536 k/De Re0,65Pr0,4
para régimen turbulento, donde De es el diámetro equivalente, que se define como

cuatro veces el cociente entre la sección de paso del fluido y el perímetro mojado:
De= 4 Sección de paso/Perímetro de mojado= 4 ab/2a+2b ≈ 2b
Donde a es la anchura de la placa plana y b es la distancia entre placas. Como esta

última es muy pequeña respecto a la primera, se puede aproximar el diámetro
equivalente a dos veces la distancia entre placas.
Para régimen laminar, se puede utilizar la correlación:
Re<400 hc= 0,742 CpGRe-0,62Pr-0,667 µ/µ0
donde G, es la velocidad másica ( G = vρ) y µ y µ0 la viscosidad del fluido

evaluada a la temperatura media y a la temperatura de la pared de la placa.
Re = (v lc ρ) /µ
lc = longitud característica del canal

ρ = densidad del fluido
v = velocidad del fluido
µ = viscosidad del fluido
Las propiedades físicas de los fluidos se determinarán para el valor medio de la

temperatura de cada uno.
ρ= 968,3 kg/m3
µ =0,335 mPa s= 0,335×10-3 kg/ms
mc = 8259,673 kg/h = 2,294 kg/s
v= 2,294 kg/s / (4×10-3m × 0,427 m × 968,3 kg/m3) = 1,387 m/s
113
Tanto para el cálculo del Reynolds como para la obtención de los coeficientes de
convección es necesario el cálculo previo de la longitud característica del canal.
Se puede obtener a partir de la expresión para el diámetro hidráulico:
lc = 4 Se / Pc = (4 b W/ 2 (b+W))( 4b W)/2W) = 2b
Se= sección transversal del canal

Pc= perímetro de la sección del canal
b= ancho del canal o distancia entre placas
W= ancho de la placa
b<<<<W
el valor de b dado por el fabricante es de b = 4 mm
Así, lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
Los valores de L y W son los valores de las dimensiones totales sino la altura y la
anchura efectiva de cada placa, es decir, tan solo se tendrán en cuenta el área de
cada placa que pone en contacto ambos fluidos. Esto es debido a que en la placa
existen diversos orificios de entrada y salida de los fluidos y orificios para el paso
de las barras
que hacen presión manteniendo juntas las placas y zonas donde los fluidos no
están puestos en contacto a través de la propia placa.
Las características de la placa las suministra el fabricante, en este caso “APV

Invensys”:
L = 1237,5 mm= 1,2375 m

W = 427 mm = 0,427 m
Superficie de intercambio de una placa (a):
a = ancho por alto= L × W = 1,2375 × 0,427 = 0,528 m2
Para el cálculo del número de Reynolds:
Re = (vlcρ) /µ = (1,387 m/s × 8×10-3 m × 968,3 kg/m3)/ 0,335×10-3 kg/ms =

32072,408
Como en valor obtenido del Re= 32072,408 → Re>400 → Régimen Turbulento

Se utiliza la siguiente expresión:
Re>400 hc= 0.2536 k/De Re0,65Pr0,4
k= 0,673 W/m ºC= 0,673 J/ sm ºC = 0,637×10-3 kJ/ sm ºC

Cp= 4,201 kJ/kg ºC
De ≈ 2b = 2×4 = 8 mm= 8×10-3 m
Pr = (2,19 + 1,96) /2 = 2,075
114
Sustituyendo los valores en la ecuación anterior:
hc= 22,96 kJ/m2s ºC
Cálculo de coeficiente individual de convección para el fluido frío (mix):
Re = (v lc ρ) /µ
lc = longitud característica del canal

ρ = densidad del fluido
v = velocidad del fluido
µ = viscosidad del fluido

ρ= 1100 kg/m3
µ =50 mPa s= 0,05 kg/ms
mc = 5500 kg/h = 1,527 kg/s
v= 1,527 kg/s / (4×10-3m × 0,427 m × 1100 kg/m3) = 0,8127 m/s
Como en el caso anterior:
lc = 4 Se / Pc = (4 b W/ 2 (b+W))( 4b W)/2W) = 2b
Se= sección transversal del canal

Pc= perímetro de la sección del canal
b= ancho del canal o distancia entre placas
W= ancho de la placa = 0,427 m
L = longitud de la placa = 1,2375 m
b<<<<W
Así, lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
Re = (v lc ρ) /µ = (0,8127 m/s × 8×10-3 m × 1100 kg/m3)/ 0,05 kg/ms = 143,035
Como el Re<400→Régimen Laminar; la expresión a emplear será:
hf= 0,742 CpfGRe-0,62Pr-0,667 µ/µ0
115
En este caso el cociente µ/µ0 , se considerará ≈ 1
Cpf= 3,34 kJ/kg ºC

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
Pr =Cp µ / k
Pr= (3,34 kJ/kg ºC × 0,05 kg/ms) / (0,504×10-3 kJ/ sm ºC) = 331,35
hf= 2,128 kJ/m2s ºC
Us= 1/ [1/ hc + e/km + 1/ hf + Rc + Rf]
hc = 22,96 kJ/m2s ºC
hf = 2,128 kJ/m2s ºC
Km= conductividad térmica del material del que están fabricadas las placas.
Para acero inoxidable calidad AISI 316-L:
Km = 13,5 W/m ºC = 13,5 J/msºC= 13,5×10-3 kJ/msºC
e = espesor de la placa = Dato del fabricante cuyo valor es de 0.6×10-3 m

Rc y Rf = resistencias térmicas originadas por los contaminantes de los
fluidos caliente y frío que se depositan en las dos superficies de las placas.
(Dichos valores se tomarán de la nata: “Tecnología de los productos lácteos”
,1998. Ralph Early, al desconocerse los valores de dichas resistencias del mix de
una composición tan concreta (40 % en sólidos totales y 20% en materia grasa)).
Rc = 1,245 ×10-4 sm2ºC/kJ

Rf = 1,904 ×10-4 sm2ºC/kJ
Por tanto:
Us= 1/ [1/ hc + e/km + 1/ hf + Rc + Rf]
Sustituyendo los valores:
Us= 1,79135 kJ/sm2ºC= 1791,35 W/m2ºC
116
Cálculo del Área de transferencia de calor (A) y del número de placas (N):
Del primer apartado se tiene que Q=312290 kJ/h= 86.747 kW
Y como:
Q= UsAΔTm → A= Q/ (UsΔTm) = 86.747×103 W/ (1791,35 W/m2ºC × 5.38 ºC)=
A= 9 m2
En cuanto al área de transferencia se debe considerar el área total de las placas

incluyendo las corrugaciones, o bien el área proyectada. La diferencia entre ambos
criterios puede ser de hasta un 50% mayor para el primer caso. Por simplicidad, se
referirá A al área proyectada de modo que vendrá definida por la expresión:
A= N a = N L W
N= Número de placas útiles para la transferencia de calor que forma el

intercambiador.
a= área proyectada de una única placa.
L= altura de las placas (en la dirección del flujo)
W= anchura de las placas
Como A es la superficie efectiva de transferencia, se tendrá en cuenta que las

placas de los extremos no son útiles para la transferencia de calor de manera que
no se incluirán a la hora de realizar el cálculo del área total.
Los valores de L y W son los valores de las dimensiones totales sino la altura y la
anchura efectiva de cada placa, es decir, tan solo se tendrán en cuenta el área de
cada placa que pone en contacto ambos fluidos. Esto es debido a que en la placa
existen diversos orificios de entrada y salida de los fluidos y orificios para el paso
de las barras que hacen presión manteniendo juntas las placas y zonas donde los
fluidos no están puestos en contacto a través de la propia placa.

Invensys”:
L = 1237,5 mm= 1,2375 m

W = 427 mm = 0,427 m
Cálculo del número de placas:
como N=A/a = 9 m2 / 0,528 m2 = 17,04 ≈ 18 placas
La sección de calentamiento tiene un área total de 9 m2 con un total

de 18 placas.
117
CÁLCULO DE LA PÉRDIDA DE CARGA
En el diseño mecánico de cualquier intercambiador la pérdida de carga es

decisiva, ya que determina la potencia y el consumo de energía en los motores de
las bombas que impulsan los fluidos y por tanto influye en la inversión, en los
costes de operación y en el mantenimiento del sistema de bombeo.
A partir (Manual Perry Chilton) se emplea la siguiente expresión para el cálculo
de la pérdida de carga en intercambiadores de calor de placas (PHE), obtenida a
partir de APVdel manual de transferencia de calor: “Design & Application of
Paraflow-Plate Heat Exchangers”:
ΔP = 2f (G2L)/ (gρDe)
donde f es el denominado factor de fricción:
f = 2,5 (GDe/µ)-0,3
m=caudal másico del fluido (kg/s)

G= vρ= velocidad másica (kg /m2s)
Cálculo de la pérdida de carga para el fluido caliente (agua)
donde f = 2,5 (GzDe/µ)-0,3
mc=8259,673 kg/h= 2,29 kg/s

G= vρ= 1,387 m/s × 968,3 kg/m3 = 2686,0642 kg/m2s
De= lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
µ= 0,335×10-3 kg/ms
L=1,2375 m
ρ=densidad del fluido caliente= 968,3 kg/m3
g=aceleración de la gravedad= 9,8 m/s2
Sustituyendo los datos en las correspondientes ecuaciones:
f = 2,5 (GDe/µ)-0,3=0,09032
ΔPc = 2f (G2L)/ (gρDe) =21245,82 kg/m2= 2,083 bar
La pérdida de carga que se produce en el fluido calefactor (agua caliente) es

de 2,083 bar
118
Cálculo de la pérdida de carga para el fluido frío (mix)
donde f = 2,5 (GDe/µ)-0,3
mf=5500 kg/h= 1,527 kg/s

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
De= lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
µ= 0,05 kg/ms
L=1,2375 m
ρ=densidad del fluido caliente= 1100 kg/m3
f = 2.5 (GDe/µ)-0.3=0,564
ΔPf = 2f (G2L)/ (gρDe) =12935,74 kg/m2= 1,268 bar
La pérdida de carga que se produce en el mix sin pasteurizar es de 1,268 bar
119
Cálculo de la Sección I (Sección de Regenerativa)
Cálculo del % de regeneración

Como se ha citado en otras ocasiones, el método que consiste en utilizar el calor
de un líquido caliente, como el mix pasteurizado, para precalentar el líquido frío
(mix sin pasteurizar), reduciéndose así el consumo de energía, se denomina
regeneración. Según el “Manual de Industrias Lácteas”, 1996. Tetra Pak:
R= ((tr – ti) × 100)/ (tp – ti)
donde:
R= eficiencia de la regeneración, %.
tr= temperatura del fluido después de la regeneración ( en este caso = 67 ºC)
ti= temperatura de entrada del fluido sin pasteurizar (en este caso = 50 ºC)
tp= temperatura de pasteurización (en este caso = 84 ºC)
así queda:
R= ((67 – 50)×100)/(84 – 50) = 50 %
Dado que el mix antes de entrar en la sección regenerativa del pasteurizador, ya

había sido previamente precalentado, el % de recuperación de calor no es del
orden del 80 al 90 %, que se podría haber conseguido.
En esta sección las propiedades del fluido y la relación de temperaturas es la

siguiente:
Propiedades del
fluido Unidades Caliente Frío
Mix Mix sin
pasteurizado pasteurizar
Densidad kg/m3 1100 1100
Conductividad
Térmica W/m ºC 0,673 0,504
Viscosidad Entrada mPa s 42 104,88
Viscosidad Salida mPa s 58 58
Tfe = Temperatura del mix sin pasteurizar a la entrada: 50 ºC

Tfs = Temperatura del mix sin pasteurizar a la salida: 67 ºC
Tce= Temperatura del mix pasteurizado a la entrada: 84 ºC
Tcs =Temperatura del mix pasteurizado a la salida: 70 ºC
120
de esta ecuación la variable mc=mf, es decir, los caudales másicos de ambos

fluidos coinciden.
Cálculo del calor intercambiado:
Q = mf Cpf (Tfs – Tfe)= 5500 kg/h ×3,34 kJ/kg ºC× (67 ºC -50 ºC)= 312290 kJ/h
Q =86,747 kW
Diferencia media de temperatura (ΔTm) se define como:
ΔTm = [(84 – 67) - (70 – 50) / ln ((84– 67) / (70 – 50))] = 18,46 ºC
Es necesario introducir un factor F, menor que la unidad, quedando que la

transferencia de calor a lo largo del intercambiador es:
Q= UAΔTm= UAFΔTlm
R= (Tce – Tcs) / (Tfs – Tfe) = (84 – 70) / (67 – 50) = 0,823
P= (Tfs – Tfe)/ (Tce – Tfe) = (67 – 50) / (84 – 50) = 0,5
ΔTm= F ΔTlm

F≈0,88, así:
ΔTm= F ΔTlm = 0,88×18,46 ºC= 16,245 ºC
Us= 1/ [1/ hc + e/km + 1/ hf + Rf + Rc]
121

(mix pasteurizado):
Cálculo del número de Reynolds:
Re = (v lc ρ) /µ
ρ= 1100 kg/m3
µ =50 mPa s= 0,05 kg/ms
v= 1,527 kg/s / (4×10-3 m × 0,427 m × 1100 kg/m3) = 0,8127 m/s
Así, lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
hc= 0,742 CpfGRe-0,62Pr-0,667 µ/µ0
En este caso el cociente µ/µ0 , se considerará ≈ 1
Cpc= 3,34 kJ/kg ºC

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
Pr =Cpc µ / k siendo k: conductividad térmica del fluido caliente (mix)
hc= 2,58 kJ/m2s ºC
122
Cálculo de coeficiente individual de convección para el fluido frío (mix

sin pasteurizar):
Re = (v lc ρ) /µ
ρ= 1100 kg/m3
µ = (104,88 mPa s + 58 mPa s) /2 = 81,44 mPa s= 0,08144 kg/ms
v= 1,527 kg/s / (4×10-3 m × 0,427 m × 1100 kg/m3) = 0,8127 m/s
lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
hf= 0,742 CpfGRe-0,62Pr-0,667 µ/µ0
Se considerará: µ/µ0 ≈ 1
Cpf= 3,34 kJ/kg ºC

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
Pr =Cpf µ / k siendo k: conductividad térmica del fluido frio (mix sin
pasteurizar)
Us= 1/ [1/ hc + e/km + 1/ hf + Rc + Rf]
123
Kw = 13,5 W/m ºC = 13,5 J/msºC= 13,5×10-3 kJ/msºC

e = espesor de la placa = 0,6×10-3 m
Rc = 5,20 ×10-5 m2k/W= 1,904 ×10-4 sm2ºC/kJ
Rf = 5,20 ×10-5 m2k/W= 1,904 ×10-4 sm2ºC/kJ
hc= 2,58 kJ/m2sºC
hf= 2,080 kJ/m2sºC
Por tanto:
Us= 1/ [1/ hc + e/km + 1/ hf + Rc + Rf]
Sustituyendo los datos en la ecuación anterior:
Us= 1,09506 kJ/sm2ºC= 1095,06 W/m2ºC
Del primer apartado se tiene que Q=312290 kJ/h= 86,747 kW
Despejando A:
Q= UsAΔTm → A= Q/ (UsΔTm) = 86,747×103 W/ (1095,06 W/m2ºC × 16,245

ºC)=
A=4,87 m2

Invensys”:
L = 1237,5 mm= 1,2375 m

W = 427 mm = 0,427 m
y como N=A/a= 4,87 m2 / 0,528 m2 = 9,23 ≈ 10 placas
La sección regenerativa tiene un área total de 4,87 m2 con un total

de 10 placas.
124
CÁLCULO DE LA PÉRDIDA DE CARGA
Cálculo de la pérdida de carga para el fluido caliente (mix pasteurizado)
donde f = 2.5 (GDe/µ)-0.3
mc=5500 kg/h= 1,527 kg/s

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
De= lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
µ= 0,05 kg/ms
L=1,2375 m
ρ=densidad del fluido= 1100 kg/m3
f = 2,5 (GDe/µ)-0.3=0,564
ΔPc= 2f (G2L)/ (gρDe) =12935,74 kg/m2= 1,268 bar
La pérdida de carga que se produce en el mix pasteurizado es de 1,268 bar
Cálculo de la pérdida de carga para el fluido frío (mix sin pasteurizar)
donde f = 2.5 (GDe/µ)-0.3
mc=5500 kg/h= 1,527 kg/s

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
De= lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
µ= 0,08144 kg/ms
L=1,2375 m
ρ=densidad del fluido= 1100 kg/m3
f = 2,5 (GDe/µ)-0.3=0,653
ΔPc= 2f (G2L)/ (gρDe) =14977,024 kg/m2= 1,468 bar
La pérdida de carga que se produce en el mix sin pasteurizar es de 1,468 bar
125
Cálculo de la Sección 5ª (II) (Sección de enfriamiento)
En esta sección las propiedades del fluido y la relación de temperaturas es la

siguiente:
Propiedades del fluido Unidades Caliente Frío

Mix Agua
pasteurizado fría
Densidad kg/m3 1100 995.9
Conductividad Térmica W/m ºC 0,47 0,612
Viscosidad Entrada mPa s 58 1,05
Viscosidad Salida mPa s 249,31 0,65
Tfe = Temperatura del agua a la entrada: 18 ºC

Tfs = Temperatura del agua a la salida: 40 ºC
Tce= Temperatura del mix pasteurizado a la entrada: 67 ºC
Tcs =Temperatura del mix pasteurizado a la salida: 28 ºC
mf: Caudal másico del fluido frío(agua de red) kg/h: desconocido
Cpf: Calor específico del fluido frío (agua de red) = 4,179 kJ/kg ºC
mc: Caudal másico del fluido caliente (mix pasteurizado) kg/h como la
ρmix= 1100 kg/m3, para una qent= 5000 l/h, se tiene que mf= 5500 kg/h
Cpc: Calor específico del fluido caliente (mix pasteurizado) = 3,34 kJ/kg ºC
de esta ecuación sólo se desconoce la variable mf:
donde mf=[ mcCpc (Tcs – Tce)]/ Cpf (Tfe – Tfs)]
mf= [5500 kg/h ×3,34 kJ/kg ºC×(28 ºC -67 ºC)]/ [4,179 kJ/kg ºC×(18 ºC- 40 ºC)]
=
mf= 7792,534 kg/h=2,16 kg/s como ρf = 995,9 kg/m3
qf= 7824,625 l/h
Cálculo del calor intercambiado:
Q = mf Cpf (Tfs – Tfe) = 7792,534 kg/h ×4,179 kJ/kg ºC× (40 ºC -18 ºC)=
Q = 716429,99 kJ/h
Q = 199 kW
126
Diferencia media de temperatura (ΔTm) se define como:
ΔTm = [(67 – 40) - (28 – 18) / ln ((67– 40) / (28 – 18))] = 17,11 ºC
Es necesario introducir un factor F, menor que la unidad, quedando que la

transferencia de calor a lo largo del intercambiador es:
Q= UAΔTm= UAFΔTlm
R= (Tce – Tcs) / (Tfs – Tfe) = (67 – 28) / (40 – 18) = 1,77
P= (Tfs – Tfe)/ (Tce – Tfe) = (40 – 18) / (67 – 18) = 0,45
ΔTm= F ΔTlm

F≈0,53, así:
ΔTm= F ΔTlm = 0,53×17,11= 9,068 ºC
Us= 1/ [1/ hc + e/km + 1/ hf + Rf + Rc]

(mix pasteurizado):
Cálculo del número de Reynolds:
Re = (v lc ρ) /µ
ρ= 1100 kg/m3
µ = (58 + 249,31)/2 = 153,65 mPas= 0,1536 kg/ms
v= 1,527 kg/s / (4×10-3 m × 0,427 m × 1100 kg/m3) = 0,8127 m/s
127
Así, lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
hc= 0,742 CpfGRe-0,62Pr-0,667 µ/µ0
Se considerará el cociente µ/µ0 ≈ 1
Cpc= 3,34 kJ/kg ºC

G= vρ= 0,8127 m/s × 1100 kg/m3 = 893,97 kg/m2s
Pr =Cpc µ / k
Se obtiene:
hc= 1,927 kJ/m2s ºC
Cálculo de coeficiente individual de convección para el fluido frío (agua

fría):
Re = (v lc ρ) /µ

mf= 7792,534 kg/h=2,16 kg/s

ρ= 995,9 kg/m3
µ = (1,05 mPa s + 0,65 mPa s) /2 = 0,85 mPa s= 8,5×10-4 kg/ms
v= 2,16 kg/s / (4×10-3 m × 0,427 m × 1100 kg/m3) = 1,149 m/s
Así, lc ≈ 2b = 2×4 = 8 mm= 8×10-3 m
Re = (v lc ρ) /µ = (1,149 m/s × 8×10-3 m × 995,9 kg/m3)/ 8,5×10-4 kg/ms =

10769,78
128
Como en valor obtenido del Re= 32072,408 → Re>400 → Régimen Turbulento

Se utiliza la siguiente expresión:
Re>400 hc= 0,2536 k/De Re0,65Pr0,4

De ≈ 2b = 2×4 = 8 mm= 8×10-3 m
G= vρ= 1,149 m/s × 995,9 kg/m3 = 1144,289 kg/m2s
Pr = (7,31+2,75) /2 = 5,03
Sustituyendo los datos en la ecuación anterior:
Us= 1/ [1/ hc + e/km + 1/ hf + Rc + Rf]
Rf = 3,40 ×10-5 m2k/W= 1,245 ×10-4 sm2ºC/kJ

Rc = 5,20 ×10-5 m2k/W= 1,904 ×10-4 sm2ºC/kJ
hc= 1,927 kJ/m2sºC
hf= 15,46 kJ/m2sºC
Us= 1,59138 kJ/sm2ºC= 1591,38 W/m2ºC
Del apartado anterior se tiene que Q= 716429,99 kJ/h= 199 kW,
Despejando A:
Q= UsAΔTm → A= Q/ (UsΔTm) = 199 ×103 W/ (1591,38 W/m2ºC ×9,068 ºC)=
A=13,79 m2

Invensys”:
L = 1237,5 mm= 1,2375 m

W = 427 mm = 0,427 m
129
y como N=A/a= 13,79 m2 / 0,528 m2 = 26,11 ≈ 27 placas
La primera sección de enfriamiento tiene un área total de 26,11 m2 con un

total de 27 placas.
130

Diseño Pasteurizador

Cargado por

Copyright:

Formatos disponibles

Diseño Pasteurizador

Cargado por

Información del documento

Descripción original:

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

Diseño Pasteurizador

Cargado por

Copyright:

Formatos disponibles

CÁLCULOS JUSTIFICATIVOS 10.

Diseño del pasteurizador

10.2. Cálculo y diseño del intercambiador de calor de

10.2.2. Hipótesis de cálculo

Para el diseño del intercambiador de calor para la pasteurización de Qmix= 5000

1. El intercambiador de calor opera en condiciones estacionarias.

como ocurre con otros productos de composición constante como la leche,

Cálculo de la Sección III (Sección de calentamiento)

La ecuación de energía, suponiendo que no existen pérdidas térmicas es la

ΔH + ΔEc + ΔEp = Q + We (Balance energético)

Para que la expresión anterior queda:

según las siguientes simplificaciones:

Formas mecánicas de energía son nulas ΔEc ≅ ΔEp ≅ We = 0

Para cualquiera que sea el tipo de equipo intercambiador empleado, si sólo se

Q= mf (h1 – h2)= mc (H2 – H1)

mfCpf (Tfe – Tfs) = mcCpc (Tcs -Tce)

Propiedades del fluido Unidades Caliente Frío

Tfe = Temperatura del mix sin pasteurizar a la entrada: 67 ºC

mf: Caudal másico del fluido frío(mix a pasteurizar) Kg/h como la

mc: Caudal másico del fluido caliente (agua caliente): desconocido

de esta ecuación sólo se desconoce la variable mc:

donde mc=[ mfCpf (Tfe – Tfs)]/ Cpc (Tcs – Tce)]

mc= 8259,673 kg/h = 2,294 kg/s como ρc= 968,3 kg/m3

qc= 8,53 m3/h

Aplicando a continuación la ecuación de Fourier a todo el equipo:

De la ecuación tanto Q como ΔTm son valores fácilmente calculables.

Cálculo del calor intercambiado (Q) se define como:

Diferencia media de temperatura (ΔTm):

ΔTm = ΔT1 - ΔT2 / (ln (ΔT1 / ΔT2))

Los intercambiadores de placas se diseñan para funcionar con el flujo de los

ΔTm = [(Tce – Tfs) – (Tcs – Tfe) / ln ((Tce - Tfs) / (Tcs – Tfe))]

Los canales de los extremos del intercambiador presentan la mitad de superficie

Para el cálculo del factor de corrección F es necesario definir la llamada relación

F: Factor de corrección (o de aproximación)

P= (Tfs – Tfe)/ (Tce – Tfe) = (84 – 67) / (90 – 67) = 0,739

El efecto de la disminución de la efectividad, debido a que la superficie de

Sin embargo si el número de placas es menor a 50, el factor F se calcula

A partir de la gráfica para el cálculo del factor de corrección para

ΔTm= F ΔTlm = 0,57×9,44 ºC= 5,38 ºC

Obtención del coeficiente global de transmisión de calor (Us sucio)

Cálculo del coeficiente de calor sucio (Us)

El coeficiente global de transferencia de calor viene dado por la expresión:

Us= 1/ [1/ hc + e/km + 1/ hf + Rf + Rc]

hc = coeficiente de transferencia de calor por convección del fluido caliente.

Su peso, en porcentaje, es significativo. Estas resistencias térmicas de

Los coeficientes de convección de los fluidos son función de sus propiedades

Obtención de los coeficientes individuales de convección a partir de

Cálculo de coeficiente individual de convección para el fluido caliente

Los coeficientes individuales se estiman a partir de correlaciones específicas para

Re>400 hc= 0,2536 k/De Re0,65Pr0,4

para régimen turbulento, donde De es el diámetro equivalente, que se define como

De= 4 Sección de paso/Perímetro de mojado= 4 ab/2a+2b ≈ 2b

Donde a es la anchura de la placa plana y b es la distancia entre placas. Como esta

Para régimen laminar, se puede utilizar la correlación:

Re<400 hc= 0,742 CpGRe-0,62Pr-0,667 µ/µ0

donde G, es la velocidad másica ( G = vρ) y µ y µ0 la viscosidad del fluido

lc = longitud característica del canal

Las propiedades físicas de los fluidos se determinarán para el valor medio de la