(259a) Modelling of Multiphasic Behavior of Biodiesel Transesterification Operating Below Critical Conditions Using CO2 as a Co-solvent with PC-SAFT EoS

Rodriguez, G., University of Pittsburgh
Beckman, E. J., University of Pittsburgh
Biodiesel is composed of fatty acid esters generated via methanolysis or ethanolysis of triglycerides. The preferred method for biodiesel production is transesterification of triglycerides with methanol or ethanol. Phase behavior plays a critical role in biodiesel transesterification. The reaction is generally carried out at high temperatures, pressures and high short alcohols to oil ratios, in order to increase the solubility of non-polar and polar phases.

The process can be performed uncatalyzed using supercritical conditions: temperatures of 250 - 400 °C and pressures of 19- 24 MPa with high methanol/ethanol to oil ratios of 40:1 molar. The severity of these conditions guarantees good reaction kinetics and the formation of a singular phase in the reactor. This process has recently received much attention in literature however it still presents challenges regarding the intensive energy requirements for the process and downstream separation of products from excess reactants.1–6

The catalyzed process has different approaches including acid and basic homogeneous catalysts and is the current commercial path to produce biodiesel.7 however ongoing challenges regarding the presence of water, free fatty acid impurities, and extra cost associated with the separation of the biodiesel from glycerol are still present.7,8

Heterogeneous catalysts have also been investigated showing economical profitability, and provide a solution for the final catalyst separation.9 They vary in nature from alkaline earth, titanium silicates, anion exchange resins and polymers,10–14

Carbon dioxide has been suggested as co-solvent with beneficial properties for biodiesel transesterification at supercritical conditions because: it is abundant in nature, it has been widely used in vegetable oil extraction, it has a rather mild critical value (31.1 °C and 7.39 MPa), and increases the solubility between alcohols and triglycerides favoring conversion15,16.

Findings from Soh et al17 showed triglycerides transesterification does not need to be done at supercritical conditions and can be coupled with heterogeneously catalyzed reactions. Carbon dioxide works as a co-solvent for both the polar and oleic phase. Expanding the two phases with CO2 increases the solubilized triglycerides on the alcohol phase, which consequently encourages the transesterification reaction. Almost complete conversion was obtained in moderate reaction times (< 2 Hr). Also, temperatures around 95 °C and pressures of 9.5 MPa were enough to achieve 98% methyl oleate yield in a tri-phasic catalyzed reaction process. Benefitial transport properties

Despite the system's simplicity, some species present on the biodiesel reaction, namely methanol/ethanol and glycerol have a high degree of molecular interaction as hydrogen bonding, and relevant polar forces. Understanding the systems phase behavior development from a modelling perspective requires models that can capture the complexity of this interactions.

Perturbed Chain Statistically Association Fluid Theory (PC-SAFT) by Gross and Sadowski,18,19 has become a widely used engineering resource due to its accurate prediction capabilities and has been successfully applied to many associative mixtures. Additionally PC-SAFT has been extended to polar compounds.20–24 In particular Nguyen et al.25,26 and Corazza et al.27 looked into biodiesel related phase equilibria, parametrizing pure methyl- and ethyl-esters at different conditions and obtaining promising results for methanol, ethanol, glycerol, and biodiesel binary and ternary systems.

Up to seven molecular parameters describe each pure component phase and are fitted to saturation data. The equation of states is written as an expansion of the Helmholtz free energy segregated into a reference fluid and the main interaction forces- i.e. dispersion, hydrogen bonding, and polar forces. Binary interaction coefficients are then calculated using binary data to mitigate accumulated errors that result from averaging approximations on mixtures.

Modelling of the biodiesel reaction in the presence of carbon dioxide has focused on supercritical conditions and most of the work done has been performed with cubic equations of state. Modelling this reaction at milder conditions of pressure and temperature where multiple phases are present using SAFT like equations has not been done before.

Modelling results achieved accurate phase behavior representation for pure, binary and ternary systems including carbon dioxide using a polar version of PC-SAFT. Group contribution methods allowed generalization for methyl-esters and ethyl-esters modelling simplifying the required components basis and minimizing the number of parameters. Small errors were obtained using very low values of binary interaction coefficients (below 7%) for binary systems. Successful modelling for liquid-liquid and vapor-liquid-liquid regimes suggest the possibility to determine optimal conditions for transesterification performed in multiphasic reactors in presence of CO2

The global aim of this study is to fully characterize biodiesel-related compounds with carbon dioxide as a co-solvent with the PC-SAFT, predict supercritical, and sub-critical phase behavior generating binary and ternary diagrams to describe interactions amongst species. These calculations could potentially suggest a new process design by estimating favorable thermodynamic conditions as pressure, temperature and reactants ratios for biodiesel production using carbon dioxide as a co-solvent in subcritical conditions based on their phase behavior equilibria.

Keywords: Biodiesel, PC-SAFT, Phase-Behavior, Carbon-Dioxide, Co-solvent


  1. Madras, G., Kolluru, C. & Kumar, R. Synthesis of biodiesel in supercritical fluids. Fuel 83, 2029–2033 (2004).
  2. Balat, M. Biodiesel Fuel Production from Vegetable Oils via Supercritical Ethanol Transesterification. Energy Sources, Part A Recover. Util. Environ. Eff. 30, 429–440 (2008).
  3. He, H., Wang, T. & Zhu, S. Continuous production of biodiesel fuel from vegetable oil using supercritical methanol process. Fuel 86, 442–447 (2007).
  4. He, H., Sun, S., Wang, T. & Zhu, S. Transesterification kinetics of soybean oil for production of biodiesel in supercritical methanol. JAOCS, J. Am. Oil Chem. Soc. 84, 399–404 (2007).
  5. Demirbas, A. Biodiesel from vegetable oils via transesterification in supercritical methanol. Energy Convers. Manag. 43, 2349–2356 (2002).
  6. Sawangkeaw, R., Bunyakiat, K. & Ngamprasertsith, S. A review of laboratory-scale research on lipid conversion to biodiesel with supercritical methanol (2001-2009). J. Supercrit. Fluids 55, 1–13 (2010).
  7. Pinnarat, T. & Savage, P. E. Assessment of noncatalytic biodiesel synthesis using supercritical reaction conditions. Ind. Eng. Chem. Res. 47, 6801–6808 (2008).
  8. Vicente, G., Martínez, M. & Aracil, J. Integrated biodiesel production: A comparison of different homogeneous catalysts systems. Bioresour. Technol. 92, 297–305 (2004).
  9. West, A. H., Posarac, D. & Ellis, N. Assessment of four biodiesel production processes using HYSYS.Plant. Bioresour. Technol. 99, 6587–6601 (2008).
  10. Schuchardt, U., Vargas, R. M. & Gelbard, G. Transesterification of soybean oil catalyzed by alkylguanidines heterogenized on different substituted polystyrenes. J. Mol. Catal. A Chem. 109, 37–44 (1996).
  11. Georgogianni, K. G., Katsoulidis, A. K., Pomonis, P. J., Manos, G. & Kontominas, M. G. Transesterification of rapeseed oil for the production of biodiesel using homogeneous and heterogeneous catalysis. Fuel Process. Technol. 90, 1016–1022 (2009).
  12. Gryglewicz, S. Rapeseed oil methyl esters preparation using heterogeneous catalysts. 70, (1999).
  13. Demirbas, a. Biodiesel from Vegetable Oils with MgO Catalytic Transesterification in Supercritical Methanol. Energy Sources, Part A Recover. Util. Environ. Eff. 30, 1645–1651 (2008).
  14. Shibasaki-Kitakawa, N. et al. Biodiesel production using anionic ion-exchange resin as heterogeneous catalyst. Bioresour. Technol. 98, 416–421 (2007).
  15. Maçaira, J., Santana, A., Recasens, F. & Angeles Larrayoz, M. Biodiesel production using supercritical methanol/carbon dioxide mixtures in a continuous reactor. Fuel 90, 2280–2288 (2011).
  16. García-Jarana, M. B. et al. Use of supercritical methanol/carbon dioxide mixtures for biodiesel production. Korean J. Chem. Eng. 33, 2342–2349 (2016).
  17. Soh, L., Curry, J., Beckman, E. J. & Zimmerman, J. B. Effect of System Conditions for Biodiesel Production via Transesterification Using Carbon Dioxide-Methanol Mixtures in the Presence of a Heterogeneous Catalyst. ACS Sustain. Chem. Eng. 387–395 (2014). doi:10.1021/sc400349g
  18. Gross, J. & Sadowski, G. Perturbed-Chain SAFT: An Equation of State Based on a Pertubation Theory of Chain Molecules. Ind. Eng. Chem. Res 1244–1260 (2001).
  19. Tumakaka, F., Gross, J. & Sadowski, G. Modeling of polymer phase equilibria using Perturbed-Chain SAFT. Fluid Phase Equilib. 194–197, 541–551 (2002).
  20. Gubbins, K. E. & Twu, C. H. Thermodynamics of polyatomic fluid mixtures—I theory. Chem. Eng. Sci. 33, 863–878 (1978).
  21. Twu, C. H. & Gubbins, K. E. Thermodynamics of Polyatomic Fluids Mixtures-II. Chem. Eng. Sci. 33, 879–887 (1977).
  22. Twu, C. H., Gubbins, K. E. & Gray, C. G. Thermodynamics of mixtures of nonspherical molecules. III. Fluid phase equilibria and critical loci. J. Chem. Phys. 64, 5186–5197 (1976).
  23. Vrabec, J. & Gross, J. Vapor - Liquid Equilibria Simulation and an Equation of State Contribution for Dipole - Quadrupole Interactions. 51–60 (2008).
  24. K. Jog, W. G. Chapman, P. Application of Wertheim’s thermodynamic perturbation theory to dipolar hard sphere chains. Mol. Phys. 97, 307–319 (1999).
  25. Nguyen-Huynh, D., Passarello, J.-P., Tobaly, P. & de Hemptinne, J.-C. Modeling Phase Equilibria of Asymmetric Mixtures Using a Group-Contribution SAFT (GC-SAFT) with a k ij Correlation Method Based on London’s Theory. 1. Application to CO 2 + n -Alkane, Methane + n -Alkane, and Ethane + n -Alkane Systems. Ind. Eng. Chem. Res. 47, 8847–8858 (2008).
  26. Nguyen Huynh, D., Nguyen Thi, T. & Van Dinh, S. T. Predicting the temperature / pressure dependent density of biodieselfuels. Petrovietnam 10, 46–58 (2012).
  27. Corazza, M. L., Fouad, W. A. & Chapman, W. G. PC-SAFT predictions of VLE and LLE of systems related to biodiesel production. Fluid Phase Equilib. 416, 130–137 (2016).


This paper has an Extended Abstract file available; you must purchase the conference proceedings to access it.


Do you already own this?



AIChE Members $150.00
AIChE Graduate Student Members Free
AIChE Undergraduate Student Members Free
Non-Members $225.00