Next Article in Journal
Rapid Truck Loading for Efficient Feedstock Logistics
Previous Article in Journal
Sen2Grass: A Cloud-Based Solution to Generate Field-Specific Grassland Information Derived from Sentinel-2 Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection

by
Merlin Simo-Tagne
1,2,*,
Ablain Tagne Tagne
3,
Macmanus Chinenye Ndukwu
4,
Lyes Bennamoun
5,
Marcel Brice Obounou Akong
3,
Maryam El Marouani
6 and
Yann Rogaume
7
1
INRAE, LERMAB, ERBE—F, 88000 Epinal, France
2
Nancy-Metz Academy, 54035, 2 rue Philippe de Gueldres, 54000 Nancy, France
3
Department of Science, University of Yaoundé I, Yaoundé P.O. Box 812, Cameroon
4
Department of Agricultural and Bioresources Engineering, Michael Okpara University of Agriculture, Umuahia P.M.B. 7267, Abia State, Nigeria
5
Department of Mechanical Engineering, University of New Brunswick, 15 Dineen Drive, Fredericton, NB E3B 5A3, Canada
6
Department of Chemistry, College of Sciences, University of Hafr Al Batin, Hafr Al Batin 39524, Saudi Arabia
7
Université de Lorraine, INRAE, LERMAB, ERBE—F, 88000 Epinal, France
*
Author to whom correspondence should be addressed.
AgriEngineering 2021, 3(1), 138-157; https://doi.org/10.3390/agriengineering3010009
Submission received: 17 February 2021 / Revised: 3 March 2021 / Accepted: 5 March 2021 / Published: 17 March 2021

Abstract

:
In this work, an indirect solar dryer for drying cassava root chips was modelled and experimentally validated using the environmental conditions of Yaoundé in Cameroon and Yamoussoukro in Ivory Coast. The dryers were operational in natural convection mode. Resolution of the equations was achieved by finite differences and the 4th order of Runge–Kutta methods. A model was proposed for performing heat and mass transfer using thermophysical properties of cassava roots, and the obtained results were satisfactory for all conditions, with moisture content difference of less than 0.2 kg/kg between the experimental and theoretical results. The model showed that the core of the product takes more time to dry, which always prolongs the drying duration. The heat and mass transfer coefficients vary during the entire process of solar drying. The drying kinetics vary during the drying with values lower than 1.2 × 10−4 kg/(kg.s). The great gradients of humidity were observed in the thickness of the sample with a regular distribution of the temperature each drying time in the thickness of the sample.

1. Introduction

The challenges facing most developing countries, especially in Africa, is not what to produce but primarily how to process and preserve what is produced [1]. In 2009, Africa alone lost 25% of food produced, which constitute about 4.5 million metric tonnes of food produced in 2009 [2], and the trend continued unabated until this day. This volume of food losses is very huge, considering that most homes in Africa live on less than one dollar a day. This massive loss of food production occurs mostly during the crop season when there is always a glut in the market, and producers or farmers struggle to sell their product. Due to a lack of processing and preservation facilities, these products, if not sold, will be allowed to rot or fed to animals [3]. Cassava root is among the agricultural products that are losses sustained by farmers because cassava (Manihot utilisima Kunz.) root is a major staple food in most African countries and across Latin America and Asia because of their high carbohydrate value [4]. In most countries, cassava roots are processed into flakes, fufu, and cassava chips for direct consumption or utilization in industries for the production of cassava flour or feeds for animals [5]. However, the major critical step employed by farmers to preserve their product is drying. The unit operations involved in cassava processing includes harvesting the roots (Figure 1a) from the farm, washing, cutting, peeling, shredding or chipping (Figure 1b), and drying [6,7,8]. The drying of the cassava chips is essential for shelf life extension, secondary conversion, reduction in bulk weight, product packaging, ease of transportation, etc. However, drying consumes a lot of energy, due to the high initial moisture content of cassava tubers immediately after harvest. Drying can be achieved through the use of artificial industrial dryers or locally by making use of the abundant sun available, as it is done in most African countries [9]. However, the use of sun-drying to dry the chips exposes the cassava chips to contamination from impurities arising from dust, animal droppings, and possible rewetting from rainfall, which prolongs the drying periods [10,11]. In most cases, to safeguard the products from the above scenario, the processor (i.e., person processing the chips) has to be around to keep watch over the product, thereby splitting their time between, going to the farm and the market, and watching over the product, which is usually tedious. Therefore, solar drying has been suggested as the most cost-effective method to solve this problem [12,13,14,15].
Solar dryer gives room for flexibility in the choice of the material and design based on the size of the enterprise and users’ income level. Various designs can range from a simple solar cabinet dryer to indirect solar dryers and hybrids solar dryers equipped with supplementary assisted heat sources, such as thermal storage material, electricity, biomass heater, etc., to continue the drying process during the night or when the weather is not clement [18,19,20,21,22,23,24]. Additionally, Ndukwu et al. [25] and Bennamoun [26] have reviewed different solar dryer designs across the African continent and concluded that different designs are now available in Africa for both direct and indirect drying chamber designs. In addition, Chaudhari and Salve [27], in their review of different solar drying technology and designs, showed that the drying kinetics of several crops based on diffusion models has been studied for solar dryers of different kinds. However, what is lacking in most of the research is simulation models that will help in the future optimization of developed dryers for different products [24]. This will reduce the bulkiness and cost of experimentations and help in the standardization of solar dryer designs. Some researchers have developed some numerical simulation models to tackle these challenges in the African continent. For example, Bahnasawy and Shenanan [28] showed the sudden effect of airflow rate on the initial moisture level when simulating an indirect solar dryer for drying dairy products in Egypt. Bennamoun and Belhamri [29] numerically simulated the solar drying of seedless grapes in Algeria with the outlet temperature dominating the airflow rate. Simo-Tagne et al. [13] presented a numerical model and simulation of solar drying of red chili under the variable external condition of the African continent. Bala and Woods [30] presented a modelling of solar drying of rough rice, while Ayadi [31] considered the three major heat transfer modes in the simulation of solar drying of a medicinal plant in Tunisia with a good result. Other numerical simulations dealing with solar drying of hygroscopic products with good results have been presented also in literature [32,33,34,35].
All the numerical simulations reviewed above showed that the results obtained for each simulation is influenced by the thermophysical properties of the product, the design of the solar dryer and the external conditions of the validation environments, which depends on the seasonal variations. Moreover, literature dealing with the numerical simulation of solar drying of cassava roots is very scarce in the literature. Therefore, the objective of this research is to develop a numerical simulation model for solar drying of cassava roots in a convective mode based on heat and mass transfer, thermophysical properties, and the sorption isotherm of cassava roots. The model will be validated under the external variable condition of Cameroon and Ivory Coast, which are two major producers of cassava chips. The outcome will aid others in producing solar dryer designs and also help designers make informed decisions on design parameters.

2. Materials and Methods

2.1. Description of the Solar Dryer

The dryer was studied without and with the solar collector as presented in Figure 2, referred to as design 1 and design 2, respectively. The solar dryer presented in design 1 was experimentally validated in the Yamoussoukro region of Ivory Coast using cassava roots and was reported in the literature by Tieu et al. [36]. The dryer is an indirect passive solar dryer that is mostly common in Africa due to a lack of electricity [25]. The sides and the bottom faces are made with 10 mm thick plywood while the drying chamber is 1.773 m long, 0.313 m high (north side), and 0.86 m wide. Cassava roots chips are placed on a spread nylon sheet on top of the clay layer inside the drying chamber and arranged to form a rectangular shape that is 1.753 m long and 0.84 m wide. The air passes into the drying chamber through the inlet air section with dimensions of 0.6 m × 0.04 m. The hot and humid air rises by buoyancy and exits through a 1.5 m high chimney with a diameter of 50 mm. The roof made up of blackened sheet metal is inclined at 10° south to capture a lot of sunlight. Figure 2 shows design 2, which has a solar collector attached to form an indirect passive solar dryer. The collector has a dimension of 2 m × 0.6 m × 0.04 m. The top of the solar collector is covered with a glass cover of 10 mm thickness. The absorber is made of a 1 mm thick steel plate painted black. Under the absorber is covered with a wooden panel that serves as thermal insulation. The dryer functions in natural convection, and therefore it is assumed that the air velocity in the solar collector, the drying chamber, and the chimney is close to the ambient air velocity given by the meteorological data of the town during the drying period. The typical weather condition of the Yaoundé environment used for the simulation is shown in Figure 3. Table 1 presents the physical properties of the experimental solar dryer used during the numerical simulations.

2.2. Thermophysical Properties of Cassava Roots Used

The sorption isotherm is very important in drying processes. It gives the equilibrium moisture content of the studied product at the given drying air relative humidity and temperature. The Guggenheim–Anderson–de Boer (GAB) model for equilibrium moisture content is highly recommended in the literature [37] because of its satisfactory results in all relative humidity ranges, allowing for a good approximation of the moisture content at the fiber saturation points. The performance of passive solar dryers is highly dependent on the constant periodic variations of the environmental conditions, and therefore it is important to have a continuous variation between the desorption isotherms and the drying air temperature as well as relative humidity. Thus, we used the GAB model to model the experimental data on desorption isotherms given by Ajala et al. [38]. The model obtained for the simulation is presented in Equation (1) with the parameters given in Table 2.
X eq = X m . G . C . RH 1 C . RH 1 + G . C . RH C . RH
X m = X mo + X mo 1 . T
C = C o exp H C R . T
G = G o exp H G R . T
R = 8.314 J/(mol·K)
From Figure 4, the moisture contents in the first layer were not influenced by the air-drying temperature when the equilibrium moisture content of the product is lower than 0.075 kg/kg. For the equilibrium moisture contents above 0.075, when the drying air temperature increases, the moisture contents at the first layer decrease. A great difference in energies needed to extract one mole of bound water in the studied products ( H C and H G ) is observed. When drying air temperature increases, equilibrium moisture contents decrease showing that the increase of temperature is favorable for the extraction of the bound water in the studied products. The desorption isotherms increase when the air relative humidity increases. When the temperature increases, the influence of the air relative humidity on the desorption isotherms of the products decreases.
Figure 4 presents some desorption isotherm curves of cassava roots. We can see a good agreement between the experimental points and the numerical curves, showing that the GAB model gave a satisfactory result. To obtain all GAB’s parameters, all the squares of the regression coefficients were higher than 0.8 ( r 2 > 0.8 ).
To obtain the specific heat and thermal conductivity of cassava roots, we used the relationships recommended for agriculture products [13]:
C p c a s = 837 + 3348 H 1 + H
λ c a s = 0.49 0.44 e x p 0.206 H
The average density of cassava roots is taken as ρ cas = 840   kg / m 3 . In effect, Baharuddin et al. [39] obtained this value in the range of 810 to 870 kg/m3 according to the maturity age of the cassava roots. The value of dry density used is given by Baharuddin et al. [39] as 670   kg / m 3 . Fernando et al. [40] found that the moisture diffusion coefficient of cassava roots is located between 1.18 × 10−9 to 6.16 × 10−9 m2/s. In this study, we used an Arrhenius variation of the mass diffusion coefficient such as the one used by Ajala et al. [41]:
D H c a s = D o e x p E a R T a
When applied on Cassava, Ajala et al. [41] found activation energy of E a = 30,300   J / mol . This value of activation energy is used for the simulation. The mass diffusivity constant used during the simulation of the drying of pretreated okra by Kuitche et al. [42] was D o = 1.8 × 10 5   m 2 / s . In this work, we used D o = 3.25 × 10 4   m 2 / s in the range of acceptable values obtained on agro-products [42].
The thermal diffusivity is given by:
D T c a s = λ c a s ρ c a s C p c a s

2.3. Heat Mass Transfer on the Products

The heat transfer by radiation between the products and the walls, roof, and the floor of the drying chamber was neglected by studies such as the one by Simo-Tagne et al. [35]. We assume that the tray does not influence heat mass transfer.

2.3.1. Roots of Cassava

We assume that the product has a thin layer configuration. Only the transfers in the thickness orientation are taken into account. Then, we have the following in the product:
T t = D T c a s 2 T x 2
H t = D H c a s 2 H x 2
We also assume that during the drying process, the contact between the clay and the inlet air is considered perfect. Thus, when the solar collector is used (design 2), we have the following:
T x = 0 = T f z = z m a x
H x = 0 = X e q T f s z = z m a x ; R H f s z = z m a x
when the solar collector is not used (design 1), we have:
T x = 0 = T aE
H x = 0 = X eq T aE ; RH aE
At the surface of the product, we have:
λ cas T x x = ep = h c T x = ep T a
D Hcas H x x = ep = h m H x = ep X eq

2.3.2. Drying Chamber

We assume that the temporal variation of the sensible energy of each component of the solar dryer is equal to the sum of all power energy exchanged between the component and all other parts of the dryer. Thus, we have the following:
For the roof:
m r C p r T r t = A r h c e x t T e x t T r + A r F r s k T s k 4 T r 4 + A r h c i n t T a T r + α r A r I G
The walls:
m w C p w T w t = A w h c e x t T e x t T w + A w h c i n t T a T w + α w A w I G
Drying air:
m a C p a T a t = A w h c i n t T w T a + A r h c i n t T r T a + A p h c p T T a + m ˙ C p a T f s T a + m ˙ C p v Y E T f s Y s T a
m a Y s t = m ˙ Y E Y s m 0 H t
The dry mass of the product used during the experiment is equal to m 0 = 14.804   kg .
P v = RH . P vsat T a
The drying kinetics is given by the expression H t .
We also assume that the humidity of the inlet and outlet air is the same in the solar collector. Thus the relative humidity of outlet air (RHfs) is given by:
RH fs = P vsatE P vsatfs RH E
RH E is the relative humidity of ambient air; and P vsatE and P vsatfs are the pressure of vapor at the saturation state of the ambient air and the outlet the solar collector, respectively.

2.3.3. Solar Collector

Glass Cover in Solar Collector

Using the same assumptions taken at the level of drying chamber, we have the following:
m g C p g A g T g c t = α g c 1 + τ g c ρ p I G + h r p g T p c T g c + h c g f T f T g c + h c e x t T a E T g c + h r g c T s k y T g c
Both h rpg and h rgc are taken from Simo-Tagne et al. [43] and given by:
h rpg = σ T pc 2 + T gc 2 T pc + T gc 1 ϵ p + 1 ϵ g 1
h rgc = ϵ g σ T sky 2 + T gc 2 T sky + T gc
where ϵ p = 0.11 , ϵ g = 0.89 and ρ p = 0.06 [13].
The absorber with insulation plate in solar collector is calculated as follows:
m p C pc A pc T pc t = α p τ gc I G 1 ρ p ρ gc + U b T aE T pc + h cpf T f T pc + h rpg T gc T pc
where ρ gc = ρ p = 0.06 [13].
U b = 1 1 h cext + e b λ b + e pc λ pc + 1 h cpf + h rpg
The airflow between absorber plate and glass cover in solar collector is as follows:
m f C pf A f T f t + m ˙ C pf l c T f z = h cpf T pc T f + h cgf T gc T f

2.3.4. Heat and Mass Transfer Coefficients

In the drying chamber, the convective heat transfer coefficients between a solid surface and the air inside are given by Churchill and Chu and taken from Incropera et al. [44]. We have the following:
h c = Nu . λ a d H
h m = Sh . D a d H
The hydraulic diameter is given by:
d H = 4 Vol a Sur a = 0.2481   m
Vola and Sura are respectively the volume and surface of the drying chamber.
The Nusselt number and the Sherwood number were used to obtain the convective heat transfer coefficient and the convective mass transfer coefficient between air drying and cassava, respectively. The Nusselt number was taken from Incropera et al. [44] as follows:
Nu Air Cass = 0.54 Ra 1 / 4
Bower and Saylor [45] gave the Sherwood number as:
Sh Air Cass = 0.23 Sc 1 / 3 Ra 0.321
Ra = Gr . Pr
With:
Sc = ν a D a
ν a = μ a ρ a
The Nusselt number was used to obtain the convective heat transfer coefficient between the wall and the air inside the drying chamber; it is given by Churchill and Chu and was taken from Incropera et al. [44]:
Nu pLa = 0.825 + 0.387 Ra L 1 / 6 [ 1 + 0.492 / Pr 9 / 16 ] 8 / 27 2
The Nusselt number for the convective heat transfer coefficient between the roof and the air inside the drying chamber is also given by Churchill and Chu and taken from Incropera et al. [44] as follows:
Nu ia = 0.677 Pr 1 / 2 Gr 1 / 4 0.95 + Pr 1 / 4
To obtain the Grashof number (Gr), the intensity of gravity g was changed to gcos α , where α is the angle between the roof and the vertical height (80° in this case).
To obtain the convective heat transfer coefficients in the solar collector due to the natural convection between the airflow and the glass cover ( h cgf ) as well as that between the airflow and the absorber ( h cpf ), the Nusselt number was obtained from Simo-Tagne et al. [43]:
N uCol = 1 + 1.44 1 1708 R a cos α + 1 1708 R a cos α 2 1 1708 sin 1.8 α 1.8 R a cos α + R a cos α 5830 1 / 3 1 + R a cos α 5830 1 / 3 1 2        
where α is the angle between the solar collector and the horizontal plane (10°), Ra number was obtained at the average temperature of the airflow and the solid surface.

2.4. Method of Resolution

To resolve Equations (6a)–(6h), the finite differences method was used because it is highly recommended for its stability. Thus, we have the following for Equation (6a,b):
T i j + 1 = C 1 T i j T i j + C 2 T i j ( T i + 1 j + T i 1 j ) H i j + 1 = C 1 H i j H i j + C 2 H i j ( H i + 1 j + H i 1 j )
C 1 T i j = 1 2 D Tcas Δ t Δ x 2 i j
C 2 T i j = D Tcas Δ t Δ x 2 i j
C 1 H i j = 1 2 D Hcas Δ t Δ x 2 i j
C 2 H i j = D Hcas Δ t Δ x 2 i j
Equations (6c) and (6d) give the following:
T 0 j = T f j z = z max
H 0 j = X eq T f j z = z max ; RH f j z = z max
Equations (6e) and (6f) give the following:
T 0 j = T aE j
H 0 j = X eq ( T aE j ; RH aE j )
Equations (6g) and (6h) give respectively the following:
T N + 1 j = ( 1 h c Δ x λ cas ) N + 1 j T N j + h c Δ x λ cas N + 1 j T a j
H N + 1 j = ( 1 h m Δ x D Hcas ) N + 1 j H N j + h m Δ x D Hcas N + 1 j X eq ( T a j ; RH a j )
To resolve Equations (7), (8), (9a), (9b), (10a) and (10d), we used the 4th order Runge–Kutta method. This method has been used with great satisfaction in literature when the partial differential equation is the only function of drying time by Simo-Tagne et al. [43]. To resolve Equation (11), we used the same assumption as Njomo [46] where the sensible heat is neglected. Thus, we have the following:
m ˙ C p f l c T f z = h c p f T p c T f + h c g f T g c T f
With:
T f z = 0 = T a E
Along with the solar collector, the heater of the fresh air coming from ambient is given by the analytical solution:
T f z = h c p f T p c + h c g f T g c h c p f + h c g f + T a E h c p f T p c + h c g f T g c h c p f + h c g f e x p h c p f + h c g f z ρ f e s c V a E C p f
With:
m ˙ = ρ f e s c l c V a E
According to the experimental solar dryer [36], the width ( l c ) and thickness ( e sc ) of the airflow in the solar collector were equal to 0.6 m and 0.04 m, respectively. Table 3 presents the thermal properties of the solar dryer used during the numerical simulations. The thermophysical parameters of the drying air used are presented in Appendix A. Prandtl and Grashop numbers are showed in Appendix B. The time step was equal to 1 s, but the results were recorded each 1 h of drying time. The 2 m length of the solar collector was divided into 250 parts, and the 15 mm thickness of each sample of cassava was divided into 50 parts. We used the Fortran 90 language to have all numerical data.

3. Results and Discussion

3.1. Model Validation in Yamoussoukro (Ivory Coast)

Figure 5 was obtained using the experimental data of Tieu et al. [36] with the solar dryer in design 1, in the town of Yamoussoukro in Ivory Coast from 18 October 2018 to 20 October 2018. All meteorological data for the days of the experiments were presented in the literature [48] and were integrated into the model. The model gave a satisfactory result between the experimental and theoretical results with the difference between experimental data and numerical data of moisture content on a dry basis with less than 0.2 kg/kg (dry basis) moisture content difference. Figure 6 shows that the design influenced the drying kinetic of the product. When the solar collector is used in design 2, the drying air temperature is more heated and increases the mass diffusion coefficients, thus increasing the drying rate. Therefore, the drying curve becomes steeper compared to design 1. However, this variation was not too high because during the night and less sunny days, the solar collector does not have a satisfactory influence. This is the reason that supplementary heat is advocated in most West African states [14].
Figure 7 presents the evolution of the temperatures within the solar drying system in comparison with the ambient temperature. During the day, the exit air from the collector changes the temperature of the drying air in the drying chamber, thereby raising its value. During the night, the temperature of the air in the drying chamber decreases to the ambient air temperature values. This is the major reason the air inlet is closed during the off-sunshine period in order to stop the renewal of the air in the drying chamber as well as to insulate the walls, roof, and floor of the drying chamber to decrease energy loss by conduction or convection.

3.2. Heat Mass Transfer Coefficients, Heat Mass Gradients and Drying Kinetic in Yaoundé (Cameroon)

Figure 8 shows the convective heat and mass transfer coefficients, the temperature of cassava, the relative humidity of drying air, and the moisture content evolutions during the drying of cassava in Yaoundé from 18 November 2018. The meteorological data of the drying periods were presented previously in Figure 3. Figure 8a shows that the heat transfer (about 25 W/m2 s) is by convection between the airflow and absorber, and between the airflow and glass cover during the sunny periods. The maximum values of the heat transfer coefficients between the drying air and the product and the walls are about 4 W/m2s. The global heat transfer coefficient below the solar collector Ub varied around 5 W/m2s while the radiative heat transfer coefficients varied around 0.5 W/m2s. The variation of the convective heat transfer coefficient between the airflow and glass cover of the solar collector was low and was found to be around 25 W/m2s. Figure 8b shows that the solar dryer decreases the value of the drying air relative humidity and increases the temperature during sunny days, which leads to a quick decrease in the moisture content of the product, thereby increasing the mass transfer coefficient (hm) during the sunny days as shown in fig Figure 8c. However, the heat transfer coefficient hm varies between 2.5 × 10−7 and 2 × 10−6 m/s.
Figure 9 shows the effects of temperature and humidity gradient on the thickness of the product (Figure 9a) and the absolute value of the drying kinetics (Figure 9b) as the product loses moisture. From the figures, the first hour of drying time showed that the part of the product in contact with the tray is near the equilibrium moisture content. This is as a result of the solar collector raising the temperature of the airflow. During this drying period, the product loses moisture from the surface as it dries, leaving the center or the core wet. Therefore, drying the core of the product takes more time, which always prolongs the drying duration. After 4 h drying, the moisture content at the core dries from 1.75 to 1 kg/kg. However, a moisture content of about 0.4 kg/kg was achieved at the center after 8 h of drying. High moisture gradients were observed on the side of the tray while on the opposite side, a low moisture gradient was observed due to the variations of the characteristics of the drying air. The temperature distribution in the thickness of the product is regular and increases during the drying process. Thus, the surfaces dried very quickly due to the low relative humidity values at the surfaces.
Figure 9b shows that after 1 h of drying, the surfaces dried very quickly compared to the center of the product. During the drying, the drying kinetics were lower at the surfaces because of low humidity. The values of absolute drying kinetics are below 1.2 × 10−4 kg/(kg·s).

4. Conclusions

An indirect solar dryer for drying cassava root chips was modelled and experimentally validated. The dryers were operational in Sub-Saharan Africa in natural convection mode with Yaoundé and Yamoussoukro chosen as the towns for validation. A resolution of equations was achieved by finite difference and the 4th order Runge–Kutta methods. Models were proposed for heat and mass transfer using thermophysical properties of cassava roots and the obtained results were satisfactory for all conditions for the environmental condition. The moisture content difference between the experimental and theoretical results was less than 0.2 kg/kg dry basis. The model showed that the core of the product takes more time to dry, which always prolongs the drying duration. However, after 4 h drying, the moisture content at the core dries from 1.75 to 1 kg/kg while it took 8 h of drying to lower it to 0.4 kg/kg. In addition, the results show that the heat and mass transfer coefficients varied during the night and the day. The variations of the convective heat transfer coefficients between the drying air and the solar dryer are higher than all the heat radiative transfer coefficients. The high and low temperatures as well as the humidity gradient were observed for each drying time.

Author Contributions

Conceptualization, M.S.-T., A.T.T., and M.C.N.; methodology, M.S.-T. and A.T.T.; software, M.S.-T. and A.T.T.; validation, M.S.-T., A.T.T., and M.C.N.; formal analysis, M.S.-T., A.T.T., and M.C.N.; investigation, M.S.-T., A.T.T., and M.C.N.; resources, M.S.-T., A.T.T., and M.C.N.; data curation, M.S.-T., A.T.T., and M.C.N.; writing—original draft preparation, M.S.-T., A.T.T., and M.C.N.; writing—review and editing, M.S.-T., A.T.T., M.C.N., L.B., M.B.O.A., M.E.M., and Y.R.; visualization, M.S.-T., A.T.T., M.C.N., L.B., M.B.O.A., M.E.M., and Y.R.; supervision, L.B., M.B.O.A., M.E.M., and Y.R.; project administration, M.S.-T. and M.C.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

M. Simo-Tagne and A. Tagne Tagne acknowledge Romain Rémond of the University of Lorraine (ENSTIB-Epinal/France) for the meteorological data obtained from Meteonorm software. LERMAB is supported by a grant overseen by the French National Research Agency (ANR) as part of the “Investissements d’Avenir” program (ANR-11-LABX-0002-01. Lab of Excellence ARBRE) and is part of ICEEL (Institut Carnot pour l’Energie et l’Environnement en Lorraine).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Aarea (m2)
C, G, Co, Goparameters of desorption isotherms (-)
Cpspecific heat capacity (J/kg·K)
D a mass diffusion coefficient of the vapor in drying air (m2/s)
DHmass diffusion coefficient (m2/s)
dHhydraulic diameter (m)
D o parameter of the mass diffusion coefficient (m2/s)
DTthermal diffusion coefficient (m2/s)
Eaenergy of activation (J/mol)
ebthe thickness of the wooden insulation in the solar collection (m)
epcthe thickness of the absorber in the solar collector (m)
e sc the thickness of the airflow in the solar collector (m)
Fa-bgeometrical factor between the areas (a) and (b) (-)
gthe intensity of the presenter (9.81 m2/s)
GrGrashof number (-)
Hmoisture content of the cassava (kg/kg, db)
hcconvective heat transfer coefficient (W/m2.K)
hmconvective mass transfer coefficient (m/s)
hrradiative heat transfer coefficient (W/m2.K)
IGglobal solar irradiation received on the solar dryer (W/m2)
l c width of the airflow in the solar collector (m)
mmass (kg)
modry mass of cassava (kg)
m ˙ mass rate of the drying air (kg/s)
NuNusselt number (-)
ShSherwood number (-)
Ppressure (Pa)
PrPrandtl number (-)
Rperfect gases universal constant (8.314 J/K·mol)
RaRayleigh number (-)
RHrelative humidity of drying air (%/100)
Surarea of the walls of the drying chamber (m2)
Ttemperature (K)
tdrying time (s)
U b global heat transfer coefficient in the bottom of the solar collector (W/m2.K)
V o l volume of the drying chamber (m3)
xaxis coordinates in the thickness direction of the sample (m)
Xmoisture content of cassava at equilibrium state (kg/kg, db)
X m o moisture content of the first layer when T = 0 K
X m o 1 moisture of bound water extracted per kelvin of temperature (kg/kg.K)
Y s humidity of the drying air (kg/kg)
zaxis coordinates of the solar collector (m)
zoaltitude (m)
H C , H G the energy needed to dry a mole of bound water (J/mol)
Greek Symbols
λ thermal conductivity (W/(m. K))
ρ density (kg/m3)
μ a dynamic viscosity of drying air (Pa·s)
ν a kinematic viscosity of the drying air (m2/s)
σ Stephan Boltzmann coefficient (5.67 × 10−8)
α the angle between the solar collector and the horizontal plane (°)
α i solar absorptivity coefficient of component i (-)
τ i solar transmissivity coefficient of component i (-)
ρ p solar reflexivity coefficient of the absorber (-)
ϵ i solar emissivity coefficient of component i (-)
Subscripts
adrying air in the drying chamber or solar collector
aE, E, extambient (exterior space of the solar drying)
airair
atmatmosphere
cascassava
epsurface
eqequilibrium
fairflow in the solar collector
fsend position of the solar collector when z = zmax (outlet point)
ggas (air)
gc, gglass cover
ithe interior of the dryer
Lalateral part
mthe first layer of bound water
pcabsorber in the solar collector
rroof
sksky
vvapor
vsatvapor in saturation state
wwall
Index
i,jorder of discretization scheme (-)

Appendix A. Thermophysical Parameters of the Drying Air

Appendices See [2,34,49].
Thermal conductivity (W/(m. K)), (Ta in °C):
λ a = 7.57 × 10 5 T a + 0.0242
Heat capacity of dry air (J/(kg. K)), (Ta in °C):
C pa = 1835 0.734 T a
Heat capacity of vapor air (J/(kg·K)), (Ta in °C):
C pv = 999.2 + 0.1434 T a + 1.101 × 10 4 T a 2 6.7581 × 10 8 T a 3
Dynamic viscosity (Pa·s), (Ta in °C):
μ a = 1.718 × 10 5 + 4.620 × 10 8 T a
Mass diffusion coefficient of vapor in air:
D a m 2 / s = 2.25 × 10 5 P atm T a 273 1.81
P atm = 101,325 12 zo + 5.2 × 10 4 z o 2
where zo in m is the altitude, T a is given in K, and Patm in Pa.
The density of the air vapor:
ρ g kg / m 3 = 353 T a
where T a is given in K.

Appendix B. Non-Dimensional Numbers

Prandtl number (-) (Ta in °C):
Pr = 2.54 × 10 4 T a + 0.7147
Grashof number (-) (Ta in K):
Gr = β g d h 3 ρ a 2 Δ T μ air 2   ;   β = 1 T a

References

  1. Ndukwu, M.C.; Akani, O.A.; Simonyan, K.J. Nigeria’s grain resource structure, and government sustainable policy: A review. Agric. Eng. Int. CIGR J. 2015, 17, 441–457. [Google Scholar]
  2. Simo-Tagne, M.; Ndukwu, M.C.; Zoulalian, A.; Bennamoun, L.; Kifani-Sahban, F.; Rogaume, Y. Numerical analysis and validation of a natural convection mix-mode solar dryer for drying red chilli under variable conditions. Renew. Energy 2020, 151, 659–673. [Google Scholar] [CrossRef]
  3. Ndukwu, M.C.; Manuwa, S.I. Impact of evaporative cooling preservation on the shelf life of fruits and vegetable in South-Western Nigeria. Res. Agric. Eng. 2015, 61, 122–128. [Google Scholar] [CrossRef] [Green Version]
  4. Ndukwu, M.C.; Asoegwu, S.N.; Ahaneku, I.E. Status of research on agricultural machinery development in Nigeria: A case study of cassava tuber processing machineries. Agric. Mech. Asia Afr. Latin Am. 2018, 49, 150–155. [Google Scholar]
  5. Ndukwu, M.C.; Onyenwigwe, D. Development of a motorized parboiled cassava tuber shredding machine. Niger. J. Technol. 2013, 32, 1–6. [Google Scholar]
  6. Egbeocha, C.C.; Asoegwu, S.N.; Okereke, N.A.A. A Review on performance of cassava, peeling machines in Nigeria. Futo J. Ser. 2016, 2, 140–168. [Google Scholar]
  7. Igboayaka, E.C.; Ndukwu, M.C.; Ernest, I.C. A Modelling approach for determining the throughput capacity and energy consumption of a cassava tuber shredder. J. Chin. Adv. Mater. Soc. 2018, 6, 801–816. [Google Scholar] [CrossRef]
  8. Ndukwu, M.C.; Afam, G.; Nwakuba, N.R. Development and optimization of a manual fed cassava root chipper for household cassava processors. Turk. J. Agric. Eng. Res. (TURKAGER) 2020, 1, 283–295. [Google Scholar] [CrossRef]
  9. Ndukwu, M.C.; Ogunlowo, A.S.; Olukunle, O.J. Cocoa bean (Theobroma Cacao L.) drying kinetics. Chil. J. Agric. Res. 2010, 70, 633–639. [Google Scholar]
  10. Sharma, A.; Chen, C.R.; Vu Lan, N. Solar-energy drying systems: A review. Renew. Sustain. Energy Rev. 2009, 13, 1185–1210. [Google Scholar] [CrossRef]
  11. Ndukwu, M.C.; Bennamoun, L. Potential of integrating Na2SO4.10H2O pellets in solar drying system. Dry. Technol. 2018, 36, 1017–1030. [Google Scholar] [CrossRef]
  12. Anyanwu, C.N.; Oparaku, O.U.; Onyegegbu, S.O.; Egwuatu, U.; Edem, N.I.; Egbuka, K.; Nwosu, P.N.; Sharma, V.K. Experimental investigation of a photo voltaic-powered solar cassava dryer. Dry. Technol. 2012, 30, 398–403. [Google Scholar] [CrossRef]
  13. Simo-Tagne, M. Simple solar dryers adapted to the social, climatic and economic context can help increase the value of tropical timber. ITTO Trop. For. Update 2018, 27, 25–27. [Google Scholar]
  14. Ndukwu, M.C.; Diemuodeke, E.O.; Abam, F.I.; Abada, U.C.; Eke-emezie, N.; Simo-Tagne, M. Development and modelling of heat and mass transfer analysis of a low-cost solar dryer integrated with biomass heater: Application for West African Region. Sci. Afr. 2020, 10, e00615. [Google Scholar] [CrossRef]
  15. Ndukwu, M.C.; Simo-Tagne, M.; Bennamoun, L. Solar drying research of medicinal and aromatic plants: An African experience with assessment of the economic and environmental impact. Afr. J. Sci. Technol. Innov. Dev. 2020, 1–14. [Google Scholar] [CrossRef]
  16. Cadyventures. Available online: https://www.cadyventures.com/product/organic-dried-cassava-chips/ (accessed on 14 February 2021).
  17. Istockphoto. Available online: https://www.istockphoto.com/fr/search/2/image?excludenudity=false&phrase=cassava%20roots (accessed on 14 February 2021).
  18. Pirasteh, G.; Saidur, R.; Rahman, S.M.A.; Rahim, N.A. A review on development of solar drying applications. Renew. Sustain. Energy Rev. 2014, 31, 133–148. [Google Scholar] [CrossRef]
  19. Fudholi, A.; Sopian, K.; Bakhtyar, B.; Gabbasa, M.; Othman, M.Y.; Ruslan, M.H. Review of solar drying systems with air-based solar collectors in Malaysia. Renew. Sustain. Energy Rev. 2015, 51, 1191–1204. [Google Scholar] [CrossRef]
  20. Hossain, M.Z.; Alam, M.M.; Hossain, M.F.; Sarker, M.S. Performance evaluation of a cabinet solar dryer for drying red pepper in Bangladesh. J. Agric. Eng. 2018, 774, 100–109. [Google Scholar] [CrossRef]
  21. Ren, H.; Ma, Z.; Lin, W.; Fan, W.; Li, W. Integrating photovoltaic thermal collectors and thermal energy storage systems using phase change materials with rotary desiccant cooling systems. Sustain. Cities Soc. 2018, 36, 131–143. [Google Scholar] [CrossRef] [Green Version]
  22. Amer, B.M.A.; Gottschalk, K.; Hossain, M.A. Integrated hybrid solar drying system and its drying kinetics of chamomile. Renew. Energy 2018, 121, 539–547. [Google Scholar] [CrossRef]
  23. Nwakuba, N.R.; Ndukwu, M.C.; Asonye, G.U.; Asoegwu, S.N.; Nwandikom, G.I. Environmental sustainability analysis of a hybrid heat source dryer. Polytechnica 2020, 9–114. [Google Scholar] [CrossRef]
  24. Simo-Tagne, M.; Ndukwu, M.C.; Ndi-Azese, M. Experimental modelling of a solar dryer for wood fuel in Epinal (France). Modelling 2020, 1, 39–52. [Google Scholar] [CrossRef]
  25. Ndukwu, M.C.; Bennamoun, L.; Abam, F.I. Experience of solar drying in Africa: Presentation of designs, operations, and models. Food Eng. Rev. 2018, 10, 211–244. [Google Scholar] [CrossRef]
  26. Bennamoun, L. Reviewing the experience of solar drying in Algeria with presentation of the different design aspects of solar dryers. Renew. Sustain. Energy Rev. 2011, 15, 3371–3379. [Google Scholar] [CrossRef]
  27. Chaudhari, A.D.; Salve, S.P. A review of solar dryer technologies. Int. J. Res. Advent Technol. 2014, 2, 218–232. [Google Scholar]
  28. Bahnasawy, A.H.; Shenana, M.E. Mathematical model of direct sun and solar drying of some fermented dairy products (Kishk). J. Food Eng. 2004, 61, 309–319. [Google Scholar] [CrossRef]
  29. Bennamoun, L.; Belhamri, A. Numerical simulation of drying under variable external conditions: Application to solar drying of seedless grapes. J. Food Eng. 2006, 76, 179–187. [Google Scholar] [CrossRef]
  30. Bala, B.K.; Woods, J.L. Simulation of the indirect natural convection solar drying of rough rice. Sol. Energy 1994, 53, 259–266. [Google Scholar] [CrossRef]
  31. Ayadi, M.; Zouari, I.; Bellagi, A. Simulation, and performance of a solar drying unit with storage for aromatic and medicinal plants. Int. J. Food Eng. 2015, 11, 597–607. [Google Scholar] [CrossRef]
  32. Lamrani, B.; Khouya, A.; Draoui, A. Energy and environmental analysis of an indirect hybrid solar dryer of wood using TRNSYS software. Sol. Energy 2019, 183, 132–145. [Google Scholar] [CrossRef]
  33. Bentayeb, F.; Bekkioui, N.; Zeghmati, B. Modeling and simulation of a wood solar dryer in a Moroccan climate. Renew. Energy 2008, 33, 501–506. [Google Scholar] [CrossRef]
  34. Simo-Tagne, M.; Bonoma, B.; Bennamoun, L.; Monkam, L.; Leonard, A.; Zoulalian, A.; Rogaume, Y. Modeling of coupled heat and mass transfer during drying of ebony wood using indirect natural convection solar dryer. Dry. Technol. 2019, 37, 1863–1878. [Google Scholar] [CrossRef]
  35. Simo-Tagne, M.; Bennamoun, L.; Léonard, A.; Rogaume, Y. Modeling, numerical simulation and validation of a convective dryer in steady conditions: Case study of tropical woods. Int. J. Model. Simul. 2020, 40, 143–161. [Google Scholar] [CrossRef]
  36. Tieu, Z.A.; Gbaha, P.; Diby, K.A. Étude expérimentale d’un grenier séchoir solaire à convection naturelle: Application au séchage du cacao et du manioc. Afr. Sci. 2019, 15, 80–95. [Google Scholar]
  37. Simo-Tagne, M.; Zoulalian, A.; Rogaume, Y.; Rémond, R.; Bonoma, B. Modélisation des isothermes de sorption, caractérisation des propriétés thermodynamiques et détermination des humidités d’équilibre d’usage des bois tropicaux. Rev. Energies Renouvelables 2016, 19, 79–96. [Google Scholar]
  38. Ajala, A.S.; Ngoddy, P.O.; Olajide, J.O. Sorption isotherms and their fitted equations for dried chips of cassava roots (Manihot Esculenta Crantz; Tme-7 variety) and the resulting isosteric heats of sorption. Food Res. 2020, 4, 703–711. [Google Scholar] [CrossRef]
  39. Baharuddin, N.H.; Mohamed, M.; Abdullah, M.M.A.B.; Muhammad, N.; Rahman, R.; Omar, M.N.; Amini, M.H.M.; Razab, M.K.A.A.; Rizman, Z.I. Potential of cassava root as a raw material for biocomposite development. ARPN J. Eng. Appl. Sci. 2016, 11, 6138–6147. [Google Scholar]
  40. Fernando, W.J.N.; Low, H.C.; Ahmad, A.L. The effect of infrared on diffusion coefficients and activation energies in convective drying: A case study for banana, cassava and pumpkin. J. Appl. Sci. 2011, 11, 3635–3639. [Google Scholar] [CrossRef] [Green Version]
  41. Ajala, A.S.; Aboiye, A.O.; Popoola, J.O.; Adeyanju, J.A. Drying characteristics and mathematical modelling of cassava chips. Chem. Process Eng. Res. 2012, 4, 1–9. [Google Scholar]
  42. Kuitche, A.; Edoun, M.; Takamte, G. Influence of pre-treatment on drying on the drying kinetic of a local okra (Hibiscus ersculentus) variety. World J. Dairy Food Sci. 2007, 2, 83–88. [Google Scholar]
  43. Simo-Tagne, M.; Zoulalian, A.; Rémond, R.; Rogaume, Y. Mathematical modelling and numerical simulation of a simple solar dryer for tropical wood using a collector. Appl. Therm. Eng. 2018, 131, 356–369. [Google Scholar] [CrossRef]
  44. Incropera, F.P.; Dewitt, D.P.; Bergman, T.L.; Lavine, A.S. Fundamentals of Heat and Mass Transfer; John Wiley and Sons: Hoboken, NJ, USA, 2007; 1070p. [Google Scholar]
  45. Bower, S.M.; Saylor, J.R. A study of the Sherwood–Rayleigh relation for water undergoing natural convection-driven evaporation. Int. J. Heat Mass Transf. 2009, 52, 3055–3063. [Google Scholar] [CrossRef]
  46. Njomo, D. Modelling the heat exchanges in a solar air heater with a cover partially transparent to infrared radiation. Energy Convers. Manag. 1991, 31, 495–503. [Google Scholar] [CrossRef]
  47. Simo-Tagne, M.; Zoulalian, A.; Remond, R.; Rogaume, Y.; Bonoma, B. Modeling and simulation of an industrial indirect solar dryer for iroko wood (Chlorophora excelsa) in a tropical environnement. Maderas Cienc. Tecnol. 2017, 19, 95–112. [Google Scholar] [CrossRef]
  48. Meteonorm. Available online: http://www.meteonorm.com/ (accessed on 24 June 2019).
  49. Simo-Tagne, M.; Rogaume, Y.; Zoulalian, A. Modeling and simulation of the drying of the beech timber (Fagus sylvatica) using oscillating regimes. Int. J. Therm. Environ. Eng. 2014, 8, 83–90. [Google Scholar]
Figure 1. Some pictures of cassava roots. (a) Cassava roots with bark from plantations [16], (b) cassava roots chips without bark ready to be dried [17].
Figure 1. Some pictures of cassava roots. (a) Cassava roots with bark from plantations [16], (b) cassava roots chips without bark ready to be dried [17].
Agriengineering 03 00009 g001
Figure 2. Schematic diagrams of the studied solar dryer without (design 1) and with (design2) a solar collector.
Figure 2. Schematic diagrams of the studied solar dryer without (design 1) and with (design2) a solar collector.
Agriengineering 03 00009 g002
Figure 3. Ambient temperature as well as global and direct solar irradiation on a horizontal plane in Yaoundé from 22 November 2018 (taken from http://www.meteonorm.com/).
Figure 3. Ambient temperature as well as global and direct solar irradiation on a horizontal plane in Yaoundé from 22 November 2018 (taken from http://www.meteonorm.com/).
Agriengineering 03 00009 g003
Figure 4. Validation of four curves of desorption isotherms of cassava roots using the parameters of Table 2 at temperatures of 53, 60, 80, and 86 °C.
Figure 4. Validation of four curves of desorption isotherms of cassava roots using the parameters of Table 2 at temperatures of 53, 60, 80, and 86 °C.
Agriengineering 03 00009 g004
Figure 5. Validation of the model using experimental data given by Tieu et al. [36] in Yamoussoukro, October 2018, design 1.
Figure 5. Validation of the model using experimental data given by Tieu et al. [36] in Yamoussoukro, October 2018, design 1.
Agriengineering 03 00009 g005
Figure 6. Influence of the design on the drying kinetic in Yamoussoukro, October 2018, designs 1 and 2.
Figure 6. Influence of the design on the drying kinetic in Yamoussoukro, October 2018, designs 1 and 2.
Agriengineering 03 00009 g006
Figure 7. Evolutions of temperatures of ambient air (TaE), drying air (Ta), of the product (T) and airflow air outlet the solar collector (Tfs) in Yamoussoukro, October 2018, design 2.
Figure 7. Evolutions of temperatures of ambient air (TaE), drying air (Ta), of the product (T) and airflow air outlet the solar collector (Tfs) in Yamoussoukro, October 2018, design 2.
Agriengineering 03 00009 g007
Figure 8. Drying kinetics in Yaoundé (Cameroon) from 22 November 2018, at noon (design 2). Heat transfer coefficients (a), Temperature, relative humidity and moisture content variations (b), mass transfer coefficient (c).
Figure 8. Drying kinetics in Yaoundé (Cameroon) from 22 November 2018, at noon (design 2). Heat transfer coefficients (a), Temperature, relative humidity and moisture content variations (b), mass transfer coefficient (c).
Agriengineering 03 00009 g008aAgriengineering 03 00009 g008b
Figure 9. The distribution of the temperature and moisture content in the thickness of the product for 18 October; initial local drying time was 8:00 am at Yaoundé. Temperature and moisture content (a) and drying kinetic (b) using design 2.
Figure 9. The distribution of the temperature and moisture content in the thickness of the product for 18 October; initial local drying time was 8:00 am at Yaoundé. Temperature and moisture content (a) and drying kinetic (b) using design 2.
Agriengineering 03 00009 g009
Table 1. Physical properties and characteristics of the components of the solar dryer [13].
Table 1. Physical properties and characteristics of the components of the solar dryer [13].
ComponentsAreas (m2)Density (kg/m3)Mass (kg)
Drying chamber
Roof1.55027205.16
Wall0.8244503.709
Cassava roots chips1.473840 (dry)
Air drying0.2386 (volume in m3) ρ a = 353 / T a 0.2386 × ρ a
Solar collector
Glass cover1.20025308.36
Absorber1.20027205.64
Wooden insulation1.2004505.4
Airflow0.6 × 0.04 (section) ρ f = 353 / T f 0.0288 × ρ f
Table 2. Parameters desorption isotherms model of the GAB model obtained using data taken from Ajala et al. [38].
Table 2. Parameters desorption isotherms model of the GAB model obtained using data taken from Ajala et al. [38].
ParametersCassava Roots
X m o (kg/kg)0.20444366
X m o 1 (kg/(kg.K))−0.00034997
C o (-)0.1809885
H C (J/mol)1917.2676
G o (-)6632.002
H G (J/mol)−11,219.687
Table 3. Thermal properties of the components of solar dryer [43,47].
Table 3. Thermal properties of the components of solar dryer [43,47].
ComponentsSpecific Heat Capacity (J/(kg·K))Thermal Conductivity (W/(m. K))Absorptivity Coefficient (-)Transmitivity Coefficient (-)
For the drying chamber
Roof9001900.910
Wall16000.130.40
Cassava roots chipsSee Equation (2)See Equation (3)(-)(-)
Air dryingSee Equations (A2) and (A3)See Equation (A1)(-)(-)
For the solar collector
Glass cover8000.13110.10.89
Absorber9001900.910
Wooden insulation12360.14(-)(-)
AirflowSee Equations (A2) and (A3)See Equation (A1)0.10.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Simo-Tagne, M.; Tagne Tagne, A.; Ndukwu, M.C.; Bennamoun, L.; Obounou Akong, M.B.; El Marouani, M.; Rogaume, Y. Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection. AgriEngineering 2021, 3, 138-157. https://doi.org/10.3390/agriengineering3010009

AMA Style

Simo-Tagne M, Tagne Tagne A, Ndukwu MC, Bennamoun L, Obounou Akong MB, El Marouani M, Rogaume Y. Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection. AgriEngineering. 2021; 3(1):138-157. https://doi.org/10.3390/agriengineering3010009

Chicago/Turabian Style

Simo-Tagne, Merlin, Ablain Tagne Tagne, Macmanus Chinenye Ndukwu, Lyes Bennamoun, Marcel Brice Obounou Akong, Maryam El Marouani, and Yann Rogaume. 2021. "Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection" AgriEngineering 3, no. 1: 138-157. https://doi.org/10.3390/agriengineering3010009

APA Style

Simo-Tagne, M., Tagne Tagne, A., Ndukwu, M. C., Bennamoun, L., Obounou Akong, M. B., El Marouani, M., & Rogaume, Y. (2021). Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection. AgriEngineering, 3(1), 138-157. https://doi.org/10.3390/agriengineering3010009

Article Metrics

Back to TopTop