Density and Viscosity of CO2 Binary Mixtures with SO2, H2S, and CH4 Impurities: Molecular Dynamics Simulations and Thermodynamic Model Validation

The aim of this study is to generate density and viscosity data for carbon capture utilization and storage (CCUS) mixtures using equilibrium molecular dynamics (EMD) simulations. Binary CO2 mixtures with SO2 and H2S impurities at mole fractions of 0.05, 0.10, and 0.20 were constructed. Simulations were performed across a temperature range of 223–323.15 K and at pressures up to 27.5 MPa using ms2 software. The simulation results were compared with predictions from established models. These included the Multi-Fluid Helmholtz Energy Approximation (MFHEA) for density, and the Lennard-Jones (LJ), Residual Entropy Scaling (ES-NIST), and Extended Corresponding States (SUPERTRAPP) models for viscosity. Available experimental data from the literature were also used for validation. Density predictions showed excellent agreement with MFHEA, especially for CO2 + SO2 mixtures, with %AARD values below 1% for 0.05 and 0.10, and 1.60% for 0.20 mole fraction SO2. For CO2 + H2S mixtures, deviations also increased with impurity concentration, reaching a maximum %AARD of 4.72% at 0.20 mole fraction. Viscosity data were validated against experimental values from the literature for a CO2 + CH4 (xCH4 = 0.25) mixture, showing strong agreement with both models and experiments. This confirms the reliability of the MD approach and the thermodynamic models, even for systems lacking experimental data. However, viscosity estimates showed higher uncertainty at lower temperatures and higher densities, a known limitation of the Green–Kubo method. This highlights the importance of selecting an appropriate correlation time to ensure the pressure correlation functions reach a plateau, avoiding inaccurate or uncertain viscosity values.