Edison Collinao PDF

Software para Estudios de Conducción de Calor
Transiente en Procesos Térmicos Industriales
Trabajo para optar al Título de:

Ingeniero Mecánico
Profesor Patrocinante:
Sr. Rogelio Moreno Muñoz
Ing. Civil Mecánico, M.Sc
Profesor Co - Patrocinante:
Sr. Jorge Morales Vilugrón
Ing. (e) en Electricidad
EDISON ADRIAN COLLINAO ROBLES

Valdivia - Chile
2003
Los Profesores Patrocinantes y los Profesores Informantes del trabajo de
titulación, comunican al Director de la Escuela de Mecánica de la Facultad de Ciencias
de la Ingeniería que el trabajo de titulación del Sr.:
EDISON ADRIAN COLLINAO ROBLES
Ha sido aprobado en el examen de defensa rendido el día , como

requisito para optar al titulo de Ingeniero Mecánico. Y, para que así conste para todos
los efectos firman:
Profesor Patrocinante:
Sr. Rogelio Moreno Muñoz, Ing. Civil Mecánico, M.Sc ______________
Profesor Co-Patrocinante:
Sr. Jorge Morales Vilugrón, Ing. (e) en Electricidad ______________
Profesor Informantes:
Sr. Juan Carlos Lehmann, Ing. Civil Mecánico ______________
Sr. Rolando Ríos Rodríguez, M.Sc en Ing. Mecánica ______________
V°B° Director de Escuela
Sr. Enrique Salinas A. ______________

AGRADECIMIENTOS
Doy gracias a Dios en primer lugar por su Don inefable, quien siempre
inmutablemente a estado proveyéndome de todas las cosas.
Extiendo mi aprecio sincero a Don Rogelio Moreno, por su apoyo y

colaboración prestada en el desarrollo del trabajo de titulación, por su alto espíritu
académico demostrado en todas las etapas de esté, así como en el transcurso de la
carrera.
Además, agradezco a Don Jorge Morales, por su aporte en la conformación

del programa computacional. También a los profesores de la Facultad quienes de una
u otra forma colaboraron en mi formación profesional.
Por ultimo, agradezco a la Corporación de Desarrollo Indígena CONADI y al

Ministerio de Educación MINEDUC, por las becas recibidas en todo el periodo de
mis estudios.
DEDICATORIA
En está, una de las etapas de mi vida,

dedico este logro, a mis queridos tíos Gloria y
Carlos, quienes desde el comienzo y hasta
ahora me entregan sus consejos y gran ayuda
permanente.
A mis queridas abuelas Matilde y

Silvia, tía Lucy y a todos mis otros familiares
quienes estuvieron y están siempre
apoyándome en todo.
A mis Padres, por el apoyo en todo lo

posible. En forma especial a mi Madre por su
amor incondicional, la obtención de esté título
es también su logro.
ÍNDICE DE CONTENIDOS
NOMENCLATURA iv
RESUMEN 1
SUMMARY 2
Capítulo 1 – INTRODUCCIÓN 3
1.1 Introducción 3
1.2 Objetivos 5
1.3 Metodología de Trabajo 6
Capítulo 2 – REVISIÓN BIBLIOGRÁFICA 7

2.1 Introducción 7
2.2 Descripción de algunos Procesos Térmicos de Calentamiento 7
y/o Enfriamiento
2.2.1 Proceso de Macerado 7
2.2.2 Proceso del Temple 9
2.2.3 Proceso de Escaldado e Hidroenfriado 11
2.3 Estudios Realizados y Nivel Actual 13
2.4 Métodos Numéricos 16
2.4.1 Método de Diferencias Finitas 16
2.5 Propiedades Termofísicas de las Maderas, Alimentos y Metales 23
2.5.1 Propiedades Térmicas de la Madera 23
2.5.2 Propiedades Térmicas de los Alimentos 25
2.5.3 Propiedades Térmicas de los Metales 27
Capítulo 3 – DESARROLLO DE LOS MODELOS MATEMÁTICOS 29

3.1 Introducción 29
3.2 Suposiciones del Modelo 29
3.2.1 Suposiciones Geométricas 29
3.2.1.1 Placas 30
3.2.1.2 Esferas 30
3.2.1.3 Cilindros 31
3.2.1.4 Paralelogramos 32
3.2.2 Suposiciones de las Propiedades Térmicas y del 32
comportamiento del material
3.2.2.1 Maderas 32
3.2.2.2 Alimentos 34
3.2.2.3 Metales 34
3.2.3 Suposiciones del Medio de Calentamiento y/o Enfriamiento 35
3.3 Modelo Matemático de Conducción de Calor 35
3.3.1 Modelo Matemático para una Placa Extensa 35
3.3.2 Modelo Matemático para una Esfera y un Cilindro Infinito 38
3.3.3 Modelo matemático para un Paralelogramo Infinito 40
3.3.4 Modelo matemático para un Cilindro Finito 42
3.4 Condiciones: Inicial y de Borde 45
3.4.1 Placa Extensa 45
3.4.2 Esfera y Cilindro Infinito 47
3.4.3 Cilindro Finito 51
3.4.4 Paralelogramo 62
3.5 Energía Almacenada /Removida al/del Producto 71
3.6 Condiciones de Estabilidad 72
3.6.1 Introducción 72
3.6.2 Placa Extensa 72
3.6.3 Esfera y Cilindro Infinito 73
3.6.4 Cilindro Finito 74
3.6.5 Paralelogramo 75
Capítulo 4 – IMPLEMENTACIÓN DEL SOFTWARE 76

4.2 Lenguaje de Programación 76
4.3 Algoritmo de Programación 77
4.4 Interfaz de Usuario 80
4.5 Especificaciones y Alcance del Software 84
Capítulo 5 – COMPARACIÓN Y DISCUSIÓN DE LOS RESULTADOS 89

5.2 Comparación de Resultados en Maderas 89
5.2.1 Resultados en el Proceso de Macerado 89
5.2.2 Resultados del Proceso de Calentamiento de Madera con Vapor 95
5.2.3 Análisis del Enfriamiento de una Placa Inclinada 97
5.3 Comparación de Resultados en Alimentos 101
5.3.1 Influencia de las Propiedades Térmicas Variables en Alimentos 104
5.3.2 Influencia del Coeficiente Convectivo en el proceso de Escaldado y 106
Enfriado de Alimentos Vegetales
5.4 Comparación de Resultados en Metales 108
5.4.1 Resultados en el Proceso del Temple 108
5.4.2 Resultados Método Gráfico versus Diferencias Finitas (MDF) 111
Capítulo 6 – CONCLUSIONES 113
REFERENCIAS 115
Referencias Bibliográficas 115
Referencias de Internet 117
ANEXOS 118
Anexo 1 119
Anexo 2 123
Anexo 3 126
NOMENCLATURA
a : Ancho de la placa extensa (m)

A : Ancho paralelogramo (m)
B : Alto paralelogramo (m)
CP : Calor especifico (W/Kg°K)
GS : Densidad básica de la madera
D : Diámetro del sólido (m)
h : Coeficiente convectivo de transferencia de calor (W/m²°K)
i, j : Posición nodal denotado por un numero entero
g : Aceleración de gravedad (m/s²)
k r (T ) : conductividad térmica radial (W/m°K)
k z (T ) : conductividad térmica axial (W/m°K)
l : Parámetro equivalente al tamaño de la malla (m)
M : Parámetro adimensional del método numérico
M1 : Parámetro adimensional del método numérico
M2 : Parámetro adimensional del método numérico
n : Posición en el tiempo denotado por un numero entero
r : Posición radial sólido (m)
R : Radio cilindro (m)
T : Temperatura del sólido (°K)
Tb : Temperatura del baño o ambiental (°K)
Td : Temperatura deseada (°K)
Ti : Temperatura inicial (°K)
Ti n : Temperatura en el nodo i para el tiempo n = 0 (°K)
Ti n +1 : Temperatura en el nodo i para el tiempo n + 1 (°K)

t : Tiempo (s)
P : Numero de nodos en la dirección r ó x (según sea el modelo)
Q : Numero de nodos en el dirección z ó y (según sea el modelo)
Q(t ) : Energía almacenada/removida al /del sólido (J)
Pr : Numero de Prandtl
Gr : Numero de Grashof
Bi : Numero de Biot
Nu : Numero de Nusselt
H : Contenido de humedad (%)
SÍMBOLOS GRIEGOS
αr(T ) : Difusividad térmica radial (m²/s)

αz(T ) : Difusividad térmica axial (m²/s)
β : Coeficiente de expansión volumétrica (1/K)
ρ : Densidad (Kg/m³)
∆Q : Incremento de energía en el sólido (J)
∆r : Tamaño de un nodo radial (m)
∆t : Incremento de tiempo (s)
∆x : Tamaño de un nodo en x (m)
∆y : Tamaño de un nodo en y (m)
∆z : Tamaño de un nodo axial (m)
ν : Viscosidad cinemática del fluido (m²/s)
1
RESUMEN
Existen numerosas aplicaciones de ingeniería, donde se encuentra el uso de

procesos térmicos de conducción de calor transiente. En Chile, la industria de la madera
debe tratar con el calentamiento de trozas en agua, en el caso de la producción de los
alimentos vegetales se recurre a los procesos de escaldado y enfriado y en la industria
metalmecánica a procesos como el temple.
El presente trabajo, consiste en elaborar un software con fines académicos, para

estudiar los procesos mencionados, y otros. Para ello se desarrollan modelos de
conducción de calor inestables lineales y no lineales, para predecir los tiempos de
calentamiento y/o enfriamiento y calor añadido/removido a/de productos de geometría
simple: cilindros, placas, esferas y paralelogramos. En la frontera de los sólidos, es
posible considerar y evaluar el coeficiente de transferencia de calor por convección, a
través de los mecánicos de convección natural y forzada.
Los modelos fueron resueltos a través del método de diferencias finitas en su

forma explicita, siendo implementados en un programa computacional, el cual fue
confeccionado en el lenguaje Visual Basic.
Los resultados obtenidos en las áreas de estudio (maderas, alimentos, metales),

son bastante satisfactorios al realizar las respectivas comparaciones. En maderas se
obtuvieron desviaciones alrededor del 2% con datos teórico de otro autores, un 1,2%
para el caso de alimentos vegetales expuestos al proceso de escaldado y enfriado al
comparar con datos experimentales. En el proceso del temple los resultados mejoran
considerablemente al modelar a propiedades térmicas variables y coeficiente convectivo
variable.
En términos generales, el software ofrece resolver una variedad de problemas, sin

errores significativos, si estos incluyen los efectos de todas las variables que afectan en
la estimación del tiempo de calentamiento y calor absorbido/removido al/del producto.
2
SUMMARY
There are many engineering applications where the use of thermal processes of
transient heat conduction is found. In Chile, the wood industry must deal with the
heating of logs in water. In the case of the production of vegetable food, it is appealed to
the processes of blanching and cooling and in the metal mechanic industry, to processes
like tempering.
The present work consists of elaborating a software with academic purposes, to

study the already mentioned processes and others. So, lineal and not lineal unstable
models of heat conduction are developed so as to predict the heating and cooling times
and also the added and removed heat to simple geometry products: cylinder, plates,
spheres and parallelograms. In the frontier of the solids, it is possible to consider and
evaluate the coefficient of heating transfer by convection through the mechanics of
natural and forced convection.
The models were solved through the method of finite differences in their explicit
form, being implemented in a computational program which was made in the Visual
Basic Language.
The results obtained in the areas of studies (wood, food, metals), are very
satisfactory when analyzing the corresponding comparisons. In wood, deviations were
obtained around 2% with another author’s theoretical data; a 1,2% in the case of
vegetable food exposed to the blanching and cooling processes when comparing them
with experimental data. In the tempering process, the results improve considerably
when modeling to thermal variable properties and variable convective coefficient.
In general terms, the software offers to solve a variety of problems, without

significant errors if they include the effects of all variables that affect when calculating
the heating times and the added and removed heat of the product.
3
Capítulo 1 – INTRODUCCIÓN
1.1 Introducción
Una de las necesidades en la industria chilena maderera, metalmecánica,

alimentos y otros, es tratar con procesos térmicos de calentamiento y/o enfriamiento en
diversos productos o materiales, con el propósito de modificar sus distintas propiedades
y que se necesitan para obtener ciertas características de mayor calidad, vida útil, o
simplemente optimizar el proceso productivo.
Como ejemplo se tienen algunos de estos procesos contenidos en los campos

industriales ya mencionados y que son el proceso del temple, macerado, escaldado e
hidroenfriado, entre otros.
Específicamente, en el caso de la industria de contrachapado, existe el proceso de

calentamiento de trozas de madera, llamado macerado, que tiene el propósito de
ablandar la madera, adquiriéndose propiedades mecánicas apropiadas para facilitar el
siguiente proceso de corte o debobinado.
El uso adecuado de todos estos procesos y otros, obedece al hecho de plantear

ciertas interrogantes que son de carácter especifico, como por ejemplo, ¿cuál es el
tiempo optimo al que debo exponer el producto o material en el proceso?; o saber
¿cuanto es el calor absorbido por el producto?, a fin de evaluar cual es la energía que se
está consumiendo en el proceso. En consecuencia se hace evidente realizar reportes de
distribuciones de temperatura versus tiempo, y lo mismo para el calor absorbido.
Existen publicaciones de tablas para determinar tiempos de calentamiento de

trozas de madera. Sin embargo, estas no son aplicables a nuestro país ya que la
información entregada solo es valida a condiciones equivalentes a las investigaciones
realizadas en América del norte por FEIGHL (1972) y FEIGHL y GODIN (1975).
4
Por otro lado en alimentos se presentan las mismas limitantes en la información,

existen variados estudios en algunos vegetales, pero no existe evidencia en otros
(espárragos, otros).
Así el objetivo de este trabajo es entregar una solución a los problemas de

conducción de calor transitorios encontrados en los procesos térmicos en unidades de
calentamiento y/o enfriamiento de productos, de procesamiento común. Se formula un
modelo matemático que aborda problemas de geometría simple, tales como cilindros,
placas, esferas y paralelogramos utilizando el método numérico explicito de diferencias
finitas a la solución de este, considerando el comportamiento del material anisotrópico y
con las propiedades térmicas constantes y variables en el tiempo.
Además se define en la interfaz sólido – fluido la transferencia de calor por

convección, a través de los mecanismos de convección natural o forzada. Así en la
primera situación el coeficiente de transferencia de calor por convección que mide la
intensidad de este fenómeno, será variable en el tiempo, el cual se determina a través de
parámetros adimensionales como el numero de Grashoft y Prandtl. Para la segunda
situación el efecto convectivo se asume uniforme y constante en el tiempo y que puede
ser calculado por medio de correlaciones que dependen de los números de Reynolds y
Prandtl.
Se desarrolla una herramienta computacional capaz de simular el

comportamiento de los productos frente al sistema térmico, es decir generar reportes de
distribuciones de temperatura en función del tiempo y calor absorbido por el producto.
Todo para el estudio posterior de la cinética del proceso.
Finalmente se realizan comparaciones entre resultados obtenidos por el software

y los encontrados en la literatura y/o los trabajos teórico – experimentales realizados
anteriormente en la Universidad Austral de Chile.
5
1.2 Objetivos
Objetivo General:
!" Elaborar un software con fines académicos, para simular procesos térmicos
industriales en unidades de calentamiento y/o enfriamiento para productos de
geometría simple, con propiedades térmicas constantes y variables.
Objetivos Específicos:
!" Desarrollar un modelo matemático para describir el comportamiento cinético que

ocurre en los procesos térmicos de calentamiento y/o enfriamiento. Aplicar el
método numérico explicito de diferencias finitas a la solución de problemas de
conducción de calor en estado transiente.
!" Confeccionar un programa computacional capaz de simular el proceso térmico

para geometrías simples tales como cilindros, placas, esferas y paralelogramos,
que permita entregar reportes de distribución de temperatura versus tiempo y calor
añadido/removido al/del material.
!" Comparar los resultados obtenidos con el programa computacional versus los
encontrados en la literatura y/o los experimentales de investigaciones de la
Universidad Austral de Chile.
6
1.3 Metodología de Trabajo
Como uno de los objetivos primarios será aplicar el método numérico a la

solución de problemas de conducción de calor, se formularán las ecuaciones pertinentes
para cada una de las geometrías en cuestión, con sus respectivas condiciones iniciales y
de contorno, posteriormente se recogerán datos de interés, como son las propiedades
térmicas de los productos metálicos, madereros y alimenticios.
Se establecerá un diagrama de flujo o algoritmo para la solución numérica

explicita del problema de conducción de calor transitorio, con esto se desarrollara el
software en ambiente Windows, a través del lenguaje de programación de visual Basic
6.0.
Finalmente se realizará la comparación de los resultados encontrados tanto en

forma experimental y/o literatura con los obtenidos por el software, estableciendo las
respectivas conclusiones.
7
Capítulo 2 - REVISON BIBLIOGRAFICA
2.1 Introducción
El Capítulo presenta una revisión bibliográfica, que comprende una descripción

de algunos procesos térmicos, como son el macerado, temple, escaldado e hidroenfriado,
los cuales ejemplifican los procesos térmicos contenidos en las áreas de la industria de la
madera, metalmecánica, alimentos y otros.
Se revisa además los estudios realizados y el nivel actual, en relación a la

aplicación de los métodos numéricos, particularmente el método de diferencias finitas,
en el estudio de los problemas de conducción de calor transiente, la importancia de estos
en el interés que le compete a la industria.
Además se revisan los métodos numéricos mas utilizados en la solución de los

problemas de transferencia de calor, realizando una descripción de estos, las ventajas y
desventajas de su aplicación.
2.2 Descripción de algunos Procesos Térmicos de Calentamiento y/o

Enfriamiento
2.2.1 Proceso de Macerado
El proceso de macerado de trozas de madera es un proceso térmico utilizado

frecuentemente en empresas de tecnología de productos forestales y cuyo propósito es
ablandar la madera, obteniéndose propiedades mecánicas apropiadas para el
subsiguiente proceso de corte ó debobinado realizado en maquinas foliadoras.
El ablandamiento se lleva a cabo en la mayoría de las industrias nacionales en

piletas de agua caliente (ver Figura1) las cuales son de temperatura variable en el
tiempo, cabe señalar que no todos los regímenes son iguales en cuanto a valores de
8
temperatura y de tiempo de calentamiento, por lo que es conveniente desarrollar un

programa computacional que permita predecir los tiempos de calentamiento y así
obtener las temperaturas óptimas para cada especie.
FIGURA N° 1 Proceso de macerado en piletas de agua

caliente.
El tiempo de calentamiento tiene que ser suficientemente largo para calentar la

madera hasta el rollizo residual, depende de la densidad de la especie y del diámetro del
rollizo. FEIHL (1972) establece que el tiempo de calentamiento aumenta con el
cuadrado del diámetro, en el Cuadro 1 se indican algunos tiempos para diferentes
diámetros y temperaturas para trozas con una densidad verde de 0,5.
CUADRO N° 1 Temperatura y tiempo de cocción según FEIHL (1972).

Diámetro Temperatura de la Temperatura final en el Tiempo
(cm) Pileta rollizo residual (Horas)
35 65 60 13
60 43 38 34
60 65 60 36
60 88 83 36
9
En la etapa final, el proceso de macerado es mas rápido si se mantiene un

gradiente de temperatura de 5°C entre la temperatura de la pileta y la madera. Por otro
lado para un buen proceso de debobinado, es necesaria una temperatura adecuada
dependiendo de la especie (gravedad especifica).
2.2.2 Proceso del Temple
El temple es un tratamiento térmico que tiene por objeto endurecer y aumentar la

resistencia de los aceros. Para ello, se calienta en general el acero a una temperatura
ligeramente más elevada que la critica superior y se enfría luego más o menos
rápidamente (según la composición y el tamaño de la pieza) en un medio conveniente,
agua, aceite, etc.
Teóricamente, en el temple, toda la masa del acero debe encontrarse en estado

austenítico (alta temperatura) en el momento de comenzar el enfriamiento brusco. Si se
enfría con una rapidez optima, todo el acero queda con una gran dureza y con una nueva
estructura llamada martensita. Dependiendo del contenido de carbono en el acero, esta
temperatura varía desde 770 a 970 °C para un contenido de carbono entre 0,9 a 0,10 %.
Para un buen proceso de temple, se debe conseguir una temperatura más elevada
que la critica superior, para que todo el acero quede en estado austenico. El temple debe
cumplir 40° - 60 °C por encima de la temperatura critica (MARTÍNEZ y GORDON,
1990; APRAÍZ, 1974).
Los factores que interviene en el temple se pueden clasificar en dos grupos:
En el primero se encuentran los que son inherentes a la calidad del material y

determinan su velocidad critica de temple. Los mas importantes son: la composición del
producto y el tamaño del grano. El segundo grupo pertenecen los factores que
determinan la velocidad de enfriamiento de los diferentes puntos de la pieza, siendo los
10
más importantes: la forma, dimensiones de las piezas, propiedades térmicas y el medio

de enfriamiento, .
Para este trabajo se orientará el análisis, solo al segundo grupo de factores

mencionados.
(i) Influencia del Tamaño de las piezas: El volumen y espesor de las piezas tiene
gran importancia en el temple de los aceros. Si se estudia el enfriamiento en agua de un
cilindro que ha sido calentado a alta temperatura, fácilmente se comprende que la capa
exterior del acero, que está en contacto con el liquido, es la que se enfría primero y más
rápidamente. La capa siguiente no se enfría con tanta rapidez, porque el calor no se
transmite directamente al liquido sino que debe atravesar antes la capa superficial del
acero, y esta capa tiene una conductibilidad bastante limitada. Fácilmente se comprende,
también, que las velocidades de enfriamiento entre la superficie y el interior de la pieza
de acero de diversos diámetros variaran mucho de unos casos a otros. En consecuencia
se concluye que en algunos casos las capas superficiales de una pieza se endurecen
perfectamente y, en cambio, otras zonas del interior, por enfriarse a velocidades
menores, no templan.
(ii) Influencia del Medio de Enfriamiento: El medio de enfriamiento junto con el

tamaño y forma de las piezas son, factores que deciden la velocidad de enfriamiento.
Con el enfriamiento en agua salada muy agitada, se consiguen las mayores
velocidades, y enfriando las piezas dentro de un horno, se obtienen las menores. Con
enfriamientos en aceite, al aire, en sales, etc., se consiguen velocidades intermedias,
pudiendo elegirse entre todos ellos el medio que mejor cumpla las condiciones de
enfriamiento deseadas en cada caso, APRAÍZ (1974).
En la Figura 2 se ve el proceso de enfriamiento de varios puntos de un cilindro de

50 mm templado en agua (a la izquierda) y aceite (a la derecha).
11
FIGURA N° 2 Curvas de enfriamiento de un cilindro en agua y en

aceite.
2.2.3 Proceso de Escaldado e Hidroenfriado
El proceso de escaldado, consiste en la aplicación de alguna forma de calor,

usualmente agua caliente, seguido de un enfriamiento que constituye una operación
básica, previas al proceso de congelado, deshidratado y enlatado de diversos productos
vegetales. El objetivo principal del escaldado es la reducción de la actividad enzimática
la cual esta en directa relación con la cinética de la transferencia de calor, (LUH y
LORENZO, 1988).
El escaldado busca varios objetivos, entre los cuales se destacan:

- Reducir la actividad enzimática, por destrucción térmica de las enzimas
responsables presentes en los productos que en caso de no aplicar el tratamiento darían a
lugar a aromas, sabores, o colores extraños y que causarían la perdida de la vitamina C,
disminuyendo así su valor comercial, (LUND, 1977; BURNETTE, 1977).
- Expulsar los gases obstruidos en los espacios intercelulares de las hortalizas

(GANTHAVORN y POWERS, 1998), evitando así un aumento de la presión
desarrollada en las conservas durante el procesamiento debido a la expansión de los
gases presentes (SHAMS y THOMPSON, 1987).
12
- Mejorar el sabor y estabilizar el color verde de los vegetales (LUND, 1977;

FENNEMA, 1993).
- Mejorar la textura (FENNEMA, 1993).
Los métodos para escaldar productos vegetales son variados, la inmersión en

agua a temperaturas entre los 80°C y 100°C seguido por inmersión en agua fría es el mas
común, la otras alternativas son a través del vapor, microondas y gas caliente (las cuales
no se describirán en este trabajo). En todo caso ninguna de estas alternativas supera al
método mas común en la industria del escaldado por agua, ya sea por razones de
eficiencia del proceso (perdida de nutrientes), ahorro energético o costos de
implementación y mantenimiento, (AGUILERA, 1996).
En el caso de producciones a escala mediana, como es el caso de nuestro país, se

ocupa frecuentemente el escaldado con agua. Por otro lado, el escaldado continuo,
utilizado para la mayoría de las industrias procesadoras de productos vegetales. En estos
procesos se esta en presencia de convección forzada como tipo de transferencia de calor,
ya que existe una velocidad de flujo, ya sea del producto contra el fluido o del fluido
contra el producto, (ORDÓÑEZ, 1996).
Los tiempos de escaldado e hidroenfriado al igual que otros procesos térmicos

dependen de ciertos factores los cuales son: la temperatura del baño, coeficiente
convectivo de transferencia de calor, temperatura inicial y final del producto, tamaño y
configuración del producto y propiedades termofísicas del producto (densidad, calor
especifico y conductividad térmica).
Los tiempos y temperaturas del proceso de escaldado mencionados en la

literatura para varios productos varían en un amplio rango, el Cuadro 2 muestra algunos
ejemplos.
13
CUADRO N° 2 Tiempos de escaldado en agua a 95°C,

fuente: LUH y LORENZO (1998).
Producto Tiempo (min.)
Apio 2
Alcachofas 5-9
Brócolis 2-3
Coles de Bruselas 4-5
Espárragos 2-5
Espinacas 2
Zanahorias 2-5
Mazorca de Maíz 6-11
2.3 Estudios Realizados y Nivel Actual
La solución de los problemas de conducción de calor transiente con condiciones de

borde o propiedades térmicas variables solo puede ser resuelta a través de los métodos
numéricos.
Trabajos realizados en las ultimas décadas por investigadores chilenos han dado
excelentes resultados en la utilización de los métodos numéricos.
Los investigadores MORENO et al (1992) desarrollaron un modelo para simular el

proceso de calentamiento de trozas de madera llamado macerado, que tiene el propósito
de ablandar la madera para el siguiente proceso de corte ó debobinado y que es utilizado
frecuentemente en industrias de contrachapados, (ver Sección 2.2.1).
El modelo se basa en la ecuación de difusión de calor, asumiendo las propiedades

termofísicas constantes y considerando la conducción de calor axial dos veces más
rápida que la radial, como resultado del comportamiento anisotropico del material. La
solución del modelo matemático pudo realizarse a través del método numérico de
diferencias finitas en su forma explicita, posteriormente se elaboró un programa
14
computacional en lenguaje programación Fortran, que fue capaz de entregar reportes de

temperatura / tiempo, así también de la energía almacenada por la troza.
Finalmente se validó el modelo propuesto, realizando una comparación entre los

resultados obtenidos por el programa y las mediciones en terreno, en ambos se tomaron
cinco especies diferentes de trozos de madera y a diferentes condiciones de operación.
Se concluyó de esto que las desviaciones de los resultados se aceptan en la practica,

además a juicio de los autores no es justificable para una mayor precisión del modelo
considerar las propiedades físicas variables excepto si se conoce previamente el la
distribución del contenido de humedad en el interior de la troza de madera.
Otros estudios realizados en el extranjero por STEINHAGEN et al (1980) muestran

excelentes resultados en la aplicación del método numérico de diferencias finitas en la
determinación de cartas de tiempos de calentamiento para trozas de madera congelados
y no congelados, los tiempos de calentamientos estimados fueron satisfactorios en
comparación con los datos experimentales, existiendo una discrepancia de menos del
10% y que se debió al considerar las propiedades térmicas constantes y además de no
tomar en cuenta el efecto de los gradientes de humedad. Además se señala que para
mayor precisión en la determinación de las relación temperatura – tiempo es posible
utilizar las propiedades variables con ayuda de herramientas computacionales.
SIMPSON (2001) presenta un trabajo reciente que refleja la importancia de

contar con una herramienta que permita predecir el comportamiento de la troza en un
tratamiento térmico dado, se busco generar en tablas, tiempos de calentamiento para
trozas de madera de sección circular y rectangular sometidas a un proceso térmico de
calentamiento en un medio como el vapor, con el fin de conseguir una esterilización de
la madera, es decir eliminar los insectos o patógenos encontrados en ella, todo como
cumplimiento de normas sanitarias en la exportación de estos productos, el método
utilizado.
15
En el área de la metalmecánica MORAGA y SÁNCHEZ (1994) desarrollaron un

modelo mediante la formulación de resistencias y capacitancias térmicas de diferencias
finitas y su comparación con valores experimentales, para la distribución de
temperaturas y velocidades de enfriamiento de una placa y un cubo de acero SAE 1045,
sometidos a un proceso térmico de temple en agua. El método numérico se utilizó para
la solución de problemas de conducción de calor en una, dos y tres dimensiones,
transientes, con propiedades térmicas variables con la temperatura y considerando la
variación en función del tiempo del coeficiente convectivo de transferencia de calor.
ROJAS (2001) desarrolló un modelo a través de la formulación de diferencias

finitas, para simular el proceso de escaldado y enfriado de alimentos vegetales
(zanahorias, espárragos, etc.) basándose en el trabajo de titulación de AGUILERA
(1996) quien a su vez trabajo en un modelo matemático, mediante métodos analíticos
para simular el proceso de transferencia de calor y a través de técnicas numéricas, para
describir la cinética del proceso de escaldado e hidroenfriado. En ambos trabajos, los
modelos siguen la tendencia de las curvas experimentales del proceso.
Otros estudios realizados en el país se han llevado a cabo en el campo de los

alimentos, tienen como común denominador el de determinar las relaciones entre las
condiciones del proceso (temperatura – tiempo) y los cambios de calidad del producto, y
que los investigadores MORAGA et al. (1996, 1997) han utilizado satisfactoriamente el
método de diferencias finitas, validando los resultados con valores experimentales.
Por otro lado a nivel internacional se han realzado diversos trabajos en este
campo, siendo uno de los más interesantes el propuesto por Q. TUAN PHAM (2002), el
cual a generando cálculos en base a los métodos numéricos de diferencia finitas,
elementos finitos y CFD (Dinámica de Fluidos Computacional) para determinar los
tiempos y carga de calor durante el proceso de refrigeración de alimentos, todo con el
objetivo de dimensionar en forma optima los equipos en los procesos de enfriamiento y
congelado (compresores, condensadores). En dicho trabajo se concluye que para la
mayoría de los casos, es una buena practica el uso de las diferencias finitas o elementos
16
finitos en el interior del producto, combinado con el coeficiente de transferencia de calor

obtenido desde ecuaciones empíricas ó simulaciones por CFD.
2.4 Métodos Numéricos
En los últimos años, con el incremento de la velocidad y capacidad

computacional, el interés por el uso de los métodos numéricos a cobrado gran
importancia en la solución de problemas gobernados por ecuaciones diferenciales.
Principalmente para problemas complicados en ingeniería, los cuales con estas
herramientas pueden ser resueltos a un menor costo y tiempo.
En transferencia de calor en la práctica se encuentran problemas que implican

geometrías y condiciones de borde complejas, o propiedades variables, que no es posible
resolver de manera analítica, pero si con la aplicación de los métodos numéricos, además
resulta mas ventajosa en términos del tiempo necesario para llegar a una solución rápida,
también se facilita el cambio de los parámetros del proceso, permitiendo al ingeniero
determinar el comportamiento del sistema térmico con mucha mayor facilidad, KREITH
y BHON (2001).
El análisis numérico se basa en proporcionar una solución solamente para puntos

discretos dentro de la región del problema, ofreciendo así una buena aproximación a la
solución exacta, simplificándose la solución, al resolver ahora un grupo de ecuaciones
recursivas.
2.4.1 Método de Diferencias Finitas
El método de diferencias finitas (MDF) tiene por objetivo resolver ecuaciones

diferenciales parciales, que gobiernan un problema físico en un dominio dado, a través
del proceso de discretización, es decir remplazar en forma aproximada una ecuación
diferencial parcial en una ecuación recursiva para un determinado numero de puntos
17
nodales de una región. Para tal efecto se recurre al uso (i) de la expansión de series de
Taylor, y (ii) a la aproximación del volumen de control
(i) Expansión de Series de Taylor: La idea de la representación de diferencia

finita de una derivada puede ser introducida por el empleo de su definición, una base
formal para el desarrollo de derivadas a la aproximación de diferencias finitas es a
través de la expansión de series de Taylor. Se considera la expansión de una función
f (x ) sobre un ponto x , en las direcciones adelanta y atrasada, respectivamente, vea la
Figura 3.
h2 h3
f (x + h ) = f (x ) + hf ′(x ) + f ′′(x ) + f ′′′(x ) + ... (2.1)
2! 3!
h2 h3
f (x − h ) = f (x ) − hf ′(x ) + f ′′(x ) − f ′′′(x ) + ... (2.2)
2! 3!
FIGURA N° 3 Representación esquemática de las

series de Taylor.
Despejando de las ecuaciones (2.1) y (2.2) se obtienen las primeras derivadas de

la aproximación de diferencia finita adelantada y atrasada, respectivamente.
f (x + h ) − f (x )
f ′(x ) = + 0(h) Adelantada (2.3)
h
18
f (x ) − f (x − h )
f ′(x ) = + 0(h) Atrasada (2.4)
h
Donde la notación 0(h) caracteriza el error de truncamiento asociado a la

aproximación de diferencia finita y que representa la diferencia entre la derivada y la
representación de diferencia finita.
Ahora restando la ecuación (2.2) a la (2.1), se obtiene la aproximación de la

diferencia central, quedando:
f (x + h ) − f (x − h )
f ′(x ) = + 0(h) 2 Centrada (2.5)
2h
Como se ve la aproximación de diferencia es de segundo orden 0(h) 2 en h, por

lo que la aproximación central resulta más exacta que las diferencias adelantada y
atrasada.
Ahora si se introduce la siguiente notación:
x = ih , x + h = (i + 1)h , x − h = (i − 1)h
f (x ) = f i , f (x + h ) = f i +1 , f (x − h ) = f i −1
Las segundas derivadas se obtienen de la misma forma que las primeras con
excepción de la aproximación central que se determina sumando las ecuaciones (2.1) y
(2.2).
En resumen, se presenta a continuación el Cuadro 3 con las primaras derivadas y
segundas ya sea en su forma adelantada, atrasada y central, para problemas transientes,
útiles en el reemplazo de las ecuaciones diferenciales parciales que gobiernan el proceso
de conducción de calor por diferencias finitas que para este trabajo son en su forma
central.
19
(ii) Aproximación de Volumen de Control: Anteriormente se menciono que a

partir de una ecuación diferencial parcial apropiada para un problema gobernado por
leyes físicas, se pueden utilizar las series de Taylor como un procedimiento puramente
matico para desarrollar la aproximación de diferencia finita.
En la alternativa de la aproximación de volumen de control, la ecuación de

diferencia finita se desarrolla por la restricción de la ecuación diferencial parcial a un
volumen de control y conservar la cantidad física tal como masa, momentum o energía
sobre el volumen de control. El concepto básico es análogo al escribir un balance
térmico sobre un pequeño volumen alrededor de un punto de una malla (ÖZISIK, 1994).
CUADRO N° 3 Primeras y segundas derivadas en aproximaciones de

diferencia finita.
Aproximación Primeras Derivadas
f i − f i −1
Atrasada f i′ =
h
f i +1 − f i
Adelantada f i′ =
h
f i +1 − f i −1
Central f i′ =
2h
Aproximación Segundas Derivadas
f i − 2 − 2 f i −1 + f i
Atrasada f i '' =
h2
f − 2 f i +1 + f i − 2
Adelantada f i '' = i
h2
f − 2 f i + f i +1
Central f i '' = i i −1
h2
Ahora bien se disponen varios esquemas de representar las derivadas como son el
método explicito, implícito, Crank-Nicolson, Combinado, siendo el utilizado en este
trabajo el esquema explicito. Sin embargo a continuación se explicaran las ideas básicas
de algunos para un mejor entendimiento, para ello se considera un problema
unidireccional de conducción de calor transiente.
20
Sea:
∂T (x, t ) ∂ 2 T (x, t )
=α 0 ≤ x ≤ L, 0 < t (2.6)
∂t ∂x 2
Para una región 0 ≤ x ≤ L dividida en M partes iguales de una malla:

L
∆x =
M
Método Explicito Simple
Un sistema numérico explicito es aquel mediante el cual se puede obtener

directamente el valor desconocido, a partir de valores conocidos en un tiempo presente,
se calcula entonces el valor de la temperatura en un tiempo futuro, Es el método que es
utilizado en el este trabajo, ya que permite variar las propiedades termofísicas de los
materiales y los coeficientes convectivos de transferencia de calor a cada paso de
tiempo.
Así entonces la ecuación 2.6 queda representada como sigue:
Ti n +1 − Ti n Ti −n1 − 2Ti n + Ti +n1

=α (2.7)
∆t ∆x 2
Ti n +1
∆t
Ti −n1 Ti n Ti +n1
∆t
∆x ∆x
FIGURA N° 4 Esquema del método de diferencias

finita en su forma explicita
21
Una característica inconveniente del método explicito es que no es

incondicionalmente estable. Es decir en un problema transitorio, la solución para la
temperatura nodal Ti n +1 debe aproximarse de forma continua a los valores finales (de
estado estable) al aumentar el tiempo. Sin embargo, con el método explicito, bajo ciertas
condiciones relativas a la malla ∆x y el incremento de tiempo ∆t esta solución se
caracteriza por oscilaciones numéricamente inducidas. Las oscilaciones se vuelven
inestables (ver Figura 5), lo que ocasiona que la solución difiera de las condiciones de
estado estable reales.
FIGURA N° 5 Efecto de la estabilidad en la

solución del método explicito.
Para evitar este tipo de errores, el valor establecido de ∆t debe mantenerse en

ciertos limites, el cual depende de ∆x y otros parámetros del sistema. Tal dependencia
se le llama criterio de estabilidad y que para los problemas tratados en este trabajo se
obtiene requiriendo el coeficiente asociado con el nodo de interés en el tiempo anterior
sea mayor o igual a cero. Así se necesita reunir todos los términos que incluyen Ti n
para obtener el cuociente, del cual se determina el valor permisible de ∆t .
Método Implícito Simple
A diferencia del método explicito el implícito, necesita resolver

simultáneamente la ecuación de diferencia, a partir del valor conocido del nodo de
interés en el tiempo presente, la ventaja de este método es que no presenta restricciones
22
en cuanto al valor de ∆t como una consideración en la estabilidad de la solución (ver

Figura 6).
La ecuación (2.6) queda como:
Ti n +1 − Ti n T n +1 − 2Ti n +1 + Ti +n1+1
= α i −1 (2.8)
∆t ∆x 2
Ti −n1+1 Ti n+1 Ti +n1+1
∆t
Ti n
∆t
∆x ∆x

finita en su forma implícita.
Método Crank-Nicolson
La idea básica del esquema propuesto por Crank-Nicolson, es modificar el lado

derecho de la ecuación (2.6) tomando la media aritmética de las ecuación (2.7) en su
esquema explicito y (2.8) en su forma implícita, entonces la aproximación de diferencia
finita de la ecuación de conducción de calor queda de la siguiente manera:
Ti n +1 − Ti n α  Ti −n1+1 − 2Ti n +1 + Ti +n1+1 Ti −n1 − 2Ti n + Ti +n1 

=  +  (2.9)
∆t 2 ∆x 2 ∆x 2 
Así entonces se requiere resolver un sistema de ecuaciones simultaneas tanto el

nodo central como los adyacentes (ver Figura 7). La ventaja al igual que el esquema
implícito es que la solución siempre converge, pero su uso necesita mayores cálculos
matemáticos.
23
Ti−n1+1 Ti n+1 Ti+n1+1
∆t
Ti−n1 Ti n Ti +n1
∆t
∆x ∆x

finita en su forma Crank - Nicolson.
2.5 Propiedades Termofísicas de las Maderas, Alimentos y Metales
Resulta imprescindible conocer con precisión las propiedades tales como la

densidad, conductividad térmica y calor especifico cuando se quiere realizar un análisis
térmico, sobre todo si los productos están frente a un proceso de transferencia de calor
transiente, donde la temperatura cambia con el tiempo y por consiguiente las
propiedades mencionadas también. A continuación se revisaran tales parámetros en
cada una de las áreas industriales que contempla este trabajo.
2.5.1 Propiedades Térmicas de la Madera
Los datos de la conductividad térmica y el calor específico dependen de varios

factores. La conductividad térmica de la madera esta afectada por un numero básico de
factores: densidad, contenido de humedad, dirección del grano (radial, tangencial,
longitudinal), estructuras irregulares y temperatura (SIMPSON, 1993; STEINHAGEN,
1977).
La conductividad térmica en la dirección radial es casi igual a la tangencial, en

tanto que la conductividad en la dirección longitudinal se reporta que es entre 1,5 a 3
veces mayor que la conductividad radial, el valor de esta ultima se puede obtener a
24
través de los datos compilados por (STEINHAGEN, 1977) y que se muestran en las
siguientes figuras.
Por otro lado, el calor especifico de la madera depende de la temperatura y el

contenido de humedad, pero es, prácticamente independiente de la densidad o tipo de
especie (Pino, Robles, Tepa,...). Al igual que la conductividad térmica STEINHAGEN
(1977) presenta valores compilados en las figuras 8, 9 y 10.
FIGURA N° 9 Factor de corrección

Ak, versus gravedad
especifica. línea1
valido para un
contenido de
humedad < 40%,
línea 2 para > 40%.
FIGURA N° 8 Conductividad térmica en la

dirección radial, versus
temperatura con el
contenido de humedad.
25
FIGURA N° 10 Calor especifico, versus

temperatura, con contenido de
humedad.
2.5.2 Propiedades Térmicas de los Alimentos
Para los alimentos, los valores de las propiedades térmicas se establecen

experimentalmente, pues dependen de muchos factores que son difíciles de predecir. Al
igual que otros materiales (maderas, metales) la conductividad térmica no es constante,
depende de la temperatura. Cuando el rango de temperaturas es limitado, esta
dependencia puede ser despreciada y es posible asumir que las propiedades no varían
con la temperatura.
Valores experimentales de conductividad térmica, calor especifico, densidad para

diferentes alimentos han sido modelados en correlaciones matemáticas. Estos modelos
han sido muy útiles en la estimación de estas propiedades. CHOI y OKOS (1986)
26
desarrollaron y presentaron modelos matemáticos para predecir dichas propiedades y

son las que se muestran en los Cuadros 4 y 5.
El valor de éstas puede ser hallado como una función de los mayores
constituyentes básicos (proteína, grasa, carbohidrato, fibra, minerales y agua)
encontrados en el alimento y de la temperatura.
La conductividad y difusividad térmica dependen de la estructura espacial del

material. Consecuentemente el modelo utiliza una fracción de volumen φ i que estima el
valor de la propiedad para la composición del material y que equivale a:
wi
ρi
φi = (2.6)
w
∑ ρi
i
Donde:
wi : Fracción de masa para cada constituyente.
ρ i : Densidad intrínseca para cada constituyente.
CUADRO N° 4 Modelo matemático para la predicción de las propiedades

termofísicas de alimentos de alto contenido de humedad,
según CHOI y OKOS (1986).
Propiedad Modelo
ρ 1 ∑ (w ρ )
i i
Cp ∑ (C pi wi )
Normal al Flujo de Calor Paralelo al Flujo de Calor
k  φi 
∑ (k φ )i i 1 ∑  k 
 i
α ∑ (α w )
i i
27
CUADRO N° 5 Propiedades térmicas según la mayor componente de los

alimentos como función de la temperatura.
Propiedad Componente Ecuación
Carbohidratos k = 0.20141 + 1.3874 ⋅ 10 −3 T − 4.3312 ⋅ 10 −6 T 2
Minerales k = 0.32962 + 1.4011 ⋅ 10 −3 T − 2.9069 ⋅ 10 −6 T 2
k (W / m° K ) Fibra k = 0.18331 + 1.2497 ⋅ 10 −3 T − 3.1683 ⋅ 10 −6 T 2
Grasa k = 0.18071 + 2.7604 ⋅ 10 −4 T − 1.7749 ⋅ 10 −7 T 2
Proteína k = 0.17881 + 1.1958 ⋅ 10 −3 T − 2.7178 ⋅ 10 −6 T 2
Agua k = 0.57109 + 1.7625 ⋅ 10 −3 T − 6.7036 ⋅ 10 −6 T 2
−3 −6 2
Carbohidratos C p = 1.5488 + 1.9625 ⋅ 10 T − 5.9399 ⋅ 10 T
Minerales C p = 1.0926 + 1.8896 ⋅ 10 −3 T − 3.6817 ⋅ 10 −6 T 2
C p = 1.8459 + 1.8306 ⋅ 10 −3 T − 4.6509 ⋅ 10 −6 T 2
C p (kJ / kg °C ) Fibra
Grasa C p = 1.9842 + 1.4733 ⋅ 10 −3 T − 4.8008 ⋅ 10 −6 T 2
Proteína C p = 2.0082 + 1.2089 ⋅ 10 −3 T − 1.3129 ⋅ 10 −6 T 2
Agua C p = 4.1762 − 9.0864 ⋅ 10 −5 T − 5.4731 ⋅ 10 −6 T 2
Carbohidratos ρ = 1.5991 ⋅ 10 3 − 0.31046 ⋅ T
Minerales ρ = 2.4238 ⋅ 10 3 − 0.28063 ⋅ T
ρ (kg / m 3 ) Fibra ρ = 1.3115 ⋅ 10 3 − 0.36589 ⋅ T
Grasa ρ = 9.2559 ⋅ 10 2 − 0.41757 ⋅ T
Proteína ρ = 1.3299 ⋅ 10 3 − 0.51840 ⋅ T
Agua ρ = 997.18 + 3.1439 ⋅ 10 −3 T − 3.7574 ⋅ 10 −3 T 2
Carbohidratos α = 8.0842 ⋅ 10 −2 + 5.3052 ⋅ 10 −4 T − 2.3218 ⋅ 10 −6 T 2
Minerales α = 1.2461 ⋅ 10 −1 + 3.7321 ⋅ 10 −4 T − 1.2244 ⋅ 10 −6 T 2
α = 7.3976 ⋅ 10 −2 + 5.1902 ⋅ 10 −4 T − 2.2202 ⋅ 10 −6 T 2
α (m 2 / s )× 10 6
Fibra
Grasa α = 9.8777 ⋅ 10 −2 + 1.2569 ⋅ 10 −4 T − 3.8286 ⋅ 10 −8 T 2
Proteína α = 6.8714 ⋅ 10 −2 + 4.7578 ⋅ 10 −4 T − 1.4646 ⋅ 10 −6 T 2
Agua α = 0.1317 + 6.2477 ⋅ 10 −4 T − 2.4022 ⋅ 10 −6 T 2
2.5.3 Propiedades Térmicas de los Metales
La selección de datos confiables de las propiedades es una parte integral de

cualquier análisis cuidadoso de ingeniería, ha de evitarse el uso ocasional de valores de
publicaciones o manuales que no hayan sido bien caracterizados o evaluados
28
(INCOPRENA y DEWITT, 1999). Al igual que los otros materiales ya revisados las
propiedades térmicas de los metales dependen de la composición misma del sólido y por
supuesto de la temperatura , tales valores pueden ser obtenidos de diferentes autores, por
ejemplo la dependencia de la temperatura de las propiedades térmicas (conductividad
térmica, calor especifico) del acero SAE 1045, de densidad igual a 7844 Kg / m 3 , según
SÁNCHEZ (1994) equivalen a:
k (T ) = 5.1572 × 10 2 + 2.7464 × 10 −4 × T − 8.1334 × 10 −5 × T 2 + 5.354 × 10 −4 × T 3

+ 1.9212 × 10 −12 × T 4 + 1.7603 × 10 −19 × T 5 (2.1)
(W / m° K )
Cp (T ) = 3.9552 × 10 2 + 1.1676 × T − 7.0682 × 10 −4 × T 2 − 2.0021× 10 −5 × T 3

+ 6.1869 × 10 −8 × T 4 + 4.582 × 10 −11 × T 5 (2.2)
(J / Kg °K )
Para otros metales INCROPENA y DEWITT (1999) presentan diferentes valores

en tablas, de propiedades en función de la temperatura, por otro lado la pagina web
www.matweb.com provee más de 27500 datos de propiedades de diversos materiales.
29
Capítulo 3 – DESARROLLO DE LOS MODELOS MATEMÁTICOS
3.1 Introducción
Como el modelamiento de los procesos térmicos de calentamiento y/o

enfriamiento es un problema complejo, que depende de factores de forma, tamaño,
propiedades del material y características del medio, se necesita desarrollar un modelo
matemático que permita describir los factores antes señalados, es decir, idealizar lo más
posible la situación real del proceso.
Así entonces, se establecerán primeramente las hipótesis o suposiciones en

relación al modelo geométrico del producto, su comportamiento y propiedades térmicas,
además de las características del medio al cual se exponen los productos. Todo con el fin
de plantear posteriormente las ecuaciones de transferencia de calor, a cada geometría en
particular.
La solución al problema del modelado matemático será resuelto a través del

método numérico de diferencias finitas en su forma explícita.
3.2 Suposiciones del Modelo
Antes de formular y desarrollar los modelos matemáticos de conducción de calor

se establecerán las suposiciones del modelo y que involucran el factor geometría del
producto, propiedades térmicas, comportamiento del material y medio circundante, todo
con el objetivo de acotar el problema complejo de simular un producto en un proceso de
calentamiento y/o enfriamiento.
3.2.1 Suposiciones Geométricas
Se consideran cuatro geometrías que comprenden, formas de procesamiento

común, existentes en diversos procesos térmicos encontrados en la industria alimenticia,
30
maderera, metalmecánica y que son las siguientes: placas, esferas, cilindros, y

paralelogramos.
La conducción de calor para cada una de ellas se supondrá de la siguiente forma:
3.2.1.1 Placas
Se considera una placa extensa, de espesor a , inicialmente a una temperatura

igual y uniforme en toda su distribución espacial, expuesta a un flujo de calor cuya
mayor incidencia está en la dirección normal a la placa, ya que la rapidez con que
alcanza la temperatura de tratamiento en el centro de ella es mucho mayor que en las
otras direcciones, o sea T (z , t ) (ver Figura 11):
FIGURA N° 11 Modelo geométrico para

una placa extensa.
3.2.1.2 Esferas:
Se tiene una esfera expuesta a flujos de calor en la dirección radial,

produciéndose en ella líneas de igual temperaturas llamadas isotermas, por lo que se
considera el flujo de calor en una dirección, es decir T (r , t ) , (ver Figura 12):
31
FIGURA N° 12 Modelo geométrico para

una esfera.
3.2.1.3 Cilindros
Es posible realizar la simulación en dos tipos de cilindros que se diferencian

según la relación existente entre el largo y el radio, es decir se supondrán cilindros
infinitos y finitos. La conducción de calor para los cilindros infinitos se considera
solamente en la dirección radial, T (r , t ) , en tanto que en el caso de los cilindros finitos
se considera radial y axial, T (r , z , t ) (vea la Figura 13), además se supone que
inicialmente el cuerpo está a una temperatura igual y uniforme en todos los puntos.
Qr Qr
r r
R R
Qz
L L
a) Cilindro Infinito b) Cilindro Finito
FIGURA N° 13 Modelo geométrico para un cilindro infinito y finito.

32
3.2.1.4 Paralelogramos:
Se simulará paralelogramos de longitud infinita ( L ≥ 5 Dmedio 2 ) que implica la
conducción de calor en las direcciones x e y, T (x, y, t ) como se ilustra en la Figura 14,

además la temperatura inicial es igual y uniforme en todos los puntos.
FIGURA N° 14 Modelo geométrico para un

paralelogramo.
3.2.2 Suposiciones de las Propiedades Térmicas y del comportamiento del

material
A continuación se realizaran las suposiciones que permitan describir las

propiedades térmicas y el comportamiento de los materiales madereros, alimenticios y
metálicos, frente a los procesos de calentamiento y/o enfriamiento de manera de poder
simular cada uno de ellos en forma, lo más real posible.
3.2.2.1 Maderas
Se considera el material totalmente homogéneo, o sea no existen

discontinuidades dentro de la madera, la humedad es constante en todo el material a
33
estudiar y tiene un mismo valor durante el proceso, además la conductividad térmica y el

calor especifico pueden ser considerados constantes o variables por el usuario. El valor
de estas propiedades en función del contenido de humedad, temperatura y densidad
básica, se obtienen a partir de las siguientes ecuaciones, basadas en los datos compilados
por (STEINHAGEN, 1977).
!"Conductividad Térmica:
kr = [k 0 + mk × (T − 273)]× Ak (W / m° K ) (3.1)
k0 = 3.042 × 10−8 × H 3 − 8.376 × 10 −6 × H 2 + 2.866 × 10 −3 × H + 0.1478 (3.2)
mk = −6.757 × 10 −11 × H 3 + 1.949 × 10−8 × H 2 + 1.184 × 10 −5 × H

(3.3)
+ 2.555 × 10 − 4
Ak = 0.62 + [1.7662 × (GS − 0.3)] (3.4)
!"Calor Especifico:
Cp = [Cp 0 + mCp × (T − 273)] (J / Kg ° K ) (3.5)
Cp0 = 7.849 × 10 −4 × H 3 − 0.2352 × H 2 + 29.55 × H + 1580 (3.6)
mCp = −3.16 × 10 −4 × H 2 + 0.0666 × H + 2.63 (3.7)
Donde:
H : Contenido de humedad de la madera (%)
mk : Pendiente ecuación (3.3) dependiente de la temperatura (W / m° K )
Ak : Factor de corrección en términos de GS
34
Según varios autores citados por (STEINHAGEN, 1977) se considera el

comportamiento anisotropico de la madera, según la relación en la conductividad
térmica K Z (Longitudinal ) = Re lación × K r (Radial ) , que se mantiene constante en todo el
proceso y que el usuario puede adoptar en promedio como (Relación = 2) .
3.2.2.2 Alimentos
El análisis para el caso de los alimentos, es complejo ya que en la realidad

existen variadas estructuras físicas, es decir en estricto rigor las propiedades termofísicas
varían en cualquier dirección, que para el caso del los alimentos están dadas por la
orientación de las fibras de este. Además la conductividad térmica es diferente si el flujo
de calor incide sobre el sólido en forma perpendicular o paralela a esta (VALENTAS et
al., 1997, HELDMAN, 1975).
En este trabajo el alimento puede ser considerado con las propiedades

térmicas constantes ó variables y como un medio isotrópico para los vegetales
(zanahorias, arvejas, espárragos, etc.), las propiedades térmicas variables se determinan
según las correlaciones de CHOI y OKOS (1986), (ver Sección 2.5.2).
3.2.2.3 Metales
Las propiedades térmicas de los metales pueden ser consideradas como

constantes ó variables con la temperatura para algunos metales, para el SAE 1045 las
correlaciones se obtienen a partir del trabajo de SÁNCHEZ (1994) (ver Sección 2.5.3),
en el Anexo 1 se muestran los polinomios de interpolación de las propiedades térmicas
como función de la temperatura, obtenidos de los datos expuestos por INCROPERA y
DEWITT (1999).
Se considera el material esencialmente homogéneo en relación a las propiedades

térmicas, además se supone el comportamiento del material isotrópico, es decir las
propiedades definidas en el estudio son las mismas en todas las direcciones.
35
3.2.3 Suposiciones del Medio de Calentamiento y/o Enfriamiento
Como los productos van a ser sometidos a un sistema térmico, en unidades de

calentamiento y/o enfriamiento, se definirán las características del medio circundante
con el fin de evaluar la transferencia de calor por convección en la interfase sólido –
fluido. Para tales efectos, se realizan ciertas suposiciones que trataran de tomar un
amplio espectro de condiciones de operación posibles durante los procesos.
- La condición en la frontera del sólido será del tipo tercera clase (h, T∞ ) .
- Se supondrá como medio de enfriamiento líquidos o gases.
- Se desprecia el mecanismo de transferencia de calor por radiación en el interior del
sólido. Es decir se supondrá un sólido no transparente.
- El modo de transferencia de calor por convección la definirá el usuario, como
resultado de un gradiente de densidad (convección natural) o un gradiente de presión
(convección forzada).
- El efecto de convección es uniforme en todo el contorno del sólido en cualquier
instante de tiempo.
- El valor del coeficiente de convección puede ser definido por el usuario como
variable o constante en el tiempo.
3.3 Modelo Matemático de Conducción de Calor
A Continuación y para cada geometría y considerando las suposiciones

realizadas, se desarrollarán los modelos matemáticos basados en la ecuación de difusión
térmica, sin generación interna de calor, en coordenadas cartesianas y cilíndricas para el
caso que corresponda.
3.3.1 Modelo Matemático para una Placa Extensa
Se realiza una modelación matemática formulada por la ecuación de conducción

de calor de Fourier con propiedades térmicas variables, ver ecuación (3.1),
36
presentándola posteriormente en coordenadas cartesianas, como se muestra en la Figura

11 de la Sección 3.2.1.1.
∂T
ρC p = ∇[k ( T )∇T ] (3.1)
∂t
La ecuación (3.1) en coordenadas cartesianas 1D T (z ,t ) , se escribe como:
∂T ∂  ∂T 
ρC p =  k( T )  (3.2)
∂t ∂z  ∂Z 
La ecuación diferencial de conducción de calor (3.2), que tiene la conductividad

térmica dependiente de la temperatura, es una ecuación no lineal cuya solución es
compleja, por lo que puede ser resuelta, utilizando el método de aproximación sucesiva,
que consiste en transformar la ecuación no lineal (3.2) en una ecuación lineal integral,
asumiendo la conductividad térmica constante del material, para un periodo de tiempo
anterior, es decir se determinará la temperatura en la primera aproximación de manera
de usarla para el cáculo de la conductividad térmica en la segunda aproximación.
Por lo tanto la ecuación (3.2) queda de la siguiente manera:
∂T ∂ 2T
ρC p = k (T ) 2 (3.3)
∂t ∂Z
Discretizando la ecuación (3.3) en un dominio Ω como se muestra en la Figura

15. Se tiene una placa de espesor = a y se considera por simetría la mitad de la placa en
la discretización.
37
FIGURA N° 15 Discretización modelo placa

extensa
El dominio ( z ,t ) se divide en intervalos ∆z , ∆t tal que:
z = (i − 1)∆z i = 2 ,3 ,...P − 1
t = n∆t n = 0 ,1,2 ,3...
La distribución de temperaturas para cualquier nodo i del dominio Ω y en

cualquier instante n + 1 , se obtiene reemplazando la ecuación diferencial (3.3) que
gobierna el proceso, por una ecuación recursiva, utilizando el método de diferencias
finitas en su forma explicita, teniendo los gradientes de temperaturas presentados en la
Sección 2.4, se escribe como.
 T n +1 − Ti n   T n − 2Ti n + Ti +n1 
ρC p  i  = k (T ) i −1  (3.4)
 ∆t   ∆z 2 
∀ i = 2,3,... P − 1 (0 < z < a / 2)
Si ∆z = l , y luego despejando Ti n +1 :
 2  1 n
Ti n +1 = Ti n  1 −  + (
Ti −1 + Ti +n1 ) (3.5)
 M  M
∀ i = 2,3,... P − 1 (0 < z < a / 2)

38
Donde:
l2 k( T )
M = , con α z =
α z ∆t ρC p
3.3.2 Modelo Matemático para una Esfera y un Cilindro Infinito
A continuación se desarrollará la modelación matemática conjuntamente para

una esfera y un cilindro infinito, que permita describir la distribución de temperaturas en
el interior del sólido. Según la Figura 12 de la Sección 3.2.1.2 y Figura 13 a) de la
Sección 3.2.1.3, la ecuación (3.1) se puede escribir en coordenadas cilíndricas y
esféricas, T (r ,t ) :
∂T ∂  ∂T  b  ∂T 
ρC p =  k (T )  +  k (T )  (3.6)
∂t ∂r  ∂r  r  ∂r 
Donde:
1 → Cilindro
b=
2 → Esfera
La solución de la ecuación (3.6) se efectúa ídem a la ecuación (3.2) del modelo

de conducción de calor de una placa extensa. Es decir la ecuación (3.6) queda de la
siguiente manera:
∂T  ∂ 2T b  ∂T 
ρC p = k (T )  2 +   (3.7)
∂t  ∂r r  ∂r 
Discretizando la ecuación (3.7) en el dominio Ω , para una esfera y un cilindro

infinito de radio R, como se muestra en la Figura 16.
39
∆r
P-1
i+1
i
Ω i-1 R
FIGURA N° 16 Discretización modelo

esfera y cilindro infinito
El dominio ( r ,t ) se divide en intervalos ∆r , ∆t tal que:
r = (i − 1)∆z i = 2 ,3 ,...P − 1
t = n∆t n = 0 ,1,2 ,3...
La distribución de temperaturas para cualquier nodo i del Ω y en cualquier

instante n + 1 , se obtiene reemplazando la ecuación diferencial (3.6) que gobierna el
proceso; por un sistema de ecuaciones algebraicas, utilizando el método de diferencias
finitas en su forma explicita, teniendo los gradientes de temperaturas presentados en la
Sección 2.4, se escribe como:
 Ti n +1 − Ti n   T n − 2Ti n + Ti +n1 b  Ti +n1 − Ti −n1 

ρC p   = k (T )  i −1 +  
(i − 1)∆r  2∆r 
 (3.8)
 ∆t   ∆r 2
∀ i = 2,3,... P − 1 (0 < r < R)
Si ∆r = l , y luego despejando Ti n +1 :
 n n 
 2  1
Ti n +1 = Ti n 1 −  +
b
(
Ti −1 − Ti +1 + 2(i − 1) Ti −1 + Ti +1 
n n
) (3.9)
 M M  
∀ i = 2,3,... P − 1 (0 < r < R)

40
Donde:
l2 k (T )
M= , con α r =
α r ∆t ρC p
3.3.3 Modelo matemático para un Paralelogramo Infinito
Se utiliza una modelación matemática formulada por la ecuación de conducción

de calor de Fourier con propiedades térmicas variables, según la Figura 14 de la Sección
3.2.1.4 la ecuación (3.1) en coordenadas cartesianas T (x , y ,t ) , se escribe como:
∂T ∂  ∂T  ∂  ∂T 
ρC p =  k x (T )  +  k y (T )  (3.10)
∂t ∂x  ∂x  ∂y  ∂y 
La solución de la ecuación (3.10) se efectúa ídem a la ecuación (3.2) del modelo

de conducción de calor de una placa extensa. Es decir la ecuación (3.10) queda de la
siguiente manera:
∂T ∂ 2T ∂ 2T
ρC p = k x (T ) 2 + k y (T ) 2 (3.11)
∂t ∂x ∂y
Discretizando la ecuación (3.11) del dominio Ω de un paralelogramo infinito

como se muestra en la Figura 17, el dominio espacial y temporal ( z ,t ) se divide en
intervalos ∆z , ∆t tal que:
x = (i − 1)∆z i = 2 ,3 ,...P − 1 (0 < x < a/2 )

y = ( j − 1)z j = 2 ,3,...Q − 1 (0 < y < b/2 )
t = n∆t n = 0 ,1,2 ,3...
41
y
∆x
∆y
Ω i,j+1
Q-1
i-1,j i,j i+1,j
b
i,j-1
1,1 P-1 P
x
a
FIGURA N° 17 Discretización modelo paralelogramo infinito
La distribución de temperaturas para cualquier nodo i , j del Ω y en cualquier

proceso; por una ecuación recursiva, utilizando el método de diferencias finitas en su
forma explicita, teniendo los gradientes de temperaturas presentados en la Sección 2.4,
se escribe como:
 Ti n +1 − Ti n   T n − 2Ti ,nj + Ti +n1, j   Ti ,nj −1 − 2Ti ,nj + Ti ,nj +1 

ρC p   = k x (T ) i −1, j  + k y (T )  (3.12)
 ∆t   ∆x 2   ∆y 2 
∀ i = 2 ,3 ,... P − 1 (0 < x < a / 2 )
j = 2,3,...Q - 1 (0 < y < b/2 )
Si ∆x = ∆y = l , y despejando Ti n, j+1 se tiene:
  1
Ti n, j+1 = Ti n, j  1 −
2
−
2
 + (
Ti −n1, j + Ti +n1 +
1
) (
Ti n, j −1 + Ti n, j +1 (3.13) )
 M1 M 2  M1 M2
∀ i = 2 ,3 ,... P − 1 (0 < x < a / 2 )
j = 2,3,...Q - 1 (0 < y < b/2 )
42
Donde:
l2 k (T )
M1 = , con α x =
α x ∆t ρC p
l2 k (T )
M2 = , con α y =
α y ∆t ρC p
3.3.4 Modelo matemático para un Cilindro Finito
Se utiliza una modelación matemática basada en la ecuación de conducción de

calor de Fourier con propiedades térmicas variables, según la Figura 13 b) de la Sección
3.2.1.3 la ecuación (3.1) en coordenadas cilíndricas T (r , z ,t ), se escribe de la siguiente
forma:
∂T 1 ∂  ∂T  ∂  ∂T 
ρC p =  k r (T )r  +  k z (T )  (3.14)
∂t r ∂r  ∂r  ∂z  ∂r 
La solución de la ecuación (3.14) se efectúa igual a la ecuación (3.2) del modelo

de conducción de calor de una placa extensa. Por lo tanto la ecuación (3.14) queda de la
siguiente manera para ser resuelta por el método de diferencias finitas.
∂T  1 ∂  ∂T  ∂T
ρC p = k r (T ) r  + k z (T ) 2 (3.15)
∂t  r ∂r  ∂r  ∂z
Discretizando la ecuación (3.15) del dominio Ω de un cilindro finito como se

muestra en la Figura 18, el dominio espacial y temporal ( r , z ,t ) se divide en intervalos
∆r , ∆z , ∆t tal que:
r = (i − 1)∆r i = 2 ,3 ,...P − 1 ( 0 < r < R )
z = ( j − 1)∆z j = 2 ,3 ,...Q − 1 ( 0 < z < L / 2 )
t = n∆t n = 0 ,1,2 ,3...
43
Z ∆r
∆z
Q-1
i,j+1
i+1,j Ω i,j i+1,j
i,j-1
1,1 P-1 P r
FIGURA N° 18 Discretización modelo cilindro
finito
La distribución de temperaturas para cualquier nodo (i , j ) del Ω y en cualquier

proceso; por una ecuación recursiva, utilizando el método de diferencias finitas en su
forma explicita, teniendo los gradientes de temperaturas presentados en la Sección 2.4,
se escribe entonces como:
 Ti ,nj+1 − Ti ,nj   T n − 2Ti ,nj + Ti +n1, j Ti +n1, j − Ti −1, j 

ρC p   = k r (T ) i −1, j +
1
+

 ∆t 
 
 (∆ r )2
(i − 1 )∆r ( 2 ∆r ) 
(3.16)
 Ti ,nj −1 − 2Ti ,nj + Ti ,nj +1 
+ k z (T )z  
 (∆z )2 
∀ i = 2 ,3,... P − 1 (0 < r < R)
j = 2,3,... Q - 1 (0 < z < L / 2 )
Si ∆r = ∆z = l , y despejando Ti n, j+1 , se tiene:

44
  1  Ti +1, j − Ti −1, j  

n n
2 2  + Ti −n1, j + Ti +n1, j  +
Ti ,nj+1 = Ti ,nj 1 − −  + 
 M1 M 2  M 2  2(i − 1) 
 
(3.17)
+
1
(
Ti ,nj −1 + Ti ,nj +1 )
M2
∀ i = 2 ,3,... P − 1 (0 < r < R)
j = 2,3,... Q - 1 (0 < z < L / 2 )
Donde:
l2 k (T )
M1 = , con α r = r
α r ∆t ρC p
l2 k (T )
M2 = , con α z = z
α z ∆t ρC p
45
3.4 Condiciones: Inicial y de Borde
A continuación se desarrollan las correspondientes ecuaciones algebraicas

basadas en las condiciones limitantes del problema, las cuales consideran
3.4.1 Placa Extensa
• Condición Inicial:
Para calentamiento la condición inicial equivale a:

Ti n = Ti ∀ i n=0
Para enfriamiento la condición inicial equivale a:
Ti n ≅ Ti n ∀ i n=n
Para el caso del Temple se considera como condición inicial en todas las
geometrías:
Ti n ≅ Ti n ∀ i n=0
• Condiciones de Borde:
- Nodo i sobre el borde adiabático:
Se considera un nodo i en el borde aislado que coincide con el eje X , en donde

∂T
se tiene la condición de borde adiabático, es decir = 0 . Como se muestra en la
∂z z =0
Figura 19, existe un nodo imaginario i − 1 que es espejo del nodo i + 1 , lo que implica
una condición de simetría, o sea Ti +n1 = Ti −n1 . Así aplicando esta condición a la ecuación
(3.5) se tiene la expresión que describe la temperatura en un tiempo n + 1 , en el centro

de la placa, y que se escribe de la siguiente manera:
46
 2  2
Ti n+1 = Ti n  1 −  + Ti +n1 (3.18)
 M M
∀i = 1
FIGURA N° 19 Condición de borde placa extensa nodo

central.
- Nodo i sobre el Borde Superficial
Se considera un nodo i en el borde superficial paralelo al eje X , en donde se

supone la transferencia de calor por convección desde un fluido contiguo con
almacenamiento de energía como se muestra en la Figura 20. Se aplica el método de
balance de energía para derivar la ecuación de diferencia finita. Por el principio de
conservación de la energía, se tiene:
Incremento de energía = Energía entrante − Energía saliente

∆Q = Q1 − Q2
∂T ∂T
ρcV = hA(T∞ − Ti ) − K z (T ) A
∂t ∂z
∆z ∂T  Ti n − Ti −n1 
ρc = h(Tb − Ti ) − K z (T ) 

2 ∂t  ∆z 
47
Si ∆z = l y T∞ = Tb (Temperatura de baño / ambiente), entonces se tiene Ti n +1 :
 2 2 ∆th  2 ∆th 2
Ti n+1 = Ti n 1 − −  + Tb + Ti −n1 (3.19)
 M ρcl  ρcl M
∀ i=P
FIGURA N° 20 Condición de borde placa extensa nodo

superficial.
3.4.2 Esfera y Cilindro Infinito
A continuación se desarrollan las correspondientes ecuaciones algebraicas de

diferencias finitas basadas en las condiciones limitantes del problema de conducción de
calor para una esfera y cilindro infinito:

Ti n = Ti ∀ i n=0
Ti n ≅ Ti n ∀ i n=n
48
- Nodo i sobre el centro adiabático de la esfera y el cilindro infinito:
Se considera en r = 0 para una esfera y un cilindro infinito la condición de

∂T
segunda clase donde ( = 0 ) de centro adiabático como se muestra en la Figura 21.
∂r r =0
0
Al introducir esta restricción a la ecuación (3.7), esta queda indeterminada por lo que
0
aplicando L’Hopital’s se tiene:
d  ∂T 
 
1  ∂T  dr  ∂r  ∂ 2T
  = = 2
r  ∂r  r =0 d
(r ) ∂r
dr
Reemplazando en la ecuación (3.7), se escribe:
∂T ∂ 2T
ρC p = k ( T )(1 + b ) 2 ,r = 0 (3.20)
∂t ∂r
Discretizando la ecuación (3.20) anterior, se tiene:

 T n +1 − Ti n   T n − 2Ti n + Ti +n1 
ρC p  i  = k ( T )(1 + b ) i −1  (3.21)
 ∆t   ∆r 2 
Por condición de simetría Ti +n1 = Ti −n1 , además si ∆r = l , y luego despejando Ti n +1 ,
tenemos la ecuación que describe la temperatura en un tiempo n + 1 para una esfera y un

cilindro infinito:
 2(1 + b )  2(1 + b ) n
Ti n+1 = Ti n  1 − + Ti +1 (3.22)
 M  M
∀i = 1
49
∆r
P-1
Ω i+1 R
i
i-1
FIGURA N° 21 Condición de borde esfera y

cilindro infinito – nodo central.
- Nodo i sobre el Borde Superficial de una esfera:
Como se muestra en la Figura 22, el flujo de calor atraviesa el cascaron de la esfera,

por lo que realizando el balance térmico, se tiene:
∆Q = Q1 − Q2
∂T ∂T
ρC pV = hA(T∞ − Ti ) − k r A
∂t ∂r
∆r   Ti n +1 − Ti n  ∆r   Ti n − Ti −n1 
3 2
4 
ρC p π  R −    = h4πR 2 ( Tb − Ti n ) − k r 4π  R −   
3  2   ∆t   2   ∆r 
Despejando Ti n +1 , se tiene:
 24 R 2 h∆t 6α∆t  24 R 2 h∆t

Ti n +1 n 
= Ti 1 − −  + Tb +
 ρC P (2 R − ∆r )
3
∆r (2 R − ∆r )  ρC p (2 R − ∆r )3
(3.23)
6α∆t
+ Ti −n1
∆r ( R − ∆r )
2
∀ i=P
50
Q1 i
/2
∆r
∆r
Q2
i-1
Ω
R
FIGURA N° 22 Condición de borde esfera – nodo

superficial.
- Nodo i sobre un borde superficial de un Cilindro Infinito:
Como se muestra en la Figura 23, el flujo de calor atraviesa el cascaron de un

cilindro de largo infinito, realizando el balance térmico, se tiene:
∆Q = Q1 − Q2
∂T ∂T
ρC pV = hA(T∞ − Ti ) − k r A
∂t ∂r
∆r  ∆r  Ti n − Ti n  ∆r  Ti n − Ti −n1 
+1
 
ρC p 2π  R −    = h 2πR(Tb − Ti ) − k r 2π  R −

n
 
 4  2  ∆t   2  ∆r 
Despejando Ti n +1 , se tiene:
 2(2 R − ∆r )k∆t 4hR∆t  4hR∆t

Ti n +1 = Ti n 1 − −  + Tb
 (4 R − ∆r )∆rρC P ρC P (4 R − ∆r )  ρC P (4 R − ∆r ) (3.24)
2(2 R − ∆r )k∆t n
+ T
(4 R − ∆r )∆rρC P i −1
∀ i=P
51
i
Q1 i-1
∆
r
R Q2
Ω
∆ r/2
FIGURA N° 23 Condición de borde cilindro infinito –

nodo superficial.
3.4.3 Cilindro Finito
A continuación se desarrollan las correspondientes ecuaciones algebraicas

basadas en las condiciones limitantes del problema de conducción de calor en un
cilindro finito:

Ti n, j = Ti ∀ i, j , n=0

Ti ,nj = Ti ,nj ∀ i, j , n=n
- Nodo i,j sobre el Borde adiabático coincidente con el eje Z:
Se considera un nodo i, j que coincide con el eje Z, en donde se tiene la

∂T
condición de borde adiabático, es decir = 0 . Como se muestra en la Figura 24,
∂r r = 0
52
existe un nodo imaginario i − 1, j que es espejo del nodo i + 1, j , lo que implica una
condición de simetría, o sea Ti +n1, j = Ti −n1, j . Así aplicando esta condición a la ecuación
0
(3.15) se genera una relación indeterminada , por lo que se requiere aplicar
0
L’Hospital’s. Entonces la ecuación, queda de la siguiente forma:
∂T ∂T ∂T
ρC p = k r (T ) 2 + k z (T ) 2 (3.25)
∂t ∂r ∂z
Discretizando la ecuación (3.25) se tiene:
 Ti ,nj+1 − Ti ,nj   T n − 2Ti ,nj + Ti +n1, j   Ti ,nj −1 − 2Ti ,nj + Ti ,nj +1 

ρC p   = k r (T ) i −1, j  + k z (T )  (3.26)
 ∆t   (∆ )2
  (∆ )2

   r   z
Luego si ∆r = ∆z = l y despejando Ti ,nj+1 :
  4 n
Ti n, j+1 = Ti n, j  1 −
4
−
2
 + Ti +1, j +
1
[
Ti n, j −1 + Ti n, j +1 ] (3.27)
 M1 M 2  M1 M2
∀ i =1
j = 2,3,...Q - 1 (0 < z < L / 2 )
Donde: M 1 y M 2 se definen en la pagina 42.

53
z ∆r
Q
T/ r = 0
∆z
Q-1
i,j+1
i+1,j
∆z
i,j
i,j-1
∆r/2
1,1 P-1 P
r
FIGURA N° 24 Condición de borde cilindro finito nodo (i,j)

sobre el borde adiabático eje z.
- Nodo i,j sobre el Borde adiabático coincidente con el eje r:
Se considera un nodo i, j que coincide con el eje r, en donde se tiene la

∂T
∂z z =0
existe un nodo imaginario i, j − 1 que es espejo del nodo i, j + 1 , lo que implica una
condición de simetría, o sea Ti ,nj +1 = Ti ,nj −1 . Así aplicar esta condición a la ecuación (3.17)
se genera la ecuación que describe la temperatura en el borde adiabático radial r, queda

se escribe de la siguiente manera:
  1  Ti +1, j − Ti −1, j  

n n
2 2  + Ti −n1, j + Ti +n1, j  + 2 Ti ,nj +1 (3.28)
Ti ,nj+1 = Ti ,nj 1 − −  + 
 M1 M 2  M 1  2(i − 1) 
  M 2
∀ i = 2,3,... P − 1 (0 < r < R)
j =1
54
z ∆r
∆z
Q-1
∆z/2
i,j+1
1,1 i,j i+1,j P

r
i-1,j
T/ z = 0 ∆r

sobre el borde adiabático eje r.
- Nodo i,j en la intersección de los Bordes Adiabáticos:
Se considera un nodo i, j en la intersección de dos bordes adiabáticos como se

ilustra en la Figura 26, donde se considera la existencia de dos nodos imaginarios i, j − 1
y i − 1, j que son espejos de los nodos i, j + 1 y i + 1, j , lo que implica una condición de
simetría, es decir Ti n, j +1 = Ti n, j −1 ; Ti +n1, j = Ti −n1, j y que aplicando estas igualdades a la
ecuación (3.17), se obtiene la ecuación que describe la temperatura en la intersección de

los bordes, la cual se escribe de la siguiente manera:
 2 2  2 n 2 n
Ti ,nj+1 = Ti ,nj 1 − −  + Ti +1, j + Ti , j +1 (3.29)
 M1 M 2  M1 M2
∀ i =1
j =1
55
z ∆r
∆z
Q-1
T/ r = 0
∆z/2
i,j+1
r
i,j i+1,j P-1 P
∆r/2 T/ z = 0
FIGURA N° 26 Condición de borde cilindro finito nodo (1,1)

en la intersección de los dos bordes
adiabáticos.
- Nodo i,j sobre el Borde superficial
Como se muestra en la Figura 27, el flujo de calor atraviesa el cascaron de un

cilindro finito, realizando el balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3 − Q4
∂T ∂T ∂T
ρcV = hA(T∞ − Ti ) − k r (T ) A − K z (T ) A
∂t ∂r ∂z
∆r  ∆r  Ti , j − Ti , j   T n − Ti −n1, j 
n +1
 = h 2πR∆z( Tb − Ti n, j ) − k r (T )2π  R − ∆r ∆z  i , j
n
 
ρc 2π  R − ∆z
 4  2  ∆t 
  2   ∆r 

∆r  ∆r  Ti , j − Ti , j +1 Ti , j − Ti , j −1 
n n n n

− K z (T )2π  R −  +
 4  2  ∆z ∆z 

Si consideramos la malla con ∆r = ∆z = l y posteriormente despejando Ti n, j+1 se
tiene:
56
 4 (2 R − l ) 2 8R∆th  8R∆th
Ti ,nj+1 = Ti ,nj 1 − − − + Tb +
 M 1 (4 R − l ) M 2 ρC p l (4 R − l ) ρC p l (4 R − l )
(3.30)
4 (2 R − l ) n
+ (
Ti −1, j +
1
)
Ti , j +1 + Ti , j −1
n n
M 1 (4 R − l ) M2
∀i = P
j = 2 ,3 ,...Q − 1
z ∆r
∆z
Q-1
i,j+1
Q3
i-1,j
∆z
i,j
Q2 Q1
i,j-1
Q4
∆r/2
1,1 P-1 P
r
FIGURA N° 27 Condición de borde cilindro finito nodo

(i,j) sobre el borde convectivo.
- Nodo i,j sobre el Borde lateral
Como se muestra en la Figura 28, el flujo de calor atraviesa la parte lateral de un

cilindro finito, realizando el balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3 − Q4
∂T ∂T ∂T
ρcV = hA(T∞ − Ti ) − k r (T ) A − K z (T ) A
∂t ∂r ∂z
57
∆z  Ti , j − TI ,J 
n +1 n
ρc 2π (i − 1)∆r ∆r  = h 2π (i − 1)∆r∆r( Tb − Ti n, j ) −
2  ∆t 
 
∆r  ∆z  Ti , j − Ti −1, j  ∆r  ∆z  Ti , j − Ti +1, j 
n n n n
 
− k r ( T )2π  (i − 1)∆r −  − k ( T )2π  (i − 1)∆r +  
 2  2  ∆r  r
  2  2  ∆r 

 Ti n, j − Ti n, j −1 
− K z ( T )2π (i − 1)∆r∆r  
 ∆ z 
 
∆r
∆z/2
z
Q1
i-1,j i,j i+1,j
Q
Q2 Q3
Q-1
Q4 i,j-1
∆z
r
1,1 P-1 P

sobre el borde convectivo.
tiene:
 2 2 2∆th  2∆th 1  (2i − 3) n (2i − 1) T n  +
Ti ,nj+1 = Ti ,nj 1 − − − + Tb +  Ti −1, j +
 M1 M 2 ρcl  ρcl M1  (2i − 2 ) (2i − 2) i +1, j 
(3.31)
2 n
+ Ti , j −1
M2
∀ i = 2 ,3,... P − 1
j =Q
58
- Nodo i,j en la intersección de los bordes lateral y superficial:
Como se muestra en la Figura 29, el nodo i, j se ve afectado por la transferencia

de calor por convección tanto en la parte lateral como superficial de un cilindro finito.
Realizando el balance térmico, se tiene:
∆Q = Q1 + Q2 − Q3 − Q4
∆r/2
z ∆r
∆z/2
Q1
i-1,j i,j
Q
Q3 Q2
Q-1 Q4 i,j-1
∆z
1,1 P-1 P
r
FIGURA N° 29 Condición de borde cilindro finito nodo

(P,Q) intersección bordes convectivos.
∂T ∂T ∂T
ρC pV = hA(T∞ − Ti ) − k r (T ) A − K z (T ) A
∂t ∂r ∂z
∆r  ∆z ∆r  Ti , j − TI ,J 
n +1
∆r  ∆r ∆z
(Tb − Ti n, j ) + h 2πR (Tb − Ti n, j ) −
n
 
ρC p 2π  R −  = h 2π  R − 
 4  2 2  ∆t 
  4  2 2
∆r  ∆z  Ti , j − Ti −1, j  ∆r  ∆r  Ti , j − Ti , j −1 
n n n n
 
− k r (T )2π  R −  − K z (T )2π  R − 
 2  2  ∆r 
  4  2  ∆z 

tiene:
59
 (2 R − l ) 4 2 2 ∆th  4 R  2 ∆th  4R 
Ti n, j+1 = Ti n, j 1 − − −  1 +  +  1 + T +
 (4 R − l ) M 1 M 2 ρcl  (4 R − l )  ρcl  (4 R − l )  b
(3.32)
+
( 2R − l ) 4 n
T +
2 n
T
(4 R − l ) M 1 i −1, j M 2 i , j −1
∀i = P
j =Q
- Nodo i,j en la Intersección de los Bordes Adiabático – Convectivo (P,1):

de calor por convección solamente en la parte superficial de un cilindro finito,
realizando el balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3
∂T ∂T ∂T
∂t ∂r ∂z
∆r  ∆z ∆r  Ti , j − Ti , j 
n +1
∆z
( )
n

ρC p 2π  R −  = h 2πR Tb − Ti n, j −
 4  2 2   ∆t  2

∆r  ∆z  Ti , j − Ti −1, j  ∆r  ∆r  Ti , j − Ti , j +1 
n n n n
 
− k r (T )2π  R −   − K z (T )2π  R − 
 2  2  ∆r   4  2  ∆z 
 
si consideramos la malla con ∆r = ∆z = l y posteriormente despejando Ti n, j+1 se
tiene:
 (2 R − l ) 4 2 8 R∆th  8R∆th
Ti ,nj+1 = Ti ,nj 1 − − − + Tb +
 (4 R − l ) M 1 M 2 ρC p l (4 R − l ) ρC p l (4 R − l )
(3.33)
+
(2 R − l ) 4 T n + 2 T n
(4 R − l ) M 1 i −1, j M 2 i, j +1
∀i = P
j =1
60
z ∆r
∆z
Q-1
∆r/2
Q3
i,j+1
Q2 Q1
r
1,1 i-1,j i,j
FIGURA N° 30 ∆z/2
Condición de borde cilindro finito nodo (P,1)
intersección borde convectivo y borde
adiabático.
- Nodo i,j en la Intersección de los Bordes Adiabático – Convectivo (1,Q):

de calor por convección solamente en la parte lateral de un cilindro finito. Realizando el
balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3
∂T ∂T ∂T
∂t ∂r ∂z
∆r ∆z ∆r  Ti , j − Ti , j  ∆r ∆z  Ti , j − Ti +1, j 
n +1
∆r ∆z
( )
n n n
ρC p 2π = h 2π Tb − Ti , j − k r (T )2π
n −
4 2 2  ∆t 
 4 2 2 2  ∆r 

∆r ∆r  Ti , j − Ti , j −1 
n n
− K z (T )2π
4 2  ∆z 

61
tiene:
 4 2 2∆th  2∆th 4 n 2 n
Ti ,nj+1 = Ti ,nj 1 − − − + Tb + Ti +1, j + Ti , j −1 (3.34)
 M 1 M 2 ρC p l  ρC p l M1 M2
∀ i =1
j=Q
z ∆r
∆z/2
Q1
i,j i-1,j
∆z
Q2
i,j-1 Q3
∆r/2
1,1 P-1 P
r
FIGURA N° 31 Condición de borde cilindro finito nodo (1,Q)

adiabático.
62
3.4.4 Paralelogramo
A continuación se desarrollaran las correspondientes ecuaciones algebraicas

basadas en las condiciones limitantes del problema:

Ti n, j = Ti ∀ i, j , n=0

Ti ,nj = Ti ,nj ∀ i, j , n=n
- Nodo i,j sobre Borde Adiabático coincidente con el eje y:
Se considera un nodo i, j que coincide con el eje y, en donde se tiene la

∂T
∂x x =0
existe un nodo imaginario i − 1, j que es espejo del nodo i + 1, j , lo que implica una
condición de simetría, o sea Ti +n1, j = Ti −n1, j . Así aplicando esta condición a la ecuación
(3.13) se obtiene la ecuación que permite calcular la temperatura en un tiempo n + 1 en

el borde adiabático, la que se escribe de la siguiente forma:
 2  2 n
Ti n, j+1 = Ti n  1 −
2
−  + Ti +1, j +
1
( )
Ti n, j −1 + Ti n, j +1 (3.35)
 M 1 M 2  M 1 M 2
∀i = 1
j = 2 ,3,...,Q − 1 (0 < y < b/2 )
63
y
∆x
T/ x = 0 Q
∆y
i,j+1
i+1,j
i,j
∆y i,j-1
∆x/2
1,1 P-1 P
x
FIGURA N° 32 Condición de borde paralelogramo nodo

(i,j) sobre el borde adiabático eje y.
- Nodo i,j sobre el Borde adiabático coincidente con el eje x:
Se considera un nodo i, j que coincide con el eje x, en donde se tiene la

∂T
∂y y =0
existe un nodo imaginario i, j − 1 que es espejo del nodo i, j + 1 , lo que implica una
condición de simetría, o sea Ti ,nj +1 = Ti ,nj −1 . Así aplicar esta condición a la ecuación (3.13)
se genera la ecuación que describe la temperatura en el borde adiabático x, la que se

escribe de la siguiente manera:
 2  1
Ti n, j+1 = Ti n  1 −
2
−  + ( )
Ti −n1, j + Ti +n1, j +
2 n
Ti , j +1 (3.36)
 M1 M 2  M1 M2
∀ i = 2 ,3 ,..., P − 1 (0 < x < a/2 )
j =1
64
y
∆x
∆y
Q-1
∆y/2
i,j+1
1 P
x
i-1,j i,j i+1,j
T/ y = 0 ∆x

(i,j) sobre el borde adiabático eje x.
- Nodo i,j en la intersección de los Bordes Adiabáticos:
Se considera un nodo i, j en la intersección de dos bordes

∂T ∂T
adiabáticos =0, = 0 como se ilustra en la Figura 34, donde se considera la
∂x x =0 ∂y y =0
existencia de dos nodos imaginarios i, j − 1 y i − 1, j que son espejos de los nodos

i, j + 1 y i + 1, j , lo que implica una condición de simetría, es decir
Ti n, j +1 = Ti ,nj −1 ; Ti +n1, j = Ti −n1, j y que aplicando estas igualdades a la ecuación (3.13), se
obtiene la ecuación que describe la temperatura en la intersección de los bordes, la cual

se escribe de la siguiente manera:
 2 2  2 n 2 n
Ti n, j+1 = Ti n  1 − −  + Ti +1, j + Ti , j +1 (3.37)
 M1 M2  M1 M2
∀i =1
j =1
65
y
∆x
∆y
Q-1
T/ x = 0
i,j+1
i+1,j
i,j P-1 P
x
∆y/2
∆x/2
T/ y = 0
FIGURA N° 34 Condición de borde paralelogramo nodo (1,1) en

la intersección de los dos bordes adiabáticos.
- Nodo i,j sobre el Borde lateral derecho:
Como se muestra en la Figura 35, el flujo de calor atraviesa el área lateral

derecha del paralelogramo infinito. Realizando el balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3 − Q4
∂T ∂T ∂T
ρC pV = hA(T∞ − Ti ) − k x A − ky A
∂t ∂x ∂y
∆x  Ti , j − Ti , j  T n −T n 
n +1
( )
n
ρC p ∆y  = h∆y Tb − Ti n, j − k x ∆y i , j i −1, j  −

2  ∆t 


 ∆x 

∆x  Ti , j − Ti , j +1 Ti , j − Ti , j −1 
n n n n
− ky +
2  ∆y ∆y 

Si consideramos una malla con ∆x = ∆y = l , posteriormente despejamos Ti n, j+1 , se
tiene:
66
 2 ∆th  2 ∆th
Ti n, j+1 = Ti n  1 −

2
−
2
−
ρ
+
 ρC l
Tb +
2 n
Ti −1, j +
1
(
Ti n, j +1 + Ti n, j −1 (3.38) )
 M 1 M 2 C l
p  p M 1 M 2
∀i = P
j = 2 ,3 ,...Q − 1
y
∆x ∆x/2
∆y
Q-1
Q3 i,j+1
i+1,j
∆y
i,j
Q2 Q1
Q4 i,j-1
1,1 P-1 P
x

(i,j) sobre el borde Convectivo.
- Nodo i,j sobre el Borde lateral superior:
Como se muestra en la Figura 36, el flujo de calor atraviesa el área lateral

superior del paralelogramo infinito. Realizando el balance térmico, se tiene:
∆Q = Q1 − Q2 − Q3 − Q4
∂T ∂T ∂T
∂t ∂x ∂y
∆y  Ti , j − Ti , j   T n − Ti +n1, j Ti n, j − Ti −n1, j 
n +1
 = h∆x Tb − Ti n, j − k x ∆y  i , j
( )
n
ρC p ∆x + −
2  ∆t 
 2 
 ∆ x ∆x 

 Ti n, j − Ti n, j −1 
− k y ∆x 
 ∆ y 
 
67
tiene:
 2 ∆th  2 ∆th
Ti n, j+1 = Ti n  1 −

2
−
2
−
ρ
+
 ρC l
Tb +
1
( )
Ti +n1, j + Ti −n1, j +
2 n
Ti , j −1 (3.39)
 M 1 M 2 C p l  p M 1 M 1
∀ i = 2 ,3 ,..., P − 1
j =Q
y ∆x
∆y/2
Q1
i-1,j i,j i+1,j
Q
Q3 Q4
Q-1
i,j-1 Q2
∆y
1,1 P-1 P
x

(i,j) sobre el borde convectivo.
- Nodo i,j en la intersección de los bordes laterales:

de calor por convección tanto en la parte lateral superior como derecha. Realizando el
balance térmico, se tiene:
∆Q = Q1 + Q2 − Q3 − Q4
∂T ∂T ∂T
∂t ∂x ∂y
68
∆x ∆y  Ti , j − Ti , j   T n − Ti −n1, j 
n +1
 = h ∆x Tb − Ti n, j + h ∆y Tb − Ti n, j − k x ∆y  i , j
( ) ( )
n
ρC p −
2 2  ∆t 
 2 2 2  ∆x 

 Ti n, j − Ti n, j −1 
− k y ∆x 
 ∆ y 
 
tiene:
 2 2 4 ∆th  4 ∆th 2 n 2 n
Ti n, j+1 = Ti n  1 − − − + Tb + Ti −1, j + Ti , j −1 (3.40)
 
M 1 M 2 ρC p l  ρC p l M1 M2

∀i = P
j =Q
∆x/2
y
∆x
∆y/2
Q1
i-1,j i,j
Q
Q3 Q2
Q-1
Q4 i,j-1
∆y
1,1 P-1 P
x

(P,Q) intersección bordes adiabáticos.
- Nodo i,j en la intersección de los Bordes Adiabático – Convectivo (P,1):

de calor por convección solamente en la parte lateral derecha. Realizando el balance
térmico, se tiene:
69
y
∆x
∆y
Q-1
∆x/2
Q3 i,j+1
Q2 Q1
1,1
x
i-1,j i,j
∆y/2
FIGURA N° 38 Condición de borde paralelogramo nodo (P,1)
adiabático.
∆Q = Q1 − Q2 − Q3
∂T ∂T ∂T
ρC pV = hA(T∞ − T ) − k x A − ky A
∂t ∂x ∂y
∆x ∆y  Ti , j − Ti , j   T n − Ti −n1, j 
n +1
 = h ∆y Tb − Ti n, j − k x ∆y  i , j
( )
n
ρC p −
2 2  ∆t 
 2 2  ∆x 

∆x  Ti , j − Ti , j +1 
n n
− ky
2  ∆y 

tiene:
 2 2 2 ∆th  2 ∆th 2 n 2 n
Ti n, j+1 = Ti n  1 − − − + Tb + Ti −1, j + Ti , j +1 (3.41)
 M M ρ C l  ρC l M M
 1 2 p  p 1 2
∀i = P
j =1
70
- Nodo i,j en la intersección de los Borde Adiabático - Convectivo (1,Q):

de calor por convección solamente en la parte superior. Realizando el balance térmico,
se tiene:
∆Q = Q1 − Q2 − Q3
∂T ∂T ∂T
∂t ∂x ∂y
∆x ∆y  Ti , j − Ti , j   T n − Ti +n1, j 
n +1
 = h ∆y Tb − Ti n, j − k x ∆y  i , j
( )
n
ρC p −
2 2  ∆t 
 2 2  ∆x 

∆x  Ti , j − Ti , j −1 
n n
− ky
2  ∆y 

tiene:
 2 2 2∆th  2∆th 2 n 2 n
Ti n, j+1 = Ti n  1 − − − + Tb + Ti +1, j + Ti , j −1 (3.42)
 
M 1 M 2 ρC p l  ρC p l M1 M2

∀i =1
j=Q
71
∆y/2
Q1
i,j i-1,j
∆y
Q2
i,j-1 Q3
∆x/2
1,1 P-1 P
x

(1,Q) intersección borde convectivo y
borde adiabático.
3.5 Energía Almacenada /Removida al/del Producto
La energía almacenada / removida de un sólido cualquiera puede ser evaluada

como:
Q = ρ × Cp × V × ∆T (3.43)
Donde V es el volumen del sólido (medio continuo), pero esta ecuación no es

valida al enfoque discreto de la solución que se implementa en el trabajo, por lo tanto la
ecuación (3.43) queda de la siguiente manera (para un modelo lineal y bi-direccional ):
( P,Q )
Q(t ) = ∑V × ρ × Cp × (T (t ) − T )
( i =1, j =1)
i i i (3.44)
Donde Vi es el volumen de control para cada nodo en la discretización y puede

ser obtenido de las ecuaciones desarrollas anteriormente.
72
3.6 Condiciones de Estabilidad
3.6.1 Introducción
Al utilizar el método numérico de diferencias finitas en su forma explicita, a la

solución de problemas de conducción de calor transitorios, se necesita asegurar la
estabilidad de la solución, ya que es una de las característica de este método, la que se
explica en detalle en el Capítulo 2, Sección 2.4.
Para los problemas que trata el presente trabajo, la condición de estabilidad se

determina agrupando el cuociente asociado al nodo de interés en el tiempo anterior y que
esté sea mayor o igual a cero. En general esto se realiza reuniendo todos los términos
que incluyen Ti n o Ti ,nj , según sea el caso.
Como se requiere que el ∆t se mantenga en ciertos limites, este depende del

tamaño de la malla y otros parámetros como las propiedades térmicas y el coeficiente de
transferencia de calor por convección, tales valores se consideran constantes ó variables
en el tiempo lo que implica encontrar la o las expresión(es) de ∆t en la situación más
desfavorable o restrictiva. A continuación se presentan dichas ecuaciones, para los
diferentes modelos geométricos.
3.6.2 Placa Extensa
De la ecuación 3.5 se obtiene el primer cuociente que acompaña a Ti n , para una

placa extensa.
Donde:
 2 
0 ≤ 1 − 
 M
Despejando ∆t , se tiene:
73
l2
∆t ≤ (3.45)
2α
La otra situación desfavorable se obtiene de la ecuación 3.19, por lo que se tiene

lo siguiente:
 2α∆t 2∆th 
0 ≤ 1 − 2 −
 l ρC p l 

Despejando ∆t , se obtiene:
1
∆t ≤
 2α 2h  (3.46)
 2 + 
l ρC P l 
3.6.3 Esfera y Cilindro Infinito
De la ecuación 3.22 se obtiene el primer cuociente que acompaña a Ti n , tanto

para la esfera como para cilindro infinito.
Donde:
 2(1 + b ) 
0 ≤ 1 − 
 M 
l2
∆t ≤ (3.47)
2α (1 + b )
De las ecuaciones 3.23 y 3.24 se obtienen respectivamente las segundas

expresiones posiblemente desfavorables en el criterio de estabilidad de una:
- Esfera:
 24 R 2 h∆t 6α∆t 

0 ≤ 1− − 

 ρC P (2 R − ∆ r )3
∆r (2 R − ∆r ) 

74
1
∆t ≤
 24 R 2 h 6α r  (3.48)
 + 
 ρC (2 R − ∆r )
 P
3
∆r (2 R − ∆r ) 

- Cilindro Infinito
 4(2 R − ∆r )α r ∆t 8 Rh∆t 
0 ≤ 1 − − 
 (4 R − ∆r )∆r 2
ρC P l (4 R − ∆r ) 
1
∆t ≤
 4(2 R − ∆r )α r 8 Rh  (3.49)
 + 
 (4 R − ∆r )∆r
2
ρC P l (4 R − ∆r ) 
3.6.4 Cilindro Finito
De la ecuación 3.32 se obtiene el cuociente que acompaña a Ti ,nj , para un
cilindro finito.
  4( 2 R − l )α r ∆t 2α z ∆t 2h(8 R − l )∆t  
0 ≤ 1 −  + +  

  ( 4 R − l )l 2
l 2
ρC P l ( 4 R − l ) 
1
∆t ≤
 2( 2 R − l )α r α z h(8 R − l )  (3.50)
2 + 2 + 
 (4 R − l )l
2
l ρC P l (4 R − l ) 
Además de la ecuación 3.34 como otra ecuación desfavorable se obtiene el

cuociente que acompaña a Ti ,nj , para un cilindro finito.
75
 4α ∆t 2α z ∆t 2∆th 
0 ≤ 1 − r2 − −
 l l2 ρC p l 

1
∆t ≤
 2α α h  (3.51)
2 2 r + 2z +
 l l ρC p l 

3.6.5 Paralelogramo
De las ecuaciones 3.13, 3.35, 3.36 o 3.37, se obtiene el cuociente que acompaña a
Ti ,nj , para un paralelogramo.
 2 2 
0 ≤ 1 − − 
 M 1 M 2 
l2
∆t ≤ (3.52)
2(α r + α z )
También de la ecuación 3.40 como otra ecuación desfavorable se obtiene el

cuociente que acompaña a Ti ,nj , para un paralelogramo.
 2α ∆t 2α ∆t 4∆th 
0 ≤ 1 − r2 − r2 −
 l l ρC p l 

1
∆t ≤
α + α 2h  (3.53)
2 r 2 z +
 l ρC p l 

76
Capítulo 4 – IMPLEMENTACIÓN DEL SOFTWARE
4.1 Introducción
El presente capítulo tiene por objetivo entregar los antecedentes de la

implementación de los modelos matemáticos desarrollados en el Capítulo 3 en un
programa computacional, elaborado con la finalidad de describir adecuadamente la
cinética de los procesos de conducción de calor transiente. Se darán las pautas
necesarias para la adecuada comprensión del funcionamiento y alcance de este,
empezando por el algoritmo de programación, utilizado para la posterior conformación
de la interfaz de usuario y rutinas de cálculo.
4.2 Lenguaje de Programación
Los modelos desarrollados, así como las correlaciones matemáticas para la

predicción de las propiedades térmicas de los productos, fueron implementados en un
programa computacional, escrito en Visual Basic versión 6.0, ejecutable en cualquier
plataforma de operación Windows®. La herramienta utilizada esta basada en objetos
(Visual) con propiedades y métodos, los cuales se usan para crear la interfaz gráfica de
usuario (GUI). En lugar de escribir numerosas líneas de código para describir la
apariencia y la ubicación de los elementos de la interfaz, simplemente se pueden
arrastrar y colocar objetos prefabricados en su lugar dentro de la pantalla.
Por otro lado, La palabra “Basic” hace referencia al lenguaje BASIC (Beginners
All-Purpose Symbolic Instruction Code), un lenguaje original creado en 1963,
utilizado por más programadores que ningún otro lenguaje en la historia de la
informática o computación. Visual Basic ha evolucionado a partir del lenguaje BASIC
original y ahora contiene centenares de instrucciones, funciones y palabras clave,
muchas de las cuales están directamente relacionadas con la interfaz gráfica de
Windows.
77
4.3 Algoritmo de Programación
Antes de elaborar cualquier aplicación computacional se hace recomendable

elaborar un diagrama de flujo o algoritmo de programación, con el fin de establecer los
lineamientos en la conformación de la interfaz de usuario y de la posterior escritura de
los códigos de programación. Al respecto, en las figuras siguientes se muestran los
diagramas correspondientes al software desarrollado y que tiene por nombre “Procesos
CCT”.
Inicio
Selección de Estudio
Estudios en Maderas Estudios en Alimentos Estudios en Metales
Selección del Modelo Geométrico: Selección del Modelo Geométrico:

Selección del Modelo Geométrico:
Placa Placa
Placa
Cilindro Infinito Cilindro Infinito
Cilindro Infinito
Cilindro Finito Cilindro Finito
Cilindro Finito
Esfera Esfera
Paralelogramo
Paralelogramo Paralelogramo
Ingreso de Datos Fijos:

Dimensiones Geometría Seleccionada [m]
Tamaño de la Malla según sea la Geometría
Temperatura Inicial (Ti, [°K])
Condición de Termino (Tdes, [°K])
Selección del Tipo de Proceso:
Calentamiento y/o Enfriamiento
FIGURA N° 40 Diagrama de flujo “Procesos CCT” (a).

78
Temperatura de
Baño/Ambiental Si
Variable
No Ingreso de:
Temp. Inicial (Tbi, [°K])
Incremento de Ti en tiempo (dt, [s])
Temp. final de baño (Tbf, [°K])
Ingreso de:
Temperatura de Baño/Ambiental
Constante (Tb, [°K])
Propiedades Térmicas
Si No
Constantes del sólido
Ingrese los Datos, de las

Propiedades Térmicas del sólido
(r [Kg/m3], k [W/m°K], [Cp J/Kg])
Coeficiente de Transferencia
No de Calor por Convección
2
(hm,[W/m °K]) según Usuario
Si
No de Calor por Convección
2
Natural (hm,[W/m °K])
Si
de Calor por Convección
Natural (hm,[W/m2°K])
2 3
FIGURA N° 41 Diagrama de flujo “Procesos CCT” (b).

79
2 3
Propiedades del
Si
Fluido Constantes
Ingresar:
No Coeficiente (hm, [W/m2°K])
Ingresar propiedades
constantes del fluido:
Número de Prandtl (Pr)
Conductividad Térmica (k,[W/m°K])
Viscosidad Cinemática (n, [m2/s])
Coef.de Expansión Volum. (b, [1/°K])
Calcular Intervalos de Tiempo (dt1)

como condición de estabilidad
en el Cálculo del Proceso
Ingresar Intervalo de Tiempo

(dt, [seg]) menoro igual a los
calculados
Inicio del Cálculo
Se inicia el tiempo
t = n= 0
Se asigna la temperatura inicial

T ni,j = Ti
(para todos los nodos del dominio
de la geometría seleccionada)
FIGURA N° 42 Diagrama de flujo “Procesos CCT” (c).

80
Se incrementa el tiempo
t = t + dt
Se calcula la nueva distribución

de temperatura con las ecuaciones
correspondientes a cada punto de
la geometria seleccionada
T n+1 i,j = f(h, k, Cp, T ni,j,...)
Cálculo de todos los parámetros
considerados variables en el proceso
es decir: Propiedades Térmicas y/o
Coeficiente de Transferencia de Calor. Actualización de la distribución
de temperaturas en n = 0
T ni,j =T n+1i ,j
Condición de Término
n
No T i,j <Tdeseada (Calentami ento)
T ni,j >T deseada (Enfriamient o)
Si
Termino
FIGURA N° 43 Diagrama de flujo “Procesos CCT” (d).
4.4 Interfaz de Usuario
La interfaz es simple y tiene la misma apariencia que diversas aplicaciones de

Windows, el estilo utilizado para este software es el de interfaz de documentos múltiples
(MDI), como se muestra en la Figura 43, contiene: menús de operación y una barra de
herramientas.
81
FIGURA N° 44 Ventana principal “Procesos CCT”
Si se hace un Clic en el botón “Continuar” se despliega la ventana o el módulo

de “Estudios en Maderas”
FIGURA N° 45 Ventana de “Estudios en Maderas”
Como se aprecia en la Figura 44, la ventana que muestra a la izquierda tres

fichas (Características del Producto, Variables del Proceso y Condiciones de Termino y
Estabilidad) destinadas a la configuración de todos los parámetros que interviene en el
calculo del proceso a estudiar. A la derecha se pueden apreciar dos imágenes (Modelo
82
Geométrico y Discretización) las cuales ayudan al usuario (Alumnos) a entender mejor

la modelación a desarrollar.
Para los datos reutilizables, como las propiedades térmicas de los productos se
dispone de ventanas que incluyen un control Data o control de origen de Datos, el cual
permite insertar, presentar y manipular datos de los registros (Campo) de controles
enlazados como el de una caja de texto (TextBox), así también es posible desplazarse de
un registro a otro (ver Figura 46).
FIGURA N° 46 Ventana “Propiedades

Térmicas” en maderas
Terminada la configuración de los parámetros, se procede a realizar el cálculo de

condición de estabilidad en la ficha “Condiciones de Termino y Estabilidad” haciendo
un Clic en el botón “Calcular dt”, se generan dos valores en donde el usuario debe
realizar una comparación para ingresar como último dato un valor menor a ellos. Así
finalmente se procede a realizar la simulación o el cálculo del proceso, en el menú
“Calcular” o en el botón de signo igual.
83
Inmediatamente aparece una ventana llamada “Calculando el Proceso” que

muestra la evolución del proceso (ver Figura 46) en términos de la temperatura vs
tiempo en diferentes puntos.
FIGURA N° 47 Ventana “Calculando el Proceso”
Finalizado el cálculo el usuario puede acceder a las graficas para realizar el

análisis pertinente, haciendo un clic en el menú “Resultados”.
FIGURA N° 48 Ventana “Resultados”

84
4.5 Especificaciones y Alcance del Software
El software “Procesos CCT” elaborado para el presente trabajo de titulación,

tiene los siguientes requerimientos de sistema:
CUADRO N° 6 Requerimientos de “Procesos CCT”
Disco Duro Se requiere de un mínimo de 100 MB

Memoria RAM 32 MB mínimo
Sistema
Windows 95 o mayor
Operativo
CPU AMD, Pentium 266 MHz
Desde el punto de vista del alcance del programa, este puede realizar lo
siguiente:
!"Análisis de Transferencia de Calor, aplicado a productos de geometría simple de

procesamiento común, como son: Placas, Cilindros, Esferas y Paralelogramos,
en estado transiente o inestable.
!"Trabajar con materiales anisotrópicos o isotrópicos, considerando sus

propiedades térmicas Constantes o Variables.
!"Trabajar en la interfaz sólido – fluido con condiciones de Convección Natural o

Forzada, además de poder establecer la temperatura de baño o ambiental
variable en el tiempo.
!"Determinar los perfiles de temperatura versus tiempos en los procesos de

calentamiento y/o enfriamiento en diferentes puntos del sólido, además del calor
absorbido/removido al/del producto.
85
!"Presenta reportes gráficos de la tendencia de la temperatura en el tiempo, calor

almacenado y del coeficiente de transferencia de calor por convección en el
tiempo.
!"Genera un archivo de Microsoft Excel en la dirección “C:\Modelación.xls”, para

el posterior análisis detallado de los perfiles de temperatura, calor, propiedades
térmicas y coeficiente de transferencia de calor por convección.
86
Capítulo 5 – COMPARACIÓN Y DISCUSIÓN DE LOS RESULTADOS
5.1 Introducción
El presente Capítulo tiene por objetivo analizar y discutir la comparación de los

resultados obtenidos por el programa computacional elaborado versus los encontrados
en la literatura y/o los experimentales de investigaciones de la Universidad Austral de
Chile.
Por otro lado se busca generar confiabilidad tanto en los modelos desarrollados
en el Capítulo 3, como en la escritura de los códigos de programación, es decir la
implementación del software.
Para ello se abordan por separado las comparaciones en los estudios realizados
en maderas, alimentos y metales.
5.2 Comparación de Resultados en Maderas
Las comparaciones en esta área se realizan, tomando como base los resultados
encontrados, en primer lugar, en el proceso de macerado de trozas de madera en agua
caliente, y en segundo lugar el proceso de calentamiento de maderas de sección circular
y rectangular, sometidas a vapor.
5.2.1 Resultados en el Proceso de Macerado
La comparación de este proceso se lleva a cabo con base un análisis teórico del
proceso de macerado, que fue estudiado por MORENO y DEVLIEGER (1990). En tal
trabajo se expone una modelación lineal de conducción de calor transiente para predecir
tiempos de calentamiento de rollizos de madera hasta una temperatura deseada para el
debobinado.
87
Un posterior trabajo de MORENO et al (1992), consistió en simular un modelo

para el proceso macerado, incorporando la evaluación del consumo de energía de la
troza. Se realiza una validación con resultados experimentales a escala 1:1, llevada a
cabo en la industria, es decir, en condiciones no optimas como son las de un laboratorio.
En base a esto, se empieza a simular el calentamiento de trozas de maderas a

través de un modelo unidireccional de conducción de calor radial, estableciendo las
hipótesis del trabajo de MORENO y DEVLIEGER (1990), que consisten en considerar
las propiedades térmicas de la madera constante para valores medios de temperatura y
contenido de humedad de las trozas durante el proceso. La temperatura inicial de la troza
es uniforme y superior a 0°C, y temperatura de baño constante.
Como podemos apreciar en los siguientes cuadros, se presentan tiempos de

macerado para diferentes temperatura inicial y contenidos de humedad de la troza con
las respectivas comparaciones entre el modelo desarrollado por MORENO y
DEVLIEGER (1990) y el software elaborado en el presente trabajo.
Se considera un Coeficiente Convectivo h = 250 W/m2°K, y una malla de 40
nodos.
CUADRO N° 7 Tiempos de macerado para diferentes temperaturas iniciales.

Características Moreno (1990) Procesos CCT
GS = 0,429 (Pino) Ti (°C) Tiempo (h) Tiempo (h) Desviación %
D = 0,4 m 5 35,0 34,9 0,3
H = 60 % 10 33,4 33,3 0,0
Tb = 326°K (53 °C) 15 31,7 31,6 0,3
Td = 321°K (48 °C) 20 29,7 29,6 0,3
CUADRO N° 8 Tiempos de macerado para diferentes contenidos de humedad.

Características Moreno (1990) Procesos CCT
GS = 0,429 (Pino) H% Tiempo (h) Tiempo (h) Desviación %
D = 0,4 m 60 33,4 33,3 0,3
Ti = 283°K (10 °C) 80 34,9 34,5 1,1
Tb = 326°K (53 °C) 100 35,9 35,8 0,3
Td = 321°K (48 °C) 120 37,5 36,8 2,0
88
Las discrepancias obtenidas entre los dos trabajos no son significativas, no

superan el 2 % (solo caso de Humedad = 120%), las cuales se atribuyen básicamente a
las diferencias en la obtención de las ecuaciones que permiten calcular las propiedades
térmicas de la madera.
La siguiente Figura muestra un campo de temperaturas típico de calentamiento

de una troza según el modelo de conducción de calor radial.
Proceso de Macerado de una Troza

Temperatura de Baño
333 Nodo Central
323
Temperatura (°K)
Nodo Rollizo Central

313 Nodo Interno a 0.1 m del centro
303 Nodo a 0.15 m del centro

Nodo Superficial
293
283 Características:
GS = 0,429 (Pino), H = 60%
273 Ti = 283 °K
0 5 10 15 20 25 30 35 40 Tb = 326 °K
Td = 321°K (rollizo central)
Tiempo (h) Numero de Nodos = 40
FIGURA N° 49 Distribución de temperatura de una troza según conducción radial.
Para validar los resultados obtenidos por el software “Procesos CCT”, se recurre
a los resultados de la validación experimental realizada por MORENO et al (1992). La
Figura 50 ilustra la temperatura en el rollizo residual provista por el modelo matemático
versus los datos tomados desde la industria para la siguientes especies (Tepa, Alamao,
Ulmo, Olivillo y Roble).
89
FIGURA N° 50 Comparaciones de temperaturas teóricas y reales en una

(Tepa, Alamo, Ulmo, Olivillo y Roble), por Moreno (1992).
Como se puede observar existen desviaciones, las que se atribuyen a las

diferencias entre el contenido de humedad real y el asumido por el modelo. El contenido
de humedad es un parámetro altamente variable desde una troza a otra y también desde
un punto a otro en el interior de ella MORENO et al (1992).
Sin embargo, se debe tener presente que las temperaturas optimas para el proceso
de corte varían en un amplio rango según LUTZ, J. F. (1972), las cuales, se muestran en
el Cuadro 9. En consecuencia se puede decir que las diferencias entre los modelos no
son groseras a nivel industrial.
CUADRO N° 9 Temperaturas optimas para debobinado,

según LUTZ, J.F. (1972).
Temperatura de
Especie GS
debobinado
Tepa 0,418 310 – 326
Alamo 0,331 281 – 300
Ulmo 0,537 340 – 356
Olivillo 0,448 318 – 334
Roble 0,492 330 – 346
Pino 0,429 315 - 331
90
El estudio de la influencia del coeficiente de transferencia de calor por

convección, considerado variable, se muestra en la Figura 51, que comprende
modelaciones obtenidas por el software y cuyas características son las siguientes.
- Diámetro de la Troza : 0,4 m

- Modelo Geométrico : Cilindro Infinito
- Temperatura inicial (Ti) : 283 °K
- Temperatura de Baño Constante (Tb) : 5 °K mayor a la Temperatura deseada
- Temperatura deseada (Td) : 321°K en el rollizo central (0.05 m del centro de
la troza).
- Propiedades térmicas constantes a una humedad H = 60%
Efecto del Coeficiente Convectivo Variable (h)

330
325
320
315
Temperatura (°K)
310
305
Modelación a h constante (250 W/m²°K) - nodo (rollizo central)
300
Modelación a h variable - nodo (rollizo central)
295
Moreno (1990) h constante
290
Modelación a h constante (250 W/m²°K) - nodo superficial
285 Modelación a h variable - nodo superficial
280
0 5 10 15 20 25 30 35 40
Tiempo (h)
FIGURA N° 51 Efecto del Coeficiente convectivo variable en los tiempos de

macerado.
Los tiempos obtenidos para cada modelación, fue el siguiente:
- Modelación a h variable : 33 hrs 38 min

- Modelación a h constante : 33 hrs 26 min
91
La influencia de un coeficiente convectivo variable es prácticamente nulo, las

curvas prácticamente se superponen en todo el proceso. Si se compara la diferencia en el
tiempo de 33,63 hrs conseguido por el pino considerando el coeficiente convectivo
variable y el tiempo de 33,3 hrs con una modelación a coeficiente constante, está es
considerada despreciable, existe un aumento del tiempo de macerado en un 0,6% al
modelar con un coeficiente variable.
La explicación a este hecho, se da simplemente por los altos coeficientes

convectivos entre la troza y el baño (Bi > 100), es decir la resistencia térmica superficial
es despreciable lo que implica que la temperatura superficial prácticamente sea igual a
la del medio en pocos instantes, tal situación se observa en la Figura 49. Se puede
plantear además en base a este fenómeno que es posible utilizar una condición de borde
de primera clase (la temperatura superficial = temperatura de baño en todo el proceso).
La evolución del coeficiente de transferencia de calor por convección en el
tiempo puede ser vista en el software, tal como se ve en la Figura 52, obtenida de la
ecuación de MCAdams (1974) (el Anexo 2 muestra otras ecuaciones para el coeficiente
convectivo).
 gβ (Tb − Tr = R ) 
0.25
hD
Nu = 0.53 Pr 
k  ν2  (5.1)
Valida para el rango 10 4 < Gr Pr < 10 9
500
400
h (W/m²°K)
300 Coeficiente Convectivo h(t)
200
100
0
0 5 10 15 20 25 30 35 40
Tiempo (h)
FIGURA N° 52 Variación del coeficiente convectivo en el tiempo.

92
La influencia de las propiedades térmicas de la madera, también puede ser

estudiada realizando diferentes modelaciones, a propiedades constantes y variables. los
resultados se describen como sigue, en el Cuadro 10:
CUADRO N° 10 Resultados de modelaciones con diferentes especies.

“Procesos CCT”
Moreno (1990) Modelación 1 Modelación 2 Modelación 3
Especie Tiempo (h) Tiempo (h) Tiempo (h) Tiempo (h)
Ulmo 38,9 38,3 37,9 38,5
Roble 38,4 38,0 37,6 38,2
Tepa 37,6 37,5 37,1 37,7
Alamo 36,8 36,6 36,2 36,8
Donde las modelaciones 1,2 y 3 consideran lo siguiente:

!" Modelación 1: propiedades térmicas de la madera constantes y propiedades
termofísicas del fluido variables.
!" Modelación 2: propiedades térmicas de la madera variables y propiedades
termofísicas del fluido variables.
!" Modelación 3: propiedades térmicas de la madera constantes y propiedades
termofísicas del fluido constantes.
Las desviaciones entre los modelos 1 y 2 que consideran las propiedades

térmicas de la madera constantes y variables respectivamente no son grandes. Se
determina una disminución en un 1% en el tiempo de macerado a propiedades variables,
lo cual se esperaba ya que la difusividad térmica de la madera aumenta con la
temperatura (ver Figuras 8, 9 y 10 sección 2.5.1), lo que implica una razón de cambio de
la temperatura en la troza mayor que al considerar propiedades térmicas constantes, la
difusividad térmica es menor en tal caso ya que se evalúa a una temperatura media
(media aritmética entre la temperatura inicial y la deseada).
Por lo tanto es posible decir que no es justificable modelar el problema a

propiedades variables, a razón de otros autores además, (MORENO et al , 1992;
SIMPSON, 2001).
93
Considerar las propiedades termofísicas del agua variables resulta no

influyente, las diferencias pueden despreciarse como se ve en le Cuadro 10. Al comparar
los tiempos, la modelación 3 considera la variación de estos parámetros y los tiempos
obtenidos son mayores en un 0,5% a la modelación 1 de igual condición pero que
considera constante tales propiedades. La explicación a tal aumento está dada por un
disminución en la viscosidad cinemática y aumento en el coeficiente de expansión
térmica del agua.
Pero desde el punto de vista práctico la incorporación de esta alternativa de

modelar a propiedades termofísicas del medio variables, resulta más rápido para el
usuario, al no consultar en tablas los datos disponibles en la literatura.
La influencia de la anisotropía de la madera, es estudiada a través del consumo

de energía, según los modelos bi-direccional (cilindro finito) y unidireccional (cilindro
infinito) versus la relación largo / diámetro de la troza, para tales efectos se considera la
siguiente relación de anisotropía α z = 2α r (recomendada por STEINHAGEN (1977)).
Calor Absorbido vs Relación (L/D)
45
Q (10^3 KJ)
Modelo Bidireccional GS = 0,418 (Tepa)

40 (Cilindro Finito) Ti = 283 °K, H = 60%
35 Tb = 323 °K
30
Modelo Unidireccional
(Cilindro Infinito) Td = 318 °K
25 L = 1,2 m
20
Anisotropía kz/kr = 2
15
10
5
0
1 2 3 4 5 6
Relación L/D
FIGURA N° 53 Energía almacenada vs L/D, L es fijo.
Las diferencias entre los dos modelos se observan en el siguiente Cuadro. Se

comete un error del -0,4% al emplear un modelo unidireccional para una relación
94
L arg o / Diametro ≥ 2 , por lo tanto se concluye que no es justificable utilizar el modelo

bi-direccional para una relación L arg o / Diametro ≥ 2 si la relación de anisotropía
equivale a (kz/kr =2).
CUADRO N° 11 Diferencias entre los modelos bi-direccional y

unidireccional.
Cilindro Finito Cilindro Infinito
L/D Q (10^3 KJ) Q (10^3 KJ) Desviación (%)
1,5 37,181 38,997 -4,88
2 21,806 21,902 -0,44
3 9,735 9,714 0,22
4 5,454 5,433 0,39
5 3,462 3,446 0,46
Otra condición de operación que el software es capaz de resolver en el proceso

de macerado, es el calentamiento a temperatura de baño variable. La Figura 54 ilustra el
efecto en la cantidad de energía almacenada por la troza de madera cuando éste se
expone a condiciones de baños de agua caliente a temperatura: constante y variable en
el tiempo.
Energía Almacenada vs Tiempo

30 333
25 323
Temperatura (°K)
20 313
Q (10^3 KJ)
Calor Almacenado a Tb Variable

15 Calor Almacenado a Tb Constante 303
Temperatura de baño (Tb) Variable
10 Tepa (0,418) H = 100% 293
D = 0,4 m Temperatura (Tb) Constante
Ti = 283 °K
5 Td = 318 °K (rollizo central) 283
0 273
0 5 10 15 20 25 30 35 40
Tiempo (h)
FIGURA N° 54 Evolución de la energia almacenada en una troza de madera en un

baño de agua caliente a temperatura variable y constante.
95
Se aprecia en la figura, la curva de calor almacenado a temperatura de baño

constante se comporta como creciente exponencial, en cambio con una condición de
temperatura de baño variable la curva de energía almacenada presenta un
comportamiento lineal hasta el cambio de régimen (tb constante) y enseguida toma un
rápido crecimiento, terminando como una exponencial.
Como se ha podido comprobar los resultados obtenidos por el software son

bastante aceptables para ser usados a nivel práctico.
Es posible afirmar que no es justificable modelar con un coeficiente convectivo

de transferencia de calor troza – agua variable. Es más, resulta perjudicial al aumentar el
tiempo de cálculo, inicialmente el coeficiente convectivo variable alcanza su valor
máximo (ver Figura 52) lo que provoca una disminución en el intervalo de tiempo ∆t a
incrementar en el cálculo (ver ecuación 3.49).
Por otro lado la influencia de las propiedades térmicas variables generan una
disminución de 1% en el tiempo de calentamiento, lo que es despreciable en la practica.
Sin embargo, queda planteada la interrogante de determinar cual es el mejor modelo que
se ajusta a resultados experimentales. Como se explico anteriormente la validación
experimental de MORENO et al (1992) fue realizada en condiciones desfavorables, se
sugiere entonces realizar un estudio experimental a nivel de laboratorio.
5.2.2 Resultados del Proceso de Calentamiento de Madera con Vapor
Para dar mayor confiabilidad en los resultados obtenidos por el software, se

realiza otra comparación en base al trabajo desarrollado por SIMPSON (2001), quien
elaboró tablas con tiempos de calentamiento para diferentes especies de madera de
sección circular y rectangular, sometidas a vapor.
96
En dicho trabajo, la solución de las ecuaciones de conducción de calor fue

realizada por MACLEAN (1930,1931) y refinadas por Ingersoll y Zobel (1948),
utilizando el método analítico de separación de variables.
Según SIMPSON (2001), para una buena solución a través de aquellas

ecuaciones, se requiere que la superficie de la madera este fría inicialmente, para asumir
inmediatamente que la temperatura de la superficie sea igual y se mantenga a la
temperatura del medio al ingresar la madera a la cámara de vapor. Es decir se usa la
condición de borde de primera clase.
En base a lo expuesto, se presupone que se genera, al momento de ingresar la

madera fría, una película de vapor condensante en la superficie, consiguiendo así altos
coeficientes de transferencia de calor por convección (según KREITH y BOHN (2001)
entre 6.000 – 120.000 W/m²°K) por lo que la resistencia de la superficie es despreciable
(Bi > 100) implicando así que la temperatura de la superficie consiga el estado de
equilibrio inmediatamente a la temperatura del medio.
Así entonces es posible realizar diferentes modelaciones en el software elaborado

(Procesos CCT), considerando los datos extraídos del trabajo de SIMPSON (2001) y
expuestos en el Cuadro 12, a una temperatura ambiente de 358 °K (85°C) y un
coeficiente convectivo de solo 350 W/m²°K estimado en base al numero de Biot (Bi >
100).
CUADRO N° 12 Datos y resultados entre “Procesos CCT” vs SIMPSON (2001)

para maderas de sección rectangular (paralelogramos).
Datos
Especies Ancho (cm) Alto (cm) Humedad % GS Temp inicial (°K)
1 Aspen (Alamo) 13,208 13,97 85 0,411 287,44
2 Aspen 8,128 8,89 95 0,432 289,11
3 Maple (Arce) 15,24 17,78 71 0,543 280,78
4 Maple 7,366 10,922 81 0,527 288,56
97
Resultados
Desviación
Tiempo Tiempo Desviación
Tiempo Simpson (2001)
Calc. (min) (min) Exp. vs
Especies Exp. vs Procesos
Simpson Procesos Procesos CCT
(min) CCT
(2001) CCT (%)
(%)
1 Aspen (Alamo) 133 146 146 -9,8 0,0
2 Aspen 46,6 59,3 56,9 -22,1 4,0
3 Maple (Arce) 214 222 221 -3,3 -2,4
4 Maple 60 58,3 59,7 0,5 0,5
Promedio Absoluto 8,7 0,5
Los tiempos de calentamiento obtenidos por el software son bastantes

satisfactorios en comparación a los tiempos conseguidos por SIMPSON (2001) con el
uso del método de separación de variables, la desviación promedio equivale a un 0,5 %,
se puede concluir de esto que el método de diferencias (FDM) ofrece una excelente
alternativa para resolver estos problemas, sin caer en errores significativos.
Sin embargo, existen diferencias en relación a los tiempos experimentales, pero

que a juicio de SIMPSON (2001) son adecuados para aplicaciones prácticas. Las
posibles explicaciones a las desviaciones observadas, se dan al probable error de no
ubicar la térmocupla en el centro geométrico, otro origen de error fueron los instantes
perdidos en ingresar las muestras a la cámara de vapor, produciéndose variaciones en las
condiciones de operación. En todo caso el autor plantea un grado de exactitud al estimar
los tiempos de calentamiento dentro de un 5% a un 15%.
5.2.3 Análisis del Enfriamiento de una Placa Inclinada
En la industria de contrachapados existen procesos de climatización que tienen

por objetivo enfriar los tableros de madera, como un proceso final para el posterior
almacenamiento de éstos. Generalmente esto se realiza como se muestra en la Figura
55, donde los tableros son enfriados a diferentes ángulos de inclinación. Se plantea el
análisis de la influencia que dicha inclinación tiene en el tiempo de enfriamiento.
98
FIGURA N° 55 Enfriamiento de una placa

inclinada.
Así entonces, se modela una placa con las siguientes características:
- Especie : Tepa
- Dimensiones : 0,04 m (espesor), 1,2 m (altura), 1,2 m (ancho).
- Numero de Nodos: 10
- Temperatura Inicial: 423 °K
- Temperatura Deseada Nodo 1 (centro): 294 °K
- Temperatura de Baño / Ambiental: 293 °K
- Se considera el mecanismo de convección natural, siendo el coeficiente de
transferencia de calor variable con temperatura.
Los tiempos de enfriamiento conseguidos al modelar a diferentes posiciones

(ángulos de inclinación) puede ser vista en la Figura 56. Los ángulos son medidos desde
la vertical en dirección horario. Se observa que el ángulo de inclinación no es influyente
hasta un ángulo de 55°. Aumentar el ángulo de inclinación hasta una posición horizontal
produce un aumento en el tiempo de enfriamiento.
La explicación a este hecho se da por una reducción en la aceleración de gravedad

( g cos θ ) lo que provoca una disminución en la fuerza de empuje, encargada de inducir
99
el movimiento del fluido, que para el caso de una placa vertical actúa exclusivamente
en la dirección ascendente (o descendente – calentamiento de una placa). Si la placa esta
inclinada con respecto a la gravedad, la fuerza de empuje tiene una componente normal
y paralela a la superficie, hay una reducción en las velocidades del fluido a lo largo de la
placa lo que implica una reducción en la transferencia de calor por convección.
Influencia del Angulo Inclinación
25,00
20,00
Tiempo (h)
15,00
Tiempo vs Inclinación de una Placa
10,00
5,00
0,00
0 15 30 45 60 75 90
Angulo (Grados)
FIGURA N° 56 Influencia del ángulo de inclinación en el enfriamiento de una placa.
La reducción del coeficiente convectivo variable en el tiempo a medida que

aumenta en el ángulo de inclinación, se ilustra en la Figura 57, donde se aprecia una
reducción del coeficiente de transferencia de calor por convección, a partir de la curva
que consideran un ángulo de 60°. La correlación utilizada e incorporada en el software
es obtenida de Fujji e Imura (extraida de KREITH y BOHN (2001)).
N u L = 0.56(GrL Pr cos θ )1 / 4 (5.2)
10 5 < GrL Pr cos θ < 1011

Valida para 0 ≤ θ ≤ 89°
100
Coeficiente convectivo variable (h) vs Tiempo
0° de Inclinación
5,5
15° de Inclinación
5
4,5 30° de Inclinación
4 45° de Inclinación
3,5
h (W/m²K)
3
2,5
2 85° de Inclinación
1,5
1
0,5
0
0 5 10 15 20 25
Tiempo (h)
FIGURA N° 57 Variación del coeficiente convectivo en el tiempo a diferentes

ángulos de inclinación.
101
5.3 Comparación de Resultados en Alimentos
El modelo de conducción de calor desarrollado, para el caso de los alimentos, es

comparado con los datos experimentales disponibles en el trabajo de titulación de
ROJAS (2001), para simular los procesos de escaldado (calentamiento) y enfriado de
alimentos vegetales de geometrías regulares.
Se presentan a continuación modelaciones para dos geometrías: Cilindro Infinito

y Esfera, en el Cuadro 13 se describen las características del producto y variables del
proceso.
CUADRO N° 13 Características de los Productos Espárrago y Arveja.

Cilindro Infinito Esfera
Alimento Espárrago Arveja
Diámetro (mm) 12 9
Densidad (Kg/m³) 1019,9 1080,8
Conductividad Térmica (W/m°K) 0,583 0,5118
Calor Especifico (J/Kg°K) 3994,5 3587,8
Temperatura Inicial (°K) 285,8 285,8
Calentamiento:
Temperatura de Calentamiento °K 368,5 368
Coeficiente Convectivo h (W/m²°K) 603 641
Temperatura Deseada Superficie °K 367,8 367,5
Numero de Biot (Bi) 6,2 5,5
Enfriamiento:
Temperatura de Enfriamiento °K 275,1 274,9
Coeficiente Convectivo h (W/m²°K) 332 581
Temperatura Deseada Superficie °K 275,2
Temperatura Deseada Centro °K 275
Numero de Biot (Bi) 3,4 5,1
Como se puede apreciar en las siguientes gráficas, la tendencia de las curvas

corresponden en forma aproximada a las curvas que entregan los datos experimentales.
102
Procesos de Escaldado e Hidroenfriado en

"Espárragos"
373
363
353
Temperatura (°K)
343 Experimental 1
333 Experimental 2
323 Experimental 3
313 Experimental 4
303 Procesos CCT
293
283
273
0 100 200 300 400 500 600
Tiempo (s)
FIGURA N ° 58 Perfil temperatura vs tiempo nodo central – cilindro infinito.
Los errores relativos entre la curva obtenida por el software “Procesos CCT” y
los datos experimentales se presentan en detalle en el siguiente cuadro.
CUADRO N° 14 Resultados de “Procesos CCT" vs datos experimentales espárrago.

Tiempo (s) Error Exp. 1 % Error Exp. 2 % Error Exp. 3 % Error Exp. 4 % Error Medio %
0 -0,2 0,1 -0,2 0,0 -0,1
32 -1,2 0,4 -0,1 1,8 0,2
64 -1,7 -1,4 -1,7 0,2 -1,2
96 -1,7 -1,5 -1,7 -1,1 -1,5
128 -1,7 -1,5 -1,8 -1,4 -1,6
160 -1,3 -1,5 -1,3 -1,6 -1,4
180 -1,1 -1,2 -1,6 -1,0 -1,2
192 -0,8 -0,8 -1,6 -1,4 -1,2
224 2,4 1,2 -1,3 0,9 0,8
256 1,1 0,7 -0,3 1,1 0,6
288 1,0 0,8 0,7 1,7 1,0
320 1,1 1,3 1,3 1,9 1,4
352 1,3 1,6 1,7 2,1 1,7
384 1,3 1,6 1,8 2,0 1,7
416 1,1 1,7 1,7 1,8 1,6
448 1,2 1,5 1,5 1,8 1,5
480 1,0 1,5 1,3 1,6 1,3
512 0,8 1,2 1,0 1,5 1,1
540 0,7 1,0 0,8 1,3 0,9
Promedio Absoluto 1,2
103
Para el caso de la Esfera, se obtiene el siguiente grafico.
Proceso de Escaldado e Hidroenfriado en

"Arvejas"
373
Experimental 1
363
Experimental 2
353
Experimental 3
Temperatura (°K)
343
333 Experimental 4
323 Procesos CCT

313
303
293
283
273
0 50 100 150 200 250 300
Tiempo (s)
FIGURA N ° 59 Perfil temperatura vs tiempo nodo central – esfera.
Al igual que la geometría anterior se tienen los siguientes errores.
CUADRO N° 15 Resultados de “Procesos CCT" vs datos experimentales arveja.

Tiempo (s) Error Exp. 1 % Error Exp. 2 % Error Exp. 3 % Error Exp. 4 % Error Medio %
0 -0,2 0,1 -0,2 0,0 -0,1
32 -1,2 0,4 -0,1 1,8 0,2
64 -1,7 -1,4 -1,7 0,2 -1,2
96 -1,7 -1,5 -1,7 -1,1 -1,5
128 -1,7 -1,5 -1,8 -1,4 -1,6
160 -1,3 -1,5 -1,3 -1,6 -1,4
180 -1,1 -1,2 -1,6 -1,0 -1,2
192 -0,8 -0,8 -1,6 -1,4 -1,2
224 2,4 1,2 -1,3 0,9 0,8
256 1,1 0,7 -0,3 1,1 0,6
288 1,0 0,8 0,7 1,7 1,0
320 1,1 1,3 1,3 1,9 1,4
352 1,3 1,6 1,7 2,1 1,7
384 1,3 1,6 1,8 2,0 1,7
416 1,1 1,7 1,7 1,8 1,6
448 1,2 1,5 1,5 1,8 1,5
480 1,0 1,5 1,3 1,6 1,3
512 0,8 1,2 1,0 1,5 1,1
540 0,7 1,0 0,8 1,3 0,9
Promedio Absoluto 1,21
104
Como se pude apreciar en ambos casos las diferencias disminuyen tendiendo a

cero, a medida que el tiempo del proceso aumenta hasta alcanzar las condiciones de
equilibrio. Los errores no son mayores que un 1,5 % en ambos casos.
5.3.1 Influencia de las Propiedades Térmicas Variables en Alimentos
Para establecer el grado de influencia o si es necesario modelar los procesos de

escaldado y enfriado de alimentos con propiedades térmicas variables, se comparan
modelaciones con las características citadas en el Cuadro 13.
La Figura 60, ilustra las diferencias obtenidas al modelar a propiedades térmicas

variables (P.T.V.) y a propiedades constantes (P.T.C.).
Proceso de Escaldado y Enfriado a

Diferentes Modelaciones en "Esparrágos"
373 Experimental
363 Experimental
353 Experimental
Temperatura (°K)
343
Experimental
333
Modelación (P.T.C.)
323
313 Modelación (P.T.V.)
303
293
283
273
0 100 200 300 400 500 600 700
Tiempo (s)
FIGURA N° 60 Curvas a diferentes modelaciones vs resultados experimentales

en “Espárragos” (0,1<Bi<40).
Como se ve la curva a propiedades variables consigue la tendencia de las otras

curvas, las desviaciones son similares a las entregadas en el Cuadro 14, el error
promedio absoluto obtenido en todo el proceso equivale a 1,2 % igual al conseguido al
modelar a propiedades constantes.
105
El siguiente grafico muestra las curvas para el caso de arvejas, se logra apreciar
que los resultados a propiedades variables no son tan satisfactorios. El error promedio
absoluto equivale a un 2,1 %, la desviación probablemente se deba a la mala
determinación de los constituyes del alimento, parámetros que intervienen sin duda en el
cálculo. Se plantea la necesidad entonces de conocer muy bien tales propiedades.
Proceso de Escaldado y Enfriado a

Diferentes Modelaciones en "Arvejas"
373
363 Experimental
353 Experimental
343 Experimental
Temperatura (°K)
Experimental
333
Modelación (P.T .C.)
323
Modelación (P.T .V.)
313
303
293
283
273
0 50 100 150 200 250 300
Tiempo (s)
FIGURA N° 61 Curvas a diferentes modelaciones vs resultados experimentales

en “Arvejas” (0,1<Bi<40).
La explicación a la poca influencia beneficiosa de tal consideración (propiedades

térmicas variables), se debe a los gradientes de temperaturas pequeños los cuales no
alcanzan a alterar en forma significativa tales parámetros.
En base a estos antecedentes, se concluye que para lograr una mayor precisión en
los resultados no es justificable utilizar modelos donde intervengan las propiedades
térmicas variables de alimentos vegetales. Sin embargo, la inclusión de esta alternativa
de modelación provista por el software, genera posibilidades de realizar otros estudios
en otros alimentos, con solo ingresar sus constituyentes básicos.
106
5.3.2 Influencia del Coeficiente Convectivo en el proceso de Escaldado y

Enfriado de Alimentos Vegetales.
El mecanismo de transferencia de calor por convección en el proceso de

escaldado y enfriado de alimentos vegetales, corresponde al de convección forzada, en la
mayoría de los casos como se menciono en la Sección 2.2.3. Existe una velocidad de
flujo del fluido adyacente al sólido o del sólido al fluido y tal parámetro es altamente
influyente en la determinación de coeficiente convectivo y por ende en los tiempos de
escaldado y enfriado, por otro lado permite diseñar y evaluar los equipos a utilizar,
específicamente determinar en forma optima la rapidez y el largo de la cinta
transportadora.
En base a estos antecedentes, el software incorpora el ingreso de los parámetros

que intervienen en la determinación del coeficiente convectivo.
El Cuadro 16, muestra los coeficientes convectivos obtenidos por ORDÓÑEZ

(1996) en forma experimental, para un cilindro infinito tanto para el calentamiento como
para el enfriamiento de productos vegetales. los datos son comparados con la correlación
de Churchill y Berstein(*) incluida en el software .
CUADRO N° 16 Valores del coeficiente convectivo (h) a diferentes

condiciones de operación vs correlación de Churchill.
Coeficientes convectivos (h)
Velocidad experimentales (W/m²°K) (h) según
Diámetro Temperatura
de flujo Correlación
(mm) medio (°K) Churchill (*)
(m/s) 1 2 3 Promedio
10 353 0,22 4191 4185 3815 4064 3772

16 353 0,18 3328 3450 3641 3473 2096
22 353 0,11 2849 2693 2780 2774 1804
10 293 0,11 2707 2751 2751 2736 2627
16 293 0,20 3691 3343 3192 3409 2877
22 293 0,11 2308 2189 2188 2240 1803
(*) Obtenida de INCROPENA y DEWITT (1999).

107
El efecto de errores medios cometidos al determinar el coeficiente convectivo a

través de la correlación, se puede observar en el siguiente Cuadro. Las desviaciones en
los tiempos no son tan grandes, si consideramos los errores cometidos en la
determinación del coeficiente convectivo en forma experimental.
CUADRO N° 17 Efecto en los tiempos de calentamiento y enfriamiento a

diferentes coeficientes convectivos (experimental / literatura).
Diámetro Temperatura h (W/m²°K) h (W/m²°K) Tiempo (s) Tiempo (s) Desviación
(mm) medio (°K) Experimental Correlación Experimental Correlación (*) (%)
10 353 4064 3772 105,0 105,6 -0,57
16 353 3473 2096 265,2 273,0 -2,94
22 353 2774 1804 498,6 510,6 -2,41
10 293 2736 2627 99,0 100,0 -1,01
16 293 3409 2877 243,0 245,0 -0,82
22 293 2240 1803 462,0 468,0 -1,30
108
5.4 Comparación de Resultados en Metales
En esta área, la Universidad Austral de Chile, no cuenta con estudios realizados a

nivel experimental, por lo que la comparación se lleva a cabo primeramente en base a
los resultados numérico – experimental de MORAGA y SÁNCHEZ (1994) del proceso
del temple en agua. Como segundo punto se comparan los resultados obtenidos por el
método gráfico y el desarrollado en el presente trabajo.
5.4.1 Resultados en el Proceso del Temple
Los autores antes mencionados, modelaron el proceso del temple en agua para
una placa de acero SAE 1045, considerando el cambio de las propiedades térmicas con
la temperatura. Las ecuaciones están escritas en la Sección 2.5.3.
Además se definieron condiciones de borde convectivas que consideran la

variación del coeficiente convectivo en función del tiempo y una condición inicial que
establece una temperatura con una distribución espacial uniforme. Sobre el coeficiente
convectivo variable, los autores presentan la ecuación que logra representar el fenómeno
de enfriamiento en las tres fases características (Conducción y radiación a través de la
capa gaseosa, transporte de vapor, conducción y convección del liquido). La ecuación
obtenida del trabajo mencionado, para fines de incorporarla en software, es la siguiente:
h (T ) = 1 ,294 × 10 -12 × T 6
- 3 ,275 × 10 - 9 × T 5
+ 3 ,095 × 10 - 6 × T 4
- 1 ,323 × 10 - 3 × T 3
+ 0 , 227 × T 2
- 1 ,087306 × T + 1 903,592 (5.3)
T( ° K) , h (T )( W / m ° K )
2
La modelación de este proceso que a continuación se desarrolla, tiene los

siguientes atributos:
- Placa de Acero SAE 1045, 6,3 mm de espesor
- Temperatura inicial: 1050 °K
- Temperatura deseada nodo central: 293.3 °K
- Temperatura del agua: 293 °K
109
Proceso del Temple (etapa de enfriamiento) en una

Placa de Acero SAE 1045
1173
1073
Temperatura (°K) 973
873
Moraga y Sanchez (1994)
773
Propiedades térmicas variables (P.T.V.)
673
Propiedades térmicas constantes (P.T.C.)
573
473
373
273
0 5 10 15 20 25 30 35
Tiempo (s)
FIGURA N° 62 Enfriamiento de una placa delgada, nodo central.
Los resultados obtenidos por el software considerando las propiedades térmicas

del acero variables son mejores que los modelados a propiedades térmicas constantes.
Las propiedades térmicas constantes fueron evaluadas a un valor medio (350°K) lo que
implica una difusividad térmica mayor al valor que verdaderamente tiene tal parámetro
en los primeros 7 segundos desde (1050°K a 350°K), sin embargo ambos convergen al
final del proceso cuando se llega a la condición de término o de equilibrio, con una
diferencia respecto al tiempo experimental de 3 segundos.
Los errores cometidos entre los datos experimentales y los conseguidos a

propiedades variables, pueden deberse a varias situaciones, la primera de ellas es que el
coeficiente convectivo utilizado para esta modelación, fue extraído a través de la lectura
de un gráfico disponible en el trabajo de MORAGA y SÁNCHEZ (1994), lo que implica
ya una incertidumbre en el valor real que éste debe tener.
El siguiente cuadro muestra las diferencias en porcentajes, entre los modelos

simulados y los datos experimentales de MORAGA y SÁNCHEZ (1994).
110
CUADRO N° 18 Errores en los resultados del temple.

Error con (P.T.C.) Error con (P.T.V.)
Tiempo (s)
(%) (%)
0 0 0
2 29,18 5,56
4 35,74 11,43
6 28,05 13,82
8 24,52 16,97
10 20,15 15,84
12 16,32 13,82
14 16,28 14,90
16 14,63 13,88
18 12,44 12,05
20 11,95 11,77
22 10,72 10,65
24 8,77 8,75
26 7,21 7,23
28 3,05 3,07
30 -0,18 -0,15
Promedio Absoluto 14,9 9,9
El siguiente gráfico, muestra claramente la inconveniencia de considerar un valor

medio como coeficiente convectivo constante.
Proceso del Temple (etapa enfriamiento) en una

Placa de Acero SAE 1045
1173
1073
Temperatura (°K)
973
873
773 Moraga y Sanchez (1994)
673 Modelación h Constante (2500 W/m²°K)
573
Modelación (h Variable)
473
373
273
0 10 20 30 40 50
Tiempo (s)
FIGURA N° 63 Enfriamiento de una placa a diferentes modelaciones

(h constante y h variable).
111
5.4.2 Resultados Método Gráfico versus Diferencias Finitas (MDF)
La segunda comparación en esta área se realiza con el interés de poder observar

las diferencias entre el método gráfico y el método numérico de diferencias finitas
empleado en el software. Para ello se plantea un ejercicio obtenido de la pagina 143 de
la publicación “Principios de Transferencia de Calor” por KREITH y BOHN (2001).
El ejercicio dice lo siguiente:
“En un proceso de fabricación, se fabrican piezas de acero en caliente y luego se

enfrían en agua. Considere un cilindro de acero de 2 m de largo y 0.2 m de diámetro (k =
40 W/m°K, α = 1.0 x 10-5 m²/s), con una temperatura inicial de 673 °K, que de pronto
se enfría en agua a 323 °K. Si el coeficiente de transferencia de calor es de 200 W/m²°K,
calcule lo siguiente 20 minutos después de la inmersión.
1.- la temperatura en el centro

2.- la temperatura superficial
3.- el calor transferido al agua durante los 20 minutos iniciales”
El desarrollo del presente ejercicio se muestra en el Anexo 3. Los resultados son

los siguientes, según método gráfico:
1.- temperatura en el centro = 445.5 °K

2.- temperatura superficial = 421 °K
3.- calor transferido al agua Q(t) = 14.9 kWhr
En base a estos datos se realiza la simulación, considerando como condición de

término la temperatura en el nodo central igual al resultado obtenido por el método
gráfico y que es igual a 445.5 °K, se calcula entonces el tiempo para alcanzar dicho
valor.
112
Los resultados obtenidos entre ambos métodos, se muestran en el Cuadro 19.
CUADRO N° 19 Comparación método gráfico versus diferencias finitas.

MDF
Resultados Método Gráfico
(5 Nodos) (10 Nodos) (15 Nodos)
Tiempo (s) 1200 1020 1154 1202
Temperatura en el centro (°K) 445,50 445,18 445,43 445,46
Temperatura superficial (°K) 421,00 424,34 422,05 421,30
Calor transferido al agua (kWhr) 14,90 11,31 13,66 14,60
Como vemos se logra comprobar la efectividad del método de diferencias finitas,

utilizado para resolver los problemas transitorios. Al aumentar el numero de nodos
(malla más fina) la precisión mejora considerablemente.
Como se ha comprobado, los modelos desarrollados e implementados en la

herramienta computacional, ofrecen resolver una variedad de problemas en maderas,
alimentos y metales, sin errores significativos, si estos incluyen los efectos de todas las
variables que afectan al tiempo de calentamiento y calor absorbido/removido al/del
sólido.
110
Capítulo 6 – CONCLUSIONES
Desde el punto de vista de los objetivos planteados, se logró dar cumplimiento a

cada uno de ellos. El software elaborado con fines académicos, posibilita realizar
análisis de conducción de calor transiente lineal y no lineal en las áreas de la industria
maderera, alimenticia y metalmecánica, además con mecanismos de convención natural
o forzada en la frontera del sólido.
El método numérico de diferencias finitas (MDF) en su forma explicita ha

permitido modelar en forma exitosa el comportamiento cinético que ocurre en los
procesos térmicos de calentamiento y/o enfriamiento.
Se validaron los resultados obtenidos por el software versus los encontrados en la

literatura y los teóricos – experimentales obtenidos en estudios previos de la Universidad
Austral de Chile, consiguiendo las siguientes afirmaciones:
Las discrepancias obtenidas entre el software y los trabajos teóricos y

experimentales de macerado no son significativas, no superan el 2 %, las cuales se
atribuyen básicamente a las diferencias en la obtención de las ecuaciones que permiten
calcular las propiedades térmicas de la madera.
La diferencia en los tiempos de macerado, considerando el coeficiente

convectivo variable y los con una modelación a coeficiente convectivo constante es
despreciable, debido a la situación de un Bi > 100. Por tanto, se puede plantear el uso
de una condición de borde de primera clase.
Elevados coeficientes convectivos, generados al ingresar madera fría a una

cámara de vapor, permiten determinar tiempos de calentamiento, también con una
incertidumbre aceptable en la practica.
111
Las desviaciones entre los modelos que consideran las propiedades térmicas de
la madera constantes y variables respectivamente no son considerables y se estima una
disminución en un 1% en el tiempo de macerado a propiedades variables. Por lo tanto es
posible decir que no es justificable modelar a propiedades variables a razón de otros
autores además (MORENO et al , 1992; SIMPSON, 2001).
Considerar las propiedades termofísicas del agua variables resulta no influyente,

se obtienen diferencias de un 0,5% al comparar los tiempos, que consideran constantes
tales propiedades.
Se comete un error del -0,4% al emplear un modelo unidireccional, por lo tanto

se concluye que no es justificable utilizar el modelo bi-direccional para una relación
L arg o / Diámetro ≥ 2 si la relación de anisotropía equivale a kz / kr = 2 .
Los errores cometidos al modelar alimentos vegetales no son mayores que un 1,5
% con respecto a resultados experimentales, lo que se considera aceptable. Para lograr
una mayor precisión en los resultados no es justificable utilizar modelos donde
intervengan las propiedades térmicas variables de alimentos vegetales. Sin embargo la
inclusión de esta alternativa de modelación provista por el software, genera
posibilidades de realizar otros estudios en otros alimentos, con solo ingresar sus
constituyentes básicos de este.
Se comparó el proceso del temple con diferentes modelaciones, siendo la que

más se ajusta a los valores experimentales, la que modela a propiedades térmicas
variables y con un coeficiente de transferencia de calor por convección variable.
El desempeño del método de diferencias finitas explicito, resultó bastante

eficiente en comparación a otros métodos (analítico, método grafico). Por lo que ofrece
resolver una variedad de problemas en maderas, alimentos y metales, sin errores
significativos, si estos incluyen los efectos de todas las variables que afectan al tiempo
de calentamiento y calor absorbido/removido al/del sólido.
115
REFERENCIAS
Referencias Bibliográficas
- AGUILERA F., “Simulación del proceso global de escaldado – hidroenfriado de

productos vegetales”, Tesis de Ingeniería en Alimentos, Universidad Austral de
Chile, 1996.
- APRAIZ B. J. "Tratamientos térmicos de los aceros". Editorial: Dossat, 1974.
- BURNETTE, F.S., “Peroxidase and its relationship to food flavor and quality: a
review, J.Food Sc, 1977.
- DÍAZ C., “Análisis teórico experimental del coeficiente de transferencia de calor
durante los procesos de hidrocalentamiento – hidroenfriamiento de alimentos”,
Tesis de Ingeniería en Alimentos, Universidad Austral de Chile, 2001.
- FEIHL, O. “Heating frozen and nofrozen veneer logs”, Forest Product Journal,
22(10): 41-50, 1972.
- FEIHL, O.; GODIN, V., “Heating veneer logs : a practical guide”, Can Dept. Env.,
Eastern Forest Prod. Lab., Otawa, Ontario, Canadá, Forestry Tech. Report 9.19
pp., 1975
- CHOI, Y. And OKOS, M.R., “Effects of temperature and composition on the
thermal properties of food, In: Food engineering and process application. Vol I,
Transport phenomena, Le Maguer, M. and Jelen, P. (eds), Elsiever Apllied
Sciencie Co., New York, 1986.
- FENNEMA, O., “Química de los alimentos, Acribia s.a., Zaragoza, España, 1993.
- GANTHAVORN, C and POWERS, J.R., “Changes in peroxidase activity,
hexanal, ascorbic acid and free sulfhydril in blanched asparagus during frozen
storage”, J Food Sc 53: 1403-1405, 1988.
- HELDMAN, D.R., “Food process engineering”, Ed. Avi publishing company,
USA, 1975.
- INCROPERA, F.P. y DEWITT D.P., “Fundamentos de transferencia de calor”, 4ª
Ed. Prentice Hall, México, 1999.
- KREITH F. & BOHN M., “Principios de Transferencia de Calor”, Ed. Thomson
Learning, Mexico, 2001.
116
- LUH, B.S. and LORENZO, M.C., “Freezing of vegetables in: commercial

vegetable processing”, Luh, B.S. and Woodroof, J.G (eds), Van Nostrand
Reinhold, New York, 1988.
- LUND, D.B., “Design of thermal process for maximizing nutrient retention.”,
Food Technol, 31(2):71-78,1977.
- McAdams, W. H., “Heat Transmission”, 3d ed., McGraw-Hill Book Company,
New York, 1954.
- MARTINEZ y GORDON, “Teoría y Tecnología del Tratamiento Térmico”, Ed.
Pueblo y Educación, Ciudad de la Habana, Cuba, 1985.
- MORAGA N. y MURILLO H., “Transferencia de calor en congelación de
alimentos cárneos con métodos de elementos y de diferencias finitas”, III
Congreso iberoamericano de ingeniería mecánica, CIDIM, 1997.
- MORAGA N. y SÁNCHEZ L., “Predicción de temperaturas y velocidades de
enfriamiento para temple de placas y cubos de acero en agua”, XII Jornadas de
transferencia de calor y materia, Universidad de Santiago de Chile, Chile, 1994.
- MORAGA N., FORTUNY P. y PAVICEVIC Y., “Modelación térmica y
validación experimental para descongelación de esfera de carne”, II Encuentro de
modelos físicos y matemáticos en ingeniería, Universidad de Santiago de Chile,
1997.
- MORAGA N., MEDINA E. y MENESES A., “Análisis de transferencia
simultanea de calor y materia en congelación de esferas de carne”, XIV Jornadas
de transferencia de calor y materia, Universidad Católica de Valparaíso, Chile,
1996.
- MORENO R., DEVLIEGER F. and INZUNZA L. – “Finite Difference Simulation
for Heat Conduction Applied to Heating Veneer Logs”, Methods numerical in
engineering and applied sciences, Edit. Cimne, Barcelona, España, 1992.
- Q. TUAN PHAM, “Calculation of processing time and heat load during food
refrigeration”, “Food for Thought – Cool” AIRAH Conference, Sydney, Australia,
2002.
- ÖZISIK M. N., “Basic Heat Transfer”, Ed. Mc GrawHill Kogakusha, Ltd., Tokio,
1977.
117
- ÖZISIK M. N., “Heat Conduction”, Ed. John Wiley & Sons, Inc., New York,
1980.
- ÖZISIK M. N., “Finite difference methods in heat transfer”, Ed. CRC Press., USA,
1994.
- ROJAS J., “Modelamiento, simulación y optimización del proceso de escaldado y
enfriado de alimentos vegetales”, Tesis de Ingeniería Civil en Informática,
Universidad Austral de Chile, 2000.
- SHAMS, M.A. and THOMPSON, D.R., “Quantitative determination of pea losses
as affected by conventional water blanching, J. Food Sc. 52:1006-1009,1987.
- SIMPSON, W.T., “Specific gravity, moisture content, and density relationships for
wood”, Gen Tech Rep FPL-GTR-76 Madison, WI: U.S. Department of
Agriculture, Forest Service, Forest Products Laboratory, 13p.
- SIMPSON, W. T., “Heating times round and rectangular cross sections of wood in
steam”, USDA Forest Service, Forest Prod. Lab. Madison, Wis., Gen. Tech. Rep.
FPL – GTR – 130, 130 pp. 2001.
- STEINHAGEN, H. P., “Thermal conductive properties of wood, green or dry from
– 40°C to + 100°C: a literature review. USDA Forest Serv., Forest Prod. Lab.,
Madison, Wis General Tech. Report FPL-9 10PP.
- STEINHAGEN, H. P., Meyers G., Kubler H., “Heating time charts for frozen and
nonfrozen veneer logs”, Forest Products Journal, FPJ Technical Note, 27 – 37 pp.,
USA, 1980.
- VALENTAS K.J., ROSTEIN E. and SING. R.P., “Handbook of Food
Engieneering Practice”, CRC Press LLc, USA, 1997.
Referencias de Internet
- www.matweb.com , Pagina web que dispone más de 27500 datos de

propiedades de diversos materiales.
- www.fs.fed.us/publications/ , Pagina web de servicio forestal, del
departamento de agricultura de los Estados Unidos, dispone de una gran
variedad de publicaciones en el área de la madera, se bajaron las
publicaciones de: STEINHAGEN, H. P, SIMPSON, W. T, entre otros.
118
ANEXOS
Anexo 1: Propiedades térmicas variables de algunos metales y del aire.

Anexo 2: Correlaciones para transferencia de calor por convección
natural y forzada.
Anexo 3: Solución ejercicio planteado en Sección 5.4.2
____________________________________________________
Se incluye 1 CD-ROM con los siguientes archivos:
!"Instalador
!"Ejecutable Software “Procesos CCT”
!"Código Fuente
!"Trabajo Titulación en formato PDF (Acrobat Reader)
119
Anexo 1
“Propiedades térmicas variables de algunos metales y del aire”
• Propiedades térmicas variables de algunos metales.
Propiedades para el SAE 1045, según MORAGA y SÁNCHEZ (1996).
Conductividad Térmica SAE 1045
k(T ) = -1E-26T 6 + 2E-19T 5 + 2E-12T 4 + 5E-08T 3 - 8E-05T 2 + 0,0003T + 51,572

60 R2 = 1
50
40
k W/m°C
30
20 Conductividad térmica k(T)
10 Polinómica (Conductividad térmica
0
0 200 400 600 800 1000 1200
Temperatura (°C)
Calor Especifico SAE 1045

y = -5E-23x6 - 5E-11x5 + 6E-08x4 - 2E-05x3 - 0,0007x2 +
1,1676x + 395,52
R2 = 1
1400
1200
1000
Cp (J/Kg°C)
800
600
Calor Especifico Cp(T)
400
Polinómica (Calor Especifico Cp(T))
200
0
0 200 400 600 800
Temperatura (°C)
120
Propiedades para el SAE 1008, según www.matweb.com
Conductividad Térmica SAE 1008
k(T) = 3,464E-08*T3 - 5,655E-05*T2 - 2,124E-02*T + 74,27

70 R2 = 0,999
60
Conductividad Térmica k(T)
50
Polinómica (Conductividad Térmica k(T))
k (W/m°K)
40
30
20
10
0
0 250 500 750 1000 1250 1500
Temperatura (°K)
Calor Especifico SAE 1008

Cp(T ) = -2,432E-13*T 6 + 9,743E-10*T 5 - 1,578E-06*T 4 + 1,321E-03*T 3 - 6,021E-
01*T 2 + 1,418E+02*T - 1,301E+04
1200 R2 = 0,95
1000
800
Cp (J/kg°K)
600
400 Calor Especifico
200 Polinómica (Calor Especifico)
0
0 200 400 600 800 1000 1200
Temperatura (°K)
121
Propiedades para el Acero Puro y Aluminio, según KREITH y BOHN (2001).
Conductividad Térmica Acero Puro

160 k(T ) = -4,202E-13*T 5 + 1,840E-09*T 4 - 2,998E-06*T 3 +
2,285E-03*T 2 - 8,755E-01*T + 2,008E+02
140
R2 = 0,998
120
K (W/m°K)
Conductividad T érmica
100
80 Polinómica (Conductividad T érmica)
60
40
20
0
0 250 500 750 1000 1250 1500 1750
Te mperatura (°K)
Calor Especifico Acero Puro

Cp(T ) = 1,181E-11*T 5 - 4,511E-08*T 4 + 6,135E-05*T 3
- 3,627E-02*T 2 + 9,550*T - 453,2
1200
R2 = 0,93
1000 Calor Especifico

Cp (J/Kg°K)
600
400
200
0
0 250 500 750 1000 1250 1500 1750
Tempe ratura (°K)
C onductividad Térm ica A lum inio

k(T)= -8,154E-11*T5 + 1,840E-07*T4 -1,557E-04*T3 + 6,118E-02*T2 -
500 1,104E+01T + 9,613E+02
450 R2 = 9,323E-01
400
350
k(W/m°K)
300
250
200
150
Conductividad Térm ica
100
50 Polinóm ica (Conductividad Térm ica)
0
0 100 200 300 400 500 600 700 800 900
Tem peratura (°K )
122
Calor Especifico Aluminio

Cp(T ) = 3,786E-11*T 5 - 9,715E-08*T 4 + 9,788E-05*T 3 -
1400 4,840E-02*T 2 + 1,217E+01*T - 3,393E+02
R2 = 1
1200
1000
Calor Especifico
Cp (J/Kg°K)
800
400
200
0
0 200 400 600 800 1000
Tempe ratura (°K)
• Propiedades Térmicas del Aire en función de la Temperatura
Según Ecket y Drake, (1959)
Densidad (Kg/m³) 353 × T −1

Viscosidad Dinamica (Kg/m s) 5.623 × 10 −7 × T 0.623
Calor Especifico (J/Kg°K) 387 × T 0.156
Conductividad Térmica (W/m°K) 4.772 × 10 −4 × T 0.717
123
Anexo 2
“Correlaciones de Transferencia de Calor por Convección Natural”
• Correlaciones - Convección Natural
Geometría Ecuación de correlación Restricciones

104 < GrL Pr < 109
c = 0.59
N u L = c (GrL Pr) n
Placa Vertical o n = 1 / 4 (flujo laminar)
109 < GrL Pr < 1013
Cilindro vertical
McAdams (1974) c = 0.1
n = 1 / 3 (flujo turbulento)
N u L = 0.56(GrL Pr cos θ )1 / 4 105 < GrL Pr cos θ < 1011

Placa Vertical o inclinada con
0 ≤ θ ≤ 89°
la superficie caliente hacia Fujji e Imura (*) (flujo laminar )
abajo.
Placas Horizontales 105 < GrL Pr < 2 × 107

c = 0.54
n = 1 / 4 (flujo laminar)
Superficie caliente 105 < GrL Pr < 2 × 107

arriba o fría abajo McAdams (1974) c = 0.14
Superficie fría arriba o 3 × 105 < GrL Pr < 3 × 1010

c = 0.27
caliente abajo
McAdams (1974) n = 1 / 4 (flujo laminar)
124
104 < GrL Pr < 109

c = 0.53
n = 1 / 4 (flujo laminar)
Cilindro horizontal largo N u D = c (GrD Pr) n
único
109 < GrL Pr < 1012
McAdams (1974)
c = 0.13
0.589(GrD Pr)1 / 4
N uD = 2 +
Esferas [1 + (0.469 / Pr) ]9 / 16 4 / 9 GrD Pr < 1011
Pr ≥ 0.7
Churchill (*)
Nota:
Las correlaciones del coeficiente convectivo natura para los Paralelogramo, equivalen a
las de un cilindro.
Las propiedades térmicas deben ser evaluadas a temperatura pelicular media

Tm = (Ti + T∞ ) / 2 .
(*) Obtenida de INCROPERA y DEWITT (1999).
• Correlaciones - Convección Forzada.
Geometría Ecuación de correlación Restricciones

Re L < 5 × 105
N u L = 0 . 664 Re L Pr 0 .33
0 .5
Pr > 0.5
(flujo laminar)
Re L > 5 × 105
Superficie Horizontal
[
N u L = 0 . 036 Pr 0 .33 Re 0L.8 − 23200 ] Pr > 0.5
(flujo turbulento)
125
Placa alargada plana

N u L = 0 . 2 Re L
2/3
Con ancho D, 1 < Re L < 4 × 105
perpendicular (flujo laminar)
al flujo en un gas
4/5
0.62 Re 0D.5 Pr 1 / 3   Re D  5 / 8 
Cilindro largo N u L = 0 .3 + 1 +   
[1 + (0.4 / Pr) ] 2 / 3 1/ 4
  282000  
perpendicular al flujo Re D Pr > 0.2
Churchill y Berstein (*)
500 < Re < 105

Cilindro de sección no
c = 0.675
circular en flujo cruzado
N u L = c Re n n = 0.092
de un gas (sección
2500 < Re < 8000
cuadrada paralelogramo)
Jakob (**) c = 0.699
n = 0.16
Esfera en un Gas
17 < Re D < 70000
flujo perpendicular N u D = c Re n
c = 037
MCAdams (1958)
al sólido n = 0.6
Esfera en un Liquido
flujo perpendicular [
NuD = 0.97+ (0.68Re0.5 ) Pr0.3 ] 1 < Re D < 2000
Kramers (1958)
al sólido
(*) Obtenida de INCROPERA y DEWITT (1999).

(**) Obtenida de ÖZISIK M. N. (1977)
Nota: El software incorpora otras correlaciones desarrolladas en la Uach, para los

procesos de hidrocalentamiento y/o hidroenfriamiento de alimentos vegetales, bajo los
mecanismos de convección natural y forzada, tales ecuaciones pueden ser revisadas en el
trabajo de titulación de Díaz (2001).
126
Anexo 3
“Solución por Método Grafico al Problema Planteado en la Sección 5.4.2”
Solución:
Como la longitud del cilindro es 10 veces su diámetro, se puede omitir los

efectos en los extremos. Para determinar si la resistencia interna es insignificante,
primero se calcula el numero de Biot.
h × r0 (200 W/m²K)(0.1 m)
Bi = = = 0.5
k 40 W/mK
Puesto que el numero de Biot es mayor que 0.1, la resistencia interna es

significativa y no se puede usar el método de capacitancia. Para utilizar la solución
gráfica se calculan los parámetros adimensionales apropiados, de acuerdo con siguiente
expresión (numero de Fourier):
α × t (1 × 10 −5 m²/s)(20 m)(60 s/min)

Fo = 2 = = 1.2
r0 0.12 m 2
y
Bi 2 × Fo = (0.5 2 )(1.2) = 0.3
La cantidad de energía interna almacenada en el cilindro por unidad de longitud es:

k 2
Qi′ = Cpρπr 2 (Ti − T∞ ) =  πr0 (Ti − T∞ )
α 
40 W / mK
= (π )(0.12 m 2 )(350 K) = 4.4 × 10 7 Ws / m
1 × 10 m / s
-5 2
De acuerdo con la siguiente Figura 1, la temperatura adimensional en el centro

que corresponde a 1/Bi = 2.0 y Fo = 1.2 es:
T (0, t ) − T∞
= 0.35
Ti − T∞
127
En vista de que Ti − T∞ se especifica como 350°C
T∞ = 50°C , T (0, t ) = (0.35)(350) + 50 = 172.5°C

La temperatura superficial con r / r0 = 1.0 y t = 1200 s se obtiene de la Figura 2
en función de la temperatura en el centro:
T (r0 , t ) − T∞
= 0.8
T (0, t ) − T∞
.
Por lo tanto, la razón de la temperatura superficial es:
T (r0 , t ) − T∞ T (0, t ) − T∞
= 0.8 = (0.8)(0.35) = 0.28
Ti − T∞ Ti − T∞
y la temperatura superficial después de 20 minutos es:
T (r0 , t ) = (0.28)(350) + 50 = 148°C

Entonces, es posible determinar la cantidad de calor transferida de la varilla de
acero a partir de la Figura 2 (b). Como Q ′´(t ) / Qi′ = 0.61 , resulta:
(2 m)(4.4 × 10 7 Ws / m)
Q(t ) = (0.61) = 14.9 kWhr
3600 s/hr
128
FIGURA N° 1 Temperaturas transitorias adimensionales como función del tiempo en un cilindro largo.
129
FIGURA N° 2 Flujo de calor en un cilindro largo como función del tiempo en un cilindro largo.

Edison Collinao PDF

Cargado por

Copyright:

Formatos disponibles

Edison Collinao PDF

Cargado por

Información del documento

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

Edison Collinao PDF

Cargado por

Copyright:

Formatos disponibles

Software para Estudios de Conducción de Calor

Transiente en Procesos Térmicos Industriales

Trabajo para optar al Título de:

EDISON ADRIAN COLLINAO ROBLES

EDISON ADRIAN COLLINAO ROBLES

Ha sido aprobado en el examen de defensa rendido el día , como

Sr. Rolando Ríos Rodríguez, M.Sc en Ing. Mecánica ______________

V°B° Director de Escuela

Sr. Enrique Salinas A. ______________

Extiendo mi aprecio sincero a Don Rogelio Moreno, por su apoyo y

Además, agradezco a Don Jorge Morales, por su aporte en la conformación

Por ultimo, agradezco a la Corporación de Desarrollo Indígena CONADI y al

En está, una de las etapas de mi vida,

A mis queridas abuelas Matilde y

A mis Padres, por el apoyo en todo lo

Capítulo 2 – REVISIÓN BIBLIOGRÁFICA 7

Capítulo 3 – DESARROLLO DE LOS MODELOS MATEMÁTICOS 29

Capítulo 4 – IMPLEMENTACIÓN DEL SOFTWARE 76

Capítulo 5 – COMPARACIÓN Y DISCUSIÓN DE LOS RESULTADOS 89

Capítulo 6 – CONCLUSIONES 113

a : Ancho de la placa extensa (m)

Ti n : Temperatura en el nodo i para el tiempo n = 0 (°K)

Ti n +1 : Temperatura en el nodo i para el tiempo n + 1 (°K)

αr(T ) : Difusividad térmica radial (m²/s)

Existen numerosas aplicaciones de ingeniería, donde se encuentra el uso de

El presente trabajo, consiste en elaborar un software con fines académicos, para

Los modelos fueron resueltos a través del método de diferencias finitas en su

Los resultados obtenidos en las áreas de estudio (maderas, alimentos, metales),

En términos generales, el software ofrece resolver una variedad de problemas, sin

The present work consists of elaborating a software with academic purposes, to

In general terms, the software offers to solve a variety of problems, without

Una de las necesidades en la industria chilena maderera, metalmecánica,

Como ejemplo se tienen algunos de estos procesos contenidos en los campos

Específicamente, en el caso de la industria de contrachapado, existe el proceso de

El uso adecuado de todos estos procesos y otros, obedece al hecho de plantear

Existen publicaciones de tablas para determinar tiempos de calentamiento de

Por otro lado en alimentos se presentan las mismas limitantes en la información,

Así el objetivo de este trabajo es entregar una solución a los problemas de

Además se define en la interfaz sólido – fluido la transferencia de calor por

Se desarrolla una herramienta computacional capaz de simular el

Finalmente se realizan comparaciones entre resultados obtenidos por el software

!" Desarrollar un modelo matemático para describir el comportamiento cinético que

!" Confeccionar un programa computacional capaz de simular el proceso térmico

1.3 Metodología de Trabajo

Como uno de los objetivos primarios será aplicar el método numérico a la

Se establecerá un diagrama de flujo o algoritmo para la solución numérica

Finalmente se realizará la comparación de los resultados encontrados tanto en

Capítulo 2 - REVISON BIBLIOGRAFICA

El Capítulo presenta una revisión bibliográfica, que comprende una descripción

Se revisa además los estudios realizados y el nivel actual, en relación a la

Además se revisan los métodos numéricos mas utilizados en la solución de los

2.2 Descripción de algunos Procesos Térmicos de Calentamiento y/o

2.2.1 Proceso de Macerado

El proceso de macerado de trozas de madera es un proceso térmico utilizado

El ablandamiento se lleva a cabo en la mayoría de las industrias nacionales en

temperatura y de tiempo de calentamiento, por lo que es conveniente desarrollar un

FIGURA N° 1 Proceso de macerado en piletas de agua

El tiempo de calentamiento tiene que ser suficientemente largo para calentar la

CUADRO N° 1 Temperatura y tiempo de cocción según FEIHL (1972).