Optimizing Uniaxial Oil Extraction of Bulk Rapeseeds: Spectrophotometric and Chemical Analyses of the Extracted Oil under Pretreatment Temperatures and Heating Intervals

: Optimizing the operating factors in edible oil extraction requires a statistical technique such as a response surface methodology for evaluating their effects on the responses. The examined input factors in this study were the diameter of pressing vessel, V D (60, 80, and 100 mm), temperature, T PR (40, 60, and 80 ◦ C), and heating time, H TM (30, 60 and 90 min). The combination of these factors generated 17 experimental runs where the mass of oil, oil yield, oil extraction efﬁciency, and deformation energy were calculated. Based on the response surface regression analysis, the combination of the optimized factors was V D : 100 (+1) mm; T PR : 80 ◦ C (+1) and H TM : 60 (0) min); V D : 60 (–1) mm; T PR : 80 ◦ C (+1) and H TM : 75 (+0.5) min and V D : 100 (+1) mm; T PR : 80 ◦ C (+1) and H TM : 90 (+1). The absorbance and transmittance values signiﬁcantly ( p < 0.05) correlated with the wavelength and temperature, but they did not correlate signiﬁcantly ( p > 0.05) with heating time. The peroxide value did not correlate signiﬁcantly with temperature, however, it correlated signiﬁcantly with heating time. Neither the acid value nor the free fatty acid value correlated with both temperature and heating time. The ﬁndings of the present study are part of our continuing research on oilseeds’ processing optimization parameters. Author Contributions: Conceptualization, C.D., A.K., D.H. and ˇC.M.; funding acquisition, D.H.; methodology, C.D., A.K., ˇC.M. and P.H.; validation, C.D.; A.K. and O.D.; formal analysis, C.D., A.K. and O.D.; data curation, C.D., A.K. and P.H.; writing—original


Introduction
Oilseed rape (Brassica napus L.) is the second-highest potential source of vegetable oil with a high nutritional value and a favorable composition of fatty acids for both food and animal feed [1,2]. Oilseed rape, also known as winter oilseed crop, is Europe's prime oilseed crop, widely grown in Germany, Poland, the Czech Republic, and France [3][4][5][6][7]. In general, oilseeds provide many nutritious and functional properties for human health, such as starch, crude protein content, oil content, fatty acids, amino acids, vitamins, phytosterols, and polyphenols [8][9][10]. Oilseed crops, such as rape, sunflower, safflower, canola, mustard, and camelina, among others, are noteworthy feedstocks for biodiesel production, being an alternative renewable energy source for reducing greenhouse gas emissions caused by fossil fuels [11]. Worldwide, the major oilseed crops grown are soya, rape, cotton, pea, sunflower, oil palm, and copra [12,13].

Determination of Moisture and Oil Content
The moisture content and percentage oil content of the rapeseeds sample were determined based on the conventional oven drying and Soxhlet extraction methods [36][37][38][39]. The hot air oven produced by MEMMERT GmbH + Co. KG, Buechenbach, Germany, was used to dry the rapeseeds sample at a temperature setting of 105 • C and a drying time of 24 h. With the Soxhlet extraction procedure, 10 g of the rapeseeds sample was ground in a mini grinder. The ground sample was loaded into a thimble and cotton wool was placed atop. The thimble was inserted into the Soxhlet extractor of 250 mL of solvent volume, which then was connected to a 500 mL round bottom flask containing 250 mL of petroleum ether. The complete setup was placed under a heating source at a temperature of 60 • C, where the solvent was heated to reflux for 24 h. The extracted oil was left in the oven without drying or heating for 3 days to remove the residual solvent. Measurements were done in triplicates. The electronic balance Kern 440-35 (Kern and Sohn GmbH, Balingen, Germany), with an accuracy of 0.001 g was used for weighing the extracted oil samples. Based on the relation given by [40], the moisture content and percentage oil content were calculated to be 6.37 ± 0.24 (% w.b.) and 41.35 ± 0.70 (%), respectively.

Box-Behnken Experimental Design of Compression Factors
When several processing factors and their interactions are likely to influence the output parameters, the response surface methodology with Box-Behnken Design (BBD) is used to find the optimum processing factors [41][42][43]. The BBD based on the input factors generated 17 experimental runs with five repetitions at the center point. The independent factors were pressing vessel diameter, heating temperature, and heating time, with each factor set at three levels. The mathematical equation defining the Box-Behnken design is given in equation (Equation (1)) as follows: where Y is the response variable; i and j are linear and quadratic coefficients; β 0 , β i , β ii , and β ij are the regression coefficients in the intercept, linear, quadratic, and interaction terms, respectively; X i and X j are the independent variables; and k is the number of factors. The factors were coded from −1 to +1 (low, center, and high) using equation (Equation (2)) [44,45] as follows: where x i is the coded value of the ith variable, X i is the uncoded value of the ith test variable, and X 0 is the uncoded value of the ith test variable at the center point.

Pretreatment of Rapeseeds Sample
Before the compression tests, the measured sample of rapeseeds was pretreated or conditioned at temperatures of 40, 60, and 80 • C and heating times of 30, 60, and 90 min using the conventional oven method (MEMMERT GmbH + Co. KG, Buechenbach, Germany).

Compression Tests of Bulk Rapeseeds Sample
The sample of rapeseeds was measured at a constant compression height of 100 mm in each of the compression chambers/vessels of diameters 60, 80, and 100 mm using their plungers. The weights of the sample were 190.33 g, 338.80 g, and 524.54 g. Based on the cross-sectional area of the pressing vessels, the volume of the sample in each vessel was calculated to be 28.27 × 10 −5 , 50.27 × 10 −5 and 78.54 × 10 −5 m 3 , respectively. Preliminary experiments were performed to determine the maximum compression force (without the serration effect or the ejection of the seedcake through the pressing holes) for each of the pressing vessels' diameters of 60, 80, and 100 mm. Maximum forces of 180, 300, and 450 kN at a compression speed of 4 mm/min were determined in an increasing order of the pressing vessels' diameters. Based on the observed allowable processing conditions stated above, each of the 17 factor combinations generated from the Box-Behnken Design were then tested using the universal compression testing equipment ( Figure 1). The data obtained were used to calculate the responses/parameters: mass of oil, oil yield, oil extraction efficiency, and deformation energy. The mass of oil was calculated gravimetrically (as the difference of mass of seed cake and initial mass of the sample M S (g)). The oil yield was calculated based on the relations reported by [46,47] as given in equation (Equation (3)).
where O Y is the oil yield (%) and M O is mass of oil (g). The oil extraction efficiency was calculated according to the relations mentioned by [46,47], as given in equation (Equation (4)).
where O S is the percentage of oil content (%) determined from the Soxhlet extraction method. The deformation energy was calculated according to the relations stated by [48][49][50], as given in equation (Equation (5)).
where E N is the deformation energy (kJ), F n+1 + F n and x n+1 − x n are the compression force (kN) and deformation (mm), n is the number of data points, and i is the number of sections in which the axis deformation was divided. The hardness, H D (kN/mm) was calculated as the ratio of maximum compression force M F (kN) [48][49][50] to that of deformation D F (mm) as given in equation (Equation (6)).
where V M is the volume of rapeseeds (m 3 ), D is the diameter of pressing vessel (×10 −6 m 2 ), and H is the pressing height of the rapeseeds sample (×10 −3 m). Using the universal

Spectrophotometric Analysis of Extracted Oil
Spectrophotometric analysis was carried out using the UV-VIS spectrophotometer (VIS V-10 Plus, Giorgio Bormac S.r.l., Carpi, Italy) to determine the absorption and transmission rates of the extracted rapeseed oil within a specified wavelength (325 and 600 nm). This information has been reported to be functional for assessing the quality of the oil for the prevention of UV radiation on human skin [51].

Chemical Analysis of Extracted Oil under Pretreatment Temperatures
The rapeseed oil extracted at pretreatment temperatures between 40 and 80 °C was analyzed in terms of peroxide value (PV), acid value (AV), and free fatty acid (FFA). The rapeseed oil at a laboratory temperature of 22 °C served as the control. The procedures

Spectrophotometric Analysis of Extracted Oil
Spectrophotometric analysis was carried out using the UV-VIS spectrophotometer (VIS V-10 Plus, Giorgio Bormac S.r.l., Carpi, Italy) to determine the absorption and transmission rates of the extracted rapeseed oil within a specified wavelength (325 and 600 nm). This information has been reported to be functional for assessing the quality of the oil for the prevention of UV radiation on human skin [51].

Chemical Analysis of Extracted Oil under Pretreatment Temperatures
The rapeseed oil extracted at pretreatment temperatures between 40 and 80 • C was analyzed in terms of peroxide value (PV), acid value (AV), and free fatty acid (FFA). The rapeseed oil at a laboratory temperature of 22 • C served as the control. The procedures reported by [8,52,53] were followed. For PV, 5 g of the oil sample was weighed into a volumetric flask, then 30 mL of chloroform and a glacial acetic acid mixture of ratio (2:3) were added for the dissolution. The mixture was shaken vigorously for exactly 1 min, followed by the addition of 30 mL of distilled water. The mixture was titrated with 0.1 M sodium thiosulphate solution until the yellow color disappeared using 1 ml of 1% starch as an indicator. Peroxide value was expressed as (meq O 2 /kg). For the determination of AV and FFA, 5 g each of the oil sample was weighed into a volumetric flask and then 100 ml of neutralized ethanol (warmed up to 60-65 • C) was added together with a 2 ml of 1% phenolphthalein and immediately titrated with ethanolic KOH (0.1 Normality) to obtain an appearance that was light pink in color. AV and FFA were expressed as (mg KOH/g oil).
The measurements were made twice, and the results averaged.

Statistical Analysis of Calculated Responses
The experimental data were statistically evaluated using STATISTICA 13 [54] by applying the following statistical techniques: basic statistics (correlation analysis) and general linear models (repeated measured ANOVA, post hoc tests, simple regression, multiple regression and response surface regression) at 5% significance level.

Determination of Maximum Compression Force
In this present study, two compression tests were conducted. The first test was the control experiment to determine the maximum compression forces for the different diameters of the pressing vessel at an initial sample pressing height of 100 mm and speed of 4 mm/min. The second test was the compression factors combination (diameter of pressing vessel, pretreatment temperature, and heating time) based on the Box-Behnken experimental design. From the control experiment, the mass of oil (g), oil yield (%), oil extraction efficiency (%), deformation energy (kJ), and hardness (kN/mm) were calculated, as presented in (Table 1). The mass of oil, deformation energy, and hardness increased with the increase in vessel diameter, whereas the deformation, oil yield, and oil extraction efficiency decreased with the increase in vessel diameter. The deformation energy is the area under the force-deformation curve ( Figure 2) [27,49,50]. The ANOVA results (Table 2) showed that the compression factors (vessel diameters and forces) had a significant effect (p < 0.05) on the calculated responses, except for deformation, which was not significant (p > 0.05). The coefficient of determination (R 2 ) values confirming the results were between 0.288 and 0.997. The results of the Box-Behnken design are explained in the succeeding sections. Here, it is important to highlight that three main responses: oil yield (%), oil extraction efficiency (%), and deformation energy (kJ) were examined in relation to the combination of the compression factors.

Oil Yield, Oil Extraction Efficiency, and Deformation Energy
Three compression factors with three levels each, namely vessel diameter (60, 80, and 100 mm), heating temperature (40, 60, and 80 °C), and heating time (30,60, and 90 min) were evaluated for the compression tests of bulk rapeseeds sample. Based on the Box-Behnken design of experiments (BBD) coupled with the response surface methodology (RSM), 17 experimental runs were obtained with the 12 factors combination and 5 repetitions at the center points (Table 3). For the different vessel diameters at each maximum compression force at a constant speed of 4 mm/min, the calculated responses were oil yield, oil extraction efficiency, and deformation energy. Oil yield and oil extraction efficiency values ranged from 16.172 to 24.783% and 39.109 to 59.934%, respectively. The deformation energy values ranged from 1.17 to 3.19 kJ. From the BBD experimental data (Table 3), the compression factors combination of vessel diameter (60 (−1), temperature (80 (+1), and heating time (60 (0)) produced the highest oil yield of 24.783%, oil extraction efficiency of 59.934%, and deformation energy of 1.255 kJ. The optimum factors for rapeseeds oil extraction in terms of oil yield, oil extraction efficiency, and deformation energy were determined based on the response surface regression analysis (Section 3.4).

Oil Yield, Oil Extraction Efficiency, and Deformation Energy
Three compression factors with three levels each, namely vessel diameter (60, 80, and 100 mm), heating temperature (40, 60, and 80 • C), and heating time (30,60, and 90 min) were evaluated for the compression tests of bulk rapeseeds sample. Based on the Box-Behnken design of experiments (BBD) coupled with the response surface methodology (RSM), 17 experimental runs were obtained with the 12 factors combination and 5 repetitions at the center points (Table 3). For the different vessel diameters at each maximum compression force at a constant speed of 4 mm/min, the calculated responses were oil yield, oil extraction efficiency, and deformation energy. Oil yield and oil extraction efficiency values ranged from 16.172 to 24.783% and 39.109 to 59.934%, respectively. The deformation energy values ranged from 1.17 to 3.19 kJ. From the BBD experimental data (Table 3), the compression factors combination of vessel diameter (60 (−1), temperature (80 (+1), and heating time (60 (0)) produced the highest oil yield of 24.783%, oil extraction efficiency of 59.934%, and deformation energy of 1.255 kJ. The optimum factors for rapeseeds oil extraction in terms of oil yield, oil extraction efficiency, and deformation energy were determined based on the response surface regression analysis (Section 3.4). Table 3. Box-Behnken design of compression factors combination, coded values, and calculated responses (oil yield, oil expression efficiency, and deformation energy).

Force-Deformation Curves of Experimental Runs (BBD)
The 17 experimental runs (R1 to R17) ( Table 3) of the factor levels combination in terms of the force-deformation curves are graphically illustrated in Figure 3. The maximum force for each vessel diameter of 60, 80, and 100 mm were determined from the initial pressing height of the bulk rapeseeds sample measured at 100 mm, which was compressed at a speed of 4 mm/min. A higher force with a bigger vessel diameter produced the maximum oil output. The curves showed a smooth pattern without the serration effect based on the control experiments (Section 3.1).

Response Surface Regression Analysis of Factor Combinations
The effect of the factor combinations on the calculated parameters or responses was statistically analyzed based on the response surface regression statistical technique. The

Response Surface Regression Analysis of Factor Combinations
The effect of the factor combinations on the calculated parameters or responses was statistically analyzed based on the response surface regression statistical technique. The results are presented in Table 4. For the mass of oil (g), the coefficients of the linear and quadratic terms, as well as the interaction terms of the vessel diameter and temperature, were significant (p < 0.05). However, the interaction terms of vessel diameter and heating time and that of temperature and heating time were not significant (p > 0.05). For oil yield (%) and oil expression efficiency (%), the coefficients of the quadratic term of the vessel diameter and interaction terms of the factors (vessel diameter, temperature, and heating time) were not significant (p > 0.05) in comparison with the coefficients of the other terms of the factors, which were significant (p < 0.05). For deformation energy, the coefficients of the linear and quadratic terms of the vessel diameter, the linear term of temperature, and the linear term of heating temperature were significant (p < 0.05), whereas the coefficients of the quadratic terms of the temperature and heating time, as well as the interaction terms of the factors, were not significant (p > 0.05). The intercept coefficients of all the models were significant (p < 0.05). The significance of the results explains the accuracy of the models for predicting the calculated responses.

Determined Regression Models for Predicting the Responses
The oil yield (%), oil extraction efficiency (%), and deformation energy (kJ) were the main responses from the compression tests of bulk rapeseeds based on the compression factor combinations, which were generated from the BBD. The linear regression models defining these responses as a function of the compression factors/predictors are expressed in equations (Equation (8)) to (Equation (10)) respectively. The intercepts and the coefficients of the predicators and their interactions were statistically significant (p < 0.05) for predicting the calculated responses (Table 4).

Determined Optimum, Predicted, and Validated Values of the Responses
The optimum, predicted, and validated values of the responses (oil yield (%), oil expression efficiency (%), and deformation energy (kJ)) are given in Table 5. Based on the response surface regression analysis [54], the optimum values of the responses in relation to the compression factor combinations: (V D : 60 (-1) mm; T PR : 80 • C (+1) and H TM : 75 (+0.5) min) and (V D : 100 (+1) mm; T PR : 80 • C (+1) and H TM : 90 (+1) min were obtained from the surface profile plots (Figure 4a-c). The predicted values were obtained from the linear regression models using equations (Equations (8)-(10)) which were validated based on additional experiments. The desirability values of the optimal responses and their factors ranged between 0.979 and 1 and the coefficient of variation and percentage error between the predicted and validated values ranged from 0.09 to 2.30%, which confirms the reliability of linear regression models (Equations (8)-(10)) for predicting the responses. The surface and area contour plots of the interaction effect of the compression factors (temperature, pretreatment, and heating time) on the responses (oil extraction efficiency and deformation energy) are illustrated in Figure 5. In Figure 5a, at a constant heating time, the increase in the diameter of pressing vessel from 60 to 100 mm (coded as −1 to +1) showed no increases in oil extraction efficiency but the increase in temperature from 40 to 80 • C (coded as −1 to +1) recorded 59%, while their combined effect decreased to 56%. In Figure 5b, at a constant temperature, the increase in the diameter of the pressing vessel neither increased nor decreased the extraction efficiency, but the increase in heating time from 30 to 90 min (coded as −1 to +1) slightly increased it to 53%, while their interaction effect decreased it to 50%. In Figure 5c, at a constant pressing vessel diameter, the increase in temperature increased the oil extraction efficiency to approximately 55% and the increase in heating time did not considerably increase the oil extraction efficiency. However, their interaction effect increased it to 57%. On the other hand, in Figure 5d, at a constant heating time, the increase in the diameter of the pressing vessel increased the deformation energy by about 2.9 kJ, while the increase in temperature did not increase the deformation energy. However, their combined effect produced 3.2 kJ of deformation energy. Furthermore, in Figure 5e, at a constant temperature, the increase in the diameter of the pressing vessel increased the deformation energy to 3 kJ, while the increase in heating time showed no increase in deformation energy but their interaction effect gave approximately 3.2 kJ. Finally, in Figure  5f, at constant pressing vessel diameter, the increase in temperature and heating time, as well as their interaction, had no significant effect on the deformation energy.

Compression and Relaxation Processes of Rapeseeds Oil Extraction
The uniaxial compression process is the dependency between compression force and deformation, which can be followed by the relaxation process, which describes the dependency between relaxation force and time at a constant strain of the bulk seeds. The relaxation process is performed to recover the residual oil in the seedcake immediately after the compression process. The compression factor combinations that produced higher oil yield (%) and oil extraction efficiency (%) were subjected to a relaxation process at a constant time of 20 min. The results in comparison with the compression process are given in Table 6

Compression and Relaxation Processes of Rapeseeds Oil Extraction
The uniaxial compression process is the dependency between compression force and deformation, which can be followed by the relaxation process, which describes the dependency between relaxation force and time at a constant strain of the bulk seeds. The relaxation process is performed to recover the residual oil in the seedcake immediately after the compression process. The compression factor combinations that produced higher oil yield (%) and oil extraction efficiency (%) were subjected to a relaxation process at a constant time of 20 min. The results in comparison with the compression process are given in Table 6 Figure 6. The thick line of the relaxation process is due to the changes in pressure of the piston (hydraulic transmission system). There was no change in pressure at the compression process, hence the thin line. The results are further explained in Section 4. Table 6. Data of the relaxation process of bulk rapeseed oil extraction with optimal compression factor combinations.

Chemical Properties of Rapeseed Oil at Pretreatment Temperatures
The chemical properties (peroxide value, acid value, and free fatty acid) of t rapeseed oil extracted at different pretreatment temperatures and heating intervals a given in Table 7

Chemical Properties of Rapeseed Oil at Pretreatment Temperatures
The chemical properties (peroxide value, acid value, and free fatty acid) of the rapeseed oil extracted at different pretreatment temperatures and heating intervals are given in Table 7. The averaged peroxide values at heating times (30, 60, and 90 min) for pretreatment temperatures (40, 60, and 80 • C) were 5.10 ± 0.77, 6.19 ± 1.61, and 5.71 ± 1.32 meq O 2 /kg. Similarly, the acid values were 1.43 ± 0.39, 1.42 ± 0.29, and 1.49 ± 0.30 mg KOH/g oil. The free fatty acid values were 0.72 ± 0.20, 0.71 ± 0.15, and 0.75 ± 0.15 mg KOH/g oil. The averaged peroxide values increased from 40 • C to 60 • C and then decreased at 80 • C with heating times between 30 and 90 min. Acid values and free fatty acid values decreased from 40 • C to 60 • C and then increased at 80 • C with heating times. In terms of varying heating times at constant temperature, peroxide values at 60 • C and 80 • C increased with an increase in heating times, whereas at 40 • C it decreased from 30 min to 60 min and then increased at 90 min. On the other hand, acid values and free fatty acid values increased from 30 min to 60 min and then decreased at 90 min at 40 • C and 60 • C, respectively. However, at 80 • C they decreased from 30 min to 60 min and then increased at 90 min. The data were further subjected to normality tests and ANOVA tests of between-subject effects to assess the significant effect of the compression factors on the calculated chemical properties. Based on the Shapiro-Wilk test (Table 8), the data showed a normal distribution function. The normal distribution function was assessed by the fact that the p-values were greater than the significance level of 5% with the corresponding high values of the coefficient of determination (R 2 ) ranging from 0.720 to 0.974 [55]. The ANOVA tests of between-subject effects on the chemical properties of the rapeseed oil (Table 9) revealed that peroxide value with temperatures and heating times was significant (p < 0.05) compared to the acid value and free fatty acid value, which showed non-significance (p > 0.05) with temperature, but were significant (p < 0.05) with heating times. The compression factor interaction effect on peroxide value indicated non-significance, whereas that of acid value and free acid value proved significant. Further statistical explanation (correlation, multiple regression, and post hoc tests) is provided in the Supplemental Materials (Section 3.10) and Discussion (Section 4).

Effect of Compression Factors on Absorbance and Transmittance of Rapeseed Oil
Spectrophotometric analysis of the extracted rapeseed oil at different pretreatment temperatures and heating times was evaluated at different wavelengths to understand the effect of the pretreatment conditions on the spectral properties. The data were subjected to various statistical analyses. The ANOVA analysis (Table 10) showed that the individual compression factors and their interactions had a significant effect (p < 0.05) on the spectral profiles (absorbance and transmittance versus wavelength). The corrected model of the spectral profiles produced high values of the coefficient of determination (R 2 ) of 0.999 and 0.995 respectively confirming the significant effect of the compression factors on the absorbance and transmittance values. The spectral profiles at different temperatures and heating times are graphically shown in Figures 7 and 8, respectively. Additional statistical evaluation (correlation, multivariate tests of significance, multiple regression, and normality tests) of the experimental data are provided in the Supplemental Materials (Section 3.10), explaining in detail the significance of the results.

Supplementary Materials
The data on chemical properties (peroxide value, acid value, and free fatty acid) were further subjected to correlation analysis (Table S1). The correlation analysis showed that peroxide value did not significantly correlate with temperature (P > 0.05). However, it did significantly correlate (P < 0.05) with heating time. The acid value and free fatty acid value also did not significantly correlate (P > 0.05) with both temperature and heating time. The

Supplementary Materials
The data on chemical properties (peroxide value, acid value, and free fatty acid) were further subjected to correlation analysis (Table S1). The correlation analysis showed that peroxide value did not significantly correlate with temperature (p > 0.05). However, it did significantly correlate (p < 0.05) with heating time. The acid value and free fatty acid value also did not significantly correlate (p > 0.05) with both temperature and heating time. The regression models for predicting the above-mentioned parameters and their coefficients of determination (R 2 ) are given in Tables S2 and S3. The coefficients of determination (R 2 ) for the regression models were found to be 0.629 and 0.039 (Table S3), respectively. The peroxide value of the rapeseed oil was significantly (p > 0.05) influenced by the heating temperatures and times. The significance of the result was further subjected to post hoc tests (Tukey HSD and Duncan) based on multiple comparisons of the mean differences of the compression factors. For both the Tukey HSD and Duncan tests, the mean difference of the temperatures of the sample of 40 • C and 80 • C showed non-significant results (p > 0.05), similar to the mean difference of the temperatures of the sample of 80 • C and 60 • C. However, the sample temperatures at 40 • C and 60 • C showed significance (p < 0.05) from their mean difference. On the other hand, the mean difference of the sample heating intervals of 30 min and 60 min with the Tukey's HSD test indicated non-significance, whereas 30 min and 90 min, as well as 60 min and 90 min, tested significantly. In comparison with the Duncan test, the mean differences between 30 and 60 min, 30 and 90 min, and 60 and 90 min were significant (Table S4). The experimental data of the spectral profiles (absorbance and transmittance versus wavelength) were tested for normality. Based on the Kolmogorov-Smirnov normality test, the absorbance values at wavelengths from 325 nm to 335 nm, at 405 nm, from 415 nm to 425 nm, at 450 nm, and from 465 nm to 600 nm were significant (p > 0.05). The transmittance values were also observed to be significant at wavelengths from 325 nm to 335 nm, from 345 nm to 360 nm, at 405 nm, and from 500 nm to 600 nm. The Shapiro-Wilk normality test, on the other hand, showed that the absorbance values at wavelengths from 325 nm to 335 nm and from 465 nm to 600 nm were significant (p < 0.05).
The transmittance values were also significant (p < 0.05) at wavelengths from 325 nm to 335 nm, at 355 nm, 360 nm, and from 500 nm to 600 nm. The Shapiro-Wilk test of normality showed high values of the coefficient of determination (R 2 ) (Tables S5 and S6) compared to the Kolmogorov-Smirnov normality test. The normal distribution of the data was assessed based on the p-values being greater than the significance level of 5% [55]. From the correlation analysis, the absorbance and transmittance values significantly correlated (p < 0.05) with wavelength and temperature but they did not correlate significantly (p > 0.05) with heating time (Table S7). The multivariate tests of significance agreed with the correlation results (Table S8). Finally, the multiple regression analysis of the effects of the factors on the absorbance and transmittance of the rapeseed oil is given in Table S9. The intercept and the coefficients of the wavelength and temperature were significant (p < 0.05) for predicting the absorbance and transmittance of the rapeseed oil at a specific heating time. The spectral profiles of the rapeseed oil under varying temperatures and heating times are illustrated in Figure S1.

Discussion
In developing countries, mechanical screw presses are suitable for oil extraction compared to advanced methods such as ultrasound-assisted extraction [12,14,15,24]. Understanding fully the uniaxial compression process under a laboratory scale is key for optimizing the mechanical screw pressing system for large scale oil production for both domestic and industrial applications.
In this present study, the applied maximum compression forces of 180, 300, and 450 kN were equivalent to the pressure values of 63.66, 59.68, and 57.3 MPa, which were calculated based on the cross-sectional area of the compression vessels/chambers of diameters 60, 80, and 100 mm respectively. The examined compression factors, among others such as screw geometry/configuration and its components, namely nozzle sizes, screws with choke worm shaft ring sizes, and press cylinders with mesh sizes, thus affect the oil output, residual oil in seedcake, input energy, and oil quality under both the uniaxial compression and mechanical oil extraction processes [56][57][58][59].
These above-mentioned compression factors and responses need to be optimized using appropriate mathematical and statistical tools. In this context, the response surface methodology (RSM) based on the Box-Behnken design (BBD) was employed to optimize the compression factors-vessel diameters, pretreatment temperatures, and heating intervals-for the responses-oil yield, oil extraction efficiency, and deformation energyfor extracting rapeseed oil under uniaxial compression. The BBD generated 17 experimental runs with 12 factor combinations and 5 repetitions at the center points. Based on the BBD experimental data (Table 3), the compression factor combinations (Run 3) with pressing vessel diameter (V D : 60 (−1)) mm, temperature (T PR : 80 • C (+1)), and heating time (H TM : 60 (0) min) produced the highest oil yield of 24.783% and oil extraction efficiency of 59.934%. The reason could be that the smaller pressing vessel of diameter 60 mm provided a much smaller space, which helped to increase the force (pressure) towards the seeds, resulting in high oil yield compared to the bigger pressing vessels of diameters 80 and 100 mm, which provided larger space and, hence, less pressure towards the seeds producing low oil yield [47]. In addition, an increase in pretreatment temperature and heating time will increase oil yield [60]. However, higher levels of heating will reduce the moisture content of the seeds, thereby impeding the cell wall of the seeds to break/crack, thus generating a low percentage of oil yield and/or oil extraction efficiency [47,60]. The experimental data were further analyzed using the response surface regression function in Statistica 13 software [54] to optimize the compression factors with the responses that were described by linear regression models ((Equations (8)- (10)). The linear regression models indicated a good fit for prediction, with high coefficient of determination values between 0.975 and 0.998. Again, the models were adequate since the F-values were greater than the p-values and the lack-of-fit p-values were non-significant (p > 0.05). [33,41,61] mentioned that a good regression model (one containing squared terms, products of two factors, linear terms, and intercepts) must be significant, the lack-of-fit p-value of the model must also be insignificant, and the coefficient of determination (R 2 ) value of the model should be closer to 1. The results of the present study agreed with the above statements or boundary conditions. Furthermore, applying the RSM approach requires that not only the optimum compression factors of the responses be determined, but also, they need to be validated through additional experiments. The optimum compression factor combinations for predicting the responses (oil yield (%), oil extraction efficiency (%) and deformation energy (kJ)) were observed as (V D : 60 (−1) mm; T PR : 80 • C (+1) and H TM : 75 (+0.5) min) and (V D : 100 (+1) mm; T PR : 80 • C (+1) and H TM : 90 (+1) min) with corresponding desirability values between 0.979 and 1. Ref. [34] explained that the desirability function can be used to determine optimum performance of the responses concerning the independent factors. The values of the coefficient of variation and percentage error of the predicated responses and their validation were between 0.09 and 2.3%, confirming the validity of the regression models ((Equations (8)-(10)).
In the industrial oil extraction process, the press cake with the residual oil content between 20% and 25% is mostly subjected to solvent extraction using n-hexane to recover the residual oil [8,20,22,35]. Under the uniaxial compression process, the residual oil can be recovered by the relaxation process, which is the dependency between the relaxation force and time at a constant strain of the bulk oilseeds. The relaxation process can be set together with the compression process. Residual oil of approximately 6% was recovered with a relaxation time of 10 min. This means that a total of 25 min was used for both the compression (15 min) and relaxation processes to extract the maximum oil yield and/or to achieve higher oil extraction efficiency. This information is useful for designing new oil extraction systems, such as mechanical screw presses, to avoid the combined use of the mechanical screw presses and the solvent extraction method to achieve high oil extraction efficiency or to minimize the residual oil in the seedcake. It is important to note that the examined compression factors (diameter of pressing vessel, temperature, and heating time) influenced the relaxation process in terms of the percentage oil yield and/or oil extraction efficiency. Lower levels of the compression factor combinations recorded lower amounts of the residual oil in the seedcake, whereas higher levels produced higher amounts of the residual oil. Besides, since the relaxation process is done at constant strain, there is no deformation energy utilization compared to the compression process where the deformation energy is characterized by the area under the force-deformation curve [27,50]. For the experimental design factor combinations (V D : 60 (−1)) mm; temperature (T PR : 80 • C (+1)) and heating time (H TM : 60 (0) min) the deformation energy of 1.255 kJ produced the highest oil yield and/or oil extraction efficiency, as stated above. In comparison with the optimal factor combinations (V D : 60 (−1) mm; T PR : 80 • C (+1) and H TM : 75 (+0.5) min), the deformation energy was 1.24 ± 0.01 kJ. The difference was almost negligible. The recorded deformation energy values indicate the threshold energy required for obtaining the maximum oil output concerning the compression factor combinations already stated above. However, the compression factor combinations (V D : 100 (+1) mm; T PR : 80 • C (+1) and H TM : 90 (+1) min) produced higher deformation energy between 3.19 and 3.27 kJ, with a lower oil yield of 22.972 ± 0.519% and oil extraction efficiency of 55.553 ± 1.255%, indicating that the compression factor combinations (V D : 60 (−1) mm; T PR : 80 • C (+1) and H TM : 60(0)/75 (+0.5) min) is the most optimal for estimating the oil yield (%), oil extraction efficiency (%), and deformation energy (kJ) of rapeseed oil extraction under uniaxial compression process.
Regarding the chemical properties of the extracted oil, the pretreatment temperatures did not significantly affect (p > 0.05) the peroxide value, acid value, and free fatty acid of the extracted rapeseed oil. The heating time at a specific pretreatment temperature had a significant effect on the peroxide value compared to the acid value and free fatty acid, which had no significant effect on the heating time. Ref. [62], cited in [24], indicated that excesses in the heating time and pretreatment temperature will result in seedcake with lower nutritional value and lower oil quality. Hence, the quality of the extracted oil can be preserved at a lower conditioning temperature of 60 • C and heating times between 50 and 60 min [24,63]. Generally, the peroxide value (PV), acid value (AV), and free fatty acid (FFA) are among the important physicochemical indicators for the quality assessment of edible oils. PV, AV, and FFA are evidence of autoxidation (free radical reaction) and hydrolytic rancidity [64][65][66]. Mathematically, FFA is half of AV [67]. The oxidation and chemical changes in oils during heating are characterized by an increase in free fatty acid content and a decrease in the total unsaturation of oils [67,68]. A high PV value may reflect either increased formation of hydroperoxides or reduced decomposition, whereas a high AV or FFA content results in increased losses during refining, poor flavor quality and stability of the finished edible oil, and rancidity of the oil [69][70][71][72]. PV, AV, and FFA values have been reported to range from 1.9 to 31.2 meq O 2 /kg oil, 0.6 to 4 mg KOH/g oil, and 1.122 to 10.261 mg KOH/g oil, respectively, which are influenced by different processing conditions and varieties of oilseeds [67,73]. In the present study, the low values of the PV, AV, and FFA achieved at the observed processing conditions meet the acceptable quality standards of edible oils [64][65][66][67][68][69][70][71][72][73].
Finally, the spectral profiles (absorbance and transmittance versus wavelength) of the extracted oil at different temperatures and heating times were described. The pretreatment temperatures and heating times increased the absorbance and transmittance values from 0.4 to 3.0 (-) and 22% to 45% with the wavelength between 325 and 600 nm. The absorbance increase occurred at wavelengths between 375 and 450 nm, whereas the transmittance increase was observed at wavelengths between 500 and 600 nm. However, at wavelengths between 525 and 600 nm, the absorbance values decreased. Similarly, the transmittance values decreased at wavelengths between 350 and 500 nm. It is important to mention that by examining the transmittance and absorbance versus wavelength of the extracted rapeseed oils, it could be also possible to analyze their oxidation stability as influenced by temperature and time pretreatment conditions [74,75]. In addition, the high absorption and low transmission rates of the extracted rapeseed oil suggest that the extracted oils at the various pretreatment conditions can be used for the prevention of ultraviolet radiation problems on human skin, as reported by [51]. However, this information needs to be studied extensively using advanced spectroscopic techniques [74][75][76][77][78].

Conclusions
The universal compression testing machine of a load capacity of 500 kN was used to evaluate the effect of compression factors, namely the diameter of the pressing vessel (V D : 60, 80, and 100 mm), pretreatment temperature (T PR : 40, 60, and 80 • C), and heating time (H TM : 30, 60, and 90 min) on the responses of oil yield (%), oil extraction efficiency (%), and deformation energy (kJ). The compression factors and the responses were optimized by applying the response surface methodology (RSM) coupled with the Box-Behnken design (BBD), which is a statistical tool for analyzing the effect of independent factors on the dependent parameters. The Box-Behnken design generated 17 experimental runs, involving 12 compression factor combinations and 5 repetitions at the center points. From the BBD experimental data, the compression factor combinations and their coded values of vessel diameter (60 (−1), temperature (60 (0), and heating time (60 (0)) produced the maximum oil yield of 24.783% and oil extraction efficiency of 59.934%, with the corresponding deformation energy of 1.255 kJ. Based on the response surface regression analysis, linear regression models ((Equations (8)-(10))) were described for predicting the responses at optimum compression factor combinations (V D : 60 (−1) mm; T PR : 80 • C (+1) and H TM : 75 (+0.5) min). The statistical lack-of-fit p-values of the regression models were non-significant (p > 0.05) and the coefficient of variation and percentage error values between the predicted and validated responses ranged between 0.08% and 2.30%, which indicate that the determined regression models were adequate for predicting the responses. The relaxation time of 10 min with the compression factor combinations (V D : 80 (0) mm, T PR : 80 • C (+1), and H TM : 90 (+1) min) recovered the maximum amount of the residual oil of approximately 6% from the seedcake. On contrary, the compression factor combinations (V D : 60 (−1) mm, T PR : 80 • C (+1), and H TM : 60 (0)/75(+0.5) min) may not require the relaxation process since the residual oil recovered was negligible.
The chemical properties (peroxide value, acid value, and free fatty acid value) and spectral properties (absorbance and transmittance versus wavelength) of the extracted rapeseed oil were not significantly affected by the studied compression factors, indicating that the rapeseed oil could be extracted at pretreatment temperatures between 40 and 80 • C and heating times between 30 and 90 min without any quality assessment problems related to domestic, industrial, and pharmaceutical applications.
The present study is part of our continuing research on the uniaxial compression process of different oilseeds to fully understand the compression factors concerning the oil extraction efficiency and deformation energy requirement to reduce the residual oil in the seedcake using the mechanical screw presses, which will be the focus of future studies.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/10 .3390/pr9101755/s1, Table S1: Correlation analysis of the chemical properties with the factors effect; Table S2: Multiple regression results of the chemical properties with the factors effect; Table S3: Test of the sum of squares whole model of the factors effect on the peroxide value, acid value, and free fatty acid value of rapeseed oil; Table S4: Post hoc tests of peroxide value of rapeseed oil with the factors effect; Table S5: Tests of normality of absorbance, A (-), of rapeseed oil with the factors effect; Table S6: Tests of normality of transmittance, T (%), of rapeseed oil with the factors effect; Table S7: Correlation analysis of the absorbance and transmittance of rapeseed oil with the factors effect; Table S8: Multivariate tests of significance of the factors effects on absorbance (-) and transmittance, T (%); Table S9: Multiple regression analysis of the factors effect on absorbance, A (-), and transmittance, T (%); Figure S1: Absorbance and transmittance versus wavelength of rapeseed oil at laboratory temperature, heating temperatures, and heating times.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.