1. Introduction
The removal of moisture can be due to simultaneous heat and mass transfer [1] . Purposely it is carried out to reduce water to the level at which microbial spoilage and deterioration reactions are greatly minimized [2] . The dried product achieved will minimize transportation, storage, packaging cost and time. Though there are many ways in which drying can be achieved, but the choice of method depend on the material and the sanitary level required [3] . The use of drying of foodstuffs is very common to improve food stability and minimize chemical and physical changes during storage. Hence, drying is the most extensively used technique for food preservation [4] . Moreover, there are various types of drying techniques that can be applied to reduce the water and which attains the purpose of food preservation as has been shown in Singh et al. [4] and Inyang, et al. [5] , these include, air drying, hot air drying, sun drying, solar drying, osmotic drying, infra-red drying, oven drying, spray drying, drum drying, microwave drying, freeze drying, vacuum drying, foam mat drying, impingement drying, acoustic drying, explosion puffing drying, hybrid dryers.
Kinetic modeling of process parameters is very useful in food processes. The processes involved are mainly (bio) chemical and physical reactions. These changes proceed at a certain rate and with certain kinetics. Kinetic modeling enables us to describe these changes and their rates quantitatively. In kinetic modeling we have a powerful tool that can help to unravel basic reaction mechanisms. The understanding of the basic mechanisms is vital for quality modeling and quality control [6] .
To understand the progress of reactions, knowledge of thermodynamics and kinetics is required. The rate with which a reaction proceeds is the resultant of the driving force and the resistance against change. There is thus an intimate link between thermodynamics and kinetics [6] .
Hence, studying drying kinetics will aid in identifying appropriate drying methods and to control the processes of drying. It is also important for engineering and process optimization. It is sometime expensive to conduct a full scale experiment to determine the suitable conditions for drying. Thus, drying kinetics is used to express the moisture removal process and its relation to the process variables and hence, a good understanding of the drying rate is important to develop a drying model [7] .
Despite the essence of kinetics, modeling of particulate or thin-layer drying of materials is necessary to understand the fundamental transport mechanism and a prerequisite to successfully simulate or scale up the whole process for optimization or control of the operating conditions. Simple models with a reasonable physical meaning are effective for engineering purposes. Mathematical modeling of dehydration process is an inevitable part of design, development and optimization of a dryer [8] . It mainly involves elaborative study of drying kinetics, which describes the mechanisms and the influence that certain process variables exert on moisture transfer [9] . In other words, it can be used to study the drying variables, evaluate the drying kinetics and to optimize the drying parameters and the conditions [10] [11] . A proper drier design requires knowledge on the characteristics of the material to be dried and the drying kinetics [12] . It is noted that higher temperature implies larger driving force for heat transfer [13] [14] . It also accelerates the drying process, as the temperature provides a larger water vapor pressure deficit [15] .
2. Mechanism of Drying
The main mechanisms of drying are surface diffusion on the pore surfaces, liquid or vapor diffusion due to moisture concentration differences and capillary action in granular and porous foods due to surface forces [16] . Generally, hygroscopic products dry at constant rate and subsequent falling rate periods and drying stops when equilibrium is established. The constant rate period of drying, the physical form of the product and external conditions such as temperature, drying air velocity, direction of air flow and relative humidity have a great influence on the surface of the product being dried in surface diffusion [17] [18] . When the surface film of the solids or particles appears to be dried and the moisture content has been reduced to its critical moisture content then the first falling rate period begins. Unlike the constant rate periods, the falling rate period is controlled by liquid diffusion as a result of moisture concentration differences and the internal conditions of the product. The internal conditions such as moisture content, the temperature and the structure of the product play an important role in the falling rate periods. This phenomena is then replaced by the second falling rate periods of drying namely vapor diffusion due to moisture concentration difference and also internal conditions of the products [19] [20] . It has been accepted that the drying phenomenon of biological products during falling rate period is controlled by the mechanism of liquid and/or vapor diffusion [19] [21] . However, only capillary and liquid diffusion theories are, generally, applicable to the drying of food materials [22] . Diffusion is the dominant physical mechanism governing moisture movement in the material [23] [24] [25] , which is dependent on the moisture content of the samples [26] . The moisture contents were expressed on dry basis, which is more convenient for modeling [27] [28] .
Figure 1 describes the drying rate and temperature as a function of time. This rate curve can also be used in identifying the dominant mechanism of a product during drying. In the initial drying period, the equilibrium air temperature is usually greater than the temperature of the product [30] . Therefore, the drying rate between A and B increases with an increase in temperature of the product until the surface temperature attains equilibrium (Corresponding to line B to C). Under constant conditions, the drying process of agricultural and biological products has been described as a number of steps consisting also of an initial constant rate period (B to C) during which drying occurs as if pure water is being evaporated, and one or several falling rate periods where the moisture movement is controlled by combined external-internal resistances or by either external or internal resistance to heat and mass transfer [31] . Mostly, many fruits and vegetables dry during the falling rate periods because the drying process is controlled by a diffusion mechanism. Drying usually stops when steady state equilibrium is reached [16] . During the constant rate period the physical form of
![]()
Figure 1. A typical drying curve of agricultural products showing constant rate and falling rate periods [29] .
the product is affected, especially the surface of the product. This period is largely controlled by capillary and gravity forces. The conditions of the drying process, like the temperature, drying air velocity, and relative humidity, also affect the product during this stage. The first falling rate period (C to D) begins when the surface film of the product appears to be dry, and the moisture content has decreased to its critical moisture content. As drying continues, the material will then experience a change from the first falling rate period to a phenomenon known as the second falling rate period (D to E) [29] .
In literatures, several models like the theoretical, empirical and semi empirical are mentioned for the analysis of the drying of hygroscopic products [21] [32] . But the most widely applied categories of thin-layer models are the semi theoretical and empirical models [29] [33] [34] [35] .
3. Thin Layer Drying Process
Thin layer refers to a layer of sufficiently small product thickness in order that one could consider that air characteristics everywhere in the layer are identically uniform without variance. Thin layer drying process also refers to the drying of individual particles or grains of material which are fully exposed to the drying air .The process is often divided into two periods of drying which are the constant drying rate period and the falling drying rate period [36] .
In the constant rate drying period, the material contain so much water that liquid surfaces exist and will dry in a manner comparable to an open faced body of water. The water and its environment, not solid will determine the rate of drying. Wet sand, soil, pigment and washed seeds are examples of materials that initially dry at constant rate. Radiation, conduction or convection provides the energy needed for drying.
In actual practice, drying of agricultural produce takes place in the falling rate period. The falling rate period is bounded by equilibrium moisture content of an equilibrium moisture curve between zero and nearly 100% relative humidity. Drying in the falling rate period involves two processes namely; the movement of moisture within the material to the surface and removal of moisture from the surface [37] . Furthermore, thin layer drying could be viewed as the drying of one layer of sample particles or slices in abundant supply of drying air with appropriate physical characteristics.
Thin layer drying is also described as the process of removal of moisture from a porous media by evaporation in which excess drying air is passed through a thin layer of a material until the equilibrium moisture content is reached [38] [39] .
Thin-layer drying equations fall into three categories namely, theoretical, semi-theoretical, and empirical models.
3.1. Theoretical Method
Theoretical models only consider the internal resistance in transfer of moisture while the other two categories are considered to study the external resistance to moisture transfer between air and product [40] [41] . The theoretical models are derived from Fick’s second law of diffusion [21] . While semi-theoretical models are generally derived from Fick’s second law and its modifications and also, from Newton’s law of cooling [16] . Fick’s second law of diffusion is the most commonly used theoretical model. Theoretical models are inadequate and tend to generate erroneous results and are complex for practical applications. The theoretical models make too many assumptions leading to a considerable number of errors [40] [42] , thus limiting their utilization in the design of dryers [29] . A theoretical equation gives a better understanding of the transport processes but an empirical equation gives a better fit to the experimental data without any understanding of the transport processes involved [19] .
3.2. Empirical Method
The empirical method is a method based on experimental data and dimensionless analysis. The empirical models of drying show a direct relationship between average moisture content and drying time [43] . However, this method omits the fundamentals of the drying process and its parameters have no physical meaning, therefore it does not provide an accurate view of the important processes that occur during the phenomenon while describing the drying curves for certain experimental conditions [44] . Empirical models fails to take into account the fundamentals of drying process and is able to explain only the drying curve for drying condition but not the processes that occurs during drying [45] . Empirical models help to understand the trend of experimental/process variables both dependent and independent. The main challenges faced by the empirical models are that they depend largely on experimental data and provide limited information about the heat and mass transfer during the drying process [16] .
3.3. Semi-Empirical Method
Semi theoretical models have been developed to ease the use and to fit the drying data of the food material to be dried [46] . Much emphasis has been given to developing semi theoretical models for achieving harmony between theory and ease of use. Such models are generally based on Newton’s Law of Cooling applied to mass transfer. When applying this law it is assumed that conditions are isothermal and resistance to moisture transfer is confined only to the surface of the product [47] . Semi theoretical models are the simplified general series solutions of the Fick’s second law. These models, however works only when temperature, relative humidity, air velocity and moisture are in the range at which the models are developed. Time requirement is less and shapes of the material dried are not considered in semi theoretical models [48] .
Among semi-theoretical models are the model of two terms, the Henderson and Pabis, Lewis, Page and Modified Page [21] .
4. Drying Constant
In the thin layer drying concept, the drying constant is the combination of drying transport properties such as moisture diffusivity, thermal conductivity, density, specific heat, interface heat, and mass coefficients [49] . Thus, knowledge of transport and material properties is necessary to apply any transport equation [50] .
5. Derivation of Thin layer Models
Research on the study of the kinetics of thin layer drying is performed with various agricultural products such as seeds, grains, fruits and some plant species with economic importance [44] . The Henderson and Pabis model, firstly used for model drying of corn is the first term of general series solution to Fick’s law. However, due to inaccuracy and high degree of temperature difference between kernel and air, the model could not be fitted during the first 1 or 2 hour of drying corn [51] . The first two terms of general series solution to Fick’s second law is the two term model. The two term model however, does not consider the shape of the product dried. The model assumes constant diffusivity and requires the temperature of the product to be constant [40] [52] .
Lewis model is a special case of Henderson and Pabis model. According to Bruce [42] the model is inaccurate; it overestimates the first period and underestimates the last period of drying.
Drying kinetics models does not take into account the effects of interactions by parameters other than the time of drying. Models that incorporate a large number of variables still do not exist but due to the complex non-linear relationship between the kinetics of drying and variables related, the development of such models is not feasible [32] [53] . The concept of thin-layer drying models for characterizing the drying behaviour was suggested, initially, by Lewis [54] who derived the semi-theoretical model for porous hygroscopic materials, which is analogous with Newton’s law of cooling. The following model was developed
(1)
where MR is moisture ratio, k is drying constant (
), t is drying time,
,
,
are moisture content at any time, equilibrium and initial, respectively. Page [55] modified the Lewis model by adding a dimensionless empirical constant (n) and used it for study the drying behavior of shelled corns.
(2)
For study the drying kinetics of soybeans, Overhults et al. [56] modified the Page model and obtained the following equation (this model is known as Modified Page-I Model)
(3)
In addition, White et al. (1978) make a little change in (3) to describe the drying kinetics of soybeans (this model is known as Modified Page-II Model)
(4)
Further, for the drying of sweet potato, Diamante and Munro [57] modified the Page model and proposed the following equation (this model is known as Modified Page Equation (2) Model)
(5)
where l is an empirical constant (dimensionless).
The drying process of apple and potato is considered as the external resistance to moisture transfer between air and product, therefore, the drying behavior of apple and potato can be characterized by semi-theoretical thin-layer drying models. The above discussed models have been tested for experimental data for apple and potato opted from literature. These models are inadequate to describe the drying behavior of apple and potato in last drying hours. Therefore, in this context we have proposed the following model by adding a linear term in the Lewis model for drying kinetics of apple and potato
(6)
where k is the drying constant (
) and a is the adjusting model constant (dimensionless) introduced for obtaining the best fit of the model to experimental data. The drying constant k depends on drying conditions, temperature, air velocity, humidity and drying methods. The moisture ratio (
) can be calculated
by (7) instead of
due to the small value of
, as compare to
, and
[58] .
(7)
Therefore, (6) reduced to
(8)
The theoretical models are derived from Fick’s second law of diffusion [21] while semi-theoretical models are generally derived from Fick’s second law by its modifications and from Newton’s law of cooling [16] . We are more concern for semi-theoretical thin-layer drying model due to better prediction of results as compare to other models [4] .
These models are just modification of one model or the other, taking the basis from Newton (Lewis) or Page model as could be seen in the literature [29] . It is to overcome the shortcomings of one another. The models are found to successfully describe the drying behavior (kinetics) of different types of fruits and vegetables. Each of the equations has its own dimensionless (model) constants (a, b, c, d, l,
) and drying constants (k, g, h, K1, K2, K0).
Table 1 below shows the mathematical models applied to thin layer drying curves of various agricultural products.
6. Classification of Model
Drying processes are usually modeled using two main models, the distributed element model and the lumped element model [16] . These can be described individually as follows:
6.1. Distributed Element Model
This model or system is based on the interaction between time and one or more
![]()
Table 1. Mathematical models applied to thin layer drying curves of various agricultural products.
spatial variables for all of its dependent variables. The distributed element model considers the simultaneous mass and heat transfer for the drying processes. It is important to note that the pressure effect is negligible compared to the temperature and moisture effect as reported by Brooker et al. [87] .
6.2. Lumped Element Model
This model or system considers the effect of time alone on the dependent variables. The lumped element model does not consider the change in temperature of a product and assumes a uniform distribution of drying air temperature within the product. The model includes assumptions from the Luikov equations, that is, the pressure variable is negligible and the temperature is constant [29] [88] . Lumped parameter model can be used as temperature is uniform in thin layer drying [19] .
7. Factors That Affects Drying Kinetics
Drying of food products like fruits and vegetables are greatly affected by the following factors: drying time, drying temperature, relative air humidity, air flow rate (or velocity), surface area, size (material thickness), volume and local or partial pressure [89] .
Amongst these conditions, the most influential in terms of drying fruits and vegetables are material thickness (size) and drying temperature [90] [91] [92] . Thus, Cerquera et al., [93] included air speed, the final moisture content, contraction or shrinkage of the product and the features that are related to the preservation and quality of the food as most influential parameters. During the drying process the air velocity and relative humidity were found to be the least significant factors that affect the drying kinetics of fruits and vegetables, while temperature and thickness were reported to be the factors that most affect thin-layer drying kinetics of fruits and vegetables [94] . But Gacula and Singh [95] showed that the temperature has more influence on the drying constant than the thickness. However, the interactions between the temperature and thickness did not significantly affect the value of the drying constant. These findings were similar to those of Wang [96] , who worked with potatoes.
8. Goodness of Fit Statistics
Thin-layer drying models are evaluated and can be compared using statistical measures. The quality of the fitted models are evaluated with the following statistical measures: Correlation coefficient (r or R), coefficient of determination (r2 or R2), Reduced chi square (χ2), mean bias error (MBE), root mean square error (RMSE), Sum square error (SSE), mean relative error root square (RRMS), modelling efficiency (EF), mean percent error (MPE), Mean square error (MSE). The model with the highest value for R or R2 will be selected to describe the drying curves when Correlation coefficient or coefficient of determination is used as the statistical measures and the reverse is the case for other statistical measures. Hence, goodness of fit is said to have occurred in a model if R2 value is high and other criteria such as χ2, RMSE, MBE, MPE, SSE, RRMS, EF, MSE values are low [27] [97] [98] [99] [100] . Nonlinear least square method based on the Levenberg-Marquardt method can also be used to fit the experimental data to select equations. The values of coefficient of determination R2 can be used to test the linear relationship between experimental and model calculated values [101] , and is one of the primary criteria for the selection of best model.
However, the selection of the most appropriate model for describing the drying behaviour of fruits and vegetables does not depend on the number of constants. Rather, it depends on various statistical indicators. The statistical indicators that have often been used to successfully select the most appropriate drying models as reported in literature [102] [103] [104] [105] .
Where N = n is the total number of observations, MRexp,i and MRpre,i are the experimental and predicted moisture ratio at any observation i. High R2 value which closer to 1 represents the best fit of the model [106] as can be seen in Table 2.
Table 2 showing various statistical measures (or indicators) which could be used for evaluating the goodness of fit.
9. Optimization
Optimization is the action of making the best or most effective use of a situation or resource(s), thus with regards to some criteria. Hence, it is a process of finding the conditions that give the maximum or minimum of a function. It is an act of obtaining the best result under given circumstances [110] . Using optimization
![]()
Table 2. Various statistical measures.
could help in predicting the best models for the processes of drying agricultural products in thin layers. Nowadays, predictive and very accurate tools such as Artificial Neural Networks (ANN) and Genetic Algorithm (GA) enjoy a special status in predicting and optimizing drying processes of various products [108] .
Genetic algorithm is a complicated optimization technique that studies a set of various defined characteristics to reach an optimal value through replication of evolutionary biological processes based on mutation and integration similar to what happens in genetics. This method has been used successfully in investigating various problems [111] . In general, many parameters including population size, number of generations, number of data items, etc. influence the optimization method. The optimal values of these parameters can be obtained through the trial and error method [108] .
Recently, interest in using Artificial Neural networks (ANN) as a modeling tool in food technology is increasing because they have found extensive utilization in solving many complex real world problems. Artificial neural networks are control systems necessary to solve problems in which the analytical methods are difficult to apply and their results have to be in a specific interval, e.g., in real time. The main objective of ANN is the development of mathematical algorithms that will enable ANNs to learn by mimicking information processing and knowledge acquisition in the human brain [112] .
ANN models contain layers of simple computing nodes that operate as nonlinear summing devices. Number of level or layers: the distribution of neurons in the neural network is done forming levels or layers of a determined number of nodes each one. There are input, output and hidden layers [113] .
The direction of the information flow is through the layers. The connectivity between the nodes of a neural network is related with the way in which the exits of neurons are directed to become into entrances other neurons. The output signal of a node can be one of the entrances another element of process, or even be an entrance to it. When any output of the neurons is input of neurons of the same level or preceding levels, the network is described as feed-forward. On the contrary, if there is at least one connected exit as entrance of neurons of previous levels or of the same level, including themselves, the network is denominated of feedback. The feedback networks that have at least a closed loop of back propagation are called recurrent [114] .
Though, both Response Surface Methodology (RSM) and ANN models provided good quality predictions, yet the ANN shows a clear superiority over RSM for both data fitting and estimation capabilities. On the other hand, ANN has the disadvantage of requiring large amounts of training data in comparison with RSM. This problem is solved by using statistical experimental design, to reduce the number of experiments [115] .
The most important stages in a biological process are modeling and optimization to improve a system and increase the efficiency of the process without increasing the cost. The classical optimization method (single variable optimization) is not only time-consuming and tedious but also does not depict the complete effects of the parameters in the process and ignores the combined interactions between physicochemical parameters. This method can also lead to misinterpretation of results. In contrast, response surface methodology (RSM) is an empirical modeling system for developing, improving, and optimizing of complex processes [115] .
Although RSM has so many advantages, and has successfully been applied to study and optimize the enzymatic processes and enzyme production from microorganisms, it is hard to say that it is applicable to all optimization and modelling studies. The past decade has seen a host of data analysis tools based on biological phenomena develop into well-established modelling techniques, such as artificial intelligence and evolutionary computing [115] .
Artificial neural networks (ANNs) are now the most popular artificial learning tool in biotechnology, with a wide applications range included optimization of bioprocesses and enzyme production from microorganisms. Indeed an ANN is a massively interconnected network structure consisting of many simple processing elements capable of performing parallel computation for data processing. The fundamental processing element of ANNs (the artificial neuron) simulates the basic functions of biological neurons [115] .
RSM is a frequently used method for empirical modeling and prediction in the processing of biological media. The ANN has recently grown to be one of the most efficient methods for empirical modeling and prediction, especially for non-linear systems [116] .
The role model for ANN is the human mind. ANN is a collection of methodologies that intend to utilize the endurance for ambiguity and uncertainty to achieve complete information and provide low cost solutions. ANN is playing an increasingly significant role in many application areas of science, engineering and technology [116] [117] .
Neural networks have the potential for solving problems in which some inputs and corresponding output values are known, but the relationship between the inputs and outputs is not well understood or is difficult to translate into a mathematical function. These conditions are commonly found in tasks involving grading, sorting and identifying agricultural products [116] .
ANNs are useful tools for food safety and quality analyses, which include modeling of microbial growth, and from this predicting food safety; interpreting spectroscopic data, predicting physical, chemical, functional and sensory properties of various agricultural products during processing, storage and distribution. ANNs hold a great deal of promise for modeling complex tasks in process control and simulation and in applications of machine perception including machine vision and electronic nose for food safety and quality control [116] .
The application of ANN for predicting the shelf life of food products in food industry is quite a new and effective approach. ANN procedure is sensitive, reliable, fast, simple and low-cost method for monitoring the authenticity of the products, which provide consumers with a safer food supply [118] .
ANN is quite a new and easy computational modeling approach used for prediction, which has become popular and accepted by food industry, researchers, scientists and students [116] .
10. Conclusion
In conclusion, the most commonly used and newly developed thin layer drying models have been discussed. The development of models and the mechanism of the drying process were explained. Various statistical approaches to determine the best model(s) for the drying processes were explained. But the most widely applied categories of thin-layer models are the semi theoretical and empirical models. Using optimization could help in predicting the best models for the processes of drying agricultural products in thin layers. Models have been found to be helpful for food industries to analyze drying kinetics of agricultural products so that the quality of these products can be maintained and assist in reducing losses incurred during bumper harvesting and processing. ANN one of the numerous modeling tools in food technology is quite a new and easy computational modeling approach used for prediction, which has become popular and accepted by the food industry, researchers, scientists and students and have found extensive utilization in solving many complex real world problems.