Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp.
8-15, ISSN 2011-642X – E-ISSN 2389-864X
“Journal of Engineering Sciences”
Original Article DOI: [Link]
Estudio numérico del proceso de cocción de ladrillos al interior de un horno tipo
MSc. Hernando Alexander Yepes Tumay1
1 Grupo de Investigación en Nuevas Tecnologías, Sostenibilidad e Innovación (GINSTI), Departamento de Ingeniería Mecánica, Universidad Francisco de Paula Santander - Ocaña, Colombia, [Link]
org/0000-0001-9399-5673, Email: hayepest@[Link]
Revista Ingenio, vol.15, n°1, pp.
8-15, 2018, doi: [Link]
Fecha de recibido:09 de julio de 2018
Fecha aprobación: 26 de octubre de 2018
RESUMEN
Palabras claves: El aumento en la demanda energética y las emisiones contaminantes, así como la dependencia de los combustibles fósiles a
transformación energética. En este estudio se evaluó numéricamente el proceso de cocción de ladrillos al interior de un horno tipo
de ladrillos
al interior del horno evidenciaron varios campos de estudio en los que se debe profundizar para lograr una mejora energética en la
ABSTRACT
Key words: The growth of energy requirements, pollutant emissions and fossil fuel dependency at the global level require strategies to achieve
results allowed to identify some topics that require a more profound study to improve the brick burning process from an energetic
1. Introducción Una de las estrategias para reducir el consumo los
El crecimiento poblacional, los nuevos estándares de vida, combustibles fósiles por sus efectos negativos y de igual
el aumento de la economía y los avances tecnológicos han forma contribuir al cumplimiento de la demanda energética,
generado un gran incremento energético a nivel mundial. De es promover la optimización de procesos que permitan
acuerdo a la Agencia Internacional de la Energía (IEA por sus un mejor aprovechamiento de los recursos energéticos.
siglas en inglés) la demanda mundial aumentará en un 14.6% En el caso Colombiano y de América Latina este aspecto
para el año 2020 y un 48% para el 2040 en comparación a es de vital consideración dado que por ser en su mayoría
la energía requerida por la población en 2012, pasando de países en vía de desarrollo, la disponibilidad y facilidad de
549 cuatrillones de BTU a 815 en dicho periodo de tiempo
[1]. América Latina no es ajena a esta tendencia e incluso en en comparación a otros países [3]. En particular el sector
los últimos 40 años ha aumentado su consumo energético
en cerca de 2.5 veces, con un crecimiento promedio anual térmicos se presenta un rezago tecnológico que deriva en un
de 8%. En el caso colombiano de acuerdo con el Banco
Interamericano de Desarrollo para el año 2040 el incremento [4-5-22]. Dentro de estos uno de los procesos de mayor
en el consumo energético será de 110% en comparación al requerimiento energético es la cocción de cerámicos [6].
año 2013 con un aumento anual promedio de 2.8%. [2].
El uso de CFD para la optimización de hornos y sistemas
Autor para correspondencia Revista Ingenio, 15(1), Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 8
Correo electrónico: hayepest@[Link] (MSc. Hernando Alexander Yepes Tumay)
La revisión por pares es responsabilidad de la Universidad Francisco de Paula Santander Ocaña
Artículo bajo la licencia CC BY-NC-ND ([Link]
MSc. Hernando Alexander Yepes Tumay
energéticos es cada vez mayor debido a las ventajas de compuerta para el ingreso de aire y la extracción de gases;
esta técnica, donde se pueden realizar análisis detallados mayores detalles pueden encontrarse en [11]. El esquema
de variables que son experimentalmente muy difíciles de del funcionamiento de este tipo de hornos se muestra en la
Figura 1. Cada cámara tiene la capacidad para almacenar
establecer mejoras dentro de diferentes procesos térmicos. 5 apiles, los cuales consisten en un arreglo organizado de
Con el objetivo de obtener un modelo para sustentar aproximadamente 600 ladrillos en dos paneles.
propuestas de mejoramiento viables para la industria
cerámica Possamai et al. [7] llevaron a cabo el modelado de
un horno para la producción de vidrio pulverizado, evaluando la parte superior y a través de 6 canales en la parte inferior.
tres grados diferentes de aislamiento. El modelo permitió Los denominados apiles se posicionan con una distancia de
separación de 60 cm.
movimiento ligero de la llama hacia uno de los costados del
horno debido a la ubicación de la salida de gases.
Ubaque et al [8] adelantaron un estudio numérico de un
óptimos de temperatura con el objetivo de incinerar residuos
de mayor turbulencia para ubicar los residuos, las cuales
dependen de la ubicación de los ladrillos. Se pudo también
comprobar de manera preliminar que la ubicación de los
desechos no afecta el proceso de cocción. En la misma
línea de optimizar el proceso de combustión Tu et al [9]
llevaron a cabo un estudio numérico sobre el efecto de la
forma de la cámara de combustión, donde encontraron que
al aumentar el ángulo de inclinación del techo y suelo de
la cámara se genera una mayor recirculación y por tanto
se obtiene una zona de reacción mucho más amplia con
menores temperaturas y bajas emisiones de NOx Liu et
al [10] utilizaron la herramienta CFD para estudiar las
características de calentamiento y uniformidad de losas
de acero en un horno de recalentamiento. Los resultados
permitieron evidenciar que durante el calentamiento de las
losas se generan altos gradientes de temperatura debido a la
Figura 1.
2.2 Dominio computacional
Bajo estas consideraciones se plantea el presente estudio,
La simulación se realizó usando el programa Ansys-Fluent
en el cual se espera evaluar numéricamente un horno industrial
19.0. Considerando el arreglo de los ladrillos, la cantidad
de apiles y que las condiciones en cada cámara se pueden
la modelación del proceso de cocción utilizando CFD,
considerar informes, el estudio de una sola de estas, es
una representación acorde de los procesos que se dan en
campos de temperatura al interior del horno.
la totalidad el horno. Bajo esta premisa el dominio de
simulación fue acotado a una longitud de 5.4 m lo cual abarca
2. Metodología una cámara completa (4.9 m) y una zona de consideración al
ladrillos ubicados en su interior.
El presente estudio se llevó a cabo tomando como dominio
el horno de la empresa “Ladrillera Ocaña”. El sistema está Con el objetivo de abordar tanto el efecto de calentamiento
constituido por dos secciones rectangulares de 44.1 m de al interior de los ladrillos como la interacción de los
largo, cada una de ellas está divida en 12 cámaras donde elementos presentes dentro de la cámara, en las simulaciones
son ubicados los ladrillos, cada cámara cuenta una puerta de desarrolladas se consideró cada uno de los apiles completos
principal para la carga y descarga de estos, además de una como toda una zona maciza, siendo un total de 4 apiles en
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 9
Estudio numérico del proceso de cocción de ladrillos al interio
la longitud de dominio determinada. Los detalles
un excesivo número de celdas y por tanto un incremento en
el costo computacional. Si bien en algunos estudios como el
desarrollado por Ubaque y colaboradores [8] implementan
estos detalles, utilizan simetrías junto con celdas tetraédricas,
las cuales incrementan la incertidumbre numérica, mientras
en el presente estudio se usaron únicamente celdas
hexaédricas.
Las puertas de carga y de ingreso de aire, se encuentran
cerradas al momento del calentamiento y la cocción, etapas
que son el enfoque del presente estudio, por lo cual son
omitidas, es decir que la pared lateral donde estarían ubicadas
se considera completamente plana. Si bien ambas compuertas
al tener un espesor menor al de las paredes del horno pueden
incidir en la transferencia de calor por conducción que tiene Figura 2. Mallado hexaédrico.
lugar a través de estas, su área de interacción es mucho Por otra parte, en las zonas cercanas a las paredes se
menor al total de las fronteras por lo cual la incidencia puede
considerarse muy pequeña. limite y los gradientes que se forman cerca de las interfaz
2.3 Mallado del dominio en parte de las zonas mencionadas.
Para la simulación se construyó un mallado con 2.062.945
de celdas hexaédricas, el cual cuenta con un Equiangle
Skewness(Eqsk) máximo de 0.67. La distribución entre
el número de celdas junto con el rango de valores de este
parámetro y su apreciación cualitativa de calidad [12] se
muestran en la Tabla 1.
Tabla 1: Distribución de celdas y calidad del mallado
De esta manera el mallado construido cuenta con un 99.16% Figura 3.
de la totalidad de sus celdas con un nivel de calidad excelente De igual forma los ladrillos al interior del horno también
o buena, donde incluso la mayoría está catalogada con el fueron incluidos dentro del dominio y por tanto fueron
más alto estándar, lo cual permite establecer una adecuada mallados. En la Figura 4 se muestra los bloques constituidos
discretización espacial que permite evitar problemas por las celdas hexaédricas.
de solución a nivel numérico durante la simulación.
Adicionalmente el incremento en el tamaño de celda usado
no supera el 20%. El dominio computacional con su mallado
se muestra en la Figura 2.
Figura 4. Bloques de ladrillos mallados.
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 10
MSc. Hernando Alexander Yepes Tumay
2.4 Mallado del dominio en la cámara fue tenida en cuenta mediante la variación
Las dimensiones generales del dominio computacional se de la temperatura. Esta aproximación es comúnmente
muestran en la Figura 3, los detalles pueden encontrarse en usada en procesos de simulación debido al requerimiento
[11]. Para las simulaciones desarrolladas se utilizaron 3 tipos computacional que implica la consideración del proceso de
combustión.
de salida y pared. En las fronteras de esta última tipología se
La composición de los gases de combustión fue
lo reportado en los estudios previos [11-13]. En la Figura 5 determinada usando los valores promedio de las mediciones
se muestra de forma esquemática las diferentes condiciones realizadas por Guerrero [11] en la chimenea del horno,
de frontera usadas durante la simulación con un esquema de estableciendo una composición de 7.12% CO2, 5.32% H2O,
6.46% O2, 4.9% CO y 76.20% N2 (% en volumen).
dominio, el color rojo señala las salidas de este y el color
blanco está asociado con las fronteras de tipo pared. Tal
como se puede apreciar existe una sola entrada y salida del
dominio considerado, sin embargo, existen diferentes tipos negativo de 10.000 Pa, emulando la presión de succión que
de pared dada su ubicación y condiciones térmicas. se genera dentro del horno producto del extractor. Tanto en
las entradas como salidas se usó el diámetro hidráulico y la
intensidad de turbulencia (IT) como condiciones de frontera
En la Tabla 2 se establecen las disposiciones térmicas
dos primeras letras hacen referencia a la pared de acuerdo
con lo expuesto en la Figura 5 y PO indica la pared opuesta
a la Pared Puertas. La tercera hace referencia al tipo de
condición, donde C indica convección, A si se considera
Tabla 2: Disposiciones térmicas en las paredes.
Figura 5. Condiciones de frontera en la simulación.
La entrada de gases corresponde a los productos de
combustión que ingresan al dominio proveniente de las
másico correspondiente y la composición de estos se utilizó
los datos reportados por Guerrero [11]. Considerando una
potencia de operación de 4310.2 kW y un factor de aireación
condición de frontera “Entrada de gases”, si bien es cierto
volátiles liberados durante la combustión del carbón y su
reportados en la literatura, congruentes para una convección
consumo mismo, esto representa una fracción muy pequeña
natural [16-17]. En PO se asumió un valor menor, dado que
proceso de combustión del carbón tampoco fue considerado
en las simulaciones realizadas ya que determinar la cantidad
aire y por tanto la velocidad es mucho menor. El piso del
exacta de volátiles, así como la tasa de liberación de forma
horno fue considerado adiabático de acuerdo a lo reportado
precisa, el análisis último para determinar la composición
en diferentes estudios [18-19], mientras que las paredes de
química y tamaño de partícula supone un alcance que va más
los tubos se considera esta misma condición dado que están
allá de lo abordado en la presente investigación. El efecto
inmersos dentro del ladrillo cerámico lo que permite un
del combustible y su descarga sobre los apiles considerados
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 11
Estudio numérico del proceso de cocción de ladrillos al interio
aislamiento bastante alto, permitiendo que la de capas límite y altos gradientes tanto de velocidad como
aproximación considerada sea apropiada. de temperatura que pudiese afectar los resultados de la
simulación.
2.5 Modelos complementarios Tabla 3: Número de celdas de los mallados usados.
El tensor de Reynolds fue calculado usando el modelo κ-ε
junto con las funciones de pared estándar. La selección de este
modelo se hizo basado en los resultados adecuados obtenidos
en diversas investigaciones ,[18]. Para la transferencia de
calor por radiación el modelo P1 fue seleccionado debido a
su buen desempeño en estudios anteriores [19] y simplicidad
computacional. Las propiedades de radiación como la
absortividad fueron determinadas usando el modelo de suma
ponderada de gases grises (wsggm) el cual es adecuado
cuando se tienen mezclas como es el caso de los productos
de combustión.
2.6 Procedimiento de simulación
La simulación se llevó a cabo en estado transitorio para Para establecer la independencia del mallado se
evaluar el proceso de calentamiento. El número de Courant compararon las temperaturas y velocidades en diferentes
convectivo en cada celda fue menor a 40 en todo el dominio.
Para tener en cuenta las etapas descritas en la Figura 1, esta región sólida. De forma cualitativa en la Figura 6 se muestran
fue dividida en tres fases. La primera abarca los efectos los contornos de temperatura en la mitad del segundo apile
iniciales hasta el acercamiento medio de la llama y por tanto de ladrillos pasada 1 hora. Es evidente que no existe una
un incremento moderado en la temperatura de los gases de variación apreciable en la distribución de esta variable entre
combustión. La segunda tiene encuentra el momento donde los tres mallados considerados.
la llama se encuentra en las inmediaciones de la cámara lo
que produce un incremento de temperatura de los gases y
encuentra en la cámara y por tanto la temperatura es máxima.
En la primera fase se hizo un aumento paulatino de los
gases de combustión desde 343K hasta 473K en un periodo de
tiempo total de 1.2 horas. Para la segunda fase la temperatura
de los gases fue incrementada a un 60% de la temperatura
de llama adiabática del combustible; el periodo de tiempo
para esta etapa fue de 5 horas. Finalmente, en la tercera fase
se considera un incremento de 200K en la temperatura de
los gases debido a la presencia de la zona de reacción en la Figura 6. Contornos de temperatura al interior de los ladrillos
cámara, la cual se realizó durante un tiempo de 8.9 horas. para los 3 mallados.
Para todos los contornos de temperatura presentados en
En cada paso de tiempo de la simulación se realizaron este trabajo, la escala de temperatura contemplada abarca de
300K a 1575K, a menos de que se indique lo contrario.
igual forma se tomó como criterio de convergencia que los
residuales para la energía y la radiación estuviesen en el En la Figura 7 se muestran las temperaturas y velocidades
orden de 10-6 o menor, en el caso de las velocidades se tomó en función del tiempo usando los tres mallados contemplados
de 10-4, mientras que para la continuidad el valor fue de 10-5 en dos localizaciones. La primera de estas es entre los dos
[21]. primeros apiles de ladrillos en el centro del tercer canal a
una altura de 30 cm, mientras que la segunda se localiza
2.7 Independencia de mallado a la misma altura y posición transversal, pero en medio
Las simulaciones de comparación se realizaron usando 3 de los apiles 2 y 3. Tal como se puede apreciar, las líneas
mallados diferentes. En la Tabla 3 se muestra el número de de temperatura en ambos puntos para los 3 mallados se
superponen casi en su totalidad, con diferencias máximas
se hizo tomando como criterio las regiones de formación punto 1 y el punto 2 respectivamente.
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 12
MSc. Hernando Alexander Yepes Tumay
Figura 8. Líneas de corriente pasadas 5.2 horas.
Se puede apreciar como los gases que viajan por la parte
inferior del horno se enfrían con mayor rapidez, producto de
una menor cantidad y probablemente una velocidad mayor
lo que aumenta la convección. Sin embargo, lo ideal sería
que la tasa de reducción de temperatura de los gases fuese
uniforme para conseguir un tratamiento similar en todos
los ladrillos. Este comportamiento sugiere la necesidad de
de los ladrillos que permita una mayor uniformidad. Entre
los espacios de separación de los apiles se evidencia
zonas de recirculación que pueden ayudar a uniformizar
la temperatura sobre los ladrillos. No obstante, como se
calentamiento uniforme de la carga.
En la Figura 9 se muestran los contornos de temperatura
En cada caso los ladrillos del extremo derecho son los
primeros en tener contacto con los productos de combustión.
Figura 7. Comparación de temperatura y velocidad en puntos
temperatura relevantes asociados con la ubicación de los
ladrillos dentro de un mismo apile. Sin embargo, se observan
Para todos los contornos de temperatura presentados en diferencias de temperatura en el interior de los mismos, lo
este trabajo, la escala de temperatura contemplada abarca de
300K a 1575K, a menos de que se indique lo contrario.
En la Figura 7 se muestran las temperaturas y velocidades
en función del tiempo usando los tres mallados contemplados
en dos localizaciones. La primera de estas es entre los dos
primeros apiles de ladrillos en el centro del tercer canal a
una altura de 30 cm, mientras que la segunda se localiza
a la misma altura y posición transversal, pero en medio
de los apiles 2 y 3. Tal como se puede apreciar, las líneas
de temperatura en ambos puntos para los 3 mallados se
superponen casi en su totalidad, con diferencias máximas
punto 1 y el punto 2 respectivamente.
Figura 9. Contorno de temperaturas ladrillos.
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 13
Estudio numérico del proceso de cocción de ladrillos al interio
efectos sobre la productividad, es claro que la simulación
Las mencionadas variaciones internas pueden observarse
con mayor detalle en la Figura 10, donde se muestra el en un sistema real, teniendo en cuenta además que no es un
contorno de temperaturas en plano transversal en la mitad horno de tipo laboratorio sino un equipo industrial.
alcanza temperaturas entre 1100 y 1200 K pasadas más de
13 horas, el interior se encuentre entre 650 y 750K. Además, 4. Conclusiones
los ladrillos que fungen como pilares del apile formando los En el presente estudio se realizó la simulación de un horno
canales inferiores alcanzan temperaturas cercanas a los 1350
K en el interior, es decir incluso mayores a las temperaturas
como dominio de cálculo una de las cámaras de cocción.
temperatura al interior de los ladrillos. Las líneas de corriente
y la temperatura de los gases el interior del horno evidencia
que el proceso de transferencia de calor ocurre de forma no
uniforme, siendo aparentemente mayor en la zona inferior
donde los gases se enfrían con mayor rapidez asociados a
un aumento en la velocidad que mejora la transferencia de
calor por convección. La distribución de temperaturas de
los ladrillos expone gradientes de temperatura importantes
sugiere que es necesario cambiar la distribución de
los ladrillos para logar una mayor uniformidad en las
Figura 10. Temperaturas al interior de los ladrillos.
temperaturas durante el proceso de cocción.
apiles en la Figura 11 se exponen los contornos de temperatura 5. Referencias
en un plano medio vertical en dirección longitudinal para [1] 2016 IEO, International Energy Outlook 2016. U.S.
los tiempos asociados con la etapa de cocción. Al cabo de Department of Energy, 2016.
las 13 horas las diferencias entre las temperaturas internas [2] L. Balza, R. Espinasa, and T. Serebrisky, Lights On?:
Energy Needs in Latin America and the Caribbean to
transversal, siendo en promedio la temperatura interna 2040. Inter-American Development Bank, 2016.
de 850K. Esta reducción en la diferencia de temperaturas [3] Agencia Internacional de Energía, “Recomendaciones
puede asociarse un espesor menor y las zonas de circulación
2014.
[4] Research and Innovation in Advanced Combustion of
los ladrillos. Industrial Use - INCOMBUSTION, “Determinación
del potencial de reducción del consumo energético
en los subsectores manufactureros códigos CIIU 10-
18.” UPME, p. 209, 2014, [Online]. Available: http://
[Link]/Estudios/2014/INFORME_III_
Caracterizacion_energetica_VerPub.pdf.
[5] A. Amell, L. Rubio, Y. Cadavid, and C. Echeverri,
Figura 11. Temperaturas plano vertical longitudinal. necesidades tecnológicas y planes de acción
tecnológica para mitigación del cambio climático.”
Ministerio de Medio Ambiente y Desarrollo
térmicos, se hace evidente que la distribución de los ladrillos Sostenible, 2013.
al interior del horno puede no ser las más adecuada en [6] J. Campos, E. Lora, L. Meriño, I. Tovar, andA. Navarro,
términos del proceso de cocción, aunque su selección puede “Ahorro de energía en la industria cerámica,” Un
Proy. la unidad planeación Min. energética Colomb. y
el Inst. Colomb. para el Desarro. la Cienc. y la Tecnol.
alternativas requiere un estudio más detallado junto con los
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 14
MSc. Hernando Alexander Yepes Tumay
Fr. Jose Caldas (Colciencias).(2014, octubre), 2010. Eng. Appl. Sci., vol. 11(2), pp. 1045–1061, 2017.
[7] T. S. Possamai, R. Oba, and V. P. Nicolau, “Numerical [19]
and experimental thermal analysis of an industrial
kiln used for frit production,” Appl. Therm. Eng., transfer in an industrial continuous furnace,” Appl.
vol. 48, pp. 414–425, Dec. 2012, doi: 10.1016/j. Therm. Eng., vol. 117, pp. 263–274, May 2017, doi:
applthermaleng.2012.05.025. 10.1016/[Link].2017.02.031.
[8] C. A. G. Ubaque, J. C. M. Piraján, and M. L. V. [20] O. Díaz-Ibarra, P. Abad, and A. Molina, “Design
Bohórquez, “Incineración de residuos sólidos of a day tank glass furnace using a transient model
Therm. Eng., vol. 52(2), pp. 555–565, Apr. 2013, doi:
16, pp. 163–172, 2012. 10.1016/[Link].2012.11.018.
[9] Y. Tu, H. Liu, S. Chen, Z. Liu, H. Zhao, and C. [21] C. Lezcano, J. Obando, and A. Amell, “Numerical
MILD combustion of natural gas,” Appl. Therm. mixtures of methane and a sub-bituminous pulverized
Eng., vol. 76(0), pp. 64–75, 2015, doi: 10.1016/j. Coal,” Combust. Sci. Technol., vol. 189(9), pp. 1591–
applthermaleng.2014.11.007. 1604, 2017, doi: 10.1080/00102202.2017.1305370.
[10] Y. J. Liu, J. D. Li, R. D. K. Misra, Z. D. Wang, and [22] G. Guerrero-Gómez, E. Espinel-Blanco, & T.
G. D. Wang, “A numerical analysis of slab heating Velásquez-Pérez, “Análisis isocinético y corrección a
characteristics in a rolling type reheating furnace condiciones de referencia en horno a cielo abierto en
with pulse combustion,” Appl. Therm. Eng., vol. el municipio de Ocaña, Norte de Santander”, Revista
107, pp. 1304–1312, Aug. 2016, doi: 10.1016/j. Ingenio, vol. 14(1), pp. 43–51, jul. 2017. Doi: https://
applthermaleng.2016.07.074. [Link]/10.22463/2011642X.2194
[11] G. Guerrero Gómez, “Hornos de combustión:
alternativas de reducción de emisiones e incremento
de rentabilidad,” Bucaramanga: Universidad De
Santander, 2017.
[12] H. Meng, G. Zhu, Y. Yu, Z. Wang, and J. Wu,
the turbulent heat transfer in the static mixer with
vol. 99, pp. 647–659, Aug. 2016, doi: 10.1016/J.
IJHEATMASSTRANSFER.2016.03.110.
[13] S. A. Jácome Manzano, “Evaluación termodinámica
del proceso de cocción y análisis de gases en hornos
2012.
[14] T. Ishii, C. Zhang, and Y. Hino, “Numerical study
of the performance of a regenerative furnace,”
Heat Transf. Eng., vol. 23(4), pp. 23–33, 2002, doi:
10.1080/01457630290090473.
[15] T. L. Bergman, F. P. Incropera, D. P. DeWitt, and A. S.
Lavine, Fundamentals of heat and mass transfer. John
Wiley & Sons, 2011.
[16] A. Tehzeeb, M. Bhuiyan, and N. Jayasuriya,
“Evaluation of brick Kiln performances using
Environ. Eng., vol. 1(2), pp. 86–93, 2012.
[17] W. Ratanathavorn, S. Charoenjai, S. Janbuala, B.
of Design Parameters for Clay Brick Kiln Using
Computational Fluid Dynamics and Experimental
Design.,” Adv. Mater. Res., vol. 1101, 2015.
[18] A. Beyene, R. Venkata, and S. Getachew, “CFD
Simulation of Biogas Fired Clay Brick Kiln,” Am. J.
Revista Ingenio, vol. 15, n°1, Enero-Diciembre 2018, pp. 8-15, ISSN 2011-642X – E-ISSN 2389-864X 15
