SciELO - Scientific Electronic Library Online

 
vol.19 número1Bondability of Black locust (Robinia pseudoacacia) and Beech wood (Fagus sylvatica) with polyvinyl acetate and polyurethane adhesivesEffects of adhesive and loading directions on the load-carrying capacity of V-nails índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


Maderas. Ciencia y tecnología

versión On-line ISSN 0718-221X

Maderas, Cienc. tecnol. vol.19 no.1 Concepción  2017

http://dx.doi.org/10.4067/S0718-221X2017005000009 

ARTÍCULO

 

Modeling and simulation of an industrial indirect solar dryer for Iroko wood (Chlorophora excelsa) in a tropical environment


Merlin Simo-Tagne1,♠, André Zoulalian2, Romain Remond3, Yann Rogaume3, Beguidé Bonoma4

1LERMaB, Post-Doctoral position, 27 rue Philippe Séguin, P.O. Box 1041, F-88051 Epinal, France.
2University of Lorraine, LERMaB, Faculty of Sciences and Technologies, Aiguillettes Campus, P.O. Box 70239-54506 VlN, Nancy, France.
3University of Lorraine, LERMaB, ENSTIB, 27 rue Philippe Séguin, P.O. Box 1041, F-88051 Epinal, France.
4Higher Teacher’s Training College, Applied Physic Laboratory, P.O.Box 47, Yaoundé, Cameroon.

Corresponding author


ABSTRACT

This paper presents a modeling of an instrumental indirect solar wood dryer less expensive functioning in a Cameroonian climate applied to the climate of Yaoundé. The dryer is easy to build and electric energy is only used for the fan. Applications are done on Iroko wood (Chlorophora excelsa), a tropical wood 50mm thick most utilized in Africa. A satisfactory agreement between experimental and numerical results was found. Influences of thickness, wood initial water content and airflow rate were studied.

Keywords: Drying curves, numerical simulation, solar greenhouse dryer, tropical climate, tropical woods.


 

INTRODUCTION

In Cameroon, forest exploitation gives provides considerable to the state. In addition, natural forests represent an important area. Iroko (Chlorophora excelsa) is prohibited exporting in the form of rough lumber. Then, it is very important to develop the first transformations in the country. For this reason, the reflections to optimize the drying are necessary in function of the local realities. Locally, after sawing, natural drying is used in the majority to dry wood. This method is known to have a lot of drying time and to destroy more boards than solar drying in an appropriate dryer. Solar drying is a method recommended where sustainable development is essential to reduce the environment impact of the industrial activities.

Cameroon has good solar potential with an average incident solar energy from 4,5 kWh/m2 in the south and 5,8kWh/m2 daily in the great north (Ayangma et al. 2008). For example, the study of Ayangma et al. (2008) using 20 years data and applied in Garoua, a city of the great north region of Cameroon gives an estimation near 5,743 kWh/(m2.day). Another works in the literature such as Lealea and Tchinda (2013) and Njomo and Wald (2006) present some estimations of solar irradiation of many Cameroonian towns in general. Figure 1 presents the position of Cameroon's incident solar energy with some others developing countries.

In the literature, we have solar dryers with the wall in glass or with energy storage (Awadalla et al. 2004, Bauer 2003, Bekkioui 2009, Bekkioui et al. 2009, Bekkioui et al. 2011, Bentayeb et al. 2008a, Bentayeb et al. 2008b, Luna 2008, Luna et al. 2010). These dryers are very expensive for the level of development of many people coming in tropical region. Mathematical modeling is an essential tool that facilitates the understanding of the physical phenomena occurring in the product to dry during solar drying. It is possible to numerically model a dryer in order to give the characteristics to impose during the construction and to satisfy the needs of a population.


Figure 1. Global solar irradiation on a horizontal plane of some developing countries (Weiss and Buchinger 2003).

 

In this work, a modeling and a numerical simulation of the instrumental solar dryer of tropical woods coming from Cameroon are done. The price of construction of this simple dryer is low. The environment of the drying is the one of Yaoundé. Experiments are done from 22th November to 12th December 2004. Before, we present a model for the desorption isotherms based on experimental values presented in the literature.

MATERIAL AND METHODS

Modeling of the desorption isotherms

We have used experimental values (Jannot et al. 2006) to propose a relation for the isotherms of desorption. For the equilibrium moisture content (Xeq,-) as a function of temperature (T, K) and fractional relative humidity (HR,-), we obtained:

(1a)

(1b)

(1c)

(1d)

 

Figure 2a illustrates the equilibrium moisture content of the desorption summarized in Eq (1a).

The fractional moisture content at the fiber saturation point (XPSF,-) is obtained using Eq (1a) with a value of 1 for relative humidity. This is presented in Eq (2) and plotted as a function of temperature in figure 2b.

(2)

(a)

(b)

Figure 2. (a) Experimental and theoretical (Dent's model) isotherms of desorption of Iroko wood and (b) evolution of the fiber saturation point versus temperature.

 

Solar data

The drying was done in Yaoundé. After using the data given in Kemajou et al. (2012), we obtained Equations (3a) and (3b) of the correlations of experimental values of diffuse and global insulation on a horizontal plane in Yaoundé respectively. Geographical coordinates of Yaoundé are: Latitude: 3,87°N; Longitude: 11,52°E; Altitude: 720m (Afungchui and Neba-Ngwa 2013).

(3a)

(3b)

Using the relationships given in Appendix A, we have deduced the global irradiations on a slope plane at 10° (3c) and on the wall at 90° (3d).

(3c)

(3d)

 

Where t is the time of the day in hour. Between 6am and 6pm, Gt* is given by (3c) and (3d) on the slope and the wall respectively. At night, Gt* is equal to zero. Figure 3 presents the diffuse and global solar irradiation on horizontal plane experimentally obtained by Kemajou et al. (2012) and estimations of theses physical parameters on inclined plane with 10° angle.



Figure 3. Global irradiations on a slope plane at 10° angle and on the wall (Kemajou et al. 2012).

 

Solar dryer design

The fan is the type SK012/4EHBWC. The electrical characteristics are 0,37kW power, 1350tr/min rate of rotation, 3,6A intensity, 240V tension and 50Hz frequency. The dryer dimensions are 3,10m length 2,40m width and 2,75m west height. The roof at the slope of 10 degrees with horizontal. The floor is well insulated to reduce the heat losses. Black body situated at the top of the dryer is made in a steel sheet (aluminum) with an inclination 10 degrees with horizontal, smooth 5e/10e. Another physical parameters of the steel sheet are 0,5mm thickness, 2m length, 1m width and 2,7kg mass. The air volume in the dryer is equal to the dryer volume with the reduction of the volume of all other component located in the dryer. The obtained air volume is equal to 13m3. Layout of the dryer is presented in figure 4.


Figure 4. Schematic representation of the indirect solar dryer.

 

The wall and the slope are made in polyethylene. We have 3,7kg and 5,28m2 respectively the mass and surface of the east side wall. 4,623kg and 6,6m2 are respectively the mass and surface of the west side wall. The north and south sides have the mass and surface respectively equal to 5,374kg and 7,673m2.The mass and surface of the slope are 5,292kg and 7,555m2 respectively.

Wood stack characteristics

We have 14 layers with 4 boards by layer. Two consecutive layers are separated by the sticks and each stick dimensions are 220cm length and 3,5cm thickness. Each board in the wood stack has 5cm thickness, 220cm length and 38,5cm width. Thus, the wood stack porosity is equal to 0,412. The total board surface and the total wood stack volume are respectively equal to 97,16m2 and 2,156m3. Using the wood stack characteristics and the dryer dimensions, we obtained a fill factor (FF) equal to 0,195.

Modeling of the Dryer

Hypotheses taken are:

*Each component inside the dryer is homogeneous;

*The thermophysical characteristics of the air are only influenced by the temperature. The thermophysical characteristics of the black body and the wall are constant;

*The entire wall has the same global solar radiation;

*We neglected the variation between the days of the global solar radiation;

*Natural convection is neglected;

*The floor of the dryer is supposed adiabatic;

*Transfers by conduction are neglected on all components of the dryer.

Model equations

Mass transfer on the wood stack:

We have used a purely diffusive mass transfer model such as Ananías et al. (2009, 2011), Bekkioui (2009), Bekkioui et al. (2009, 2011) and Bentayeb et al. (2008a, 2008b):

(4)

with: m00(1-ε)Vp=1108,26Kg and Sb=2LpNc(lNpc+et)=97,16m2. We have used ρ0=514,1Kg/m3 (Gérard et al. 1998).

 

Applied on temperate and Chilean woods, the literature proposes the variation of the global mass transfer coefficient K in function of the air temperature, air velocity, air relative humidity, fiber saturation point and equilibrium water content given on Equation (5) (Alvear et al. 2003, Ananías et al. 2009, Ananías et al. 2011, Bekkioui 2009, Bekkioui et al. 2009, Bekkioui et al. 2011, Bentayeb et al. 2008a, Bentayeb et al. 2008b):

(5)

a0=0,2265ms/kg; b0=268,9m2/kg; c0=2543,6K; p=2,7158

In this relation, e is used expressed in mm (Chrusciel et al. 1999). Simo-Tagne et al. (2016) show that it is possible to use the same influences of the global mass transfer coefficient in the case of tropical woods.

To estimate the external air temperature, we have used the relation (Benlahmidi 2013):

(6)

We have Tamax =303,75K and Tamin =292,95K. Initially, we have taken the temperature of all the component of the dryer equal to the value .

Mass transfer of air

The global mass balance of water in the dryer is given by Eq (7) taken from Bekkioui (2009), Bekkioui et al. (2009, 2011) and Bentayeb et al. (2008a, 2008b).

(7)

With

(8)

G=0,05kg/s.

Air relative humidity (HR,-) is given by:

(9a)

(9b)

with z=720m for Yaoundé (Afungchui and Neba-Ngwa 2013).

(9c)

 

Thermal transfer on the wood stack

We have neglected the exchange between the wood stack and the floor, also between the wood stack and the cover. Thus the rate of accumulation of thermal energy in the product is equal to the sum of the rate of thermal energy gain from the product due to sensible and latent heat, the rate of thermal energy received from air by the product due to convection and the rate of radiation between black body and wood stack. Thus we have:

(10)

 

Thermal transfer of the air

The rate of accumulation of thermal energy in the air inside the dryer is written as the sum of the rate of solar energy accumulated inside dryer from solar radiation through the roof and the wall respectively, the rate of thermal energy change in the air chamber due to inflow and outflow of the air in the chamber, the rate of thermal energy loss from air inside due to latent heat, the rates of thermal energy exchanged by convection between inside air and the roof, between inside air and the wall, between inside air and the black body and between inside air and wood stack respectively. Thus we have:

(11)

 

Thermal transfer on the black body

The rate of accumulation of thermal energy in the black body situated on the top of the drying chamber is equal to the sum of the rate of solar energy transferred by the roof and accumulated on the black body from solar radiation, the rate of thermal energy transfer by convection between the inside air and the black body, the rates of thermal energy exchanged by radiation between the roof and the black body, between the black body and the wall and between the black body and the wood stack respectively. Thus we have:

(12)

 

Thermal transfer of the roof slope in polyethylene

The rate of accumulation of thermal energy in the roof is equal to the sum of the rate of solar energy accumulated in the roof from solar radiation, the rate of thermal energy exchanged by radiation between the roof and the black body, and between the roof and the sky, the rate of thermal energy transferred by convection between the inside air and the roof, also between the outside air and the roof. Thus we have:

(13)

Thermal transfer on the wall of the dryer

The rate of accumulation of thermal energy in the wall (polyethylene cover) is equal to the sum of the rate of solar energy accumulated in the wall from solar radiation, the rate of thermal energy transfer by convection between the inside air and the wall, also between the outside air and the wall, the rate of thermal energy exchanged by radiation between the wall and the black body, and between the wall and the sky. Thus we have:

(14)

The relations below are used in the program. Form factors are obtained after calculations with the dimensions of the components of the dryer.

(15)

(16)

Vext=1,3m/s (meteorological data), Vint=1,5m/s;hciL=7w/(m2K);hci=8w/(m2K); hcto=8w/(m2K);

(17a)

We have Re=5687,36; Dh=0,07m; Pr=0,708; λair 0, 026248 w/(mK) ; hcb=8w/(m2K). Thus we have hb =7,77w/(m2K).

The thermophysical correlations on the wood and on the component of the dryer are taken in the literature Simo-Tagne (2014) and Simpson and TenWolde (1999). Those of air, black body (in aluminum) and wall (in polyethylene) are taken in the literature (Jannot 2011, Lienhard IV and Lienhard V 2011, Nadeau and Puiggali 1995). We have:

Cpto=900J/(kg.K); Cpp=CppL=2300J/(kg.K); αppL=0,05;αto=0,91;τp=0,95;

(17b)

(17c)

(17d)

(17e)

(17f)

 

Second term of second member of (17e) is equal to zero in a non-hygroscopic domain, domain that is limited by the water content at the fiber saturation points. The geometric factor is obtained using the method developed in Clark and Korybalsky (1974). This method is also used by Bekkioui (2009), Bekkioui et al. (2009, 2011) and Bentayeb et al. (2008). We have obtained:

Ftop=0,537; Fpciel=0,8; Fpto=0,81; FpLto=0,75; Ftob=0,3; Fbbb=0,323.

The parameter used to compare experimental and numerical results is given by the average relative error:

(18)

Method of resolution and experimental drying process

Method of resolution

Finite difference method (Gonçalves 2005) is used to resolve equation relative to the mass transfer on the air (Eq. 7). We obtained:

(19)

Known air humidity of the air at the time t0, (Ys,to) we deducted the air humidity at the time t1=t0+Δt(Ys,t1) . For all the others relations, we have used Runge Kutta method in the order 4 (Gonçalves 2005). For example, Eq.10 is treated such as the following:

(20a)

(20b)


Given Tbo we have:

(20c)

(20d)

(20e)

(20f)

(20g)

 

Fortran language in the version 77 was used to generate our results and the step time used was 15s. We recorded our results each 1hr drying time.

Experimental drying process

The mission of promoting local materials of Cameroon (MIPROMALO) had built and experienced the solar dryer studied in this paper. After having positioned the wood stack, the fan is commanded ON. The air relative humidity inside is controlled after each two hours during the day. If the air relative humidity is more greater than 0,7 (near 0,9), the air inside is replaced by fresh outside air. When the air relative humidity inside is lower than 0,7 (near 0,3), the fresh air replaces the air inside the dryer. At night, inside air is exchanged each 1hr in order to avoid atmospheric saturation.

RESULTS AND DISCUSSION

Known that initial drying time 0 corresponds to 00pm (local time in Yaoundé). Figure 5 shows a satisfactory agreement between theoretical and experimental measurements. We have a average relative error Er equal to 4,49%. From initial water content equal to 0,4kg/kg with the thickness equal to 50mm, 20days are necessary to dry iroko wood to 0,15kg/kg, figure 5. Figure 6 shows the influence of the environment on the water content of our wood. During the night, equilibrium water content is great and it is weak at 12 o'clock. Then, the drying of wood stack is weak during the night. It is clear that in this condition, the wood stack dried until the equilibrium value of 0,11kg/kg after 768h drying time, figure 6.

Figure 5. Experimental and theoretical water content evolution versus drying time, 50mm.




Figure 6. Theoretical water content and equilibrium water content evolutions versus drying time, 50mm.

 

If we want to use this solar dryer to dry our wood in order to use in a house in Yaoundé, knowing that the means values in the Yaoundé house are HR=0,727 and T=24,8°C (Kameni et al. 2014), relation (1a) gives equilibrium water content equal to 0,144kg/kg for our wood. In the same experimental process and climatic conditions, our indirect solar dryer gives these equilibrium values after 18 days. Figure 7 shows the time variations of temperatures relative to the wood stack (red), to the interior air (blue) and to the black body (green) during 528hrs drying time. Curves are presented 528hrs of the drying process. We notice that the curves are similar evolutions than the temperature of the air outside the dryer. Black body absorbs most energy necessary to facilitate drying process during the night. Figure 8 shows that air relative humidity vary with the time. When air temperature increases, air relative humidity decreases. Also, air temperature inside the dryer is greater than air exterior with a difference equal to 10°C.


Figure 7. Predicted temperatures of wood stack, black body and air interior during the process, 50mm.


Figure 8. Theoretical relative humidity, outside air temperature and inside air temperature, 50mm.

 

Figure 9 presents the variation of average water content of wood stack with the drying time in function to the wood thickness numerically obtained. When wood thickness decreases, drying process is rapid. When the thickness changes, final water content is near than equilibrium water content.


Figure 9. Predicted water content versus drying time and board thickness, 50mm.

 

Figure 10a shows that the drying kinetic is fast when initial water content is important because it is easy to dry free water. Also, final water content of all curves at different initial water content is near to equilibrium water content. It is clear that, if initial water content or thicknesses of wood are different, final water content of each plank is different. Thus, it is important to do a pretreatment to the wood stack in order to homogenize initial water content of the stack. Also, it is possible to do this homogenization at the end of the drying process with an increase of interior air relative humidity. Figure 10b shows that when air flow rate increases, drying kinetic decreases. Increase of air flow rate helps to decrease the value of air temperature because communication between air drying and the components of the dryer is reduced. Figure 11 presents the consequence of two types of ventilation of the dryer. Open night conduct (in the night, we remove moist air in the dryer) helps to increase drying kinetic, compared to the case where in the night fan is stopped. When moist air is not remove in the night, it is possible that interior air humidity becomes more than the board humidity at the surface. Thus, the drying kinetic becomes very lowest comparing by the case where ventilation is not stopped. It is clear that the drying duration is great in the case where ventilation is stopped during the night.



(a)

(b)

Figure 10. Predicted water content versus drying time with 50mm of thickness. (a) Influence of initial water content, (b) influence of air flow rate.


Figure 11. Predicted water content in function of the ventilation, 50mm.

 

CONCLUSIONS

The model developed gives a satisfactory agreement in comparison with experimental results. The influences of initial wood water content, airflow velocity and board thickness on the drying process are conformed. Mathematical model proposed can be used to design the performing dryer in order to build in tropical region many indirect solar dryers less expensive. During drying process, it is important to homogenize initial water content and the board thickness in order to obtain a same drying kinetic of each board. At night, it is economically important to stop the fan and open the drying chamber when inside air is near saturation. In the future, it will interesting to discuss which terms are important to carefully characterize in this model.

ACKNOWLEDGEMENTS

The principal author acknowledges the International Tropical Timber Organization (ITTO) for financially supporting a part of this work (ITTO Ref. Number: 012/15A). The authors are grateful to the administration and all staff of the mission of promoting local materials of Cameroon (MIPROMALO) for the experimental data given and all explanations obtained in order to put it in our program. We acknowledge the two anonymous reviewers for all suggestions given to ameliorate the presentation of this paper.

REFERENCES

Afungchui, D.; Neba-Ngwa, R. 2013. Global solar radiation of some regions of Cameroon using the linear Angstrom and non-linear polynomial relations (part I) model development. International Journal of Renewable Energy Research 3(4):984-992.         [ Links ]

Alvear, M.; Broche, W.; Salinas, C.; Ananias, R.A. 2003. Drying kinetic of Chilean coigüe: Study of the global drying coefficient. 8th International IUFRO Wood Drying Conference: 383-387.         [ Links ]

Ananías, R.A.; Chrusciel, L.; Zoulalian, A.; Salinas-Lira, C.; Mougel, E. 2011. Overall mass transfer coefficient for wood drying curves predictions. Mass Transfer in Multiphase Systems and its Applications. [on line] Prof. Mohamed El-Amin (Ed.), ISBN: 978-953-307- 215-9, InTech. Available from: <http://www.intechopen.com/books/mass-transfer-in-multiphase-systems-and-its-applications/overall-mass-transfer-coefficient-for-wood-drying-curves-predictions>         [ Links ]

Ananías, R.A.; Mougel, E.; Zoulalian, A. 2009. Introducing an overall mass-transfer coefficient for prediction of drying curves at low temperature drying rates. Wood Science and Technology 43(1):43-56.         [ Links ]

Awadalla, H.S.F.; El-Dib, A.F.; Mohamad, M.A.; Reuss, M.; Hussein, H.M.S. 2004. Mathematical modelling and experimental verification of wood drying process. Energy Conversion and Management (45):197-207.         [ Links ]

Ayangma, F.; Nkeng, G.E. ; Bonoma, D.B.; Nganhou, J. 2008. Evaluation du potentiel en énergie solaire au Cameroun : Cas du Nord Cameroun. African Journal of Science and Technology, Science and Engineering series 9(2):32-40.         [ Links ]

Bauer, K. 2003. Development and optimization of a low temperature drying schedule for Eucalyptus grandis (Hill) ex maiden in a solar assisted timber dryer. Ph.D. Thesis, Institut für Agrartechnik in den Tropen und Subtropen, Universität Hohenheim, Germany.         [ Links ]

Bekkioui, N. 2009. Séchage solaire du bois: modélisation simplifiée du séchage solaire d'une pile de bois dans un séchoir solaire à parois vitrées. Ph.D.Thesis, University of Mohammed V-AGDAL, Morocco.         [ Links ]

Bekkioui, N.; Hakam, A.; Zoulalian, A.; Sesbou, A.; Kortbi, M.E. 2011. Solar drying of pin lumber: Verification of a mathematical model. Maderas. Ciencia y Tecnologia 13(1):29-40.         [ Links ]

Bekkioui, N.; Zoulalian, A.; Hakam, A.; Bentayeb, F.; Sesbou, F. 2009. Modelling of a solar wood dryer with glazed walls. Maderas. Ciencia y Tecnologia 11(3):191-205.         [ Links ]

Benlahmidi, S. 2013. Etude du séchage convectif par l'énergie solaire des produits rouges. Ph.D.Thesis. University of Mohamed Khider-Biskra. Algeria.         [ Links ]

Bentayeb, F.; Bekkioui, N.; Camacho, E. F. 2008b. Simulation of a solar dryer functioning in a Moroccan climate. VI Minsk International Seminar. Heat Pipes, Heat Pumps, Refrigerators:214-219.         [ Links ]

Bentayeb, F.; Bekkioui, N.; Zeghmati, B. 2008a. Modelling and simulation of a wood solar dryer in a Moroccan climate. Renewable Energy (33):501-506.         [ Links ]

Chrusciel, L.; Mougel, E.; Zoulalian, A.; Meunier, T. 1999. Characterisation of water transfer in a low temperature convective wood drier: influence of the operating parameters on the mass transfer coefficient. Holz Roh Werkstoff (57):439-445.         [ Links ]

Clark, J.A.; Korybalski, M.E. 1974. Algebraic methods for the calculation of radiation exchange in an enclosure. Warme-und Stoffubertragung 7:31-44.         [ Links ]

Galilée, L.P. 2005. L'air humide Cours de climatisation, (Chapter 1). BTS Cours FEE 1ère Année.         [ Links ]

Gérard, J.; Kouassi, A.E.; Daigremont, C.; Détienne, P.; Fouquet, D.; Vernay, M. 1998. Synthèse sur les caractéristiques technologiques de référence des principaux bois commerciaux africains. Série FORAFRI, Document 11, CF, CIRAD, CIFOR.         [ Links ]

Gonçalves, E. 2005. Résolution numérique, discrétisation des EDP et EDO. Institut National Polytechnique de Grenoble.         [ Links ]

Jannot, Y. 2011. Thermique solaire. Cours de transfert thermique.         [ Links ]

Jannot, Y.; Kanmogne, A.; Talla, A.; Monkam, L. 2006. Experimental determination and modelling of water desorption isotherms of tropical woods: afzelia, ebony, iroko, moabi and obeche. Holz als Roh-und Werkstoff (64):121-124.         [ Links ]

Kameni, N.M.; Tchinda, R.; Orosa, J.A.; Roshan, G. 2014. Study of dioxide carbon concentration and indoor àir quality in some buildings in the equatorial region of Cameroon (Yaoundé). Iranian Journal of Health Sciences 2(2):1-15.         [ Links ]

Kemajou, A.; Mba, L.; Pako-Mbou, G. 2012. Energy efficiency in air-conditioned buildings of the tropical humid climate. IJRRAS 11(2):235-240.         [ Links ]

Lealea, T.; Tchinda, R. 2013. Estimation of diffuse solar radiation in the South of Cameroon. Journal of Energy Technologies and Policy 3(6):32-42.         [ Links ]

Lienhard IV, J.H.; Lienhard V, J.H. 2011. A heat transfer textbook. Fourth edition, Cambridge Massachusetts.         [ Links ]

Luna, D.L. 2008. Modélisation et conception préliminaire d'un séchoir solaire pour bois de pin avec stockage d'énergie. Ph.D. Thesis, ENSAM, France.         [ Links ]

Luna, D.L.; Nadeau, J.P.; Jannot, Y. 2010. Model and simulation of a solar kiln with energy storage. Renewable Energy 36(11):2533-2542.         [ Links ]

Nadeau, J.P.; Puiggali, J.R. 1995. Séchage, des processus physiques aux procédés industriels. Paris, New York, Londres, Tec and Doc.         [ Links ]

Njomo, D.; Wald, L. 2006. Solar irradiation retrieval in Cameroon from meteosat satellite imagery using Helio 2 method. ISESCO Science and Technology Vision 2(1):19-24.         [ Links ]

Simo-Tagne, M. 2014. Numerical study of heat and mass transfer during the thermal drying of tropical woods. International Journal of Thermal and Environment Engineering 8(2):9-15.         [ Links ]

Simo-Tagne, M.; Monkam, L.; Rémond, R.; Zoulalian, A.; Rogaume, Y.; Beguide-Bonoma. 2016. Experimental determination of the global mass transfer coefficient of the tropical woods in order to deduce the drying curves at the lower temperature. International Journal of Thermal and Environment Engineering 12(1):9-14.         [ Links ]

Simpson, W.T.; TenWolde, A. 1999. Physical properties and Moisture Relations of Wood. Chapter 3 from Forest Products Laboratory. Wood handbook.         [ Links ]

Weiss, W.; Buchinger, J. 2003. Solar drying. AEEE INTEC.         [ Links ]


Appendix

(Jannot 2011)

St=Gt-D; Gt, St and D are global, direct and diffuse insolation on horizontal plane. For inclined plane i, we have:

(A1)

(A2)

(A3)

(A4)

ρ is the albedo

h is height of the sun:

(A5)

(A6)

a is azimuth; Lat is latitude; i is inclination angle; δ is declination; w is horar angle giving by:

(A7)

TS is Solar time expressed in h.

(A8)

j day number in the year;

(A9)

TL is the legal time (h), ET is Equation of the time (h) lref and long are reference of longitude and longitude of the site (°) respectively. We have used C=1h and lref=0° for Yaoundé.

Nomenclature

Cpa:mass heat of the air (J/(kg.K));

Cpto:mass heat of the black body (J/(kg.K));

Cpb:mass heat of the wood (J/(kg.K));

Cpp:mass heat of the roof slope (polyethylene) (J/(kg.K));

Cppl:mass heat of the wall (polyethylene) (J/(kg.K));

D:diffuse solar irradiation on a horizontal plane (w/m2);

:derivative with respect of the drying time (s-1);

e:thickness of the plank (mm);

Eb:desorption heat (J/kg);

Er:average relative error (%);

et:stick thickness (m);

Fbb-b:geometric factor black body-wood(-);

FF:volume of the wood stack divided by the volume of the dryer chamber (-);

Fpl-ciel:geometric factor between wall and sky (-);

Fpl-to:geometric factor between wall and black body (-);

Fp-ciel:geometric factor between roof slope and sky (-);

Fto-b:geometric factor between black body and wood (-);

G:mass flow rate of dry air (kg/s);

Gt:global solar radiation on a horizontal plane (w/m2),

Gt*(10°):global solar radiation on the roof (w/m2);

Gt*(Wall):global solar radiation on the wall (w/m2);

hb:convective transfer between wood and air (w/(m2K));

hcb:interior convective coefficient between the wood stack and the air (w/(m2K));

hcext:convective coefficient between exterior air and wall (w/(m2K));

hci:interior convective coefficient between the roof slope and the air (w/(m2K));

hcil:interior convective coefficient between the wall and the air (w/(m2K));

hcto:interior convective coefficient between the black body and the air (w/(m2K));

HR:air relative humidity (-);

hvext:convective coefficient between exterior air and the roof slope (w/(m2K));

hvint:convective coefficient between interior air and the roof slope (w/(m2K));

K:mass global transfer coefficient (kg/(m2s));

l:width of each plank (m);

L:latent heat of evaporation (J/kg);

Lp:length of each plank (m);

ma:air mass (kg);

mo:wood dry mass(kg);

mp:mass of the roof slope (kg);

mto:mass of the black body (kg);

N:number of experimental wood water content in dry basis (-)

Nc:number of layers;

Npc:number of boards by layer;

Patm:pressure of the atmosphere (Pa);

Pr:Prandtl number (-);

Pvsat:pressure of vapor in saturation (Pa);

R:perfect gas constant (8,314J/(mol.K));

r:regression of determination (-);

Sb:total exchange surface of the wood (m2);

Sbb:black body surface (m2);

Sp:surface of the roof slope (m2);

Spl:surface of the wall (m2);

t:drying time (s);

Ta:interior air temperature (K);

Taext:temperature of the air exterior of the dryer (K);

Tamax:maximum air temperature (K);

Tamin:minimum air temperature (K);

Tb:wood temperature (K);

Tciel:temperature of the sky (K);

Tpl:temperature of the wall (K);

Tto:black body temperature (K);

Vext:ambient air velocity (m/s);

Vint:air velocity inside dryer(m/s);

Vdryer:volume of the dryer (m3);

Vp:wood stack volume (m3);

X:water content of the wood stack in dry basis (kg/kg);

Xeq:equilibrium water content in dry basis (kg/kg);

Xfsp:fractional moisture content in dry basis at the fiber saturation point(kg/kg);

Xm:fractional moisture content in dry basis at the monolayer saturation point (kg/kg);

Xexpi:experimental water content of the wood stack in dry basic number i (kg/kg);

Xthi:theoretical water content of the wood stack in dry basic number i (kg/kg);

YE:inlet air humidity of the dryer (kg/kg);

YS:inside air humidity of the dryer (kg/kg);

z:altitude of the dryer (the one of the town)(m);

αp:absorptivity of the roof (-);

αpL:absorptivity of the wall (-);

αto:absorptivity of the black body (-);

τp:transmittivity of the roof and wall (-);

λair:thermal conductivity of the air (W/(m.K));

μair:dynamical viscosity of the air (Pa.s);

ε:wood stack porosity (-);

ρa:air density (kg/m3);

ρo:wood dry density (kg/m3);

σ:Stefan-Boltzmann constant value (5, 67x10-8W/(m2K4);

∆t:time step (s);

:symbol of the absolute value.


Corresponding author: simotagne2002@yahoo.fr

Received: 17.01.2016 Accepted: 20.11.2016

Creative Commons License Todo el contenido de esta revista, excepto dónde está identificado, está bajo una Licencia Creative Commons