Mixed-Gas Selectivity Based on Pure Gas Permeation Measurements: An Approximate Model

An approximate model based on friction-coefficient formalism is developed to predict the mixed-gas permeability and selectivity of polymeric membranes. More specifically, the model is a modification of Kedem’s approach to flux coupling. The crucial assumption of the developed model is the division of the inverse local permeability of the mixture component into two terms: the inverse local permeability of the corresponding pure gas and the term proportional to the friction between penetrants. Analytical expressions for permeability and selectivity of polymeric membranes in mixed-gas conditions were obtained within the model. The input parameters for the model are ideal selectivity and solubility coefficients for pure gases. Calculations have shown that, depending on the input parameters and the value of the membrane Peclét number (the measure of coupling), there can be both a reduction and an enhancement of selectivity compared to the ideal selectivity. The deviation between real and ideal selectivity increases at higher Peclét numbers; in the limit of large Peclét numbers, the mixed-gas selectivity tends to the value of the ideal solubility selectivity. The model has been validated using literature data on mixed-gas separation of n-butane/methane and propylene/propane through polymeric membranes.


Introduction
Polymeric membranes have been widely used in various gas separation applications mainly due to their high performance regarding permeation and selectivity, as well as the easy processability of polymers. Gas permeation through a dense polymeric membrane is governed by a solution-diffusion mechanism where gas dissolution on the feed side and diffusion across the membranes determine an overall gas separation process. Although single gas permeation properties are more often reported, the selectivity of membranes for mixed gas can be different compared to the ideal selectivity based on single gas measurements [1][2][3][4][5][6][7][8][9]. In some cases, ideal selectivities are higher than mixed gas selectivities. For example, CO 2 /N 2 selectivities based on single gas measurement might be higher than mixed-gas selectivities due to the swelling of the polymer in CO 2 environment and plasticization effect [10]. The same tendency was observed for hydrocarbons separation through polyalkylmethylsiloxanes composite membranes [9]. The separation of an eightcomponent mixture of saturated and unsaturated hydrocarbons C1-C4 was studied. It was shown that the values of ideal selectivities for C 3 H 8 /CH 4 and n-C 4 H 10 /CH 4 gas pairs were higher than mixed-gas selectivities. This effect was explained by significant swelling of the membrane material in the hydrocarbon mixture, first of all, due to the presence of n-butane.
On the other hand, an opposite effect of the increase of the separation selectivity for gas mixtures is a character feature of high free volume glassy polymers (so-called polymers of intrinsic microporosity) [1,7]. This unique effect was first reported for poly(1trimethylsilyl-1-propyne) (PTMSP) by Srinivasan et al. [11]. It was shown that helium and nitrogen permeation was drastically reduced in the presence of a more strongly sorbing and/or more condensable SF 6 in the feed; the phenomenon was called "pore blocking" or "light gas rejection" effect. Thus, the ideal selectivity α 12 , calculated by the ratio of the permeabilities for pure gases 1 and 2 can differ rather significantly from the actual, mixed-gas selectivity α mix 12 . From the theoretical point of view, the deviation of the actual selectivity from the ideal selectivity is due to the complex, interdependent transport of mixture components through the polymer membrane. The reason for the marked reduction or enhancement of membrane selectivity can be designated as the coupling effect. This effect has sorption and diffusion constituents. Sorption coupling is associated with the intermolecular interaction of sorbed species and competition between unlike molecules for the limited number of available "sorption sites" in the polymer. Diffusion coupling is related to friction interaction between the mixture components in the membrane and changes the mobility of the components, slowing the transfer of lighter (faster) gas and accelerating the transfer of heavier (slower) gas.
There is a longstanding, perhaps not entirely correct, opinion that the gas mixture components permeate across rubbery polymers essentially independently of each other. For this reason, pure gas solubilities, diffusivities, and permeabilities can be used in mixed-gas calculations [12]. In the framework of linear irreversible thermodynamics, this means that the cross-terms in the flux differential equations can be neglected, in other words, the coupling of component fluxes can be ignored. Gas permeation properties of rubbery polymers were studied by a combination of the Maxwell-Stefan and Flory-Huggins theories in the works [13][14][15]. It has been obtained that the experimental data on mass transfer of CO 2 /C 2 H 6 [13,14] and CO 2 /CH 4 [15] mixtures in polyethylene oxide and polydimethylsiloxane, respectively, can be described with acceptable accuracy without regard to diffusional coupling. On the other hand, for pervaporation and gas separation by rigid glassy polymers, consideration of cross-terms in the flux equations seems necessary and, accordingly, the coupling effect cannot be neglected. It has been found that the friction interactions of penetrant-glassy polymer and penetrant-penetrant are of the same order [16,17].
It should be particularly emphasized that mixed-gas measurements need to modify and improve conventional and widely used experimental techniques [18,19]. Another problem is related to the accurate measurement of gas mixture sorption in polymer materials. This is a non-trivial and time-consuming procedure that has only been implemented for a few polymer membranes [20,21]. According to Genduso et al. [15], this problem represents a "bottleneck in gas transport analysis" in polymer membranes. The limited experimental data on diffusion and sorption of gas mixtures in polymeric membranes, on the one hand, restrain the formulation of rigorous theoretical approaches and, on the other hand, increases the importance of developing engineering predictive models to analyze the permeation properties of polymeric membranes.
The seminal work on the analysis of coupling in pervaporation (the description of pervaporation and gas separation are mathematically identical) was performed by Kedem [22] based on the friction-coefficient formalism. Kedem derived the flux expressions and considered two limiting cases: weak and strong coupling between the penetrant fluxes. However, the application of Kedem's model to the description of mixed-gas permeation revealed difficulties associated with the interpretation of component permeability coefficients and obtaining the correct expression for the separation factor.
In this work, an attempt was made to develop a simplified analytical model for the evaluation and possibly prediction of mixed-gas transport properties of polymer membranes based on available experimental data on single gases. The basic idea was to split the local permeability into two contributions, one of which depends on the friction between the penetrants, and the other contribution can be identified with pure gas permeability.

Basic Equations
The Maxwell-Stefan approach (equivalent to the friction-coefficient formalism) based on the concepts of linear irreversible thermodynamics is the classical method for the description of mass transfer [23][24][25]. If the isothermal process in the three-component system, penetrant 1/penetrant 2/polymer 3, is considered, then the Maxwell-Stefan equation (in the absence of external forces and shear forces) for each of the components can be described as follows: This equation defines the balance of the thermodynamic force, namely, the negative gradient of the component chemical potential along the z direction perpendicular to the membrane surface, exerted on the molecule of type i, and friction forces acting between this component and the molecules of other species. The friction force is the product of the friction coefficient f ij and the difference in velocity between species. The friction coefficient f ij is proportional to the concentration of component j, which for polymer systems is conveniently expressed in terms of the volume fraction in the membrane, i.e., Here the factor RT is the product of the gas constant and absolute temperature, ζ ij are the mutual friction (resistant) coefficients inversely proportional to the Maxwell-Stefan diffusion coefficients D ij . These coefficients are non-symmetric and satisfy the relation Note that one can also use the more traditional molar fraction-based diffusivities D ij . The relation between D ij and D ij is given by D ij = D ij V/V j where V is the mean molar volume of the system components). The average velocity of molecules u i is defined through the molar flux of the component N i (moles of species i per membrane unit area per unit time) and its molar concentration c i With the account for Equations (2) and (3) and assuming u 3 = 0 (stagnant polymer matrix), Equation (1) can be rearranged to yield N i as a function of the thermodynamic forces: where Here, D i is the local diffusion coefficient of species 1 or 2 in the membrane, φ i = c i V i is the volume fraction of component i. Note that D 12 /V 2 = D 21 /V 1 . In the absence of intermolecular friction between the components 1 and 2, i.e., at D 12 , D 21 → ∞ , Equations (4) and (5) are reduced to where A superscript "0" indicates the absence of diffusional coupling between fluid fluxes. The diffusion coefficient D 0 i referred to a polymer-fixed frame of reference.
Assuming ideal behavior in the gas phase which would be in equilibrium with the sorbed gas at position z, i.e., dµ i = RT d(ln p i ), Equations (4) and (5) can be rewritten as where is the local permeability coefficient, K i = c i /p i is the local solubility coefficient.

Kedem's Solution
Equations (9) and (10) are identical to those obtained by Kedem [22] (in Kedem's paper the coupling coefficient was written as Q = V 1 /D 21 = V 2 /D 12 ). In general, permeability coefficients depend on local partial pressures: Π i = Π i (p 1 , p 2 ). Assuming that the permeability coefficients are not partially pressure dependent, Equations (9) and (10) can be integrated to give where l is the actual thickness of a membrane and Without coupling ( B 1 , B 2 → 0 ) Equations (12) and (13) are reduced to where p f i and p p i are the partial pressures of the type i component in feed and permeate, respectively. In a strongly coupled system (B 1 , B 2 ≥ 5) Taking into account Equation (14), from the ratio of Equations (16) and (17) it follows that (y 1 /y 2 where i /p f is the molar fraction of component i in feed gas mixture, p f is the total feed pressure. A basic quantity of interest is the separation factor (SF) for a binary gas mixture defined as In terms of β 12 , Equation (18) reduces to β 12 = Π 1 x 2 /(Π 2 x 1 ). This relation includes permeability coefficients, which, according to Kedem's assumption, are constant. For the model to be practically usable, permeability Π i should preferably be identified with a parameter available from the experiment, such as the permeability of the pure component, and the permeability ratio Π 1 /Π 2 should be identified with the ideal selectivity (α 12 ). With this interpretation, the SF is expressed as (20) However, this expression, apart from the fact that it diverges at x 1 → 0 and tends to zero at x 1 → 1 , is incompatible with Kedem's important result for a strongly coupled system. This result (for negligible pressures on the permeate side) was indicated by Kedem herself and is as follows This conclusion can be obtained from the differential flux equation for one of the components. Let us consider Equation (4). In the strong coupling limit ( D 12 → 0 ) the first term in the right-hand side of this equation can be neglected. Moreover, as follows from Equation (6), the diffusion coefficient D 1 → D 12 /(c 2 V 2 ) in this limit, which leads to relation (21). According to the relation, in the limit of strong coupling, the SF tends to the ratio of the solubility coefficients. This inconsistency is evident in Equations (20) and (21). The reason for the discrepancy in the results for β 12 lies in Kedem's assumption concerning the constancy of the permeability coefficients when integrating Equations (9) and (10). A proposed approach, applying a partitioning of the inverse permeability coefficients into a term explicitly dependent on the intermolecular friction of the penetrants and a term determined by frictional interaction of the penetrant with the polymer matrix, is presented below.

Approximations and Flux Equations
Now, let us return to the differential flux Equations (9) and (10), and consider the permeability coefficients Π i (Equation (11)). Substituting Formulas (6) into Equation (11) gives the following expressions for the local permeability coefficients: Since the volume fraction of type i penetrant φ i = c i V i = K i p i V i , the above equations become: 1 where In the absence of coupling ( D 12 , D 21 → ∞ ), the first term in the right-hand side of Equations (22) and (23) vanishes, and the permeability coefficients Π i are reduced to the coefficients Π 0 i , which, as a first approximation, can be assumed to be independent of concentration and identified with the corresponding values for single gases. The relations (22) and (23) play a crucial role in the model development.
Substituting Equations (22) and (23) into Equations (9) and (10), respectively, give where ξ = z/l is dimensionless membrane thickness, the parameters B i are defined in Equation (15). These equations include the ratio of the local solubility coefficients. We make a simplifying assumption K 1 /K 2 = const that allows an analytical integration of the above differential equations. The assumption of constant K 1 /K 2 is fulfilled when the sorption of the mixture is consistent with the Henry model, as well as with the multicomponent Langmuir model. It seems reasonable to equate the ratio K 1 /K 2 with the ratio of solubility coefficients for pure gases α 12, S ≡ α S , i.e., Identification (27) means that, within the framework of this work, the deviation of the membrane selectivity from ideal selectivity is due only to frictional interaction of the gas penetrants 1 and 2. Namely, the coupling effect is based on the diffusional coupling only. In practice, there may be noticeable deviations K 1 /K 2 from the ideal solubility selectivity (see, e.g., [20,26]). Given that mixed-gas solubility data are far from always available and their experimental measurement is a time-consuming procedure, the approximation (27) is reasonably adequate for design calculations.
Using Equation (27) and the relation B i = y i B, Equations (25) and (26) can be rewritten as is the quantity that includes the total molar flux N of the gas mixture and 1-2 diffusivity D 21 , which is a measure of coupling. In the membrane literature, a dimensionless quantity of this type is often referred to as the "membrane Peclét number". Integration of the differential equation system (28) and (29) leads to the result (see Appendix A for some details): where Here r is the pressure ratio r = p p /p f (p f and p p denote the feed and permeate pressure, respectively). One can see that besides the partial pressure drop, the flux Equations (31) and (32) include an additional contribution ("coupling pressure") due to diffusion coupling. Note that this additional contribution has different signs for components 1 and 2, i.e., the coupling results in lower flux (compared to pure gas flux) for one component and higher flux for the other component. This is the apparent difference between the developed model and Kedem's approach, in which the additional "driven force" (complementary to ∆p i ) has a positive sign for both components (it is easy to show from the flux equations (10)). In the limit B → 0 (no coupling), Equations (31) and (32) are reduced in form to Equation (15). However, unlike Equation (15), these equations include permeability coefficients which are identified with corresponding quantities for individual gases. It follows from Equations (31) and (32) that Using the equality N i = y i N, one gets an expression for the total flux where In the above equations α 12 is the ratio of permeability coefficients for pure gases (the ideal selectivity) and α D is the ideal diffusivity selectivity. (By the way, note that if the penetrants have the same mobility, i.e., α D = 1, then, as seen from (35), the total flux coincides with the total flux of the gas mixture without coupling. This fact is a direct consequence of assumption (27)).
Any two of the three Equations (31), (32) and (35) can be used to obtain the equation for calculating the permeate composition. So, dividing Equation (31) by Equation (32) gives After solving this equation, the SF can be calculated.

Separation Factor
First, consider the trivial case of no coupling, i.e., B = 0. In this situation d 1 = 0 and Equation (37) reduces to the quadratic equation with the familiar result for the permeate composition [27]. Using the definition (19), the SF for a binary gas mixture with no coupling is given by It is evident that β 12 ≤ α 12 (equal sign for r = 0). Next, consider the situation of coupling between component fluxes (B > 0). If the downstream pressure is very small, then Equation (37) reduces formally to the quadratic equation with the following solution for the SF: It should be noted that this equation is indeed quadratic if one considers the Peclét number B as the parameter that does not depend on the permeate composition. In general, this is not the case (the Peclét number is a function of permeate composition, see Section 4.4 below). Nevertheless, Equation (39) is useful as an approximate starting point for further model development. In the case of nonzero downstream pressure (0 < r < 1), Equation (37) is reduced to a cubic equation for the SF, which is more convenient to solve numerically. Equation (39) implies that: (1) the SF is finite over the entire concentration range of the feed mixture (in distinction to expression (20)), in the limit x 1 → 0 one has β 12 | x 1 =0 = α 12 penetrants with the same diffusion mobility the SF is equal to the solubility selectivity; (3) β 12 = α S also if the Peclét number is much greater than one; (4) in the absence of coupling (B = 0) the SF is equal to the ideal selectivity. The last two facts determine the range of SF values: An important qualitative conclusion follows from (40): the coupling effect is the negative for the diffusivity selectivity membranes (α D > 1), i.e., real selectivity (mixture selectivity) is less than the ideal one. On the contrary, the coupling effect is positive for the solubility selectivity membranes (α D < 1), i.e., real selectivity is greater than the ideal one.

Component Permeabilities
The permeability coefficient of the mixture component i and mixed-gas selectivity are defined as Π mix From Equation (35) for the total flux, the component fluxes are With the relations y 1 = β 12 x 1 /(x 2 + β 12 x 1 ) and y 2 = 1 − y 1 , one gets As pointed out above in Equation (40), the SF in the limit of large Peclét numbers is equal to the solubility selectivity. Using this result in Equation (43) and for simplicity, setting r = 0, the limiting permeability coefficients are This suggests that α mix 12 ≡ Π mix 1 /Π mix 2 = α S in the case of a strongly coupled system. It can be shown (see Appendix B) that the expressions (44) are equivalent to the relations previously obtained by Krishna and van Baten [28] in the model description of gas mixture permeation across micro-and mesoporous membranes.
The implications of the expressions (44) are obvious. If the preferentially permeating component (component 1) of the gas mixture is more mobile, i.e., α D > 1, then its permeability Π mix 1 will be lower compared to the pure gas permeability Π 0 1 ; the mixed-gas permeability of component 2 will, on the contrary, be higher than in pure gas conditions. If the preferentially permeating component of the gas mixture is less mobile, i.e., α D < 1, then under conditions of a gas mixture it penetrates through the membrane faster than in a pure state; the second, more diffusion-mobile component, penetrates slower compared to corresponding pure gas. The situation described is reminiscent of Le Chatelier's principle in chemical thermodynamics. To summarize, one can write Note that conditions (45) are obtained for the situation of strong coupling (large Peclét numbers) and assuming downstream pressure is very small. Calculations show that inequalities (45) are also valid under less stringent constraints.
Examples of the SF and component permeability calculations are shown in Figures 1 and 2, respectively. Consider first the case of a sorption-selective membrane, for which α D < 1 and α S > α 12 > 1. The SF increases with the degree of coupling (i.e., the Peclét number) (Figure 1a) and tends to the value of the ideal solubility selectivity. If the pressure ratio r = 0, then the SF exceeds the value of ideal selectivity α 12 (=15 in this example). At non-zero permeate pressure, the SF also grows monotonically with increasing the Peclét number. In this case, two intervals can be distinguished: at small Peclét numbers (B less than about 1) β 12 < α 12 , and at appreciable Peclét numbers (B more than about 1) the SF even exceeds the corresponding parameter at zero permeate pressure. In the case of a diffusion-selective membrane (α D ∼ α 12 > 1 and α S ∼ 1), the SF is less than the ideal selectivity and decreases with an increase in the Peclét number (Figure 1b). In addition, an increase in downstream pressure (from 0 to 0.1p f ) leads to an additional decrease in the SF.
case of a sorption-selective membrane: the permeability of fast component 1 decreases, and the permeability of slow component 2 increases with Peclét's number (Figure 2b). At the same time, an increase in downstream pressure has almost no effect on the mixed-gas permeabilities values.

Explicit Form of the Peclét Number
Above, the SF was estimated as a function of the Peclét number using the ideal selectivity and solubility selectivity as input parameters. The effect of coupling was treated formally, in the sense that the Peclét number was assumed to be a parameter independent of feed and permeate composition. A more rigorous approach, not limited to the description of experimental data, should take this dependence into account and provide an explicit expression for the Peclét number. In addition, it is important to evaluate the case of a sorption-selective membrane: the permeability of fast component 1 decreases, and the permeability of slow component 2 increases with Peclét's number (Figure 2b). At the same time, an increase in downstream pressure has almost no effect on the mixed-gas permeabilities values.

Explicit Form of the Peclét Number
Above, the SF was estimated as a function of the Peclét number using the ideal selectivity and solubility selectivity as input parameters. The effect of coupling was treated formally, in the sense that the Peclét number was assumed to be a parameter independent of feed and permeate composition. A more rigorous approach, not limited to the description of experimental data, should take this dependence into account and provide an explicit expression for the Peclét number. In addition, it is important to evaluate the The behavior of the SF will become clearer when considering the dependence of relative permeability coefficients Π i = Π mix i /Π 0 i on Peclét numbers (Figure 2). For a sorption-selective membrane, component 1 is less mobile than component 2 (α D < 1). As evident from Figure 2a, the mixed-gas permeability of component 1 exceeds pure gas permeability, while the mixed-gas permeability of component 2, on the contrary, is lower than its permeability in pure gas conditions. Namely, slow gas 1 moves faster in the mixture, and fast gas 2, on the contrary, reduces its diffusion mobility. Everything occurs according to the above inequality (45) at α D < 1. Calculations have shown (dashed lines in Figure 2a) that an increase in downstream pressure leads to a larger deviation of the permeability coefficients from the ideal (pure gas) values. In the case of a diffusion-selective membrane, mixed-gas permeabilities behave in the opposite way to the case of a sorption-selective membrane: the permeability of fast component 1 decreases, and the permeability of slow component 2 increases with Peclét's number (Figure 2b). At the same time, an increase in downstream pressure has almost no effect on the mixed-gas permeabilities values.

Explicit Form of the Peclét Number
Above, the SF was estimated as a function of the Peclét number using the ideal selectivity and solubility selectivity as input parameters. The effect of coupling was treated formally, in the sense that the Peclét number was assumed to be a parameter independent of feed and permeate composition. A more rigorous approach, not limited to the description of experimental data, should take this dependence into account and provide an explicit expression for the Peclét number. In addition, it is important to evaluate the range of Peclét numbers characteristic of mixed-gas transport through sorption-or diffusion-selective membranes.
Let us now rewrite Equation (30) for the Peclét number as the following: where N is the dimensionless total permeate flux. The explicit expression of N follows from Equation (35): The expression for B contains the unknown Maxwell-Stefan diffusivity D 21 which is an inverse friction coefficient between penetrants 1 and 2. Given that Maxwell-Stefan diffusivities are not accessible from the experiment, various interpolation relations are used to predict those [29]. In this paper the simple (geometric mean) approximation is used: Note that this expression for 1-2 diffusivity is similar (and even simpler) to the interpolation formula used by Krishna [14] when modeling water/ethanol pervaporation across a polyimide membrane. Substituting of Equation (48) into Equation (46) gives where S 2 is the experimental solubility coefficient of component 2 in units cm 3 (STP)/(cm 3 ·atm), the feed pressure expressed in atm, a correction factor k ∼ 1 was introduced into the equation to compensate for the approximate nature of equality (48). Let us estimate the value of B depending on the composition of the mixture to be separated. For the sake of simplicity, let the pressure ratio r = 0. It follows from Equation (47) that the total dimensionless flux varies from 0 to α 12 when the mixture composition changes from x 1 = 0 (the pure component 2) to x 1 = 1 (the pure component 1). According to Equation (49), the value of B increases with the upstream pressure and molar volumes of penetrants. At fixed upstream pressure Equation (49) becomes For instance, if to take p f = 10 atm, √ V 1 V 2 ≈ 85 cm 3 /mol, and k = 1, then const ≈ 0.04. Since component 1 is the preferentially permeating component (α 12 > 1), the coupling effect is maximum in the limit of infinite dilution of component 2 in the feed mixture. In addition, as can be seen from Equation (50), the value of B decreases with increasing diffusivity selectivity. Therefore, one would expect the coupling effect to be weaker for membranes with diffusion selectivity than for membranes with sorption selectivity.
Diffusivity selectivity glassy polymers (polyimides, microporous polymers), which are used to separate olefins and paraffins C 2+ , have almost no sorption selectivity (α S ∼ 1), and the ideal selectivity values are α 12 ∼ α D = 2 − 30 [8,30]. If one accepts α 12 = α D = 10 and S 2 = 10 cm 3 (STP)/(cm 3 ·atm) (as the average solubility of propane in aromatic polyimides) for a rough estimate, then from Equation (50) one obtains that the Peclét numbers vary from 0.2 to 1.9. For solubility selectivity glassy polymers (polyacetylenes, polynorbornenes, PIM-1) used for C 2+ /methane separations, the diffusion selectivity is less than 1, and the sorption selectivity is in a wide range from 5 to~10 3 [8,31]. For typical values α 12 = 10, α D = 0.2 and S 2 = 4 cm 3 (STP)/(cm 3 ·atm) (as the solubility of methane in PTMSP at 25 • C) one gets from Equation (50) the range of Peclét numbers from 0.4 to 3.6. Thus, it can be stated that in the hydrocarbon separation by glassy polymers, the Peclét number reaches values of the order of unity (notice that the above estimates were made at a total feed pressure of 10 atm).
Let us list the basic equations for the mixed-gas selectivity and permeability calculation within the framework of the developed model. Permeate composition is calculated using Equation (37) with the expression (49) for the Peclét number. The relationship of permeate composition with the SF is given by Equation (19). Penetrant fluxes are then calculated using Equations (31) and (32) or, equivalently, via the expressions (42). For negligible pressures of the permeate side, computations are simplified. In this case, the SF can be found by Formula (39). In doing so, the Peclét number can be considered in the first approximation as a model parameter, or it can be estimated using the explicit expression (49). Input parameters for the proposed model are the pure gas permeabilities and solubility coefficients.

Comparison of the Model with Experimental Data
The results of the model calculations were validated against the experimental data for two examples of hydrocarbon separation by the glassy polymers. The first example is nbutane/methane separation through high free volume disubstituted polyacetylene-poly(1trimethylsilyl-1-propyne) (PTMSP). The second example is propylene/propane separation through polyimide 6FDA-TrMPD. These examples illustrate the transport properties of sorption-and diffusion-selective membranes, respectively. Pure gas values of selectivity and solubility coefficients were used as input parameters for calculations ( Table 1). Results of the model calculations and their comparison with mixed-gas experimental data are presented in Figure 3. Table 1. Permeabilities and solubility coefficients (Henry's constants), and their ratio for n-C 4 H 10 /CH 4 (2 mol %/98 mol %) and C 3 H 6 /C 3 H 8 (equimolar mixture) in PTMSP and 6FDA-TrMPD membranes, respectively. Diffusivity selectivities are shown in the last column.

Polymer
Π 1 (barrer) light gas CH4 diffusion in free volume elements of "microporous" PTMSP is hindered, blocked by the sorbed molecules of heavier, condensable n-C4H10 [11,35]. Along with this basic reason, one should also take into account the effect of competitive sorption, which can lead to a change (increase) of solubility selectivity [20,33]. In terms of the developed model, the blocking of methane transport by n-butane is a consequence of frictional interaction between penetrants, in other words, the blocking mechanism is a manifestation of coupling between methane and n-butane fluxes. through PTMSP membrane versus n-C4H10 feed molar fraction; (b) Mixed-gas C3H6/C3H8 selectivity and relative permeabilities of the mixture components through polyimide 6FDA-TrMPD versus C3H6 feed molar fraction. The results of selectivity and permeability calculations are given for pressure ratio r = 0 (solid lines) and r = 0.05 (dashed lines). Experimental selectivities and permeabilities for 2 mol % n-C4H10/98 mol % CH4/PTMSP and 50 mol % C3H6/50 mol % C3H8/6FDA-TrMPD systems are shown by the symbols.

Conclusions
An approximate model for the mixed-gas selectivity and permeability of non-porous and microporous polymeric membranes was developed. This model is the modification of Kedem's friction-coefficient approach [22] to flux coupling in pervaporation. The central argument of the model is to decompose the inverse local permeability of penetrants into two terms. The first one is the inverse permeability of pure gas, and the second one is proportional to the friction between penetrants. As a result, under certain simplifying assumptions (constant values of local pure gas permeability and solubility selectivity), analytical expressions for permeability and selectivity of polymeric membranes in mixed-gas conditions were derived. It is shown that the use of simple (geometrical mean) approximation for the Maxwell-Stefan penetrant 1/penetrant 2 diffusivity allows considering the developed model as a predictive one if pure gas permeability and solubility coefficients are available as input. It was found that the coupling between penetrant fluxes causes the membrane selectivity to deviate from pure gas selectivity. For diffusion-selective membranes, there is a reduction in selectivity, while for sorption-selective membranes there is an enhancement in selectivity with respect to pure gas selectivity. A quantitative measure of the coupling effect in the model is the membrane Peclét number. As estimates have shown, the Peclét number is less or of the order of unity in hydrocarbon separation by glassy polymers. In the limit of large Peclét numbers, the derived expressions have the same final form as those obtained earlier by Krishna and van Baten [28] in modeling the permeability of gas mixtures across the inorganic microporous ma-  (a) Mixed-gas n-C 4 H 10 /CH 4 selectivity and relative permeabilities Π i = Π mix i /Π 0 i of the mixture components through PTMSP membrane versus n-C 4 H 10 feed molar fraction; (b) Mixed-gas C 3 H 6 /C 3 H 8 selectivity and relative permeabilities of the mixture components through polyimide 6FDA-TrMPD versus C 3 H 6 feed molar fraction. The results of selectivity and permeability calculations are given for pressure ratio r = 0 (solid lines) and r = 0.05 (dashed lines). Experimental selectivities and permeabilities for 2 mol % n-C 4 H 10 /98 mol % CH 4 /PTMSP and 50 mol % C 3 H 6 /50 mol % C 3 H 8 /6FDA-TrMPD systems are shown by the symbols.
In the case of PTMSP, the n-C 4 H 10 /CH 4 selectivity is noticeably greater than the ideal value (=3.5) and grows with n-butane concentration in feed (Figure 3a (top)). The component permeabilities (in units of pure gas values), on the contrary, decrease with the increasing concentration of n-butane (Figure 3a (bottom)). Light gas CH 4 permeability in mixed-gas conditions decreases as compared to the pure gas value (Π mix CH 4 < Π 0 CH 4 ), whereas the situation is the opposite for heavier n-C 4 H 10 (Π mix C 4 H 10 > Π 0 C 4 H 10 ). This behavior is consistent with the second inequality (45).
In the case of polyimide 6FDA-TrMPD, the mixed-gas C 3 H 6 /C 3 H 8 selectivity is almost half that of pure gas selectivity and decreases with increasing concentration of preferentially permeating component (propylene) in feed (Figure 3b (top)). The permeability of both components increases with the propylene concentration in the feed (Figure 3b (bottom)). Permeability behavior in mixed-gas conditions follows the first inequality (45): for more mobile gas C 3 H 6 permeability decreases compared to pure gas one, while for less mobile gas C 3 H 8 there is an increase in permeability compared to a pure gas one.
Generally, one can state a semi-quantitative agreement with the experiment for both penetrants/polymer systems. As can be seen in Figure 3, the agreement between the calculated and experimental selectivities improves at non-zero permeate pressure. The calculation results depend on the Peclét number. It was evaluated using Equation (49) with fixed parameters p f = 10 atm, √ V 1 V 2 = 85 cm 3 /mol, and k = 1. Given the crude approximation (48) for 1-2 diffusivity, the developed model shows more than satisfactory performance.
In closing, it should be pointed out the following. It has been known since [35] that the mixed-gas n-C 4 H 10 /CH 4 selectivity of PTMSP membranes significantly exceeds the ideal selectivity, by a factor of 6-18 according to various experimental data [31]. The increase in selectivity is associated with a marked (3-9 times) decrease in CH 4 permeability under mixed-gas conditions compared to pure gas conditions [31]. The physical reason for depressing the CH 4 permeability is commonly attributed to a blocking mechanism: light gas CH 4 diffusion in free volume elements of "microporous" PTMSP is hindered, blocked by the sorbed molecules of heavier, condensable n-C 4 H 10 [11,35]. Along with this basic reason, one should also take into account the effect of competitive sorption, which can lead to a change (increase) of solubility selectivity [20,33]. In terms of the developed model, the blocking of methane transport by n-butane is a consequence of frictional interaction between penetrants, in other words, the blocking mechanism is a manifestation of coupling between methane and n-butane fluxes.

Conclusions
An approximate model for the mixed-gas selectivity and permeability of non-porous and microporous polymeric membranes was developed. This model is the modification of Kedem's friction-coefficient approach [22] to flux coupling in pervaporation. The central argument of the model is to decompose the inverse local permeability of penetrants into two terms. The first one is the inverse permeability of pure gas, and the second one is proportional to the friction between penetrants. As a result, under certain simplifying assumptions (constant values of local pure gas permeability and solubility selectivity), analytical expressions for permeability and selectivity of polymeric membranes in mixed-gas conditions were derived. It is shown that the use of simple (geometrical mean) approximation for the Maxwell-Stefan penetrant 1/penetrant 2 diffusivity allows considering the developed model as a predictive one if pure gas permeability and solubility coefficients are available as input. It was found that the coupling between penetrant fluxes causes the membrane selectivity to deviate from pure gas selectivity. For diffusion-selective membranes, there is a reduction in selectivity, while for sorption-selective membranes there is an enhancement in selectivity with respect to pure gas selectivity. A quantitative measure of the coupling effect in the model is the membrane Peclét number. As estimates have shown, the Peclét number is less or of the order of unity in hydrocarbon separation by glassy polymers. In the limit of large Peclét numbers, the derived expressions have the same final form as those obtained earlier by Krishna and van Baten [28] in modeling the permeability of gas mixtures across the inorganic microporous materials. In this limiting case, the mixed-gas selectivity is equal to the value of the ideal solubility selectivity.
The model was tested for the separation of n-butane/methane and propylene/propane mixtures through PTMSP and polyimide 6FDA-TrMPD membranes, respectively. These two examples illustrate the hydrocarbon separation with sorption-selective (PTMSP) and diffusion-selective (6FDA-TrMPD) membranes. A rough quantitative agreement of the model results with the experimental data was obtained. Thus, the model evaluation of mixed-gas selectivity using pure gas permeation data might be considered as a first estimate of the applicability of membranes to separate a target gas mixture.

Conflicts of Interest:
The authors declare no conflict of interest.