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

Merlin Simo-Tagne, Ablain Tagne Tagne, Macmanus Chinenye Ndukwu, Lyes Bennamoun, Marcel Brice Obounou Akong, Maryam El Marouani, Yann Rogaume

► To cite this version:

Merlin Simo-Tagne, Ablain Tagne Tagne, Macmanus Chinenye Ndukwu, Lyes Bennamoun, Marcel Brice Obounou Akong, et al.. Numerical Study of the Drying of Cassava Roots Chips Using an Indirect Solar Dryer in Natural Convection. AgriEngineering, 2021, 3 (1), pp.138-157. 10.3390/agriengineering3010009. hal-03533349

HAL Id: hal-03533349

https://hal.univ-lorraine.fr/hal-03533349

Submitted on 10 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

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

Merlin Simo-Tagne<sup>1,2,*</sup> (ID), Ablain Tagne Tagne<sup>3</sup>, Macmanus Chinenye Ndukwu<sup>4 ID</sup>, Lyes Bennamoun<sup>5 ID</sup>, Marcel Brice Obounou Akong<sup>3</sup>, Maryam El Marouani<sup>6 ID</sup> and Yann Rogaume<sup>7 ID</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** 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–19. https://doi.org/10.3390/agriengineering3010009

Received: 17 February 2021

Accepted: 5 March 2021

Published: 17 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**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 \times 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.

**Keywords:** cassava roots; desorption isotherms; drying kinetics; heat and mass transfer; natural convection; solar dryer; Sub-Saharan Africa

## 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

<footnote n="1">1 INRAE, LERMAB, ERBE—F, 88000 Epinal, France</footnote>

<footnote n="2">2 Nancy-Metz Academy, 54035, 2 rue Philippe de Gueldres, 54000 Nancy, France</footnote>

<footnote n="3">3 Department of Science, University of Yaoundé I, Yaoundé P.O. Box 812, Cameroon; tagneablain97@gmail.com (A.T.T.); marcelobounou@yahoo.fr (M.B.O.A.)</footnote>

<footnote n="4">4 Department of Agricultural and Bioresources Engineering, Michael Okpara University of Agriculture, Umuahia P.M.B. 7267, Abia State, Nigeria; ndukwumcu@mouau.edu.ng</footnote>

<footnote n="5">5 Department of Mechanical Engineering, University of New Brunswick, 15 Dineen Drive, Fredericton, NB E3B 5A3, Canada; lyes.bennamoun@unb.ca</footnote>

<footnote n="6">6 Department of Chemistry, College of Sciences, University of Hafr Al Batin, Hafr Al Batin 39524, Saudi Arabia; melmarouani@uhb.edu.sa</footnote>

<footnote n="7">7 Université de Lorraine, INRAE, LERMAB, ERBE—F, 88000 Epinal, France; yann.rogaume@univ-lorraine.fr</footnote>

<footnote>* Correspondence: simotagne2002@yahoo.fr or Merlin.Simo-Tagne@ac-nancy-metz.fr; Tel.: +33-644-911-921</footnote>
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–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–15].

<figure><img src="image_4.png" /><figcaption>**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].</figcaption></figure>

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–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–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.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Schematic diagrams of the studied solar dryer without (design 1) and with (design2) a solar collector.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**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/).</figcaption></figure>
## 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_{\text{eq}} = \frac{X_{\text{m}} \cdot \text{G.C.RH}}{(1 - \text{C.RH})(1 + \text{G.C.RH} - \text{C.RH})} \quad (1\text{a})$$

$$X_{\text{m}} = X_{\text{mo}} + X_{\text{mo1}} \cdot T \quad (1\text{b})$$

$$C = C_o \exp\left(\frac{\Delta H_{\text{C}}}{R \cdot T}\right) \quad (1\text{c})$$

$$G = G_o \exp\left(\frac{\Delta H_{\text{G}}}{R \cdot T}\right) \quad (1\text{d})$$

$$R = 8.314 \text{ J/(mol}\cdot\text{K)} \quad (1\text{e})$$

<table><caption>**Table 1.** Physical properties and characteristics of the components of the solar dryer [13].</caption><thead><tr><th>Components</th><th>Areas (m<sup>2</sup>)</th><th>Density (kg/m<sup>3</sup>)</th><th>Mass (kg)</th></tr></thead><tbody><tr><td colspan="4">Drying chamber</td></tr><tr><td>Roof</td><td>1.550</td><td>2720</td><td>5.16</td></tr><tr><td>Wall</td><td>0.824</td><td>450</td><td>3.709</td></tr><tr><td>Cassava roots chips</td><td>1.473</td><td>840 (dry)</td><td></td></tr><tr><td>Air drying</td><td>0.2386 (volume in m<sup>3</sup>)</td><td>ρ<sub>a</sub> = 353/T<sub>a</sub></td><td>0.2386 × ρ<sub>a</sub></td></tr><tr><td colspan="4">Solar collector</td></tr><tr><td>Glass cover</td><td>1.200</td><td>2530</td><td>8.36</td></tr><tr><td>Absorber</td><td>1.200</td><td>2720</td><td>5.64</td></tr><tr><td>Wooden insulation</td><td>1.200</td><td>450</td><td>5.4</td></tr><tr><td>Airflow</td><td>0.6 × 0.04 (section)</td><td>ρ<sub>f</sub> = 353/T<sub>f</sub></td><td>0.0288 × ρ<sub>f</sub></td></tr></tbody></table>

<table><caption>**Table 2.** Parameters desorption isotherms model of the GAB model obtained using data taken from Ajala et al. [38].</caption><thead><tr><th>Parameters</th><th>Cassava Roots</th></tr></thead><tbody><tr><td>X<sub>mo</sub>(kg/kg)</td><td>0.20444366</td></tr><tr><td>X<sub>mo1</sub>(kg/(kg.K))</td><td>-0.00034997</td></tr><tr><td>C<sub>o</sub>(-)</td><td>0.1809885</td></tr><tr><td>ΔH<sub>C</sub>(J/mol)</td><td>1917.2676</td></tr><tr><td>G<sub>o</sub>(-)</td><td>6632.002</td></tr><tr><td>ΔH<sub>G</sub>(J/mol)</td><td>-11,219.687</td></tr></tbody></table>
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<sub>C</sub> and ΔH<sub>G</sub>) 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² > 0.8).

<figure><img src="image_7.png" /><figcaption>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.</figcaption></figure>

To obtain the specific heat and thermal conductivity of cassava roots, we used the relationships recommended for agriculture products [13]:

$$C_{\text{pcas}} = 837 + 3348 \left( \frac{H}{1+H} \right) \qquad (2)$$

$$\lambda_{\text{cas}} = 0.49 - 0.44\exp(-0.206H) \qquad (3)$$

The average density of cassava roots is taken as ρ<sub>cas</sub> = 840 kg/m³. In effect, Baharuddin et al. [39] obtained this value in the range of 810 to 870 kg/m³ 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³. Fernando et al. [40] found that the moisture diffusion coefficient of cassava roots is located between 1.18 × 10<sup>-9</sup> to 6.16 × 10<sup>-9</sup> m²/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_{\text{Hcas}} = D_{\text{o}}\exp\left(-\frac{E_{\text{a}}}{RT_{\text{a}}}\right) \qquad (4)$$

When applied on Cassava, Ajala et al. [41] found activation energy of E<sub>a</sub> = 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 \times 10^{-5} \text{ m}^2/\text{s}$. In this work, we used $D_o = 3.25 \times 10^{-4} \text{ m}^2/\text{s}$ in the range of acceptable values obtained on agro-products [42].

The thermal diffusivity is given by:

$$D_{Tcas} = \frac{\lambda_{cas}}{\rho_{cas}C_{pcas}} \qquad (5)$$

## 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:

$$\frac{\partial T}{\partial t} = D_{Tcas} \frac{\partial^2 T}{\partial x^2} \qquad (6a)$$

$$\frac{\partial H}{\partial t} = D_{Hcas} \frac{\partial^2 H}{\partial x^2} \qquad (6b)$$

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_{max}) \qquad (6c)$$

$$H(x = 0) = X_{eq}(T_{fs}(z = z_{max}); RH_{fs}(z = z_{max})) \qquad (6d)$$

when the solar collector is not used (design 1), we have:

$$T(x = 0) = T_{aE} \qquad (6e)$$

$$H(x = 0) = X_{eq}(T_{aE}; RH_{aE}) \qquad (6f)$$

At the surface of the product, we have:

$$- \lambda_{cas} \left. \frac{\partial T}{\partial x} \right|_{x=ep} = h_c (T_{x=ep} - T_a) \qquad (6g)$$

$$- D_{Hcas} \left. \frac{\partial H}{\partial x} \right|_{x=ep} = h_m (H_{x=ep} - X_{eq}) \qquad (6h)$$

### 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_{pr} \frac{\partial T_r}{\partial t} = A_r h_{cext} (T_{ext} - T_r) + A_r F_{r-sk} (T_{sk}^4 - T_r^4) + A_r h_{cint} (T_a - T_r) + \alpha_r A_r I_G \qquad (7)$$

The walls:

$$m_w C_{pw} \frac{\partial T_w}{\partial t} = A_w h_{cext} (T_{ext} - T_w) + A_w h_{cint} (T_a - T_w) + \alpha_w A_w I_G \qquad (8)$$
Drying air:

$$m_a C_{pa} \frac{\partial T_a}{\partial t} = A_w h_{cint} (T_w - T_a) + A_r h_{cint} (T_r - T_a) + A_p h_{cp} (T - T_a) + \dot{m} C_{pa} (T_{fs} - T_a) + \dot{m} C_{pv} (Y_E T_{fs} - Y_s T_a) \quad (9a)$$

$$m_a \frac{\partial Y_s}{\partial t} = \dot{m} (Y_E - Y_s) - m_0 \frac{\partial H}{\partial t} \quad (9b)$$

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) \quad (9c)$$

The drying kinetics is given by the expression $\frac{\partial H}{\partial 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 ($RH_{fs}$) is given by:

$$RH_{fs} = \frac{P_{vsatE}}{P_{vsatfs}} RH_E \quad (9d)$$

$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:

$$\frac{m_g C_{pg} \partial T_{gc}}{A_g \partial t} = \alpha_{gc} (1 + \tau_{gc} \rho_p) I_G + h_{rpg} (T_{pc} - T_{gc}) + h_{cgf} (T_f - T_{gc}) + h_{cext} (T_{aE} - T_{gc}) + h_{rgc} (T_{sky} - T_{gc}) \quad (10a)$$

Both $h_{rpg}$ and $h_{rgc}$ are taken from Simo-Tagne et al. [43] and given by:

$$h_{rpg} = \frac{\sigma(T_{pc}^2 + T_{gc}^2)(T_{pc} + T_{gc})}{\frac{1}{\epsilon_p} + \frac{1}{\epsilon_g} - 1} \quad (10b)$$

$$h_{rgc} = \epsilon_g \sigma (T_{sky}^2 + T_{gc}^2) (T_{sky} + T_{gc}) \quad (10c)$$

where $\epsilon_p = 0.11$, $\epsilon_g = 0.89$ and $\rho_p = 0.06$ [13].

The absorber with insulation plate in solar collector is calculated as follows:

$$\frac{m_p C_{pc} \partial T_{pc}}{A_{pc} \partial t} = \frac{\alpha_p \tau_{gc} I_G}{1 - \rho_p \rho_{gc}} + U_b (T_{aE} - T_{pc}) + h_{cpf} (T_f - T_{pc}) + h_{rpg} (T_{gc} - T_{pc}) \quad (10d)$$

where $\rho_{gc} = \rho_p = 0.06$ [13].

$$U_b = \frac{1}{\frac{1}{h_{cext}} + \frac{e_b}{\lambda_b} + \frac{e_{pc}}{\lambda_{pc}} + \frac{1}{h_{cpf} + h_{rpg}}} \quad (10e)$$

The airflow between absorber plate and glass cover in solar collector is as follows:

$$\frac{m_f C_{pf} \partial T_f}{A_f \partial t} + \frac{\dot{m} C_{pf} \partial T_f}{l_c \partial z} = h_{cpf} (T_{pc} - T_f) + h_{cgf} (T_{gc} - T_f) \quad (11)$$
### 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 = \frac{Nu.\lambda_a}{d_H} \qquad (12a)$$

$$h_m = \frac{Sh.D_a}{d_H} \qquad (12b)$$

The hydraulic diameter is given by:

$$d_H = \frac{4Vol_a}{Sur_a} = 0.2481 \text{ m} \qquad (12c)$$

$Vol_a$ and $Sur_a$ 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.54Ra^{1/4} \qquad (12d)$$

Bower and Saylor [45] gave the Sherwood number as:

$$Sh_{Air-Cass} = 0.23Sc^{1/3}Ra^{0.321} \qquad (12e)$$

$$Ra = Gr.Pr \qquad (12f)$$

With:

$$Sc = \frac{\nu_a}{D_a} \qquad (12g)$$

$$\nu_a = \frac{\mu_a}{\rho_a} \qquad (12h)$$

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} = \left[ 0.825 + \frac{0.387Ra_L^{1/6}}{[1 + (0.492/Pr)^{9/16}]^{8/27}} \right]^2 \qquad (12i)$$

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} = \frac{0.677Pr^{1/2}Gr^{1/4}}{[0.95 + Pr]^{1/4}} \qquad (12j)$$

To obtain the Grashof number (Gr), the intensity of gravity $g$ was changed to $g\cos\alpha$, where $\alpha$ 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_{cfg}$) as well as that between the airflow and the absorber ($h_{cpf}$), the Nusselt number was obtained from Simo-Tagne et al. [43]:

$$Nu_{Col} = 1 + 1.44 \left[ \frac{1 - \frac{1708}{R_a \cos \alpha} + 1 - \frac{1708}{R_a \cos \alpha}}{2} \right] \left[ 1 - \frac{1708(\sin(1.8\alpha))^{1.8}}{R_a \cos \alpha} \right] + \frac{\left( \frac{R_a \cos \alpha}{5830} \right)^{1/3} - 1}{2} \left[ \frac{\left( \frac{R_a \cos \alpha}{5830} \right)^{1/3} - 1}{2} \right] \qquad (12k)$$
where $\alpha$ 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):

$$ \left\{ \begin{aligned} T_i^{j+1} &= C_{1T_i^j} T_i^j + C_{2T_i^j} (T_{i+1}^j + T_{i-1}^j) \\ H_i^{j+1} &= C_{1H_i^j} H_i^j + C_{2H_i^j} (H_{i+1}^j + H_{i-1}^j) \end{aligned} \right. \qquad (13) $$

$$ C_{1T_i^j} = 1 - \frac{2D_{\text{Tcas}}\Delta t}{(\Delta x)^2} \Big|_i^j \qquad (14a) $$

$$ C_{2T_i^j} = \frac{D_{\text{Tcas}}\Delta t}{(\Delta x)^2} \Big|_i^j \qquad (14b) $$

$$ C_{1H_i^j} = 1 - \frac{2D_{\text{Hcas}}\Delta t}{(\Delta x)^2} \Big|_i^j \qquad (14c) $$

$$ C_{2H_i^j} = \frac{D_{\text{Hcas}}\Delta t}{(\Delta x)^2} \Big|_i^j \qquad (14d) $$

Equations (6c) and (6d) give the following:

$$ T_0^j = T_f^j(z = z_{\text{max}}) \qquad (15a) $$

$$ H_0^j = X_{\text{eq}}(T_f^j(z = z_{\text{max}}); RH_f^j(z = z_{\text{max}})) \qquad (15b) $$

Equations (6e) and (6f) give the following:

$$ T_0^j = T_{aE}^j \qquad (15c) $$

$$ H_0^j = X_{\text{eq}}(T_{aE}^j; RH_{aE}^j) \qquad (15d) $$

Equations (6g) and (6h) give respectively the following:

$$ T_{N+1}^j = \left(1 - \frac{h_c \Delta x}{\lambda_{\text{cas}}}\right) \Big|_{N+1}^j T_N^j + \frac{h_c \Delta x}{\lambda_{\text{cas}}} \Big|_{N+1}^j T_a^j \qquad (16a) $$

$$ H_{N+1}^j = \left(1 - \frac{h_m \Delta x}{D_{\text{Hcas}}}\right) \Big|_{N+1}^j H_N^j + \frac{h_m \Delta x}{D_{\text{Hcas}}} \Big|_{N+1}^j X_{\text{eq}}(T_a^j; RH_a^j) \qquad (16b) $$

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:

$$ \frac{\dot{m}C_{\text{pf}}}{l_c} \frac{\partial T_f}{\partial z} = h_{\text{cpf}}(T_{\text{pc}} - T_f) + h_{\text{cfg}}(T_{\text{gc}} - T_f) \qquad (17a) $$

With:

$$ T_f(z = 0) = T_{aE} \qquad (17b) $$
Along with the solar collector, the heater of the fresh air coming from ambient is given by the analytical solution:

$$T_f(z) = \frac{h_{\text{cfp}} T_{\text{pc}} + h_{\text{cgf}} T_{\text{gc}}}{h_{\text{cfp}} + h_{\text{cgf}}} + \left( T_{\text{aE}} - \frac{h_{\text{cfp}} T_{\text{pc}} + h_{\text{cgf}} T_{\text{gc}}}{h_{\text{cfp}} + h_{\text{cgf}}} \right) \exp \left( - \frac{(h_{\text{cfp}} + h_{\text{cgf}})z}{\rho_f e_{\text{sc}} V_{\text{aE}} C_{\text{pf}}} \right) \quad (18a)$$

With:

$$\dot{m} = \rho_f e_{sc} l_c V_{aE} \quad (18b)$$

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.

<table><caption>Table 3. Thermal properties of the components of solar dryer [43,47].</caption><thead><tr><th>Components</th><th>Specific Heat Capacity (J/(kg·K))</th><th>Thermal Conductivity (W/(m·K))</th><th>Absorptivity Coefficient (-)</th><th>Transmitivity Coefficient (-)</th></tr></thead><tbody><tr><td colspan="5" style="text-align: center;">For the drying chamber</td></tr><tr><td>Roof</td><td>900</td><td>190</td><td>0.91</td><td>0</td></tr><tr><td>Wall</td><td>1600</td><td>0.13</td><td>0.4</td><td>0</td></tr><tr><td>Cassava roots chips</td><td>See Equation (2)</td><td>See Equation (3)</td><td>(-)</td><td>(-)</td></tr><tr><td>Air drying</td><td>See Equations (A2) and (A3)</td><td>See Equation (A1)</td><td>(-)</td><td>(-)</td></tr><tr><td colspan="5" style="text-align: center;">For the solar collector</td></tr><tr><td>Glass cover</td><td>800</td><td>0.1311</td><td>0.1</td><td>0.89</td></tr><tr><td>Absorber</td><td>900</td><td>190</td><td>0.91</td><td>0</td></tr><tr><td>Wooden insulation</td><td>1236</td><td>0.14</td><td>(-)</td><td>(-)</td></tr><tr><td>Airflow</td><td>See Equations (A2) and (A3)</td><td>See Equation (A1)</td><td>0.1</td><td>0.9</td></tr></tbody></table>

## 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.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Validation of the model using experimental data given by Tieu et al. [36] in Yamoussoukro, October 2018, design 1.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Influence of the design on the drying kinetic in Yamoussoukro, October 2018, designs 1 and 2.</figcaption></figure>

## 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/m² 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/m²s.
The global heat transfer coefficient below the solar collector $U_b$ varied around 5 W/m²s while the radiative heat transfer coefficients varied around 0.5 W/m²s. 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/m²s. 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 ($h_m$) during the sunny days as shown in fig Figure 8c. However, the heat transfer coefficient $h_m$ varies between $2.5 \times 10^{-7}$ and $2 \times 10^{-6}$ m/s.

<figure><img src="image_10.png" /><figcaption>**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.</figcaption></figure>

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 core 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 \times 10^{-4}$ kg/(kg·s).
<figure><img src="image_11.png" /><figcaption>**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).</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**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.</figcaption></figure>

**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**

<table><tr><td>A</td><td>area (m<sup>2</sup>)</td></tr><tr><td>C, G, C<sub>o</sub>, G<sub>o</sub></td><td>parameters of desorption isotherms (-)</td></tr><tr><td>Cp</td><td>specific heat capacity (J/kg·K)</td></tr><tr><td>D<sub>a</sub></td><td>mass diffusion coefficient of the vapor in drying air (m<sup>2</sup>/s)</td></tr><tr><td>D<sub>H</sub></td><td>mass diffusion coefficient (m<sup>2</sup>/s)</td></tr><tr><td>d<sub>H</sub></td><td>hydraulic diameter (m)</td></tr><tr><td>D<sub>o</sub></td><td>parameter of the mass diffusion coefficient (m<sup>2</sup>/s)</td></tr><tr><td>D<sub>T</sub></td><td>thermal diffusion coefficient (m<sup>2</sup>/s)</td></tr><tr><td>E<sub>a</sub></td><td>energy of activation (J/mol)</td></tr><tr><td>e<sub>b</sub></td><td>the thickness of the wooden insulation in the solar collection (m)</td></tr><tr><td>e<sub>pc</sub></td><td>the thickness of the absorber in the solar collector (m)</td></tr><tr><td>e<sub>sc</sub></td><td>the thickness of the airflow in the solar collector (m)</td></tr><tr><td>F<sub>a-b</sub></td><td>geometrical factor between the areas (a) and (b) (-)</td></tr><tr><td>g</td><td>the intensity of the presenter (9.81 m<sup>2</sup>/s)</td></tr><tr><td>Gr</td><td>Grashof number (-)</td></tr><tr><td>H</td><td>moisture content of the cassava (kg/kg, db)</td></tr><tr><td>h<sub>c</sub></td><td>convective heat transfer coefficient (W/m<sup>2</sup>·K)</td></tr><tr><td>h<sub>m</sub></td><td>convective mass transfer coefficient (m/s)</td></tr><tr><td>h<sub>r</sub></td><td>radiative heat transfer coefficient (W/m<sup>2</sup>·K)</td></tr><tr><td>I<sub>G</sub></td><td>global solar irradiation received on the solar dryer (W/m<sup>2</sup>)</td></tr><tr><td>l<sub>c</sub></td><td>width of the airflow in the solar collector (m)</td></tr><tr><td>m</td><td>mass (kg)</td></tr><tr><td>mo</td><td>dry mass of cassava (kg)</td></tr><tr><td>ṁ</td><td>mass rate of the drying air (kg/s)</td></tr><tr><td>Nu</td><td>Nusselt number (-)</td></tr><tr><td>Sh</td><td>Sherwood number (-)</td></tr><tr><td>P</td><td>pressure (Pa)</td></tr><tr><td>Pr</td><td>Prandtl number (-)</td></tr><tr><td>R</td><td>perfect gases universal constant (8.314 J/K·mol)</td></tr><tr><td>Ra</td><td>Rayleigh number (-)</td></tr><tr><td>RH</td><td>relative humidity of drying air (%/100)</td></tr><tr><td>Sur</td><td>area of the walls of the drying chamber (m<sup>2</sup>)</td></tr><tr><td>T</td><td>temperature (K)</td></tr><tr><td>t</td><td>drying time (s)</td></tr></table>
<table><tr><td>U<sub>b</sub></td><td>global heat transfer coefficient in the bottom of the solar collector (W/m<sup>2</sup>.K)</td></tr><tr><td>Vol</td><td>volume of the drying chamber (m<sup>3</sup>)</td></tr><tr><td>x</td><td>axis coordinates in the thickness direction of the sample (m)</td></tr><tr><td>X</td><td>moisture content of cassava at equilibrium state (kg/kg, db)</td></tr><tr><td>X<sub>mo</sub></td><td>moisture content of the first layer when T = 0 K</td></tr><tr><td>X<sub>mo1</sub></td><td>moisture of bound water extracted per kelvin of temperature (kg/kg.K)</td></tr><tr><td>Y<sub>s</sub></td><td>humidity of the drying air (kg/kg)</td></tr><tr><td>z</td><td>axis coordinates of the solar collector (m)</td></tr><tr><td>zo</td><td>altitude (m)</td></tr><tr><td>ΔH<sub>C</sub>, ΔH<sub>G</sub></td><td>the energy needed to dry a mole of bound water (J/mol)</td></tr></table>

**Greek Symbols**

<table><tr><td>λ</td><td>thermal conductivity (W/(m. K))</td></tr><tr><td>ρ</td><td>density (kg/m<sup>3</sup>)</td></tr><tr><td>μ<sub>a</sub></td><td>dynamic viscosity of drying air (Pa·s)</td></tr><tr><td>ν<sub>a</sub></td><td>kinematic viscosity of the drying air (m<sup>2</sup>/s)</td></tr><tr><td>σ</td><td>Stephan Boltzmann coefficient (5.67 × 10<sup>-8</sup>)</td></tr><tr><td>α</td><td>the angle between the solar collector and the horizontal plane (°)</td></tr><tr><td>α<sub>i</sub></td><td>solar absorptivity coefficient of component i (-)</td></tr><tr><td>τ<sub>i</sub></td><td>solar transmissivity coefficient of component i (-)</td></tr><tr><td>ρ<sub>p</sub></td><td>solar reflexivity coefficient of the absorber (-)</td></tr><tr><td>ε<sub>i</sub></td><td>solar emissivity coefficient of component i (-)</td></tr></table>

**Subscripts**

<table><tr><td>a</td><td>drying air in the drying chamber or solar collector</td></tr><tr><td>aE, E, ext</td><td>ambient (exterior space of the solar drying)</td></tr><tr><td>air</td><td>air</td></tr><tr><td>atm</td><td>atmosphere</td></tr><tr><td>cas</td><td>cassava</td></tr><tr><td>ep</td><td>surface</td></tr><tr><td>eq</td><td>equilibrium</td></tr><tr><td>f</td><td>airflow in the solar collector</td></tr><tr><td>fs</td><td>end position of the solar collector when z = z<sub>max</sub> (outlet point)</td></tr><tr><td>g</td><td>gas (air)</td></tr><tr><td>gc, g</td><td>glass cover</td></tr><tr><td>i</td><td>the interior of the dryer</td></tr><tr><td>La</td><td>lateral part</td></tr><tr><td>m</td><td>the first layer of bound water</td></tr><tr><td>pc</td><td>absorber in the solar collector</td></tr><tr><td>r</td><td>roof</td></tr><tr><td>sk</td><td>sky</td></tr><tr><td>v</td><td>vapor</td></tr><tr><td>vsat</td><td>vapor in saturation state</td></tr><tr><td>w</td><td>wall</td></tr></table>

**Index**

<table><tr><td>i,j</td><td>order of discretization scheme (-)</td></tr></table>

**Appendix A. Thermophysical Parameters of the Drying Air**

Appendices See [2,34,49].

Thermal conductivity (W/(m. K)), (T<sub>a</sub> in °C):

$$ \lambda_a = 7.57 \times 10^{-5} T_a + 0.0242 \quad (A1) $$

Heat capacity of dry air (J/(kg. K)), (T<sub>a</sub> in °C):

$$ C_{pa} = 1835 - 0.734 T_a \quad (A2) $$
Heat capacity of vapor air (J/(kg·K)), (T_a in °C):

$$C_{pv} = 999.2 + 0.1434T_a + 1.101 \times 10^{-4}T_a^2 - 6.7581 \times 10^{-8}T_a^3 \quad (A3)$$

Dynamic viscosity (Pa·s), (T_a in °C):

$$\mu_a = 1.718 \times 10^{-5} + 4.620 \times 10^{-8}T_a \quad (A4)$$

Mass diffusion coefficient of vapor in air:

$$D_a(\text{m}^2/\text{s}) = \frac{2.25 \times 10^{-5}}{P_{\text{atm}}} \left(\frac{T_a}{273}\right)^{1.81} \quad (A5a)$$

$$P_{\text{atm}} = 101,325 - 12zo + 5.2 \times 10^{-4}zo^2 \quad (A5b)$$

where zo in m is the altitude, T_a is given in K, and P_atm in Pa.

The density of the air vapor:

$$\rho_g(\text{kg/m}^3) = \frac{353}{T_a} \quad (A6)$$

where T_a is given in K.

## Appendix B. Non-Dimensional Numbers

Prandtl number (-) (T_a in °C):

$$Pr = -2.54 \times 10^{-4}T_a + 0.7147 \quad (A7)$$

Grashof number (-) (T_a in K):

$$Gr = \frac{\beta g d_h^3 \rho_a^2 \Delta T}{\mu_{\text{air}}^2} ; \quad \beta = \frac{1}{T_a} \quad (A8)$$

## 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.

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. [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. [CrossRef]

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.

5. Ndukwu, M.C.; Onyenwigwe, D. Development of a motorized parboiled cassava tuber shredding machine. *Niger. J. Technol.* **2013**, *32*, 1–6.

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.

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. [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. [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.

10. Sharma, A.; Chen, C.R.; Vu Lan, N. Solar-energy drying systems: A review. *Renew. Sustain. Energy Rev.* **2009**, *13*, 1185–1210. [CrossRef]

11. Ndukwu, M.C.; Bennamoun, L. Potential of integrating Na₂SO₄·10H₂O pellets in solar drying system. *Dry. Technol.* **2018**, *36*, 1017–1030. [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. [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.

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. [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*. [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?excludensity=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. [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. [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**, *77*, 100–109. [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. [CrossRef]

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. [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. [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. [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. [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. [CrossRef]

27. Chaudhari, A.D.; Salve, S.P. A review of solar dryer technologies. *Int. J. Res. Advent Technol.* **2014**, *2*, 218–232.

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. [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. [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. [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. [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. [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. [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. [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. [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.

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.

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. [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.
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. [CrossRef]

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.

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.

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. [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.

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. [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. [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. [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.