Theoretical Model for Predicting Moisture Ratio during Drying of Spherical Particles in a Rotary Dryer

Amathematical model was developed for predicting the drying kinetics of spherical particles in a rotary dryer. Drying experiments were carried out by drying fermented ground cassava particles in a bench scale rotary dryer at inlet air temperatures of 115–230C, air velocities of 0.83m/s–1.55m/s, feedmass of 50–500 g, drum drive speed of 8 rpm, and feed drive speed of 100 rpm to validate the model. The data obtained from the experiments were used to calculate the experimental moisture ratio which compared well with the theoretical moisture ratio calculated from the newly developed Abowei-Ademiluyi model. The comparisons and correlations of the results indicate that validation and performance of the established model are rather reasonable.


Introduction
Rotary drying is a very complicated process that can be applied not only to thermal drying but also movement of particles within the dryer.Several authors have carried out investigations on the steady state modeling of the rotary drying process.Static models are in general differential equations, and they are suitable for investigation of static distributions.Myklestad [1] was the first to obtain an expression to predict product moisture content throughout a rotary dryer based on drying air temperature, initial moisture content, and product feed rate.Thin layer drying equations contribute to the understanding of the heat and mass transfer phenomena in agricultural products and computer simulations for designing new and improving existing commercial drying processes [2].They are used to estimate drying times of several products and also to generalize drying curves.In thin layer drying model, the rate of change in material moisture content in the falling rate drying period is proportional to the instantaneous difference between material moisture content and the expected material moisture content when it comes into equilibrium with the drying air [3].
Many authors have developed semiempirical models based on the diffusion theory to predict the drying kinetics of moist substances in thin layer as shown in Table 1 (where MR is the moisture ratio).The constants , , , ,   ,  1 , and  in eight models by most authors have been found to be functions of inlet air temperature, inlet air velocity, humidity, and so forth, the mass of feed was not accounted for by all the authors and in the drying of substances with high moisture content like fermented ground cassava, dairy products, and some pharmaceutical product in rotary dryer, and the mass of feed should be accounted for in the thin layer drying equation.It was observed that although several models have been proposed, there is not a general theory to describe the mechanism of rotary drying, and it seems that specific models for an equipment and material are more useful than general models [4].
Therefore the objective of this study is to develop a theoretical model for predicting the drying kinetics of spherical particles in a rotary dryer accounting for the quantity of materials to be dried in the model.

Theoretical Development of Thin Layer Drying Equation.
Fick's diffusion equation ( 1) has been accepted for describing  the drying characteristics of biological and chemical products in the falling rate period [16] as follows: where  is the diffusion coefficient,  is moisture content (dry basis) at any time , and  is drying time.The equation of diffusion for a spherical particle at constant diffusivity and radial (as shown in Figure 1) flux takes the following form: In order to solve (2) the following assumptions were adopted: (1) moisture movement is only diffusion and unidirectional; (2) diffusion coefficient  is independent of moisture concentration; (3) drying process is isothermal, that is, adiabatic dryer; (4) material to be dried is spherical in shape; (5) shrinkage is neglected.
Using − 2 as a separation constant we obtain from (2) Integrating (3) using separation of variables gives Equation ( 4) is of the form Equation ( 6) is a Bessel equation of order zero, the solution of which is [17] But Combining ( 5) and (9) gives ( ⋅ ) = ()(), so that and applying the boundary conditions in (11).The solution to (10) in the case of a sphere is expressed as where  is the radius of sphere, MR is moisture ratio,   is initial moisture content (% db),   is equilibrium moisture content (% db),   is moisture content at time  (% db), and  is drying time (hr).
From the work of Abowei [18], the mass of hydrocarbon  was accounted for when modeling one-dimensional diffusion of oil spill in water and obtain a general solution (see (13)) to predict the diffusion of known quantity of crude oil in water.This equation is analogous to the diffusion equation ( 12) describing diffusion of moisture in porous spherical particles as follows: where   is the quantity of oil spilled and   is the concentration of oil spilled at any time. is the area where oil is spill;   is the diffusion coefficient, and  is the time.
Comparing ( 12) with ( 13), the term 13), and hence ( 12) can be rewritten as where MR is the moisture ratio. is surface area available for moisture transfer which for the rotary dryer is ( 2 + ),   in equation ( 14) is now the mass of fermented ground cassava which is analogous to the   in equation ( 13),  is the average density of sample to be dried, and the density is added   equation ( 14) to make the equation dimensionless, since moisture ratio MR is dimensionless, so that ( 14) becomes where  is the average radius of particle to be dried  is the radius of rotary dryer drum,  is the length of the rotary dryer, and where   is the activation energy.Equation ( 15) is the new theoretical Abowei-Ademiluyi model for predicting the drying of any spherical particles in a rotary dryer.Equation ( 15) was simulated to obtain the theoretically determined moisture ratio.

Dimensional Analysis Approach.
In order to remove moisture from a moist material in a rotary dryer, the moisture ratio (MR) can be taken to be a function of the change in temperature Δ, the quantity of fermented ground cassava to be dried   , the latent heat , diameter of particle  to be dried, inlet air velocity V, and drum speed , so that mathematically the moisture ratio MR is dimensionless as where  is a correction factor.
Applying dimensional analysis we have Applying the Buckingham  method gives ∑  : 0 = .

Experimental Work.
In order to validate the model, fermented ground cassava particles was dried in a bench scale rotary dryer (Figure 2).The developed theoretical model (15) was simulated with Microsoft Excel 2007 using the following data: (i) average density product (kg/m 3 ) = 400, [  Development Project farm (ADP) at Rumuokoro, Port Harcourt.The choice of this cassava cultivar TMS 30572 was based on its preference by farmers, because of its high yield and suitability for gari processing [20].The cassava cultivar was peeled, washed, grated, and packed in sack for pressing.The dewatered mash was allowed to ferment naturally for 72 hrs; sieved with a mesh of 3.5 mm, and then dried in a bench rotary dryer (Figure 2).

Experimental Procedure.
At the beginning of each experiment, the dryer was allowed to reach steady state at the desired airflow rate, inlet air temperature, feed drive speed, and drum drive speed.When steady state condition had been attained, the fermented ground cassava mash of known moisture content was introduced into the dryer feed hopper.The drying conditions used in the experiments are inlet air temperatures of 115 ∘ C, 140 ∘ C, 190 ∘ C, and 230 ∘ C, air velocities of 0.83, 1.02, 1.397, and 1.55 m/s, mass of feed of 50 g, 100 g, 200 g, and 500 g, feed drive speeds 100 rpm, and drum drive speeds of 8 rpm.The decrease in mass of fermented mash was monitored with time per pass.The initial moisture content of samples was determined separately before start of experiment.The weight loss during drying was used to calculate the moisture content.The drying data obtained were used to calculate the experimental moisture ratio (MR) to predict the kinetics of drying fermented ground cassava.

Simulated Theoretical Results
. The theoretical moisture ratio results are presented in Figures 3, 5, and 7.The theoretical moisture ratio decreases with drying time as inlet air temperature, inlet air velocity, and mass of feed increases.A similar profile is also exhibited in Figures 4, 6, and 8 for the experimental moisture ratio.The theoretical moisture ratio plots show a typical drying curve generally obtained during drying of moist materials [3,21].
It can be observed from the theoretical moisture ratio plots that the Abowei-Ademiluyi model does not give values for moisture ratio at  = 0, and this will not be a problem, since the initial moisture content from which the moisture ratio at  = 0 (i.e.,   ) was calculated is always known at start of drying.Hence the Abowei-Ademiluyi model can be used to predict the drying kinetics of spherical particles at any Modelling and Simulation in Engineering

Comparison of Theoretical and Experimented
Results.The simulated theoretical result compared favorably with those of the experimental results.The similarity is shown from the high value ( close to 1) obtained for the coefficient of multiple determinations  2 at different inlet air temperature and inlet air velocity as shown in Figures 9 and 10.However, better fit could be obtained if the average density particle is correctly chosen.The theoretical (Abowei-Ademiluyi model) moisture ratio also compared well with experimentally moisture ratio at different mass of feed as shown in Figure 11.

Conclusion
The new theoretical, Abowei-Ademiluyi model has been developed for predicting drying kinetics of spherical particles at any known particle diameter, rotary drum diameter, and dryer length.The new model also account for the mass of feed.Model validation was carried out by drying fermented ground cassava particles in a bench scale rotary dryer at inlet air temperatures of 115-230 ∘ C, air velocities of 0.83 m/s-1.55m/s, feed mass of 50-500 g, drum drive speed 8 rpm,   and feed drive speed of 100 rpm.The theoretical moisture ratio calculated from the model compared favorably with experimental moisture ratio.

Figure 1 :
Figure 1: Showing hypothetical profile of moisture diffusion from a spherical particle in a rotary dryer.

Figure 3 :
Figure 3: Variation of moisture ratio with time at different air inlet temperatures using the theoretical model.

Figure 7 :
Figure 7: Variation of moisture ratio with time at different mass of feed using the theoretical model.

Figure 8 :
Figure 8: Experimental moisture ratio at different mass of feed.

Figure 11 :
Figure 11: Theoretical (Abowei-Ademiluyi model) and experimental moisture ratios at different mass of feed.

Table 1 :
Mathematical models given by various authors for the drying curves.