Abstract

Ionic liquid mixed matrix membranes (IL3Ms) were synthesized using polyethersulfone (PES) as the base polymer and silica-aluminophosphate (SAPO-34) as the dispersed particles, and their CO2 permeation was investigated. Three of the most widely used models for gas separation—the Maxwell, Lewis–Nielson, and Maxwell–Wagner–Sillar (MWS) models—were then applied to the membranes. Large deviations were found between the model predictions and experimental data. FESEM images suggested that local agglomeration and disorientation of the SAPO-34 particles within the membrane afforded substantial changes in the morphology. The MWS model, which considers the shape factor, was modified to incorporate the volume fraction of the wetted dispersed phase and the ideal shape factor. A direct relationship was found between the filler concentration and the shape factor. The modified model was shown to produce absolute and relative errors of less than 3%. When validated against data from the literature, the deviation remained within 5%. The modified model can be used to estimate the gas permeance of an IL3M.

1. Introduction

Membranes, both polymeric and inorganic, have several specific advantages over conventional separation processes and are finding increasing applications in both gas and liquid separation. This has accelerated research in this area, particularly in the identification and development of new membrane materials. In gas separation, the major challenge is to synthesize membranes that can touch, if not cross, the Robeson plot [1] in terms of flux and selectivity or separation factor. Mixed matrix membranes (MMMs) are promising novel materials, and their high potential in both gas and liquid separation has already been recognized.

Transport of small gas molecules through a polymeric dense membrane occurs when random molecular motions are projected by the polymer matrix. The driving force in this process of sorption, diffusion, and permeation is the concentration gradient or partial pressure difference across the membrane [2, 3]. In sieving materials, gas transportation is driven by the difference in diffusion rates between permeating molecules of different shapes and sizes. More specifically, the rate of diffusion is determined by the gas penetrant and the sieving pore size [4]. Molecular sieving materials are known to offer high separation factors since their pore dimensions are similar to the dimensions of the gas molecules. Diffusion occurs when the molecules jump through a narrow opening between one sorption cavity and the next [5]. Adsorption and surface diffusion are further factors in gas permeation through small-pore-diameter membranes [6].

MMMs are heterogeneous membranes that generally constitute inorganic fillers dispersed within a polymer matrix, although the inclusion of organic or metal-organic framework materials has also shown excellent potential [7]. In the ideal MMM morphology, a defect-free boundary is assumed between the dispersed particles and polymer interface [8]. The strong adhesion between the polymer and filler reduces the free volume in the polymer [9]—a process known as polymer rigidification [10]. In most cases, this rigidification [11] improves selectivity while reducing permeability [12].

With proper selection of the filler, many MMMs are excellent separation media for use in gas separation. However, there remains a need to improve both the flux and selectivity of MMMs. This can be achieved by appropriate manipulation of the membrane. In the proposed approach, an ionic liquid is incorporated into the membrane matrix. Ionic liquids are a large family of green solvents that demonstrate excellent solvent properties in addition to being nonvolatile and stable even at high temperatures. For membrane gas separation, many studies have investigated the use of ionic liquids in the form of supported liquid membranes (SLMs) and have yielded encouraging results [1315]. Ionic liquids have also been polymerized to produce ionic liquid polymeric membranes (ILPMs) [1618]. However, the permeation separation performance of these membranes is still lower than that of SLM and is yet to exceed the 2008 Robeson upper bound. A study mentioned by Ilyas et al. [19] has successfully synthesized supported ionic liquid membranes with 41 CO2/CH4 selectivity with no observation of ionic liquid leaching. It was a very good finding but it was limited to 10 bar, and the group has not focusing the effects of inorganic materials. Another work from Bhattacharya and Mandal [20], which was also not considering the effects of inorganic materials, has successfully tested their MMM at a maximum pressure of ∼7 bar and found the membranes work better for H2S/CH4 separation. Based on these studies of using ionic liquid, it can be concluded that the ionic liquid can be the way forward of commercializing the membranes for gas separation purposes.

In this study, an ionic liquid was homogeneously distributed in the casting suspension and used to prepare an MMM comprising both fine adsorbent particles (SAPO-34) and an ionic liquid. The goal of this study was to achieve improvements in both flux and selectivity. The study also attempted to develop models suitable for predicting gas transportation in the membrane that considers the membrane’s composition, characteristics, and morphology. A proper theoretical understanding of gas transportation in an MMM is required to support the design of research and development studies and to allow its implementation in industrial separation [21].

Gas transportation across MMMs has been a major focus of recent membrane research, mainly owing to their increasing use in gas-sweetening applications. Most existing theoretical models rely on old equations that were developed to address the quite different problem of estimating the electrical conductivity of a nonhomogeneous substance comprising a dispersed solid phase. By assuming that an analogy can be drawn with gas permeation, these conductivity models have been used to evaluate gas separation [22]. Similar electrical and thermal conductivity models have been applied to the analysis of MMMs [23].

Different models have used different approaches to predict the relative permeability of gases. Based on a review from Aroon et al. [8, 24], there are few permeation models used to predict effective permeability of a gaseous penetrant for mixed matrix membranes comprising porous particles. Another study from Hashemifard et al. [23] was also considering the models reviewed in [24] for their prediction of the relative permeability. In this study, the Maxwell, Lewis–Nielson, and Maxwell–Wagner–Sillar (MWS) models were used. However, these models cannot predict the permeability of a gas species through an MMM comprising a third component. One objective of the study was therefore to evaluate the existing models and to establish the adjustments needed to make them applicable for predicting gas transport across an IL3M. The MMM used in this study was prepared using polyethersulfone (PES), SAPO-34, and an ionic liquid. The analysis focused on the effect of the ionic liquid on the dispersion geometry of the SAPO-34 particles. This reflected our reported finding that incorporation of an ionic liquid triggers agglomeration of an inorganic filler within the membrane matrix [25].

2. Evaluation of Existing MMM Transport Models

Predicting the permeability of an MMM has proved a challenge. Interestingly, the theoretical models developed to address electrical conductivity in a nonhomogeneous medium are applicable to this entirely different problem, including the models proposed by Maxwell, Bruggeman, Pal, and Lewis–Nielson [8, 21, 23, 26, 27]. These models have received attention primarily because of their simplicity—the permeability of the MMM can be assumed to depend on those of the base polymer and the dispersed phase and the loading of the filler in terms of volume fraction [10]. The MWS model also incorporates the particle shape factor. Hashemifard et al. [15] and Nasir et al. [28] reviewed the applicability of the models at different particle loadings. Table 1 lists the available models and the parameters they require to estimate permeability.

The Maxwell model (equation (1) in Table 1) was first developed in 1873 to predict the permittivity of dielectric conduction through a heterogeneous medium [34]. The model involves a constitutive equation in which the heterogeneous medium is assumed to be an analogue of the membrane flux. This allows the model to be applied for gas transportation in an MMM. However, the Maxwell model fails when a high loading of filler in the MMM produces the maximum packing of the volume fraction. Furthermore, this model is not suitable for predicting permeation when parameters such as the shape, size distribution, and aggregation of the particles must be considered.

An alternative model was originally proposed to derive the elastic modulus of particulate composites. This is known as the Lewis–Nielson model and is given by equation (4) in Table 1 [32]. For random close packing of uniform spheres, is normally set at 0.64 [22]. When , the permeability ratio λdm ⟶ ∞ as the relative permeability diverges. In this case, is very sensitive as it is affected by the shape, size distribution, and aggregation of the particles. However, the Lewis–Nielson model does not consider these parameters. As , the Lewis–Nielson model reduces to the Maxwell model.

The Pal model (equation (3) in Table 1) was developed for calculating the thermal conductivity of a composite particulate and has also been used to calculate permeability [35]. This model addresses the filler particle packing problem by applying a differential effective medium approach. The Pal model correctly predicts the behavior as by considering the effect of morphology on permeability through parameter . However, it involves an implicit relationship that must be solved numerically for the relative permeability, Pr.

Bouma et al. used the MWS model to predict the gas transportation characteristics of an MMM [33]. This model is given by equation (6) in Table 1, where n is the shape factor of the particle. Notably, the maximum and minimum values of the shape factor are used to derive the effective permeability of gas diffusion in an MMM, and these can be given by series or parallel models, respectively. However, no current model allows the effect of a third component (e.g., an ionic liquid) homogeneously distributed in the continuous membrane matrix to be considered when calculating the gas permeability of an IL3M. This study addressed the effect of adding an ionic liquid to an MMM and the optimization of the parameter affected by the ionic liquid in the theoretical calculation of CO2 permeability.

3. Membrane Synthesis and Permeation Measurement

An MMM is a combination of polymer and molecular sieves. By combining the process characteristics of the polymer matrix and the superior transportation properties of the molecular sieve, MMMs combine the benefits of both materials [36]. In recent research, MMMs have been combined with ionic liquids to produce IL3Ms.

3.1. Synthesis of the Ionic Liquid Mixed Matrix Membranes

In this study, PES was chosen to provide the polymer support required for fabricating the enhanced MMM. This was the commercial Ultrason E 6020 P used in industrial processing and was purchased from BASF, Germany, in flake form. Commercial SAPO-34 with a kinetic diameter of 0.38 nm was purchased from ACS Material in powder form. Two ionic liquids of 99% purity were used: 1-ethyl-3-methyl-imidazolium bis(trifluoromethyl-sulfonyl)imide (emim[Tf2N]) and 1-ethyl-3-methylimidazolium trifluoromethanesulfonate (emim[CF3SO3]). These were purchased from Sigma Aldrich and used without further purification. The IL3M synthesized using emim[Tf2N] (IL3M-emim[Tf2N]) was used for model correlation, while that synthesized using emim[CF3SO3] (IL3M-emim[CF3SO3]) was used for validation.

The PES flakes were dried overnight prior to use to remove any trapped moisture. The ionic liquid was mixed with the N-methyl-2-pyrrolidone (NMP) solvent, and the SAPO-34 particles were dispersed in the mixture. The PES was then dried and weighed, before being dissolved in the mixture. The solution was stirred continuously to ensure homogeneity and degassed to remove bubbles prior to casting. The cast membranes were dried in an oven at 160°C for 24 h. The details of the materials used, method of membrane preparation, and membrane characterization are presented in more detail in our recent publication [25]. For use as a benchmark, a PES-based polymeric membrane was synthesized. The IL3Ms were synthesized using a fixed ionic-liquid concentration of 20 wt.% and SAPO-34 filler concentrations of 10, 20, and 30 wt.%.

3.2. Membrane Morphology Characterization

Membrane morphology is critical for understanding transport mechanisms. This was characterized in the present study to clarify the effect of the added ionic liquid on the SAPO-34 orientation in the membrane matrix. Field emission scanning electron microscopy (FESEM, ZEIS SUPRA TM 55VP) was used for this characterization. Before scanning, the membranes were cut using liquid nitrogen to provide a smooth surface. Figure 1 shows the FESEM membrane sample of PES-SAPO-34 mixed matrix membrane with 10 wt.% ionic liquids. The average membrane thickness of all synthesized membranes is ∼10 µm.

3.3. Gas Separation Performance Analysis

Gas permeance analyses were conducted following the procedure reported in our previous study [25]. Gas permeation was then calculated from the time taken by the gas mixture to pass across the membrane at constant volume. This was recorded using a bubble flow meter, after a steady-state condition had been reached. The separation performance study was conducted at a pressure range of 10–30 bar. To ensure the reliability of the results, three membranes were fabricated using each composition, and all the fabricated membranes were tested using the same gas permeation method.

4. Present Model

Three of the models listed in Table 1 were selected for validation purposes: the Maxwell, Lewis–Nielson, and MWS models. Statistical analysis was performed to assess the level of agreement between the experimental data and the predictions made by the models. When deriving the gas permeance using the theoretical models, a defect-free PES membrane was assumed. The experimental data on CO2 permeance are given in Table 2. The CO2 permeance of SAPO-34 was taken from the literature [37]. Since the only data available on CO2 permeance through an inorganic SAPO-34 membrane were obtained at a test pressure of 20 bar, the experimental data were also recorded at 20 bar. Two types of imidazolium-based ionic liquids were used to synthesize the IL3Ms. The gas permeation data for the IL3M-emim[Tf2N] were used for model development.

The MWS model required the shape factor of the SAPO-34 particles to be specified. Based on [38], a value of n = 0.33 was used because SAPO-34 has a cubic form, as shown in the FESEM image in Figure 2. Figure 3 compares the relative permeance of CO2 against SAPO-34 loading predicted by the above-mentioned three models with the experimentally determined values.

The absolute average error (AARE) was calculated using equation (1), and Table 3 summarizes the deviations between the models and the experimental data:

From Table 3, the relative permeability predicted by the existing models showed a major deviation from the experimental data; the AARE values were high (>50%), suggesting the need for model improvement and modification. Figure 4 shows FESEM images of the SAPO-34 orientation in the PES polymer matrix at different loadings. The FESEM analysis allowed the effective shape of the fillers to be determined as the loading was increased. The addition of ionic liquid caused the SAPO-34 to agglomerate, thereby changing the overall orientation of the inorganic filler. The analyses conducted using the theoretical models highlighted the importance of considering the morphology of the particles. As noted above, MWS is the only model that allows this. This model was therefore selected for modification, to bring it into line with the experimental data. A corrected parameter is needed to enable use of the model for accurate prediction of gas separation in an IL3M.

The FESEM cross sections revealed that the orientation of the SAPO-34 depended on the filler and the ionic liquid loading. The addition of ionic liquid changed the overall orientation of the SAPO-34 such that the resulting shape was no longer cubic, as shown in Figure 2. This resulted from the SAPO-34 particles agglomerating and combining with each other. This agglomeration became more pronounced as the SAPO-34 loading was increased. The FESEM results established that the shape factor is an important parameter in determining the Pr theoretical value since it may also favor enhancement of permeation via the hopping mechanism. As shown in Figure 4, the orientation of the filler within the membranes changed as the SAPO-34 loading in the membrane was increased. The SAPO-34 agglomerated locally. This became more pronounced as the loadings reached 20 and 30 w/w%, suggesting that the overall shape of the combined filler was also changing. To consider the shape factor, the MWS model required the optimization and improvement discussed in the following section.

5. Model Improvement and Optimization

To improve the MWS model, an additional parameter was incorporated. This modeled the effect that addition of the ionic liquid had on the orientation of the particle shapes. It was assumed that the addition of ionic liquid wetted the inorganic fillers, producing an undefined geometrical orientation of the fillers in the polymer matrix. Figure 5 shows these parameters.

The ionic liquid was assumed to wet the inorganic particles, making and n sensitive to their distribution, size, and orientation. The presence of ionic liquid also caused the SAPO-34 to become wetted, thereby increasing the volume fraction of the fillers. In practice, these parameters are difficult to determine directly via FESEM analysis. The model prediction can be made accurate if these parameters are fitted for each filler loading, as the morphology is a function of the filler loading. The MWS model developed in this study for use with an IL3M was modified to incorporate a corrected shape factor, ni, and the volume fraction of the wetted dispersed phase, :where is the ideal shape factor (corrected shape factor) and is the wetted volume fraction of the dispersed phase through the IL3M. The volume fraction of the wetted dispersed phase, , was then redefined as follows:where is the ratio of the ionic liquid distance that has surrounded the particles as shown in Figure 5, determined as follows:where is the radius of the ionic liquid surrounding the dispersed phase, is the radius of the dispersed phase itself, and is determined from the correlation between and obtained from the fitting analysis. If there is no wetting effect from the ionic liquid, will reduce to .

Physically, the ideal shape factor, , is the shape factor of the SAPO-34 particles after incorporation of ionic liquid since the ionic liquid affects the SAPO-34 orientation in the polymer matrix. This is defined as follows:where is the optimized shape factor from the fitting exercise and is the maximum packing volume fraction of the inorganic filler (). If , is reduced to . Sikander (2015) obtained an expression similar to equation (5) when optimizing the shape factor [39]. Any intermediate point value should therefore be directly obtainable from equation (5).

5.1. Determination of Volume Fraction of Wetted Dispersed Phase

The objective of this approach is to optimize the volume fraction of the wetted dispersed phase and to find the correlation between and , as given by equation (3). The optimization procedures conducted against the experimental data were used to generate the values shown in Table 4. These values demonstrate that the incorporation of ionic liquid was sufficient to wet the SAPO-34, eventually increasing the volume fraction of the wetted dispersed phase, .

5.2. Determination of Ideal Shape Factor, ni

The ideal shape factor, ni, used in equation (5) was determined via optimization analysis. Figure 6 shows the optimization curve obtained by applying the optimization procedure to the original MWS model against the experimental data. The evaluation of relative CO2 permeability was repeated using the optimized shape factor no = 0.020, and a lower AARE (<25%) was returned. This analysis confirmed that the n parameter can be used to improve the predictions made by the model if the optimized values are known. Even the use of a single optimized value of no improved the MWS model predictions, and the accuracy could be further improved if the shape factor could be fitted at a specific SAPO-34 loading. This reflects the fact that the membrane morphology is a function of the distribution and aggregation of the particles.

If the shape factor at a specific SAPO-34 loading is known, the gas relative permeability at any operating pressure can be predicted since the shape factor is independent of pressure. The ni factor was introduced to reflect the orientation and distribution of the particles. Optimization of the ni value produced the results shown in Table 5.

As equation (5) shows, the ideal shape factor depends on the value of . The ni value is therefore expected to increase in line with the SAPO-34 loadings.

5.3. Determination of the Corrected Dispersed Phase Permeability, Pd

The actual dispersed phase permeability, Pd, of CO2 across SAPO-34 is critical to determine Pr(eff) in equation (6) in Table 1. Since the Pd value given in [37] was for SAPO-34 loading at 20-bar pressure, Pd (CO2) = 1731.4 GPU was used in the present study to evaluate the theoretical performance only at this pressure. In addition, the effect of pressure on relative permeability was investigated, and the Pd values for CO2 at the other testing pressures needed to be optimized accordingly.

The fitting and optimization for Pd was performed using an indirect method. The least squares method was applied to the variable parameter values to reduce the predicted errors from the model. The corrected Pd was obtained when the least error resulted. Since the fitting method worked for 10 wt.% of SAPO-34, it was also used to determine the corrected Pd at 20 wt.% and 30 wt.%. Table 6 summarizes the fitted CO2Pd values for all SAPO-34 loadings at pressures of 10, 15, 25, and 30 bar. The Pd value at 20 bar was as stated and maintained accordingly.

As can be seen from Table 6, the optimized CO2Pd values decreased with increasing pressure. This might be due to the concentration polarization of the intrinsic SAPO-34 inorganic membrane. Van de Graaf et al. reported that CH4 became enriched at the support pores, eventually depleting the driving force of the quickly permeating gas. The depletion of CO2 at increased pressures was also attributed to the high diffusion resistance caused by the enrichment of CH4 at the support [40]. Avila et al. reported that the concentration polarization became more pronounced at a pressure of 5.8 MPa, owing to the growth of separation layers at the tubular support [41].

5.4. Results of Model Modification

After applying the above modifications, the theoretical predictions of CO2 relative permeability (Pr(eff)) for the PES-SAPO-34 IL3Ms were calculated using equation (2). Table 7 shows the relationship between Pr-experimental and Pr-calculated for CO2. The AARE between the experimental and calculated Pr values was less than 2%.

For clarity, the relationship between Pr-experimental and Pr-calculated is plotted in Figure 7. The results confirm that the modified model successfully predicted the CO2 relative permeability of the IL3M.

6. Model Validation

6.1. Validation Using Different Data

As noted above, IL3M-emim[CF3SO3] was used for validation. The value of the dispersed phase permeability, Pd, was that used in the previous section. However, the volume fraction of the wetted dispersed phase, , and the ideal shape factor, ni, needed correction since a different ionic liquid was used. To determine these values, the correlations obtained from equations (3) and (5) were applied to and ni, respectively. Table 8 shows the and ni values for the IL3M-emim[CF3SO3].

The CO2 relative permeability (Pr(eff)) of the PES-SAPO-34-emim[CF3SO3] IL3M was predicted theoretically using equation (12). Table 9 shows the relationship between the Pr-experimental and Pr-calculated values. The AARE between the experimental and calculated Pr values was less than 3%. The approach was shown to achieve better agreement between the experimental and calculated values of Pr(eff). The experimental and calculated values of Pr are shown in Figure 8, with R2 = 0.95.

6.2. Validation Using Data from the Literature

In this section, we review data from the literature on membrane fabrication incorporating ionic liquids. Hudiono et al. investigated CO2/CH4 separation using MMMs comprising a styrene-based polymer, SAPO-34 as the zeolite filler, and emim[Tf2N] as the ionic liquid [17]. The CO2 permeability reported in that study is shown in Table 10. The study confirmed that the incorporation of a sufficient volume of ionic liquid into the styrene-based-SAPO-34 MMM increased CO2 permeability and selectivity. They concluded that the addition of ionic liquid had wetted the surface of the SAPO-34, making it compatible with the styrene-based poly(RTIL). The high CO2 diffusivity of the emim[Tf2N] produced an overall enhancement of gas permeability.

This study suggested that the presence of ionic liquid changed the orientation of the SAPO-34. The ideal shape factor, ni, was therefore introduced by incorporating the ionic liquid concentration into the model. The SAPO-34 fillers were also believed to have been wetted by the RTIL, and was also calculated using equation (3). The ideal shape factor, ni, was calculated using equation (5). The CO2 relative permeability was recalculated using equation (12), which incorporated both and ni. The calculated data are shown in Table 11, and the relationship between the experimental and theoretical data is presented in Figure 9.

As can be seen from Table 11, the AARE fell below 5% when the MWS model was modified by incorporating and ni.

7. Conclusions

Three of the most widely used models for estimating the permeability of a gas through an MMM were applied to membranes modified by adding an ionic liquid. All three models showed high rates of error. The effective shape factor was found to be a critical parameter in the MWS model. The addition of ionic liquid had produced local agglomeration of SAPO-34 within the membrane matrix, thereby changing the overall orientation of the filler. Therefore, the shape factor was modified to include both the ionic liquid and filler loading effects. The ideal shape factor, ni, and volume fraction of wetted dispersed phase, , were determined. Calculating relative permeance using a modified MWS model incorporating and ni showed better agreement between the experimental and predicted values. The modified model was validated against a different set of experimental data and against values reported in the literature. Both showed good agreement between the predicted and experimental values, with an AARE of less than 5%. Our results confirm that the modified MWS model can be used to predict the performance of IL3Ms.

Data Availability

Previously reported CO2 and CH4 permeability data were used to support this study and are available at https://doi.org/10.1016/j.seppur.2014.08.019. These prior studies (and datasets) are cited at relevant places within the text as reference [22].

Additional Points

(i) CO2 permeation prediction was successfully modified through the Maxwell–Wagner–Sillar model. (ii) The effects of ionic liquid loading were embedded in the modified model. (iii) The AARE% between the experimental and model values is less than 5.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to acknowledge the fund provided by the Ministry of Education (MOE), Malaysia (ERGS-0153AB-C67), and the assistance and support provided by Universiti Teknologi Petronas.