Computation of partial molar properties using continuous fractional component Monte Carlo

Journal Article (2018)
Author(s)

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)

Research Group
Engineering Thermodynamics
Copyright
© 2018 A. Rahbari, R. Hens, I.K. Nikolaidis, A. Poursaeidesfahani, M. Ramdin, I. G. Economou, O. Moultos, D. Dubbeldam, T.J.H. Vlugt
DOI related publication
https://doi.org/10.1080/00268976.2018.1451663
More Info
expand_more
Publication Year
2018
Language
English
Copyright
© 2018 A. Rahbari, R. Hens, I.K. Nikolaidis, A. Poursaeidesfahani, M. Ramdin, I. G. Economou, O. Moultos, D. Dubbeldam, T.J.H. Vlugt
Research Group
Engineering Thermodynamics
Issue number
21-22
Volume number
116
Pages (from-to)
3331-3344
Reuse Rights

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.