Beta Fulltext view is in preview — article structure may vary. Browse all articles
Contents
Food Science & Nutrition Technology Research Article 17 min read

Mathematical Modelling of the Drying Characteristics of Indian Spinach Using Machine Learning Approach

Jegede AO* and Isa J*
* Corresponding author
ISSN: 2574-2701  10.23880/fsnt-16000308  Received: August 08, 2023  Published: September 11, 2023
  views
 27 references
 4 figures
 8 tables
PDF
Keywords
Drying Machine Learning Models Moisture Ratio
Abstract

This study employs machine learning techniques to investigate the drying characteristics of Indian spinach leaves. Drying experiments were conducted at varying air temperatures (50°C, 60°C, 70°C) and velocities (1.2m/s and 1.4m/s) while monitoring drying time, moisture content, and important parameters. Artificial neural networks and support vector machine were used to predict drying rate and final moisture content based on conditions. Higher temperatures and air velocities lead to shorter drying times and lower moisture contents, with higher air velocities enhancing moisture diffusion and reducing energy intensity. Nine drying models were fitted, with the Page and Hii, et al. models showing strong performance. The Multi- Layer Perceptron (MLP) 3-8-1 configuration achieves high accuracy, and Support Vector Machine (SVM) models with the Radial Basis Function (RBF) kernel perform well in predicting moisture ratio. These findings offer insights into Indian spinach drying behavior and highlight the effectiveness of machine learning models. The obtained values aid in optimizing drying conditions, and the robust performance of MLP and SVM models holds potential for applications beyond drying, including diabetic retinopathy diagnosis.

Introduction

India spinach, scientifically known as Basella alba L., is a nutrient-rich leafy vegetable indigenous to India. Alongside varieties like Malabar spinach (Basella rubra L.) and Ceylon spinach (Basella cordifolia L.), it belongs to the Basellaceae family [1]. Renowned for its high vitamin content, including A, C, and K, as well as minerals like iron and calcium, India spinach is a versatile crop adaptable to various climates. It can be grown from seeds or cuttings, with a harvest period of approximately six weeks. Notably, it contains mild okra-like sap that prompts some to cook the leaves [2]. Malabar spinach, another member of the Basella species, offers significant health benefits. Its leaves are rich in vitamins A and C, calcium, potassium, and iron, surpassing regular spinach in vitamin C content and even outpacing kale in vitamin A content. The plant’s shoots are known for their febrifuge, diuretic, and laxative effects, while its roots are used to treat diarrhea. Poultices made from Malabar spinach leaves serve as applications for boils and sores. Its protein, amino acids, minerals, vitamins, and flavones make it a vital dietary and therapeutic resource [3]. Food preservation is vital for spoilage prevention and nutritional retention. Preservation methods like canning, freezing, salting, and food irradiation have emerged to combat spoilage caused by chemical and microbial agents. Thermal dehydration, an ancient method, is cost-effective and widely used. Although it can compromise quality slightly, the goal is to extend shelf life economically [4]. Drying is a cost-efficient approach to extend the shelf life of various food items. Moisture reduction is necessary to control microbial growth, measured by water activity. Drying methods, dryers, and pre- and post-processing interventions play vital roles in achieving optimal dried product quality. Techniques like osmotic dehydration and blanching enhance quality, while post-processing activities like coating and blending maintain final product quality [4]. Machine learning has revolutionized various sectors, including agriculture and food preservation. Machine learning algorithms optimize drying processes by analyzing variables like temperature, humidity, and airflow to enhance speed and product quality. Accurate prediction of drying characteristics is crucial for quality optimization, especially for India spinach, which exhibits unique drying behavior determinants [1]. Traditional mathematical modeling oversimplifies complex interactions in drying kinetics. Machine learning, using extensive datasets of drying conditions and spinach properties, offers a robust model for precise drying behavior prediction [3]. However, machine learning in drying processes, particularly for India spinach, remains underexplored due to scarce datasets and specific drying characteristics. A model merging machine learning techniques with domain expertise is needed to gain insights into Indian spinach’s drying behavior, revolutionizing drying protocols [5]. Innovative drying technologies combining techniques with physical fields like microwaves and ultrasonics are reshaping fruit and vegetable preservation, addressing challenges of energy consumption and nutrient loss. Microwave-assisted drying, known for faster drying and versatile product generation, stands out among these techniques [6]. AI’s role in the food industry has expanded, aiding sorting, parameter prediction, quality control, and safety through techniques like expert systems, artificial neural networks (ANN), adaptive neuro-fuzzy inference systems (ANFIS), and machine learning (ML). These tools enhance classification, sorting, and quality standard. The objective of this study is to model the drying characteristics of Indian spinach using machine learning approach.

Materials and Methods

Sample Preparation and Materials used

Fresh Indian spinach leaves were procured from Shasha market in Akure, Southwestern region of Nigeria before the experiment. The healthy leaves were selected and washed with tap water and the surface moisture on the sample surface was gently removed with blotting paper. The materials and equipment that were used for the study include: Indian Spinach leaves, drying system (mechanical dyer), knife, Tray, Small bowl, Stopwatch, Salt, Beakers, measuring cylinder, weighing balance, Distilled water, Sample holder, Blotting paper and Hand towel.

Drying Experiment Procedure

The experiment was carried out with the mechanical dryer in Processing and Storage Laboratory of the department of Agricultural and Environmental Engineering, Federal University of Technology Akure, Nigeria. The initial moisture contents of the leave sample of Indian Spinach were determined using the standard oven drying method. The dryer was allowed to reach steady-state conditions. The drying experiments were carried out at 50, 60, and 70 drying air temperatures 1.2 m/s and 1.4 m/s. The leaves were uniformly spread on the stainless steel trays of the dryer and the variation in the weight of the sample were monitored at 30 minutes intervals until final moisture content were attained. During the experiment, important parameters affecting the dryer performance such as air temperature, relative humidity, pressure, moisture content, and air velocity, ambient temperature and the inlet and outlet temperatures of drying air in the dryer chamber were measured and recorded. The dried sample was packed in an air tight container for further experiment or analysis. Determination of Drying Parameters Determination of Moisture Content: The moisture content of the sample during drying was determined using equation 1 for wet.

100 w d w w W W M W − = × (1)

Where d M is the moisture content of the sample on dry basis (%), w M is the moisture content of the sample on wet basis (%), while w W is the weight of wet sample (g) and d W is the weight of dry sample (g) Determination Moisture Ratio: Drying ratio was the ratio of the final weight of dried samples to its initial weight. It was dimensionless.

Determination of Drying Rate: The determination of the drying rate was carried out by measuring the weight change during the drying process. During drying, there is loss of moisture with time hence the drying rate indicates the rate of loss of this moisture.

The drying rate was calculated as follows:

o F M M dM dt t − = (2)

1 Where Mo is the initial moisture content, MF is the final moisture content, % d.b and t1 is the time.

Determination of Effective Moisture Diffusivity: Effective moisture diffusivity, Deff is defined to describe the rate of moisture movement, no matter which mechanism is involved. A complete drying profile consists of the first stage of drying, a constant-rate period, and a falling rate period. The unsteady state diffusion of moisture by Fick’s second law was used to determine it;

2 eff M D M t ∂ = ∇ ∂

(3) Calculation of Deff using the Fick’s second law is a tool for describing the drying process and possible mechanisms for the transfer of moisture within food products.

Where M is the moisture content (kg water/kg dry matter), t is the drying time and, Deff is the effective diffusivity (m2/s).

Determination of activation energy: The Activation Energy was calculated by using an Arrhenius type equation (Lopez et al., 2000):

exp o Ea D D RTa   = −    

(4) Where Ea is the energy of activation (kJ/mol), R is universal gas constant (8.3143 kJ/mol), Ta is absolute air temperature (K), and Do is the pre-exponential factor of the Arrhenius equation (m2/s)

Machine Language Approach

To estimate the drying rate and final moisture content of the spinach based on the drying circumstances, ANN (artificial neural network) and SVM (support vector machine) machine learning approaches were created using MATLAB 2020b. The model was trained, validated, and tested using, respectively, 70%, 15%, and 15%. To find the model’s ideal hyper parameter, trial and error methods were employed. To confirm that the model can accurately predict the drying properties of Indian spinach, the coefficient of determination, mean squared error, and mean absolute error were used to measure the model’s accuracy.

Modeling of the Drying Process

For fitting the drying characteristics of Indian spinach, some models was used to estimate the parameters associated with each model. The best fitted model was determined by the highest coefficient of determination (R2) and the lowest values of the root mean square error (RMSE) and chi-square (χ2) amongst the predicted and experimental results. Then these data were fitted to the selected models including the Page model, Midilli, et al. model, Logistic model, Two-term model, and Logarithmic model (Table 1). The results of nonlinear regression analyses and statistical analyses of the different models including the constants of the models and the comparison criteria will be used to evaluate goodness of fit namely, R2, χ2, RMSE and SSE.

ModelsEquationReferences
Midilli, et al.MR = a exp(-ktn) + btMidilli A and Kucuk H [7]
PageMR = exp(-ktn)Kaleemullah S and Kailappan R [8]
Wang and SinghMR = 1 + bt + at2Wang CY and Singh RP [9]
LogarithmicMR = a exp(-kt) + cOnwude, et al. [10]
Henderson and PabisMR = a.exp(-kt)Kashaninejad, et al. [11]
PelegMR= 1−𝑡(𝑎+𝑏𝑡)Da Silva, et al. [12,13]
Hii, et al.𝑀𝑅= 𝑎℮−𝑘𝑡𝑛+𝑐𝐸𝑥𝑝(−𝑔𝑡𝑛)Hii, et al. [14]
Diffusion approach𝑀𝑅= 𝑎Exp (−𝑘𝑡) +(1−𝑎)𝐸𝑥𝑝(−𝑔𝑡)Wang, et al. (2007)
Verma, et al.𝑀𝑅= 𝑎Exp (−𝑘𝑡) +(1+𝑎)𝐸𝑥𝑝(−𝑔𝑡)Yaldiz O and Ertekin C [15]

Table 1: Thin layer model equations.

Results and Discussions

Effect of Temperature on Drying Behaviour

The drying behavior of Indian spinach at various temperatures (70°C, 60°C, and 50°C) with an airflow velocity of 1.2 m/s was analyzed in Figure 1. The impact of temperature on drying time and moisture content was assessed. At 70°C, the drying process began with a high initial moisture content of 93.27%, and as time progressed, drying time increased while moisture content steadily decreased to 27.15% in 147.65 minutes. Beyond this point, drying time stabilized, and moisture content continued to decrease.

Figure 1: The impact of temperature on drying time and moisture content was assessed. At 70°C, the drying process began with a high initial moisture content of 93.27%, and as time progressed, drying time increased while moisture content steadily decreased to 27.15% in 147.65 minutes. Beyond this point, drying time stabilized, and moisture content continued to decrease.
Click to enlarge
Figure 1: The impact of temperature on drying time and moisture content was assessed. At 70°C, the drying process began with a high initial moisture content of 93.27%, and as time progressed, drying time increased while moisture content steadily decreased to 27.15% in 147.65 minutes. Beyond this point, drying time stabilized, and moisture content continued to decrease.
Figure 2: Effect of temperature on the drying behaviour of Indian spinach at 1.4 m/s.
Click to enlarge
Figure 2: Effect of temperature on the drying behaviour of Indian spinach at 1.4 m/s.

Comparing the temperatures, higher values led to shorter drying times and lower moisture contents. At 60°C, both drying time and moisture content decreased compared to 70°C. Similarly, at 50°C, drying time decreased further, and moisture content reached its lowest point, implying that higher temperatures expedited moisture removal. Drying rates, representing moisture removal speed, mirrored this trend at all temperatures. The rates decreased as moisture content reduced. This observation aligned with prior research by Premi, et al. [16], Jun-Wen, et al. [17], and Rodríguez, et al. [18], confirming the influence of higher temperatures on shorter drying times and reduced moisture content. Figure 2 demonstrated the drying behavior at temperatures of 70°C, 60°C, and 50°C, with a constant air velocity of 1.4 m/s. A consistent decline in moisture content was observed as temperature decreased, highlighting the efficacy of higher temperatures in accelerating moisture removal. This phenomenon resulted from enhanced kinetic energy of water molecules, elevating evaporation rates. Additionally, drying time decreased with higher temperatures due to increased heat transfer rates, promoting faster moisture evaporation. The moisture ratio, reflecting moisture content relative to the initial amount, decreased over time at each temperature, with higher temperatures leading to swifter reduction. This trend echoed the drying time results. Drying rates increased with higher temperatures due to greater temperature differentials between spinach and air, improving moisture transfer efficiency. This finding aligned with Jun-Wen, et al. [17]. Although drying rates initially rose, they eventually declined due to decreasing moisture content, reducing the driving force for moisture transfer. This study’s outcomes were consistent with research by Ningning, et al. [19], Krokida, et al. [20], and Clemente, et al. [21], reinforcing the correlation between higher temperatures and reduced moisture content, shorter drying times, increased drying rates, and accelerated moisture removal.

Effect of Air Velocity on Drying Behaviour

The data in Figure 3 demonstrates the drying behavior of Indian spinach under different air velocities (1.4 m/s and 1.2 m/s). As drying time increases, moisture ratio decreases for both velocities, indicating moisture loss. Higher air velocity results in lower moisture ratio, indicating faster moisture evaporation and more efficient drying. Drying time is shorter at 1.4 m/s due to enhanced heat and mass transfer, accelerating moisture removal. Drying rate initially increases, then decreases as moisture content reduces. Higher air velocity leads to higher drying rate due to better heat and mass transfer. Moisture content decreases as drying time increases for both velocities. Higher air velocity leads to lower moisture content, reflecting more effective moisture evaporation.

Figure 3: Effect of air velocity on the drying behaviour of Indian spinach at 50 °C.
Click to enlarge
Figure 3: Effect of air velocity on the drying behaviour of Indian spinach at 50 °C.

Figure 4 explores air velocity’s impact on Indian spinach drying at a constant 60°C. Higher air velocity (1.4 m/s) facilitates greater moisture removal, evidenced by faster moisture content reduction. Drying time is shorter with increased air velocity, indicating accelerated drying. Moisture ratio steadily decreases at both velocities, with higher air velocity causing more rapid decrease due to enhanced evaporation. Drying rate rises over time, consistently higher at 1.4 m/s, indicating faster moisture evaporation. This aligns with a study by Luna et al. [22], indicating higher air velocity leads to quicker moisture removal with increased temperature.

Figure 4: Effect of air velocity on the drying behaviour of Indian spinach at 60 °C.
Click to enlarge
Figure 4: Effect of air velocity on the drying behaviour of Indian spinach at 60 °C.

Mathematical Modelling

The drying modeling of Indian spinach was investigated at various air velocities (1.2 m/s) and temperatures (70°C, 60°C, 50°C) using different models. Among them, the Page model demonstrated excellent performance with high R² values (0.9885 to 0.9987), indicating strong fit across temperatures. RMSE values (0.0119 to 0.0353) supported its accuracy, and low chi-square values affirmed the goodness of fit. The Hii, et al. model also showed impressive results (R²: 0.9875 to 0.9990) with low RMSE and chi-square values. Midili, et al. model had a reasonable fit (R²: 0.9538 to 0.9886) with moderate RMSE and satisfactory chi-square values. The Lewis, Verma, et al. and Diffusion models had moderate fits (R²: 0.8091 to 0.9358) with higher RMSE, suggesting room for improvement. Logarithmic, Henderson & Pabis, and Wang and Singh models also showed moderate fits (R²: 0.8240 to 0.9570), indicating reasonable average deviations, but chi-square values indicated room for enhancement. The Peleg model had high R² values (0.9461 to 0.9928) with low RMSE and chi-square values, indicating a good fit. Table 3 presented drying modeling results at different temperatures (1.4 m/s). At 70°C, the Page model excelled (R²: 0.9982, RMSE: 0.0198, Χ²: 0.0005). Similarly, at 60°C, the Page model performed well (R²: 0.9924, RMSE: 0.0367, Χ²: 0.0016). At 50°C, the Midili, et al. model was the top performer (R²: 0.9810, RMSE: 0.0536, Χ²: 0.0041). These models provide insights for optimizing Indian spinach drying under specific temperature conditions.

Air
velocity
TemperatureModelConstantRMSEΧ²
1.2 m/s70 °CLewisk=0.0060.8330.18060.0373
Pagek=1.02E-07, n=3.22030.99170.03530.0017
Midili, et al.a=0.9891, k=-0.0615, n=0.4788, b=-0.01040.96510.06820.0093
Hii, et al.a=0.5, k=1.14E-07, n=3.2078, c=0.5,
g=1.14E-07
0.99160.03360.003
Verma, et al.a=0.0831, k=0.0072, g=-0.00720.80910.16320.0426
Diffusiona=-86.8508, k=0.0181, g=-0.01780.93580.09620.0148
Wang and Singha=-0.0011, b=-1.85E-050.9570.07640.0078
Henderson & Pabisa=1.1662, k=0.00720.80910.16320.0355
Logarithmica=125.3997, k=4.09E-05, C=-124.260.92480.10010.016
Pelega=397.5017, b=-0.91290.94610.090.0108
60 °CLewisk=0.00310.88940.12790.0178
Pagek=8.96E-07, n=2.50850.99870.01190.0002
Midili, et al.a=0.9808, k=-0.0367, n=0.5065, b=-0.00550.98860.03160.0015
Hii, et al.a=0.4958, k=7.48E-07, n=2.5379,
c=0.4936, g=7.32E-07
0.9990.01050.0002
Verma, et al.a=0.075, k=0.0038, g=-0.00380.86930.10870.0158
Diffusiona=-84.0404, k=0.0095, g=-0.00930.97090.05170.0036
Wang and Singha=-0.0007, b=-6.21E-060.98940.0310.0012
Henderson & Pabisa=1.15, k=0.00380.86930.10870.0142
Logarithmica=107.9925, k=2.60E-05, C=-106.87710.95770.06080.0049
Pelega=733.0774, b=-1.11320.9820.04380.0023
50 °CLewisk=0.00120.84380.10850.0125
Pagek=1.00E-07, n=2.5950.98850.02420.0007
Midili, et al.a=0.9444, k=-0.0377, n=0.4586, b=-0.00280.95380.04740.0029
Hii, et al.a=0.4984, k=1.46E-07, n=2.531, c=0.4986,
g=1.41E-07
0.98750.02620.0009
Verma, et al.a=0.0611, k=0.0016, g=-0.00160.8260.09260.0103
Diffusiona=-61.1444, k=0.0045, g=-0.00440.94280.05330.0034
Wang and Singha=-0.0003, b=-2.45E-060.98750.06070.0041
Henderson & Pabisa=1.1221, k=0.00160.8260.09260.0097
Logarithmica=90.5617, k=1.49E-05, C=-89.45230.8990.07010.0059
Pelega=2256.81, b=-3.08040.99280.02060.0005

Table 2: Goodness of fit and constant parameter for the selected models for 1.2 m/s.

Air
velocity
TemperatureModelConstantRMSEΧ²
1.4 m/s70 °CLewisk=0.0090.87840.17340.0343
Pagek=1.74E-07, n=3.30750.99820.01980.0005
Midili, et al.a=0.9735, k=3.32E-08, n=3.6575,
b=6.85E-05
0.99860.01520.0005
Hii, et al.a=0.4932, k=2.17E-07, n=3.2543, c=0.4908,
g=2.19E-07
0.9980.01850.0009
Verma, et al.a=0.0825, k=0.0104, g=-0.01040.85970.15780.0398
Diffusiona=-102.2737, k=0.0249, g=-0.02450.96250.0830.011
Wang and Singha=-0.0051, b=-6.65E-070.92670.1150.0176
Henderson & Pabisa=1.165, k=0.01040.85970.15780.0332
Logarithmica=3.8986, k=0.0017, C=-2.79790.9310.10620.018
Pelega=195.0655, b=-0.01780.92690.1150.0176
60 °CLewisk=0.00360.81950.18220.0365
Pagek=4.65E-08, n=3.10890.99130.0390.0019
Midili, et al.a=0.9687, k=-0.0658, n=0.4421, b=-0.00710.96130.0680.0073
Hii, et al.a=0.0087, k=6.5757, n=3.4364, c=0.9926,
g=-8.59E-09
0.99410.02650.0013
Verma, et al.a=0.0942, k=0.0046, g=-0.00460.79120.16090.0356
Diffusiona=-100.889, k=0.0117, g=-0.01150.92610.09720.013
Wang and Singha=-7.76E-06, b=-1.14E-050.96820.06180.0047
Henderson & Pabisa=1.1885, k=0.00460.79120.16090.0316
Logarithmica=110.4495, k=3.15E-05, C=-109.28380.90880.10430.015
Pelega=712.8187, b=-1.38580.94960.0840.0086
50 °CLewisk=0.00330.83610.1830.0363
Pagek=1.02E-07, n=2.89940.99240.03670.0016
Midili, et al.a=1.0352, k=1.02E-06, n=2.4904, b=0.00E000.9810.05360.0041
Hii, et al.a=0.4995, k=3.81E-08, n=3.0746, c=0.4999,
g=3.69E-08
0.99520.03270.0017
Verma, et al.a=0.1001, k=0.0042, g=-0.00420.80630.16050.0335
Diffusiona=-103.7812, k=0.0103, g=-0.01010.93170.09730.0123
Wang and Singha=-0.0003, b=-7.64E-060.97620.05530.0036
Henderson & Pabisa=1.2002, k=0.00420.80630.16050.0305
Logarithmica=142.3656, k=2.16E-05, C=-141.19860.92910.09490.0117
Pelega=764.637, b=-1.18940.96190.07620.0069

Table 3: Goodness of fit and constant parameter for the selected models for 1.4 m/s.

Effective Moisture Diffusivity and Activation Energy

Table 4 shows the drying of India spinach and examined its moisture diffusivity and activation energy. The effective moisture diffusivity (Deff) was measured at velocities of 1.2 m/s and 1.4 m/s, along with temperatures of 70°C, 60°C, and 50°C. Higher velocities led to increased Deff values, indicating faster moisture removal. Activation energy (Ea) was also studied, with values of 85.81 kJ/mol and 43.05 kJ/ mol at velocities of 1.2 m/s and 1.4 m/s, respectively. The decrease in Ea at higher velocities suggests less energy is required for drying. This is supported by the enhanced heat and mass transfer between spinach and air. These findings align with previous work by Fernando and Amarasinghe [23], while Thirupathihalli, et al. [24] reported a similar trend in mango drying, contrasting with Luna, et al. [22] on Piper umbellatum.

VelocityTemperatureabD (x10-12 m2/s)
eff
rσEa (kJ/mol)
1.2 m/s700.5353-2.22E-0422.45-0.91361.324285.81
600.3239-9.07E-059.190.98780.565
500.2061-3.45E-053.490.99420.444
1.4 m/s700.6945-3.61E-0436.53-0.94162.260943.05
600.5072-1.35E-0413.650.85041.6155
500.6592-1.40E-0414.210.97283.3417

Table 4: Effective moisture diffusivity and activation energy. a and b are the intercept and slope of lnMR vs time, r is correlat

Machine Learning Modeling

Moisture Ratio: The fitness statistics for Moisture Ratio (MR) prediction were analyzed using Artificial Neural Network (ANN) and Support Vector Machine (SVM) models during the drying process of Indian spinach as show in Table 5 & 6. Among the ANN configurations, MLP 3-8-1 achieved high accuracy (R=0.9992, RMSE=0.0069) for testing. MLP 3-7-1 and MLP 3-4-1 also showed strong performance. SVM models, including Linear, Polynomial, RBF, and Sigmoid kernels, were evaluated. RBF kernel stood out with R=0.9846, RMSE=0.027, and MSE=0.0035 for overall stages. Polynomial kernel also performed well. Linear kernel had decent results, while Sigmoid kernel showed the weakest performance. Hyper parameters were provided for each kernel. The RBF kernel demonstrated the best correlation, accuracy, and lowest errors. Proper SVM model selection, particularly RBF, enhances MR prediction accuracy in Indian spinach drying [24, 25, 26, 27].

NetworkTrainingValidationTestingHidden
activation
Output
activation
RRMSEMSERRMSEMSERRMSEMSE
MLP 3-8-10.99920.01080.00010.99780.02090.00040.99920.00690ExponentialExponential
MLP 3-8-10.99930.00930.00010.99660.01820.00030.99870.00770.0001TanhLogistic
MLP 3-7-10.99920.010.00010.99660.01830.00030.99880.00710TanhLogistic
MLP 3-4-10.99860.01350.00020.99520.02190.00050.99920.00590TanhLogistic
MLP 3-8-10.99780.01670.00030.99610.02110.00040.99750.01270.0002ExponentialTanh
MLP 3-3-10.99840.0140.00020.99450.02390.00060.99750.00990.0001TanhExponential
MLP 3-8-10.99770.01740.00030.99340.02660.00070.99890.00790.0001TanhLogistic
MLP 3-4-10.99690.01980.00040.99710.0170.00030.99470.01480.0002ExponentialLogistic
MLP 3-4-10.99690.01960.00040.99540.02480.00060.99640.01420.0002TanhTanh
MLP 3-7-10.99720.01880.00040.99610.01980.00040.99320.01960.0004ExponentialExponential
MLP 3-4-10.99620.02220.00050.99680.01870.00040.99170.02110.0004ExponentialExponential
MLP 3-6-10.99550.0240.00060.99280.02620.00070.9950.01460.0002TanhLogistic
MLP 3-7-10.99510.02510.00060.99170.02880.00080.99690.0140.0002TanhTanh
MLP 3-3-10.99650.02160.00050.98920.03740.00140.99720.01030.0001LogisticTanh
MLP 3-8-10.99480.0260.00070.99170.030.00090.99670.01630.0003ExponentialTanh
MLP 3-6-10.99630.02250.00050.99010.03290.00110.99330.0190.0004ExponentialTanh
MLP 3-7-10.99070.03470.00120.99190.02790.00080.99020.02020.0004LogisticTanh
MLP 3-5-10.99480.0260.00070.99140.0320.0010.98750.02220.0005LogisticExponential
MLP 3-5-10.99430.02670.00070.99070.03260.00110.9830.02790.0008LogisticIdentity
MLP 3-5-10.98610.04150.00170.99110.03140.0010.9870.03060.0009ExponentialLogistic

Table 5: The Fitness Statistics of ANN (MLP) model for MR prediction.

StatisticsLinearPolynomialRBFSigmoid
TrainingR0.8720.97760.98880.8415
RMSE0.1720.08210.05230.192
MSE0.02960.00670.00270.0369
TestingR0.73840.95190.94990.5994
RMSE0.19260.09130.07520.208
MSE0.03710.00830.00570.0433
OverallR0.85020.97670.98460.8137
RMSE0.17840.08460.05910.1962
MSE0.03180.00710.00350.0385
HyperparameterC=10.000,
nu=0.100
C=10.000, nu=0.400,
degree=3.000, gamma=0.680,
coefficient=0.000
C=10.000, nu=0.300,
gamma=2.870
C=7.000, nu=0.400,
gamma=0.221,
coefficient=0.000

Table 6: The fitness statistics of SVM model for MR prediction.

Drying Rate: Table 7 shows the multi-Layer Perceptron (MLP) models for India spinach drying prediction, uncovering varying performance. The most effective MLP 3-5-1 model with Exponential hidden and Logistic output activation yielded R values of 0.9398 (validation) and 0.9873 (testing). Similarly, the MLP 3-5-1 with Logistic hidden and Exponential output achieved high scores (R = 0.9388 for validation, 0.9756 for testing). Respectable results were seen with the MLP 3-6-1 model (R = 0.9254 for validation, 0.9720 for testing) and the MLP 3-2-1 (R = 0.9187 for validation, 0.9913 for testing). These models sufficiently predicted drying rate behavior. Table 8 shows the Support Vector Machine (SVM) analysis which indicated the RBF kernel outperformed others. It showed the highest correlation (R = 0.8836) and lowest RMSE (0.1212) in overall predictive accuracy, surpassing Linear, Polynomial, and Sigmoid kernels across training and testing. Hyper parameters for each kernel configuration were also detailed [6].

NetworkTrainingValidationTestingHidden
activation
Output
activation
RRMSEMSERRMSEMSERRMSEMSE
MLP 3-5-10.61470.15010.02250.93980.03980.00160.98730.0270.0007LogisticExponential
MLP 3-5-10.62670.14870.02210.93880.03890.00150.97560.0280.0008ExponentialLogistic
MLP 3-5-10.61420.15020.02260.92180.04470.0020.9910.02040.0004LogisticExponential
MLP 3-6-10.62380.14890.02220.92540.04240.00180.9720.02720.0007ExponentialLogistic
MLP 3-2-10.61280.15070.02270.91870.04630.00210.99130.02730.0007TanhLogistic
MLP 3-10-10.60980.15080.02270.97690.0350.00120.92510.0410.0017TanhLogistic
MLP 3-2-10.6110.15180.0230.91430.05260.00280.99110.02150.0005LogisticLogistic
MLP 3-2-10.6110.15180.0230.91430.05260.00280.99110.02150.0005LogisticLogistic
MLP 3-10-10.61880.15020.02260.90320.04710.00220.96430.03870.0015ExponentialLogistic
MLP 3-10-10.61880.15020.02260.90320.04710.00220.96430.03870.0015ExponentialLogistic
MLP 3-9-10.60790.15140.02290.89850.0480.00230.98020.03470.0012ExponentialLogistic
MLP 3-8-10.57430.15720.02470.93170.03980.00160.91270.04830.0023ExponentialTanh
MLP 3-4-10.58720.15460.02390.87130.05630.00320.97910.03290.0011ExponentialLogistic
MLP 3-7-10.5660.15660.02450.89830.05330.00280.93480.03920.0015TanhIdentity
MLP 3-6-10.54180.15970.02550.93150.04770.00230.87020.05110.0026ExponentialIdentity
MLP 3-2-10.53030.1610.02590.93680.04580.00210.80360.0610.0037LogisticTanh
MLP 3-7-10.5650.15690.02460.88840.05120.00260.85520.05470.003TanhIdentity
MLP 3-5-10.5440.15950.02540.90320.05160.00270.84340.05610.0031LogisticTanh
MLP 3-10-10.57060.15620.02440.8870.05770.00330.86620.05690.0032TanhTanh
MLP 3-9-10.53510.1610.02590.87690.05920.00350.91550.04220.0018LogisticIdentity

Table 7: The fitness statistics of ANN (MLP) model for DR prediction.

StatisticsLinearPolynomialRBFSigmoid
TrainingR0.32090.34280.85750.344
RMSE0.24180.23250.1340.2307
MSE0.05840.05410.0180.0532
TestingR0.69950.91260.94990.5708
RMSE0.21590.19790.0840.202
MSE0.04660.03920.00710.0408
OverallR0.43890.4530.88360.4116
RMSE0.23430.22270.12120.2225
MSE0.05490.04960.01470.0495
HyperparameterC=10.000,
nu=0.100
C=7.000, nu=0.500,
degree=2.000, gamma=0.215,
coefficient=0.000
C=7.000,
epsilon=0.100,
gamma=9.300
C=9.000, nu=0.500,
gamma=0.222,
coefficient=0.000

Table 8: The fitness statistics of SVM model for DR prediction.

Conclusion

From this study, it can be inferred that higher temperatures speed up Indian spinach drying by reducing moisture content and drying time. Increased air velocities also enhance moisture evaporation. Efficient convective heat and mass transfer were recorded. The Page and Hii, et al. models accurately represent drying behavior. Higher air velocities increase effective moisture diffusivity, quickening moisture removal and lowering activation energy. A 3-8-1 MLP setup with Tanh activation in hidden layers and Logistic activation in the output is optimal for Moisture Ratio prediction. In SVM models, the RBF kernel consistently outperforms others in correlation, RMSE, and MSE for drying prediction. Polynomial kernel is viable, Linear and Sigmoid kernels perform less well. SVM kernel-specific hyper parameters give insights. MLP models using Exponential hidden and Logistic output activations excel in predicting Indian spinach drying.

References

  1. Hanumappa M (2019) Malabar spinach (Basella alba) is a Nutritious and Ornamental Plant. University of the District of Columbia, USA.
  2. Garavand AT, Rafiee S, Keyhani A (2011) Study on Effective moisture diffusivity, activation energy and mathematical modeling of thin layer drying kinetics of bell pepper. Australian Journal of Crop Science AJCS 5(2): 128-131.
  3. Touil A, Chemkhi Saber, Zagrouba F (2014) Moisture Diffusivity and Shrinkage of Fruit and Cladode of Opuntia ficus-indica during Infrared Drying. Journal of Food Processing 2014: 1-9.
  4. Mujumdar AS, Devahastin S (2008) Fundamental principles of drying in Guide to Industrial Drying.
  5. Deshmukh SA, Gaikwad DK (2014) A review of the taxonomy, ethnobotany, phytochemistry diffusivity, activation energy and mathematical modeling of thin layer drying kinetics of bell pepper. Australian Journal of Crop Science 5(2): 128-131.
  6. Caywood MS, Roberts DM, Colombe JB, Greenwald HS, Weiland MZ (2017) Gaussian process regression for predictive but interpretable machine learning models: An example of predicting mental workload across tasks. Frontiers in Human Neuroscience 10: 647.
  7. Midilli A, Kucuk H (2003) Mathematical modeling of thin layer drying of pistachio by using solar energy. Energy Conversion Management 44(7): 1111-1122.
  8. Kaleemullah S, Kailappan R (2006) Modelling of thin- layer drying kinetics of red chilies. Journal of Food Engineering 76(4): 531-537.
  9. Wang CY, Singh RP (1978) A single layer drying equation for rough rice. American society Agricultural Engineering, Paper No. 78-3001, USA.
  10. Onwude DI, Hashim N, Abdan K, Janius R, Chen G (2018) Modelling the mid-infrared drying of sweet potato: Kinetics, mass and heat transfer parameters, and energy consumption. Heat Mass Transfer 54: 2917-2933.
  11. Kashaninejad M, Mortazavi A, Safekordi A, Tabil LG (2007) Thin-layer drying characteristics and modelling of pistachio nuts. Journal of Food Engineering, 78(1): 98- 108.
  12. Da Silva WP, eSilva CMDPS, Gama FJA, Gomes JP (2014) Mathematical models to describe thin-layer drying and to determine drying rate of whole bananas. Journal of the Saudi Society of Agricultural Sciences 13(1): 67-74.
  13. Derossi A, Severini C, Cassi D (2011) Mass Transfer Mechanisms during Dehydration of Vegetable Food: Traditional and Innovative Approach. In: El-Amin M (Ed.), Advanced Topics in Mass Transfer. Intech, pp: 305- 354.
  14. Hii CL, Law CL, Cloke M, Suzannah M (2009) Thin layer drying kinetics of cocoa and dried product quality. Biosystems Engineering 102(2): 153-161.
  15. Yaldiz O, Ertekin C (2001) Thin layer solar drying of some vegetables. Drying Technology 19(3-4): 583-597.
  16. Premi M, Sharma H, Upadhyay A (2012) Effect of Air Velocity and Temperature on the Drying Kinetics of Drumstick Leaves (Moringa Oleifera). International Journal of Food Engineering 8(4): 1-21.
  17. Jun-Wen B, Hong-Wei X, Hai-Le M, Cun-Shan Z (2018) Artificial Neural Network Modeling of Drying Kinetics and Color Changes of Ginkgo Biloba Seeds during Microwave Drying Process. Journal of Food Quality 2018: 1-8.
  18. Rodríguez J, Clemente G, Sanjuán N, Bon J (2014) Modelling drying kinetics of thyme (Thymus vulgaris L.): Theoretical and empirical models, and neural networks. Food Science and Technology International 20(1): 13- 22.
  19. Ningning J, Jiyang MA, Rongge MA, Yang Z, Panyu C, et al. (2023) Effect of slice thickness and hot-air temperature on the kinetics of hot-air drying of Crabapple slices. Food Science and Technology, pp: 1-8.
  20. Krokida MK, Karathanos VT, Maroulis ZB, Marinos- Kouris D (2003) Drying kinetics of some vegetables. Journal of Food Engineering 59(4): 391-403.
  21. Clemente G, Sanjuán N, Andres Carcel J, Mulet A (2014) Influence of Temperature, Air Velocity, and Ultrasound Application on Drying Kinetics of Grape Seeds. Drying Technology 32(1): 68-76.
  22. Dorneles LNS, Goneli ALD, Cardoso CAL, da Silva CB, Hauth MR, et al. (2019) Effect of air temperature and velocity on drying kinetics and essential oil composition of Piper umbellatum L. leaves. Industrial Crops & Products 142, pp: 111846.
  23. Fernando JAKM, Amarasinghe ADUS (2016) Drying kinetics and mathematical modeling of hot air drying of coconut coir pith. SpringerPlus 5: 807.
  24. Murthy TPK, Manohar B (2013) Hot air-drying characteristics of mango ginger: Prediction of drying kinetics by mathematical modeling and artificial neural network. J Food Sci Technol 51(12): 3712-3721.
  25. Qadri OS, Osama K, Srivastava AK (2020) Foam mat drying of papaya using microwaves: Machine learning modeling. Journal of Food Process Engineering 43(6).
  26. Lertworasirikul S, Saetan S (2010) Artificial neural network modeling of mass transfer during osmotic dehydration of kaffir lime peel. Journal of Food Engineering 98(2): 214-223.
  27. Chaurasia P, Younis K, Qadri OS, Srivastava G, Osama K (2019) Comparison of Gaussian process regression, artificial neural network, and response surface methodology modeling approaches for predicting drying time of mosambi (Citrus limetta) peel. Journal of Food Process Engineering 42(2).

Cite this article

BibTeX
APA
RIS
@article{jegede2023,
  title   = {Mathematical Modelling of the Drying Characteristics of Indian Spinach
Using Machine Learning Approach},
  author  = {Jegede AO* and Isa J},
  journal = {Food Science & Nutrition Technology},
  year    = {2023},
  volume  = {8},
  number  = {3},
  doi     = {10.23880/fsnt-16000308}
}
Jegede AO* and Isa J (2023). Mathematical Modelling of the Drying Characteristics of Indian Spinach
Using Machine Learning Approach. Food Science & Nutrition Technology, 8(3). https://doi.org/10.23880/fsnt-16000308
TY  - JOUR
TI  - Mathematical Modelling of the Drying Characteristics of Indian Spinach
Using Machine Learning Approach
AU  - Jegede AO* and Isa J
JO  - Food Science & Nutrition Technology
PY  - 2023
VL  - 8
IS  - 3
DO  - 10.23880/fsnt-16000308
ER  -