ANALISIS DE TRANSFERENCIA DE CALOR MEDIANTE

ELEMENTOS FINITOS EN UN HORNO DE COMBUSTIÓN

PARA LA PRODUCCIÓN DE PILONCILLO
ANALYSIS OF HEAT TRANSFER BY MEANS OF FINITE ELEMENTS IN
A COMBUSTION OVEN FOR PILONCILLO PRODUCTION
Adriana Getsemaní Mora Canseco*, Anselmo Osorio Mirón, Esthefanía de la Cruz Reyes

Universidad Veracruzana, México, getzady@[Link]; Prolongación Oriente 6 No. 1009, Col. Rafael Alvarado,
Orizaba, Veracruz, México.

RESUMEN: La mayoría de los equipos utilizados en el proceso de producción de piloncillo carecen de

requerimientos necesarios para proporcionar las características deseables del producto, lo que dificulta
su posicionamiento en el mercado, provocando baja demanda del producto. En el proceso de piloncillo, el
mecanismo principal de transferencia de calor se lleva a cabo en las hornillas paneleras, en donde
existen pérdidas por chimeneas, paredes y una ineficiencia térmica debido al material de construcción y a
un diseño no apropiado. En este trabajo, se realiza un análisis de transferencia de calor mediante
elementos finitos en un horno de combustión para la producción de piloncillo, con el objeto mejorar la
fase de combustión e incrementar la eficiencia térmica. La metodología que se llevó a cabo, es la
siguiente: primero se diseñó la geometría en el software SolidWorks 2010. Enseguida se realizó un
análisis de transferencia de calor al horno mediante simulación numérica utilizando el software de
Dinámica de Fluidos Computacionales (CFD) Comsol Multiphysics 3.5, para ello se definió el modelo y la
dimensión para la simulación, se ingresó el elemento del método numérico de resolución, se importó la
geometría, se procedió a definir las variables involucradas en el proceso, posteriormente se definió el tipo
de malla, finalmente se seleccionaron los gráficos para el postprocesado. Con la simulación se obtuvieron
resultados gráficos y numéricos de las temperaturas alcanzadas en la fase de combustión así como el
flujo de calor convectivo, que permiten mejorar el diseño de la fuente de energía en el proceso.

Palabras Clave: Piloncillo, Horno de combustión, Transferencia de calor, Elementos finitos.

ABSTRACT: Most of the equipment used in the brown sugar (piloncillo) production process no
requirements for providing the desirable characteristics of the product, making it difficult to market
position, causing low product demand. In the process of brown sugar, the primary mechanism of heat
transfer takes place in ovens, where there are losses by smokestacks, walls and thermal inefficiency
because of the material of construction and design misused. In this sense, an analysis of heat transfer
through finite elements in a combustion furnace for the production of brown sugar is realized, with the
object of improving the combustion phase and increase thermal efficiency. The methodology for the
development of the paper is as follows: first, geometry was designed in SolidWorks 2010 software. After,
an analysis of heat transfer to the furnace was realized by a numerical simulation software like Comsol
Multiphysics 3.5 that using the Computational Fluid Dynamics (CFD); here it defined the model and
dimension for simulation, was entered the element of resolution numerical method, the geometry was
imported, then we defined the variables involved in the process, and defined the mesh type, finally the

“XI Simposio Internacional de Automatización”

 Mora, A. G.; Osorio A., De la Cruz, E. |”Análisis de transferencia de calor mediante elementos finitos en un horno de combustión
para la producción de piloncillo”

graphics for postprocessing were selected. The graphical and numerical results of the temperatures
reached in the combustion phase as well as convective heat flux, obtained by numerical simulation can
improved the design of the power source in the process.

KeyWords: Piloncillo, Furnace combustion, Heat transfer, Finite elements.

1. INTRODUCCION canal de evacuación y después suben por la

chimenea, de esta manera se tiene un tiro vertical
(TV) y un tiro horizontal (TH).
El piloncillo es considerado un producto alimenticio,
el cual es obtenido de la evaporación y Para calcular la altura de la chimenea (H) se aplica
cristalización de los jugos de la caña de azúcar; se la fórmula:
distingue por tratarse de un producto natural, sin
refinación, con mínimo procesamiento y con la
mayor cantidad de constituyentes originales del H 3 ∗ TV (1)
jugo de la caña de azúcar. Comúnmente, el
piloncillo se obtiene en pequeñas factorías rurales
denominadas “trapiches”, la evaporación del jugo Las dimensiones y características del hogar
de caña se lleva a cabo en una paila soportada en dependerán del tipo y capacidad del horno y del
un horno de combustión, alimentado por bagazo de recorrido que tenga el fuego en el interior de la
caña de azúcar previamente secado.[1] cámara; sin embargo, por grande que sea el horno,
el hogar siempre deberá tener medidas mínimas.
La mejora del proceso de producción puede incluir
cada una de las etapas, desde la molienda hasta la Para el diseño la cámara de combustión se
concentración, enfriamiento y moldeo, por lo que al tomaron en cuenta los siguientes parámetros:
observar detenidamente el proceso tradicional de
producción de piloncillo, se puede encontrar como  Humedad y tamaño del bagazo
inconveniente la existencia de pérdidas de energía  Temperaturas requeridas
en el horno de combustión originadas por el  Capacidad de la hornilla
material de construcción así como por el mal  Aspectos económicos
dimensionamiento del mismo.
Mientras que para el ducto se consideraron dos
Este trabajo hace referencia a una parte del secciones: en la primera, denominada zona caliente
proceso tecnificado en la elaboración de piloncillo a debido a su proximidad a la cámara, la temperatura
partir del jugo de la caña de azúcar, conservando de los gases de combustión es elevada; en la
las características del proceso tradicional con la segunda, la temperatura de los gases de
finalidad de mantener la identidad del producto, y combustión es relativamente baja y se denomina
considerarlo como un alimento natural zona fría.
mínimamente procesado. Para el análisis térmico
del horno de combustión, se ha tomando en cuenta
el combustible empleado y el material de A partir de estos datos, se realizó la geometría en el
construcción. software SolidWorks 2010.

2. CONTENIDO Para el análisis de transferencia de calor mediante

En este apartado se definieron las condiciones de elementos finitos, empleando el software Comsol
diseño empleadas para el horno de combustión, se Multiphysics versión 3.5, se consideraron los
describe el horno como, un intercambiador de calor resultados de los balances de materia en el área de
en el que el fluido de proceso se calienta por combustión y el balance de energético del proceso
conducción, con la energía procedente de una de producción.
llama de combustión, y por convección desde los
gases calientes de la misma combustión. En la figura 1, se muestra el diagrama utilizado para
el balance de materia en el área de combustión en
El diseño del horno analizado es el de tipo cruzado, el proceso de producción de piloncillo. [2]
en este tipo de horno los gases primero suben,
bajan, recorren el horno horizontalmente hacia el

“XI Simposio Internacional de Automatización”

 Mora, A. G.; Osorio A., De la Cruz, E. |”Análisis de transferencia de calor mediante elementos finitos en un horno de combustión para la
producción de piloncillo”

Gases de
Bagazo
Sistema o combustión
proceso
1 productivo 3

Humedad
del Residuo
bagazo s
Aire 2
(T,H

Figura 1 Balance de materia en el área de

combustión

En la figura 2, se esquematiza el sistema que se

utilizó para el balance energético del proceso de
producción de piloncillo.

Bagazo E5
Figura 3. Navegador de modelos Comsol Multiphysics
Gases de 3.5
Bagazo E1
combustión E3 Se importó la geometría en 3D generada en el
software SolidWorks 2010 y se procedió a definir
Hornilla Residuos E4 las variables involucradas en el proceso. La
Aire E2 ecuación 2 modela la distribución de la temperatura
en el horno de combustión para la transferencia
Pérdidas por las conductiva y convectiva de calor.
paredes
E6
Figura 2 Balance energético en el área de combustión
∙ T Q ρC u ∙ T (2)

Donde T es temperatura,  es la densidad (Kg/ m3),

C es la capacidad de térmica a presión constante
2.1 Materiales y métodos
(J/kgk), Q fuente de calor (W/m3) y u el campo de
velocidad (m/s). [3]
Para el análisis de transferencia de calor en el
horno de combustión, se inició definiendo el
modelo, seleccionando el modelo físico y el tipo de Posteriormente, se definieron las condiciones de
dimensión para la simulación, de igual manera se frontera: temperatura teórica de entrada y salida de
ingresaron las variables dependientes y el elemento los gases, temperatura ambiente, conductividad
del método numérico de resolución, a través del térmica del gas (aire), densidad y coeficiente de
Navegador de Modelos como se muestra en la transferencia de calor. El flujo por las paredes se
Figura 3. despreció (no sale calor por éstas en el modelo
computacional). Los valores teóricos del flujo de
gases se muestran en la Tabla I.

“XI Simposio Internacional de Automatización”

 Mora, A. G.; Osorio A., De la Cruz, E. |”Análisis de transferencia de calor mediante elementos finitos en un horno de combustión para la
producción de piloncillo”

Tabla I: Valores teóricos del flujo de gases

Variable térmica Unidades Valor

Temperatura en la cámara de
K 1300
combustión en la entrada
Temperatura en la chimenea
K 401
en la salida
Temperatura ambiente K 300
-3
Densidad gas @1300k Kg m 0.2679
-3
Densidad gas @744k Kg m 0.4643

Calor especifico gases de -1 -1

J kg m K 1189
combustión Figura 4. Mallado libre de horno de combustión
-2 -1
Coeficiente de convección Wm K 24.4
1500x1 En la Figura 5, se muestra la solución de geometría
Flujo de calor de entrada W 3
0 de la hornilla, el análisis se lleva a cabo mediante el
Presión de entrada Pa 0 elemento Lagrange-Cuadrático en estado
-1 estacionario, con un resolvedor de sistema lineal
Velocidad de entrada del aire ms 14.1
Directo (PARDISO).
Los valores teóricos de los gases mostrados en la Tabla I
indican las variables aplicadas en la simulación numérica.

Se definido el tipo de mallado libre, por el método

de triángulo.
Para los parámetros del resolvedor se estableció un
análisis estacionario, con tolerancia relativa de
1.0e-6, con un número máximo de 25 iteraciones.

Para el postprocesado se seleccionaron los gráficos

subdominio, arista y línea de flujo. En los gráficos
subdominio y arista se estableció la Temperatura
como cantidad predefinida, indicando además en el
grafico subdominio una selección de elementos Figura 5. Simulación de horno de combustión,
aleatoria y un valor de 1 para la fracción de indicando las temperaturas
elementos. Las líneas de flujo se predefinieron por
el flujo de calor convectivo. La Figura 6, representa la solución del
postprocesado, se aprecia el comportamiento de
las líneas de flujo de calor convectivo alcanzando
una temperatura máxima de 1363.938 K
2.2 Resultados equivalente a 1090.788 °C, el área de combustión
de gases presenta las temperaturas más elevadas,
En la Figura 4, se muestra el resultado del mallado mientras el área de la chimenea muestra las zonas
de la hornilla mediante elementos finitos, con un más frías.
tamaño normal y con una resolución geométrica de
20. Cabe señalar que no se tomó en cuenta ningún
incremento de tiempo ni de tamaño. El mallado de
la geometría consta de 3227 elementos.

“XI Simposio Internacional de Automatización”

 Mora, A. G.; Osorio A., De la Cruz, E. |”Análisis de transferencia de calor mediante elementos finitos en un horno de combustión para la
producción de piloncillo”

cámara de combustión y manifestando una mejor

transferencia de calor hacia la paila de evaporación.

Figura 6. Comportamiento de las líneas de flujo de

calor convectivo en horno de combustión Figura 8 Presentación de horno de combustión en 3D

La Figura 7 muestra como el diseño creado permite Obtenidos los datos de diseño de un horno de
concentrar una temperatura máxima de 1044.786 combustión para la producción de piloncillo, se
°C en la zona de retorno de los gases de pone de manifiesto que el material de construcción
combustión, en esta zona se encuentra la última del mismo, influye en la eficiencia de distribución de
sección de la paila de evaporación, y en donde se temperatura. En la tabla II se condensan los datos
requieren temperaturas elevadas para que el jugo de diseño para un horno de combustión para la
de caña alcance temperaturas superiores a los producción de piloncillo, a partir de caña de azúcar.
100°C y se formen las mieles.
Datos generales
Posición Horizontal
Longitud (m) 3.2 m
Ancho 0.8 m
Altura 1.10 m
3 3
Capacidad (m ) 0.54 m
Composición
A Hornilla
B Ducto de humos
C Chimenea
Datos de diseño Hornilla
Material Ladrillo refractario

Figura 7. Zona más caliente en el horno de Consumo energético Bagazo

combustión
Temperatura de diseño 1650 °C
La figura 8 muestra el prototipo final, diseñado en el (°C)
software SolidWorks 2010, a partir de datos Temperatura de 1027 °C
recabados en visitas realizadas a diversos trapiches trabajo (°C)
de la zona centro del estado de Veracruz y del Conductividad térmica 0.025 W/(m*K)
análisis numérico efectuado al horno de
3
combustión. Densidad 1 Kg/m

Capacidad térmica a 1189 J/(kg*K)

El horno de combustión que se presenta es de tiro
Fuente de calor 1500e3 W/ m3
cruzado, en donde los gases tienen una mejor
circulación, manteniéndose más tiempo en la

“XI Simposio Internacional de Automatización”

 Mora, A. G.; Osorio A., De la Cruz, E. |”Análisis de transferencia de calor mediante elementos finitos en un horno de combustión para la
producción de piloncillo”

Datos de diseño Ducto de humos

5. REFERENCIAS BIBLIOGRÁFICAS
Material Ladrillo refractario

Temperatura de diseño (°C) 1650 °C

1.[Link]
Temperatura de trabajo (°C) 950 °C E/COVECAINICIO/IMAGENES/ARCHIVOSPDF/AR
CHIVOSDIFUSION/MONOGRAF%CDA%20DE%2
Conductividad térmica 0.025 W/(m*K) [Link]

Densidad 1 Kg/m
3 2 “Propuesta de mejora en la eficiencia energética y
de producción de piloncillo en trapiches de la zona
Capacidad térmica a 1189 J/(kg*K)
presión constante centro del estado de Veracruz” Tesis Profesional,
3 Universidad Veracruzana Trujillo V. A., Guerra S.Y.
Fuente de calor 1500e3 W/ m
Orizaba, Ver. 2011.
Datos de diseño Chimenea
3 Torres R. y Grau J. Introducción a la mecánica
Material Ladrillo refractario de fluidos y transferencia de calor con Comsol
Mutiphysics. Editorial Addlink media 2007.
Altura 3.05 m

Se obtienen las especificaciones y condiciones de operación del

6. SÍNTESIS CURRICULARES DE LOS
horno, condensadas en la tabla II. AUTORES
3. CONCLUSIONES Adriana Getsemaní Mora Canseco
8 de Septiembre de 1989. Orizaba, Veracruz
Este trabajo, permitió el análisis de transferencia de México.
calor en la fase de combustión de la producción de Estudiante de 9 Semestre de la Carrera de
piloncillo a partir de caña de azúcar, considerando Ingeniería Química, en la FCQ-UV, Orizaba, Ver.
la metodología de trabajo planteada inicialmente. Asistente de Proyecto PROMEP Convocatoria
2011,
Los resultados numéricos muestran que el horno
Estancia de Verano en la Investigación UV, 2012,
diseñado para el proceso de piloncillo, es eficiente,
FIME-Xalapa, Ver.
dado que presenta distribuciones de variables o
Asistente de Proyecto FOMIX CONACYT-Veracruz
parámetros mas uniformes en cuanto a
Convocatoria 2009-3, Etapa III.
temperaturas que los diseños tradicionales
Asistencia a Cursos de Formación y Capacitación
empleados en la mayoría de los trapiches de la
en Simulación y Automatización de Procesos, 2011-
zona centro del estado de Veracruz, tomando en
2012.
cuenta que en la mayoría de estos trapiches, se
construyen con ladrillo rojo, que soporta menor
DATOS DE CONTACTO
temperatura comparado con el ladrillo refractario, el
Dirección:
cual se empleó en el presente análisis.
Facultad de ciencias químicas.
4. AGRADECIMIENTOS (Orizaba).
Se agradece al FOMIX CONACYT-VERACRUZ, Universidad veracruzana.
MEXICO, por el apoyo brindado en la realización de Prol. Oriente 6 no. 1009.
este trabajo (PROYECTO 128573). Col. Rafael alvarado.
C. P. 94340 orizaba, ver., méxico.
Agradecimiento especial al Ing. Sarahí pacheco Teléfono:
Espinoza, Ing. Gerardo Evaristo González Landero, 045 271 1394504
M. C. Ignacio Sánchez Bazán y M. C. Luis Alberto Correo electrónico:
Sánchez Bazán, por su apoyo y contribución a la getzady@[Link]
elaboración de este trabajo.

“XI Simposio Internacional de Automatización”