Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium in Food and Fuel Systems

*Jacqueline M. Ortega Bacicheti, Guilherme D. Machado, Fábio Nishiyama, Vladimir F. Cabral and Donato Aranda*

#### **Abstract**

Ethanol has been employed as a solvent in biodiesel production and vegetable oil refining since it is more economically attractive and less toxic than methanol and hexane. Furthermore, ethanol has demonstrated easy recovery, good selectivity, and distribution coefficient for free fatty acids (FFA), which is the primary target in the refining process since high acidity oil can lead to the formation of side products. As the knowledge of phase equilibrium behavior of fatty systems is essential to design and optimize the extraction of FFA, this chapter will present two new UNIFAC subgroups for ethanol: EtOH-B, focused on biodiesel production; and EtOH-D, focused on the deacidification process. Except for ethanol and water subgroups fitted in this study, all remaining UNIFAC parameters were taken from the literature. The new EtOH-B and EtOH-D parameters provide a considerably lower mean square error (1.20% and 0.87%) than the other works available in the literature. The results show that new ethanol subgroups and the developed methodology are valuable tools in predicting liquid-liquid phase equilibrium for ethyl biodiesel and vegetable oil deacidification systems considered, resulting in reduced computational calculations and a relatively small split with the complex dataset established by the UNIFAC-LL model.

**Keywords:** biodiesel, deacidification, ethanol, liquid-liquid equilibrium, vegetable oil, UNIFAC

#### **1. Introduction**

Several studies have introduced ethanol as a solvent for liquid-liquid extraction, resulting in a satisfactory reduction of free fatty acid (FFA) content in oils [1–4]. High acidity level vegetable oils need to be refined either for human consumption or for fuel production. In this way, ethanol can act both as a solvent for oil deacidification by liquid-liquid extraction and as a reagent for transesterification reaction.

Considering edible vegetable oil, high levels of FFA have been associated with digestive problems for consumers, autoimmune disorders, and various types of cancer [5–8]. Besides health problems, the content of pigments, phosphatides, odoriferous molecules, and free fatty acids in edible vegetal oil is detrimental to the oil quality, oxidative stability, and consumer acceptance. Among these impurities, FFA is the most detrimental because they increase the acidity of the oil, inducing an undesirable rancid flavor [9]. In order to reduce neutral oil losses and production of undesired compounds, caused by unwanted chemical reactions [10], liquid-liquid extraction using ethanol as a solvent can replace the conventional physical and chemical refining processes.

Besides deacidification of edible vegetable oils, ethanol can be employed in biodiesel production. The conventional biodiesel production process requires the use of a feedstock with reduced acidity; however, residual oils are characterized by a high FFA content, which promotes soap production with the alkali homogeneous catalysts used in transesterification [11]. Thus, feedstock purification processes, such as oil deacidification, through solvent extraction are needed. Using residual oils for biodiesel production can provide several advantages when compared to the conventional process with refined oils, such as the minimization of environmental impacts related to the disposal of waste oils, a noncompetitive relationship with the food crops, and availability of the feedstock supply, which accounts for over 70% of the refined oil biodiesel production costs [12].

In order to correctly design, operate and optimize extraction columns and the subsequent additional purification or solvent recovery units, liquid-liquid phase equilibrium (LLE) data are necessary to determine the thermodynamic properties of the systems composed of oil, solvent, and FFA [3, 13]. In general, most of the phase equilibrium data for biodiesel systems use methanol as the reacting alcohol [14], studies regarding systems with ethanol are more recent [15]. Some advantages of using ethanol as a solvent are its high affinity for FFAs at ambient temperature, low toxicity when compared to methanol, and its easy removal under mild conditions [16, 17].

Thermodynamic modeling can be used to predict or correlate the experimental behavior of fatty acid systems described through LLE data. The approach using excess Gibbs energy models is widely applied to describe biodiesel and vegetable oil systems through thermodynamic models, such as NRTL [18–22], UNIQUAC [18, 22–24], and UNIFAC [25–28]. Although thermodynamic models, such as UNIQUAC and NRTL, are capable of accurately representing the experimental data, both are restricted to the specificity of the type of vegetable oil used, as the phase equilibrium parameters are adjusted solely for the system described in the experimental data used for the LLE data fitting. Numerous types of refined or residual feedstocks are available for biodiesel production or require deacidification to be safe for human consumption. Thus, the use of a predictive group-contribution method, such as UNIFAC, to describe the phase behavior of the variety of biodiesel reaction systems or deacidification of edible vegetable oil systems is motivated, as larger number of systems would be useful.

In the UNIFAC model, alcohols, such as methanol and propanol, used to be represented by a specific group. As stated by Magnussen et al. [29] in their work, 1 propanol (CH3CH2CH2OH) can be represented, in principle, by the CH3, CH2, and OH subgroups. However, the fitting of these parameters in the UNIFAC parameter table depends on several other molecules that present these subgroups. An approach that represents molecules by a subgroup of their own, as in the case of methanol (CH3OH), can provide more accurate property calculations.

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

On account of the growing tendency of substitution of environmentally damaging chemicals for less impactful compounds, added to its availability and price in the context of biofuels in Brazil and other countries, the use of ethanol as a solvent or reactant is expected to vastly increase. Therefore, the study of the possibility of using ethanol as UNIFAC subgroup must also be evaluated.

In the literature, the thermodynamic modeling of these systems focuses on the UNIQUAC and NRTL models. Reipert et al. [30] correlated the LLE experimental data of refined babassu oil, lauric acid, and hydrated ethanol using the NRTL model. The mass fraction root mean square error (RMSD) between the observed and estimated compositions was 0.85%. Gonçalves et al. [31] applied NRTL and UNIQUAC models to correlate LLE experimental data for corn oil with a hydrated ethanol solvent system at room temperature and atmospheric pressure. The RMSD of 0.89% and 0.92%, respectively. Rodrigues et al. [32] investigated the LLE of Brazil nut oil or Macadamia nut oil + commercial oleic acid + ethanol + water, at 298.2 K. They reported RMSD not higher than 1.5% using the NRTL and UNIQUAC models. Basso et al. [19] determined LLE data for glycerol + ethanol + fatty acid ethyl ester from crambe oil system and adjusted parameters for NRTL model. The RMSD between experimental and calculated values by the NRTL model was less than 0.82%. The authors verified the prediction capacity of the UNIFAC model by testing two different sets of UNIFAC binary interaction parameters and obtained a RMSD of 2.27% and 3.97%, respectively. Da Silva et al. [20] established experimental data for Jatropha curcas oil + oleic acid + ethanol + water systems at (288.15, 298.15, 308.15, and 318.15) K and correlated the experimental data by the nonrandom two-liquid (NRTL) model with temperature-dependent binary parameters. For all systems, the RMSD was lower than 0.96%. Basso et al. [33] obtained LLE experimental data of systems containing glycerol + ethanol + ethyl biodiesel from macauba pulp oil, performed thermodynamic modeling, and simulated the settling step of this biodiesel using simulation software. Binary interaction parameters were adjusted for NRTL and UNIQUAC models and the RMSD between experimental data and calculated values were 0.44%, 1.07%, 3.52%, and 2.82%, respectively, using the NRTL, UNIQUAC, UNIFAC-LLE, and UNIFAC-Dortmund models.

All the aforementioned studies indicate that the NRTL and UNIQUAC models can reproduce the experimental liquid-liquid behavior of the systems involving vegetable oil + fatty acids + ethanol + water and biodiesel + glycerol + ethanol. However, although these thermodynamic models represent satisfactorily the experimental data, the NRTL and UNIQUAC parameters obtained in these studies are specific to these systems, consequently, they are restricted to the specific vegetable oil used in the parameter correlation. Considering the numerous types of vegetable oils that need deacidification, and different triacylglycerides matrices that can be used for biodiesel production, new parameter-fitting procedures are required for other systems, and new experimental databanks are required [9, 17, 30]. However, the new parameterfitting procedures required for new systems are not mandatory if we consider the group contribution method. UNIFAC (UNIQUAC Functional-group Activity Coefficient) model [34] is a group contribution method established by Fredenslund, Jones e Prausnitz in 1975, that can predict the liquid-liquid behavior of systems containing fatty acids using the activity coefficient calculation, which considers the interaction between the subgroups that form the molecules. Fredenslund et al. [34] proposed a group contribution method that could use experimental data available in the literature to predict the LLE of systems with no disposable experimental data.

Models based in group contribution concept estimate the properties of a mixture considering it as a solution of different functional groups that forms the molecule in the

**Figure 1.** *Soybean oil (left) and coconut oil (right) molecular structures.*

mixture. As can be seen in **Figure 1**, soybean and coconut oils are different molecules, but the same chemical compounds are encountered in these molecules. Considering these few functional groups as building blocks, a vast array of vegetable oil can be represented, as all vegetable oils present the same functional groups only differing in quantities.

This type of predictive model applicable to a larger number of systems would be useful and important in the food industry and biofuel industry, respectively, given the diversity of raw materials and the fact that most systems are multicomponent, and the variety of generated biodiesel systems.

Physical and chemical properties are considered as a sum of all contributions made by each one of the functional groups. These contributions are quantified according to the interaction parameters, adjusted by experimental data.

This model was extensively applied to study the liquid-liquid equilibrium of vegetable oils.

Batista et al. [2] applied UNIFAC and ASOG (analytical solution of groups) group contribution methods to correlate interaction parameters of triolein + oleic acid + ethanol and triolein + stearic acid + ethanol systems. The parameters were validated for canola, corn, and Spanish oil and the root mean square deviation between experimental and calculated molar fractions were 1.31% and 1.32% for the UNIFAC and ASOG models, respectively.

Bessa et al. [35] tested the predictive capability of the original UNIFAC model parameters and then modified them in terms of new readjusted binary interaction parameters. Due to inadequate results obtained by UNIFAC model without any changes in its subgroups, the authors introduced a new group ("OHgly") and two matrices of parameters were adjusted. The authors obtained good predictions and a significant improvement in the performance of this group contribution model has been achieved.

Noriega and Narvaez [28] proposed a new set of UNIFAC group interaction parameters to describe the LLE for all the systems involved in biodiesel production. The parameters presented a RMSD up to 2.07%.

In the UNIFAC model, the activity coefficient is given in terms of a combinatorial contribution taking into account entropy effects arising from differences in molecular size and shape, and a residual contribution taking into account energetic interactions between the functional groups in the mixture. To achieve a better description of the experimental data composed by the molecules presented in **Table 1**, the chosen subgroups to represent the studied fatty systems were as follows: CH3, CH2, CH, CH=CH, COOH, CH2COO, OH, H2O, and the new proposed group EtOH.

Therefore, this chapter proposes the adjustment of thermodynamic parameters of the UNIFAC model considering a new ethanol subgroup in order to predict the LLE for systems containing food or fuel. The fact that vegetable oils present different types *Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*


**Table 1.**

*Chosen subgroups to represent the studied fatty systems.*

of constituent fatty acids but with similarity in their structure makes it possible to use a predictive approach that represents the interactions between groups based on UNIFAC thermodynamic model in the proper representation of the LLE for industrial separation processes. Thus, this study aims to present such methodology in the prediction of food systems containing FFA to be separated as well as systems containing fuels, focusing on biodiesel produced from vegetable oils.

#### **2. Methodology**

The isofugacity criterion for phase equilibrium is conventionally used to describe a condition at which the chemical potential of each component is equal in the phases among which this component can distribute. If fugacities are expressed in terms of activity coefficients, the LLE using an excess Gibbs energy model is represented by Eq. (1) [34].

$$\mathbf{x}\_i^{\text{AP}} \mathbf{y}\_i^{\text{AP}} = \mathbf{x}\_i^{\text{OP}} \mathbf{y}\_i^{\text{OP}} \tag{1}$$

where *xAP <sup>i</sup>* and *xOP <sup>i</sup>* are the molar fractions of component *i* in the alcohol and oil phases, respectively, and *γAP <sup>i</sup>* and *γOP <sup>i</sup>* are the activity coefficients of component *i* in the alcohol and oil phases, respectively.

In this chapter, the activity coefficients are calculated through UNIFAC model [34], which considers the combinatorial and residual contributions, Eq. (2).

$$\ln \gamma\_i = \ln \gamma\_i^C + \ln \gamma\_i^R \tag{2}$$

The combinatorial contribution *γ<sup>C</sup> <sup>i</sup>* (Eq. (3)) is related to the difference in size and shape of molecules. The volume fraction ∅*<sup>i</sup>* and surface fraction *θ<sup>i</sup>* of each molecule *i* are obtained using Eq. (5) and (6), respectively.

$$\ln \chi\_i^{\mathbb{C}}(\mathbf{x}) = \ln \frac{\mathfrak{Q}\_i}{\mathfrak{x}\_i} + \frac{Z}{2} q\_i \ln \frac{\theta\_i}{\mathfrak{Q}\_i} + l\_i - \frac{\mathfrak{Q}\_i}{\mathfrak{x}\_i} \sum\_{j=1}^{n\_c} \mathfrak{x}\_j l\_j,\tag{3}$$

$$l\_i = \frac{z}{2} \left( r\_i - q\_i \right) - (r\_i - 1); Z = 10. \tag{4}$$

*Ethanol and Glycerol Chemistry – Production, Modelling, Applications, and Technological Aspects*

$$\mathcal{Q}\_{i} = \frac{r\_{i}\mathbf{x}\_{i}}{\sum\_{j=1}^{n\_{c}} r\_{j}\mathbf{x}\_{j}} \tag{5}$$

$$\theta\_i = \frac{q\_i \mathbf{x}\_i}{\sum\_{j=1}^{n\_c} q\_j \mathbf{x}\_j} \tag{6}$$

where *xi* is the mole fraction of component *i* and *ri* and *qi* are the measures of molecular van der Waals volumes and molecular surface areas, respectively. They are calculated as the sum of the group volume and group area parameters, *Rk* and *Qk*.

$$r\_i = \sum\_{k=1}^{n\_\sharp} \nu\_{ki} R\_k \tag{7}$$

$$q\_i = \sum\_{k=1}^{n\_g} \nu\_{ki} \mathbf{Q}\_k \tag{8}$$

where *νki* is the number of groups of type *k* in molecule *i*. Here, we use the values of *Rk and Qk* reported by Magnussen *et al.* [29]. Os valores *Rk* e *Qk* são calculados a partir do volume e área superficial dos grupos de Van der Waals (*Vwk* e *Awk*, respectivamente), retratado por Bondi (1968):

$$R\_k = \frac{V\_{wk}}{15,17} \tag{9}$$

$$Q\_k = \frac{A\_{wk}}{2.5 \times 10^9} \tag{10}$$

The residual contribution *γ<sup>R</sup> <sup>i</sup>* , described through Eq. (11), is due to group areas and group interaction parameters. Γ*<sup>k</sup>* and Γð Þ*<sup>i</sup> <sup>k</sup>* are the residual activity coefficient from group k in the solution and in a solution containing only molecules of type *i*, respectively.

$$\ln \gamma\_i^R(T, \varkappa) = \sum\_k \nu\_{ki} \left( \ln \Gamma\_k - \ln \Gamma\_k^{(i)} \right) \tag{11}$$

By Eq. 12 it is possible to calculate the residual activity coefficient.

$$\ln \Gamma\_k = Q\_k \left[ 1 - \ln \left( \sum\_m^k \theta\_m \boldsymbol{\nu}\_{mk} \right) - \sum\_m \left( \frac{\theta\_{dm} \boldsymbol{\nu}\_{km}}{\sum\_n \theta\_n \boldsymbol{\nu}\_{nm}} \right) \right] \tag{12}$$

where *θ<sup>m</sup>* is the area fraction of group m given by:

$$\theta\_m = \frac{Q\_m X\_m}{\sum\_n Q\_n X\_n},\# \tag{13}$$

and *Xm* is the molar fraction of the group m in the mixture.

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

$$X\_m = \frac{\sum\_j^n Q\_m X\_m}{\sum\_j \sum\_n \nu\_{nj} \alpha\_j} \neq \tag{14}$$

The group interaction parameter *ψnm* between groups n and m is given by:

$$\psi\_{mn} = \exp\left(-\frac{a\_{mn}}{T}\right), \neq \tag{15}$$

where *amn* is an adjustable parameter of the binary interaction between groups *m* and *n*, and it has units of kelvins. Each group-group interaction possesses two parameters: *anm*, *amn*, and *anm* 6¼ *amn*. These parameters were obtained from a database using a wide range of experimental phase-equilibrium results.

The parameter fitting procedure was applied separately to fit the EtOH interaction parameters for biodiesel separation and deacidification systems. In deacidification process, the presence of water in ethanol solvent increases the polarity of the solution. This feature and the presence of glycerol in biodiesel system makes the behavior of EtOH subgroup different in deacidification systems from that in biodiesel systems, even assuming that the subgroups of both systems were similar (the CH3, CH2, CH, CH ¼ CH, COOH, CH2COO), and accepting that UNIFAC model is a group contribution method that can predict the liquid-liquid behavior considering the interaction between the subgroups that form the molecules.

As the ethanol parameters calculated for the biodiesel system are not adequate for deacidification process, the UNIFAC parameters for new ethanol subgroup were fitted for the biodiesel process separately from deacidification process. As result, this study will present UNIFAC interaction parameters for ethanol subgroup in biodiesel systems (EtOH-B) and in deacidification systems (EtOH-D). The UNIFAC parameters were fitted considering the phase compositions in molar fractions and binary interaction parameters in 1/K.

A data bank containing 56 systems was compiled. We used a total of 88 tie-lines from eight biodiesel types from vegetable oils with ethanol, and 246 tie-lines for the 14 types of vegetable oils with hydrated ethanol, at temperatures ranging from 288.15 to 333.15 K, all systems studied at atmospheric pressure. **Table 2** shows a summary of the equilibrium systems used in the parameter fitting and validation procedure.

For biodiesel separation fitted ethanol-related parameters, the LLE database of seven ethyl systems of vegetable oils (soybean, canola, palm, *jatropha curcas*, cottonseed, crambe, and sunflower) at different temperatures were used to fit the interaction parameters of the ethanol subgroups. The fitted ethanol-related parameters were validated against data for the macauba (*Acrocomia aculeata*) biodiesel system. Lastly, the results obtained are compared with those obtained by the UNIFAC-LL parameters [29] and parameters available from literature work [35].

For deacidification fitted ethanol-related parameters, the LLE database of 13 systems of vegetable oils (brazil nut, corn, cottonseed, garlic, grapeseed, *jatropha curcas*, palm, peanut, rice bran, sesame, soybean, and sunflower) at different temperatures were used to fit the interaction parameters of the ethanol subgroups. The validation procedure was performed using canola oil system and a total of five tie-lines at 303.15 K, and it is worth mentioning that canola oil was not included in the parameter estimation process. Lastly, the results obtained are compared with those obtained by the UNIFAC-LL parameters [29] and parameters available from literature work [26, 28, 43].


#### **Table 2.**

*LLE database of the systems involved in the correlation process.*

The software Microsoft Excel was used in the parameter fitting procedure coupled to XSEOS [44] and SOLVER® add-ins. The XSEOS add-in, an open-source code programmed in visual basic for applications (VBA) with several excess Gibbs energy models and equations of state, was employed to evaluate the activity coefficient of the UNIFAC model, while the SOLVER® add-in with the generalized reduced gradient (GRG) method [45] was used as the numerical calculation tool.

The parameter fitting was performed by minimizing the objective function (Eq. 16) using the generalized reduced gradient (GRG) nonlinear solving method.

$$\text{OF} = \sum\_{n}^{\text{NLa}} \sum\_{i}^{\text{NCr}} \left( \mathbf{W}\_{i.n}^{\text{OPexp}} - \mathbf{W}\_{i.n}^{\text{OPcalc}} \right)^{2} + \left( \mathbf{W}\_{i.n}^{\text{APexp}} - \mathbf{W}\_{i.n}^{\text{APcalc}} \right)^{2} \tag{16}$$

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

where *NLa* and *NCn* represent the total number of tie-lines in each group and the total number of components or pseudo components in tie-line *n*, respectively. *Wi*,*<sup>n</sup>* correspond to the mass fraction of the component or pseudo component *i* in the tieline *n* in the oil phase (*OP*) or alcohol phase (*AP*), while the superscripts *exp* and *calc* the experimental or calculated values, respectively. During the parameter fitting procedure, the minimization of eq. (16) is constrained by the isofugacity condition expressed by Eq. (1), and the calculated mass fraction sum, which must be equal to 1.

In order to evaluate the capability of the ethanol interaction parameters [27, 43] to describe LLE in systems containing FFA accurately, the percentual mass fraction root mean square deviation between the experimental and calculated phase composition values RMSD ð Þ, Eq. (17), was compared to the results obtained using parameters available in the literature (UNIFAC-LL from Magnussen et al. [46], Bessa et al. [35], Noriega and Narváez [28], Hirata et al. [26]).

$$\text{RMSD} = 100. \sqrt{\frac{\sum\_{n}^{\text{Nla}} \sum\_{i}^{\text{NCla}} \left(W\_{i,n}^{\text{OPexp}} - W\_{i,n}^{\text{OPcalc}}\right)^2 + \left(W\_{i,n}^{\text{APexp}} - W\_{i,n}^{\text{APcalc}}\right)^2}{2 \text{NLA} \cdot \text{NCla}} \text{\*} \tag{17}$$

The subgroups chosen to represent the fatty acids systems were as follows: CH3, CH2, CH, CH ¼ CH, COOH, CH2COO, H2O, and ethanol, EtOH. Only the binary interaction parameters corresponding to water and ethanol were estimated in this study.

In order to reduce the number of components considered in the interaction parameter estimation methodology, the pseudo-component approach can be used to characterize a complex mixture as a single component [9]. Several authors had successfully applied this approach assuming the vegetable oil as a single pseudo component with an average molar mass and average physical-chemical properties [4, 9, 17, 19, 20, 26, 28, 30–32, 35, 39–43].

Based on the results of the reported studies, it is assumed that the use of the described methodology does not cause expressive deviations in the thermodynamic modeling of deacidification and biodiesel separation systems. Therefore, each vegetable oil was replaced by a pseudo component, and the same approach was used for commercial fatty acids. The pseudo component is a thermodynamic tool applied to represent the edible oil as a single TAG and FFAs as a unique FFA. A weighted average of the vegetable oil and FFA molar masses and subgroup numbers was applied to each pseudo component, considering the fatty acid profile of each vegetable oil, [2–4, 20, 31, 32, 39–42]. The molar mass data of the components were obtained from the NIST Chemistry WebBook database.

#### **3. Results and discussion**

For biodiesel separation system, the methodology was initially applied to fit parameters for the new interaction subgroup representing the ethanol molecule (EtOH-B), while the UNIFAC-LL for the subgroups forming the other components of the system (ethyl biodiesel and glycerol) remained unchanged. The fitting process used seven ethyl biodiesel from soybean, canola, palm, *jatropha curcas*, cottonseed, crambe, and sunflower oils available in the literature. **Table 3** summarizes the UNIFAC-LL parameters [29] used and the interaction parameters fitted for the new EtOH-B subgroup proposed [43].


#### **Table 3.**

*UNIFAC interaction parameters for* CH3*,* CH2*,* CH*,* CH ¼ CH*,* COOH*,* CH2COO*,* H2O*, and EtOH subgroups.*

The bold numbers are the fitted one.

**Figures 2** and **3** show the experimental points and calculated tie-lines obtained using the new interaction parameters summarized in **Table 3** for different types of biodiesel. For all diagrams presented in this chapter, the following classification is valid:

<sup>Δ</sup> and ▲ for experimental and calculated points, respectively;

□ and ■ for experimental and calculated overall composition, respectively;



The tie-lines overlapping seen in **Figures 2** and **3** suggest that the fitted UNIFAC parameters for the proposed EtOH subgroup can predict the behavior of the systems considered with high accuracy.

In order to validate these new fitted EtOH parameters, we perform a liquid-liquid equilibrium prediction for ethyl biodiesel from macauba pulp, which was not used during the parameter fitting process.

#### **Figure 2.**

*Ternary LLE diagram for the ethyl biodiesel from soybean oil system at 293.15 K (left) and* jatropha curcas *oil at 318.15 K (right). Experimental data from [14, 36], respectively.*

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

**Figure 3.**

*Ternary LLE diagram for the ethyl biodiesel from cottonseed oil system at 333.15 K (left) and crambe oil at 298.2 K (right). Experimental data from [19, 38], respectively.*

**Figure 4** shows the results of this prediction in the form of a ternary diagram, which exhibit small deviations from the experimental data, thus validating the methodology used.

The results obtained in the fitting and validation processes were then compared with UNIFAC parameters proposed by Magnussen et al. [29] and Bessa et al. [35]. The UNIFAC-LL database fitted by Magnussen et al. [29] in 1981 is extensively widespread and applied to describe fatty systems equilibrium [47–49]. A more recent research by Bessa and collaborators [35] refitted all interaction parameters of the UNIFAC-LL and proposed a new OH subgroup used to represent uniquely this subgroup present in the glycerol molecule, thus having to fit 42 interaction parameters.

**Table 4** shows the percentage mean square error (MSE%). The results using the parameters proposed by Machado et al. [43] are always better than those using the Magnussen et al. [29] and Bessa et al. [35] parameters.

**Figure 5** shows ternary LLE diagrams for the ethyl biodiesel from macauba pulp oil system comparison with Bessa et al. [35] fitted parameters (left) and with UNIFAC-LL parameters (right).

#### **Figure 4.**

*Ternary LLE diagram for the ethyl biodiesel from macauba pulp oil system at 298.15 K for validation procedure. Experimental data from [33].*


#### **Table 4.**

*Comparative percentage mean square error (MSE%) for this study's fitting and the literature (EtOH).*

#### **Figure 5.**

*Ternary LLE diagram for the ethyl biodiesel from macauba pulp oil system at 298.15 K comparison with the Bessa et al. [35] (left, red tie-lines); and with the UNIFAC-LL [29] data (right, green tie-lines). Experimental data from ref. [19].*

**Figure 5** (right) shows that the UNIFAC-LL results have large deviations and the use of such parameters to predict equilibrium between the glycerol-rich phase and the biodiesel-rich phase for the systems considered by Machado et al. [43] is not

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

recommended. Although the parameters fitted by Bessa et al. [35] could better describe the LLE than the UNIFAC-LL [29] parameters, still there is considerable discrepancy between the predictions and the experimental data. Such observations corroborate the results shown in **Table 4**. **Figure 4** shows that ethanol subgroup fitting tie-lines represent considerably better the equilibrium of the biodiesel-rich phase with the glycerol-rich phase. It is worth mentioning that, only eight parameters were fitted, while Bessa et al. [35] considered 42 adjustable parameters in their study.

For deacidification system, the methodology was initially applied to fit parameters for the new interaction subgroup representing the ethanol molecule (EtOH-D) and H2O subgroups, while the UNIFAC-LL for the subgroups forming the other components of the system (vegetable oil and free fatty acids—FFA) remained unchanged. The fitting process used experimental LLE data of 17 (Brazil nut, corn, cottonseed, garlic, grapeseed, *jatropha curcas*, macadamia, palm, peanut, rice bran, sesame, soybean, and sunflower) different vegetable oils available in the literature. **Table 5** shows the original UNIFAC-LL [29] parameters used and the interaction parameters fitted for the water (H2O) and EtOH-D subgroups.

**Figures 6** and **7** show the experimental points and calculated tie-lines obtained using the new interaction parameters summarized in **Table 5** for corn and soybean oil. Due to the systems presenting different water contents added to the ethanol solvent, ethanol and water were used as a mixed solvent to represent the pseudo-quaternary systems in a triangular diagram.

Analyzing **Figures 6** and **7** it can be noted an inversion in the tie-line slopes as water is added to the solvent. This phenomenon occurs because water decreases the solubility between oil and ethanol; hence, it increases the two-liquid phase regions [30, 32, 39, 42].

The ternary diagrams presented in **Figures 6** and **7** show that the calculated results are very close to the experimental data. Therefore, the H2O and EtOH-D adjusted interaction parameters using the UNIFAC model correlated with the high accuracy of the LLE behavior of the considered system containing fatty acids.

**Table 6** compares the RMSD from the experimental mass fraction data from the calculated data obtained by UNIFAC modeling using the parameters fitted by Bacicheti et al. [27] for EtOH-D subgroup with those obtained using parameters available in the literature (UNIFAC-LL from Magnussen et al. [29]), Noriega and Narváez [28], Hirata et al. [26]). **Table 6** still presents the RMSD between experimental and calculated data obtained using Machado et al. [43] parameter set.


**Table 5.**

*UNIFAC-LL [29] parameters and 18 interaction parameters fitted for* H2O *and EtOH-D subgroups.*

**Figure 6.**

*LLE diagram for corn oil + commercial oleic acid 1 + and solvent (ethanol + %water) system at 298.15 K [31]. All compositions are on a mass basis.*

**Figure 7.**

*LLE diagram for soybean oil + commercial linoleic acid 1 + and solvent (ethanol + water) system at 323.15 K [42]. All compositions are on a mass basis.*


*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

#### **Table 6.**

*Root mean square deviation (*RMSD *for fittings in this chapter and the literature.*

As commented previously, Noriega and Narvaez [28] fitted 30 UNIFAC parameters of the subgroups related to pseudo-ternary and pseudo-quaternary systems of biodiesel + alcohol + glycerol and oil + fatty acid + alcohol + water. Hirata et al. [26] used plenty of pseudo-quaternary data available in the literature to fit all the 30 UNIFAC interaction parameters of interest for fatty systems. The present study adjusted only 16 interaction parameters for pseudo-ternary biodiesel systems and adjusted only 18 interaction parameters for pseudo-quaternary deacidification systems.

The global RMSD displayed in **Table 6** shows that the parameters set proposed by UNIFAC-LL [29] and Noriega and Narváez [28] are not suitable to describe the liquidliquid equilibrium of deacidification and biodiesel-related systems. The global RMSD using UNIFAC-LL and Noriega and Narváez [28] were 9.24 and 13.43, respectively, while using the parameters fitted by Bacicheti et al. [27] resulted in a global RSMD of 0.87.

The considerable difference in RMSD between experimental and calculated data applying Machado et al. [43] parameter set (10.81%) and Bacicheti et al. [27] parameter set (0.87%) subgroup parameters emphasize the difference in ethanol subgroup in biodiesel separation and deacidification systems due to water and glycerol molecules in the system.

It was expected that Noriega and Narváez [28] parameters and Machado et al. [43] parameter set resulted in slightly high deviations (**Table 6**) since they were adjusted to biodiesel systems. But we were expecting that the other UNIFAC subgroups (CH3, CH2, CH, CH ¼ CH, COOH, CH2COO) would be able to represent the liquid-liquid equilibrium behavior, giving a fair and enough performance of UNIFAC model. However, the presence of water is too sharp, turning the hole systems distinct and the deviations significant, despite having similar UNIFAC subgroups.

In contrast to results obtained by the other UNIFAC parameter matrices, the RMSD obtained by using Hirata et al. [26] proposed parameters were closer to this study's parameters performance. It is worthy to quote that, although using a similar databank, this chapter ensues in lower RMSD between experimental and calculated data and, moreover, it has the upgrade of adjusting only 18 interaction parameters, resulting in a considerable faster procedure and relatively lower computational calculations than Hirata et al. [26] that fitted 30 interaction parameters.

The validation procedure of the parameters was taken for canola oil system, which was not used in the parameter fitting process. **Table 7** lists the RMSD between the experimental and calculated molar fractions using Bacicheti et al. [27] parameter set with those obtained using UNIFAC-LL from Magnussen et al. [29]), Noriega and Narváez [28], and Hirata et al. [26]). As previously stated, canola oil [18] data did not take part in this study's parameter fitting, but it was used by Hirata et al. [26] in their parameter fitting procedure, resulting in alike RMSD.

According to **Table 7**, the Bacicheti et al. [27] parameter set exhibits lower RMSD values than those obtained by applying the parameter matrix from the previously cited authors; hence the EtOH-D subgroup can successfully describe the deacidification systems of vegetable oils.

**Figure 8** exhibits the ternary diagrams for the validation with canola oil data. Bacicheti et al. [27] parameter set exhibited small deviations from the experimental data than the other parameter matrices. In contrast, the results using the Machado et al. [43] parameter set and the parameters of Noriega and Narváez [28] have a considerable deviation from the experimental data, and thus, they are not recommended for quantitatively predicting the equilibrium of deacidification of vegetable oils.

As shown in **Figure 8**, Machado et al. [43] and Noriega and Narváez [28] parameter set, and UNIFAC-LL [29] present the highest average deviation from the experimental points. The average deviation values between the experimental and calculated data obtained using the interaction parameters from Hirata et al. [26] were close to


**Table 7.**

*RMSD between experimental and calculated data mass fraction for the validation system composed of canola oil + commercial oleic acid 3 + ethanol + water at 303.15 K.*

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

**Figure 8.**

*LLE diagram for validation procedure with canola oil system [18]. Comparison between literature data (traced line) with predictions using parameters from Bacicheti et al. [27] (black lines), Noriega and Narváez [28] (yellow lines), Machado et al. [43] (blue lines), UNIFAC-LL [29] (purple lines), and Hirata et al. [26] (green lines). All compositions are on a mass basis.*

those obtained by Bacicheti et al. [27]. However, Bacicheti et al. [27] only refitted 18 parameters, taking fewer computational calculations than those required by Hirata et al. [26], who readjusted all of them.

#### **4. Conclusions**

The new ethanol subgroups for the UNIFAC model presented in this chapter could correctly and accurately describe the liquid-liquid phase equilibrium of deacidification and biodiesel systems. The ethanol binary interaction parameters were fitted using a data bank of LLE available in the literature for biodiesel vegetable oil systems.

Using relatively few computational calculations and aiming to result in a relatively small split with the complex data set established by the UNIFAC-LL model, Bacicheti et al. [27] maintain CH3, CH2, CH, CH ¼ CH, COOH and CH2COO subgroups parameters, and refitted just water and ethanol UNIFAC interaction parameters. Ethanol-fitted parameters were then validated with very small deviations in the prediction of macauba pulp oil biodiesel for EtOH-B and canola oil for EtOH-D. The overall deviation calculated in the validation of these new UNIFAC parameters was 1.20% for biodiesel systems and 0.87% for deacidification systems.

When compared to prediction results for biodiesel systems using UNIFAC-LL parameters and those fitted by Bessa et al. [35], Machado et al. [43] had better results for the tie-lines, despite the small number of parameters fitted. The same results were obtained considering deacidification systems of vegetable oils, which ethanol interaction parameters proposed by Bacicheti et al. [27] exhibited small deviations from the experimental data and lower deviations than Noriega and Narváez [28], UNIFAC-LL [29], and Hirata et al. [26].

The method applied can contribute to a better description of the phase behavior of fatty systems involved in the deacidification of vegetable oil using liquid-liquid extraction and biodiesel separation process, as this biofuel is subject to strict composition and purity regulations. Moreover, taking into account the lower root mean square deviation between the experimental and calculated molar fractions obtained here and the relatively low computational calculations for the parameter fitting procedure, this chapter presents an easier and faster alternative approach instead to refitting all UNIFAC parameters.

The results show that the methodology employed is consistent and may be useful in predicting equilibrium when experimental equilibrium data are not available.

#### **Funding**

This work was supported by the Brazilian agencies CAPES [88887.486092/2020– 2100] and CNPq [grant number 421930/2016–0].

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

### **Author details**

Jacqueline M. Ortega Bacicheti1,2, Guilherme D. Machado2,3\*, Fábio Nishiyama2,3, Vladimir F. Cabral1,2 and Donato Aranda<sup>3</sup>

1 Food Engineering Program, State University of Maringá, Maringá, PR, Brazil

2 Chemical Engineering Program, State University of Maringá, Maringá, PR, Brazil

3 Department of Chemical Engineering, Federal Tecnologic University of Paraná-UTFPR, Londrina, PR, Brazil

\*Address all correspondence to: guilhermed@utfpr.edu.br

© 2023 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### **References**

[1] Sérgio MBO, Soletti BJI, Carvalho SHV, Queimada AJ, Coutinho JAP. Liquid– liquid equilibria for the canola oil biodiesel +ethanol+glycerol system. Fuel. 2011;**90**: 2738-2745. DOI: 10.1016/j.fuel.2011.03.017

[2] Batista E, Monnerat S, Stragevitch L, Pina CG, Gonçalves CB, Meirelles AJA. Prediction of liquidliquid equilibrium for Systems of Vegetable Oils, fatty acids, and ethanol. Journal of Chemical & Engineering Data. 1999;**44**:1365-1369. DOI: 10.1021/je9900169

[3] Cuevas MS, Rodrigues CEC, Gomes GB, Meirelles AJA. Vegetable oils Deacidification by solvent extraction: Liquidliquid equilibrium data for systems containing sunflower seed oil at 298.2 K. Journal of Chemical & Engineering Data. 2010;**55**:3859-3862. DOI: 10.1021/je900791w

[4] Rodrigues CEC, Reipert ÉCD, de Souza AF, Filho PAP, Meirelles AJA. Equilibrium data for systems composed by cottonseed oil+commercial linoleic acid+ethanol+water+tocopherols at 298.2K. Fluid Phase Equilibria. 2005;**238**: 193-203. DOI: 10.1016/j.fluid.2005. 09.027

[5] Tripathy D, Mohanty P, Dhindsa S, Syed T, Ghanim H, Aljada A, et al. Elevation of free fatty acids induces inflammation and impairs vascular reactivity in healthy subjects. Diabetes. 2003;**52**:2882-2887. DOI: 10.2337/ diabetes.52.12.2882

[6] Röhrig F, Schulze A. The multifaceted roles of fatty acid synthesis in cancer. Nature Reviews. Cancer. 2016;**16**: 732-749. DOI: 10.1038/nrc.2016.89

[7] Scoville EA, Allaman MM, Adams DW, Motley AK, Peyton SC, Ferguson SL, et al. Serum

polyunsaturated fatty acids correlate with serum cytokines and clinical disease activity in Crohn's disease. Scientific Reports. 2019;**9**:2882. DOI: 10.1038/ s41598-019-39232-z

[8] Tsoukalas D, Fragoulakis V, Sarandi E, Docea AO, Papakonstaninou E, Tsilimidos G, et al. Targeted Metabolomic analysis of serum fatty acids for the prediction of autoimmune diseases. Frontiers in Molecular Biosciences. 2019;**6**:120. DOI: 10.3389/fmolb.2019.00120

[9] Rodrigues CEC, Meirelles AJA. Extraction of free fatty acids from Peanut oil and avocado seed oil: Liquidliquid equilibrium data at 298.2 K. Journal of Chemical & Engineering Data. 2008;**53**:1698-1704. DOI: 10.1021/ je7007186

[10] Rodrigues CEC, Gonçalves CB, Marcon EC, Batista EAC, Meirelles AJA. Deacidification of rice bran oil by liquid– liquid extraction using a renewable solvent. Separation and Purification Technology. 2014;**132**:84-92. DOI: 10.1016/j.seppur.2014.05.009

[11] Lu W, Xu Z, Li M, Ma Y, Xiao Z. Deep eutectic solvent (DES) assisted deacidification of acidic oil and retaining catalyst activity: Variables optimization and catalyst characterization. Industrial Crops and Products. 2022;**183**:114990. DOI: 10.1016/j.indcrop.2022.114990

[12] Hajjari M, Tabatabaei M, Aghbashlo M, Ghanavati H. A review on the prospects of sustainable biodiesel production: A global scenario with an emphasis on waste-oil biodiesel utilization. Renewable and Sustainable Energy Reviews. 2017;**72**:445-464. DOI: 10.1016/j.rser.2017.01.034

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

[13] Mu T, Rarey J, Gmehling J. Group contribution prediction of surface charge density profiles for COSMO-RS(Ol). AIChE Journal. 2007;**53**:3231-3240. DOI: 10.1002/aic.11338

[14] Silva JRF, Mazutti MA, Voll FAP, Cardozo-Filho L, Corazza ML, Lanza M, et al. Thermophysical properties of biodiesel and related systems: (liquid +liquid) equilibrium data for Jatropha curcas biodiesel. The Journal of Chemical Thermodynamics. 2013;**58**:467-475. DOI: 10.1016/j.jct.2012.10.006

[15] Mazutti MA, Voll FAP, Cardozo-Filho L, Corazza ML, Lanza M, Priamo WL, et al. Thermophysical properties of biodiesel and related systems: (liquid+liquid) equilibrium data for soybean biodiesel. The Journal of Chemical Thermodynamics. 2013;**58**: 83-94. DOI: 10.1016/j.jct.2012.10.026

[16] Bou Orm R, Citeau M, Comitis A, Savoire R, Harscoat-Schiavo C, Subra-Paternault P, et al. Walnut oil deacidification by liquid–liquid extraction with ethanol in a single- and multistage crossflow process. OCL. 2020;**27**:35. DOI: 10.1051/ocl/2020029

[17] Mohsen-Nia M, Modarress H, Nabavi HR. Measuring and modeling liquid–liquid equilibria for a soybean oil, oleic acid, ethanol, and water system. Journal of the American Oil Chemists' Society. 2008;**85**:973-978. DOI: 10.1007/ s11746-008-1288-9

[18] Batista E, Monnerat S, Kato K, Stragevitch L, Meirelles AJA. Liquidliquid equilibrium for Systems of Canola oil, oleic acid, and short-chain alcohols. Journal of Chemical & Engineering Data. 1999;**44**:1360-1364. DOI: 10.1021/je990015g

[19] Basso RC, de Almeida Meirelles AJ, Batista EAC. Liquid–liquid equilibrium

of pseudoternary systems containing glycerol+ethanol+ethylic biodiesel from crambe oil (Crambe abyssinica) at T/K= (298.2, 318.2, 338.2) and thermodynamic modeling. Fluid Phase Equilibria. 2012; **333**:55-62. DOI: 10.1016/j.fluid.2012. 07.018

[20] da Silva CAS, Sanaiotti G, Lanza M, Meirelles AJA, Batista EAC. Liquidliquid equilibrium data for systems containing Jatropha curcas oil + oleic acid + anhydrous ethanol + water at (288.15 to 318.15) K. Journal of Chemical & Engineering Data. 2010;**55**:2416-2423. DOI: 10.1021/je900831p

[21] da Silva AE, Lanza M, Batista EAC, Rodrigues AMC, Meirelles AJA, da Silva LHM. Liquidliquid equilibrium data for systems containing palm oil fractions + fatty acids + ethanol + water. Journal of Chemical & Engineering Data. 2011;**56**:1892-1898. DOI: 10.1021/je100 9015

[22] Mohsen-Nia M, Dargahi M. Liquidliquid equilibrium for systems of (corn oil + oleic acid + methanol or ethanol) at (303.15 and 313.15) K, Journal of Chemical & Engineering Data 52 (2007) 910–914. doi:10.1021/ je700051j

[23] Serres JDS, Soares D, Corazza ML, Krieger N, Mitchell DA. Liquid–liquid equilibrium data and thermodynamic modeling for systems related to the production of ethyl esters of fatty acids from soybean soapstock acid oil. Fuel. 2015;**147**:147-154. DOI: 10.1016/j. fuel.2015.01.059

[24] Mohsen-Nia M, Khodayari A. Deacidification of sunflower oil by solvent extraction: (liquid+liquid) equilibrium data at T=(303.15 and 313.15) K. The Journal of Chemical Thermodynamics. 2008;**40**:1325-1329. DOI: 10.1016/ j.jct.2008.01.029

[25] Ferreira MC, Bessa LCBA, Abreu CRA, Meirelles AJA, Caldas Batista EA. Liquid-liquid equilibrium of systems containing triolein + (fatty acid/ partial acylglycerols/ester) + ethanol: Experimental data and UNIFAC modeling. Fluid Phase Equilibria. 2018; **476**:186-192. DOI: 10.1016/j.fluid.2018. 07.030

[26] Hirata GF, Abreu CRA, Bessa LCBA, Ferreira MC, Batista EAC, Meirelles AJA. Liquid–liquid equilibrium of fatty systems: A new approach for adjusting UNIFAC interaction parameters. Fluid Phase Equilibria. 2013;**360**:379-391. DOI: 10.1016/j.fluid.2013.10.004

[27] Bacicheti JMO, Machado GD, Cabral VF. Liquid–liquid equilibrium calculations of systems containing vegetable oil + fatty acids + ethanol + water using new parameters for the unifac subgroups of ethanol and water. Fluid Phase Equilibria. 2021; **548**:113182. DOI: 10.1016/j.fluid.2021. 113182

[28] Noriega MA, Narváez PC. UNIFAC correlated parameters for liquid-liquid equilibrium prediction of ternary systems related to biodiesel production process. Fuel. 2019;**249**:365-378. DOI: 10.1016/j.fuel.2019.03.124

[29] Magnussen T, Rasmussen P, Fredenslund A. UNIFAC parameter table for prediction of liquid-liquid equilibriums. Industrial & Engineering Chemistry Process Design and Development. 1981;**20**:331-339. DOI: 10.1021/i200013a024

[30] Reipert ÉCD, Rodrigues CEC, Meirelles AJA. Phase equilibria study of systems composed of refined babassu oil, lauric acid, ethanol, and water at 303.2K. The Journal of Chemical Thermodynamics. 2011;**43**:1784-1790. DOI: 10.1016/j.jct.2011.05.039

[31] Gonçalves CB, Batista E, Meirelles AJA. Liquidliquid equilibrium data for the system corn oil + oleic acid + ethanol + water at 298.15 K. Journal of Chemical & Engineering Data. 2002;**47**: 416-420. DOI: 10.1021/je010273p

[32] Rodrigues CEC, Silva FA, Marsaioli A, Meirelles AJA. Deacidification of Brazil nut and macadamia nut oils by solvent extraction: Liquidliquid equilibrium data at 298.2 K. Journal of Chemical & Engineering Data. 2005;**50**:517-523. DOI: 10.1021/je049687j

[33] Basso RC, da Silva CAS, de Sousa CO, de Meirelles AJA, Batista EAC. LLE experimental data, thermodynamic modeling and sensitivity analysis in the ethyl biodiesel from macauba pulp oil settling step. Bioresource Technology. 2013;**131**:468-475. DOI: 10.1016/j.biortech. 2012.12.190

[34] Fredenslund A, Jones RL, Prausnitz JM. Group-contribution estimation of activity coefficients in nonideal liquid mixtures. AICHE Journal. 1975;**21**:1086-1099. DOI: 10.1002/aic.690210607

[35] Bessa LCBA, Ferreira MC, Abreu CRA, Batista EAC, Meirelles AJA. A new UNIFAC parameterization for the prediction of liquid-liquid equilibrium of biodiesel systems. Fluid Phase Equilibria. 2016;**425**:98-107. DOI: 10.1016/j.fluid. 2016.05.020

[36] Mesquita FMR, Feitosa FX, Sombra NE, de Santiago-Aguiar RS, de Sant'Ana HB. Liquid–liquid equilibrium for ternary mixtures of biodiesel (soybean or sunflower) + glycerol + ethanol at different temperatures. Journal of Chemical & Engineering Data. 2011;**56**:4061-4067. DOI: 10.1021/ je200340x

*Ethanol as a Subgroup of the UNIFAC Model in the Prediction of Liquid-Liquid Equilibrium… DOI: http://dx.doi.org/10.5772/intechopen.108420*

[37] de Rocha EGA, Follegatti-Romero LA, Duvoisin S, Aznar M. Liquid–liquid equilibria for ternary systems containing ethylic palm oil biodiesel+ethanol +glycerol/water: Experimental data at 298.15 and 323.15K and thermodynamic modeling. Fuel. 2014;**128**:356-365. DOI: 10.1016/j.fuel.2014.01.074

[38] Mesquita FMR, Bessa AMM, de Lima DD, de Sant'Ana HB, de Santiago-Aguiar RS. Liquid–liquid equilibria of systems containing cottonseed biodiesel +glycerol+ethanol at 293.15, 313.15 and 333.15K. Fluid Phase Equilibria. 2012;**318**: 51-55. DOI: 10.1016/j.fluid.2012.01.016

[39] Rodrigues CEC, Filipini A, Meirelles AJA. Phase equilibrium for systems composed by high unsaturated vegetable oils + linoleic acid + ethanol + water at 298.2 K. Journal of Chemical & Engineering Data. 2006;**51**:15-21. DOI: 10.1021/je0495841

[40] Gonçalves CB, Meirelles AJA. Liquid–liquid equilibrium data for the system palm oil + fatty acids + ethanol + water at 318.2K. Fluid Phase Equilibria. 2004;**221**:139-150. DOI: 10.1016/ j.fluid.2004.05.002

[41] Rodrigues CEC, Antoniassi R, Meirelles AJA. Equilibrium data for the system Rice bran oil + fatty acids + ethanol + water at 298.2 K. Journal of Chemical & Engineering Data. 2003;**48**: 367-373. DOI: 10.1021/je0255978

[42] Rodrigues CEC, Peixoto ECD, Meirelles AJA. Phase equilibrium for systems composed by refined soybean oil +commercial linoleic acid+ethanol +water, at 323.2K. Fluid Phase Equilibria. 2007;**261**:122-128. DOI: 10.1016/ j.fluid.2007.07.021

[43] Machado GD, Castier M, Voll AP, Cabral VF, Cardozo-Filho L, Aranda DAG. Ethanol and methanol

Unifac subgroup parameter estimation in the prediction of the liquid-liquid equilibrium of biodiesel systems. Fluid Phase Equilibria. 2019;**488**:79-86. DOI: 10.1016/j.fluid.2019.01.012

[44] Castier M, Amer MM. XSEOS: An evolving tool for teaching chemical engineering thermodynamics. Education for Chemical Engineers. 2011;**6**:e62-e70. DOI: 10.1016/j.ece.2010.12.002

[45] Lasdon LS, Waren AD, Jain A, Ratner M. Design and testing of a generalized reduced gradient code for nonlinear programming. ACM Transactions on Mathematical Software. 1978;**4**:34-50. DOI: 10.1145/355769.355773

[46] Magnussen T, Rasmussen P, Fredenslund A. UNIFAC parameter table for prediction of liquid-liquid equilibria. Industrial & Engineering Chemistry Process Design and Development. 1981; **20**:331-339. DOI: 10.1021/i200013a024

[47] Homrich POB, Ceriani R. Phase equilibria for systems containing refined soybean oil plus Cosolvents at different temperatures. Journal of Chemical & Engineering Data. 2018;**63**:1937-1945. DOI: 10.1021/acs.jced.7b01051

[48] Homrich POB, Dias LG, Mariutti LRB, Bragagnolo N, Ceriani R. Liquid-liquid equilibria and density data for pseudoternary systems of refined soybean oil + (hexanal, or heptanal, or butyric acid, or valeric acid, or caproic acid, or caprylic acid) + dimethyl sulfoxide at 298.15 K. The Journal of Chemical Thermodynamics. 2019;**131**: 149-158. DOI: 10.1016/j.jct.2018.10.030

[49] Kanda LR, Voll FAP, Corazza ML. LLE for the systems ethyl palmitate (palmitic acid)(1)+ethanol(2)+glycerol (water)(3). Fluid Phase Equilibria. 2013; **354**:147-155. DOI: 10.1016/j.fluid.2013. 06.027

*Edited by Rampal Pandey, Israel Pala-Rosas, José L. Contreras and José Salmones*

Due to their unique physicochemical properties, low cost, and wide availability, ethanol and glycerol have gained attention for their use as alternative feedstocks in the sustainable production of several commodity and specialty products. As a result, during the last decades, there has been intense research aimed at developing the potential applications of these biomass-derived compounds. *Ethanol and Glycerol Chemistry - Production, Modelling, Applications, and Technological Aspects* discusses recent advances and different aspects of the production, direct applications, and processing of ethanol and glycerol from a multidisciplinary perspective that includes the medical field, fuels, and chemical synthesis.

Published in London, UK © 2023 IntechOpen © Perception7 / iStock

Ethanol and Glycerol Chemistry - Production, Modelling, Applications,

and Technological Aspects

Ethanol and Glycerol

Chemistry

Production, Modelling, Applications,

and Technological Aspects

*Edited by Rampal Pandey, Israel Pala-Rosas,* 

*José L. Contreras and José Salmones*