See discussions, stats, and author profiles for this publication at: [Link]

net/publication/276276204

Comparación entre modelos matemáticos y estimación de coeficientes de

transferencia de calor para un intercambiador de calor.

Conference Paper · May 2015

CITATIONS READS

0 816

5 authors, including:

D. A. Carbot-Rojas Marcos Alfredo López Castro

Centro Nacional de Metrologia Centro Nacional de Investigación y Desarrollo Tecnológico
6 PUBLICATIONS   50 CITATIONS    4 PUBLICATIONS   9 CITATIONS   

SEE PROFILE SEE PROFILE

R. F. Escobar Jiménez J.F. Gómez-Aguilar

Centro Nacional de Investigación y Desarrollo Tecnológico CONACyT-Tecnológico Nacional de México/CENIDET
109 PUBLICATIONS   1,220 CITATIONS    275 PUBLICATIONS   3,761 CITATIONS   

SEE PROFILE SEE PROFILE

Some of the authors of this publication are also working on these related projects:

PEMFC Supercapacitor Battery View project

Calculo fraccionario View project

All content following this page was uploaded by R. F. Escobar Jiménez on 11 June 2015.

The user has requested enhancement of the downloaded file.

 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

COMPARACIÓN ENTRE MODELOS MATEMÁTICOS Y ESTIMACIÓN DE

COEFICIENTES DE TRANSFERENCIA DE CALOR PARA UN INTERCAMBIADOR DE
CALOR
Diego Alessis Carbot Rojasa, Marcos Alfredo López Castroa, R. F. Escobarb, J .F. Gómez-Aguilarb, J. A. Hernández Pérezc
a
Posgrado en Ingeniería Electrónica, Centro Nacional de Investigación y Desarrollo Tecnológico, Interior Internado Palmira
S/N, Palmira C.P. 62490, Cuernavaca, Mor., México.
b
Centro Nacional de Investigación y Desarrollo Tecnológico - Tecnológico Nacional de México, Interior Internado Palmira
S/N, Palmira C.P. 62490, Cuernavaca, Mor., México.
c
Centro de Investigación en Ingeniería y Ciencias Aplicadas, Universidad Autónoma del Estado de Morelos, Av. Universidad
1001, Col. Chamilpa, C.P. 62209 Cuernavaca, Morelos, México
Resumen
En este trabajo se presenta la comparación de dos modelos matemáticos para la estimación de
temperaturas en un intercambiador de calor de tubos concéntricos con flujos a contra corriente y sin
cambio de fase. Uno de los modelos es representado por ecuaciones diferenciales ordinarias y el otro
modelo por ecuaciones diferenciales parciales. Para realizar la estimación de las temperaturas en el
modelo de ecuaciones diferenciales ordinarias, se consideraron dos diferentes métodos para el cálculo
del coeficiente global de temperatura en modelo termodinámico, el primero es descrito por una sola
ecuación que describe de forma global la transferencia de calor entre ambas corrientes y el segundo
método considera el cálculo del coeficiente global de transferencia de calor para cada una de las
corrientes. Por otra parte, en el modelo de ecuaciones diferenciales parciales se empleó un modelo
termodinámico donde se calcula el coeficiente de transferencia de calor convectivo para cada una de las
corrientes del intercambiador de calor (interna y externa).
Palabras clave: Estimación de temperatura, Intercambiador de calor, modelado matemático.
Introducción
Los modelos de parámetros agrupados son ampliamente usados para describir la dinámica de un sistema
en función del tiempo, así como, para realizar el diseño de un controlador, su representación es a través
de ecuaciones diferenciales ordinarias (EDO). Por otro lado, los modelos de parámetros distribuidos son
aquellos que toman en cuenta los cambios de las variables en tiempo y en espacio y son representados
por un conjunto de ecuaciones diferenciales parciales (EDP).
En los trabajos presentados en [1-3] se propone el uso del modelo matemático que representa la
transferencia de calor mediante el uso de un par de ecuaciones diferenciales ordinarias, considerando el
coeficiente global de transferencia de calor (𝑈𝑈) constante, así como las propiedades físicas del agua (𝜌𝜌 y
𝐶𝐶𝐶𝐶). Sin embargo, en [4] y [5] se mostró que para que el modelo del intercambiador de calor fuera
válido en un amplio rango de operación, el coeficiente de transferencia de calor tenía que ser estimado y
que éste dependía del flujo del fluido y de la temperatura. Por lo tanto, si se trabaja con un modelo que
considere el coeficiente global de transferencia de calor como constante, el modelo únicamente será
válido en una vecindad de valores cercanos a los valores considerados para el cálculo del coeficiente
global de transferencia de calor. En [6] se utiliza un modelo matemático que representa la transferencia
de calor mediante ecuaciones diferenciales parciales en el cual los coeficientes convectivos (ℎ𝑖𝑖, ℎ𝑒𝑒) se
mantienen constantes, y al usar los coeficientes convectivos como constantes el modelo tienen un error
mayor en la estimación de las temperaturas de salida. Para desarrollar la comparación entre los dos
modelos matemáticos se consideraron los dos propuestos en la literatura [1] y [6], el modelo empleado
en [1], considera el coeficiente global de transferencia de calor constante, mientras que para el caso del
modelo propuesto en [6] se emplean los coeficientes convectivos de transferencia de calor para cada
corriente, en donde los coeficientes convectivos son considerados como constantes.
En este trabajo se realizará la comparación de estos dos modelos dinámicos [1] y [6], con una
modificación en el modelo termodinámico de ambos modelos dinámicos, esta modificación consiste en

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
3293
 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

estimar en línea tanto el coeficiente global de transferencia de calor para el caso del modelo presentado
en [1], así como realizar la estimación en línea de los coeficientes convectivos para el caso del modelo
presentado en [6]. Para el caso del modelo que emplea el coeficiente global de transferencia de calor (U)
se proponen dos diferentes enfoques, el primero considera el cálculo de un sólo coeficiente global de
transferencia de calor para el modelo dinámico y el segundo considera un coeficiente global para cada
corriente. Para realizar este trabajo se consideran las propiedades físicas de fluido (densidad 𝜌𝜌(𝑇𝑇, 𝑃𝑃) y
calor específico 𝐶𝐶𝐶𝐶(𝑇𝑇, 𝑃𝑃)) dependientes de la presión y de la temperatura, y el coeficiente global de
transferencia de calor y los coeficientes convectivos se consideran dependientes del flujo y de la
temperatura [7].
Tabla 1. Nomenclatura. Subíndices

𝑽𝑽 Volumen, 𝑚𝑚3 𝒇𝒇 Frio

𝑭𝑭 Flujo volumétrico, 𝑚𝑚3 /𝑠𝑠 𝒄𝒄 Caliente
𝑪𝑪𝒑𝒑 Calor especifico del fluido,
𝐽𝐽
°𝐶𝐶 𝒆𝒆𝒆𝒆𝒆𝒆 Entrada
𝐾𝐾𝐾𝐾 𝒔𝒔𝒔𝒔𝒔𝒔 Salida
𝑨𝑨 Área de transferencia de calor, 𝑚𝑚2 𝒆𝒆 Externo
𝑼𝑼 Coeficiente global de transferencia 𝒊𝒊 Interno
𝑊𝑊
de calor, 2 °𝐶𝐶
𝑚𝑚 𝒑𝒑 Pared
𝐓𝐓 Temperatura, °𝐶𝐶
𝒓𝒓 Radio , 𝑚𝑚 𝒙𝒙 Posición axial
𝒉𝒉 Coeficiente convectivo, 𝑊𝑊/𝑚𝑚2 𝐾𝐾 Números adimensionales
𝑳𝑳 Longitud del tubo, 𝑚𝑚
𝑽𝑽𝑽𝑽𝑽𝑽 Velocidad del fluido, 𝑚𝑚/𝑠𝑠 𝑵𝑵𝑵𝑵 Número de Nusselt
𝑫𝑫 Diámetro, 𝑚𝑚 𝑹𝑹𝑹𝑹 Número de Reynolds
𝒕𝒕 Tiempo, 𝑠𝑠 𝑷𝑷𝑷𝑷 Número de Prandtl
𝑬𝑬 Espesor de la pared, 𝑚𝑚

Letras griegas

𝝆𝝆 Densidad del fluido, 𝑘𝑘𝑘𝑘/𝑚𝑚3

𝝀𝝀 Coeficiente de conductividad térmica, 𝑊𝑊/𝑚𝑚𝑚𝑚
𝝁𝝁 Viscosidad dinámica, 𝑘𝑘𝑘𝑘/𝑚𝑚𝑚𝑚

Metodología
La Figura 1 muestra el esquema del intercambiador de calor bajo estudio, la Tabla 1 muestra las
mediciones físicas de los tubos del intercambiador de calor. Para efectos de cálculo del volumen y área
de transferencia de calor del intercambiador.
Tabla 2. Dimensiones del intercambiador de calor.
Constantes Valores Unidades
𝐿𝐿 790 mm
𝐷𝐷𝑖𝑖𝑖𝑖 6 mm
𝐷𝐷𝑒𝑒𝑒𝑒 10 mm
𝐷𝐷𝑖𝑖𝑖𝑖 16 mm
𝐷𝐷𝑒𝑒𝑒𝑒 20 mm
Figura 1. Dimensiones de los tubos del intercambiador de calor.

Modelo dinámico de parámetros agrupados para el intercambiador de calor

La dinámica del intercambiador de calor se obtiene por un balance de energía para cada uno de los lados
del intercambiador de calor, dada en la ecuación (1)

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
3294
 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

𝑑𝑑𝑑𝑑𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 𝐹𝐹𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 𝑈𝑈𝑈𝑈

= �𝑇𝑇𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 − 𝑇𝑇𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 (𝑡𝑡)� + �𝑇𝑇 (𝑡𝑡) − 𝑇𝑇𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓 (𝑡𝑡)�
𝑑𝑑𝑑𝑑 𝑉𝑉𝑓𝑓 𝜌𝜌𝑓𝑓 𝐶𝐶𝑝𝑝𝑝𝑝 𝑉𝑉𝑓𝑓 𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 (1)
𝑑𝑑𝑑𝑑𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝐹𝐹𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 𝑈𝑈𝑈𝑈
= (𝑇𝑇𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 − 𝑇𝑇𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 (𝑡𝑡)) + �𝑇𝑇 (𝑡𝑡) − 𝑇𝑇𝑐𝑐𝑐𝑐𝑐𝑐𝑐𝑐 (𝑡𝑡)�
𝑑𝑑𝑑𝑑 𝑉𝑉𝑐𝑐 𝜌𝜌𝑐𝑐 𝐶𝐶𝑝𝑝𝑝𝑝 𝑉𝑉𝑐𝑐 𝑓𝑓𝑓𝑓𝑓𝑓𝑓𝑓
El modelo matemático, que se muestra en la ecuación (1), tiene en cuenta las siguientes suposiciones:
• Volumen constante en ambos tubos.
• El coeficiente de transferencia de calor está relacionado con las temperaturas de los fluidos y con
el flujo y no es considerado constante.
• No existe transferencia de calor entre el tubo externo y el medio ambiente.
• No hay almacenamiento de energía calorífica en las paredes de los tubos.
• Las temperaturas de entrada son medibles.
Cálculo del coeficiente de transferencia de calor (U)
Para estimar las temperaturas en el intercambiador de calor con precisión, en este trabajo se propone
calcular el coeficiente de transferencia de calor (𝑈𝑈) en línea (con datos experimentales en cada tiempo),
a través de los coeficientes convectivos (ℎ𝑖𝑖, ℎ𝑒𝑒), de aquí en adelante, el subíndice 𝜑𝜑 se utiliza para
referirse al lado externo o interno.
El coeficiente convectivo (ℎ𝜑𝜑 ) puede ser calculado para cada lado del intercambiador de calor usando:
𝑁𝑁𝑁𝑁𝜑𝜑
ℎ𝜑𝜑 = (2)
𝐷𝐷/𝜆𝜆𝜑𝜑
Para una sola fase, el coeficiente convectivo se determina utilizando el número de Nusselt (𝑁𝑁𝑁𝑁), que se
calcula para cada lado del intercambiador de calor usando las siguientes formulas (correlación Dittus &
Boelter):
Para flujos laminares: Para flujos turbulentos:
𝑁𝑁𝑁𝑁𝜑𝜑 = 0.664𝑅𝑅𝑅𝑅𝜑𝜑 𝑃𝑃𝑃𝑃𝜑𝜑 0.5 0.3
(3) 𝑁𝑁𝑁𝑁𝜑𝜑 = 0.02𝑅𝑅𝑅𝑅𝜑𝜑 0.8 𝑃𝑃𝑃𝑃𝜑𝜑 0.4 (4)
Donde 𝑃𝑃𝑃𝑃 es el número de Prandtl que se calcula con la ecuación (5) y 𝑅𝑅𝑅𝑅 es el número de Reynolds que
se calcula con la ecuación (6),
𝐶𝐶𝐶𝐶𝐶𝐶 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌
𝑃𝑃𝑃𝑃 = (5) 𝑅𝑅𝑅𝑅 = (6)
𝜆𝜆 𝜇𝜇
El coeficiente global de transferencia de calor (𝑈𝑈) se puede calcular mediante dos métodos, el primer
método es descrito en (7), el cual es representado por una sola ecuación que describe de forma global la
transferencia de calor que existe entre ambos fluidos [8], y el segundo método es el descrito en (8) donde
el cálculo del coeficiente global de transferencia de calor se considera para cada uno de los fluidos [8].
En este trabajo se realiza la estimación de las temperaturas empleando ambos métodos.
1
𝑈𝑈 = 𝑟𝑟
𝑟𝑟𝑒𝑒 ln � 𝑟𝑟𝑒𝑒 � (7)
1 𝑖𝑖 1
+ +
ℎ𝑒𝑒 𝜆𝜆 ℎ𝑖𝑖
1 1
𝑈𝑈𝑒𝑒 = , 𝑈𝑈𝑖𝑖 =
𝐴𝐴𝑒𝑒 𝐴𝐴𝑒𝑒 𝐴𝐴 1 𝐴𝐴𝑖𝑖 𝐴𝐴𝑖𝑖 𝐴𝐴 1 (8)
+ ∗ ln � 𝐴𝐴𝑒𝑒 � + + ∗ ln �𝐴𝐴 𝑖𝑖 � +
𝐴𝐴𝑖𝑖 ∗ ℎ𝑖𝑖 𝜆𝜆 ∗ 2𝜋𝜋 ∗ 𝐿𝐿 𝑖𝑖 ℎ𝑒𝑒 𝐴𝐴𝑒𝑒 ∗ ℎ𝑒𝑒 𝜆𝜆 ∗ 2𝜋𝜋 ∗ 𝐿𝐿 𝑒𝑒 ℎ𝑖𝑖

Modelo dinámico de parámetros distribuidos para el intercambiador de calor

Otro modelo propuesto es el presentado en [6] en el cual se asumen las siguientes
suposiciones:
(i) Se determina el tipo de flujo dependiendo el numero de Reynolds.

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
3302
3295
 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

(ii) las propiedades térmicas de las paredes de las tuberías y los fluidos se suponen constantes,
(iii) con el fin de obtener un modelo del intercambiador de calor se considera un volumen
elemental diferencial,
(iv) Las temperaturas de los fluidos y la pared de separación se suponen constantes en esta
sección diferencial de los tubos.
Basado en los supuestos anteriores, la transferencia de calor de transición entre los dos fluidos y la pared
de separación se rige por un modelo puramente advectivo como se indica a continuación:
𝜕𝜕𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡) 𝜕𝜕𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡)
− 𝑉𝑉𝑐𝑐 = 𝐻𝐻𝑐𝑐 �𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡) − 𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡)�
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡) 𝜕𝜕𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡)
+ 𝑉𝑉𝑓𝑓 = 𝐻𝐻𝑓𝑓 �𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡) − 𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡)� (9)
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡)
= 𝐻𝐻𝑐𝑐𝑐𝑐 �𝑇𝑇ℎ (𝑥𝑥, 𝑡𝑡) − 𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡)� + 𝐻𝐻𝑓𝑓𝑓𝑓 �𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡) − 𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡)�
𝜕𝜕𝜕𝜕
Con condiciones iniciales: (0, 𝐿𝐿 ) = y (0,0) = .

Aplicando el método de diferencias finitas para el sistema de ecuaciones (9), el modelo

del intercambiador de calor se puede escribir como sigue:
𝑑𝑑𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡) 𝑉𝑉𝑐𝑐 𝑉𝑉𝑐𝑐
= �− − 𝐻𝐻𝑐𝑐 � 𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡) + 𝐻𝐻𝑐𝑐 𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡) + 𝑇𝑇 (𝑥𝑥 + ∆𝑥𝑥, 𝑡𝑡)
𝑑𝑑𝑑𝑑 ∆𝑥𝑥 ∆𝑥𝑥 𝑐𝑐
𝑑𝑑𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡) 𝑉𝑉𝑓𝑓 𝑉𝑉𝑓𝑓
= �− − 𝐻𝐻𝑓𝑓 � 𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡) + 𝐻𝐻𝑓𝑓 𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡) + 𝑇𝑇 (𝑥𝑥 + ∆𝑥𝑥, 𝑡𝑡)
𝑑𝑑𝑑𝑑 ∆𝑥𝑥 ∆𝑥𝑥 𝑓𝑓 (10)
𝑑𝑑𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡)
= 𝐻𝐻𝑐𝑐𝑐𝑐 𝑇𝑇𝑐𝑐 (𝑥𝑥, 𝑡𝑡) + 𝐻𝐻𝑓𝑓𝑓𝑓 𝑇𝑇𝑓𝑓 (𝑥𝑥, 𝑡𝑡) − �𝐻𝐻𝑐𝑐𝑐𝑐 + 𝐻𝐻𝑓𝑓𝑓𝑓 �𝑇𝑇𝑝𝑝 (𝑥𝑥, 𝑡𝑡)
𝑑𝑑𝑑𝑑
𝐻𝐻𝑐𝑐 , 𝐻𝐻𝑓𝑓 , 𝐻𝐻𝑐𝑐𝑝𝑝 y 𝐻𝐻𝑓𝑓𝑓𝑓 son parámetros de acoplamiento, que están en función de los coeficientes de
transferencia de calor por convección y se calculan como se muestra a continuación.
2ℎ𝑖𝑖 𝑟𝑟𝑖𝑖 2ℎ𝑒𝑒 (𝑟𝑟𝑖𝑖 + 𝐸𝐸𝑖𝑖 )
𝐻𝐻𝑐𝑐𝑐𝑐 = 𝐻𝐻 𝑓𝑓𝑓𝑓 =
𝜌𝜌𝑝𝑝 𝐶𝐶𝐶𝐶𝑝𝑝 ((𝑟𝑟𝑖𝑖 + 𝐸𝐸𝐸𝐸𝑖𝑖 )2 − 𝑟𝑟𝑖𝑖 2 ) 𝜌𝜌𝑝𝑝 𝐶𝐶𝐶𝐶𝑝𝑝 ((𝑟𝑟𝑖𝑖 + 𝑒𝑒𝑖𝑖 )2 − 𝑟𝑟𝑖𝑖 2 )

2ℎ𝑖𝑖 2ℎ𝑒𝑒 (𝑟𝑟𝑖𝑖 + 𝐸𝐸𝑖𝑖 )

𝐻𝐻𝑐𝑐 = 𝐻𝐻𝑓𝑓 =
𝜌𝜌𝑐𝑐 𝐶𝐶𝐶𝐶𝑐𝑐 𝑟𝑟𝑖𝑖 𝜌𝜌𝑓𝑓 𝐶𝐶𝐶𝐶𝑓𝑓 (𝑟𝑟𝑒𝑒2 − (𝑟𝑟𝑖𝑖 + 𝐸𝐸𝑖𝑖 )2 )

Cálculo de los coeficientes convectivos

Como se mencionó anteriormente para estimar las temperaturas de salida con una mejor precisión, se
calculan en línea los coeficientes convectivos (ℎ𝑖𝑖, ℎ𝑒𝑒).
En la práctica, la mayoría de las predicciones del coeficiente de transferencia de calor por convección
ℎ𝑖𝑖, ℎ𝑒𝑒 se basan en correlaciones que han sido derivados de experimentos llevados a cabo en un
intercambiador de calor específico.
Para una sola fase, el coeficiente convectivo se determina utilizando el número de Nusselt (𝑁𝑁𝑁𝑁), en este
caso para calcular en número de Nusselt utilizamos las siguientes fórmulas:
Para flujos laminares: Para flujos turbulentos:
𝑁𝑁𝑁𝑁𝜑𝜑 = 0.05𝑅𝑅𝑅𝑅𝜑𝜑 0.67 𝑃𝑃𝑃𝑃𝜑𝜑 0.3 (11) 𝑁𝑁𝑁𝑁𝜑𝜑 = 0.023𝑅𝑅𝑅𝑅𝜑𝜑 0.577 𝑃𝑃𝑃𝑃𝜑𝜑 0.4 (12)

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
3296
 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

Donde 𝑅𝑅𝑅𝑅 se refiere al número de Reynolds el cual determina si un flujo es laminar o turbuelento y 𝑃𝑃𝑃𝑃
es el número de Prandtl.
𝐶𝐶𝐶𝐶𝐶𝐶 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌
𝑃𝑃𝑃𝑃 = (13) 𝑅𝑅𝑅𝑅 = (14)
𝜆𝜆 𝜇𝜇
El coeficiente convectivo (ℎ𝜑𝜑 ) puede ser calculado para cada lado del intercambiador de calor usando:

𝑁𝑁𝑁𝑁𝜑𝜑
ℎ𝜑𝜑 = (11)
𝐷𝐷/𝜆𝜆𝜑𝜑

Resultados
Para desarrollar la comparación de los modelos bajo las mismas condiciones de operación se emplearon
datos experimentales tomados de una planta piloto de un intercambiador de calor los cuales han sido
reportados en [9] y [10].
80
Tcsal est Tfsal est Tcsal real Tfsal real 80 Tcsal est Tfsal est Tcsal real Tfsal real
70
70
Temperatura (°C)

Temperatura (°C)
60
60

50
50

40
40

30
0 2000 4000 6000 8000 10000 12000 14000 16000 30
Tiempo (s)
0 2000 4000 6000 8000 10000 12000 14000 16000
Tiempo (s)
Figura 2. Comparación entre las temperaturas de salida
Figura 3. Comparación entre las temperaturas de salida
estimadas por el modelo de parámetros agrupados y los
estimadas por el modelo de parámetros agrupados y los
datos obtenidos experimentalmente.
datos obtenidos experimentalmente.

Las características de operación son las siguientes: el flujo de agua caliente se mantiene constante en
1.667𝑥𝑥10−6 𝑚𝑚3 /𝑠𝑠 y el flujo de agua fría se varía en un intervalo de 6.667x10-6 m3/s y 8.922x10-7 m3/s.
La temperatura inicial del agua fría se varía en un intervalo entre 26°C y 29°C mientras que la
temperatura inicial del agua caliente se varía en un intervalo entre 49°C y 83°C.
En la Figura 2 se presenta la comparación entre el modelo de parámetros agrupados utilizando la
ecuación (7) (para el cálculo del coeficiente de transferencia de calor), y los datos experimentales del
intercambiador de calor. Para la temperatura de agua fría el error de estimación del modelo fue de 1.87 y
para la corriente de agua caliente de 0.63. Este error fue calculado considerando el segmento de datos
entre 0 y 6000, donde se cree hubo un sesgo debido a un error de medición ocasionado por bajo flujo en
la corriente de agua fría. Sin considerar los datos del 0 al 6000, para la temperatura de agua fría el error
de estimación del modelo fue de 0.3914y para la corriente de agua caliente de 0.4672.
En la Figura 3 se presenta la comparación entre el modelo de parámetros agrupados y los datos
experimentales obtenidos del intercambiador de calor, utilizando la ecuación (8) para el cálculo de los
coeficientes globales de transferencia de calor calculados para cada corriente. Para la temperatura de
agua fría el error de estimación del modelo fue de 1.62 y para la corriente de agua caliente de 0.53. Este
error fue calculado considerando el segmento de datos entre 0 y 6000, donde se cree hubo un sesgo
debido a un error de medición ocasionado por bajo flujo en la corriente de agua fría. Sin considerar los
datos del 0 al 6000, para la temperatura de agua fría el error de estimación del modelo fue de 0.0086 y
para la corriente de agua caliente de 0.3641.

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
3297
 Memorias del XXXVI Encuentro Nacional de la AMIDIQ
5 al 8 de Mayo de 2015, Cancún, Quintana Roo, México

En la Figura 4 se presenta la comparación entre

80
Tfsal est Tcsal est Tcsal real Tfsal real
el modelo de parámetros distribuidos utilizando
la ecuación (15) para calcular los coeficientes
70
convectivos (hi, ho) para cada una de las
Temperatura (°C)

60
corrientes. Para la temperatura de agua fría el
error de estimación del modelo fue de 1.6236 y
50 para la corriente de agua caliente de 1.12. Este
error fue calculado considerando el segmento de
40
datos entre 0 y 6000, donde se cree hubo un
sesgo debido a un error de medición ocasionado
30
0 2000 4000 6000 8000
Tiempo (s)
10000 12000 14000 16000 por bajo flujo en la corriente de agua fría. Sin
Figura 4. Comparación entre las temperaturas de salida considerar los datos del 0 al 6000, para la
estimadas por el modelo de parámetros distribuidos y los temperatura de agua fría el error de estimación
datos obtenidos experimentalmente. del modelo fue de 0.84 y para la corriente de
agua caliente de 0.72.

Conclusión
De acuerdo con los resultados obtenidos de la comparación de los modelos y de los datos
experimentales, el modelo de parámetros agrupados (modelo de EDO) cuyo modelo termodinámico
emplea el coeficiente global de transferencia de calor para cada uno de los fluidos fue el que mostró
mejor desempeño en la estimación de las temperaturas. Por lo que podemos considerar este modelo
adecuado para el diseño de controladores u observadores de estado.

Bibliografía
1. Astorga-Zaragoza, C.M., A. Zavala-Río, V. M. Alvarado, R. M. Méndez, y J. Reyes-Reyes. Performance monitoring of
heat exchangers via adaptive observers. Measurement 40 (Mayo 2007): 392-405.
2. Hangos, K., J. Bokor, y G. Szederkényi. Analysis and control of nonlinear process systems. Budapest, Hungary: Springer,
2004.
3. Steiner, M. Low order dynamic models for heat exchangers.» Proceedings of the international Sympos on District Heat
Simulation, 1989.
4. Escobar, R. F., C. M. Astorga-Zaragoza, A. C. Telléz-Anguiano, D. Juárez-Romero, J. A. Hernández, y G. V. Guerrero-
Ramírez. Sensor fault detection and isolation via high-gain observers: Application to a double-pipe heat exchanger. ISA
Transactions, nº 50 (2011): 480-486.
5. Escobar, R., C. Astorga Zaragoza, D. Juárez Romero, D. Hernández, y M. Adam Medina. Estimación en Línea del
Coeficiente de Transferencia de Calor en un Intercambiador de Calor. Asociación Mexicana de Control Automático., 2009
6. Bagui, F., M. A. Abdelghani Idrissi, y H. Chafouk. Heat exchanger Kalman filtering with process dynamic
acknowledgement. Computers & Chemical Engineering., nº 28 (2004): 1465-1473.
7. López Zapata, Betty Y., Manuel Adam Medina, R. F. Escobar, y C. M. Astorga Zaragoza. Diagnóstico de fallas en
intercambiadores de calor: enfoque Multi-Modelos. Asociación Mexicana de Control Automático., 2009: 1-6.
8. Holman, J. P. Heat Transfer. Primera Edición. Traducido por Raul Valenzuela M. McGraw Hill, 1999.
9. Garcia Morales, J. «Diagnóstico de fallas para componentes de un intercambiador de calor empleando el enfoque de
observadores en modos deslizantes super-twisting. CENIDET, Cuernavaca, Morelos, s.f.
10. López Zapata, Betty Yolanda. Desarrollo de un sistema de detección de fallas utilizando un enfoque multi-modelos LTI
aplicado a un intercambiador de calor.» CENIDET, Cuernavaca, Morelos, 2009.

© 2015 Academia Mexicana de Investigación y Docencia en Ingeniería Química

ISBN 978-607-95593-3-5
View publication stats
3298