Modelado de La Transferencia de Calor en El Tratam
Modelado de La Transferencia de Calor en El Tratam
Modelado de La Transferencia de Calor en El Tratam
net/publication/311166068
CITATIONS READS
0 1,275
3 authors:
Luis Medina-Torres
Universidad Nacional Autónoma de México
118 PUBLICATIONS 1,548 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Obtaining and evaluation of liquid crystalline and other nanostructured systems applicated in nanomedicine. Obtención y evaluación de cristales líquidos y otras
nanoestructuras de aplicación en nanomedicina. View project
All content following this page was uploaded by José Javier Valencia-López on 03 December 2016.
Recibido Dic. 11, 2015; Aceptado Feb. 16, 2016; Versión final Abr. 16, 2016, Publicado Dic. 2016
Resumen
Palabras clave: dinámica de fluidos computacional; transferencia de calor; producto enlatado; leche
condensada
A mathematical model and its numerical solution in Computational Fluid Dynamics software (CFD, for its
acronym in English) for the determination of heat transfer on a foodstuff canned, with non-Newtonian
rheological properties. The heat transfer inside the can is analyzed in order to obtain temperature profiles
that help to determine the optimum conditions of heat transfer in the process of sterilization of canned
products. In addition, the influence of heat transfer on the kinetics of destruction of pathogenic
microorganisms such as the spore-producing bacteria Clostridium botulinum, is analyzed.
Keywords: computational fluid dynamics; heat transfer; canned product; condensed milk
INTRODUCCIÓN
En la actualidad los tratamientos térmicos siguen siendo la manera más utilizada de preservar los alimentos,
desde que fue propuesta por Nicolás Appert en 1810, independientemente del desarrollo de tecnologías
alternativas denominadas tratamientos no térmicos (campo de pulsos eléctricos, alta presión hidrostática,
UV, entre otros) (Farid y Ghani, 2002). Existen dos métodos diferentes de procesamiento térmicos
convencionales aplicados en la industria alimentaria: el procesamiento aséptico, en donde el producto
alimenticio es esterilizado antes de ser envasado y aquel en el que el producto alimenticio se envasa y
después se esteriliza (Erdogdu et al., 2010). Siendo uno de los procesos industriales de mayor consumo
energético e inversión de capital representando entre el 10 y 15 % de las manufacturas totales (Dincer,
1998), razones por las que continúan abiertas líneas de investigación enfocadas a la mejora continua de los
procesos térmicos con la finalidad de reducir costos de producción, optimizar el uso de recursos energéticos
y maximizar la retención de nutrientes. El objetivo final es procesar alimentos seguros, de alta calidad y a un
precio que el consumidor esté dispuesto a erogar (Ghani, 2006).
El diseño de un proceso térmico, requiere la identificación del microorganismo a inactivar y que guarda
relación con el producto alimenticio. Para alimentos de baja acidez (pH > 4.6) como es el caso de la leche
condensada, se da una atención especial al Clostridium botulinum, microorganismo que se toma como base
para el diseño de procesos de esterilización, ya que una bacteria anaerobia formadora de esporas
resistentes al calor y productor de una toxina letal para el hombre aún en concentraciones de 10 ng/kg de
peso corporal (Alvarado et al., 2009). El tratamiento térmico de alimentos no es selectiva por lo que son
inherentes al proceso pérdidas del valor nutricional y sensorial del alimento. El diseño efectivo de un
proceso térmico, requiere del conocimiento de cinéticas de destrucción de microorganismos, enzimas y
nutrientes asociados a la calidad del alimento, además de perfiles de temperatura de la zona de
calentamiento más lento (ZCL) del envase (Alvarado et al., 2009).
Debido a la cinética de primer orden que sigue la inactivación microbiana es prácticamente imposible la
obtención de un producto estéril. La destrucción completa de formas viables, representaría la destrucción
del alimento (Valle y Merson, 1981). Los tratamientos térmicos, además de inactivar microorganismos,
también destruyen algunas vitaminas termolábiles como la tiamina y promueven la oxidación de lípidos
ocasionando un sabor desagradable al producto alimenticio. Los efectos del calor sobre proteínas y lípidos
han sido descritos con detalle por Mauron en 1982, Witting y Dimick en 1982 (Fellows, 2000). Motivo por el
cual la industria alimenticia se ha dado a la tarea de investigar mejoras en la tecnología del procesamiento
térmico como la agitación (Singh y Ramaswamy, 2015) o el procesamiento aséptico con el fin de
proporcionar una buena calidad del alimento sin comprometer la seguridad de la salud pública (Singh et al.,
2015).
Los mecanismos más importantes de transmisión de calor en el procesamiento térmico de alimentos son:
conducción, convección natural y/o forzada. Otro factor importante son las propiedades termofísicas de los
alimentos en los procesos térmicos ya que éstas determinan la velocidad de transferencia de calor en el
interior del producto. Así que en los procesos térmicos de calentamiento y enfriamiento se ha establecido
que las propiedades primarias comprenden: la conductividad térmica, el calor específico y la densidad del
alimento (Ashrae, 1977). Si la transferencia de calor es controlada por conducción, la ZCL se localizará en el
centro geométrico del envase, teniendo poco efecto sobre su ubicación, la dependencia de las propiedades
termodinámicas con la temperatura (Lucatero, 2000). La convección natural por su parte, provoca que la
ZCL migre hacia el fondo del envase durante el proceso térmico, desplazándose hacia el eje axial y luego
hacia el centro de la lata (Jiménez et al., 2005).
Los modelos matemáticos proporcionan una de las herramientas básicas para describir procesos físicos y
predecir su comportamiento en condiciones variadas. A través de su uso es posible el diseño y mejorar el
entendimiento de proceso en el área de ingeniería alimentaria. Una simulación del proceso utilizando un
modelo matemático es esencial para la optimización y el control de procesos. (Tirado et al., 2007). Asi que
debido a la dificultad de ubicar la ZCL para alimentos enlatados y los errores asociados al uso de
termopares, han generado interés en el entendimiento de los fenómenos de transporte asociados a dicho
proceso, lo que implica la solución numérica de las ecuaciones de conservación de momentum, energía y
masa de manera acoplada, convirtiéndolo en un problema de convergencia difícil. La solución no hubiera
sido posible, sin la evolución de los sistemas de cómputo de alta velocidad, que permitieron el desarrollo de
software de Dinámica de Fluidos Computacional (CFD). La simulación numérica provee información más
detallada y precisa sobre la distribución de temperatura y patrones de flujo dentro del envase. Por lo
anterior, es de gran importancia plantear un modelo que describa la transferencia de conservación de
cantidad de movimiento, masa y energía.
Varios modelos matemáticos predictivos han sido empleados para la estimación de la distribución de
temperatura para un alimento enlatado. En cuanto a los alimentos líquidos enlatados muchos investigadores
han evaluado experimentalmente el coeficiente de transferencia de calor durante el tratamiento térmico y
han registrado el efecto de diferentes procesos y variables de productos dependiendo el coeficiente de
transferencia de calor (Singh et al., 2015).
METODOLOGÍA
Se plantea el problema en función de las leyes de conservación de cantidad de movimiento y energía que
constituyen la base de estudio del fenómeno convectivo, las cuales describen las variaciones de velocidad y
temperatura con respecto al tiempo y posición en el sistema (Moraga et al., 2011). El problema planteado se
resuelve tomando en cuenta una lata individual de radio R y altura L (ver Figura 1), la cual se encuentra a
una temperatura inicial To, que contiene un alimento constituido por un fluido homogéneo. Se resuelve un
caso axisimétrico que permite modelar un problema originalmente en 3-D como uno en 2-D, debido a la
geometría cilíndrica del envase a modelar. Fluido en régimen laminar, no-Newtoniano, isotrópico. El modelo
planteado se resuelve y se valida con resultados publicados por (Moraga et al., 2011) una vez validado el
modelo, se resuelve el mismo caso de estudio pero con los datos fisicoquímicos y reológicos que describen
a la leche condensada.
Se dibuja una geometría rectangular con dimensiones de 0.12 m de altura y 0.05 m de radio, se construyó
una malla en forma tetraédrica
Ecuación de continuidad
r r z (1)
0
r r z
r r p 2 1 r 2 r (2)
r z r 2r
t r r z r r r z 2
z z p 2 1 z 2 z (3)
r z z 2z g ref 1 T - Tref
t r r z r r r z2
Ecuación de Energía
T T T 2T 1 T 2T (4)
Cp z r k 2
t z r r r r z 2
dv y
n (5)
xy m
n
m
dx
Condición inicial
t 0s T 25C
Condiciones de Frontera
T v r
En r 0; 0 0
r r
En r R; v r 0 v z 0; T 121C
En z 0; v r 0 v z 0; T 121C
En z L; v r 0 v z 0; T 121C
En la Figura 2 se muestra el modelo en 2-D indicando las condiciones de frontera para la geometría
construida.
Para la determinación de la inactivación microbiana (U) se utiliza la ecuación propuesta por Bhuvaneswari y
Anandharamakrishnan. (Bhuvaneswari et al., 2014)
T Tref (6)
U 10 ZT
RESULTADOS Y DISCUSIÓN
Se presentan los resultados en dos etapas: los resultados que validan el modelo y luego una parte donde se
muestra un caso de estudio, utilizando el modelo ya validado.
Antes de resolver las ecuaciones para el modelo de la leche condensada, se realiza una validación del
modelo comparándolo con los resultados obtenidos en un artículo previo publicado por (Moraga et al.,
2011), en el cual se presenta un modelo sobre la transferencia de calor por convección natural durante la
esterilización de un fluido no-Newtoniano enlatado, utilizando como fluido modelo Carboximetil Celulosa
(CMC), los parámetros fisicoquímicos utilizados se muestran en la Tabla 1.
La simulación y los resultados obtenidos en la validación del modelo se presentan a continuación, donde se
hace una malla fina en forma tetraédrica con 614 elementos a computar en un tiempo de 1500 s (Moraga N.
et al., 2011) (Ver Figura 3).
Fig. 3 Resultados obtenidos en la validación del modelo con una malla de 614 y 7892 elementos
En la Figura 3 a) se puede observar que los perfiles de temperatura obtenidos con una malla de 614
elementos son muy similares sin embargo se puede apreciar que existe un perfil en el cual se presenta una
perturbación debido al número de elementos resueltos en la malla, es por ello que se hace una nueva
simulación tomando los mismos datos de la Tabla 1 pero tomando una malla con 7892 elementos a
resolver, obteniendo los resultados que se muestran en la Figura 3 b). En la Figura 3 b) se puede observar
que al resolver el modelo una malla extremadamente fina con 7892 elementos con un tiempo de
computación de 1000 s la perturbación que se muestra en la Figura 3 a) en los perfiles de temperatura se
corrigió de tal manera que los resultados obtenidos en este trabajo son iguales a los reportados (Moraga et
al., 2011).
Es importante señalar que el tipo de malla y el número de nodos a resolver es importante en el modelo, ya
que a mayor número de elementos o nodos, los resultados obtenidos serán más confiables aunque el
tiempo de resolución será mayor. Una vez validados los perfiles de temperatura y debido a que está provoca
un cambio en la densidad y por consiguiente un movimiento del fluido. A continuación se presenta la
comparación de las líneas de corriente que describen el movimiento del fluido dentro de la lata para ambos
modelos (Ver Fig. 4). En la Figura. 4 se puede observar el movimiento del fluido durante el periodo del
tiempo en alcanzar la temperatura de esterilización, donde las líneas de corriente son de tipo oscilatorias
posicionándose casi al fondo de la lata debido a las fuerzas de flotación involucradas dentro de la lata.
Una vez validado el modelo se puede tomar un líquido alimenticio con comportamiento no- Newtoniano para
evaluar el tiempo y la temperatura que ayuden a determinar las condiciones óptimas en el proceso de
esterilización de productos enlatados.La leche condensada es un tipo de leche evaporada por ser el
producto obtenido de la evaporación de agua de la leche entera o leche total o parcialmente descremada
con o sin la adición de azúcar. Aunque el producto se suministra a granel como una materia prima para la
industria alimentaria la cual se prepara principalmente para el consumo directo donde se empaca en
cantidades relativamente pequeñas en latas o en botellas selladas, las cuales se someterán a un proceso
térmico. Los parámetros fisicoquímicos y reológicos utilizados para la evaluación del modelo caso de
estudio se describen en la Tabla 2.
La simulación y los resultados obtenidos en la validación del modelo se presentan a continuación, donde se
realiza un malleo fino en forma tetraédrica con 7892 elementos a computar en un tiempo de 1000 s (Ver
Fig. 5). El periodo de esterilización debe variar de 15 a 20 minutos a temperaturas de 115 -116° C o
superiores (Van Den Berg, 2001). En la Figura 5 se muestra la formación de los perfiles de temperatura
dentro de una lata a tomados en diferente tiempo del proceso de esterilización, resultado de la solución
numérica del problema planteado.
Fig. 5. Capturas de imágenes sobre los perfiles de temperatura dentro de la lata tomadas en diferente tiempo de
esterilización, a) 0 min, b) 1 min, c) 3 min, d) 5 min, e) 10 min, f) 15 min
En la Figura 5 se puede observar la formación de los perfiles de temperatura, donde a los 3 minutos del
proceso de esterilización se comienzan a formar perfiles tipo herradura característicos de los procesos de
tratamiento térmico de productos enlatados (Galotto et al., 2011), cuando el tiempo de esterilización es de
15 minutos se puede observar la clara formación de los perfiles de temperatura de tipo herradura los cuales
se obtienen cuando la lata alcanza el estado estacionario (~121°C), donde la zona de calentamiento más
lenta está indicada en color azul la cual se posiciona al fondo de la misma debido a las propiedades del
fluido y las fuerzas de flotación y gravimétricas que están actuando en el proceso de convección. En el
artículo de (Erdogdu et al., 2010) analizan la transferencia de calor de una lata que contiene agua (fluido
Newtoniano) en el cual se menciona que el del fluido con mayor temperatura comienza a desplazarse hacia
arriba hasta que es desviado por la superficie de la lata generando un movimiento en la dirección radial y el
fluido con menor temperatura es desplazado hacia la parte inferior ocasionando una recirculación del flujo y
aumentando la transferencia de calor ubicando así la ZCL en el centro de la lata. En la Figura 5 se puede
observar que el efecto del movimiento del fluido descrito por (Erdogdu et al., 2010) es similar al movimiento
que se presenta en la leche condensada (fluido no-Newtoniano), ya que los perfiles de temperatura de tipo
herradura que se forman pueden ser diferentes según sean las variaciones de la densidad del fluido
ocasionadas por los gradientes de temperatura que se tienen en los diferentes puntos dentro de la lata.
El tratamiento térmico de un alimento líquido como es el caso de la leche condensada debe llevarse a cabo
después del enlatado para lograr una destrucción mínima de nutrientes y aspectos físicos como color y
textura. En este trabajo se planteó un modelo de trasferencia de calor con convección natural libre donde se
transfiere el calor de forma natural hacía todo el alimento enlatado debido a las fuerzas de flotación que se
originan en el envase tal y como se muestra en la Figura 5. Por otro lado la convección forzada se transfiere
el calor por aplicación externa de energía como la agitación giratoria o lateral de la lata con lo cual se mejora
la transferencia de calor (Singh et al., 2015).
CONCLUSIONES
Se cuenta con un modelo en 2D para la transferencia de calor y de cantidad de movimiento para un fluido
no-Newtoniano, dicho modelo resuelve las ecuaciones diferenciales parciales de manera simultánea
mediante el método de elemento finito con un software comercial de Dinámica de Fluidos Computacional,
con el cual se puede estimar: los perfiles de temperatura a través del tiempo, el movimiento del fluido
contenido en la lata ocasionado por los efectos de calentamiento convectivo y efectos viscosos, el tiempo
que tarda en llegar a la temperatura de esterilización y conocer la cinética de inactivación de
microorganismos patógenos.
Es importante conocer los perfiles de temperatura que gobiernan el proceso de esterilización dentro de la
lata ya que con estos es posible conocer la cinética de destrucción de microorganismos patógenos, aunque
cabe señalar que dicha cinética presentada en este trabajo fue calculada para una bacteria en especial
Clostridium botulinum, esta bacteria es tomada como base para diseñar procesos de esterilización de
alimentos enlatados debido a que es una bacteria termo resistente, entonces si se inactiva esta bacteria
podemos asegurar que los demás microorganismos patógenos también han sido inactivados y así poder
tener un alimento inocuo. Con la solución del modelo plateado en este trabajo es posible observar las
condiciones de operación del proceso de esterilización y así poder optimizar el proceso.
REFERENCIAS
Alvarado J., Botello E., Calderón M., Jiménez H., Martínez G., Navarrete J. Fenomenología de la
esterilización de alimentos líquidos enlatados. Revista de la Facultad Ingeniería Universidad de Antioquia:
50, 87-98 (2009)
Ashrae, Fundamentals Handbook. Thermal Properties of Foods, American Society of Heating. Refrigerating
and Air Conditioning Engineers Inc. New York, N.Y. USA. (1977)
Dincer, I. Heat transfer during heat sterilization and cooling processes of canned products. International
Communications of Heat and Mass Transfer: 34, 101-105 (1998)
Erdogdu F., Uyar R., Palazoglu T.K., Experimental comparison of natural convection and conduction heat
transfer, doi: 10.1111/j.1745-4530.2008.00309.x, Journal of Food Process Engineering 33, 85–100, (2010)
Fellows, P. J. Food Processing Technology. 2nd Edition, Woodhead Publishing Limited, England (2000)
Tirado F., Acevedo D., Puello P., Determinación Computacional del Coeficiente de Transferencia de Calor y
Propiedades Termofísicas de Alimentos, doi: 10.4067/S0718-07642014000300008, Información
Tecnológica, 25 (3), 53-58, (2014)
Ghani, A. G. A computer simulation of heating and cooling liquid food during sterilization process using
computational fluid dynamics. Association for Computing Machinery New Zealand Bulletin: 3, ISSN 1176-
9998 (2006) (En línea) http://journaldatabase.info/download/pdf/computer_simulation_heating_cooling.
(2013)
Ghani, A. G., Farid, M. M. y Chen, X. D. Numerical simulation of transient temperature and velocity profiles in
a horizontal can during sterilization using computational fluid dynamics, Journal of Food Engineering: 51,
77-83 (2002)
Hosrt Dieter Tscheuschner, Fundamentos de tecnología de los alimentos, Editorial Zaragoza, España, 300-
303, (2001)
Jiménez H., González L. M., Botello, J.E., Navarrete J. L., Estudio Numérico de la Esterilización Térmica de
Alimentos Líquidos Enlatados que Contienen Partículas, Empleando el Enfoque de Medio Poroso. Revista
Mexicana de Ingeniería Química, 4, 1-23 (2005)
Moraga N., Torres A., Galotto M. J., Guarda A., Non-newtonian canned liquid food, unsteady fluid Mechanics
and heat transfer prediction for Pasteurization and sterilization. Journal of Food Process Engineering: 34,
2000–2025 (2001)
Singh A., Ramaswamy H., Effect of product related parameters on heat-transfer rates to canned particulate
non-Newtonian fluids (CMC) during reciprocation agitation thermal processing. doi:
10.1016/j.jfoodeng.2015.05.004. Journal of Food Engineering: 165, 1-12, (2015)
Singh A., Pratap A., Ramaswamy H., Computational techniques used in heat transfer studies on canned
liquid particulate mixtures. doi: 10.1016/j.tifs.2015.02.001 Trends in Food Science & Technology: 43, 83-103,
(2015)
Valle P., Merson R. L., Cálculo del tiempo de tratamiento térmico en botes, Método General y Gráfico.
Tecnología de Alimentos 16 (3), 10-22 (1981)
Varma, M.N. y Kannan, A. CFD studies on natural convective heating of canned food in conical and
cylindrical containers. Journal of Food Engineering: 4, 1024-1036 (2006)
Van Den Berg J.C.T., Evaporated and condensed milk, Government dairy an Adviser for milk products, The
Hague, Netherlands, (2001)