## Abstract

Monte Carlo (MC) simulations in ensembles with a fixed chemical potential or fugacity, for example the grand-canonical or the osmotic ensemble, are often used to compute phase equilibria. Chemical potentials can be computed either with an equation of state (EoS) or from molecular simulations. The accuracy of the computed chemical potentials depends on the quality of the (critical) parameters used in the EoS and the applied force field in the simulations. We investigated the consistency of both approaches for computing fugacities of the industrially relevant gases CO_{2}, CH_{4}, CO, H_{2}, N_{2}, and H_{2}S. The critical temperature (T_{c}), pressure (P_{c}), and acentric factors (ω) of these gases are computed from MC simulations in the Gibbs ensemble. The effect of cutoff radius and tail corrections on the computed values of T_{c}, P_{c}, and ω is investigated. In addition, MC simulations in the Gibbs ensemble are used to compute the VLE of the 15 possible binary systems comprising the gases CO_{2}, CH_{4}, CO, H_{2}, N_{2}, and H_{2}S, and the ternary systems CO_{2}/CH_{4}/H_{2}S and CO_{2}/CO/H_{2}. Binary interaction parameters (k_{ij}) of these natural/synthesis gas mixtures are obtained by fitting the Peng-Robinson (PR) EoS to the binary VLE data from the MC simulations. The computed properties from the MC simulations are compared with the PR EoS, the GERG EoS, and experimental results. The MC results show that including tail corrections in the simulations is crucial to obtain accurate critical properties. The force fields used for the gases can reproduce the fugacities of the gases within 5% of the experimental data. The dew-point curves of all the 15 binaries were predicted correctly by the MC simulations, but the bubble-point curves for the systems H_{2}/CO, CH_{4}/H_{2}, H_{2}S/N_{2}, and H_{2}S/CO significantly deviate from the experiments.

Original language | English |
---|---|

Pages (from-to) | 174-181 |

Journal | Fluid Phase Equilibria |

Volume | 428 |

DOIs | |

Publication status | Published - 2016 |

## Keywords

- Binary interaction parameter
- Equation of state
- Molecular simulation
- Natural gas
- Synthesis gas