Download

RESEARCH ARTICLE

Mathematical modelling of lemon verbena leaves drying in a continuous flow dryer equipped with a solar pre-heating system

S. Soodmand-Moghaddam1, M. Sharifi1, H. Zareiforoush2, and H. Mobli1

1Department of Agricultural Machinery Engineering, Faculty of Agricultural Engineering and Technology, University of Tehran, Karaj, Iran; 2Department of Mechanization Engineering, Faculty of Agricultural Sciences, University of Guilan, Rasht, Iran

Corresponding Author: M. Sharifi, m.sharifi@ut.ac.ir

Received: 18 August 2019 / Accepted: 24 December 2019 / Published: 14 January 2020

Abstract

The modelling of the drying process of lemon verbena leaves in a continuous flow dryer equipped with a solar pre-heating system was performed at three levels of drying temperature (50, 40, and 30 °C) and three levels of air velocity (2, 1.5 and 1 m/s). During the experiments, lemon verbena leaves were dried to the final moisture content of 10 from 76% wet basis in the continuous flow dryer. Drying kinetic showed to drying temperature and air velocity exerted significant influence on the drying time. Also, the dried lemon verbena leaves quality was obtained by determining the essential oil content of the product after drying in different conditions in the dryer. Finally, it was observed that the highest essential oil content was maintained at a temperature of 40 °C and air velocity of 1 m/s. Suitability of 10 different mathematical drying models was used to describe drying lemon verbena leaves in this dryer. The results have shown that Midilli and Kucuk’s model can successfully predict the experimental data in all air temperatures and air velocities. In Midilli and Kucuk’s model, the amounts of R2 were above 999 × 10–3 and the amounts of root- mean-square error (RMSE) and chi-square (χ2) were less than 174 × 10–4 and 19 × 10–4.

Keywords: air velocity, air temperatures, drying kinetic, Midilli and Kucuk’s model

 

1. Introduction

Lemon verbena (Lippia citriodora Kunth.) is one of the well-known medicinal plants that has many properties, including antispasmodic, carminative, diuretic, digestive and sedative agent (da Cunha et al., 2012). The most important constituents of lemon verbena essential oil (EO) include neral, limonene and geranial (Ishkeh et al., 2019). The most important process for keeping medicinal plants is drying. Drying is one of the methods of food storage (Javed et al., 2019). The drying methods of medicinal plants have a great influence on their quality (Asekun et al., 2007). Therefore, it is important to choose the suitable drying methods for sensitive products, such as medicinal plants (Motevali et al., 2014).

Reducing moisture content (MC) reduces or stops the enzymatic or microbial activity and reduces the weight of products and facilitates the processing of products. Factors such as temperature and velocity of inlet air, initial and final humidity of the product, relative humidity and ambient air temperature affect the amount of moisture variation relative to time (Movagharnejad and Nikzad, 2007). In a research, the drying of lemon balm leaves in an indirect mode solar dryer with a double-pass air heater was investigated. According to the tests results, the temperature difference between the collector outlet and ambient declined significantly with increasing volumetric flow rate, whereas increasing the flow rate improved the collector thermal efficiency by ~20% (Shamekhi-Amiri et al., 2018).

A common method for predicting moisture loss during drying time is the empirical modelling or regression modelling, which normally measures different variables by performing empirical experiments, and the best algebraic equation between variables is selected based on the fitting of the measured data with known algebraic equations related to the kinetics of drying. Modelling the actual process of drying in terms of mathematical relations is necessary to define optimum drying conditions (Panchariya et al., 2002). Akpinar (2010) fitted the thin-layer drying data in a solar dryer with forced convection and under the open sun with natural convection of mint leaves to 10 different mathematical models. Wang and Singh’s model was found as the best model describing the thin layer drying behaviour of mint leaves (Akpinar, 2010).

In a study, 14 mathematical models of mint and thymus thin layer drying were tested to specify the suitable model in an indirect-mode forced convection solar dryer. The Midilli and Kucuk’s model and Page and modified Page models were found for describing the drying curves of the mint and thymus, respectively (El-Sebaii and Shalaby, 2013).

Aral and Beşe (2016) investigated the thin-layer drying characteristics of hawthorn fruit in a convective dryer at different air temperatures and air velocities. The results of this study showed that drying time decreases with increasing air temperature and air velocity. Also, Midilli et al.’s model was obtained as the best model for drying behaviour in this study (Aral and Beşe, 2016).

Jayatunga and Amarasinghe (2018) examined the drying kinetics of black pepper in a conventional spouted bed dryer and a spouted bed dryer fitted with a non-porous draft tube. Their results showed that the highest drying rate was achieved when black pepper was dried without draft tube at the highest drying temperature. Also, a new thin layer drying model was proposed and compared with five commonly used models and the drying behaviour of black pepper in spouted bed dryer was best fitted to (Verma et al., 1985) model and the proposed model (Jayatunga and Amarasinghe, 2018). In a study, seven different mathematical drying models were compared for describing thin layer drying of cashew kernels in a batch dryer. The results showed that the Page model was found to be the best for describing the drying behaviour of cashew kernels (Asiru et al., 2013).

Knowing the drying behaviour is important in designing, simulating and optimising the drying process. Several studies have been conducted on the drying of agricultural products, including vegetables, such as drying savory leaves (Arslan and Özcan, 2012), drying turmeric rhizomes in a solar conduction dryer (Borah et al., 2015) and drying of thyme (Thymus Vulgaris L.) (Doymaz, 2011). Drying kinetic study is necessary to obtain information about the time required for drying and select the appropriate drying model (Baini and Langrish, 2008).

Due to the possibility of undesirable changes in the quality of food because of drying, its control is of particular importance; therefore, in order to maintain the safety of products, their moisture should be reduced to a certain amount. To achieve this goal, it is necessary to model the drying of various agricultural products so that the model of drying out the product can be predicted based on the pattern obtained from the model.

As conventional methods of drying leaves plants have many problems, such as lower product quality, in this article we used a continuous flow dryer. In this type of dryer, two sources of solar energy and natural gas were used to provide heat. Due to the importance and sensitivity of the lemon verbena medicinal plant to drying, the aim of this study was to develop a continuous dryer equipped with a solar pre-heating system, to find out the suitable model and to describe the drying characteristics of lemon verbena leaves in different air temperatures and air velocities.

2. Materials and methods

Material

Lemon verbena leaves were directly bought from local producers in Rasht, Iran. Leaves were packed in polyethylene packaging and stored in a refrigerator at a temperature of 6 °C. The initial MC of the samples was determined by placing three samples of leaves in a hot air oven at a temperature of 103 °C for 24 h (Zheng et al., 2005). The average initial MC of the lemon verbena leaves was equal to 76% wet basis (w.b.).

Experimental procedure

Drying experiments were performed in a continuous dryer equipped with a solar pre-heating system at the Renewable Energy Site of the University of Guilan, Rasht, Iran (49°35′19″E longitude and 37°16′33″N latitude). The dryer consists of air relative humidity/temperature sensors, solar water heater, gas water heater, continuous flow dryer, heat exchanger, water pump and centrifugal air blower (Figure 1). The air temperature and relative humidity in the drying chamber were measured by an air relative humidity/temperature sensor (Model AM2301, China).



Fig 1

Figure 1. Schematic view of the continuous dryer equipped with solar pre-heating system.

The variation of the recorded temperatures from the set points during the experiments was ±1 °C. The air velocity of the input air was measured by a turbine-type anemometer (Lutron Model, 80 AM, Taiwan). The dryer was a continuous flow type with four sets of conveyor belts placed above each other. The speed of the conveyor belt was adjusted using a single-phase inverter (LG Model IC5, Korea). Based on some pretests, the speed of the conveyor was adjusted to 43.1 × 10–2 m/s.

The movement of the drive pulley of the conveyor belts was performed by an electric motor (0. 75 kW). For the production of hot air and drying, the product is used in conventional dryers from the thermal energy of a gas burner. A solar water heater combined with a gas water heater (Butane Model 3318IF, Iran) was used to heat the water. The heated water was then entered into a heat exchanger (Hepaco Model HP-1500, Iran) using pipelines. The airflow was blown into the water pipe by a blower to direct the heat of the pipes into the drying chamber. A portion of the required thermal energy was provided by a thermosiphon type of solar water heater having the capacity of 150 L. The four rows of conveyor belts moved in the opposite direction in the drying chamber (with dimensions of 150 × 150 × 80 cm3). The centrifugal blower was able to change the amount of inlet air by a single-phase inverter (LG Model IC5, Korea).

After the dryer reached desired conditions, about 150 g of the samples was put into the dryer chamber for drying. Eventually, the MC of the lemon verbena reached 10% (w.b.) for each experiment. In this research, drying of the product was performed at 30, 40 and 50 °C temperatures and air velocities of 2, 1.5 and 1 m/s. The leaves were weighed (every 10 min) by a digital balance (A&D Model GF6000, Japan).

Drying calculation

The MC and moisture ratio of lemon verbena leaves samples during drying experiments was calculated using the following equations (Omolola et al., 2018; Erenturk et al., 2004):

QAS-12-658-E1.jpg

QAS-12-658-E2.jpg

where MC, Wo and Wd are the MC (%), weight (g) of the sample and dried sample, respectively, and Mt and Mo are the MC at any drying time and the initial MC, respectively. This equation can be simplified to Mt/Mo because the equilibrium MC value (Me) is small compared with that of Mt and Mo (Essalhi et al., 2018; Rabha et al., 2017). Also, Equation (3) was used to calculate the drying rate during leaves drying experiments (Labed et al., 2016).

QAS-12-658-E3.jpg

Extraction of the EO

Medicinal plants have unique and valuable properties. EO derived from medicinal plants had long been used in the food, pharmaceutical, perfume and cosmetics industries (Sari, 2018). One of the most important methods of preserving medicinal plants is drying (Bhardwaj et al., 2019). The method and temperature of drying of these plants have a great influence on their EO content and quality. For this reason, in the present study, the amount of EO of lemon verbena leaves after drying was obtained for each sample at different drying temperatures and speeds.

For EO extraction by the hydro-distillation method, 30 g of each sample was placed into a Clevenger with 450 ml of water for 3 h (Abd El-Gaber et al., 2018). Isolation of EO from water due to density difference was easy (Saleh et al., 2019). The EO content of each sample after weighing by a digital balance was calculated using Equation (4).

QAS-12-658-E4.jpg

Mathematical modelling

Ten drying equations given in Table 1 were evaluated to select the most suitable model for describing the drying curves of lemon verbena leaves. Regression analysis was performed using MATLAB and SPSS 11.5.1 software package. In order for the superior model to be chosen, the value of the coefficient of determination must be high and the value of RMSE and χ2 must be low (Torki-Harchegani et al., 2016). The calculations were performed using Equations (5), (6) and (7) (Akpinar, 2010):

Table 1. Mathematical models applied to the moisture ratio curves.
Models Equations References
Two-term exponential MR = a · e(−kt) + (1− a)·e(−kat) Badaoui et al. (2019)
Henderson and Pabis MR = a · e(−kt) Mghazli et al. (2017)
Logarithmic MR = a · e(−kt) + c Motevali et al. (2016)
Diffusion approach MR = a · e(−kt) + (1− a)·e(−kbt) Vijayan et al. (2016)
Midilli and Kucuk MR = a · e(−kt)n + bt El-Sebaii and Shalaby (2013)
Newton MR = e(−kt) Yogendrasasidhar and Setty (2018)
Modified Page MR e(−(kt)n) Duffie and Beckman (2013)
Two-term MR = a · e(−kt) + b · (−kt) (Torki-Harchegani et al. (2016)
Page MR = e(−ktn) Vijayan et al. (2016)
Wang and Singh MR = 1+ at + bt2 Aregbesola et al. (2015)

QAS-12-658-E5.jpg

QAS-12-658-E6.jpg

QAS-12-658-E7.jpg

where MRpre,i and MRexp,i are the predicted and experimental moisture ratio, respectively, and n and N are the number of the drying constants and observations, respectively.

3. Results and discussion

Drying kinetics

Lemon verbena leaves were dried to in a continuous dryer equipped with a solar pre-heating system at different temperatures (30, 40 and 50 °C) and different air velocities (2, 1.5 and 1 m/s). Experimental data obtained from the drying process of lemon verbena leaves were converted into MC. The MC diagram of the lemon verbena leaves during the drying period is shown in Figure 2. From the figure, it is clear that with an increase in the air temperature from 30 to 50 °C, the slope of the MC curve increases.



Fig 2

Figure 2. Variation of moisture content as a function of time for different temperatures and air velocities of (A) 2, (B) 1.5 and (C) 1 m/s.

Air temperature is an important parameter in drying the products. Therefore, with increasing air temperature, the moisture removal increased and, as shown in Figure 3, the drying time also significantly decreased (Akhondi et al., 2011; Mujaffar and John, 2018). It was also observed that with an increase in the air velocity from 1 to 2 m/s, the drying time decreased. As shown in Figure 3, the minimum drying time of leaves (110 min) was obtained at 50 °C and an air velocity of 2 m/s. The temperature and air velocity of the drying have a great influence on the drying time of the product (Ertekin and Yaldiz, 2004). With increasing air velocity, the moisture transfer rate increases (Motevali et al., 2010). At air velocity of 2 m/s and with an increase in the air temperature from 30 to 50 °C, the drying time reduced by 26.66%, while at air velocities of 1.5 and 1 m/s, the drying time was reduced by 25 and 23.52%, respectively. With an increase in air velocity from 1 to 1.5 and 2 m/s, the drying time was reduced by 5.888 and 11.76%, respectively. Drying time is one of the important factors in preserving the quality of medicinal plants.



Fig 3

Figure 3. A comparison of the drying time of lemon verbena leaves at different air velocities.

The resulted MC curves showed that the drying time and the slope of the curved surface are inversely related. At the beginning of the drying process, the speed of moisture loss is higher due to increased moisture release from the inside to the product surface. At the end of the drying period, the leaves are dampened with a milder gradient and ultimately the sample is weighed down the constant (Doymaz et al., 2006).

Figure 4 shows the changes in the drying rate versus MC at temperatures of 30, 40 and 50 °C and air velocities of 1, 1.5 and 2 m/s. The process of leaves drying occurred in the falling rate period.



Fig 4

Figure 4. Variation in the drying rate versus moisture content of lemon verbena leaves at air velocities of (A) 2, (B) 1.5 and (C) 1 m/s.

As shown in the figure, the drying rate increased when the drying air temperature increased, and at air temperature of 50 ° C at all air velocities, the highest drying rate was obtained. The reason for this increase could be the increased heat transfer potential between the drying air and product leaves.

Quality evaluation

Table 2 shows the effect of three air velocities of 1, 1.5 and 2 m/s and temperatures of 30, 40 and 50 °C on the EO content of lemon verbena leaves. As can be seen from the data in the table, due to the sensitivity of the EO to high temperature, the EO content decreases with increasing temperature.

Table 2. The essential oil content (%) of lemon verbena leaves at different drying temperatures and air velocities.
Air velocity (m/s) Air temperature (°C)
  30 40 50
1 0.47 0.43 0.33
1.5 0.53 0.50 0.37
2 0.57 0.60 0.40

According to the results, the EO content at a temperature of 30 °C was less than that of 40 °C. This may be due to the sudden accumulation of EO in its storage cells. The lowest EO content was obtained at drying temperature of 50 °C and air velocity of 1 m/s (0.33%) and the highest EO content was obtained at drying temperature of 40 °C and air velocity of 2 m/s (0.60%). According to the results, drying air temperature and air velocity affected EO content (Pirbalouti et al., 2013; Torki-Harchegani et al., 2016).

Mathematical modelling of drying curves

In this study, different regression models were evaluated with respect to R2, RMSE and χ2 values. The criterion for choosing the superior model was a greater amount of R2 and lower values of χ2 and RMSE. Table 3 presents the results of statistical indices obtained from different models for lemon verbena drying at different levels of temperature and air velocity. As can be seen from the data, the highest matching with experimental data was obtained in Midilli and Kucuk’s model for drying of lemon verbena leaves. The arbitrary constant of this model for all drying conditions (temperature and velocity) is given in Table 4.

Table 3. Accuracy of fitting models based on different statistical indices.
Temperature 30 °C 40 °C 50 °C
Model name RMSE χ2 R2 RMSE χ2 R2 RMSE χ2 R2
1 m/s                  
Two-term exponential 78 × 10-4 132 × 10-4 9,950 × 10-4 3,078 × 10-4 120 × 10-4 9,959 × 10-4 3,003 × 10-4 90 × 10-4 9,920 × 10-4
Henderson and Pabis 3,001 × 10-4 4,679 × 10-4 5,890 × 10-4 3,012 × 10-4 2,004 × 10-4 8,907 × 10-4 3,129 × 10-4 2,348 × 10-4 7,987 × 10-4
Logarithmic 2,104 × 10-4 450 × 10-4 9,567 × 10-4 2,978 × 10-4 658 × 10-4 9,485 × 10-4 3,036 × 10-4 246 × 10-4 9,568 × 10-4
Diffusion approach 0 126 × 10-4 0 2,584 × 10-4 15 × 10-4 6,539 × 10-4 3,002 × 10-4 17 × 10-4 5,377 × 10-4
Midilli and Kucuk 87 × 10-4 19 × 10-4 9,994 × 10-4 71 × 10-4 18 × 10-4 9,990 × 10-4 173 × 10-4 10 × 10-4 9,995 × 10-4
Newton 2,348 × 10-4 2.8976 5,468 × 10-4 2,786 × 10-4 7,655 × 10-4 8,134 × 10-4 3,100 × 10-4 6,754 × 10-4 8,076 × 10-4
Modified Page 165 × 10-4 18 × 10-4 9,896 × 10-4 3,187 × 10-4 16 × 10-4 9,980 × 10-4 3,658 × 10-4 38 × 10-4 9,879 × 10-4
Two-term 3,015 × 10-4 126 × 10-4 9,979 × 10-4 3,006 × 10-4 66 × 10-4 9,965 × 10-4 3,457 × 10-4 21 × 10-4 9,891 × 10-4
Page 1,116 × 10-4 392 × 10-4 8,746 × 10-4 1,053 × 10-4 164 × 10-4 9,757 × 10-4 2,256 × 10-4 244 × 10-4 9,431 × 10-4
Wang and Singh 0 2.998 0 3,012 × 10-4 4,358 × 10-4 5,678 × 10-4 3,014 × 10-4 3,178 × 10-4 7,034 × 10-4
1.5 m/s                  
Two-term exponential 158 × 10-4 37 × 10-4 9,974 × 10-4 3,175 × 10-4 66 × 10-4 9,965 × 10-4 3,149 × 10-4 96 × 10-4 9,915 × 10-4
Henderson and Pabis 2,627 × 10-4 3,737 × 10-4 6,912 × 10-4 2,981 × 10-4 1,381 × 10-4 8,896 × 10-4 2,849 × 10-4 1,986 × 10-4 8,147 × 10-4
Logarithmic 3,103 × 10-4 490 × 10-4 9,649 × 10-4 3,064 × 10-4 689 × 10-4 9,385 × 10-4 3,127 × 10-4 303 × 10-4 9,756 × 10-4
Diffusion approach 0 26 × 10-4 0 2,634 × 10-4 4 × 10-4 6,939 × 10-4 2,383 × 10-4 9 × 10-4 5,678 × 10-4
Midilli and Kucuk 59 × 10-4 18 × 10-4 9,996 × 10-4 61 × 10-4 16 × 10-4 9,993 × 10-4 162 × 10-4 4 × 10-4 9,996 × 10-4
Newton 1,902 × 10-4 2.7311 3,617 × 10-4 2,652 × 10-4 7,785 × 10-4 7,034 × 10-4 2,652 × 10-4 7,785 × 10-4 7,034 × 10-4
Modified Page 145 × 10-4 12 × 10-4 9,931 × 10-4 3,162 × 10-4 14 × 10-4 9,976 × 10-4 3,159 × 10-4 26 × 10-4 9,979 × 10-4
Two-term 3,158 × 10-4 27 × 10-4 9,964 × 10-4 3,156 × 10-4 76 × 10-4 9,958 × 10-4 3,161 × 10-4 17 × 10-4 9,906 × 10-4
Page 2,108 × 10-4 291 × 10-4 9,632 × 10-4 2,153 × 10-4 65 × 10-4 9,963 × 10-4 2,136 × 10-4 146 × 10-4 9,830 × 10-4
Wang and Singh 0 2.9576 0 2,734 × 10-4 3,594 × 10-4 6,999 × 10-4 2,686 × 10-4 2,735 × 10-4 7,206 × 10-4
2 m/s                  
Two term exponential 169 × 10-4 40 × 10-4 9,886 × 10-4 3,195 × 10-4 70 × 10-4 9,884 × 10-4 3,019 × 10-4 86 × 10-4 9,930 × 10-4
Henderson and Pabis 3,101 × 10-4 3,897 × 10-4 6,012 × 10-4 3,283 × 10-4 2,581 × 10-4 7,895 × 10-4 3,485 × 10-4 2,316 × 10-4 7,645 × 10-4
Logarithmic 3,870 × 10-4 679 × 10-4 9,438 × 10-4 3,184 × 10-4 729 × 10-4 9,168 × 10-4 3,547 × 10-4 403 × 10-4 9,470 × 10-4
Diffusion approach 0 126 × 10-4 0 2,587 × 10-4 12 × 10-4 7,839 × 10-4 2,012 × 10-4 11 × 10-4 6,675 × 10-4
Midilli and Kucuk 51 × 10-4 15 × 10-4 9,998 × 10-4 60 × 10-4 16 × 10-4 9,995 × 10-4 174 × 10-4 9 × 10-4 9,994 × 10-4
Newton 987 × 10-4 1.8759 5,648 × 10-4 2,027 × 10-4 6,548 × 10-4 8,237 × 10-4 1,627 × 10-4 6,775 × 10-4 8,036 × 10-4
Modified Page 120 × 10-4 10 × 10-4 9,945 × 10-4 2,154 × 10-4 12 × 10-4 9,980 × 10-4 2,127 × 10-4 16 × 10-4 9,989 × 10-4
Two-term 2,984 × 10-4 23 × 10-4 9,970 × 10-4 3,026 × 10-4 56 × 10-4 9,968 × 10-4 2,261 × 10-4 15 × 10-4 9,936 × 10-4
Page 3,024 × 10-4 1,247 × 10-4 9,012 × 10-4 2,587 × 10-4 146 × 10-4 9,087 × 10-4 2,547 × 10-4 268 × 10-4 9,580 × 10-4
Wang and Singh 0 2.7547 0 3,346 × 10-4 4,337 × 10-4 5,698 × 10-4 3,623 × 10-4 3,765 × 10-4 6,206 × 10-4

 

Table 4. Arbitrary constant of the Midilli and Kucuk’s model for different temperatures and air velocities.
V (m/s) T (°C) n b k a
  30 1.126 −0.00021 0.00576 0.9913
1 40 1.129 −0.00001 0.00904 0.9954
  50 1.193 −0.00076 0.01845 0.9864
  30 1.034 −0.00011 0.00875 0.9924
1.5 40 1.127 −0.00000 0.00975 0.9962
  50 1.184 −0.00043 0.01743 0.9878
  30 1.128 −0.00012 0.00764 0.9986
2 40 1.124 −0.00001 0.00986 0.9972
  50 1.192 −0.00102 0.02573 0.9789

The variation of moisture ratio of lemon verbena leaves with time for 2, 1.5 and 1 m/s air velocities is shown in Figure 5. The results have shown that Midilli and Kucuk’s model can successfully predict the experimental data. From the plots in the figure, it is observed that Midilli and Kucuk’s model exhibits a good fit to the experimental data. Also other researchers have confirmed Midilli and Kucuk’s model as the suitable model to describe the drying behaviour of different products, such as drying kinetics of pumpkin slice (Benseddik et al., 2018), yacon slices (Shi et al., 2013), sage leaves in a cabinet dryer (Doymaz and Karasu, 2018) and apple slices (Meisami-Asl et al., 2010).



Fig 5

Figure 5. Variation of moisture ratio (MR) experimental and MR predicted with drying time at air velocities of (A) 2, (B) 1.5 and (C) 1 m/s and different air temperatures.

4. Conclusions

In this study, a continuous flow dryer equipped with a solar pre-heating system was used for drying the lemon verbena leaves. The developed system has the proper ability to provide the primary thermal energy for drying lemon verbena leaves. The benefits of using solar energy for drying include reducing pollution emission and fossil fuel consumption. The temperature of the air for drying the lemon verbena was 30, 40 and 50 °C during the experiments, which was sufficient for leaves drying. The increase in drying temperature and air velocity reduced the drying time of lemon verbena leaves. As the air velocity in the dryer increased, the drying rate increased at all three drying temperatures. Also the MC of lemon verbena leaves was reduced to 10% w.b. in this dryer. The resulted MC curves showed that the time of drying and the slope of the curved surface were inversely related.

The dried lemon verbena leaves quality was obtained by determining the EO content of the product after drying in different dryer conditions. Finally, it was observed that the highest EO content was maintained at a temperature of 40 °C and air velocity of 1 m/s. A number of mathematical models were tested for the prediction of drying lemon verbena leaves in this dryer. According to statistical analysis, Midilli and Kucuk’s model was selected due to the high coefficient of determination, low RMSE and low χ2and also due to the simplicity of the relationship between kinetics and prediction of drying process for lemon verbena in different air velocities and temperatures. In these models, the amount of R2 was above 999 × 10–3 and the amounts of RMSE and χ2 were less than 174 × 10–4 and 19 × 10–4, respectively. Therefore, among the models examined in this article, Midilli and Kucuk’s model was selected as the appropriate model for the process of lemon verbena drying in a dryer equipped with a solar pre-heating system for all drying conditions.

Conflict of interest

The authors declare no conflicts of interest with respect to research, authorship and/or publication of this article.

Compliance with ethical standards

This article followed all ethical standards for a research without direct contact with human or animal subjects.

Funding

This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.

References

Abd El-Gaber, A.S., El Gendy, A.N.G., Elkhateeb, A., Saleh, I.A. and El-Seedi, H.R., 2018. Microwave extraction of essential oil from Anastatica hierochuntica (L): comparison with conventional hydro-distillation and steam distillation. Journal of Essential Oil Bearing Plants 21: 1003–1010. https://doi.org/10.1080/0972060X.2018.1504695
Akhondi, E., Kazemi, A. and Maghsoodi, V., 2011. Determination of a suitable thin layer drying curve model for saffron (Crocus sativus L) stigmas in an infrared dryer. Scientia Iranica 18: 1397–1401. https://doi.org/10.1016/j.scient.2011.08.014
Akpinar, E.K., 2010. Drying of mint leaves in a solar dryer and under open sun: modelling, performance analyses. Energy Conversion and Management 51: 2407–2418. https://doi.org/10.1016/j.enconman.2010.05.005
Aral, S. and Beşe, A.V., 2016. Convective drying of hawthorn fruit (Crataegus spp.): effect of experimental parameters on drying kinetics, color, shrinkage, and rehydration capacity. Food Chemistry 210: 577–584. https://doi.org/10.1016/j.foodchem.2016.04.128
Aregbesola, O.A., Ogunsina, B.S., Sofolahan, A.E. and Chime, N.N., 2015. Mathematical modeling of thin layer drying characteristics of dika (Irvingia gabonensis) nuts and kernels. Nigerian Food Journal 33: 83–89. https://doi.org/10.1016/j.nifoj.2015.04.012
Arslan, D. and Özcan, M.M., 2012. Evaluation of drying methods with respect to drying kinetics, mineral content, and color characteristics of savory leaves. Food and Bioprocess Technology 5: 983–991. https://doi.org/10.1007/s11947-010-0498-y
Asekun, O.T., Grierson, D.S. and Afolayan, A.J., 2007. Effects of drying methods on the quality and quantity of the essential oil of Mentha longifolia L. subsp. Capensis. Food Chemistry 101: 995–998. https://doi.org/10.1016/j.foodchem.2006.02.052
Asiru, W.B., Raji, A.O., Igbeka, J.C. and Elemo, G.N., 2013. Mathematical modelling of thin layer dried cashew kernels. Nigerian Food Journal 31: 106–112. https://doi.org/10.1016/S0189-7241(15)30083-7
Badaoui, O., Hanini, S., Djebli, A., Haddad, B. and Benhamou, A., 2019. Experimental and modelling study of tomato pomace waste drying in a new solar greenhouse: evaluation of new drying models. Renewable Energy 133: 144–155. https://doi.org/10.1016/j.renene.2018.10.020
Baini, R. and Langrish, T.A.G., 2008. An assessment of the mechanisms for diffusion in the drying of bananas. Journal of Food Engineering 85: 201–214. https://doi.org/10.1016/j.jfoodeng.2007.06.035
Benseddik, A., Azzi, A., Zidoune, M.N. and Allaf, K., 2018. Mathematical empirical models of thin-layer airflow drying kinetics of pumpkin slice. Engineering in Agriculture, Environment and Food 11: 220–231. https://doi.org/10.1016/j.eaef.2018.07.003
Bhardwaj, A.K., Kumar, R. and Chauhan, R., 2019. Experimental investigation of the performance of a novel solar dryer for drying medicinal plants in Western Himalayan region. Solar Energy 177: 395–407. https://doi.org/10.1016/j.solener.2018.11.007
Borah, A., Hazarika, K. and Khayer, S.M., 2015. Drying kinetics of whole and sliced turmeric rhizomes (Curcuma longa L.) in a solar conduction dryer. Inf. Process. Agric. 2: 85–92. https://doi.org/10.1016/j.inpa.2015.06.002
da Cunha, A.P., Nogueira, M.T., Roque, O.R. and Barroso, J.M.G., 2012. Plantas aromáticas e óleos essenciais: composiçã;o e apli¬cações, Fundaçã;o Calouste Gulbenkian publisher, Português, pp. 35–39.
Doymaz, I., 2011. Drying of thyme (Thymus Vulgaris L.) and selection of a suitable thin-layer drying model. Journal of Food Processing and Preservation 35: 458–465. https://doi.org/10.1111/j.1745-4549.2010.00488.x
Doymaz, İ. and Karasu, S., 2018. Effect of air temperature on drying kinetics, colour changes and total phenolic content of sage leaves (Salvia officinalis). Quality Assurance and Safety of Crops & Foods 10: 269–276. https://doi.org/10.3920/QAS2017.1257
Doymaz, I., Tugrul, N. and Pala, M., 2006. Drying characteristics of dill and parsley leaves. Journal of Food Engineering 77: 559–565. https://doi.org/10.1016/j.jfoodeng.2005.06.070
Duffie, J.A. and Beckman, W.A., 2013. Solar engineering of thermal processes, John Wiley & Sons, pp. 98–100 https://doi.org/10.1002/9781118671603
El-Sebaii, A.A. and Shalaby, S.M., 2013. Experimental investigation of an indirect-mode forced convection solar dryer for drying thymus and mint. Energy Conversion and Management 74: 109–116. https://doi.org/10.1016/j.enconman.2013.05.006
Erenturk, S., Gulaboglu, M.S. and Gultekin, S., 2004. The thin-layer drying characteristics of rosehip. Biosystems Engineering 89: 159–166. https://doi.org/10.1016/j.biosystemseng.2004.06.002
Ertekin, C. and Yaldiz, O., 2004. Drying of eggplant and selection of a suitable thin layer drying model. Journal of Food Engineering 63: 349–359. https://doi.org/10.1016/j.jfoodeng.2003.08.007
Essalhi, H., Benchrifa, M., Tadili, R. and Bargach, M.N., 2018. Experimental and theoretical analysis of drying grapes under an indirect solar dryer and in open sun. Innovative Food Science and Emerging Technologies 49: 58–64. https://doi.org/10.1016/j.ifset.2018.08.002
Ishkeh, S.R., Asghari, M., Shirzad, H., Alirezalu, A. and Ghasemi, G., 2019. Lemon verbena (Lippia citrodora) essential oil effects on antioxidant capacity and phytochemical content of raspberry (Rubus ulmifolius subsp. sanctus). Scientia Horticulturae 248: 297–304 (Amsterdam).
Javed, H.U., Wang, D., Wu, G.-F., Kaleem, Q.M., Duan, C.-Q. and Shi, Y., 2019. Post-storage changes of volatile compounds in air-and sun-dried raisins with different packaging materials using HS-SPME with GC/MS. Food Research International 119: 23–33. https://doi.org/10.1016/j.foodres.2019.01.007
Jayatunga, G.K. and Amarasinghe, B., 2018. Mathematical modeling of drying kinetics of black pepper in a spouted bed dryer with and without non porous draft tubes. 2nd International Conference on Trends in Multidisciplinary Business and Economic Research (TMBER), March 23–24, 2018, Sydney Australia, pp. 118–134.
Labed, A., Moummi, N., Aoues, K. and Benchabane, A., 2016. Solar drying of henna (Lawsonia inermis) using different models of solar flat plate collectors: an experimental investigation in the region of Biskra (Algeria). Journal of Cleaner Production 112: 2545–2552. https://doi.org/10.1016/j.jclepro.2015.10.058
Meisami-Asl, E., Rafiee, S., Keyhani, A. and Tabatabaeefar, A., 2010. Determination of suitable thin layer drying curve model for apple slices (variety-Golab). Plant Omics 3: 103.
Mghazli, S., Ouhammou, M., Hidar, N., Lahnine, L., Idlimam, A. and Mahrouz, M., 2017. Drying characteristics and kinetics solar drying of Moroccan rosemary leaves. Renewable Energy 108: 303–310. https://doi.org/10.1016/j.renene.2017.02.022
Motevali, A., Minaei, S., Banakar, A., Ghobadian, B. and Darvishi, H., 2016. Energy analyses and drying kinetics of chamomile leaves in microwave-convective dryer. Journal of the Saudi Society of Agricultural Sciences 15: 179–187. https://doi.org/10.1016/j.jssas.2014.11.003
Motevali, A., Minaei, S., Banakar, A., Ghobadian, B. and Khoshtaghaza, M.H., 2014. Comparison of energy parameters in various dryers. Energy Conversion and Management 87: 711–725. https://doi.org/10.1016/j.enconman.2014.07.012
Motevali, A., Minaei, S., Khoshtaghaza, M.H., Kazemi, M. and Nikbakht, A.M., 2010. Drying of pomegranate arils: comparison of predictions from mathematical models and neural networks. International Journal of Food Engineering 6(3): 1556–3758. https://doi.org/10.2202/1556-3758.1889
Movagharnejad, K. and Nikzad, M., 2007. Modeling of tomato drying using artificial neural network. Computers and Electronics in Agriculture 59: 78–85. https://doi.org/10.1016/j.compag.2007.05.003
Mujaffar, S. and John, S., 2018. Thin-layer drying behavior of West Indian lemongrass (Cymbopogan citratus) leaves. Food Science & Nutrition 6: 1085–1099. https://doi.org/10.1002/fsn3.642
Omolola, A.O., Kapila, P.F. and Silungwe, H.M., 2018. Mathematical modeling of drying characteristics of Jew’s mallow (Corchorus olitorius) leaves. Information Processing in Agriculture 6(1): 109–115. https://doi.org/10.1016/j.inpa.2018.08.003.
Panchariya, P.C., Popovic, D. and Sharma, A.L., 2002. Thin-layer modelling of black tea drying process. Journal of Food Engineering 52: 349–357. https://doi.org/10.1016/S0260-8774(01)00126-1
Pirbalouti, A.G., Mahdad, E. and Craker, L., 2013. Effects of drying methods on qualitative and quantitative properties of essential oil of two basil landraces. Food Chemistry 141: 2440–2449. https://doi.org/10.1016/j.foodchem.2013.05.098
Rabha, D.K., Muthukumar, P. and Somayaji, C., 2017. Experimental investigation of thin layer drying kinetics of ghost chilli pepper (Capsicum Chinense Jacq.) dried in a forced convection solar tunnel dryer. Renewable Energy 105: 583–589. https://doi.org/10.1016/j.renene.2016.12.091
Saleh, I.A., El Gendy, A.N.G., Afifi, M.A. and El-Seedi, H.R., 2019. Microwave extraction of essential oils from senecio serpens GD Rowly and comparison with conventional hydro-distillation method. Journal of Essential Oil Bearing Plants 22(4): 955–961. https:// doi.org/10.1080/0972060X.2019.1649202
Sari, I., 2018. Evaluation of industrial cosmetic products and replacing them with natural products according to modern, traditional and islamic medicine. Health Biotechnology and Biopharma 1(4): 19–31. DOI: 10.22034/HBB.2018.02
Shamekhi-Amiri, S., Gorji, T.B., Gorji-Bandpy, M. and Jahanshahi, M., 2018. Drying behaviour of lemon balm leaves in an indirect double-pass packed bed forced convection solar dryer system. Case Studies in Thermal Engineering 12: 677–686. https://doi.org/10.1016/j.csite.2018.08.007
Shi, Q., Zheng, Y. and Zhao, Y., 2013. Mathematical modeling on thin-layer heat pump drying of yacon (Smallanthus sonchifolius) slices. Energy Conversion and Management 71: 208–216. https://doi.org/10.1016/j.enconman.2013.03.032
Torki-Harchegani, M., Ghanbarian, D., Pirbalouti, A.G. and Sadeghi, M., 2016. Dehydration behaviour, mathematical modelling, energy efficiency and essential oil yield of peppermint leaves undergoing microwave and hot air treatments. Renewable and Sustainable Energy Reviews 58: 407–418. https://doi.org/10.1016/j.rser.2015.12.078
Verma, L.R., Bucklin, R.A., Endan, J.B. and Wratten, F.T., 1985. Effects of drying air parameters on rice drying models. Trans. ASAE 28: 296–301.
Vijayan, S., Arjunan, T.V. and Kumar, A., 2016. Mathematical modeling and performance analysis of thin layer drying of bitter gourd in sensible storage based indirect solar dryer. Innovative Food Science and Emerging Technologies 36: 59–67. https://doi.org/10.1016/j.ifset.2016.05.014
Yogendrasasidhar, D. and Setty, Y.P., 2018. Experimental studies and thin layer modeling of pearl millet using continuous multistage fluidized bed dryer staged externally. Engineering Science and Technology, an International Journal 22(2): 428–438. https://doi.org/10.1016/j.jestch.2018.10.010
Zheng, X., Jiang, Y. and Pan, Z., 2005. Drying and quality characteristics of different components of alfalfa, In: 2005 ASAE Annual Meeting. American Society of Agricultural and Biological Engineers, p. 1. August 13–15, 2005, Hong Kong, pp, 1157–1164. https://doi.org/10.1142/9789812771957_0168.