Revision as of 10:45, 24 January 2018 (view source) Rimni (talk | contribs) ← Older edit		Revision as of 10:47, 24 January 2018 (view source) Rimni (talk | contribs) Newer edit →
Line 314:		Line 314:
	<div class="auto" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;">		<div class="auto" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;">
	'''Tabla 3.''' Valores numéricos para los parámetros fijos, para el caso de la variación de la reflectancia.</div>		'''Tabla 3.''' Valores numéricos para los parámetros fijos, para el caso de la variación de la reflectancia.</div>
		+
	{| style="width: 100%;border-collapse: collapse;font-size: 85%;"		{| style="width: 100%;border-collapse: collapse;font-size: 85%;"
Line 345:		Line 346:
	<div class="auto" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;">		<div class="auto" style="width: auto; margin-left: auto; margin-right: auto;font-size: 75%;">
	'''Tabla  4.''' Valores numéricos para los parámetros fijos, para el caso de la variación de la carga de fluido.</div>		'''Tabla  4.''' Valores numéricos para los parámetros fijos, para el caso de la variación de la carga de fluido.</div>
		+
	{| style="width: 100%;border-collapse: collapse;font-size: 85%;"		{| style="width: 100%;border-collapse: collapse;font-size: 85%;"

---

Revision as of 10:47, 24 January 2018

Thermal model of a solar concentration cooker: a parameter variation study to estimate thermal properties

Resumen. Se realiza una variación de los parámetros de diseño para estimar propiedades térmicas que caracterizan una estufa solar, tomando de base un modelo que reproduce experimentalmente el comportamiento de una estufa en base a un concentrador parabólico compuesto (CPC) de revolución. El modelo se presenta en términos de un sistema acoplado de tres ecuaciones diferenciales no lineales, las cuales se resuelven numéricamente. Se calculan las temperaturas para el fluido, los reflectores y el contenedor de la estufa solar, en función del tiempo, a partir de eso, se infieren los valores para la potencia de cocción estándar y el rendimiento térmico.

Palabras clave: Estufa solar, Modelo térmico, Potencia de cocción, Rendimiento térmico.

Abstract

A variation of the design parameters is made to estimate some thermal properties that characterize a solar cooker, based on a thermal model that experimentally reproduces the behavior of a solar cooker based on a compound parabolic concentrator (CPC) of revolution. The model is presented in terms of a coupled system of three nonlinear differential equations, which are solved numerically. The temperatures for the fluid, the reflectors and the container of the solar cooker are calculated as a function of time. From this, the values for standard cooking power and thermal performance are inferred.

Keywords: Thermal model, cooking power, solar cooker.

1 Introducción

Se presenta un estudio de propiedades térmicas de una estufa solar con colector en base a un CPC de revolución. El estudio se realiza tomando de base un modelo térmico, en términos de un sistema de ecuaciones ordinarias acoplado, no lineales, que aproxima muy bien el comportamiento experimental de la temperatura del agua en el recipiente contenedor de la estufa solar [1]. En base a ello es posible estimar algunos de los parámetros térmicos, obtenidos con estándares internacionales, que caracterizan el desempeño térmico de las estufas solares.

Trabajos similares se han desarrollado, por ejemplo en [2], se muestra el rendimiento y la distribución de temperaturas obtenidas para una estufa solar de tipo caja, considerando tanto un modelo matemático de tipo transitorio, así como el trabajo experimental correspondiente. En [3] se presenta un modelo paramétrico de operación para una estufa solar para predecir su potencia de cocción basado en parámetros controlados y variables no controladas. En [4] se obtuvo un modelo matemático para una estufa solar de tipo caja que ayuda a establecer los elementos de los parámetros que determinan e influyen de manera relevante en el proceso de calentamiento de la estufa. En [5], se muestran los resultados de emplear un reflector plano en una estufa solar de tipo caja, estableciéndose los logros en el incremento de la temperatura de calentamiento y la eficiencia de la misma. Se han realizado otros modelos similares para estufas solares tipo caja y de concentración [6, 7, 8 y 9].

El modelo térmico se desarrolla mediante en el análisis de la física térmica de la estufa solar "Jorhejpatarnskua" [1]. El modelo simula las temperaturas promedio del agua en el recipiente absorbedor, de la superficie de la superficie del concentrador y de la superficie del recipiente absorbedor. A partir de la temperatura promedio del agua en función del tiempo, se pueden estimar la potencia de cocción estándar y el rendimiento térmico de la estufa solar, en términos de algunos parámetros de diseño: el área del colector, la carga de fluido, la reflectancia del material del concentrador y el cociente absortancia-emitancia del recubrimiento del absorbedor (olla contenedora de alimentos).

2 DESCRIPCIÓN DEL MODELO TÉRMICO

2.1. Ecuaciones del balance de energía en los componentes de la estufa solar

Considérese el caso de una estufa solar en operación, con radiación incidiendo en ella y con carga de fluido (agua) en su olla contenedora. Se plantea un modelo de transferencia de calor, que consiste en considerar el balance de energía para el recipiente absorbedor, las láminas reflectoras y el fluido, en términos de un sistema acoplado de ecuaciones diferenciales ordinarias.

2.1.1. Ecuación de transferencia de calor en el recipiente absorbedor

Realizando el balance de energía en para la superficie del recipiente absorbedor de la estufa solar, se puede expresar como:

Donde

 Es el flujo de calor de la radiación incidente sobre la olla contenedora,

 .

 Es el flujo de calor por radiación entre la olla y el cielo,

 .

 Es el flujo de calor por convección entre la olla y el ambiente,

Es el flujo de calor por convección entre la olla y el interior de la estufa.

Es el flujo de calor por radiación entre la olla y las láminas reflectoras,

Además,

 Es la temperatura ambiente.

 Es la temperatura del firmamento. Y se relaciona con la temperatura ambiente como [10]:

 .

 Es la temperatura promedio del aire entre el absorbedor y la superficie reflectora

 .

 Es la temperatura promedio del aire dentro del recipiente,

 .

 Es la temperatura promedio del agua.

 Es la temperatura promedio de la superficie de los reflectores.

 Es la temperatura promedio de la superficial de la olla contenedora.

 Es la emitancia de la película selectiva del absorbedor.

 Es la constante de Stefan Boltzmann.

 Es el coeficiente de calor por convección entre la olla y el ambiente.

 Es el coeficiente de calor por convección entre el absorbedor y el aire dentro de la estufa solar.

 Es el coeficiente de calor por convección entre la olla y el fluido.

 Es el calor específico del aluminio.

 Es la masa del recipiente de la olla.

 Es la masa de agua.

 Es la masa de las láminas reflectoras.

 Es la absortancia de los reflectores.

 Es la absortancia del absorbedor.

 Es el área de la olla contenedora.

 Es el tiempo.

 Es el calor específico de las láminas reflectoras de aluminio.

 Es el calor específico del agua.

 Es el rendimiento térmico.

 Es la reflectancia de la superficie reflectora.

 Es el número promedio de reflexiones en el CPC.

 Es la potencia de cocción.

 Es la potencia de cocción estándar.

 Es el coeficiente de calor por convección entre la superficie reflectora y el ambiente.

 Es la irradiancia directa.

 Es la irradiancia reflejada.

 Es el flujo de calor por la radiacón incidente.

 Es el flujo de calor  por convección entre las láminas reflectoras y el ambiente.

2.1.2. Ecuación de transferencia de calor en la superficie reflectora

Aplicando el balance de energía en la superficie de la superficie reflectora de la estufa, se obtiene la siguiente ecuación:

 .

Donde

 Es el flujo de calor por radiación entre la superficie reflectora y el firmamento.

 .

Donde

 Es el área del colector solar.

 Es la emitancia de las láminas reflectoras.

2.1.3. Ecuación de transferencia de calor al fluido

Finalmente, se plantea la ecuación de balance de energía en el fluido contenido en la olla:

 .

Donde

 Es el flujo de calor por convección entre la olla y el interior de la estufa,

Es el flujo de calor por radiación entre el fluido y la olla contenedora,

 .

Las ecuaciones (1), (2) y (3) representan un sistema de ecuaciones diferenciales acopladas, no lineales para las variables

 ,   y   .

Se han tomado los valores para los coeficientes de transferencia de calor por convección de [11 y 12].

3 SOLUCIÓN DEL SISTEMA DE ECUACIONES Y CÁLCULO DE PARÁMETROS

Para resolver el conjunto de ecuaciones diferenciales acopladas y no lineales, (1), (2) y (3), se ha usado el método de Runge-Kutta de cuarto orden, a través de un código programado en fortran. Las condiciones iniciales que se han tomado como la temperatura ambiente. Es decir;

 ;    ; y   

De acuerdo a los resultados obtenidos al comparar la solución numérica de la temperatura promedio del agua, con los datos experimentales, se encuentra un error menor al 5% en los casos de prueba considerados [1].

3.1 Cálculo de la potencia de cocción y del rendimiento térmico

Al resolver el sistema de ecuaciones, se obtuvieron los valores numéricos de las temperaturas promedio, de la superficie del recipiente, de la superficie reflectora y del agua, en función del tiempo. Con esta última se puede estimar la potencia de cocción mediante el uso del término de lado derecho de la ecuación (3), la potencia de cocción estándar, se define cuando la diferencia de temperatura entre la temperatura del fluido y la temperatura ambiente, es decir, [12].

Donde , es la diferencia entre la temperatura del fluido y del ambiente ( ).

La potencia de cocción estándar [12]:

.

(7)

El rendimiento térmico se calcula usando la expresión (8) ver por ejemplo [13].

(8)

4 RESULTADOS DE LAS SIMULACIONES NUMÉRICAS VARIANDO PARÁMETROS DE DISEÑO

A continuación se presenta la variación de algunos de los parámetros de diseño de la estufa solar, obteniendo las simulaciones numéricas correspondientes, para la potencia de cocción estándar y el rendimiento térmico.

El conjunto de parámetros considerados, son los siguientes: cociente absortancia-emitancia de la película selectiva absorbedora, la reflectancia del reflector del concentrador, el tamaño del colector y la cantidad fluido en la olla.

A continuación se muestran las gráficas correspondientes, considerando la variación de los parámetros.

Se consideran fijos los otros parámetros, los valores numéricos tomados para la simulación, se muestran en la tabla 1.

PARÁMETRO	VALOR
COCIENTE EMITANCIA-ABSORTANCIA	0.41
REFLECTANCIA	0.93
MASA DE AGUA (kg)	5.00

En la figura 1, se muestra la potencia de cocción estándar y el rendimiento térmico, en función del área de captación del colector de la estufa solar, la cual se ha tomado en el intervalo (0.1, 1.5) ${\textstyle m^{2}}$ , correspondiente al caso de estufas solares de pequeñas dimensiones.

Figura 1. Muestra la potencia de cocción estándar y el rendimiento térmico, en términos del área del colector de la estufa solar.

En ambos casos, se considera una curva de aproximación a los datos de la simulación numérica, relacionando la potencia de cocción estándar y el rendimiento térmico, en función del área de captación de la estufa solar.

En la figura 2 se muestran, en el eje vertical la potencia de cocción estándar y el rendimiento térmico y en el eje horizontal el cociente absortancia-emitancia, que se ha tomado en el intervalo (1, 10). El límite inferior corresponde superficies que emiten lo mismo que absorben y el superior a las que absorben 10 veces más con respecto a lo que emiten. Este límite corresponde al de películas selectivas solares optimizadas. Los valores numéricos de los parámetros fijos se muestran en la tabla 2.

PARÁMETRO	VALOR
ÁREA DEL COLECTOR ( ${\textstyle m^{2}}$ )	0.60
REFLECTANCIA	0.93
MASA DE AGUA (kg)	5.00

Figura 2. Muestra la potencia de cocción estándar y el rendimiento térmico, en términos del cociente emitancia-absortancia.

De manera similar, en la figura 3, se muestran la potencia de cocción estándar y el rendimiento térmico en función de la reflectancia de la superficie reflectora, se ha tomado en el intervalo de (0.5, 0.93). El límite inferior corresponde a la reflectancia, en el espectro solar de una superficie de acero acabado espejo y el superior a una superficie optimizada para aplicaciones solares. Los valores numéricos de los otros parámetros se muestran en la tabla 3:

PARÁMETRO	VALOR
COCIENTE EMITANCIA-ABSORTANCIA	2.20
ÁREA DEL COLECTOR ( ${\textstyle m^{2}}$ )	0.60
MASA DE AGUA (kg)	5.00

Figura 3. Muestra la potencia de cocción estándar y el rendimiento térmico, en función de la reflectancia.

Finalmente en la gráfica 4, se muestran la potencia de cocción estándar y el rendimiento térmico, en función de la cantidad de agua en el recipiente contenedor. Los valores para los demás parámetros se muestran en la tabla 4.

PARÁMETRO	VALOR
COCIENTE EMITANCIA-ABSORTANCIA	2.20
REFLECTANCIA	0.93
ÁREA DEL COLECTOR ( ${\textstyle m^{2}}$ )	0.60

Figura 4. Muestra la potencia de cocción estándar y el rendimiento térmico, en función de la masa de agua en el contenedor de la estufa solar.

En las gráficas que muestra las variación de la potencia de cocción estándar y del rendimiento térmico, en función de las diferentes variables, se han considerado una línea de ajuste a los datos obtenidos en la simulación numérica, para la potencia de cocción estándar y el rendimiento térmico, obteniendo algunas correlaciones marcadas.

Finalmente, para estudiar las variaciones de la potencia de cocción estándar y el rendimiento térmico en función de los parámetros de diseño, se ha considera la razón de cambio. A continuación se muestra en la figura 5, la razón de cambio de potencia de cocción con respecto a cada uno de los parámetros que han sido adimensionalizados y mapeados el intervalo (0,1).

Figura 5. Razón de cambio de la potencia de cocción estándar con respecto a las variables de diseño de la estufa solar.

De la misma manera, en la figura 6, se muestra la razón de cambio del rendimiento térmico de la estufa solar, en función de cada uno de los parámetros considerados.

Figura 6. Razón de cambio del rendimiento térmico con respecto a las variables de diseño de la estufa solar.

En ambos casos es notorio que la razón de cambio mayor es con respecto a la emitancia de la película selectiva absorbedora, en rangos pequeños de la emitancia. En segunda instancia, la razón de cambio con respecto a la reflectancia implica cambio significativos, sobre todo cuando la reflectancia en mayor al 90%. En el caso del área del colector y de la masa en el contenedor, la razón cambio se mantiene aproximadamente constante.

5 Conclusiones

Este estudio propone un modelo térmico para simular el comportamiento térmico de la cocina de concentración solar de Jorhejpataranskua en términos de un sistema no lineal acoplado de ecuaciones diferenciales ordinarias. Las incógnitas en estas ecuaciones diferenciales son: temperaturas promedio del agua, de la superficie del recipiente absorbente y de la superficie reflectora. Se utilizó un análisis de variación de los parámetros de diseño basado en simulaciones numéricas de, específicamente, la temperatura promedio del fluido, para calcular la potencia de cocción estándar y el rendimiento térmica. El estudio consideró las variaciones en cuatro parámetros de diseño: el área del colector, la masa de fluido en la olla, el cociente absorbancia-emitancia y la reflectancia. Se obtuvieron curvas de ajuste para los datos de las simulaciones para los cálculos de potencia de cocción estándar y rendimiento térmico, con una expresión funcional de estos en términos de las variables de diseño.

Utilizando la técnica del análisis de la razón de cambio, se encontró que la mayor variación para los parámetros de diseño correspondía al caso de la emitancia de la película absorbente sobre la superficie del recipiente que contiene el alimento a cocinar con la olla de la estufa solar. También se encontró que la tasa de cambio relacionada con la reflectancia del concentrador es significativa en los casos en que es mayor que 0,9. Finalmente, las tasas de cambio obtenidas para el tamaño del colector y la cantidad de agua en la olla permanecieron prácticamente constantes.

Agradecimientos

Agradecemos el financiamiento del proyecto CONACyT número 247719 a la Universidad Intercultural Indígena de Michoacán para la realización del presente trabajo.

Referencias

1. M. González-Avilés & J. J. González Avilés, Thermal Model of a Solar Cooker Jorhejpatarnskua. Energy Procedia 57 (2014) 1623-1631.

2. A. A. El-Sebaii, R. Domański, and M. Jaworski. Experimental and theoretical investigation of a box-type solar cooker with multi-step inner reflectors. Energy 19.10 (1994) 1011-1021.

3. P. A. Funk & D.L. Larson. Parametric model of solar cooker performance, Solar Energy, V 62.1 (1998) 63-68.

4. T. C. Thulasi Das, S. Karmakar & D. P. Rao, Solar Box Cooker: part I–modeling and part II–analysis andsimulation. Solar Energy 52.3 (1994) 274.

5. A.H. Algifri & H.A. Al-Towaie. Efficient orientation impacts of box-type solar cooker on the cooker performance. Solar Energy 70.2 (200) 165-170.

6. A.K. Bidotnark & N. Türkmen. Modelling of a hot box solar cooker. Energy Conversion and Management 37.3 (1996) 303-310

7. E.R. Pejack. Mathematical model of the thermal performance of box-type solar cookers. Renewable Energy 1.5–6 (1991) 609-615

8. K. Klemens Schwarzer & M.E. Vieira da Silva. Characterisation and design methods of solar cookers. Solar Energy 82.2 (2008) 157-163

9. J.R. Morales, M. Gordon, M., H. Terrés, A. Ortega & A. Lizardi, A. Evaluación de reflectores internos en una estufa solar de tipo caja. Científica 12.4 (2008) 175-183

10. J.A. Duffie & W.A. Beckman Solar engineering of thermal processes. Third ed. New York,

Wiley, 2013.

11. S.A. Channiwala & N.I. Doshi. Heat loss coefficients for box-type solar cookers. Solar Energy. 42.6 (1989) 495-501.

12. ASAE S580. Testing and Resporting Solar Cooker Performance. American Society of Agricultural Engineers, ASAE S580 JAN03 (2003) 824-826.

13. A. Kundapur & C.V. Sudhir. Proposal for new world standard for testing solar cookers. Journal of Engineering Science and Technology 4 (2009) 272-281.