Index

Abstract

The conventional solvent extraction (CSE) method is commonly used to extract oil from oilseeds. Even though the extraction of oil from soybean seeds using the CSE method has been studied extensively, the application of response surface methodology (RSM) has not yet been widely explored. In this paper, the CSE method was employed using three different types of solvents – ethanol (Eth), n-hexane (Hex), and a mixture of Hex:Eth in a ratio of 2:1 – and optimized using RSM to extract the oil from Cambodian soybean seeds. For each solvent, a central composite design (CCD) was applied to optimize the extraction conditions, such as the extraction times (15, 30, and 45 min) and solvent-to-solid ratios (5:1, 10:1, and 15:1 mL/g). This CCD generated 11 experimental runs for each type of solvent. The optimum oil yields when using ethanol, Hex:Eth, and n-hexane were 20.53%, 18.78%, and 17.39%, respectively. The optimum condition was obtained at 15 min extraction time and 15:1 (mL/g) solvent-to-solid ratio for all solvent types. The coefficients of determination of ethanol, Hex:Eth, and n-hexane were 0.9710, 0.9954, and 0.9263, respectively, which mathematically indicated a good model for the prediction. Therefore, the model was considered accurate and reliable for predicting oil yields in this study.

Keywords: Cambodian soybeans, Central composite design, Ethanol, Extraction time, Hexane, Mixture of solvents, Model reduction, Oil yield, Optimization, Solvent-to-solid ratio.

Received: 31 October 2022 / Revised: 7 December 2022 / Accepted: 19 December 2022 / Published: 9 January 2023

Contribution/ Originality

Cambodian soybean oil has not yet been introduced to the market even though Cambodia is an agricultural country. This is the first report on the use of RSM to optimize the interactive effect between solvent-to-solid ratio and extraction time on the Cambodian soybean oil extraction process using different types of solvents.

1. INTRODUCTION

Cambodia is an agricultural country where soybeans are grown for local consumption. In Cambodia, 92% of soybean growers are located in the two provinces of Ratanakiri (52%) and Preah Vihear (40%), with the remainder in the provinces of Kratie (4%), Stung Treng (3%), and Mondulkiri (1%) (Nget et al., 2021). However, Cambodian soybean oil has not yet been introduced to the market. Soybean oil is a vegetable oil extracted from soybean seeds. It is generally obtained through mechanical pressing or solvent extraction methods (Lima et al., 2008). Among these methods, solvent extraction is the method commonly used to extract oil from seed materials with lower oil content, such as soybean seeds (Keneni, Bahiru, & Marchetti, 2021). The conventional solvent extraction (CSE) method is the process by which a molecule transfers from one solvent to another due to the difference in solubility between two immiscible (or slightly soluble) solvents (Chen & Wang, 2017). The main advantages of CSE are that it provides a higher oil yield, larger processing capacity, and lower refinery losses. In the CSE method, n-hexane is the solvent commonly used in the oil processing industry because it is inexpensive and has high solubility. Previous studies (Mohammed & AI-Saddi, 2003; Nikolic, Stankovic, Cakic, & Mitic, 2009; Shittu, Mari, & Dangora, 2019) have explored the extraction of soybean oil using n-hexane. Their results proved that n-hexane is the solvent that provides the highest soybean oil extraction yield. However, the use of n-hexane as an extraction solvent has many drawbacks due to its non-renewable fossil origin and high flammability, leading to environmental and public health concerns (Nde & Foncha, 2020; Oliveira, Garavazo, & Rodrigues, 2012; Russin, Boye, Arcand, & Rajamohamed, 2011; Tabtabaei & Diosady, 2013) . Normally, ethanol can be used instead of n-hexane to extract oil from soybean seeds. Certain studies (Oliveira et al., 2012; Rodrigues, Aracava, & Abreu, 2010; Sawada, Venâncio, Toda, & Rodrigues, 2014; Toda, Sawada, & Rodrigues, 2016) have shown that ethanol is a solvent that can serve as a substitute for n-hexane in oil production as it has a proven ability to extract oil from the seeds and offers environmental benefits. However, although many researchers have studied the extraction of soybean oil using the CSE method, response surface methodology (RSM) has not yet been widely explored. RSM is one of the most used tools for optimizing the influence process variables in experimental designs (Khuri, 2017; Kumari & Gupta, 2019). RSM mainly identifies the response by using given process variables following experimental designs, such as central composite design (CCD), that fit an empirical full second-order polynomial model. In previous studies (Aydar, 2018; Bokhari, Yusup, & Ahmad, 2012; Nde & Foncha, 2020; Popoola, Akinoso, & Raji, 2016; Yolmeh & Jafari, 2017) , CCD has been used for the optimization of process parameters during oil extraction and provided successfully optimized conditions. Therefore, the current study aims to determine the effect of different solvents on soybean oil yield and to optimize the extraction conditions, including solvent-solid ratio and extraction time, using RSM.

2. METHODOLOGY

2.1. Chemicals and Reagents

N-hexane of 99% purity (RCI Labscane, Thailand) and ethanol of above 99.9% purity (DAEJUNG, South Korea) were used to extract the oil from Cambodian soybean seeds.

2.2. Sample Preparation

Soybean seeds were collected from Ratanakiri province, Cambodia. The soybean seeds were cracked and then manually dehulled to remove the hulls and impurities from the seeds. After that, the dehulled soybeans were ground and sieved to obtain soybean powder with a particle size lower than 0.425 mm. The moisture content of the dried soybean seeds was 8.16 ± 0.3% as measured by a moisture analyzer (MOC63u, Shimadzu).

2.3. Experimental Design

RSM with CCD was used for the analysis and optimization of the interaction effects between variables for Cambodian soybean oil extraction using the CSE method. A CCD was applied to two independent variables with three coded levels (-1, 0, +1) and two un-coded levels (-a, +a), in which a was 1.4142. The independent variables were extraction time (A), ranging from 15 to 45 (min), and solvent-to-solid ratio (B), ranging from 5:1 to 15:1 (mL/g). The three types of solvents used in this study were ethanol, n-hexane, and a mixture of n-hexane and ethanol (Hex:Eth) in a ratio of 2:1. The independent variables with their coded and uncoded values are shown in Table 1.

Table 1. Coded and uncoded factor levels generated from the CCD for each solvent.
Variables
Levels
-1.4142
-1
0
1
+1.4142
A: Extraction time (Min)
8.7868
15
30
45
51.2131
B: Solvent-to-solid ratio (mL/g)
2.9289
5:1
10:1
15:1
17.071

This CCD comprised three center points with a one-star point, and the response for CCD optimization was soybean oil yield. The CCD was based on the second-order polynomial model as in the following equation:

2.4. Oil Extraction Procedure

In this study, soybean oil was extracted using the CSE method with three different solvent types: ethanol, n-hexane, and a mixture of n-hexane and ethanol (Hex:Eth) in a ratio of 2:1. The soybean powder was added to the solvent at the defined solvent-to-solid ratios (5:1 to 15:1 mL/g). During the extraction, the mixture of sample and solvent was heated at 60 ± 2°C and stirred at 500 rpm for the defined extraction times (15 to 45 min). After the extraction, the mixtures were centrifuged at 4,000 rpm for 15 min. After centrifugation, the liquid extracts were filtered with 25 µm filter paper. Next, the miscella was evaporated at 40°C, 72 rpm, and 30 min. The collected oil was dried at 105°C for 30 min in a hot-air drying oven and then transferred into a desiccator for 15 min.

2.5. Determination of Soybean Oil Yield

The extracted oil was calculated and expressed as a percentage, which was defined as the weight of the extracted oil over the wet weight of the sample taken:

Where %Yield is the percentage by weight of extraction yield in wet basis (% wet weight basis), mf is the weight of round bottom flask containing distilled soybean oil (g), mi is the tare weight of round bottom flask (g), and m is the weight of the sample used for extraction (g).

2.6. Statistical Analysis

In this study, Design-Expert software version 13.05.0 (State-Ease Inc., USA) was used for experimental design and statistical analysis. Analysis of variance (ANOVA) was used to estimate the suitability and validity of the model with an above 95% confidence level. The optimal condition was then chosen based on the desirability (numerical optimization) of the response variable.

3. RESULTS AND DISCUSSION

3.1. Soybean Oil Yield Using Different Types of Solvent

The results of the experiments in terms of un-coded values with independent variables and responses are depicted in Table 2. The variation in the oil values indicated that the extraction factors affected the extracted oil yields.

3.2. Effect of Solvent Types on Extracted Oil Yield

The extraction solvent is one of the most significant factors affecting the extraction efficiency of oil from oilseeds. Oil can be extracted from plant seeds using different solvents (Bardone et al., 2012; Shivani, Khushbu, Faldu, Thakkar, & Shubramanian, 2011). The extraction of soybean oil was performed using the CSE method and three different extraction solvents (ethanol, n-hexane, and Hex:Eth). The results showed that ethanol was the best solvent, providing the highest crude oil yield, followed by the Hex:Eth mixture and, finally, n-hexane. This result could be partially explained by the difference in solvent polarities (Li et al., 2014). Ethanol is polar, while n-hexane is a non-polar solvent. Therefore, ethanol can extract both non-polar and polar lipids, such as triglycerides, phospholipids, tocopherols, waxes, and sugars, while n-hexane can extract only non-polar lipids. This explains the higher yield obtained when using ethanol as an extraction solvent.

Table 2. Results of soybean oil yields extracted using various types of solvents.
Std.
Runs
Independent variables
Oil yields (%)
A
B
Ethanol
N-hexane
Hex:Eth (2:1)
7
1
30
2.929:1
11.65
14.81
15.01
4
2
45
15:1
21.74
17.63
18.69
6
3
51.213
10:1
19.8
17.19
17.97
9
4
30
10:1
20.09
16.8
17.29
1
5
15
5:1
13.52
13.8
15.38
8
6
30
17.071:1
20.14
18.01
19.77
3
7
15
15:1
20.47
17.41
18.93
11
8
30
10:1
19.12
16.64
18.11
5
9
8.787
10:1
18.7
15.36
16.64
10
10
30
10:1
18.85
16.88
17.01
2
11
45
5:1
16.41
16.2
15.40
Note: A: Extraction time (Min), B: Solvent-to-solid ratio (mL/g).

3.3. Regression Model for Soybean Oil Yield

The linear terms, interaction term, and quadratic terms were employed in multiple regression analysis to obtain predicted responses. The second-order polynomial equations in terms of coded and actual factors for the various types of solvent were presented as follows:

For ethanol:

The coded equation was used to determine the relative effect of the factors on the response by comparing the coefficients of factors. Therefore, the influence of the factors on the oil yields extracted using ethanol, n-hexane, and the Hex:Eth mixture were B > B2 > A > AB > B2, B > A > AB > A2 > B2, and B > A > A2 > B2 > AB, respectively. Additionally, the equation in terms of actual factors can be used to predict the response for given levels of each factor. The given levels should be specified in the original units for each factor. Moreover, the sign of the regression coefficients can be used to show how the terms influenced oil yield. A positive coefficient sign indicates that the term's effect favors the extracted oil yield, while a negative sign indicates that the term's effect will be antagonistic to the extracted oil yield. The results of the statistical analysis of the regression models are indicated in Table 3, with F-test and p-values calculated using ANOVA. Various statistical parameters, such as the probability value (p-value), the regression coefficient (R2 value), the Fisher model value (F-value), and the lack of fit F value, were used to modify the responses to fit the mathematical model.

The result of ANOVA for soybean oil extraction using ethanol, n-hexane, and Hex:Eth is thus indicated in Table 3. The F-values of the models Ycoded(1), Ycoded(2), and Ycoded(3) were 44.47, 215.4, and 17.7, respectively. Therefore, there was only a 0.04%, 0.01%, and 0.34% chance that a strong F-value of the models Ycoded(1), Ycoded(2), and Ycoded(3) could occur due to noise. Moreover, the p-values of the regression model of the oil extracted using ethanol (p = 0.0004), n-hexane (p < 0.0001), and Hex:Eth (p = 0.0034) were all lower than 0.05, indicating that the model was significant. In addition, the p-values of lack of fit of oil extraction using ethanol (p = 0.5332), n-hexane (p = 0.5229), and Hex:Eth (p = 0.6376) were all higher than the significant level.

Table 3. Multiple regression analysis results for the significant variables and interactions on the response of soybean oil yield using various types of solvents.
SM
SS
DF
MS
F-value
P-value
Remark
Ethanol (Yethanol)
Model
94.96
5
18.99
44.47
0.0004
Significant
A
4.08
1
4.08
9.56
0.0271
 
B
73.73
1
73.73
172.63
< 0.0001
 
AB
0.6561
1
0.6561
1.54
0.2702
 
0.0231
1
0.0231
0.0541
0.8253
 
14.7
1
14.7
34.42
0.002
 
Residual
2.14
5
0.4271
 
Lack of fit
1.28
3
0.4283
1.01
0.5332
Not significant
Pure error
0.8505
2
0.4252
Cor total
97.1
10
Std. dev.
0.6535
R² = 0.978
Mean
18.23
Adjusted R² = 0.956
C.V. %
3.59
Predicted R² = 0.886
Press
11.05
Adeq precision = 19.301
n-hexane (Yn-hexane)
Model
16.52
5
3.3
215.4
< 0.0001
Significant
A
3.39
1
3.39
221.03
< 0.0001
B
11.44
1
11.44
745.63
< 0.0001
AB
1.19
1
1.19
77.46
0.0003
0.411
1
0.411
26.8
0.0035
0.2311
1
0.2311
15.07
0.0116
Residual
0.0767
5
0.0153
Lack of fit
0.0468
3
0.0156
1.05
0.5229
Not significant
Pure error
0.0299
2
0.0149
 
Cor total
16.6
10
 
Std. dev.
0.1239
R² = 0.9954
 
Mean
16.43
Adjusted R² = 0.9908
 
C.V. %
0.7538
Predicted R² = 0.9759
 
Press
0.4002
Adeq precision = 45.3716
 
Hex:Eth (YHex:Eth)
Model
23.52
5
4.7
17.7
0.0034
Significant
A
0.3448
1
0.3448
1.3
0.3064
 
B
23.02
1
23.02
86.6
0.0002
 
AB
0.0169
1
0.0169
0.0636
0.811
 
0.1177
1
0.1177
0.4427
0.5353
 
0.0586
1
0.0586
0.2204
0.6585
 
Residual
1.33
5
0.2659
 
Lack of fit
0.6757
3
0.2252
0.6892
0.6376
Not significant
Pure error
0.6536
2
0.3268
Cor total
24.85
10
Std. dev.
0.5156
R² = 0.9465
Mean
17.29
Adjusted R² = 0.893
C.V. %
2.98
Predicted R² 0.7475
Press
6.28
Adeq precision = 12.6003
Note: A: Extraction time (Min), B: Solvent-to-solid ratio (mL/g), DF: Degree of freedom, SM: Source of model, MS: Mean of squares, SS: Sum of squares, Adeq precision: Adequate precision.

Thus, the lack of fit was insignificant, and the insignificant lack of fit was good. The fitness of the polynomial model was expressed by the coefficient of determination R², coefficient of adjusted R² and predicted R². It has been suggested that these values should be at least 0.80 for a good fit of the model, and the difference between adjusted R² and predicted R² should be less than 0.2 (Yujie, 2011). Therefore, the value of R², adjusted R², and predicted R² of the model indicated that this regression model was a good fit for the mathematical model. The signal-to-noise ratio was measured with adequate precision. It is preferable to have a ratio of more than 4. The adequate precisions of the regression model using ethanol, n-hexane, and Hex:Eth were 19.30, 45.37, and 12.60, respectively. Therefore, this model represented a good sign and had a strong enough signal to be used for optimization.

3.4. Effect of Extraction Time on Soybean Oil Yield

The extraction time is an important parameter of the oil extraction process (Panchal, Deshmukh, & Sharma, 2014). The effect of extraction time on soybean oil yield varied depending on the type of solvent used. Based on the results of the statistical analysis in Table 3, the extraction time had a highly significant effect on oil yield when n-hexane was used to extract the oil (p < 0.0001), a significant effect for ethanol (p = 0.0271), and an insignificant effect for the Hex:Eth mixture (p = 0.3396). Figure 1 shows the effect of extraction time and solvent-to-solid ratio on oil yield extracted using the various types of solvents. It was observed that when n-hexane and ethanol were used in the extraction, the oil yield increased slowly due to the increase in the extraction time. However, when the mixture of Hex:Eth was used, the oil yield did not change in response to an increase in extraction time. The effect of extraction time on oil yield can be explained by the reduction of residual oil in the sample. At the beginning of the extraction process, the solvent starts to dissolve the oil in the materials. The higher (total) oil content in the seed materials in contact with the solvent results in the higher oil molecules being diffused into the solvent. Thus, the extraction time should be long enough to let the oil molecules diffuse into the solvent (Abed, Kurji, & Abdul-Majeed, 2015). However, when the equilibrium time was reached, the oil yield remained stable in response to an increase in extraction time (Shao et al., 2012). This is due to the low solvent density and oil content remaining in the materials after the extraction time has been extended (Saxena, Sharma, & Sambi, 2011).

3.5. Effect of Solvent-to-Solid Ratio on Soybean Oil Yield

The solvent-to-solid ratio is also an important parameter of the oil extraction process (Lawson, Oyewumi, Ologunagba, & Ojomo, 2010; Panchal et al., 2014). According to the statistical analysis, as shown in Table 3, the solvent-to-solid ratio had a strong significant effect on the oil yield extracted from the soybean seeds. The effects of the solvent-to-solid ratio on the oil yield extracted using ethanol, n-hexane, and the Hex:Eth mixture are shown in Figure 1. It was observed that an increase in solvent-to-solid ratio led to a significant increase in the extraction yield of soybean oil. The increase in oil yield due to the increase in the solvent-to-solid ratio can be explained by the fact that at a lower solvent-to-solid ratio, the amount of solvent in contact with the sample decreases. Then the oil is less soluble in the solvent (Jisieike & Betiku, 2020). However, at a higher solvent-to-solid ratio, the amount of solvent contacting the sample increases, which can improve the rate of mass transfer, resulting in higher extraction productivity (Bokhari et al., 2012). Additionally, the increased solvent-to-solid ratio results in a faster diffusion rate of oil from the extracted materials to the solvent (Predescu et al., 2016).

Figure 1. 3D response surface plots showing the effect of solid-to-solvent ratio and extraction time on the oil yield using different types of solvents: (a). Ethanol, (b). N-hexane, and (c). Mixture of Hex:Eth

3.6. Model Reduction

Model reduction is the elimination of the not significant terms from the model, such as the term for a predictor variable or the interaction between predictor variables. Model reduction helps to increase the precision of the predictions by enhancing the fit between data and model. The regression model of the soybean oil extraction using n-hexane was not reduced because all terms of the model were significant. The reduced regression models of the soybean oil extraction in terms of coded and actual factors are presented in equations as follows:

For ethanol:

Yactual (4) = + 5.38351 + 0.047630A + 1.91305B – 0.065294B2
Ycoded (4) = + 19.41 + 0.7145A + 3.04B – 1.63B2
For the mixture of (Hex:Eth):
Yactual (5) = + 13.89799 + 0.339292B
Ycoded (5) = + 17.29 + 1.70B

The indicators used to assess the fit and quality of the reduced models are presented in Table 4.

Table 4. Indicators to assess the fit and quality of reduced models.
Indicators
Reduced Yethanol
Reduced YHex:Eth
R2
0.9710
0.9263
Adjusted R2
0.9586
0.9182
P-value of the model
<0.0001
<0.0001

The results showed that the reduced models had good accuracy as the R2 values were close to 1 and the values of the adjusted R2 were high enough to suggest a satisfactory correlation between the values predicted by the model and the values observed in the experiments. Additionally, the p-values of the reduced models generated by ANOVA on the reduced regression models were less than 0.0001, indicating that the models obtained were acceptable.

3.7. Optimization of Extraction Conditions Using RSM with CCD

Evaluation of the optimum oil yield was based on numerical optimization using Design-Expert software. The optimization was achieved by setting the optimization goals as follows: minimize the extraction time, maximize the solvent-to-solid ratio, and maximize the oil yield, and the importance of all goals was 3 pluses (+++). The optimal formulation was chosen based on the desirability of the response variable (Pal & Gauri, 2018). The resulting optimum oil yields of soybean oil extraction using ethanol, n-hexane, and the Hex:Eth mixture were 20.10%, 17.39%, and 18.99%, with a desirability of 0.943, 0.930, and 0.942, respectively. The optimum oil yield of all solvents used was obtained at 15 min of extraction time and a 15:1 (mL/g) solvent-to-solid ratio.

3.8. Model Validation

To verify and confirm the agreement between the optimum results obtained from the model and the actual values, three new confirmation experiments were carried out under the suggested optimal conditions, using the same extraction procedure, and the results are presented in Table 5 (Choudhary & Pramanik, 2021).

Table 5. Comparison of predicted and experimental oil yield values from soybean seeds extracted by CSE using different types of solvents.
Solvent
A: Extraction time (Min)
B: Solvent-to-solid ratio (mL/g)
Predicted value
Actual
value
RSE (%)
Ethanol
15
15:1
20.10
20.52 ± 0.65
2.039
N-hexane
15
15:1
17.39
17.23 ± 0.12
0.925
Hex:Eth
15
15:1
18.99
18.02 ± 0.52
5.093
Noted: The data are expressed as mean ± standard deviation for experimental value, RSE: Residue standard error.

The average oil yields that were obtained from the extraction using ethanol, n-hexane, and the mixture were 20.52 ± 0.65%, 17.23 ± 0.12%, and 18.02 ± 0.52%, respectively, while the predicted oil yields were 20.10%, 17.39%, and 18.99%, respectively. The differences in maximum oil yield between the predicted values and the experimental values using ethanol, n-hexane, and the Hex:Eth mixture were 2.039%, 0.925%, and 5.093%, respectively. These errors were within an acceptable range of 10% and confirmed that the predicted model was sufficiently accurate within the 99% prediction interval. Therefore, the model based on the CCD was considered to be accurate and reliable for predicting oil yield extracted from soybean seeds.

4. CONCLUSION

The application of RSM with CCD in this study provided a good fit for the model based on the predicted and experimental data. Among the three solvent types, the highest crude soybean oil yield was obtained using ethanol, followed by the Hex:Eth mixture and then n-hexane. According to the statistical analysis, the solvent-to-solid ratio had a highly significant effect on the extraction process for all solvent types. Similarly, the extraction time had a highly significant effect on the extraction when n-hexane was used. However, when ethanol and Hex:Eth were used in the extraction process, the extraction time had a slightly significant and insignificant effect on the extraction process, respectively. Based on the numerical optimization, the optimum oil yields when using ethanol, the mixture, and n-hexane were 20.53%, 18.78%, and 17.39%, respectively. The optimum condition was obtained at 15 min extraction time and 15:1 (mL/g) solvent-to-solid ratio for all solvent types. In future research, the quality of oils extracted from Cambodian soybean seeds using n-hexane, ethanol, and the mixture of n-hexane and ethanol should be determined, and the refining process should be optimized.

Funding: This research is supported by Cambodia Higher Education Improvement Project (Grant number: 6221-KH).

Competing Interests: The authors declare that they have no competing interests.

Authors’ Contributions: All authors contributed equally to the conception and design of the study.

REFERENCES

Abed, K. M., Kurji, B. M., & Abdul-Majeed, B. A. (2015). Extraction and modelling of oil from Eucalyptus camadulensis by organic solvent. Journal of Materials Science and Chemical Engineering, 3(8), 35–42. https://doi.org/10.4236/msce.2015.38006

Aydar, A. Y. (2018). Utilization of response surface methodology in optimization of extraction of plant materials. In: Statistical Approaches With Emphasis on Design of Experiments Applied to Chemical Processes. In (pp. 158–169). London: IntechOpen.

Bardone, E., Brucato, A., Keshavarz, T., Straccia, M. C., Siano, F., Coppola, R., . . . Volpe, M. G. (2012). Extraction and characterization of vegetable oils from cherry seed by different extraction processes. Chemical Engineering Transactions, 27, 391–396. https://doi.org/10.3303/CET1227066

Bokhari, A., Yusup, S., & Ahmad, M. M. (2012). Optimization of the parameters that affects the solvent extraction of crude rubber seed oil using response surface methodology. Recent Advances in Engineering March 2015, 23–33. https://doi.org/10.1002/er.7126

Chen, H., & Wang, L. (2017). Posttreatment strategies for biomass conversion. Technologies for Biochemical Conversion of Biomass, 197-217. https://doi.org/10.1016/b978-0-12-802417-1.00008-9

Choudhary, A. K., & Pramanik, H. (2021). Optimization and validation of process parameters via RSM for minimizing use of resources to generate electricity from a DEFC. International Journal of Energy Research, 45(14), 20413-20429. https://doi.org/10.1002/er.7126

Jisieike, C. F., & Betiku, E. (2020). Rubber seed oil extraction: Effects of solvent polarity extraction time and solid-solvent ratio on its yield and quality. Biocatalysis and Agricultural Biotechnology, 24, 101522. https://doi.org/10.1016/j.bcab.2020.101522

Keneni, Y. G., Bahiru, L. A., & Marchetti, J. M. (2021). Effects of different extraction solvents on oil extracted from jatropha seeds and the potential of seed residues as a heat provider. Bioenergy Research, 14(4), 1207–1222. https://doi.org/10.1007/s12155-020-10217-5

Khuri, A. I. (2017). Response surface methodology and its applications in agricultural and food sciences. Biometrics & Biostatistics International Journal, 5(5), 155–163. https://doi.org/10.1007/s11947-016-1855-2

Kumari, M., & Gupta, S. K. (2019). Response surface methodological (RSM) approach for optimizing the removal of trihalomethanes (THMs) and its precursor’s by surfactant modified magnetic nanoadsorbents (sMNP)-An endeavor to diminish probable cancer risk. Scientific Reports, 9(1), 1-11. https://doi.org/10.1038/s41598-019-54902-8

Lawson, O. S., Oyewumi, A., Ologunagba, F., & Ojomo, A. (2010). Evaluation of the parameters affecting the solvent extraction of soybean oil. Journal of Engineering and Applied Sciences, 5(10), 51-55.

Li, Y., Fine, F., Fabiano-Tixier, A. S., Abert-Vian, M., Carre, P., Pages, X., & Chemat, F. (2014). Evaluation of alternative solvents for improvement of oil extraction from rapeseeds. Comptes Rendus Chimie, 17(3), 242–251. https://doi.org/10.1016/j.crci.2013.09.002

Lima, M. S. A., Rocha, L. A., Molina, E. F., Caetano, B. L., Marçal, L., Mello, C., . . . Nassar, E. J. (2008). Thermoanalysis of soybean oil extracted by two methods. Quimica Nova, 31(3), 527–529. https://doi.org/10.1590/S0100-40422008000300012

Mohammed, A. H. A.-K., & AI-Saddi, D. A. (2003). Factors affecting the extraction process. Chemical and Petroleum Engineering, 4(2), 1–13.

Nde, D. B., & Foncha, A. C. (2020). Optimization methods for the extraction of vegetable oils. A Review Processes, 8(2), 1-21. https://doi.org/10.3390/pr8020209

Nget, R., Aguilar, E. A., Cruz, P. C. S., Reaño, C. E., Sanchez, P. B., Reyes, M. R., & Prasad, P. V. (2021). Overview of farmers’ perceptions of current status and constraints to soybean production in ratanakiri province of Cambodia. Sustainability, 13(8), 1-16. https://doi.org/10.3390/su13084433

Nikolic, N., Stankovic, M. Z., Cakic, S. M., & Mitic, S. (2009). Effect of extraction techniques on yield and composition of soybean oil. Macedonian Journal of Chemistry and Chemical Engineering, 28(2), 173–179. https://doi.org/10.20450/mjcce.2009.208

Oliveira, C. M., Garavazo, B. R., & Rodrigues, C. E. (2012). Liquid–liquid equilibria for systems composed of rice bran oil and alcohol-rich solvents: Application to extraction and deacidification of oil. Journal of Food Engineering, 110(3), 418-427. https://doi.org/10.1016/j.jfoodeng.2011.12.027

Pal, S., & Gauri, S. K. (2018). A desirability functions-based approach for simultaneous optimization of quantitative and ordinal response variables in industrial processes. International Journal of Engineering Science and Technology, 10(1), 76-87. https://doi.org/10.4314/ijest.v10i1.6

Panchal, B., Deshmukh, S., & Sharma, M. (2014). Optimization of oil extraction and characterization from Tamarindus indica Linn seed oil. International Journal of Oil, Gas and Coal Engineering, 2(1), 1-6. https://doi.org/10.11648/j.ogce.20140201.11

Popoola, Y. Y., Akinoso, R., & Raji, A. O. (2016). Optimization of oil extraction from giant bushel gourd seeds using response surface methodology. Food Science & Nutrition, 4(5), 759-765. https://doi.org/10.1002/fsn3.341

Predescu, N. C., Papuc, C., Nicorescu, V., Gajaila, I., Goran, G. V., Petcu, C. D., & Stefan, G. (2016). The influence of solid-to-solvent ratio and extraction method on total phenolic content flavonoid content and antioxidant properties of some ethanolic plant extracts. Revista de Chimie, 67(10), 1922–1927.

Rodrigues, C. E. C., Aracava, K. K., & Abreu, F. N. (2010). Thermodynamic and statistical analysis of soybean oil extraction process using renewable solvent. International Journal of Food Science and Technology, 45(11), 2407–2414. https://doi.org/10.1111/j.1365-2621.2010.02417.x

Russin, T. A., Boye, J. I., Arcand, Y., & Rajamohamed, S. H. (2011). Alternative techniques for defatting soy: A practical review. Food and Bioprocess Technology, 4(2), 200-223. https://doi.org/10.1007/s11947-010-0367-8

Sawada, M. M., Venâncio, L. L., Toda, T. A., & Rodrigues, C. E. C. (2014). Effects of different alcoholic extraction conditions on soybean oil yield fatty acid composition and protein solubility of defatted meal. Food Research International, 62, 662–670. https://doi.org/10.1016/j.foodres.2014.04.039

Saxena, D. K., Sharma, S. K., & Sambi, S. S. (2011). Comparative extraction of cottonseed oil. ARPN Journal of Engineering and Applied Sciences, 6(1), 84-89.

Shao, D., Atungulu, G. G., Pan, Z., Yue, T., Zhang, A., & Li, X. (2012). Study of optimal extraction conditions for achieving high yield and antioxidant activity of tomato seed oil. Journal of Food Science, 77(8), 1–7. https://doi.org/10.1111/j.1750-3841.2012.02804.x

Shittu, S. K., Mari, H. H., & Dangora, N. D. (2019). Statistical model for solvent oil extraction from soybean (Glycine max (L.)). Food Research, 3(2), 182-187. https://doi.org/10.26656/fr.2017.3(2).111

Shivani, P., Khushbu, P., Faldu, N., Thakkar, V., & Shubramanian, R. B. (2011). Extraction and analysis of Jatropha curcas L. seed oil. African Journal of Biotechnology, 10(79), 18210-18213. https://doi.org/10.5897/ajb11.776

Tabtabaei, S., & Diosady, L. L. (2013). Aqueous and enzymatic extraction processes for the production of food-grade proteins and industrial oil from dehulled yellow mustard flour. Food Research International, 52(2), 547-556. https://doi.org/10.1016/j.foodres.2013.03.005

Toda, T. A., Sawada, M. M., & Rodrigues, C. E. C. (2016). Kinetics of soybean oil extraction using ethanol as solvent: Experimental data and modeling. Food and Bioproducts Processing, 98, 1-10. https://doi.org/10.1016/j.fbp.2015.12.003

Yolmeh, M., & Jafari, S. M. (2017). Applications of response surface methodology in the food industry processes. Food and Bioprocess Technology, 10(3), 413–433. https://doi.org/10.1007/s11947-016-1855-2

Yujie, C. (2011). Optimization of the technology for preparing soluble dietary fiber from extruded soybean residue. In: Soybean-Applications and Technology, Prof. Tzi-Bun Ng (Ed.). In (pp. 56-66). London: IntechOpen.

Views and opinions expressed in this article are the views and opinions of the author(s), Journal of Food Technology Research shall not be responsible or answerable for any loss, damage or liability etc. caused in relation to/arising out of the use of the content.