Monte Carlo Simulations for the Estimation of the Effective Permeability of Mixed-Matrix Membranes

Recent years have seen the explosive development of mixed-matrix membranes (MMMs) for a myriad of applications. In gas separation, it is desired to concurrently enhance the permeability, selectivity and physicochemical properties of the membrane. To help achieving these objectives, experimental characterization and predictive models can be used synergistically. In this investigation, a Monte Carlo (MC) algorithm is proposed to rapidly and accurately estimate the relative permeability of ideal MMMs over a wide range of conditions. The difference in diffusivity coefficients between the polymer matrix and the filler particle is used to adjust the random progression of the migrating species inside each phase. The solubility coefficients of both phases at the polymer-filler interface are used to control the migration of molecules from one phase to the other in a way to achieve progressively phase equilibrium at the interface. Results for various MMMs were compared with the results obtained with the finite difference method under identical conditions, where the results from the finite difference method are used in this investigation as the benchmark method to test the accuracy of the Monte Carlo algorithm. Results were found to be very accurate (in general, <1% error) over a wide range of polymer and filler characteristics. The MC algorithm is simple and swift to implement and provides an accurate estimation of the relative permeability of ideal MMMs. The MC method can easily be extended to investigate more readily non-ideal MMMs with particle agglomeration, interfacial void, polymer-chain rigidification and/or pore blockage, and MMMs with any filler geometry.