Sistemas Automáticos de Control 2015

Índice

1. Introducción 3
2. Planteamiento del problema 5
3. Objetivos 6
4. Marco teórico 7
4.1 Procedimientos de secado 7
4.2 Funcionamiento del proyecto propuesto 8
4.2.1 Estudio analítico del colector 9
4.2.2 Estudio analítico de la cámara 10
5. Estudio de campo 11
6. Modelamiento matemático del proceso de secado: Identificación de
variables 15
6.1 Análisis térmico del colector 15
6.1.1 Ecuaciones a considerar en el colector 16
6.2 Análisis térmico y diseño de la cámara de secado 20
6.2.1 Análisis térmico de la cámara con el fruto 20
6.2.2 Diseño de la cámara 21
6.2.3 Simulación en Solidworks 22
7. Simulación del modelo en Matlab 24
7.1 Irradiación solar 25
7.2 Colector Solar 26
7.2.1 Número de Reynolds 26
7.2.2 Número de Nusselt 27
7.2.3 Coeficientes de convección 27
7.2.4 Valores de temperatura 28
7.3 Cámara de Secado 29
7.3.1 Humedad del cacao 29
7.3.2 Constante ‘b’ de los sistemas concentrados 30
7.3.3 Variación de la masa total del cacao 31
7.4 Resultados finales 32

1
 Sistemas Automáticos de Control 2015

8. Función de transferencia del proceso 33

8.1 Salida con el modelo 34
8.2 Salida con la función transferencia 35
9. Conclusiones 37
10.Bibliografía 38

2
 Sistemas Automáticos de Control 2015

1. INTRODUCCIÓN

Este inicio de año la MEF espera que la economía peruana crezca un 4.2% lo cual
equivale a un incremento del 1.4% del PBI, tales números son esperanzadores ya que
se creía una crisis económica en el país.

Esta mejora del PBI se debe a las diferentes tecnologías y conocimientos comerciales
introducidos por empresas extranjeras en diferentes sectores de la industria peruana.
Gracias a sus inversiones mejoran la calidad de vida de cada trabajador, esto en
resumen es una disminución de pobreza, lo cual es una motivación a seguir creciendo y
apostando por los recursos naturales producidos y su respectiva exportación.

Esto se demuestra claramente en la agroindustria peruana con la producción del

cacao. En el año 2014 se exportaron 234 millones de dólares americanos y la
producción nacional superó las 76 mil toneladas de dicho producto, lo que generó
alrededor de 7.5 millones de jornales anuales y benefició de manera directa a más de
90 mil familias, según informó Marco Vinelli, director general de Negocios del
Ministerio de Agricultura y Riego.

Una de las regiones responsables de esta producción se encuentra en Piura. Esta

región produce cacao orgánico de manera artesanal y la Universidad de Piura está al
tanto de este medio de producción, lo cual no está mal sin embargo como universidad
siempre se buscarán medios para mejorar la vía de producción y que beneficien a la
población.

Y a consecuencia de esto el presente trabajó plantea una solución a uno de los tantos
procesos en la producción del cacao orgánico: el secado del cacao, esta solución busca
la mejora en el secado con uso de energías limpias, es decir un secado más rápido y
con mejores índices que lo calificarían como cacao óptimo para la exportación, y con
ayuda de conocimientos en ingeniería se podrán intuir métodos que mejoren la
rentabilidad en este proceso. Pues el aroma, textura y fermentación del cacao
dependen mucho del proceso de secado y si se mejoran estos parámetros que definen
un buen fruto del cacao, se tendrá una agroindustria mucho más competente en el
amplio mercado del cacao en el mundo, y si a esto le agregamos la valoración del
cacao peruano que tiene en productoras de chocolates y otros productos derivados de
este fruto, productoras de origen Chino, Alemán, Holandés, Italiano y otros. Será una
gran innovación que mejore la economía mucho más a lo que se cita en el inicio del
texto.

Se plantea una mejora en su proceso de secado, teniendo en cuenta que el objetivo

del presente trabajo es encontrar el comportamiento del secado de cacao mediante
colector y una cámara diseñada. Con el propósito de reducir principalmente el tiempo

3
 Sistemas Automáticos de Control 2015

de secado del cacao en base al control de la humedad, permitiéndose también la

mejora en la rentabilidad de todo el proceso.

Para poder almacenar de manera segura, es decir, mantener el valor nutritivo de los
alimentos orgánicos, especialmente de semillas o granos como el cacao, es necesario
disminuir su humedad (propia de la cosecha) hasta niveles que no permitan la
proliferación de microorganismos (hongos y bacterias) que puedan afectar en su
comercialización y consumo posterior.

En este informe final, teniendo en cuenta las desventajas del secado al aire libre
mostradas en el estudio de campo realizado y adjuntado en este mismo, tales como:
necesidad de aproximadamente 4 días de secado a sol pleno, incertidumbre sobre las
condiciones climáticas, poco control sobre la operación de secado, una gran área
requerida, exposición a la contaminación por el polvo, microorganismos u otros
materiales; lo que se reduce en una degradación de las propiedades organolépticas del
cacao.

Así mismo, también se detalla la obtención de las distintas curvas que se necesitan
para entender la dinámica del secado, empezando con el proceso del colector y su
respectivo modelo en el simulink; la entrada de calor a la cámara y su modelo en el
mismo software. Para así obtener las relaciones de temperatura en el tiempo.
Relaciones que son importantes para obtener el nivel de humedad del cacao que es el
parámetro influyente en el proceso de secado.

Se señalará también la construcción de una cámara que sirva como instrumento para
validar el modelamiento planteado por el grupo.

2. PLANTEAMIENTO DEL PROBLEMA

4
 Sistemas Automáticos de Control 2015

En este punto se mencionarán inconvenientes tanto por error humano como errores
que no se pueden evitar, o como se llaman comúnmente perturbaciones de un sistema
incontrolables, como lo son la temperatura, viento, disminución de la radiación solar y
otros.

En cuestión al error humano se derivan dos grandes hechos: el primero es que el

equipo encargado del secado de manera artesanal no lo haga de manera óptima lo
cual trae como resultado una disminución del producto tanto en cantidad como en
calidad. El segundo hecho es consecuencia del primero, si se tiene un déficit en la
producción esta no abastecerá el lote pedido para la exportación. Este problema trae
un retraso si se quiere ser un digno productor de cacao y apuntar a ser un país con el
mejor índice de calidad en cacao.

Entonces debemos plantear que se está haciendo mal y como corregirlo, es conocido
que el método para secar cacao es: dejarlo tendido en cierta área para que la radiación
reduzca la humedad, se recuerda que es la explicación más simplificada.

Esta práctica trae como consecuencia un cacao menos limpio ya que no se lleva un
control de salubridad, no se consideran las impurezas que posee el aire del ambiente,
entonces si se preselecciona el cacao “contaminado” y se aleja del mejor cacao ya
seco, es obvio que habrá una disminución con respecto a la cosecha total.

Si a este inconveniente se le agrega que no todo el cacao cosechado seca a la misma

velocidad y en el mismo intervalo de tiempo por causa de la naturaleza, la radiación no
es la misma, o simplemente un día nublado retrasa la producción. Que se puede hacer
para reducir este daño ajeno al hecho humano.

Pues con el error humano se puede mejorar notablemente con capacitaciones para el
cuidado del cacao, pero que propuesta se puede plantear para la naturaleza. Entonces
como salida inmediata se debería atenuar el repentino cambio de temperatura
aprovechando sus máximos valores, es decir cuando se da la máxima radiación y con
qué mecanismos se puede aprovechar esta.

Teniendo en cuenta este problema se tiene que plantear un estudio teórico y

experimental para así modelar un sistema de control que mejore la producción tanto
en cantidad como en calidad, estudios basados en lo aprendido en el curso de Sistemas
automáticos de control y demás cursos basados en el estudio de la energía.

5
 Sistemas Automáticos de Control 2015

3. OBJETIVOS

Se plantearon varios objetivos en el transcurso de los informes, en este último se hará

una recapitulación.

 En principio se debe identificar los parámetros que influyen en el secado del

cacao, distinguir las variables de entrada y de salida, para así hacer uso de los
conceptos teóricos. Hacer uso de fórmulas referenciadas a la ingeniería
energética, ecuaciones como: balance de masa, balance de energía y
transferencia de calor, y demás ecuaciones termodinámicas.
 Definidas las ecuaciones que nos ayuden a relacionar los parámetros
influyentes en el proceso de secado de cacao, se plantean las ecuaciones de
salida para tener una medida de lo que se desea obtener. Esto es: encontrar los
parámetros de salida como la temperatura y humedad del cacao como objetivo
principal.
 Hacer un modelo matemático en simulink para referenciar la entrada de calor,
es decir se harán dos modelos uno para el colector y otro para la cámara, esto
en función de las ecuaciones obtenidas como objetivo previo. Es importante
plasmar las ecuaciones en dicho software para así poder analizar una salida
esperada y proceder a diseñar un sistema de control.
 La construcción de una cámara de secado que junto al colector solar sirvan
como instrumentación para corroborar lo analizado en el modelo. Estos
instrumentos cumplen con cierto diseño explicado posteriormente. Esta
cámara también será diseñada en solidworks.
 Con estos objetivos que se plantean como etapas del trabajo se resumen en un
objetivo final que es el de desarrollar un control sobre el proceso de secado
para así manipular las variables de salida en todo el proceso con ayuda de
tecnologías, y así aplicar los conocimientos del curso y hacer el análisis de un
sistema automático de control.

6
 Sistemas Automáticos de Control 2015

4. MARCO TEÓRICO
4.1. PROCEDIMIENTOS DE SECADO
El secado es un proceso por el cual se le quita humedad o agua a un material, en
nuestro caso consistirá en la substracción de estos valores al cacao, esto se logra
con transferencias de calor, cabe resaltar que es el único medio, no se admiten
medios mecánicos. Esta transferencia de calor se puede dar en tres formas:
conducción, convección y radiación. Las que se usarán en el presente informe serán
las dos últimas ya que se aprovecha el viento (flujo de aire) y la radiación solar como
medios de extracción de humedad y agua.
Cabe recalcar que se puede obtener secado por otros recursos energéticos, pero el
informe hace referencia a energías renovables o limpias, específicamente hablando
un secado producido por la energía solar.
El viento transfiere calor al sistema en forma de convección y el sol por sus rayos UV
emite radiación al sistema. Se debe hacer aclarar que el viento generará una tasa de
calor al cacao siempre y cuando se encuentre a una temperatura mayor que la del
cacao, en el caso del sol esto siempre se dará pero el sol es un recurso que no está
disponible las 24 horas.
En base a estas descripciones breves de nuestros recursos para realizar el secado se
derivan dos tipos de secado:

 Secado abierto: Es un secado de manera natural cuyas tasas de calor no

siguen una dirección específica, esto retarda el proceso y genera un secado
de menor calidad, es el proceso que se practica en la región a estudiar.
 Secado cerrado: Es un secado que consiste en establecer cierta distribución
de la tasa de calor en un entorno cerrado y aislado del medio ambiente al
producto que se desea secar, es mucho mejor que el abierto ya que no se
desperdicia la tasa de calor y esto trae como consecuencia un proceso
mucho más rápido como una mejora en la calidad del producto.

En el presente informe se usará un sistema de secado del tipo cerrado y con energía
solar por las siguientes ventajas:

 Ahorro en el uso de fuentes externas como combustibles.

 Un secado de manera higiénica sin mucho costo en los materiales.
 La calidad del secado es mejor en comparación al secado solar de tipo
abierto.
 Se reducen las pérdidas que se daban por putrefacción en el secado abierto
convencional que se suele usar en la región.
 Aumenta la productividad para los agricultores lo cual es un incentivo
económico.

7
 Sistemas Automáticos de Control 2015

4.2. FUNCIONAMIENTO DEL PROYECTO PROPUESTO

En el proceso de secado de un producto cualquiera a una atmósfera determinada (aire

a cierta humedad relativa), se debe tener presente que se pueden dar las siguientes
posibilidades:

a) El producto pierda agua (se seca)

b) El producto absorbe agua (se hidrata)
c) El producto no pierde ni gana agua (punto de equilibrio)

De acuerdo con lo anterior, se deduce que en el proceso de secado de un producto

intervienen variables como la humedad relativa del aire ( ), que mientras menor es
éste será capaz de contener mayor cantidad de vapor de agua, la temperatura (T) y
velocidad del aire (v). Pero se da la influencia solo de dichas variables si el secado se
realiza al aire libre, como el secado tradicional del cacao realizado en nuestro país. De
manera que si se desea realizar un secado controlado y con medios que direccionan al
aire y además lo calientan, El sistema se hará más complejo ya que se usarán análisis
de fluidos y termodinámicos tanto en el colector y cámara. (Figura 1)

Figura 1

8
 Sistemas Automáticos de Control 2015

4.2.1. Estudio analítico del colector.

Un diagrama esquemático del colector solar es mostrado (Figura 2). Consiste en dos cubiertas
de vidrio y una placa de absorción. El mortero es el material de almacenamiento y es colocado
debajo de la placa. El colector está orientado al sol de mediodía y está inclinado un ángulo β
respecto al plano horizontal.

Figura 2
El aire fluye entre las cubiertas de vidrio y la placa de absorción, adquiriendo la energía térmica
de esta última. Durante el día el mortero almacena energía térmica, misma que durante la
noche sirve para mantener el aire caliente.

Análisis Térmico:

La radiación solar es transmitida por las cubiertas de vidrio y absorbida por la placa. EL aire
fluye entre ambos vidrios y por encima de la placa de absorción. Para el balance de energía de
los componentes que conforman el sistema se van a asumir las siguientes hipótesis:

- Las capacidades caloríficas del aire, el vidrio, la placa de aluminio y el aislante son
despreciables.
- No hay gradiente de temperatura a lo largo del espesor de las cubiertas de vidrio.
- El material de aislamiento tiene una temperatura promedio (Ts) en un tiempo (t).
- No existe estratificación perpendicular al flujo de aire en los ductos.
- El sistema está perfectamente aislado y no hay fugas de aire.

9
 Sistemas Automáticos de Control 2015

4.2.2. Estudio analítico de la cámara.

En el sistema de cámara de secado, el aire de secado ya calentado por el colector

es introducido completamente en la parte inferior del interior de la cámara. El
potencial de secado del aire caliente se reduce gradualmente a medida que pasa a
través de las bandejas superiores.
Con el fin de superar el problema de la falta de uniformidad en la distribución del
flujo de aire en los secadores convencionales del gabinete, se propone un nuevo
diseño de cámara con ayuda de las fuentes bibliográficas. Este diseño de cámara
cuenta con una sección transversal trapezoidal (Figura 3.a) a la estructura
convencional de un secador (rectangular). Esto hace que parte del aire de secado
que se introduce entre las bandejas superiores mantenga el secado a alto potencial
distribuido uniformemente dentro de las bandejas (Figura 3.b), es decir se cree que
este cambio daría lugar a una distribución más uniforme del flujo de aire y la
eliminación humedades más uniforme desde todas las bandejas simultáneamente
en la cámara.
En conclusión: Se supone que la geometría de montaje lateral de la cámara podría
tener algún efecto crucial en la distribución del aire caliente dentro de la cámara de
secado.
Es evidente que la geometría de trapezoide de la cámara (tamaño de entrada de
aire, el ángulo de cámara lado con respecto a la horizontal, ángulo de inclinación,
α, y la anchura de la cámara) afectan al patrón de flujo de aire en la cámara
secadora. En muchos secadores convencionales la instalación de techo de tipo
campana y el deflector de aire son comunes, por lo tanto, también se consideran
los efectos de estos dos parámetros. Algunas otras dimensiones importantes del
gabinete de la secadora pueden ser: tamaño de salida de aire y su ubicación,
distancia de las bandejas del producto a secar desde el techo de la secadora.

Figura 3.a Figura 3.b

10
 Sistemas Automáticos de Control 2015

5. ESTUDIO DE CAMPO

Con ayuda de alumnos egresados pudimos obtener información relevante sobre

nuestro problema, información detallada sobre el proceso de cosecha del cacao, ya
que nuestros ex compañeros tuvieron visitas a pueblos aledaños a la ciudad de
Canchaque por motivos semejantes a los nuestros.

Nos contaron que llegando al lugar donde se cosecha el cacao, éste luego de un
proceso llamado “horneado” desprendía un olor característico en la zona.

Cabe recalcar que la definición de “horneado” que nos dieron no es más que un
proceso de fermentación, en el cual los granos de cacao recién recogidos se colocan y
guardan en sacos de color oscuro y se exponen a radiación solar en este estado
durante 3 o 4 días, para que la miel propia del grano se vaya, y así obtenga este
aroma tan característico.

Este proceso de “horneado” no solo sirve para que el cacao coja aroma, sino que
también logra que el cacao se hinche y tome un color característico, lo cual aumenta
su calidad de exportación.

Nos aportaron datos como el periodo de fermentación (3 a 4 días), pero estos días
tienen que ser con sol a pleno, en el intervalo de tiempo de 8 am a 5 pm. Es por eso
uno de los motivos del presente trabajo

Seguido a este proceso, viene el proceso de secado al sol, en el cual los granos de
cacao se esparcen en una carpa y se dejan secar, moviendo los granos cada hora, así
durante aproximadamente 4 días.

Apenas termina este proceso, se tiene que almacenar el cacao nuevamente y de

manera muy ajustada en los sacos oscuros para que siga fermentando con el calor
residual en los granos y agarre un mejor aroma.

11
 Sistemas Automáticos de Control 2015

Cosecha de cacao1

1
Imágenes tomadas el año 2012 por alumnos de la Universidad de Piura.

12
 Sistemas Automáticos de Control 2015

Fermentación

Proceso de fermentación de cacao2

2
Imágenes tomadas el año 2012 por alumnos de la Universidad de Piura.

13
 Sistemas Automáticos de Control 2015

Secado

Secado de manera convencional3

3
Imágenes tomadas el año 2012 por alumnos de la Universidad de Piura.

14
 Sistemas Automáticos de Control 2015

6. MODELAMIENTO MATEMÁTICO DEL PROCESO DE SECADO:

IDENTIFICACIÓN DE VARIABLES
6.1. Análisis térmico del colector

Nomenclatura:

Símbolos Descripción Unidades

A Área m2
b ancho de la placa m
C calor específico a presión constante J kg-1 K-1
Dh diámetro hidráulico m
d ancho del ducto m
f factor de fricción
hc coeficiente de transferencia de calor por convección W m-2 k-1
hr coeficiente de transferencia de calor por radiación W m-2 k-1
It radiación solar W m-2
K conductividad térmica W m-1 k-1
L longitud del colector m
l Espesor m
ṁa flujo másico de aire kg s-1
Nu número de Nusselt
p Perímetro m
ℜ número de Reynolds
T Temperatura K
t Tiempo S
U coeficiente de transferencia de calor W m-2 k-1
V velocidad del aire m s-1
v viscosidad cinemática m2 s-1
Letras griegas
α Absortividad
β ángulo de inclinación del colector
ε Emisividad
ρ Densidad kg m-3
σ constante Stefan-Boltzmann W m-2 K-4
τ Transmisividad
Subíndices
a Aire
f fluido (aire)
f1 flujo de aire 1
f2 flujo de aire 1
m Mortero
p placa de absorción
v1 cubierta de vidrio 1
v2 cubierta de vidrio 2
Tabla n°1

15
 Sistemas Automáticos de Control 2015

6.1.2 Ecuaciones a considerar en el colector.

Se va a trabajar con las temperaturas promedios de las cubiertas de vidrio y de la placa. Así
también para modelar el intercambio de calor entre estos componentes y el flujo de aire
vamos a h uso de las temperaturas promedio en ambas corrientes de aire:

T a+ T 1 T 1+T 2
T f 1= …(6.1 .2.1) T f 2= … (6.1.2 .2)
2 2

 Para la primera cubierta de vidrio:

d Tv1
mv 1 C v 1 =α I t A v ,1 +hr v 2 , v1 ( T v 2−T v 1) A v 1−hc v 1 ,f 1 ( T v 1−T f 1 ) A v 1−hc v1 , a ( T v 1−T a ) A v 1−hr v 1 ,cielo ( T v1−
dt v1

 Para la segunda cubierta de vidrio:

d Tv2
mv 2 C v 2 =τ α I t A v ,2 +hr p , v 2 ( T p−T v 2 ) A v 2 +hc f 1 , v 2 ( T f 1 −T v2 ) A v2 −hr v 2 , v1 ( T v 2−T v 1) A v 2−hc v 2 ,f 2 ( T v 2
dt v1 v2

… (6.1.2 .4)

 Para la placa de absorción:

dTP
mP C P =τ τ v 2 α p I t A p−hr p ,v 2 ( T p−T v 2 ) A p−hc p ,f 2 ( T p−T f 2 ) A p −U p ,m ( T p−T m ) A p
dt v1

… (6.1.2 .5)

 Para el material de almacenamiento:

dTm
U p , m ( T p −T m ) A p =mm C m … (6.1.2 .6)
dt

Temperaturas del flujo de aire

Para ambos flujos de aire planteamos las ecuaciones de conservación de la energía:

d Tf1
Mf 1 Ca =hc v1 , f 1 ( T v 1−T f 1 ) A v 1−hc f 1, v 2 ( T f 1−T v2 ) A v2 …( 6.1.2 .7)
dt

d Tf2
Mf 2 Ca =hc v 2, f 2 ( T v 2−T f 2 ) A v 2−hc p , f 2 ( T f 2−T p ) A p …(6.1 .2.8)
dt

16
 Sistemas Automáticos de Control 2015

Parámetros

Comenzaremos definiendo los valores en la tabla n°2 de aquellos parámetros que conocemos,
para luego dar paso al cálculo del resto de parámetros.

Parámetros Valores Unidades

Af 0.022 m2
Ap 0.810 m2
A v ,1 0.782 m2
A v2 0.782 m2
b 0.531 m
Ca 1007 J kg-1 K-1
Cm 1050 J kg-1 K-1
d 0.041 m
Km 1.4 W m-1 k-1
L 1.526 m
lm 0.082 m
lp 0.001 m
ṁa 0.078 kg s-1
p 1.144 m
V 3 m s-1
V ext 4 m s-1
αp 0.9
αv 0.1
β 13°
εp 0.95
εv 0.8
τv 0.85
ρm 2000 kg m-3
σ 5.67x10-8 W m-2 K-4
Tabla n°2

Propiedades del aire en función de la temperatura

En base a los valores mostrados vamos a interpolar ecuaciones que me permitan definir las
propiedades del aire a cualquier temperatura dentro del rango de funcionamiento con ayuda
de la tabla n°3.4

El calor específico se mantiene constante por lo que se le ha considerado un parámetro fijo y

está incluido en la tabla n°2.

4
Tablas de propiedades de gases, Termodinámica, Yunus Cengel.

17
 Sistemas Automáticos de Control 2015

Tabla n°3

K a =0.02369+0.00007 ∙ T f …(6.1 .2.9)

−5
v=( 1.31486+0.00972 ∙ T f ) ∙10 …(6.1 .2.10)

Pr=0.7359−0.00026 ∙ T f …(6.1 .2.11)

ρa =1.26297−0.00335∙ T f …(6.1.2 .12)

−5
μ= ( 1.73786+0.00449 ∙T f ) ∙10 … (6.1.2 .13)

Caracterización del flujo de aire

Para el estudio de un flujo interno forzado es necesario determinar el número de Reynolds y el

diámetro hidráulico:

4 Af V Dh
D h= =0.076923 m… (6.1.2 .14) ℜ= …(6.1 .2.15)
p v

Con los parámetros del sistema vamos a obtener un ℜ mayor a 10000, es decir, estamos
tratando con un flujo turbulento.

Vamos a asumir que el flujo turbulento está completamente desarrollado con superficies lisas.
Esto me daría un resultado más conservador a que si se trabajara con superficies ásperas.
−2
f =( 0.79 ln ℜ−1.64 ) … (6.1 .2.16)

Cálculo de los coeficientes de transferencia de calor por convección

Primero obtenemos el número de Nusselt, que es el coeficiente de transferencia de calor

adimensional, con la ecuación de Gnielinski:

( f /8 ) ( ℜ−1000 ) Pr
Nu= … (6.1.2 .17)
1+12.7 ( f /8 ) ( Pr 2/3 −1 )
0.5

Luego se puede determinar el coeficiente de transferencia de calor por convección:

18
 Sistemas Automáticos de Control 2015

Nu K f
h= …(6.1 .2 .18)
d

Para hallar el coeficiente de transferencia de la cubierta superior de vidrio al viendo se utiliza la

siguiente relación empírica:

hc v 1 ,a =2.8+3 V ext …(6.1 .2.19)

Cálculo de los coeficientes de transferencia de calor por radiación

El coeficiente de transferencia de calor radiación entre dos superficies es:

( )
−1
1 1
hr ij =ε eq σ ( T i2 +T j2 ) ( T i +T j ) … (6.1.2 .20) ε eq= + −1 …(6.1 .2.21)
εi ε j

Medidas del colector apreciadas en la figura 4:

Figura 4.

Estas medidas ya están definidas por ser el colector propuesto por la universidad.

19
 Sistemas Automáticos de Control 2015

6.2 Análisis térmico y diseño de la cámara de secado.

6.2.1 Análisis térmico de la cámara con el fruto.

Nomenclatura dada en la tabla n°4:

Símbolos Descripción Unidades

As área superficial del grano m2
C calor específico a presión constante J kg-1 K-1
d diámetro equivalente del cacao m
h coeficiente de transferencia de calor por convección W m-2 k-1
K conductividad térmica W m-1 k-1
ṁa flujo másico de aire kg s-1
Nu número de Nusselt
ℜ número de Reynolds
T temperatura K
t tiempo S
v viscosidad cinemática m2 s-1
Letras griegas
ρ densidad kg m-3
μ viscosidad dinámica
γ humedad en base húmeda
Subíndices
a aire
c cacao
f fluido (aire)
Tabla n°4

A pesar de que el grano de cacao tiene una forma elipsoide, el considerarla una esfera no nos
va a desviar demasiado del comportamiento real del proceso. Con esta hipótesis podemos
hacer uso de la relación de Whitaker para flujos sobre una esfera:

( )
1/ 4
μ∞
Nu=( 0.4 ℜ1 /2 +0.06 ℜ2/ 3 ) Pr 0.4 …(6.2.1 .1)
μs

Luego se puede determinar el coeficiente de transferencia de calor por convección:

Nu K f
h= …(6.2 .1 .2)
d

El grano de cacao podemos considerarlo un sistema concentrado, es decir, su temperatura

interna permanece uniforme durante un proceso de transferencia de calor. Esto me da la
siguiente ecuación:

d ( T−T ∞ ) −h A s
= dt …(6.2 .1.3)
T i −T ∞ ρV cp

20
 Sistemas Automáticos de Control 2015

Humedad del cacao

Se ha trabajo con la humedad en base húmeda, esto es, la relación de humedad de obtiene al
dividir la masa de humedad entre la masa total:

mhum
γ= … (6.2.1 .4)
mhum+ mseca

Si queremos determinar la masa del cacao en función de la masa seca que es constante a lo
largo de todo el proceso y la humedad que varía en función del tiempo usaremos la siguiente
relación:

1
mtotal = ∗mseca … (6.2.1 .5)
1−γ

6.2.2 Diseño de la cámara.

Se obtuvo un criterio detallado en el informe anterior sobre la cámara construida por el grupo,
a continuación en la figura 5 y figura 6 se muestran los esquemas de la cámara hechos por el
software Solidworks.

Figura 5

21
 Sistemas Automáticos de Control 2015

Figura 6

Esta cámara será donada por el grupo para siguientes pruebas con la finalidad de validar lo
obtenido. (Imágenes de la cámara)

22
 Sistemas Automáticos de Control 2015

6.2.3 Simulación en Solidworks.

Teniendo en cuenta nuestros conocimientos en solidworks se ha tratado de representar el

comportamiento de nuestro modelo teniendo en consideraciones parte de los flujos y la
disipación de calor a través de las paredes de la cámara

- temperatura en la entrada inferior de la cámara: 60 °C

- temperatura en la chimenea superior: 56 °C

- material: madera prensada

h=20
[ ]
W
2
m K

Figura __

Figura__

23
 Sistemas Automáticos de Control 2015

Figura__
Esta simulación en el solidworks representa el comportamiento térmico de la cámara
en diferentes puntos de esta. Se debe recalcar que aunque el comportamiento es el
esperado, no es una simulación exacta, la cual se haría con un software dedicado a los
fluidos.
Como la fuente de calor proviene del flujo de aire calentado por el colector, el
solidworks solo lo representa como una fuente externa, una transferencia de calor por
convección, en la geometría determinada.
Con temperaturas dadas anteriormente, se puede apreciar que la parte con mayor
temperatura es la entrada de la cámara lo cual se aprecia en la leyenda de
temperaturas, la sección que nos importa calentar es la central la cual referencia una
temperatura promedio de 56.652 °C (ver gráfica).
Esta temperatura es función de la entrada, sin embargo se llega a la conclusión que el
diseño de cámara no hace que la temperatura baje significativamente, debido a la
dinámica que trae consigo el recorrido del fluido dentro de la cámara.

SI AGREGAMOS IMÁGENES DE LA CÁMARA ENUMERAR LAS

IMÁGENES DEL SOLIDWORKS

24
 Sistemas Automáticos de Control 2015

7. SIMULACIÓN DEL MODELO EN MATLAB.

Se ha decidido usar el Simulink de Matlab porque nos da la ventaja de poder trabajar nuestro
sistema en forma de diagrama de bloques. El cuál se aprecia en la figura __:

Figura __

Las salidas de interés en nuestro sistema son la temperatura a la salida del colector, la
temperatura de los granos de cacao y la humedad del cacao. Conviene recordar que el objetivo
de esta investigación es optimizar el secado del grano de cacao para que logre un valor de
óptimo de humedad del 8% en un tiempo de 3 días, menos de la mitad del tiempo que se
necesita en el secado a la intemperie.

Para una mejor lectura de los resultados hemos hecho uso del Workspace, el cual me permite
visualizar las temperaturas y la humedad en el transcurso de las horas, al dividirse el tiempo
(Clock) entre 3600. Para ellos en la ventana de comandos escribimos:

>> plot(t,’función’)

Todos los resultados que se van a mostrar son para un tiempo de simulación de 259200
segundos, es decir 72 horas (3 días).

25
 Sistemas Automáticos de Control 2015

7.1 IRRADIACIÓN SOLAR

La irradiación solar tiene una naturaleza senoidal, cuyo pico va a ser igual al doble de la
amplitud de la onda. Además no se da a lo largo de todo día, durante la noche la irradiación
solar es nula, es por eso que hacemos uso de una función pulso que al multiplicarse por la
función de onda logra simular este aspecto de la irradiación.

Figura__

Hay irradiación desde la 6 am que sale el sol hasta las 6 pm que este se oculta. Otro dato
importante es que en Piura la irradiación máxima en promedio alcanza los 1000 W/m2.

Figura__

26
 Sistemas Automáticos de Control 2015

7.2 COLECTOR SOLAR

Debido a que el comportamiento del colector solar es gobernado por un sistema de seis
ecuaciones diferenciales el modelo final (figura __) cuenta con una cantidad basta de
relaciones.

Figura __

A pesar de lo enorme que resulta este subsistema hubiera sido más grande y complejo si no
hubieras usado otros subsistemas dentro del mismo subsistema. A continuación se mostrarán
los más importantes.

7.2.1 Número de Reynolds

Se representa el número de Reynolds en el siguiente diagrama de bloques (Figura __)

Figura__

27
 Sistemas Automáticos de Control 2015

7.2.2 Número de Nusselt

Figura__

7.2.3 Coeficiente de Convección

Figura__

28
 Sistemas Automáticos de Control 2015

En la práctica es usual encontrar valores de coeficientes de convección que van desde los 20
hasta los 30 W/m2K. Se muestra (Figura__) qué valores nos arroja nuestro sistema para el flujo
de aire que recorre encima de la placa de aluminio:

Figura__

En nuestro caso el coeficiente de convección está oscilando entre 26.1 y 26.4 W/m 2K, lo cual es
coherente con la realidad.

7.2.4 Valores de Temperatura

Se muestra en la Figura__ los valores de temperatura de la placa de aluminio, del mortero que
almacena energía en forma de calor y del flujo a la salida del colector. Con fines ilustrativos se
va a fijar un tiempo de simulación de 24 horas para que se visualice mejor el comportamiento
de estas tres temperaturas.

Figura__

29
 Sistemas Automáticos de Control 2015

7.3 CÁMARA DE SECADO

Se aprecia en la Figura__ el diagrama de boques de la cámara.

Figura__

La salida del colector solar va a ser entrada a la cámara de secado. Este flujo de aire caliente va
a elevar la temperatura de los granos de cacaos causando que la humedad contenida en estos
vaya decreciendo.

7.3.1 Humedad del cacao

Figura__

La humedad del cacao, así como de cualquier grano en general, va decreciendo

exponencialmente, dependiendo de la temperatura del proceso de secado. En la mayoría de
procesos la temperatura de del flujo caliente es constante y para estos casos se ha
desarrollado relaciones empíricas. A pesar de esto hemos podido adaptar las relaciones
experimentales a nuestro modelo (Figura__) para que acepten una temperatura que varía con
el tiempo y que además tenga en cuenta las condiciones iniciales del ambiente. Al final se ha
llegado a un resultado que es coherente con la realidad. Para un día ver Figura__

30
 Sistemas Automáticos de Control 2015

Figura__

Al comienzo del día la humedad disminuye muy poco. Ya cerca del mediodía es cuando el
cacao empieza a perder humedad con mayor rapidez. Es por eso que la curva tiene una mayor
pendiente.

7.3.2 Constante ‘b’ de los sistemas concentrados

Con este diagrama de bloques se obtuvo la constane ver Figura__

Figura__

31
 Sistemas Automáticos de Control 2015

7.3.3 Variación de la masa total del cacao

Siendo la masa inicial del cacao de 5 kg observamos que el producto final va a pesar
aproximadamente la mitad. Esto se explica porque al inicio más de la mitad del grano era
humedad5. Después del secado sólo un pequeño porcentaje es humedad, por lo que la mayoría
del peso es de la masa seca ver Figura__.

Figura__.

Se comprueba la reducción de masa en el tiempo, cabe recalcar que se toma una muestra
genérica de 5 Kg en la cámara repartidas en las bandejas.

5
Modelling of thin layer drying kinetics of cocoa beans during artificial and natural drying.

32
 Sistemas Automáticos de Control 2015

7.4 Resultados Finales

Para un tiempo de simulación de 3 días vamos a tener lo siguiente ver Figura__.

Figura__.

Esta gráfica muestra 8verFigura__) una temperatura parecida a la simulación en el solidworks.

Figura__

Observamos que al final del proceso se obtiene un valor aproximado del 10% de

33
 Sistemas Automáticos de Control 2015

8. FUNCIÓN DE TRANSFERENCIA DEL PROCESO.

Para determinar la naturaleza del proceso necesitamos introducir una entrada escalón en el
siguiente diagrama de bloques (ver Figura__)., pero dada la característica del proceso no va a
hacer unitaria sino de magnitud 500 para que físicamente sea coherente. De haber introducido
una entrada escalón unitario se habría obtenido resultados que no reflejan el comportamiento
del sistema

Figura__

En lugar de determinar el P(s) directamente para la humedad, es preferible controlar la

temperatura del grano porque la naturaleza de la función sería de primer orden ver Figura__.

Figura__.

34
 Sistemas Automáticos de Control 2015

De la gráfica se obtiene la siguiente información:

T ss =48.84

A=T ss −20=28.84

Entonces:

20+0.632∗A=38.227

T ( t )=38.227 → t=4.269 horas=15368 seg

Ya que la magnitud de la entrada es 500:

A
K= =0.05768
500

Finalmente tenemos:

0.05768
P ( s )=
15368 s +1

Para comprobar que en efecto esa es la función transferencia compararemos la salida

obtenida por el modelo que hemos hecho y la salida del proceso que acabamos de hallar:

Salida con el modelo:

Figura__

35
 Sistemas Automáticos de Control 2015

Figura__

8.1 Salida con la función transferencia.

Se simula el proceso mediante diagrama de bloques (ver Figura__) con referencia en la

irradiación.

Figura__

36
 Sistemas Automáticos de Control 2015

La respuesta a este proceso se muestra a continuación (ver Figura__.)

Figura__.

Una forma de comprobar P(s) y el modelo es ver que tanto se asemeja a la salida del modelo
matemático propuesto. Ver Figura__.

Figura__.

Se puede comprobar que la función de transferencia calculada se aproxima muy bien al

comportamiento del sistema que tenemos. Lo que sí hay que tener en cuenta al trabajar con el
P(s) es que éste me arroja la diferencia entre la temperatura de salida y la temperatura inicial,
es por eso que comienza desde cero su respuesta.

37
 Sistemas Automáticos de Control 2015

9. CONCLUSIONES

 El análisis de sistemas concentrados utilizado para determinar la razón de transferencia de

calor por convección hacia el grano de cacao en la cámara de secado, se aproxima
bastante a la realidad de nuestro proceso porque asume que la temperatura varía con el
tiempo, pero permanece uniforme a través del sistema en todo momento.

 Debido a que no se llegó hacer ninguna prueba experimental del sistema de secado, es
decir, se desconoce la temperatura de salida se tiene que asumir un factor de
compensación o de incertidumbre en el cálculo de la temperatura en la cámara. Se asume
0.1 porque éste factor siempre debe ser menor que 1.

 El sistema de secado por radiación solar presentado en este trabajo posee una mayor
eficiencia energética comparado con otros sistemas, los cuales necesitan otra fuente de
energía como combustible o eléctrica que en los lugares donde se cosecha el cacao es
difícil de conseguir.

 La intervención de una persona es mínima para el proceso de secado con el sistema

colector-cámara. Lo que representa un ahorro de tiempo y espacio que favorece al
productor.

 Para mejorar la capacidad de recibir radiación solar del colector solar se podría idear un
control giratorio o un control de inclinación, de manera que conforme transcurre el día se
obtenga la mayor radiación posible, mediante el movimiento del colector.

 Tanto el colector como la cámara pueden mejor su poder calorífico, si se mejora la calidad
del aislante térmico y/o su constitución física.

 Para reducir el tiempo de secado lo recomendable seria aumentar el flujo de aire entrante,
es decir, subir la potencia de los ventiladores en la entrada del sistema ideado. Lo que no
se recomienda es aumentar la temperatura de una manera artificial en el colector y/o en la
cámara, debido a que se podría llegar a tostar el grano y perdería sus propiedades.

 La cámara de secado fabricada en éste trabajo tiene las dimensiones para sólo 3 bandejas,
lo que significaría pocos kilogramos de cacao por cada tres días. Debido a que se ha
buscado una relación con la exportación de cacao de Morropón, donde la cosecha del
cacao es en pocas cantidades, porque se basa en calidad más que en cantidad.

38
 Sistemas Automáticos de Control 2015

10. BIBLIOGRAFÍA.

 “Transferencia de Calor y Masa” Yunus A. Cengel

 “Performance evaluation of an inclined multi-pass solar air heater with in-built thermal
storage on deep-bed drying application” Dilip Jain y Rajeev Kuman Jain
 “Modelación y Simulación en Ambiente EcosimPro de una Cámara de Secado para
Cacao” Parra Rosero P., Saavedra R. y Ipanaqué Alama W.
 “Caracterización física de la semilla de cacao criollo, forastero amazónico y trinitario de
la localidad de Cumboto, estado Aragua” Johanna Angulo, Lucía Graziani de Fariñas,
Ligia Ortiz de Bertorelli y Pablo Parra
 “La ingeniería en el desarrollo – Manejo y tratamiento de granos poscosecha” M. De
Lucia y D. Assennato - Consultores en la FAO
 “Modelling of thin layer drying kinetics of cocoa beans during artificial and natural
drying” C.L. Hii, C. L. Law y M. Cloke

39