Next Article in Journal
Effects of Owinema Bio-Preparation on Vermicomposting in Earthworm Ecological Boxes
Next Article in Special Issue
Glassy PEEK-WC vs. Rubbery Pebax®1657 Polymers: Effect on the Gas Transport in CuNi-MOF Based Mixed Matrix Membranes
Previous Article in Journal
ECE/SAE Dual Functional SuperPin Plus Curved Reflex Reflector by Use of New Structured Corner Cubes
Previous Article in Special Issue
Crosslinked Facilitated Transport Membranes Based on Carboxymethylated NFC and Amine-Based Fixed Carriers for Carbon Capture, Utilization, and Storage Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mass Transfer Through Graphene-Based Membranes

1
Institute of Chemical Process Fundamentals of the CAS, v.v.i., Rozvojova 135, 165 02 Prague 6, Czech Republic
2
Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 166 28 Prague 6, Czech Republic
3
Institute of Environmental Technology VSB-TUO, 17. listopadu 15/2172, 708 33 Ostrava-Poruba, Czech Republic
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(2), 455; https://doi.org/10.3390/app10020455
Submission received: 18 November 2019 / Revised: 29 December 2019 / Accepted: 3 January 2020 / Published: 8 January 2020
(This article belongs to the Special Issue Preparation, Characterization and Modelling of Advanced Membranes)

Abstract

:
The problems related to the transport of gases through nanoporous graphene (NG) and graphene oxide (GO) membranes are considered. The influence of surface processes on the transport of gas molecules through the aforementioned membranes is studied theoretically. The obtained regularities allow finding the dependence of the flux of the gas molecules passing through the membrane on the kinetic parameters which describe the interaction of the gas molecules with the graphene sheets. This allows to take into account the influence of external fields (e.g., resonance radiation), affecting the aforementioned kinetic parameters, on the transport of gas molecules through the membranes. The proposed approach makes it possible to explain some experimental results related to mass transfer in the GO membranes. The possibility of the management of mass transfer through the NG and GO membranes using resonance radiation is discussed.

1. Introduction

The scope of our work is to describe the transport properties with existing models without any additional data. Moreover, the gas separation mechanism for GO membranes is unclear due to inconsistent permeation and separation results reported in the literature.
It is known that graphene-based materials are very efficient in many areas of chemical engineering in particular for the elaboration of new type of batteries for modern electronics [1] and for creation of new types of membranes for the separation of gas mixtures [2]. Graphene is a two-dimensional one-atom-thick sheet of carbon atoms packed in the honeycomb crystal lattices, it is characterized by unique physicochemical properties in particular the high chemical stability, mechanical strength, thermal conductivity and charge mobility [1,2]. It was shown that pristine graphene is impermeable in regard to smallest molecules [2].
A graphene sheet with nano and sub-nano pores, known as nanoporous graphene (NG), is used for the size-selective separation of gas mixtures based on the molecular sieving effect [2]. The permeability and selectivity of NG membranes are higher than the aforementioned values for polymer membranes. Ion bombardment and electron beam exposure are often used to prepare pores in the graphene sheets. As indicated in [3], molecules adsorbed on the graphene sheet and participating in surface diffusion [4] are more likely to find the pores in NG membranes than the gas molecules. Thus, surface diffusion can be one of the main processes affecting the transport of gas molecules through NG membranes [5,6,7]. It should be noted that the problem related to the management of the transport of gas molecules using the porous graphene is very actual (e.g., for removing carbon dioxide from the gas mixtures) [8]. The core factor affecting permeance of NG membranes is the pore size. However, as indicated in the aforementioned paper, it is difficult to fabricate defect-free porous membranes in a large scale that have pores matching the kinetic diameters of common gas molecules, ranging from 2.5 to 4.0 Å. The combined use of the ionic liquid and porous graphene for more effective separation of the gas mixtures is considered in [8]. It is demonstrated from classical molecular dynamics simulations that a monolayer of ionic liquid on the NG membrane can be used to control the pore size of porous graphene. However, this method requires modification of the original membrane. In the present paper is discussed the new method to control the gas transport through NG membranes using external fields (e.g., resonance laser radiation) leading to excitation of gas molecules and changing their interaction with the graphene sheet. In this case do not need any modifications of the original membrane. The effect of resonance radiation on the separation of the gas mixtures passing through the dense membranes was considered in [9].
Graphene oxide (GO) can be considered as the oxidized form of graphene [2]. GO membranes contain stacked GO sheets as well as vacancy defects and functional groups distributed in carbon lattice, these membranes can be cheaply produced by oxidation and exfoliation of graphite [2]. The problems related to preparation of GO membranes and mass transfer in them were discussed in particular in [10,11,12]. In [10] it was shown that self-standing GO membranes demonstrate outstanding performance for hydrogen separations, including mixtures of H2/CO2 and H2/hydrocarbons. An efficient spray-evaporation method was explored in [11] to optimize performance of GO membranes. The very interesting experimental results related to the transport of gases such as CH4, N2, CO2, H2 and He through GO membranes are given in [12]. These results can be used in a study of the separation of gas mixtures and in a development of new models for mass transfer through GO membranes. The model for the gas transport through GO membranes proposed in the present paper allows explaining some experimental results given in the above-mentioned paper.
Below, we consider the possible influence of the surface processes (including the features of the gas molecule scattering by graphene sheet and surface diffusion of adsorbed molecules) on the transport of gases through the NG and GO membranes. This allows for the dependence to be found of the flux of gas molecules passing through the membranes on the kinetic parameters which describe the interaction of the gas molecules with the graphene sheets. The influence of external fields (e.g., resonance radiation), affecting the aforementioned kinetic parameters, on mass transfer in the graphene-based membranes is discussed. The core idea of this paper consists in the search of new ways to manage mass transfer in the porous graphene and graphene oxide membranes by means of external fields (e.g., resonance radiation) affecting the above-mentioned kinetic parameters.

2. The Theoretical Description of Mass Transfer Through NG Membranes

Firstly, let us consider the transport of molecules of the gas mixture through the NG membrane. The flux density of the gas molecules of component i passing through the NG membrane jir can be presented as [5].
j ir = j is + j id ,
where the value of jis is related to adsorption of the gas molecules of component i on the graphene sheet, their diffusion along the graphene sheet to reach the pores and the penetration through the two-dimensional pores to the other membrane side, the value of jid is related to the direct penetration of gas molecules of component i through the pores of the NG membrane. We assume, similarly to [3], that the rates of adsorption and desorption of gas molecules are much higher than the diffusion rate of adsorbed molecules along the graphene sheet and the penetration rate of the molecules through the pores to the other side of the NG membrane, as well as that the adsorbed phase on the graphene sheet is in equilibrium with the free gas near the considered side of NG membrane. We suppose also that the number density of adsorbed gas molecules on the graphene sheet is sufficiently small and we can use the Henry law for adsorption of gas molecules, as well as that the size of the pores is greater than the size of the gas molecules and these molecules can penetrate through the two-dimensional pores in the NG membrane.
The value of jis is related to the frequency of encounters of adsorbed molecules of component i with the two-dimensional pores on the unit surface of the graphene sheet Qip and the probability of transfer (translocation) of the molecules after encounters with the two-dimensional pores onto the other side of the NG membrane (for the sake of simplicity, we assume that the two-dimensional pores are round and have the same diameter). In the case of localized mobile adsorption, for the molecule movement from one site to an adjacent one, the molecule must possess sufficient energy to pass over the energy barrier between the two sites [4]. If this barrier is low compared to the average energy of the molecule, surface diffusion of the adsorbed gas molecules on the graphene sheet can be considered on the basis of the model of the two-dimensional gas [6]. The value of Qip can be found similarly to the method used in the theory of heterogeneous nucleation [13,14]. In the case of localized mobile adsorption (when surface diffusion is realized by jumping mechanism), the value of Qip in view of [13] can be written as
Q ip = π d n p 1 4 n ia a 0 ν exp E is k T = ε p d n ia a 0 ν exp E is k T .
where d is the diameter of the two-dimensional pore, np is the number density of the two-dimensional pores on the graphene sheet, nia is the number density of the adsorbed gas molecules of component i on the graphene sheet, a0 is the jump distance for the adsorbed molecule (lattice parameter), νʹ is the vibration frequency of the adsorbed molecule (the values a0 and νʹ are assumed usually to be constant), Eis is the activation energy for surface diffusion of adsorbed molecules of component i, k is the Boltzmann constant, T is the temperature (we assume further the system to be isothermal), εp = npsp is the porosity of the NG membrane, where sp = πd2/4 is the area occupied by the two-dimensional pore. The flux density of component i passing from one side of the graphene sheet onto its another side jia can be written as
j ia = Q ip exp E it k T ,
where Eit is the effective activation energy for the molecule translocation from one side of the graphene sheet onto its another side. The rate constants for different steps in translocation of molecules from one side of the NG membrane onto its another side are considered in [7]).
In the case of the model of the two-dimensional ideal gas for the adsorbed molecules, the value of Qip in view of [14] is given by
Q ip = π d n p n ia π π k T 2 m i 1 / 2 = 4 ε p n ia d k T 2 π m i 1 / 2 ,
where mi is the mass of molecules of component i.
It is seen from Equations (2) and (4) that, at the constant εp, the value of Qip increases with a decrease in the pore diameter d. The flux density of molecules of the component i related to the direct flux of molecules passing from the gas phase near one side of the NG membrane into the gas phase near another membrane side (the gas phase pathway) jid according to [7] can be presented as
j id = ε p p i 2 π m i k T 1 / 2 G i , gas ,
where pi is the partial pressure of component i of the gas mixture, the value of Gi,gas describes the correction to the Knudsen equation for the gas flux of molecules of component i passing through the sub-nanometer two-dimensional pore (this correction is related to the sufficiently small value of d/dm, where dm is the kinetic diameter of gas molecules). The value of Gi,gas according to [7] can be written as
G i , gas = exp Δ E i , gas k T ,
where ΔEi,gas is the energy barrier height relative to the gas phase in the transport of gas molecules of component i through the two-dimensional pore. It is seen from Equations (5) and (6) that the value of jid decreases with an increase in ΔEi,gas.
A decrease in the flux density of gas molecules jid passing through the NG membrane on the basis of the steric consideration is discussed in [5]. The reduction in the value of Gi,gas in Equation (5) and accordingly in the value of jid is related in this case to the fraction of molecule trajectories that intersect the NG membrane within the pore diameter but cannot permeate the pore. It should be noted that Gi,gas → 1 with an increase in d/dm. The resulting flux density of gas molecules passing through NG membrane jir for the case of localized mobile adsorption in view of Equations (2), (3), and (5) is given by
j ir = ε p 1 d a 0 ν n iaf n iap exp E is + E it k T + p i f 2 π m i k T 1 / 2 p i p 2 π m i k T 1 / 2 G i , gas ,
where pif and pip are the partial pressures of component i at the feed and permeate sides of the NG membrane, niaf is the number density of adsorbed molecules of component i on the feed side of the membrane, niap is the number density of adsorbed molecules of component i on the permeate side of the NG membrane.
The number density of adsorbed molecules of component i nia on the graphene sheet at the aforementioned assumptions can be written as [4]
n ia = α i τ i p i 2 π m i k T 1 / 2 = α i τ i n i v i 4 .
where αi is the sticking coefficient of molecules of component i, τi is the adsorption time for the molecules of the component i (τi = τ0 exp(Eia/(kT)), where the value of the pre-exponential factor τ0 is assumed usually to be constant, Eia is the adsorption energy), vi = (8kT/(πmi))1/2 is the mean thermal velocity of gas molecules of component i, and ni is the number density of molecules of component i in the gas phase.
In the case when the adsorbed phase on the graphene sheet can be modeled by the ideal two-dimensional gas the resulting flux density of gas molecules passing through the NG membrane jir in view of Equations (3)–(5) can be written as:
j ir = ε p 4 d ( n iaf n iap ) k T 2 π m i 1 / 2 exp E it k T + p i f 2 π m i k T 1 / 2 p i p 2 π m i k T 1 / 2 G i , gas .
Equations (7)–(9) show that the influence of resonance (e.g., laser) radiation on the porous graphene membrane allows the transport of gas molecules through the membrane and separation of gas mixtures to be managed due to a change in the radiation field of the values of αi and τi (these values can be different for excited by radiation and unexcited molecules [9]). A decrease in the values of αi and τi for molecules of foreign gases under the effect of resonance radiation can lead to the purification of the graphene sheets from the adsorbed molecules of foreign gases and to a decrease in their transport through NG membranes.
It is seen from Equations (7) and (9) that heating of the gas-membrane system by radiation leading to an increase in the system temperature also affects the transport of gas molecules through the NG membrane and accordingly the separation of gas mixtures passing through the membrane.

3. Developing of the Model for Mass Transfer Through GO Membranes

Graphene oxide is used for the production of highly effective membranes for the purification of water and separation of gas mixtures [2]. The GO membrane can be considered as a composite, comprised of pristine and oxidized parts of the system consisting of graphene sheets. The GO membranes have two types of pores: the slit-like pores between the GO sheets and the small holes in the sheets related to defects in the GO sheets [2]. The empty space between the non-oxidized regions of GO monolayered sheets can be assumed lying in the range 0.4–0.9 nm. Distance between the point-like defects in the graphene sheets can be estimated as 10.0–13.3 nm [11]. The proposed in the present paper model for the description of mass transfer in GO membranes is related to following main points (which can serve as a kind of the scheme to describe the theoretical model and its novelty): the transport of gas molecules in GO membranes is realized mainly through the slit-like pores formed by graphene sheets. Mass transfer in these pores can be affected by the processes such as adsorption and desorption of gas molecules, surface diffusion of adsorbed molecules on the pore walls, and the gas molecule scattering by the pore walls. The external fields (e.g., resonance laser radiation) can affect the parameters characterizing the aforementioned processes that allows to develop new methods controlling mass transfer in GO membranes without modification of the membranes. The proposed model gives also the ability to explain some experimental results given in [12].
The effective diffusion coefficient for the activated Knudsen diffusion (or in other words translational diffusion) of component i of a gas mixture in the nanoscale slit-like pore of the GO membrane with the flat parallel pore walls Die according to [15] can be written as
D ie = 1 3 H ie v i exp E id k T ,
where Hie = H − dim is the effective accessible distance between the pore walls in relation to molecules of component i (the accessible pore width [15]), dim is the kinetic diameter of molecule of component i, H is the real width of the slit-like pore, Eid is the energy barrier for the activated Knudsen diffusion, which depends on the pore width. It is worth noting that the pre-exponential factor in Equation (10) in the case of the cylindrical pore is equal to dievi/3, where die is the effective diameter of the cylindrical pore in relation to molecules of the component i (die = dc − dim, where dc is the real diameter of the cylindrical pore [16]). It follows from the aforementioned that the greater the molecule size (the kinetic diameter dim) the smaller the effective pore diameter and respectively the diffusion coefficient. The value of Eid in Equation (10) decreases with a decrease in the value of dim/H. As indicated in [12], the activation energy for the diffusion of gas molecules in the cylindrical pore is small (close to zero) at dim/dc < 0.6 (in the case of the slit-like pore we should to use the value of dim/H instead of dim/dc). It is worth to note that Equation (10) corresponds to the case of diffuse molecule scattering by the pore walls. For the sake of simplicity, further we assume that the value of Eid in Equation (10) is small (close to zero) and Hie = H. Under these assumptions and Equation (10), we can estimate the diffusion coefficient of the classic Knudsen diffusion for the slit-like pore as
D i , Kn = 1 3 H v .
The experimentally observed excess of the measured value of the flux of light gas molecules passing through the carbon nanotubes in relation to the value of the flux of gas molecules according to the classic Knudsen diffusion was discussed in [17]. As indicated in [17], this effect can be related to specular reflections of gas molecules from the pore walls. The Maxwell diffuse-specular model of gas molecule scattering leads to the factor ψi = (2 − fi)/fi in the value of the diffusion coefficient in the classic Knudsen diffusion [18]. Here fi is the coefficient characterizing the probability of diffuse scattering of gas molecules of component i. The value 1 − fi = fi,sp describes the probability of specular reflection of gas molecules of component i. The probability of specular reflections of molecules from the pore walls depends in particular on the surface roughness, the relationship between the masses and sizes of molecule (atom) of the gas phase and the dense phase (pore walls) and on the adsorption energy of gas molecules.
Let us consider the dependence of the value of fi,sp on the mass of molecules of a gas phase in the simplest case of a smooth and clean (without adsorbed molecules) surface. The probability of elastic (which is assumed to be specular) reflection fi,sp of the gas molecule of component i according to [19] can be presented as
f i , sp = exp E ic k T .
where Eic is the specific energy of the interaction of the gas molecules of component i with the pore walls, which depends on the adsorption energy Eia and can be estimated for light gas particles as
E ic 4 m i m w E ia ,
where mw is the mass of molecules (atoms) of the dense phase (material of the pore walls).
It is seen from Equations (12) and (13) that the smaller the values of mi and Eia at the constant value of mw, the greater the probability of elastic reflections of gas molecules. In particular, this fact can lead to a sufficiently large probability of specular reflections of light gas particles such as H2 and He from the walls of the carbon nanotubes. The similar regularities may apparently take place in the transport of gases through the slit-like pores of GO membranes.
In the backscattering scheme, it is supposed that a fraction fi of molecules of component i is scattered according to the cosine law (diffuse scattering) and a fraction 1 − fi is scattered in the direction opposite to the velocity of incident molecules and the value of ψi = fi/(2 − fi) [20]. The backscattering may occur in the case of the roughness of the pore walls and adsorption of molecules of a foreign gas on the pore walls. It was shown in [21] that adsorption of one component of a gas mixture on the pore walls can influence the scattering of the second component of a gas mixture that affects the flux of the gas molecules of the second component passing through the channel (pore). The problems related to the effect of adsorption, surface diffusion, specular reflection, and backscattering of the gas molecules on the gas flux passing through a long cylindrical channel (pore) in the free molecular regime of gas flow are considered in [20,21,22,23].
It is worth noting that surface diffusion of molecules adsorbed on the walls of the slit-like pores of the GO membrane can also increase the flux of molecules passing through the pores. Let us estimate the influence of surface diffusion of molecules adsorbed on the walls of the slit-like pore of the GO membrane on the total flux of molecules passing through the pore. For the sake of simplicity, we assume that the total flux of molecules of component i passing through the slit-like pore can be presented as a sum of the flux of the gas molecules and the surface flux of adsorbed gas molecule on both walls of the pore (the additive approximation). Taking into account Equation (8) and the aforementioned assumptions, the total flux of molecules passing through some effective cross-section of the slit-like pore σe = Hl, where l is the length of the line on the flat pore wall which is directed perpendicularly to the direction of mass transfer, can be written as
J it = H l ψ i D i , Kn d n i d x 2 l D is d n ia d x = H l D i , Kn ψ i + 1 2 H D is v i τ i α i D i , Kn d n i d x ,
where x is the longitudinal coordinate directed along the slit-like pore in the direction of mass transfer and Dis is the coefficient of surface diffusion which, for the model of localized surface diffusion on the graphene sheets, is given by [4]:
D is = D i 0 exp E is k T ,
where the value of Di0 is usually assumed to be constant.
Let us introduce the dimensionless value Ji0= Jit/Ji,Kn, where Ji,Kn is the flux of molecules passing through the effective cross-section of the slit-like pore σe by the classic Knudsen diffusion
J i , Kn = H l D i , Kn d n i d x .
Taking account Equations (14) and (16), we can write the following equation for Ji0
J i 0 = ψ i + γ i ,
where γi = Disviτiαi/(2HDi,Kn). We assume that parameters related to the interaction of gas molecules with the slit-like pore walls (in particular ψi, αi, τi and Dis) are the surface-average values in relation to graphitic areas and oxidized regions of the GO sheets. Equation (17) gives the possibility of analyzing the influence of different factors, related to the interaction of gas molecules with the walls of the slit-like pore, on mass transfer through the pore.
Let us consider the joint influence of surface diffusion and the Maxwell diffuse-specular scheme of the molecule scattering (when ψi = (2 − fi)/fi) or the backscattering scheme for molecules (when ψi = fi/(2 − fi)) on mass transfer in the slit-like pore with the flat parallel walls. Figure 1 shows the dependence of the value of Ji0 on fi for the different values of parameter γi in Equation (17) for the Maxwell diffuse-specular scheme of the gas molecule scattering and the backscattering scheme.
It is seen from Figure 1 that in the case of a diffuse-specular scheme of molecule scattering the value of Ji0 increases with an increase in γi and decreases with an increase in fi and in the case of a backscattering scheme Ji0 increases with an increase in both γi and fi.

4. Validation of Model; Discussing and Comparing the Model With Experimental Data

The very interesting experimental results related to gas transport through GO membrane are given in [12]. These results can be used in the study and optimization of processes such as the separation of gas mixtures and purification of gas phase on the impurity component by means of GO membranes. The authors of the aforementioned paper note that for gas permeation, different transport models were proposed, and these models negate each other and cannot be used to explain the results of the aforementioned paper. The explanation of these experimental results is necessary for understanding the mechanism of mass transfer in GO membranes. Let us consider some results presented in the above-mentioned paper. Figure 1 in this paper shows that the permeance for the three larger molecules CH4, N2 and CO2 corresponds to the data that follow from the classic theory of the Knudsen diffusion. The permeance for the two smaller molecules He and H2 are relatively far above these data. To explain these results authors of the above-mentioned paper, propose a two-pathway transport model for GO membranes: the inter-sheet pathway A and the inner-sheet pathway B. The inter-graphene sheet pathway (pathway A) is composed primarily of inter-galleries between stacked GO sheets which can be modeled by slit-like pores. The inner-sheet pathway (pathway B) represents the structural defects of the GO sheets and it is modeled in this paper by a straight channel. It is also assumed that the tortuosity factor for the pathway A is much greater than for the pathway B (the tortuosity factor for pathway B is assumed to be unity). The obtained experimental results are explained by a difference in the permeance for the larger molecules (CH4, N2 and CO2) and for the smaller molecules (H2 and He) for the pathways A and B. The greater permeance for smaller molecules in GO membrane is explained by their greater permeance by the pathway B. The transport of gas molecules through GO membranes is considered on the basis of the solution-diffusion model. The total permeance for the GO membrane in relation to gas molecules is written in the above-mentioned paper as a sum of the permeances by the pathway A and the pathway B. It should be noted that in a great quantity of graphene sheets of GO membranes, defects of individual sheets can be overlapped by adjacent graphene sheets [24,25] that leads to an increase in the tortuosity factor in diffusion of gas molecules through the membrane by the pathway B and accordingly to a significant decrease in the permeance by the pathway B. As indicated in [25], the average area percentage of holes within a GO sheet is only approximately 2%, which is too small to make a significant impact on the permeability of GO membranes and only the interlayer galleries are assumed to be the major mass-transport pathways. It is worth to note that the model proposed in [12] does not include the parameters characterizing the features of the gas molecule scattering by the walls of pores in GO membranes and surface diffusion of molecules adsorbed on the walls of slit-like pores.
The model involving the parameters describing features of scattering of gas molecules by the walls of the slit-like pores in GO membranes allows explaining some experimental results which are given in [12] without using the concept of two pathways for the transport of gas molecules in GO membranes. It is worth to note that presented in aforementioned paper experimental data show that the presence of molecules of CO2 in the system can lead to a decrease in the permeance of H2 molecules in the GO membrane in comparison with diffusion only H2 molecules (without the presence of CO2 molecules). This effect can be explained by the influence of strongly adsorbed CO2 molecules on the scattering of weakly adsorbed H2 molecules. The fraction of diffusely scattered H2 molecules fH can be estimated as fH = f1(1 − θC) + f2θC = θC(f2f1) + f1 [20], where f1 is the fraction of diffusely scattered H2 molecules on a part of the surface not occupied by adsorbed CO2 molecules, f2 is the fraction of diffusely scattered H2 molecules on a part of the surface occupied by adsorbed molecules of CO2, θC is the surface coverage by molecules of CO2 of the pore walls.
Let us consider for simplicity the case when the parameter γi in Equation (17) is much less than the value of ψi and the effect of surface diffusion on mass transfer through the slit-like pore can be neglected (in particular this can take place for molecules H2 due to small values of the adsorption energy and sticking coefficient for these molecules). In this case, the value of WH = JH/JH,Kn, where JH is the flux of molecules H2 passing through the slit-like pore, JH,Kn is the flux of H2 molecules passing through the pore in accordance with the classic Knudsen equation, is given by
W H = ψ H = 2 f H f H = 2 [ θ C ( f 2 f 1 ) + f 1 ] θ C ( f 2 f 1 ) + f 1 .
Figure 2 shows the dependence of WH on θC at f1 = 0.2 for different values of f2.
It is seen from Figure 2 that the value of JH decreases with an increase in f2 and θC. Figure 3 shows the dependence of WH on f1 for different values of θC at f2 = 1.
It is seen from Figure 3 that the value of WH increases with a decrease in f1. Taking into account [20,26] and aforementioned assumptions, we can write the following equation for the permeance FH for H2 molecules passing through the GO membrane by aforementioned pathway A in the presence of adsorbed molecules of CO2 on the walls of the slit-like pores:
F H = 1 L H v H ε [ 2 θ C ( f 2 f 1 ) f 1 ] 3 R T ξ [ θ C ( f 2 f 1 ) + f 1 ) ] ,
where L is the GO membrane thickness, H′ is the mean value of H for slit-like pores, vH is the mean thermal velocity of molecules of H2, R is the gas constant, while ε and ξ are the porosity and tortuosity for the aforementioned pathway A in the GO membrane.
Let us consider the correlation between the experimental results for permeance of H2 in the GO membrane from reference [12] and the results obtained on the basis of the proposed model. It is worth to note that the aforementioned correlation can be made only qualitatively. This is related to the fact that in the description of experimental results in the above-mentioned paper do not participate parameters characterizing the interaction of gas molecules with the pore walls but in the approach used in the present paper the aforementioned parameters play an essential (in some cases the main) role. It is seen from Figure 2 in the above-mentioned paper that the permeance of weakly adsorbed molecules of pure H2 can be relatively far above the data which correspond to the classic Knudsen diffusion. The similar dependence for the permeance of pure H2 (when θC = 0) follows from Equation (19) of the present paper in f1 < 1. Figure 3 in the aforementioned paper shows the experimental results for the transport of H2 in the presence in the system of strongly adsorbed CO2 molecules. It is seen from this figure that the permeance of H2 in the presence of CO2 is smaller than the permeance of pure H2. This result follows also from Equation (19) in the case θC > 0 and f2 > f1.
It is worth to note that the influence of resonance radiation on the sticking coefficient and the adsorption time of the molecules of CO2 due to excitation of these molecules can lead to a decrease in the values of above-mentioned parameters for molecules of CO2 and accordingly to a decrease in the value of θC. This in turn leads to an increase in the permeance of H2. The aforementioned effect can be used in particular in the separation of H2/CO2. mixtures. Heating of the membrane by radiation also leads to a decrease in θC on the walls of slit-like pores of the GO membrane and accordingly to an increase in the permeance of H2. The aforementioned effects can be used to manage processes such as the membrane separation of gas mixtures and cleaning of gases from impurity molecules. We hope to perform in the near future experiments associated with the verification of the presented theoretical results.
Note that the shape of gas molecules can also affect the transport of gases through GO and NG membranes [10,27]. The influence of the molecule shape on the kinetic diameter of molecules, and accordingly on the transport of gas through GO membrane is discussed in [10]. In [27] it was shown by means of molecular dynamics simulations how the size, shape, and functionalization of pores in NG membranes affect the permeance of CO2, H2S, and CH4 molecules.

5. Conclusions

The transport of gas molecules through NG and GO membranes, taking into account the influence of the surface processes on graphene sheets, was studied theoretically. The combined effect of the features of molecule scattering by graphene sheets and surface diffusion of adsorbed molecules on these sheets on the transport of gas molecules through GO membranes was analyzed. It was shown that some experimental results given in [12] can be explained on the basis of the one-pathway model for mass transfer in GO membranes (diffusion of gas through the slit-like pores of the membranes). In particular, it was shown that the presence in the system of CO2 molecules can decrease the permeance of H2 molecules that corresponds to the experimental results given in the aforementioned paper. It was shown that resonance radiation affecting parameters such as the sticking coefficient and adsorption time of gas molecules can be used for the management of mass transfer through NG and GO membranes. The presented results can be used in the intensification and optimization of processes such as the separation of gas mixtures and natural gas purification using graphene-based membranes. Of course, much more future measurements where we can consider these and possibly other models need to be performed, but that is out of the scope of the present manuscript.

Author Contributions

V.L. is the author of the modeling and contributed in writing of the article; P.I., supervised this work and contributed in writing of the article, K.F. was the author of the idea to do this research and contributed in writing of the article; O.Š. financed this work from her project and contributed in writing of the article. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Czech Science Foundation grant No. 17-05421S, and by EU structural funding in the frame of Operational Programme Research, Development and Education, project No. CZ.02.1.01./0.0/0.0/17_049/0008419 “COOPERATION”.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

DKndiffusion coefficient in classic theory of Knudsen diffusion
Dscoefficient of surface diffusion
Eaadsorption energy
Esactivation energy of surface diffusion
Etactivation energy of molecule translocation
Edenergy barrier of activated Knudsen diffusion
ddiameter of two dimensional pores of porous graphene membrane
dcdiameter of cylindrical pore
dmkinetic diameter of molecule
FHpermeance for H2
fprobability of diffuse scattering of gas molecules
fspprobability of elastic (which is assumed to be specular) reflection
Hwidth of slit-like pore
Jflux density of molecules passing through slit-like channel (pore)
jflux density of molecules passing through porous graphene
kBoltzmann constant
Lthickness of graphene oxide membrane
mmass of gas molecule
mwmass of molecules (atoms) of dense phase (material of pore walls).
nnumber density of gas molecules
nanumber density of adsorbed molecules
npnumber density of the two-dimensional pores on the graphene sheet
pgas pressure
xcoordinate directed along slit-like pore in direction of mass transfer
Qpfrequency of encounters of adsorbed molecules with the two-dimensional pores
vmean thermal velocity of gas molecules
Ttemperature
Greek symbols
αsticking coefficient of molecules
εpporosity of porous graphene membrane
θCsurface coverage by molecules of CO2 of pore walls in graphene oxide membrane
σeeffective cross-section
τadsorption time of molecules
Subscripts
aadsorption
CCO2
ddirect
eeffective
ffeed
HH2
icomponent i
mmolecule
ppermeate
ssurface
spspecular

References

  1. Shi, L.; Zhao, T. Recent advances in inorganic 2D materials and their applications in lithium and sodium batteries. J. Mater. Chem. A 2017, 5, 3735–3758. [Google Scholar] [CrossRef]
  2. Huang, L.; Zhang, M.; Li, C.; Shi, G. Graphene-based membranes for molecular separation. J. Phys. Chem. Lett. 2015, 6, 2806–2815. [Google Scholar] [CrossRef]
  3. Schrier, J. Carbon dioxide separation with a two-dimensional polymer membrane. ACS Appl. Mater. Interfaces 2012, 4, 3745–3752. [Google Scholar] [CrossRef]
  4. Do, D.D. Adsorption Analysis: Equilibria and Kinetics; Imperial College Press: London, UK, 1998. [Google Scholar]
  5. Sun, C.; Bai, B. Fast mass transport across two-dimensional graphene nanopores: Nonlinear pressure-dependent gas permeation flux. Chem. Eng. Sci. 2017, 165, 186–191. [Google Scholar] [CrossRef]
  6. Sun, C.; Bai, B. Diffusion of gas molecules on multilayer graphene surfaces: Dependence on the number of graphene layers. Appl. Therm. Eng. 2017, 116, 724–730. [Google Scholar] [CrossRef]
  7. Drahushuk, L.W.; Strano, M.S. Mechanisms of gas permeation through single layer graphene membranes. Langmuir 2012, 28, 16671–16678. [Google Scholar] [CrossRef] [PubMed]
  8. Tian, Z.; Mahurin, S.M.; Dai, S.; Jiang, D. Ion-gated gas separation through porous grapheme. Nano Lett. 2017, 17, 1802–1807. [Google Scholar] [CrossRef] [PubMed]
  9. Levdansky, V.; Izák, P. Membrane separation of gas mixtures under the influence of resonance radiation. Sep. Purif. Technol. 2017, 173, 93–98. [Google Scholar] [CrossRef]
  10. Bouša, D.; Friess, K.; Pilnáček, K.; Vopička, O.; Lanč, M.; Fónod, K.; Pumera, M.; Sedmidubský, D.; Luxa, J.; Sofer, Z. Thin, high-flux, self-standing, graphene oxide membranes for efficient hydrogen separation from gas mixtures. Chem. Eur. J. 2017, 23, 11416–11422. [Google Scholar] [CrossRef] [PubMed]
  11. Guan, K.; Shen, J.; Liu, G.; Zhao, J.; Zhou, H.; Jin, W. Spray-evaporation assembled graphene oxide membranes for selective hydrogen transport. Sep. Purif. Technol. 2017, 174, 126–135. [Google Scholar] [CrossRef]
  12. Ibrahim, A.; Lin, Y.S. Gas permeation and separation properties of large-sheet stacked graphene oxide membranes. J. Membr. Sci. 2018, 550, 238–245. [Google Scholar] [CrossRef]
  13. Pound, G.M.; Simnad, M.T.; Yang, L. Heterogeneous nucleation of crystals from vapor. J. Chem. Phys. 1954, 22, 1215–1219. [Google Scholar] [CrossRef]
  14. Hirth, J.P.; Pound, G.M. Condensation and Evaporation: Nucleation and Growth Kinetics; Pergamon Press: Oxford, UK, 1963. [Google Scholar]
  15. Lim, Y.; Bhatia, S.K.; Nguyen, T.X.; Nicholson, D. Prediction of carbon dioxide permeability in carbon slit pores. J. Membr. Sci. 2010, 355, 186–199. [Google Scholar] [CrossRef]
  16. Gilron, J.; Soffer, A.J. Knudsen diffusion in microporous carbon membranes with molecular sieving character. J. Membr. Sci. 2002, 209, 339–352. [Google Scholar] [CrossRef]
  17. Holt, J.K.; Park, H.G.; Stadermann, Y.W.M.; Artyukhin, A.B.; Grigoropoulos, C.P.; Noy, A.; Bakajin, O. Fast mass transport through sub–2–nanometer carbon nanotubes. Science 2006, 312, 1034–1037. [Google Scholar] [CrossRef] [PubMed]
  18. Nicholson, D.; Bhatia, S.K. Scattering and tangential momentum accommodation at a 2D adsorbate—Solid interface. J. Membr. Sci. 2006, 275, 244–254. [Google Scholar] [CrossRef]
  19. Karlov, N.V.; Prokhorov, A.M. Selective processes induced by resonance laser radiation at the phase boundary of two media. Sov. Phys. Usp. 1977, 20, 721–735. [Google Scholar] [CrossRef]
  20. Berman, A.S.; Maegley, W.J. Internal rarefied gas flows with backscattering. Phys. Fluids 1972, 15, 772–779. [Google Scholar] [CrossRef]
  21. Levdansky, V.V.; Pavlyukevich, N.V.; Izak, P.; Zdimal, V. Influence of adsorption of foreign gases on mass transfer in fine (nanoscale) capillaries. Heat Transf. Res. 2016, 47, 691–699. [Google Scholar] [CrossRef]
  22. Levdansky, V.V.; Smolik, J.; Moravec, P. Effect of surface diffusion on transfer processes in heterogeneous systems. Int. J. Heat Mass Transf. 2008, 51, 2471–2481. [Google Scholar] [CrossRef]
  23. Levdansky, V.V.; Smolik, J.; Moravec, P. Influence of surface phenomena on free-molecule gas flow in fine channels. Int. Commun. Heat Mass Transf. 2007, 34, 796–800. [Google Scholar] [CrossRef]
  24. Boutilier, M.S.H.; Sun, C.; O’Hern, S.C.; Au, H.; Hadjiconstantinou, N.G.; Karnik, R. Implications of permeation through intrinsic defects in graphene on the design of defect-tolerant membranes for gas separation. ACS Nano 2014, 8, 841–849. [Google Scholar] [CrossRef]
  25. Sun, P.; Wang, K.; Zhu, H. Recent developments in graphene-based membranes: Structure, mass-transport mechanism and potential applications. Adv. Mater. 2016, 28, 2287–2310. [Google Scholar] [CrossRef]
  26. Shirasaka, H.; Shimonosono, T.; Hirata, Y.; Sameshima, S. Analysis of gas permeability of porous alumina powder compacts. J. Asian Ceram. Soc. 2013, 1, 368–373. [Google Scholar] [CrossRef] [Green Version]
  27. Tronci, G.; Raffone, F.; Cicero, G. Theoretical study of nanoporous graphene membranes for natural gas purification. Appl. Sci. 2018, 8, 1547. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Dependence of Ji0 on fi for different values of parameter γi = 1: γi = 5, 2: γi = 3, 3: γi = 1; (a) A diffusion-specular scheme of gas molecule scattering; (b) A backscattering scheme for gas molecules.
Figure 1. Dependence of Ji0 on fi for different values of parameter γi = 1: γi = 5, 2: γi = 3, 3: γi = 1; (a) A diffusion-specular scheme of gas molecule scattering; (b) A backscattering scheme for gas molecules.
Applsci 10 00455 g001
Figure 2. Dependence of WH on θC at f1 = 0.2; 1: f2 = 0.4, 2: f2 = 0.6, 3: f2 = 0.8.
Figure 2. Dependence of WH on θC at f1 = 0.2; 1: f2 = 0.4, 2: f2 = 0.6, 3: f2 = 0.8.
Applsci 10 00455 g002
Figure 3. Dependence of WH on f1 for different values of θC at f2 =1; 1: θC = 0.1, 2: θC = 0.2, 3: θC = 0.3.
Figure 3. Dependence of WH on f1 for different values of θC at f2 =1; 1: θC = 0.1, 2: θC = 0.2, 3: θC = 0.3.
Applsci 10 00455 g003

Share and Cite

MDPI and ACS Style

Levdansky, V.; Šolcová, O.; Friess, K.; Izák, P. Mass Transfer Through Graphene-Based Membranes. Appl. Sci. 2020, 10, 455. https://doi.org/10.3390/app10020455

AMA Style

Levdansky V, Šolcová O, Friess K, Izák P. Mass Transfer Through Graphene-Based Membranes. Applied Sciences. 2020; 10(2):455. https://doi.org/10.3390/app10020455

Chicago/Turabian Style

Levdansky, Valerij, Olina Šolcová, Karel Friess, and Pavel Izák. 2020. "Mass Transfer Through Graphene-Based Membranes" Applied Sciences 10, no. 2: 455. https://doi.org/10.3390/app10020455

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop