Diffusion of hydrocarbons diluted in supercritical carbon dioxide

Abstract

Abstract Mutual diffusion of six hydrocarbons (methane, ethane, isobutane, benzene, toluene or naphthalene) diluted in supercritical carbon dioxide ( $${\hbox {CO}}_{2}$$ CO 2 ) is studied by molecular dynamics simulation near the Widom line, i.e., in the temperature range from 290 to 345 K along the isobar 9 MPa. The $${\hbox {CO}}_{2}$$ CO 2 + aromatics mixtures are additionally sampled at 10 and 12 MPa and an experimental database with Fick diffusion coefficient data for those systems is provided. Taylor dispersion experiments of $${\hbox {CO}}_{2}$$ CO 2 with benzene, toluene, n-dodecane and 1,2,3,4-tetrahydronaphthalene are conducted along the $$p =$$ p = 10 MPa isobar. Maxwell–Stefan and Fick diffusion coefficients are analyzed, together with the thermodynamic factor that relates them. It is found that the peculiar behavior of the Fick diffusion coefficient of some $${\hbox {CO}}_{2}$$ CO 2 mixtures in the extended critical region is a consequence of the thermodynamic factor minimum due to pronounced clustering on the molecular scale. Further, the strong dependence of the Fick diffusion coefficient on the molecular mass of the solute as well as the breakdown of the Stokes–Einstein relation near the Widom line are confirmed. Eleven correlations for the prediction of the Fick diffusion coefficient of $${\hbox {CO}}_{2}$$ CO 2 mixtures are assessed. An alternative two-step approach for the prediction of the infinite dilution Fick diffusion coefficient of supercritical $${\hbox {CO}}_{2}$$ CO 2 mixtures is proposed. It requires only the state point in terms of temperature and pressure (or density) as well as the molecular solute mass as input parameters. First, entropy scaling is applied to estimate the self-diffusion coefficient of $${\hbox {CO}}_{2}$$ CO 2 . Subsequently, this coefficient is used to determine the infinite dilution Fick diffusion coefficient of the mixture, based on the finding that these two diffusion coefficients exhibit a linear relationship, where the slope depends only on the molecular solute mass.