Highly Efficient Deacidification Process for Camelina sativa Crude Oil by Molecular Distillation

Recovery and reuse of high-acidity vegetable oil waste (higher content of free fatty acids) is a major concern for reducing their effect on the environment. Moreover, the conventional deacidification processes are known to show drawbacks, such as oil losses or higher costs of wastewater treatment, for which it requires great attention, especially at the industrial scale. This work presents the design of a highly efficient and sustainable process for Camelina sativa oil deacidification by using an ecofriendly method, namely molecular distillation. Experimental studies were performed to identify operating conditions for removing of free fatty acids (FFA) by molecular distillation which involves the oil evaporation in high vacuum conditions. The experimental studies were supported by statistical analysis and technical-economic analysis. Response surface methodology (RSM) was employed to formulate and validate second-order models to predict deacidification efficiency, FFA concentration, and triacylglyceride (TAG) concentration in deodorized oil based on three parameters effects, validated by statistical p-value < 0.05. For a desirability function value of 0.9826, the optimal parameters of evaporator temperature at 173.5 °C, wiper speed at 350 rpm, and feed flowrate at 2 mL/min were selected. The results for process design at optimal conditions (using conventional and molecular distillation methods) showed an efficiency over 92%, a significant reduction in FFA (up to 1%), and an increase in TAG (up to 93%) in refined oil for both methods. From an economical point of view, the deacidification by molecular distillation of Camelina sativa oil is a sustainable process: no wastewater generation, no solvents and water consumption, and lower production costs, obtaining a valuable by-product (FFA).


Introduction
Mono-, di-and triglycerides are the main components of crude oils, in addition to tocopherols, squalene, and sterols in traces (200-800 ppm), whose presence improve oil oxidative stability. Besides these desired components, crude oils contain undesirable impurities as phospholipids (100-500 ppm), free fatty acids (FFA) (5-20%), metal ions and metal complexes (2-15 mg·kg −1 ), oxidized products (2-6 meq·kg −1 ), and moisture (1-3%) [1]. FFA content is one of the most important aspects that influences edible oil quality or their usage as raw material in biodiesel production because they are more susceptible to autoxidation than the esterified fatty acids (FFA act like pro-oxidants, initiating the oxidation mechanism in lipids, leading to rancidity) or produce soaps (by mixing with methanol for transesterification reaction) [1]. Long storage time and temperature over 24 • C, contribute to the degradation and quick rancidity of oils [2]. Moreover, some fresh crude oil extracted by cold pressing can have high acidity and high FFA content.

Experimental Procedures
Two experimental procedures were used for deacidification of Camelina sativa oil: chemical (reaction of FFA with alkali) and physical (molecular distillation) methods.

Camelina sativa Oil Deacidification Using Alkali
This method consisted in treating the oil with an aqueous solution of sodium hydroxide which reacted with FFA, forming soap.
Experiments were performed using the method proposed elsewhere [17]. One gram of alkali was dissolved in 1 mL distilled water and further added to a sample of 150 g Camelina sativa oil. The sample was stirred 30 min, and then 15 mL hexane was added and was stirred again for another 30 min. Finally, the soap was separated at 3500 rpm for 30 min to evaporate the solvent. To remove the soap traces, three successive washes with distilled water at 95 • C were performed, and washing water was removed by centrifugation. Samples were prepared to determine the acidity value and to identify the content of FFA and triacylglyceride (TAG).

Camelina sativa Oil Deacidification by Molecular Distillation
Samples of 500 g Camelina sativa oil were subjected to a wiped film molecular distillation plant (KDL5, UIC GmbH, Alzenau, Germany) to separate FFA based on migration of molecules with different molecular volumes at vacuum conditions. The plant was heated by thermal oil. Working procedure is presented in [18]. The operating conditions were varied in the specific ranges imposed by equipment limits: the evaporator temperature (maximum 250 • C), the rotational wiping speed (250-350 rpm), feed flow rate (2-6 mL/min), while pressure (0.001-100 mbar) and condenser temperature (less than 30 • C). Light and heavy samples were collected in different devices and prepared to determine the acidity value and to identify the content of free fatty acids and triglycerides by chromatographic analysis. Separation process efficiency was calculated with Equation (3): Deacidification efficiency (%) = feed FFA content (g) − bottom product FFA content (g) feed FFA content (g) ·100 Seventeen experiments were performed by varying operating conditions, such as evaporator temperature, the rotational wiping speed, and feed flow rate in equipment range, to evaluate the influence on deacidification efficiency and TAG and FFA compositions.

Samples Characterization 2.3.1. Acidity Value
Camelina sativa oil sample acidity was determined by a method according to the European Standard SR EN No. 660 [19]. The sample was diluted in a mixture of ethanol:diethyl ether 1:1 (v/v) in the presence of phenolphthalein as an indicator and neutralized with 0.5M KOH solution in methanol. The acidity (I A ) was calculated using Equation (4), where V is volume of potassium hydroxide solution used in the titration in cm 3 , n is the normality of solution, m is analyzed sample mass in g, and M KOH is molar weight of potassium hydroxide in g/mol.

FFA and Glycerides Content
Monoacylglycerides (MAG), diacylglycerides (DAG), triacylglycerides (TAG), and FFA compositions were determined by GC-FID using a gas-chromatograph Bruker Scion 436 according to SR EN Standard No. 14105 [20]. For sample analyses, a fused capillary column ZB-5HT type with a nitrogen flow of 5 mL/min was used. Because high temperature is necessary to vaporize the glycerides the injector temperature started at 50 • C and increased with the rate of 50 • C·min −1 until 380 • C with a 2 min hold. Electronic flow controller (EFC) was settled at 2.5 mL/min and the following oven temperature program was used: 50 • C for 1 min, rate 15 • C·min −1 until 180 • C, rate 7 • C·min −1 until 230 • C, rate 10 • C·min −1 until 370 • C. Detector temperature (FID) was set at 380 • C and the gas flows at 15 mL/min for the make-up gas (nitrogen), 15 mL/min for hydrogen, and 300 mL/min for air. The injection volume used was 1 µL. All samples were analyzed in five replicates and the results are presented as the average value ± standard deviation (SD).

Statistical Analysis
To find optimal operating conditions for the deacidification process of Camelina sativa oil, a design of experiments (DOE) was used. DOE main steps followed for this analysis were: setting up experiments objectives, factors and levels selection (number of independent variables manipulated to determine their effect on the measured responses), selection of dependent variables (measured during experiments, their value depends on independent variable effect), design experimental matrix, select ranges of independent variables and number of replicates of experimental data, collect experimental data, regression model formulation and factor effect evaluation, and validation applying statistical ANOVA and regression methods. Different kinds of dependent or independent variables can be selected for the oil deacidification process, and different optimal parameters can be identified for particular cases [21][22][23][24][25][26]. Based on this statistical technique, a second-degree model was formulated to describe the system, including effects or interactions of factors, solved by minimizing/maximizing the response variables. As a DOE technique, the Box-Behnken surface experimental design was used in this study to investigate the effects of three factors (evaporator temperature, wiper rolling speed, and feed flow rate), for some responses of the oil deacidification process (minimum FFA content, maximum TGA content). The independent variable limits depended on equipment limits and raw material composition.
The Box-Behnken design consists of three levels (minimum, center, and maximum values, K = 3) selected for each independent variable, with coded labels (X 1 , X 2 , and X 3 ) and values as presented in Table 1. For the molecular distillation process, three dependent variables were explored: the FFA composition in heavy product (Y 1 ), deacidification efficiency (Y 2 ), and the TAG composition in heavy product (Y 3 ), as presented in Table 2. The Box-Behnken design matrix contains combination of two variables as for 2 2 factorial design, the third variable being maintained at level zero and some central point values. A second-order model (in coded variables) was formulated to provide optimal surfaces for each response Y i using Equation (5).
Considering three X i factors, the model has M coefficients β j,i for each i responses (i = 1 . . . K, j = 0 . . . M) which represent the change of responses Y i per unit of effect of independent variables X i , and K number of factors.
The regression models for responses considered main effects (X 1 , X 2 , X 3 ) determined by the influence of each variable independent from the others, interaction effects (X 1 X 2 , X 1 X 3 , X 2 X 3 ) described as a cumulative effect of two dependent variables and the quadratic effect of each variable (X 2 1 , X 2 2 , X 2 3 ). The model coefficient values were generated by the least squares regression followed by analysis of variance (ANOVA) for model validation. The total number of experiments (N) for second-order model was calculated using Equation (6): where K is the factors number and C p is the replicates at the central point. Thus, for Box-Behnken model design, 17 experimental runs were performed with five replications at the central point. For the Camelina sativa oil deacidification process, the Box-Behnken matrix with three factors in coded form and three responses is shown in Table 3. The optimal values of dependent variables were predicted using desirability function. The model lack of fit was calculated in terms of coefficient of determination (R 2 ) and p-values for 95% confidence interval considered for analysis of variance (ANOVA). The optimal values of factors (minimum/maximum) were identified on their response surfaces for acceptable values of fitted responses evaluating desirability functions.

Process Design, Economic and Environmental Analysis
Molecular distillation processes were designed using Aspen Plus v.10 (AspenTech Technology Inc., Bedford, MA, USA, 2017) simulator to perform mass and energy balances, equipment design, and total costs. Oil sample was considered as a mixture of mono-, di-and triglycerides; tocopherols; and fatty acids identified as majority in composition of Camelina sativa oil. All compounds are in the simulator database and their properties were calculated using RK-ASPEN (Redlich-Kwong Soave with Wong-Sandler mixture rules) properties package [27]. The molecular distillation unit was modeled by combination of three simulator models: HeatX-Flash2-HeatX, for light molecule evaporation, vaporized molecule migration, and condensation, with mixture temperature and pressure as specifications. Economic analysis was performed using Aspen Plus v.10 simulator.

Results and Discussions
According to the GC-FID analysis, the fresh press extraction of Camelina sativa oil, purchased from a regional farm, contained 3.23 ± 0.33% FFA, 0.75 ± 0.10% MAG, 5.95 ± 0.24% DAG, and 90.07 ± 0.32% TAG and has an acidity value of 3.35 ± 0.10 mg KOH/g oil. After six months of storage, chromatographic analysis of the same oil showed a modified composition: 7.67 ± 0.33% FFA, 0.84 ± 0.09% MAG, 6.01 ± 0.30% DAG, and 85.49 ± 0.42% TAG and acidity value of 7.17 ± 0.07 mg KOH/g oil. These results show the great tendency of this oil for oxidation, due to the presence of components with many double bonds.

Camelina Sativa Oil Deacidification Using Alkali
Six experiments were performed to remove FFA from three samples of Camelina sativa oil with initial acidity 3.35 mg KOH/g oil and three samples with 7.17 mg KOH/g oil. After deacidification with sodium hydroxide, the samples' acid value decreased to ≈1 mg KOH/g, respectively, to ≈2.1 mg KOH/g., while the losses were about 18% for all the samples. Deacidification efficiency was around 82% for all samples, and 0.6 g/g oil was collected as wastewater.

Camelina Sativa Oil Deacidification by Molecular Distillation
Camelina sativa oil with 7.17 ± 0.07% mg KOH/g oil acidity was used. Thirty-nine experimental runs were performed to investigate the effects of evaporator temperature (seven experiments), wiper rolling speed (eight experiments), and feed flow rate (20 experi-Sustainability 2021, 13, 2818 7 of 18 ments) on deacidification efficiency. The vacuum conditions of 0.1 Pa were maintained in the system, and condenser temperature at 19 • C for all experiments. FFA, tocopherol, and MAG molecules, having smaller volume than glycerides, are first vaporized at specified pressure conditions and migrate to condenser surface. The evaporation temperature for Camelina sativa oil (as a mixture of FFA, TAG, DAG, MAG, and other components) was selected at 150 • by many experimental trials when the first distillate was collected (evaporation degree between 1% and 10%) at pressure 0.1 Pa. This temperature, much lower than reported in literature [11][12][13][14][15], allowed us to reduce the heat consumption required by the evaporator, even if a high vacuum must be ensured. The wiper rolling speed and feed flowrates were chosen according to the limits of the equipment. Deacidification efficiency was greater than 90% for some operating conditions, no wastewater generated, and the oil losses are around 1-2%. Separated FFA can be valorized as valuable products.

Effect of Evaporator Temperature
Seven molecular distillation experiments were performed at different evaporator temperatures between 150 • C to 180 • C, fixed feed flow rate (2 mL/min), wiper rolling speed (300 rpm), pressure (0.1 Pa), and condenser temperature (19 • C). At this lower pressure, the mixture boiling temperature was not too high and TAG could not degrade (for a separation at 200 • C, it was observed that compounds with lower number of carbons appeared due to the glycerides degradation). During the separation process, light product rich in FFA and heavy product rich in TAG were separated. Increasing the evaporator temperature, the FFA concentration in oil decreased from 3.29% at 150 • C to 1.61% at 180 • C (Figure 1a), 90% for all samples. At temperature greater than 170 • C, the composition of FFA had a smaller change. The composition of TAG in oil was also affected by the evaporation temperature, increasing by 2% from 91.37% at 150 • C to 93.34% at 180 • C ( Figure 1a). The deacidification efficiency increased with the temperature, having a maximum value of 80.39% at 180 • C (Figure 1b).

Camelina Sativa Oil Deacidification by Molecular Distillation
Camelina sativa oil with 7.17 ± 0.07% mg KOH/g oil acidity was used. Thirty-nine experimental runs were performed to investigate the effects of evaporator temperature (seven experiments), wiper rolling speed (eight experiments), and feed flow rate (20 experiments) on deacidification efficiency. The vacuum conditions of 0.1 Pa were maintained in the system, and condenser temperature at 19 °C for all experiments. FFA, tocopherol, and MAG molecules, having smaller volume than glycerides, are first vaporized at specified pressure conditions and migrate to condenser surface. The evaporation temperature for Camelina sativa oil (as a mixture of FFA, TAG, DAG, MAG, and other components) was selected at 150° by many experimental trials when the first distillate was collected (evaporation degree between 1% and 10%) at pressure 0.1 Pa. This temperature, much lower than reported in literature [11][12][13][14][15], allowed us to reduce the heat consumption required by the evaporator, even if a high vacuum must be ensured. The wiper rolling speed and feed flowrates were chosen according to the limits of the equipment. Deacidification efficiency was greater than 90% for some operating conditions, no wastewater generated, and the oil losses are around 1-2%. Separated FFA can be valorized as valuable products.

Effect of Evaporator Temperature
Seven molecular distillation experiments were performed at different evaporator temperatures between 150 °C to 180 °C, fixed feed flow rate (2 mL/min), wiper rolling speed (300 rpm), pressure (0.1 Pa), and condenser temperature (19 °C). At this lower pressure, the mixture boiling temperature was not too high and TAG could not degrade (for a separation at 200 °C, it was observed that compounds with lower number of carbons appeared due to the glycerides degradation). During the separation process, light product rich in FFA and heavy product rich in TAG were separated. Increasing the evaporator temperature, the FFA concentration in oil decreased from 3.29% at 150 °C to 1.61% at 180 °C (Figure 1a), 90% for all samples. At temperature greater than 170 °C, the composition of FFA had a smaller change. The composition of TAG in oil was also affected by the evaporation temperature, increasing by 2% from 91.37% at 150 °C to 93.34% at 180 °C ( Figure 1a). The deacidification efficiency increased with the temperature, having a maximum value of 80.39% at 180 °C ( Figure 1b).

Effect of Wiper Rolling Speed
Experimental studies were performed at five values between 250 and 350 rpm at two evaporator temperatures (160 • C and 170 • C), keeping the fix flow rate of 2 mL/min at pressure 0.1 Pa and condenser temperature 19 • C. As can be seen in Figure 2, the wiper rolling element (model UIC GmbH, Alzenau, Germany) had a great influence in the deacidification process by molecular distillation due to its role of dispersing liquid on the evaporator surface in a thin film and liquid mixing.

Effect of Wiper Rolling Speed
Experimental studies were performed at five values between 250 and 350 rpm at two evaporator temperatures (160 °C and 170 °C), keeping the fix flow rate of 2 mL/min at pressure 0.1 Pa and condenser temperature 19 °C. As can be seen in Figure 2, the wiper rolling element (model UIC GmbH, Alzenau, Germany) had a great influence in the deacidification process by molecular distillation due to its role of dispersing liquid on the evaporator surface in a thin film and liquid mixing. The increasing speed of wiper rolling decreased oil acidity from 2.46% to 0.88% at 170 °C (results presented in Figure 2a) and from 3.12% to 0.89% at 160 °C. Additionally, deacidification efficiency showed a significant increase at high speeds (88.83% at 160 °C and 91.75% at 170 °C). TAG concentration showed a moderate increase from almost 91.5% at 250 rpm to 93% at 350 rpm with small differences at the two evaporation temperatures studied. Under the same conditions, calculated deacidification efficiency showed an increase with 28% (for 160 °C) and 20% (for 170 °C), for an increase in the wiper rolling rotations with 100 rpm, due to a uniform oil distribution on the evaporator walls facilitating the evaporation (Figure 2b).

Effect of Feed Flow Rate
These studies were performed at an evaporator temperature of 170 °C and different values of feed flowrates (2 mL/min, 4 mL/min and 6 mL/min) with wiper speed at 250 rpm, 300 rpm, and 350 rpm at constant pressure (0.1 Pa) and condenser temperature (19 °C). The TAG composition decreased and FFA composition in bottom product increased when the flow rates increased. At 6 mL/min, the lowest FFA could be greater than 2.5% ( Figure 3a). The increasing speed of wiper rolling decreased oil acidity from 2.46% to 0.88% at 170 • C (results presented in Figure 2a) and from 3.12% to 0.89% at 160 • C. Additionally, deacidification efficiency showed a significant increase at high speeds (88.83% at 160 • C and 91.75% at 170 • C). TAG concentration showed a moderate increase from almost 91.5% at 250 rpm to 93% at 350 rpm with small differences at the two evaporation temperatures studied. Under the same conditions, calculated deacidification efficiency showed an increase with 28% (for 160 • C) and 20% (for 170 • C), for an increase in the wiper rolling rotations with 100 rpm, due to a uniform oil distribution on the evaporator walls facilitating the evaporation (Figure 2b).

Effect of Feed Flow Rate
These studies were performed at an evaporator temperature of 170 • C and different values of feed flowrates (2 mL/min, 4 mL/min and 6 mL/min) with wiper speed at 250 rpm, 300 rpm, and 350 rpm at constant pressure (0.1 Pa) and condenser temperature (19 • C). The TAG composition decreased and FFA composition in bottom product increased when the flow rates increased. At 6 mL/min, the lowest FFA could be greater than 2.5% (Figure 3a). At low wiper speed, FFA content did not vary too much, but at high values, this difference became significant. This was determined by the fact that the feed flow rate influenced the system heat transfer: more energy for evaporation was needed if flow rate was increasing. In Figure 3b, the influence of feed flow rate versus deacidification effi- At low wiper speed, FFA content did not vary too much, but at high values, this difference became significant. This was determined by the fact that the feed flow rate influenced the system heat transfer: more energy for evaporation was needed if flow rate was increasing. In Figure 3b, the influence of feed flow rate versus deacidification efficiency is presented at different wiper rolling speeds. It was noticed that, at low speeds, the decrease was less than 5%, but at high speed, the deacidification efficiency decreased by 10% (Figure 3b).

Fitting Models
Responses as FFA content in the bottom product (Y 1 ), deacidification efficiency (Y 2 ), and TAG concentration in the bottom product (Y 3 ) were fitted by a second-order model with regression coefficients presented in Table 4, using a Box-Behnken design of experimental data. The fitted response surfaces were calculated with Equations (7)-(9): For each response, all coefficients, and their p-value with 95% confidence interval were computed to validate the significance of regression models. The regression model is significant at 5% level if the p-value is smaller than 0.05. From Table 4, it can be observed that some terms did not have a significant effect: quadratic X 2 and interaction X 1 X 3 for Y 1 , only quadratic X 3 for Y 2 and quadratic X 1 , and interaction X 1 X 2 for Y 3 . In Table 5, the experimental and data predicted with proposed second-order models are presented. It was observed that close correlations existed between the data. Also, the good correlation can be noticed in Figure 4 when predicted data are plotted versus experimental data for all three responses, all points being placed near a 45 • line, especially for Y 2 .  For each response, all coefficients, and their p-value with 95% confidence interval were computed to validate the significance of regression models. The regression model is significant at 5% level if the p-value is smaller than 0.05. From Table 4, it can be observed that some terms did not have a significant effect: quadratic X2 and interaction X1X3 for Y1, only quadratic X3 for Y2 and quadratic X1, and interaction X1X2 for Y3. In Table 5, the experimental and data predicted with proposed second-order models are presented. It was observed that close correlations existed between the data. Also, the good correlation can be noticed in Figure 4 when predicted data are plotted versus experimental data for all three responses, all points being placed near a 45° line, especially for Y2.     Table 6 shows the analysis of variance (ANOVA) for all the responses investigated in this study. The models, linear, quadratic, and interaction terms of the three independent variables (temperature, speed, and flow rate), lack of fit (model error), and pure error (replicates error) were evaluated according to their degree of freedom, sum of squares, mean square, F-value, p-value, and model coefficient of determination R 2 . A model is good if it has a high F-value, low p-value, insignificant lack of fit, and high R 2 . The results of this analysis are presented in Table 6. For the FFA content (Y 1 ), the ANOVA analysis of the quadratic model showed that the F-value of 21.156 implies that the model is significant at p < 0.001. The coefficient of determination R 2 is 0.965, measuring how properly experimental data are approximated the regression model. The lack of fit was not significant (p > 0.05) and indicated that the model prediction is adequate. Examination of the coefficients for the first response indicated that the linear terms as evaporator temperature, wiper speed, and flow rate and the interaction terms between wiper speed and flow rate were highly significant with higher F-value and lower p-value. Quadratic terms of evaporator temperature and flow rate and interaction term between temperature and speed were significant, while quadratic term of speed and interaction term between temperature and flow rate were not significant with low F-value and p > 0.05. For the deacidification efficiency (Y 2 ), the ANOVA analysis showed that a greater F-value of 2123.15 implies that the model is significant at p < 0.001. The coefficient of determination R 2 was 0.999 and all the linear, quadratic, and interaction terms were significant (p < 0.05), except the quadratic effect of flow which was insignificant at p > 0.05. To determine if the models were compatible with experimental data, lack of fit test was calculated from residuals of fitting. The model was adequate, as indicated by the error analysis that showed insignificant lack of fit with p > 0.05. In case of the third response, TAG concentration (Y 3 ), the obtained results showed that the model was significant with a F-value of 29.522 and p < 0.001. The coefficient of determination R 2 was 0.974. All the linear terms were highly significant with high F-value and low p-value, while the quadratic term of evaporator temperature and the interaction term between temperature and speed were insignificant with low F-value and p > 0.05.

Optimization by Using Response Surface Methodology
This optimization method was used for Camelina sativa oil deacidification by molecular distillation process to locate one point or region/regions on factor surfaces, where fitted responses were acceptable (optimal) values. Response surfaces were computed for factor ranges between minimum and maximum values of fitted responses obtained with secondorder model. The shape of these surfaces gave information about regions where responses had optimal values and the ranges of optimal factors that affected these responses. As can be seen from Y 1 surfaces (Figure 5a), FFA composition should be at a minimum value in regions of higher X 1 and X 2 and lower X 3 .

Figure 5.
Response surfaces for Y 1 , Y 2 , and Y 3 based on: (a) X 1 X 2 effects, X 3 = ct, (b) X 1 X 3 effects, X 2 = ct, (c) X 2 X 3 effects, X 1 = ct. These factors will have a positive impact on FFA composition (Y 1 ≈ 1% FFA). Y 2 surfaces had maximum regions (90-95%), identified in Figure 5b also in the same regions, for higher X 1 and X 2 and lower X 3 . In Figure 5c, Y 3 response surfaces in the same regions for maximum value of 93.5% are presented. First response (Y 1 -FFA) must be minimized while second (Y 2 -DE) and third (Y 3 -TAG concentration) responses must be maximized.
2D plots were computed to identify optimal values of responses when each factor had a variation between −1 and 1, and the other factors had constant values (from ranges identified on response surfaces). In Figure 6a, three profiles of response Y 1 affected by X 1 , X 2 and X 3 are presented, these profiles showing minimum values.  Computing the individual desirability functions with STATISTICA© for factor values between −1 and 1; their profiles are presented in Figure 7. Profiles presented in Figure 6b,c represented the responses of Y 2 and Y 3 affected by the same factors. These profiles presented maximum values. All three response profiles presented the optimal values at X 1 = 0.34629, X 2 = 1, and X 3 = −1. These values were in the regions identified in the regions from response surfaces.
For all responses, individual desirability functions were computed with values between 0 and 1, where 0 was for an undesirable response (response is greater than a target value for minimum point or less than the smallest value for maximum) and 1 indicated a completely desirable response (response is less than the smallest value for minimum point or greater than a target value for maximum point). Intermediate values of responses between smallest value and target value represented desirable responses. Smallest, medium, and target values for desirability functions are presented in Table 7, considering a linear variation for intermediate values. Computing the individual desirability functions with STATISTICA© for factor values between −1 and 1; their profiles are presented in Figure 7. In Table 8, the optimal values of responses and factors are centralized. Final desirability function, defined as geometric mean of individual desirability functions, was 0.9826 for all factors combined. These optimal values of coded factors corresponded to optimal operating conditions such as: evaporator temperature of 173.5 °C, wiper rolling speed of 350 rpm, and a feed flow rate of 2 mL/min for the deacidification process of Camelina sativa oil by molecular distillation. At these conditions, separation parameters for this process were minimum composition of FFA in bottom product of 1.007% FFA, maximum deacidification efficiency of 92.167%, and maximum concentration of TAG in bottom product of 93.360% TAG, with a desirability of 0.9826.
To validate the model prediction, two parallel trials were performed in the optimal operating conditions. FFA concentration of 1.15 ± 0.2%, a TAG concentration of 92.85 ± 0.05%, with a deacidification efficiency of 91.30 ± 0.1% were obtained, which is in good agreement with theoretical prediction.

Process Design, Economic and Environmental Analysis
In Figure 8, the process flowsheets for two deacidification processes of Camelina sativa oil are presented: ODAW-conventional method with alkali ( Figure 8a) and ODMDmolecular distillation (Figure 8b). A 100 kg/batch of crude oil was fed in each process and 92.5 kg/batch of refined oil was obtained from ODMD process with 12.8% more than in ODAW process. In Table 8, the optimal values of responses and factors are centralized. Final desirability function, defined as geometric mean of individual desirability functions, was 0.9826 for all factors combined. These optimal values of coded factors corresponded to optimal operating conditions such as: evaporator temperature of 173.5 • C, wiper rolling speed of 350 rpm, and a feed flow rate of 2 mL/min for the deacidification process of Camelina sativa oil by molecular distillation. At these conditions, separation parameters for this process were minimum composition of FFA in bottom product of 1.007% FFA, maximum deacidification efficiency of 92.167%, and maximum concentration of TAG in bottom product of 93.360% TAG, with a desirability of 0.9826. To validate the model prediction, two parallel trials were performed in the optimal operating conditions. FFA concentration of 1.15 ± 0.2%, a TAG concentration of 92.85 ± 0.05%, with a deacidification efficiency of 91.30 ± 0.1% were obtained, which is in good agreement with theoretical prediction.

Process Design, Economic and Environmental Analysis
In Figure 8, the process flowsheets for two deacidification processes of Camelina sativa oil are presented: ODAW-conventional method with alkali ( Figure 8a) and ODMDmolecular distillation (Figure 8b). A 100 kg/batch of crude oil was fed in each process and 92.5 kg/batch of refined oil was obtained from ODMD process with 12.8% more than in ODAW process.
The annual number of batches was different for these two processes, due to the different operational time, 15 h for ODMD, and 50 h for the other process. ODAW process involved a reaction step, free fatty acids reacted with alkali to form soap, which retained 5% of crude vegetable oil. Three filtering and washing steps were needed to removed soaps and the 60 kg/batch of resulting wastewater needed to be treated. ODMD process involved a vacuum step and heating and cooling agents to assure evaporator temperature at 173.5 • C (to vaporize FFA) and 19 • C for cooling zone (to condense FFA molecules). Different tanks were used for raw materials and products. The annual number of batches was different for these two processes, due to the different operational time, 15 h for ODMD, and 50 h for the other process. ODAW process involved a reaction step, free fatty acids reacted with alkali to form soap, which retained 5% of crude vegetable oil. Three filtering and washing steps were needed to removed soaps and the 60 kg/batch of resulting wastewater needed to be treated. ODMD process involved a vacuum step and heating and cooling agents to assure evaporator temperature at 173.5 °C (to vaporize FFA) and 19 °C for cooling zone (to condense FFA molecules). Different tanks were used for raw materials and products.
The main advantages of the ODMD process are the reduction of crude oil losses from 18% to 1.5%, and no wastewater generation. For this process, heating/cooling/vacuum systems costs are increased. In Table 9, raw materials, utility consumptions, and the production cost estimation for both processes are presented.   The main advantages of the ODMD process are the reduction of crude oil losses from 18% to 1.5%, and no wastewater generation. For this process, heating/cooling/vacuum systems costs are increased. In Table 9, raw materials, utility consumptions, and the production cost estimation for both processes are presented. Production cost was evaluated as a contribution of variable costs (equipment, raw materials, utilities, and treatment costs), fixed costs (operating, non-operating and administration labor costs, land, and others) and depreciation costs (Figure 9). Raw material specific costs ($/kg deacidified oil) were 40% higher for the alkali process due to the consumption of alkali solutions, wastewaters, and solvents. Also, the variable costs for the alkali process were influenced by treatment. For the molecular distillation process, the higher contribution of production costs was given by raw materials and utility costs.
Depreciation, $/kg deacidified oil 0.15 0.09 Cost of production, $/kg deacidified oil 3.39 2.34 Production cost was evaluated as a contribution of variable costs (equipment, raw materials, utilities, and treatment costs), fixed costs (operating, non-operating and administration labor costs, land, and others) and depreciation costs (Figure 9). Raw material specific costs ($/kg deacidified oil) were 40% higher for the alkali process due to the consumption of alkali solutions, wastewaters, and solvents. Also, the variable costs for the alkali process were influenced by treatment. For the molecular distillation process, the higher contribution of production costs was given by raw materials and utility costs.

Conclusions
The efficiency of the Camelina sativa oil deacidification process by molecular distillation was determined by combining experimental studies (to identify operating conditions) with statistical analysis (for optimal conditions) and technical-economic analysis (to evaluate production costs). The results were compared with a classical process, using alkali solutions to remove FFA. Experimental studies for deacidification of Camelina sativa oil with 7.67% FFA using both methods led to a process efficiency over 80% and final FFA content lower than 1%, but the conventional process generated wastewater with soap and 18% oil loss.
Molecular distillation involves the mixture evaporation in vacuum conditions, and trial experiments were needed to obtain operation parameters. Changing values of some operating conditions, such as evaporator temperature, wiper rolling speed, and feed flow rate in a specified range, led to a modification of deacidification efficiency with 35.44%. Moreover, 88% of FFA were removed from oil, while TAG composition in refined oil varies by 8%. Optimal operation conditions were performed using DOE, formulating a Box-Behnken design with three factors and three responses: evaporator temperature was 173.5 °C, wiper speed was 350 rpm, and feed flow rate was 2 mL/min, evaluated for a process desirability of 0.9826. In these conditions, the deacidification efficiency was 92.167%, FFA composition in deacidified oil was 1.007%, and TAG concentration was 93.36%. The technical and economic analysis accomplished for both deacidification processes (ODAW and ODMD), for 100 kg/batch Camelina sativa oil in optimal conditions, allows the evaluation of processes sustainability: raw materials and utility consumption per kg of deacidified Camelina sativa oil (1.22 kg vs. 1.08 kg and 1 kWh vs. 0.25 kWh), wastewater generation

Conclusions
The efficiency of the Camelina sativa oil deacidification process by molecular distillation was determined by combining experimental studies (to identify operating conditions) with statistical analysis (for optimal conditions) and technical-economic analysis (to evaluate production costs). The results were compared with a classical process, using alkali solutions to remove FFA. Experimental studies for deacidification of Camelina sativa oil with 7.67% FFA using both methods led to a process efficiency over 80% and final FFA content lower than 1%, but the conventional process generated wastewater with soap and 18% oil loss.
Molecular distillation involves the mixture evaporation in vacuum conditions, and trial experiments were needed to obtain operation parameters. Changing values of some operating conditions, such as evaporator temperature, wiper rolling speed, and feed flow rate in a specified range, led to a modification of deacidification efficiency with 35.44%. Moreover, 88% of FFA were removed from oil, while TAG composition in refined oil varies by 8%. Optimal operation conditions were performed using DOE, formulating a Box-Behnken design with three factors and three responses: evaporator temperature was 173.5 • C, wiper speed was 350 rpm, and feed flow rate was 2 mL/min, evaluated for a process desirability of 0.9826. In these conditions, the deacidification efficiency was 92.167%, FFA composition in deacidified oil was 1.007%, and TAG concentration was 93.36%. The technical and economic analysis accomplished for both deacidification processes (ODAW and ODMD), for 100 kg/batch Camelina sativa oil in optimal conditions, allows the evaluation of processes sustainability: raw materials and utility consumption per kg of deacidified Camelina sativa oil (1.22 kg vs. 1.08 kg and 1 kWh vs. 0.25 kWh), wastewater generation (0.73 kg vs. 0 kg), and production costs ratio (45% higher for ODAW process). The results showed that using ODMD method some drawbacks of oil deacidification processes can be removed: the method can be applied for oils with different acidities and different compositions, oil losses are minimized, and wastewater treatment are reduced.