Development of Molecular Distillation Based Simulation and Optimization of Refined Palm Oil Process Based on Response Surface Methodology

The deodorization of the refined palm oil process is simulated here using ASPEN HYSYS. In the absence of a library molecular distillation (MD) process in ASPEN HYSYS, first, a single flash vessel is considered to represent a falling film MD process which is simulated for a binary system taken from the literature and the model predictions are compared with the published work based on ASPEN PLUS and DISMOL. Second, the developed MD process is extended to simulate the deodorization process. Parameter estimation technique is used to estimate the Antoine’s parameters based on literature data to calculate the pure component vapor pressure. The model predictions are then validated against the patented results of refining edible oil rich in natural carotenes and vitamin E and simulation results were found to be in good agreement, within a ±2% error of the patented results. Third, Response Surface Methodology (RSM) is employed to develop non-linear second-order polynomial equations based model for the deodorization process and the effects of various operating parameters on the performance of the process are studied. Finally, an optimization framework is developed to maximize the concentration of beta-carotene, tocopherol and free fatty acid while optimizing the feed flow rate, temperature and pressure subject to process constrains. The optimum results of feed flow rate, temperature, and pressure were determined as 1291 kg/h, 147 ◦C and 0.0007 kPa respectively, and the concentration responses of betacarotene, tocopherol and free fatty acid were found to be 0.000575, 0.000937 and 0.999840 respectively.


Introduction
Palm oil is consumed worldwide.It is known that edible palm oil is full of nutrients that contribute to daily energy needs.The general compositions of crude palm oil are 95% of triglycerides, 3%-5% of fatty acid and 1% of minor trace component (phytonutrient).Palm oil is rich in phytonutrients; it provides a valuable source of beta-carotene and tocopherol (vitamin E).Crude palm oil has carotenes concentration of 1000 ppm-3000 ppm [1].The carotenes from palm oil are a source of pro-vitamin A and recent findings have demonstrated that beta-carotene is used to inhibit tumor progression and hence, reduces cancer formation in humans.The purpose of the refined palm oil process is to remove unwanted fatty acid and impurities.It can be classified as both chemical and physical refining method.However, the physical method is proven to be more effective compared to the chemical method in terms of food safety and environment [2].The physical method consists of three processes; degumming, bleaching and deodorization.The pre-process of deodorization involves operation at a temperature higher than 200 • C at normal pressure.It is performed to remove odor from the oil and additional fat.It generally results in not only removing fatty acid but also destroying all carotenes.Most of the carotenes in palm oil are destroyed in the conventional refining process to produce a lighter color oil [3].
The molecular distillation (MD) technique is a viable process to recover carotene from the palm oil and is suitable for such high a value-added product [4].It is an important thermal separation technology.Typical application involves high boiling substances mainly of temperature-sensitive compounds.The compounds are distilled at a pressure range from 1 atmosphere to high vacuum [5].At this pressure the evaporation temperature will fall during the short residence time, hence, no thermal decomposition can occur.Due to pressure being low as vacuum sometimes, it is known as a vacuum distillation column.Mean free path has an effect on this separation process.The mean free path of MD column is longer than the distance of evaporator and condenser.As a result, the molecules reach the condensing surface without being swerved on collision of others foreign gas.The larger molecule size has a shorter mean free path.Therefore, MD separates substances by its molecular weight.The main objectives of this research are divided into four sections: (a) MD process development and validation for a binary system from literature; (b) Application of the MD process for palm oil deodorization and validation of the process against patented experimental data; (c) Design of experiment by Response Surface Methodology (RSM) leading to the development of non-linear equations based MD model for palm oil deodorization; and (d) Optimization of the deodorization process based on the non-linear equations developed in (c).Note, there are many prevalent applications of RSM in food science, chemical and biological engineering related to molecular distillation [6][7][8][9].
Optimization in engineering shows an important methodology to enhance the process rather than making decisions intuitively.Many techniques and algorithms, such as neural network, genetic algorithm, Monte Carlo, Excel solver and others have been successfully applied in various processes.Here, Excel solver has been used for optimizing the operation of the refined palm oil process.The simulation, design of experiment and optimization are presented and discussed in later sections.

Molecular Distillation Process: Modeling and Validation
In this work, MD is used for the deodorization of palm oil.The objective of this section is to develop and simulate a molecular distillation process.Since there was no library model available in ASPEN HYSYS to simulate a falling film molecular distillation therefore, a simulation of single flash vessel was proposed to represent MD, as considered by other authors in previous studies [10].In order to ensure that the simulation of MD in ASPEN HYSYS was valid, we referred to the study by the authors of [10] for the development of a computational tool for simulating MD process in oil refining using ASPEN PLUS and validation of the results against those predicted by DISMOL (a dedicated molecular distillation process simulator).In their study, an equi-molar mixture of Dibutyl phthalate (DBP) and Dibutyl sebacate (DBS) was fed to the flash process at pressure 0.133 Pa with an initial feed flow rate of 50 kg/h.Temperature was manipulated to achieve the same distillation mass ratio with DISMOL simulator (0.2120).Their ASPEN PLUS results are shown in Figures 1 and 2 together with the results of this work, which used ASPEN HYSYS with Redlich Kwong Soave with Wong Sandler mixture rules (RKSWS) for the calculation of thermo-physical properties.(  1.It is clear that the predictions of ASPEN HYSYS are close to that of ASPEN PLUS in terms of temperature and close to DISMOL in terms of Residue DBP molar fraction.In terms of distillation of DBP molar fraction, both ASPEN PLUS and ASPEN HYSYS predictions are within ±2% of that predicted by DISMOL.Thus, the developed flash process within ASPEN HYSYS can be applied to any MD process, such as the palm oil deodorization process considered in this work.(  1.It is clear that the predictions of ASPEN HYSYS are close to that of ASPEN PLUS in terms of temperature and close to DISMOL in terms of Residue DBP molar fraction.In terms of distillation of DBP molar fraction, both ASPEN PLUS and ASPEN HYSYS predictions are within ±2% of that predicted by DISMOL.Thus, the developed flash process within ASPEN HYSYS can be applied to any MD process, such as the palm oil deodorization process considered in this work.The plotted data between ASPEN PLUS and ASPEN HYSYS seem to be in very good agreement.In order to achieve the same distillation mass ratio obtained by DISMOL, the flash equipment for ASPEN PLUS and ASPEN HYSYS had to operate at a temperature of 63.1 • C and 63.5 • C respectively.Both temperatures result in higher molar fraction of DBP in distillate and lower DBP at residue.Equation 1 was proposed by the authors of [10] to transform the flash result into MD.The efficiency factor is multiplied with the molar fraction for DBP at distillate and DBP molar fraction at residue before mass balance equations are solved.1.It is clear that the predictions of ASPEN HYSYS are close to that of ASPEN PLUS in terms of temperature and close to DISMOL in terms of Residue DBP molar fraction.In terms of distillation of DBP molar fraction, both ASPEN PLUS and ASPEN HYSYS predictions are within ±2% of that predicted by DISMOL.Thus, the developed flash process within ASPEN HYSYS can be applied to any MD process, such as the palm oil deodorization process considered in this work.

Refined Palm Oil Deodorization Process Simulation by ASPEN HYSYS
The MD process model developed and validated in Section 2 is now used for the simulation of the deodorization of refined palm oil process in ASPEN HYSYS. Figure 3 shows the MD process for deodorization in refined palm oil.The operation is carried out at a high vacuum pressure between 5 × 10 −4 and 8 × 10 −4 kPa, at a temperature of 175 • C-180 • C and a feed flow rate between 250 and 2000 kg/h [11].The vapor pressure estimations and the simulation of refined palm oil are discussed below.The MD process model developed and validated in Section 2 is now used for the simulation of the deodorization of refined palm oil process in ASPEN HYSYS. Figure 3 shows the MD process for deodorization in refined palm oil.The operation is carried out at a high vacuum pressure between 5 × 10 −4 and 8 × 10 −4 kPa, at a temperature of 175 °C-180 °C and a feed flow rate between 250 and 2000 kg/h [11].The vapor pressure estimations and the simulation of refined palm oil are discussed below.

Vapor Pressure Estimation
To design and simulate a molecular distillation, thermo-physical properties, vapor pressure of each component and the thermodynamic model are all required.This process used Redlich Kwong ASPEN (RK-Aspen) as a thermodynamic model or property method to calculate physical properties.The RK-Aspen is an equation of state model (EOS) that is good for vapor phase modeling with fewer binary parameters required and is consistent in critical region [12].
The feed stream of the process consists of five main components: tripalmitin triolein, oleic acid, tocopherol and beta-carotene.The mass fraction of each component is shown in Table 2.For a multicomponent system, it was not an easy task to calculate molar fraction by use of mass and component balances in an underdetermined situation.Consequently, thermodynamic properties and vapor pressure or pure component parameters for Antoine equations are required, and must be established.The thermo-physical properties are introduced in the simulation in terms of physical properties such as boiling point (TB), critical temperature (TC), critical pressure (PC), acentric factor (ω) and liquid molar volume (Vl20) for each component, and can be seen in Table A1 in the Appendix, as accomplished by Lim et al. [13].Equation ( 2) refers to an extension of the Antoine equation from Lim et al.'s work, and its parameters are shown in Table 3.However, the equation available in ASPEN HYSYS is shown in Equation ( 3).Therefore, the parameters of the Antoine equation used in this study were obtained by minimizing the least square error of Equation ( 4) through Solver Add-in, in Microsoft Excel and Matlab optimization.

Vapor Pressure Estimation
To design and simulate a molecular distillation, thermo-physical properties, vapor pressure of each component and the thermodynamic model are all required.This process used Redlich Kwong ASPEN (RK-Aspen) as a thermodynamic model or property method to calculate physical properties.The RK-Aspen is an equation of state model (EOS) that is good for vapor phase modeling with fewer binary parameters required and is consistent in critical region [12].
The feed stream of the process consists of five main components: tripalmitin triolein, oleic acid, tocopherol and beta-carotene.The mass fraction of each component is shown in Table 2.For a multicomponent system, it was not an easy task to calculate molar fraction by use of mass and component balances in an underdetermined situation.Consequently, thermodynamic properties and vapor pressure or pure component parameters for Antoine equations are required, and must be established.The thermo-physical properties are introduced in the simulation in terms of physical properties such as boiling point (T B ), critical temperature (T C ), critical pressure (P C ), acentric factor (ω) and liquid molar volume (V l20 ) for each component, and can be seen in Table A1 in the Appendix, as accomplished by Lim et al. [13].Equation (2) refers to an extension of the Antoine equation from Lim et al.'s work, and its parameters are shown in Table 3.However, the equation available in ASPEN HYSYS is shown in Equation (3).Therefore, the parameters of the Antoine equation used in this study were obtained by minimizing the least square error of Equation ( 4) through Solver Add-in, in Microsoft Excel and Matlab optimization.
where A, B, C, D and E are the Antoine Parameters.P is the vapor pressure in kPa, and T is temperature in K.
where S is sum of the squared error, m is the number of experimental data point within the temperature range of 373 K-503 K, P calc i is the calculated vapour pressure obtained from Lim et al.'s work [13] using Equation ( 2) and P exp i is the expected vapour pressure in Equation (3) (as available in ASPEN HYSYS).The pure component parameters for modified Antoine equation (Equation ( 3)) and coefficient of determination (R 2 ) of each component are summarized in Table 4.The coefficient of determination was calculated by regression data analysis in Microsoft Excel between the correlation of lnP calc and lnP exp as Equation (5).
where ln P

Refined Palm Oil Deodorization Process Simulation and Validation
The main task of this part is to simulate, validate, and compare the performance of molecular distillation with 4 experiments reported in the patent of refining of edible oil rich in natural carotene and vitamin E [11].The distillation mass ratio for the process was 0.036.The MD process was simulated, and prediction error of the simulations is shown in Table 5.
Table 5. Simulations results of refining of edible oil in study [11]  The carotene recovery of the patent recorded in [11] and the current work is demonstrated in Figure 4.It can be clearly seen that the refined palm oil deodorization process simulation results show a very good agreement with the patent, with a prediction error of less than 3%.[11] and this work.

RSM-Based Process Model
While we could use the ASPEN HYSYS based process model developed in Section 3 for further simulation and optimisation, here, we have decided to develop a Response Surface Methodology (RSM)-based process model and we resort to Design of Experiments (DoE).Like many authors [14,15] who have considered process simulator as an alternative to real process experiments and applied RSM on the simulated data in the past, we assume the ASPEN HYSYS based MD process described in Section 3 as the real process for the deodorization of refined palm oil.
The design of experiment (DoE) consists of three stages which are (i) screening methodology (preliminary research) to identify the vital factors, (ii) response surface methodology and (iii) model validation.The screening method is to detect factors that impact the process.In preliminary research for molecular distillation [6][7][8][9], it is found that temperatures at feed container, temperature at condenser and rotation speed of the design are the factors that had effects on the concentration of beta-carotene, tocopherol and free fatty acid.Therefore, these parameters are kept constant.Both temperatures at feed container and at condenser are fixed at 150 °C.In this work, the feed flow rate, column pressure and column temperature are manipulated for further investigation of the process.
RSM with Central Composite Design (CCD) is applied to develop, and obtain a quadratic model.The CCD consists of three portions which are the fraction of 2 k factorial design (level are coded as 1 and −1), 2k axial portion and the number of center point.In this case, the least number of 3 center points is used.Therefore, the design contained a total of 17 simulation runs according to Equation ( 6) [16]: where N is the number total runs, k is independent variable, n0 is number of center points.
The feed flow rate (x1), distillation temperature (x2), and pressure (x3) were independent variables.These three variables are used and assumed to influence the responses of beta-carotene concentration (y1), tocopherol concentration (y2) at residue stream and free fatty acid (oleic acid) concentration (y3) in the distillate stream.
In CCD, the axial value is set as 0.59.High value and low value are coded as 1 and -1

RSM-Based Process Model
While we could use the ASPEN HYSYS based process model developed in Section 3 for further simulation and optimisation, here, we have decided to develop a Response Surface Methodology (RSM)-based process model and we resort to Design of Experiments (DoE).Like many authors [14,15] who have considered process simulator as an alternative to real process experiments and applied RSM on the simulated data in the past, we assume the ASPEN HYSYS based MD process described in Section 3 as the real process for the deodorization of refined palm oil.
The design of experiment (DoE) consists of three stages which are (i) screening methodology (preliminary research) to identify the vital factors, (ii) response surface methodology and (iii) model validation.The screening method is to detect factors that impact the process.In preliminary research for molecular distillation [6][7][8][9], it is found that temperatures at feed container, temperature at condenser and rotation speed of the design are the factors that had effects on the concentration of beta-carotene, tocopherol and free fatty acid.Therefore, these parameters are kept constant.Both temperatures at feed container and at condenser are fixed at 150 • C. In this work, the feed flow rate, column pressure and column temperature are manipulated for further investigation of the process.
RSM with Central Composite Design (CCD) is applied to develop, and obtain a quadratic model.The CCD consists of three portions which are the fraction of 2 k factorial design (level are coded as 1 and −1), 2k axial portion and the number of center point.In this case, the least number of 3 center points is used.Therefore, the design contained a total of 17 simulation runs according to Equation (6) [16]: where N is the number total runs, k is independent variable, n 0 is number of center points.The feed flow rate (x 1 ), distillation temperature (x 2 ), and pressure (x 3 ) were independent variables.These three variables are used and assumed to influence the responses of beta-carotene concentration (y 1 ), tocopherol concentration (y 2 ) at residue stream and free fatty acid (oleic acid) concentration (y 3 ) in the distillate stream.
In CCD, the axial value is set as 0.59.High value and low value are coded as 1 and -1 respectively while 0 was the center point code [17].The coded values and its design range are shown in Tables 6  and 7 respectively.Those values are coded according to the following Equations ( 7)-( 9) [17].
where X i and x i are coded and actual values of variable i, respectively; and subscript i are 1 = feed flow rate, 2 = column temperature and 3 = column pressure, and H, L are highest and lowest values.The regression analysis of the responses is fitted with non-linear second-order polynomial equation represented by Equation ( 10): Processes 2017, 5, 40 where y is the responses, β 0 is constant number, β i , β ii and β ij are coefficients for the linear, quadratic and interaction, respectively.Table 8 shows the experimental (simulated by MD process described in Section 3) results for 17 runs of independent variables with respect to responses: (y 1 ), (y 2 ) and (y 3 ).The regression coefficient value for linear, quadratic and interaction are shown in Table 9.The effect of each independent variable on the responses according to 17 runs is described below.

The Effect of Feed Flow Rate on Responses
The results of feed flow rate on the responses of beta-carotene (y 1 ), tocopherol (y 2 ) and FFA (y 3 ) are investigated at a temperature of 150 • C and pressure of 5.05 × 10 −4 kPa (run No. 2 and 3).Decreasing the feed flow rate from 1500 to 1000 kg/h at constant temperature and pressure does not make any changes to the responses.The contour plots and response surface curves in Figures 5-7, reveal that the feed flow rate does not show significant deviation on responses.

The Effect of Feed Flow Rate on Responses
The results of feed flow rate on the responses of beta-carotene (y1), tocopherol (y2) and FFA (y3) are investigated at a temperature of 150 °C and pressure of 5.05 × 10 −4 kPa (run No. 2 and 3).Decreasing t h e feed flow rate from 1500 to 1000 kg/h at constant temperature and pressure does not make any changes to the responses.The contour plots and response surface curves in Figures 5-7, reveal that the feed flow rate does not show significant deviation on responses.

The Effect of Temperature on Responses
Concentration of beta-carotene, tocopherol and FFA as a function of temperature are investigated at constant feed flow rate and pressure at 1500 kg/h and 5.05 × 10 −4 kPa respectively (run no. 4 and 5).In Figures 8-10, the temperature shows significant deviation in the responses of beta-carotene (y1), tocopherol (y2) and FFA (y3).It can be clearly seen that, the concentrations of the responses decrease with increasing the distilling temperature.According to the Clauses Clapeyron relation, the vapor pressure of substance increases non-linearly with temperature.Therefore, increasing temperature (Equation ( 11)) will increase mean free path which will lengthen the average distance of traveling molecule reaching the condensation board.The longer the distance of molecule travels the higher the rate of vaporization of molecule.As a result, low concentration of beta-carotene and tocopherol are found in the residue and low concentration of free fatty acid is found in the distillate.

The Effect of Temperature on Responses
Concentration of beta-carotene, tocopherol and FFA as a function of temperature are investigated at constant feed flow rate and pressure at 1500 kg/h and 5.05 × 10 −4 kPa respectively (run no. 4 and 5).In Figures 8-10, the temperature shows significant deviation in the responses of beta-carotene (y1), tocopherol (y2) and FFA (y3).It can be clearly seen that, the concentrations of the responses decrease with increasing the distilling temperature.According to the Clauses Clapeyron relation, the vapor pressure of substance increases non-linearly with temperature.Therefore, increasing temperature (Equation ( 11)) will increase mean free path which will lengthen the average distance of traveling molecule reaching the condensation board.The longer the distance of molecule travels the higher the rate of vaporization of molecule.As a result, low concentration of beta-carotene and tocopherol are found in the residue and low concentration of free fatty acid is found in the distillate.

The Effect of Temperature on Responses
Concentration of beta-carotene, tocopherol and FFA as a function of temperature are investigated at constant feed flow rate and pressure at 1500 kg/h and 5.05 × 10 −4 kPa respectively (run no. 4 and 5).In Figures 8-10, the temperature shows significant deviation in the responses of beta-carotene (y 1 ), tocopherol (y 2 ) and FFA (y 3 ).It can be clearly seen that, the concentrations of the responses decrease with increasing the distilling temperature.According to the Clauses Clapeyron relation, the vapor pressure of substance increases non-linearly with temperature.Therefore, increasing temperature (Equation ( 11)) will increase mean free path which will lengthen the average distance of traveling molecule reaching the condensation board.The longer the distance of molecule travels the higher the rate of vaporization of molecule.As a result, low concentration of beta-carotene and tocopherol are found in the residue and low concentration of free fatty acid is found in the distillate.
no. 4 and 5).In Figures 8-10, the temperature shows significant deviation in the responses of beta-carotene (y1), tocopherol (y2) and FFA (y3).It can be clearly seen that, the concentrations of the responses decrease with increasing the distilling temperature.According to the Clauses Clapeyron relation, the vapor pressure of substance increases non-linearly with temperature.Therefore, increasing temperature (Equation ( 11)) will increase mean free path which will lengthen the average distance of traveling molecule reaching the condensation board.The longer the distance of molecule travels the higher the rate of vaporization of molecule.As a result, low concentration of beta-carotene and tocopherol are found in the residue and low concentration of free fatty acid is found in the distillate.

The Effect of Pressure on Responses
The effect of pressure on responses (composition of beta-carotene, tocopherol and FFA) is investigated at feed flow 1500 kg/h and at temperature of 150 °C (run no. 1, 2, 7 and 15).The responses are found to increase with the increase of pressure, as can be seen in Figures 8-10.In agreement with the theory of mean free path for molecular distillation as shown in Equation ( 11), increasing the pressure will reduce the mean free path, which shortens the average distance of molecules reaching the condensation board.The shorter the average distance of traveling molecule the lesser the vaporization of molecule.Consequently, high concentration of beta-carotene and tocopherol are achieved at residue stream which consequently increase the free fatty acid concentration at the distillate.

√2
(11 where  is the mean free path, R is gas constant, d is molecular diameter NA is Avogadro's number

The Effect of Pressure on Responses
The effect of pressure on responses (composition of beta-carotene, tocopherol and FFA) is investigated at feed flow 1500 kg/h and at temperature of 150 °C (run no. 1, 2, 7 and 15).The responses are found to increase with the increase of pressure, as can be seen in Figures 8-10.In agreement with the theory of mean free path for molecular distillation as shown in Equation ( 11), increasing the pressure will reduce the mean free path, which shortens the average distance of molecules reaching the condensation board.The shorter the average distance of traveling molecule the lesser the vaporization of molecule.Consequently, high concentration of beta-carotene and tocopherol are achieved at residue stream which consequently increase the free fatty acid concentration at the distillate.

√2
where  is the mean free path, R is gas constant, d is molecular diameter NA is Avogadro's number

The Effect of Pressure on Responses
The effect of pressure on responses (composition of beta-carotene, tocopherol and FFA) is investigated at feed flow 1500 kg/h and at temperature of 150 • C (run no. 1, 2, 7 and 15).The responses are found to increase with the increase of pressure, as can be seen in Figures 8-10.In agreement with the theory of mean free path for molecular distillation as shown in Equation ( 11), increasing the pressure will reduce the mean free path, which shortens the average distance of molecules reaching the condensation board.The shorter the average distance of traveling molecule the lesser the vaporization of molecule.Consequently, high concentration of beta-carotene and tocopherol are achieved at residue stream which consequently increase the free fatty acid concentration at the distillate.
where λ is the mean free path, R is gas constant, d is molecular diameter N A is Avogadro's number and P is the pressure.The independent variables and responses are fitted by the following second order polynomial Equations ( 12)-( 14) which show complex relationships between the independent variables and the responses.In addition, term constants of x 1 , x 1 x 2 , and x 1 x 3 are not considered in this case due to the insignificant effect to the responses.
The R 2 value of responses, concentration of beta-carotene, tocopherol and FFA are 0.959, 0.961 and 0.951 respectively, show a very good fit.The nonlinear Equations ( 12)-( 14) now represent the MD process model which will be used for optimizing the process in the next section.
It is observed that the response surface curves of beta-carotene, tocopherol and FFA concentration display a similar relationship with respect to independent variables.The concavity of the curves shows the optimum point of each of the responses.However, those concavities are not high enough to indicate the visually maximum point of responses.

Optimization of MD for Refined Palm Oil Process
In this section the aim is to find the optimum operating parameters that will maximize the concentration of beta-carotene, tocopherol at the residue and concentration of FFA at the distillate stream.The mathematical model (Equations ( 12)-( 14)) and the objective function have been written in Equation (15).This work is more concerned with the enhancement of concentration of beta-carotene and tocopherol as 70%, than it is with reducing the concentration of free fatty acid.Therefore, in the objective function, different weighting was applied by a factor of 0.7 in y 1 , y 2 and 0.3 for y 3 .The optimization variables are feed flow rate, temperature and pressure.Lower and upper bounds for variables are specified between a valid lower and upper bound.The actual (re-scaled) upper bound and lower bound is shown in the Table 10.
Maximize: f(x) = 0.7(y 1 + y 2 ) + 0.3y 3 (15)  Solver in Microsoft Excel is an efficient tool to solve both equality and inequality constrained non-linear optimization problem [18].An Excel worksheet for optimization is prepared as shown in Figure A1.The objective function is set to target cell for maximization by changing the independent variables (refer to Figure A2).Initial value for feed flow rate, temperature and pressure were set as 1500 kg/h, 150 • C and 0.00055 kPa, respectively.The solution method used for this work is evolutionary; it is for solving a non-smooth non-linear equation.The optimum results of x 1 , x 2 and x 3 were found to be 1291 kg/h, 147 • C and 0.0007 kPa, respectively, and the corresponding responses of y 1 , y 2 and y 3 were 0.000575, 0.000937 and 0.999840, respectively.
Lastly, model validation is performed to ensure that data inserted into simulation satisfied the prediction.ASPEN HYSYS is simulated using the optimum operating conditions.The responses are shown in Table 11 and regression statistic of model and simulation is shown in Table 12.Clearly, the model equations developed from RSM can be used for predicting responses of the deodorization process.

Conclusions
Simulation and optimization of the refined palm oil process has been considered in this work.A single flash vessel model is developed in ASPEN HYSYS to represent Molecular Distillation (MD) of the refined palm oil process, which has been validated with experimental results from the literature with prediction efficiency errors of less than 3%.Response surface methodology with central composite design was then applied to design of experiments with the validated ASPEN HYSYS model to investigate the effects of feed flow rate, temperature and pressure on concentration of beta-carotene, tocopherol and free fatty acid (FFA) in the MD.The results show that response surface curves of beta-carotene, tocopherol and FFA displayed a similar relationship with respect to temperature and pressure.The result suggested that temperature and pressure and their interaction significantly influences the responses of beta-carotene, tocopherol and FFA, while feed flow rate and their interaction less significantly affected by the responses of beta-carotene, tocopherol and FFA.Finally, the process optimization through Excel solver was carried out using RSM equations, and the optimum operating parameters of MD for refined palm oil process were obtained.These parameters are feasibly suitable for refined palm oil in industrial manufacturing.

Figure 1 .
Figure 1.Comparison of distillated mass flow rate in study [10] and this work as a function of flash vessel operating temperature (°C).

Figure 2 .
Figure 2. Comparison of distillate mole fraction in study [10] and this work as a function of flash vessel operating temperature (°C).

Figure 1 . 15 Figure 1 .
Figure 1.Comparison of distillated mass flow rate in study [10] and this work as a function of flash vessel operating temperature ( • C).

Figure 2 .
Figure 2. Comparison of distillate mole fraction in study [10] and this work as a function of flash vessel operating temperature (°C).

Figure 2 .
Figure 2. Comparison of distillate mole fraction in study [10] and this work as a function of flash vessel operating temperature ( • C).
values of the dependent variable, mean or average and fitted values, respectively.

Figure 4 .
Figure 4. Carotene recovery of refining of edible oil in study[11]  and this work.

Figure 5 .
Figure 5. Response surface curve (a) and contour plot (b) of beta-carotene concentration as a function of temperature and feed flow rate.

Figure 5 .Figure 6 .Figure 7 .
Figure 5. Response surface curve (a) and contour plot (b) of beta-carotene concentration as a function of temperature and feed flow rate.Processes 2017, 5, 40 9 of 15

Figure 6 .
Figure 6.Response surface curve (a) and contour plot (b) of tocopherol concentration as a function of temperature and feed flow.

Figure 6 .Figure 7 .
Figure 6.Response surface curve (a) and contour plot (b) of tocopherol concentration as a function of temperature and feed flow.

Figure 7 .
Figure 7. Response surface curve (a) and contour plot (b) of FFA concentration as a function of temperature and feed flow.

Figure 8 . 15 Figure 8 .Figure 9 .Figure 10 .
Figure 8. Response surface curve (a) and contour plots (b) of beta-carotene concentration as a function of temperature and pressure.Processes 2017, 5, 40 10 of 15 Figure 8. Response surface curve (a) and contour plots (b) of beta-carotene concentration as a function of temperature and pressure.

Figure 9 . 15 Figure 8 .Figure 9 .Figure 10 .
Figure 9. Response surface curve (a) and contour plot (b) of tocopherol concentration as a function of temperature and pressure.

Figure 10 .
Figure 10.Response surface curve (a) and contour plot (b) of free fatty acid (FFA) concentration as a function of temperature and pressure.

Table 1 .
Results of the DISMOL simulator, ASPEN PLUS and that developed in ASPEN HYSYS.

Table 1 .
Results of the DISMOL simulator, ASPEN PLUS and that developed in ASPEN HYSYS.

Table 1 .
Results of the DISMOL simulator, ASPEN PLUS and that developed in ASPEN HYSYS.

Table 2 .
Composition of refined palm oil and its properties.

Table 2 .
Composition of refined palm oil and its properties.
and this work.

Table 6 .
Central composite circumscribed (Min and Max = Star Points) design with three independent variables.

Table 7 .
Coded variables and design range.

Table 8 .
Experimental design using Central Composite Design (CCD) of 17 runs.

Table 9 .
Estimated regression coefficient for refined palm oil.

Table 11 .
Predicted result from the Response Surface Methodology (RSM) model and ASPEN HYSYS Simulation.

Table 12 .
The regression statistics of model and simulation.