Computation of partial molar properties using continuous fractional component Monte Carlo
A. Rahbari (TU Delft - Engineering Thermodynamics)
Remco Hens (TU Delft - Engineering Thermodynamics)
I.K. Nikolaidis (National Center for Scientific Research “Demokritos”)
A. Poursaeidesfahani (TU Delft - Engineering Thermodynamics)
M. Ramdin (TU Delft - Engineering Thermodynamics)
Ioannis G. Economou (Texas A&M University at Qatar, National Center for Scientific Research “Demokritos”)
Othon Moultos (TU Delft - Engineering Thermodynamics)
D. Dubbeldam (Universiteit van Amsterdam)
T.J.H. J. H. Vlugt (TU Delft - Engineering Thermodynamics)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
An alternative method for calculating partial molar excess enthalpies and partial molar volumes of components in Monte Carlo (MC) simulations is developed. This method combines the original idea of Frenkel, Ciccotti, and co-workers with the recent continuous fractional component Monte Carlo (CFCMC) technique. The method is tested for a system of Lennard–Jones particles at different densities. As an example of a realistic system, partial molar properties of a [NH3, N2, H2] mixture at chemical equilibrium are computed at different pressures ranging from P = 10 to 80 MPa. Results obtained from MC simulations are compared to those obtained from the PC-SAFT Equation of State (EoS) and the Peng–Robinson EoS. Excellent agreement is found between the results obtained from MC simulations and PC-SAFT EoS, and significant differences were found for PR EoS modelling. We find that the reaction is much more exothermic at higher pressures.