Molecular dynamics study of the thermodynamic properties of triglyceride/methanol mixtures in the presence of cosolvent

Abstract

Biodiesel is a renewable source of energy that has the potential to replace fossil fuel-based resources. It is produced via the transesterification reaction in which a triglyceride reacts with methanol in the presence of a catalyst. Cosolvents are usually added to the reactants in transesterification to improve the mutual solubility of methanol and triglycerides and to create a single phase. The commonly used cosolvents in transesterification include tetrahydrofuran (THF) and diethyl ether (DEE). Understanding the thermodynamic properties of this system is important in the selection of cosolvent. In this study, molecular dynamics simulations were used to study the mixing behavior between the components of ternary mixtures of triolein/methanol/THF and triolein/methanol/DEE via the application of Flory-Huggins theory of solutions. A simple approach to calculation of the Flory-Huggins interaction parameters is used in which the mixing energy is derived from averaged non-bonded pair interaction energies. Apart from pure triolein, these ternary systems resemble the experimentally studied ternary systems employing various plant-based oils. The binary interaction parameters, Gibbs Free energy and chemical potentials of the components were calculated to study the mixing behavior. Our results show that along the experimental phase boundary the binary interaction parameters remain relatively unchanged with a deviation observed at the highest volume fractions of triolein where the highest binary interactions were observed indicating the possibility of phase separation. The Gibbs energy of mixing for the systems changed from positive to negative as one moves from two phase to single phase part of the ternary system as expected and in agreement with experimental data. The chemical potentials of the ternary mixture components were calculated from the partial derivatives of the Gibbs free energy of mixing. The activities calculated from Flory-Huggins chemical potentials were compared to those calculated using the UNIFAC model. There is good agreement between the two, with both able to predict a single and biphasic mixture. However, disagreements are noted at low and high concentrations of triolein.