Non-equilibrium Molecular Dynamics Simulation on Pure Gas Permeability Through Carbon Membranes

Use of Boundary-Driven Nonequilibrium Molecular Dynamics for Determining Transport Diffusivities of Multicomponent Mixtures in Nanoporous Materials

The boundary-driven molecular modeling strategy to evaluate mass transport coefficients of fluids in nanoconfined media is revisited and expanded to multicomponent mixtures. The method requires setting up a simulation with bulk fluid reservoirs upstream and downstream of a porous media. A fluid flow is induced by applying an external force at the periodic boundary between the upstream and downstream reservoirs. The relationship between the resulting flow and the density gradient of the adsorbed fluid at the entrance/exit of the porous media provides for a direct path for the calculation of the transport diffusivities. It is shown how the transport diffusivities found this way relate to the collective, Onsager, and self-diffusion coefficients, typically used in other contexts to describe fluid transport in porous media. Examples are provided by calculating the diffusion coefficients of a Lennard-Jones (LJ) fluid and mixtures of differently sized LJ particles in slit pores, a realistic model of methane in carbon-based slit pores, and binary mixtures of methane with hypothetical counterparts having different attractions to the solid. The method is seen to be robust and particularly suited for the study of study of transport of dense fluids and liquids in nanoconfined media.

Structure-Property Relationship on the Example of Gas Separation Characteristics of Poly(Arylene Ether Ketone)s and Poly(Diphenylene Phtalide)

Three poly(arylene ether ketone)s (PAEKs) with propylidene (C1, C2) and phtalide (C3) fragments, and one phtalide-containing polyarylene (C4), were synthesized. Their chemical structures were confirmed via 1H NMR, 13C NMR and 19F NMR spectroscopy. The polymers have shown a high glass transition temperature (>155 °C), excellent film-forming properties, and a high free volume for this polymer type. The influence of various functional groups in the structure of PAEKs was evaluated. Expectedly, due to higher free volume the introduction of hexafluoropropylidene group to PAEK resulted in higher increase of gas permeability in comparison with propylidene group. The substitution of the fluorine-containing group on a rigid phtalide moiety (C3) significantly increases glass transition temperature of the polymer while gas permeation slightly decreases. Finally, the removal of two ether groups from PAEK structure (C4) leads to a rigid polymer chain that is characterized by highest free volume, gas permeability and diffusion coefficients among the PAEKs under investigation. Methods of modified atomic (MAC) and bond (BC) contributions were applied to estimate gas permeation and diffusion. Both techniques showed reasonable predicted parameters for three polymers while a significant underestimation of gas transport parameters was observed for C4. Gas solubility coefficients for PAEKs were forecasted by "Short polymer chain surface based pre-diction" (SPCSBP) method. Results for all three prediction methods were compared with the ex-perimental data obtained in this work. Predicted parameters were in good agreement with ex-perimental data for phtalide-containing polymers (C3 and C4) while for propylidene-containing poly(arylene ether ketone)s they were overestimated due to a possible influence of propylidene fragment on indices of oligomeric chains. MAC and BC methods demonstrated better prediction power than SPCSBP method.

Flow of methane in shale nanopores at low and high pressure by molecular dynamics simulations

Flow in shale nanopores may be vastly different from that in the conventional permeable media. In large pores and fractures, flow is governed by viscosity and pressure-driven. Convection describes the process. Pores in some shale media are in nanometer range. At this scale, continuum flow mechanism may not apply. Knudsen diffusion and hydrodynamic expressions such as the Hagen-Poiseuille equation and their modifications have been used to compute flow in nanopores. Both approaches may have drawbacks and can significantly underestimate molecular flux in nanopores. In this work, we use the dual control volume-grand canonical molecular dynamics simulations to investigate methane flow in carbon nanopores at low and high pressure conditions. Our simulations reveal that methane flow in a slit pore width of 1-4 nm can be more than one order of magnitude greater than that from Knudsen diffusion at low pressure and the Hagen-Poiseuille equation at high pressure. Knudsen diffusion and Hagen-Poiseuille equations do not account for surface adsorption and mobility of the adsorbed molecules, and inhomogeneous fluid density distributions. Mobility of molecules in the adsorbed layers significantly increases molecular flux. Molecular velocity profiles in nanopores deviate significantly from the Navier-Stokes hydrodynamic predictions. Our molecular simulation results are in agreement with the enhanced flow measurements in carbon nanotubes.