Table Of ContentUNIVERSIDAD NACIONAL AGRARIA
LA MOLINA
ESCUELA DE POSGRADO
MAESTRÍA EN TECNOLOGÍA DE ALIMENTOS
“MODELAMIENTO DE LA TRANSFERENCIA DE CALOR
DURANTE EL ESCALDADO DEL LOCHE (Cucurbita moschata)
Y PAPA (Solanum tuberosum)”
Presentada por:
JULIO MAURICIO VIDAURRE RUIZ
TESIS PARA OPTAR EL GRADO DE MAGISTER SCIENTIAE EN
TECNOLOGÍA DE ALIMENTOS
Lima - Perú
2015
Dedico este trabajo a todas las personas interesadas en la ingeniería de alimentos, a mi
familia y a mis mentores de la UNALM.
"Si he logrado ver más lejos, ha sido porque he subido a hombros de gigantes"
Issac Newton (1643-1727)
AGRADECIMIENTOS
- Quiero agradecer a todas las personas involucradas en este trabajo, por sus opiniones,
sugerencias, apoyo y sobre todo por la confianza brindada y por creer en mí.
- A mis padres y hermana por darme la oportunidad de continuar estudiando e
investigando.
- A mi asesor el M.Sc. Walter Francisco Salas Valerio, por encaminarme a la
investigación en ingeniería de alimentos.
- A los miembros del comité consejero por la paciencia en leer mi trabajo y por las
correcciones realizadas.
- Al M.Sc Miguel Ángel Solano Cornejo e Ing. William Sánchez Chávez, por
convertirse en críticos de este trabajo.
- A la familia Laca Reyes, por acogerme en su hogar durante mi estancia en la
Universidad Nacional Agraria La Molina y por contribuir en mi educación.
- A la Ing. Lady Judith Cabrera De la Cruz, por motivarme a la investigación en
ingeniería y ser mi motivo para continuar mejorando.
- A todos los profesores que me enseñaron durante mi estancia en la Universidad
Nacional Agraria La Molina.
ÍNDICE GENERAL
Pág.
RESUMEN
ABSTRACT
I. INTRODUCCIÓN 1
II. REVISIÓN DE LITERATURA 3
2.1 EL LOCHE (Cucurbita moschata Duch.) 3
2.1.1 Generalidades 3
2.1.2 Composición química y usos 4
2.2 LA PAPA (Solanum tuberosum L.) 5
2.2.1 Generalidades 5
2.2.2 Composición química y usos 6
2.3 EL ESCALDADO 7
2.3.1 Generalidades 7
2.4 TRANSFERENCIA DE CALOR DURANTE EL ESCALDADO 10
2.4.1 Solución analítica de la ecuación de difusión de calor 10
2.4.2 Solución numérica de la ecuación de difusión de calor 17
2.4.3 Método de diferencias finitas 18
2.5CONDICIONES INICIALES Y DE FRONTERA EN LA ECUACIÓN DE 23
DIFUSIÓN DE CALOR
2.5.1 Generalidades 21
2.5.2 Métodos para determinar (h) 25
2.6 DIFUSIVIDAD TÉRMICA (α) 27
2.6.1 Generalidades 27
2.6.2 Métodos para determinar (α) 28
III. MATERIALES Y MÉTODOS 31
3.1 LUGAR DE EJECUCIÓN 31
3.2 MATERIA PRIMA 31
3.3 EQUIPOS, MATERIALES Y REACTIVOS 31
3.3.1 Equipos 31
3.3.2 Materiales 32
3.3.3 Reactivos 33
3.4 MÉTODOS DE CONTROL 32
3.4.1 Análisis físicos y químicos de las materias primas 33
3.4.2 Análisis proximal de las materias primas 34
3.5 ESCALDADO DEL LOCHE Y PAPA 34
3.6 MEDICIÓN DE LA TEMPERATURA 35
3.7 DETERMINACIÓN DEL COEFICIENTE DE TRANSFERENCIA DE 36
CALOR (h)
3.8 DETERMINACIÓN EXPERIMENTAL DE LA DIFUSIVIDAD 37
TÉRMICA (α)
3.9 DETERMINACIÓN DE LA CONDUCTIVIDAD TÉRMICA (k) 38
3.10 SIMULACIÓN DE LA TRANSFERENCIA DE CALOR MEDIANTE LA 39
SOLUCIÓN ANALÍTICA EN 3D
3.11 SIMULACIÓN DE LA TRANSFERENCIA DE CALOR MEDIANTE 40
DIFERENCIAS FINITAS EN 3D
3.12 SIMULACIÓN DE LA TRANSFERENCIA DE CALOR MEDIANTE 48
DIFERENCIAS FINITAS EN 3D CON PROPIEDADES TÉRMICAS
VARIABLES
3.13 ANÁLISIS ESTADÍSTICO DE LOS DATOS 49
IV. RESULTADOS Y DISCUSIONES 50
4.1 CARACTERIZACIÓN DE LAS MATERIAS PRIMAS 50
4.1.1 Características físico químicas y composición proximal del loche y la papa 50
4.2 DETERMINACIÓN DEL COEFICIENTE DE TRANSFERENCIA DE 51
CALOR POR CONVECCIÓN (h)
4.3 DETERMINACIÓN DE LA DIFUSIVIDAD TÉRMICA (α) DEL LOCHE 53
Y PAPA A DIFERENTES TEMPERATURAS
4.4 DETERMINACIÓN LA CONDUCTIVIDAD TÉRMICA (k) DEL LOCHE 56
Y PAPA A DIFERENTES TEMPERATURAS
4.5 SIMULACIÓN DEL ESCALDADO DEL LOCHE Y LA PAPA 57
4.5.1 Método de separación de variable (solución analítica) 57
4.5.2 Método de diferencias finitas explícitas en 3D 58
4.5.3 Validación del método diferencias finitas explícitas en 3D 60
4.5.4 Inclusión de la difusividad térmica y conductividad térmica variables con 63
la temperatura en la simulación por diferencias finitas
4.6 VALIDACIÓN DE LA SIMULACIÓN POR DIFERENCIAS FINITAS EN 64
3D CON PROPIEDADES TÉRMICAS VARIABLES CON LOS DATOS
EXPERIMENTALES DEL ESCALDADO DEL LOCHE Y LA PAPA
4.6.1 Escaldado del loche 64
4.6.2 Escaldado de la papa 68
4.6.3 Verificación de método durante el escaldado del loche y la papa en la 73
coordenada (10,1,1)
V. CONCLUSIONES 76
VI. RECOMENDACIONES 77
VIII. REFERENCIAS BIBLIOGRÁFICAS 78
IX. ANEXOS 92
ÍNDICE DE TABLAS
Tabla 1: Composición proximal del zapallo loche por cada 100g de materia 4
húmeda
Tabla 2: Composición proximal de la papa por cada 100g de materia húmeda 7
Tabla 3: Resumen de la ecuación de difusión de calor en diferentes 11
coordenadas geométricas
Tabla 4: Resumen de las soluciones para la conducción transitoria 13
unidimensional en una placa infinita de espesor 2L, un cilindro
infinito de radio ro y una esfera de radio ro, sujetos a convección
desde todas las superficies*
Tabla 5: Condiciones de frontera para la ecuación de difusión de calor en la 24
superficie (x=0)
Tabla 6: Valores típicos del coeficiente de transferencia de calor (h) en 26
operaciones de procesamiento de alimentos
Tabla 7: Valores típicos de difusividad térmica de algunas hortalizas 29
Tabla 8: Propiedades térmicas y físicas de la papa según la literatura 48
Tabla 9: Características físico químicas del loche y la papa 50
Tabla 10: Composición proximal del loche y la papa por cada 100 g 51
Tabla 11: Coeficientes de transferencia de calor (h) de diferentes tamaños de 52
cubos de aluminio, sometidos a diferentes temperaturas de
calentamiento
Tabla 12: Difusividad térmica (α) del loche y papa a diferentes temperaturas de 53
calentamiento
Tabla 13: Variaciones de tiempo en segundo máximas que se deben utilizan en 60
diferentes redes nodales para que la simulación por diferencias
finitas sea estable
Tabla 14: Parámetros estadísticos de ajuste de los datos experimentales con los 65
datos simulados para cubos de loche de 1x1x1 cm
Tabla 15: Parámetros estadísticos de ajuste de los datos experimentales con los 66
datos simulados para cubos de loche de 2x2x2 cm
Tabla 16: Parámetros estadísticos de ajuste de los datos experimentales con los 67
datos simulados para cubos de loche de 3x3x3 cm
Tabla 17: Parámetros estadísticos de ajuste de los datos experimentales con los 69
datos simulados para cubos de papa de 1x1x1 cm
Tabla 18: Parámetros estadísticos de ajuste de los datos experimentales con los 70
datos simulados para cubos de papa de 2x2x2 cm
Tabla 19: Parámetros estadísticos de ajuste de los datos experimentales con los 71
datos simulados para cubos de papa de 3x3x3 cm
Tabla 20: Parámetros estadísticos de ajuste de los datos experimentales con los 75
datos simulados para cubos de loche y papa de 3x3x3 en la
coordenada (10, 1, 1)
ÍNDICE DE FIGURAS
Figura 1: Efecto del escaldado en los tejidos celulares: S, almidón 8
gelatinizado; CM, membrana citoplasmática alterada; CW, pared
celular poco alterada; P, pectinas modificadas; N, núcleo y las
proteínas citoplasmáticas desnaturalizadas; C, cloroplastos y
cromoplastos distorsionados
Figura 2: Intersección de dos placas infinitas 15
Figura 3: Intersección de un cilindro infinito y una placa infinita 15
Figura 4.1: Posición central (1, 1, 1) de inserción de la termocupla 35
Figura 4.2: Posición (10, 1, 1) de inserción de la termocupla 35
Figura 5: Sistema de adquisición de datos de temperatura durante el escaldado 35
de los vegetales
Figura 6: Sistema de adquisición de datos de transferencia de calor en un 37
cilindro infinito hueco de aluminio, para determinar difusividad
térmica.
Figura 7: Intersección de tres placas infinitas con el mismo espesor, para 39
formar un cubo finito en 3D
Figura 8: Generación de la red nodal en la octava parte de un cubo en 3D 41
Figura 9: Nodo interno conectado con 6 nodos, en coordenadas cartesianas 42
Figura 10: Nodo externo conectado con 5 nodos, en coordenadas cartesianas 44
Figura 11: Nodo externo conectado con 4 nodos, en coordenadas cartesianas 45
Figura 12: Nodo externo conectado con 3 nodos, en coordenadas cartesianas 46
Figura 13: Difusividad térmica (m2/s) del loche y la papa a diferentes 54
temperaturas
Figura 14: Interfaz de la aplicación computacional para simular la transferencia 58
de calor mediante la solución analítica
Figura 15: Interfaz de la aplicación computacional para simular la transferencia 59
de calor mediante diferencias finitas explícitas
Figura 16: Ventana de reporte de criterios de estabilidad y convergencia de la 59
simulación por diferencias finitas explícitas en 3D
Figura 17: Comparación entre la simulación analítica v/s la simulación 61
numérica con una red 5x5x5 nodos para cada eje
Figura 18: Comparación entre la simulación analítica v/s la simulación 61
numérica con una red 10x10x10 nodos para cada eje
Figura 19: Comparación entre la simulación analítica v/s la simulación 62
numérica con una red 15x15x15 nodos para cada eje
Figura 20: Comparación entre la simulación analítica v/s la simulación 62
numérica con una red 20x20x20 nodos para cada eje
Figura 21: Comparación entre las simulaciones del escaldado de cubos de loche 65
de 1x1x1 cm con los datos experimentales
Figura 22: Comparación entre las simulaciones del escaldado de cubos de loche 66
de 2x2x2 cm con los datos experimentales
Figura 23: Comparación entre las simulaciones del escaldado de cubos de loche 67
de 3x3x3 cm con los datos experimentales
Figura 24: Comparación entre las simulaciones del escaldado de cubos de papa 69
de 1x1x1 cm con los datos experimentales
Figura 25: Comparación entre las simulaciones del escaldado de cubos de papa 70
de 2x2x2 cm con los datos experimentales
Figura 26: Comparación entre las simulaciones del escaldado de cubos de papa 71
de 3x3x3 cm con los datos experimentales
Figura 27: Comparación entre las simulaciones del escaldado de cubos de loche 74
y papa de 3x3x3 cm en la coordenada (10, 1, 1) con los datos
experimentales
Description:2.4.1 Solución analítica de la ecuación de difusión de calor analítica, esta se logra para la cubos de 3x3x3 cm cuando se genera una red nodal de