Volume 6 - Year 2019 - Pages 9-38
DOI: 10.11159/jffhmt.2019.002

Separation of Butanol Using Pervaporation: A Review of Mass Transfer Models

Hoda Azimi, Jules Thibault, F. Handan Tezel

Department of Chemical and Biological Engineering, University of Ottawa
Ottawa, Ontario, Canada K1N 6N5

Abstract - Pervaporation is a suitable technique for butanol recovery from aqueous solutions especially ABE fermentation broths. The performance of the membrane, in terms of permeation flux and modelling of the mass transfer through membranes provides a deeper understanding, which may assist to orient the research and the development of pervaporation processes. Modelling of the mass transport through the membrane is based on sorption and diffusion of the components into and across the membrane. In this study, an overview of the different models used for the pervaporation separation of butanol is presented. Up to now, the solutiondiffusion based models were the main methods used and the Maxwell-Stefan theory was very limited. In addition, to our knowledge, the pore-flow model was not applied for the modelling of butanol pervaporation separation. Moreover, the Maxwell-Stefan theory seems to be the most accurate mass transfer model to be used for the pervaporation separation method in comparison to the other models.

Keywords: Pervaporation, Butanol, Mass Transfer, Solution-Diffusion, Maxwell-Stefan, Pore-Flow.

© Copyright 2018 Authors - This is an Open Access article published under the Creative Commons Attribution License terms Creative Commons Attribution License terms. Unrestricted use, distribution, and reproduction in any medium are permitted, provided the original work is properly cited.

Date Received: 2017-10-26
Date Accepted: 2018-04-17
Date Published: 2019-02-11

1. Introduction

The depletion of fossil fuels and world concerns about climate change have motivated researchers to find replacement fuels that would be renewable and near carbon neutral such as biodiesel, bioethanol and biobutanol. Many studies have been performed to propose ways to make the production of biofuels, and specially biobutanol, economically viable [1–6]. Biobutanol, a four-carbon alcohol, has enviable properties in comparison to the other biofuels. Some advantages of using biobutanol as a biofuel are its low volatility, low hydroscopicity and lower corrosiveness. The biological production of n-butanol is achieved via Acetone-Butanol-Ethanol (ABE) fermentation. ABE fermentation using Clostridium acetobutylicum yields acetone, butanol and ethanol in a typical ratio of 3:6:1, respectively. However, to make butanol economically viable as a biofuel, the bioconversion efficiency and product purity via the ABE fermentation process should be improved to compensate for the currently very low product yield, product toxicity to microorganisms and multiple end-products [4,7–17].

Many studies have been undertaken in an attempt to increase the efficiency of butanol produced from ABE fermentation [2,3,18–23]. One method to improve the productivity of the ABE fermentation, the in-situ recovery of the solvents, has attracted considerable attention. The main separation methods used for butanol separation from the ABE fermentation process are vacuum fermentation [24] adsorption [25], gas-stripping [26], liquid-liquid extraction [27], perstraction [28,29], reverse osmosis [30] and pervaporation [1]. Among these methods, pervaporation is a highly regarded separation technique because of its low energy consumption, high selectivity and absence of harmful effect on microorganisms [1,31,32]. To mitigate the higher cost associated with the incorporation of a pervaporation membrane for the in situ recovery of ABE solvents from fermentation broths, it is necessary to resort to highly efficient membranes.

The performance of a pervaporation membrane is assessed by the permeation flux and selectivity. To improve membrane performance, a judicious synergistic combination of theory and laboratory work is necessary. To explain the migration of species across pervaporation membranes, a few models have been proposed: solution-diffusion model [33], Maxwell-Stefan theory [34], and pore-flow model [35,36]. With these models, it is possible to infer the performance of pervaporation membranes, but pervaporation remains a complex process to be modelled.

Until now, many reviews have considered the pervaporation fundamentals and membranes; however, there are only few studies published in the literature to review the mathematical models for pervaporation processes [37–45]. In this review paper, a brief summary of the mass transfer models used in pervaporation separation is presented followed by the discussion on the models applied for butanol recovery via this process. The main focus of this manuscript was the models applied to pervaporation of butanol separation.

1. 1. Introduction to Pervaporation

Pervaporation is a partial pressure or more generally concentration driven process, which is the combination of two mass transfer mechanisms: permeation and evaporation [46–48]. During pervaporation, a phase change from liquid to vapour occurs. A liquid feed solution to be separated is in contact with one side of the membrane surface and the permeating product leaves at a low vapor pressure from the other side that is kept under vacuum [49] or swept with a purge gas [50]. The permeate is then condensed or released depending on the objective of the separation [51]. Because of the presence of vacuum or the sweeping of an inert purge gas on the permeate side, a driving force across the membrane prevails. The required heat of vaporization comes from the liquid feed such that the temperature of the feed must be set accordingly. In recent years, several studies have been undertaken to better understand the pervaporation process and use it for numerous separation processes including the separation of water and alcohol mixtures using hydrophobic membranes [41,46,51–54]. Generally, pervaporation applications fall under three categories: (i) removal of water from organic solvents, (ii) removal of organic compounds from aqueous solutions, such as the recovery of the aromatic components and biofuels from fermentation broths, and (iii) separation of anhydrous organic mixtures. In this process, the separation depends on the chemical nature of the macromolecules that comprise the membrane, the physical structure of the membrane, the physicochemical properties of the mixtures to be separated, and the permeant-permeant and permeant-membrane interactions. Figure 1 shows the simplified schematic diagram of a typical pervaporation separation experimental system used to test pervaporation membranes.

Figure 1. Simplified schematic diagram for a typical pervaporation separation setup.

1. 2. Pervaporation Membranes

Membranes, which have been used for pervaporation separation of butanol, are either zeolite membranes such as silicalite zeolite membranes and ultrathin zeolite X films [55–58] or polymeric membranes. Different kinds of polymers have been reported in the literature for butanol separation such as styrene butadiene rubber (SBR) [59], ethylene propylene diene rubber (EPDM) [60], polytetrafluoroethylene (PTFE) [61], polypropylene (PP) [62], polyurethane (polyether based) (PUR) [63], polyether block-amide (PEBA) [64], poly (vinylidinedifluoride) (PVDF) [65], poly (methoxy siloxane) (PMS) [66], poly (dimethylsiloxane) (PDMS) [19], poly (1-(trimethylsilyl)-1-propyne) (PTMSP) [67] and polyamide-imide (PAI) containing cyclodextrin (CD)[68]. Polymeric membranes are less expensive and more flexible in comparison to zeolite membranes [46,69]. However, a composite membrane incorporating these two materials have also been used by some researchers [70–73] to benefit from the outstanding characteristics of each medium to enhance membrane performance.

2. Models Used for Mass Transfer in Pervaporation

Description of mass transport through membranes is essential to enhance the design efficiency and consequently obtaining better separation performance. Up to now, different models have been introduced to study the pervaporation separation process by considering the species behaviour, membrane properties and operating conditions. Similar to the proposed models for membrane-based gas separation applications, a resistance-in-series theory can be used to describe the mass transfer of components through membranes in pervaporation separation system. According to this theory, the sorption of species in the liquid feed takes place at the membrane surface and species then diffuse through the membrane prior to desorbing on the other side of the membrane. However, the mass transfer resistance at the permeate side has been reported to be negligible in the overall mass transfer due to considerably lower pressure in the permeate side [45]. Therefore, desorption step can be modelled in a similar way to the sorption step such that it will not be discussed individually in this study. Figure 2 provides an overview of the different mass transfer models, which have been used in pervaporation separation processes.

In addition, Figure 3 shows the chemical potential (𝜇), total pressure (p) and concentration (C) profiles of a migrating species across the membrane thickness in the pervaporation separation process assuming a high vacuum on the permeate side and linear concentration profile within the membrane.

Figure 2. Summary of different mass transfer models used for pervaporation separation processes (star refers to models not yet used for butanol mixtures).
Figure 3. Gradient profiles across the membrane and the two boundary layers prevailing for the pervaporation system.

2. 1. Solution-Diffusion Model

The solution-diffusion model is used for membranes under the assumption that the mass transport across the membrane proceeds in three steps: (1) the sorption of the penetrants into the membrane, (2) the diffusion of the sorbed components through the membrane as a result of concentration gradient, and (3) the desorption of the diffused components on the permeate side. Moreover, it is assumed that the total pressure is constant within the membrane and the thermodynamic equilibrium has been reached at the two interfaces.

Considering the Fick’s first law of diffusion, the permeation flux of component i can be obtained using Eq. (1) [45].


where Ji is the permeation flux of component i (kmol.m-2.s-1), DM,i is the diffusion coefficient of component i through the membrane (m2.s-1) and CM,i is the concentration of component i within the membrane (kmol.m-3). This equation can be expressed in terms of the bulk concentration or activity of component on both sides of the membrane [43,45] as the mass transfer driving force (Eq. (2)).



In Eq.(2), PM,i is the permeability of component i, which is the product of solubility coefficient (SM,i) and diffusion coefficient (DM,i) [45]. The solution-diffusion model has been widely used for the development and optimization of membranes. However, for a more accurate prediction of membrane performances, it is necessary to determine the diffusion coefficient and the solubility coefficent or their combined product (permeability). The experimental determination of these parameters comes with a significant level of uncertainties, which makes it difficult to use the solution-diffusion model with confidence. Moreover, the permeability by itself does not provide any information that could be used for membrane development. In addition, the original solution-diffusion theory cannot consider the coupling effect between the migrating species which occurs due to the strong interaction between the penetrant molecules and, consequently, this model is only applicable when the coupling flux is negligible in comparison with individual fluxes of penetrants. However, the influence of the coupling flux may be a significant factor which should be taken into account in modelling the mass transfer in pervaporation separation [74]. The interactions between the dissolved components that have high affinity with each other result in the coupling of the species through the membrane penetration. In other words, the presence of one species in a mixture affects the Gibbs free energy of the other components and consequently results in changing the penetration behaviour through the membrane [75–77]. However, the coupling effect consideration in a model would increase the level of complexity of the model and needs additional experimental parameters.

A semi-empirical approach was proposed [45] for the solution-diffusion model where the driving force in the original solution-diffusion model was replaced by a pressure-based fugacity and the equation was derived from Fick’s first law of diffusion for a one-directional mass transfer through the membrane as given in Eq.(3).


where  is the thermodynamic diffusion coefficient (m2.s-1), is the fugacity (bar) and gi is the activity coefficient of component i. The coupling effect of the species was considered in this model by introducing coupling coefficients Bij and Bji to establish a relationship between the local activity coefficients and the local fugacity parameters for a binary feed solution as given in Eq. (4).



Eq.(4) can be further extended to consider ternary feed solutions as discussed by Lipnizki and Tragardh [45].

The integration of Eq. (3) and taking a geometric average (gM,i) between the activity coefficient of the components result in Eq. (5) for the calculation of the permeation flux.



where is the modified thermodynamic diffusion coefficient. According to Franke [78], this model showed a fairly good agreement with the experimental data for a ternary feed mixture of water, ethanol, and butanol using a PAN/PAV pervaporation membrane.

The modified solution-diffusion model would be very useful when coupling effect impacts the mass transfer through the membrane. However, if coupling effects are negligible, this equation would only add unnecessary complexity to the model. Akin to the main solution-diffusion model, semi-empirical models based on the solution-diffusion model can also be used for process and module design.

In the following section, the models which have been used to estimate the sorption peoperties and the diffusivity of the components through membranes are dicussed.

2.1.1. Predictive Models for the Sorption Properties

One of the major mass transfer steps in the pervaporation process is the transfer of a penetrant from the boundary layer to the membrane surface where it can be sorbed. The boundary layer is the thin layer in the vicinity of the membrane interface in which the diffusive mass transfer is the dominating transport mechanism as a result of the decrease of the flow velocity in comparison to the velocity of the bulk liquid [44]. The diffusive flux at the boundary layer of the membrane is a function of the bulk phase velocity, the membrane module geometry, the viscosity of the feed solution and the properties of the penetrants [79]. The diffusion rate of a species in the boundary layer is dictated by the difference of concentration between the bulk liquid and the liquid-membrane interface. The concentration at the interface depends on the solubility of a species in the membrane and the rate of diffusion within the membrane. However, for pervaporation, the rate of diffusion through the membrane in relatively small and the concentration at the surface of the membrane is usually in equilibrium with the liquid bulk concentration. At the interface, there exists a competition for sorption between the different components. The more soluble components with the highest affinity with the membrane material will sorb preferentially into the membrane whereas the other components will sorb less into the polymer and mostly remain in the bulk of the fluid.

Sorption isotherms are used to represent the sorption capacity of a membrane for each species, which is related to the activity or volume fraction of the penetrant in the membrane. Modelling of sorption behaviour is relatively well established and will be briefly reviewed in the following sections. Langmuir and Henry’s Law Isotherms

The most common way to represent sorption capacity is through isotherms, which are most often determined experimentally. Henry’s law isotherm (Eq. (6)) is typically used for rubbery membrane, while Langmuir isotherm (Eq. (7)) is more accurate in the case of glassy polymers [80].



Where CH,i is the sorption capacity (kmol.m-3), kD,i is the Henry’s law constant referring to component i (bar-1), pi is the partial pressure (bar),  is the Langmuir maximum sorption capacity in the polymeric membrane (kmol.m-3) and bi is the adsorption or hole affinity constant (bar-1). Langmuir’s and Henry’s law isotherms are theoretical models to predict the equilibrium sorption of component i inside the membrane. Even though the parameters of these isotherms must usually be obtained experimentally, some attempts have been made with molecular dynamic simulations to predict species solubility [81,82]. However, these isotherms are representative for the prediction of sorption capacity of pure feed in the membrane; therefore, these models are only applicable for dilute solutions where the competition among species sorption is negligible. In addition, more complex models such as the extended Langmuir and dual-mode sorption models might be more appropriate to predict the sorption behaviour of multicomponent mixtures which will be discussed in more details in a subsequent section [45,83,84]. Solubility Parameter Theory

Solubility parameter theory is a semi-empirical model that has generally been used for the selection of a membrane material to separate a specific component from mixtures. This parameter represents the affinity between the migrating species and the membrane by taking into account the interactions between the solute and the polymer [85]. A high affinity between the penetrating component and the membrane prevails when the solubility parameter for them are similar. However, excessive similarity between them would result in the immobilization of the components in the membrane.

The solubility parameter is defined by Eq. (8) where di is the solubility parameter (J1/2.m-3/2), DEvap,i (J.mol-1) is the total energy required to remove a molecule from its adjacent molecules, and Vi is the molar volume of component i (m3.mol-1) [86].


To improve the estimation accuracy of the solubility parameter, Hansen [87] proposed a three-dimensional solubility parameter by considering the total energy (DEvap,i) as a summation of energies required to overcome dispersion forces (DEdf,i), dispersion polar interactions (DEdi,i) and breaking hydrogen bonds (DEhb,i), i.e. . Therefore, the three-dimensional solubility parameter is defined as follows (Eq. (9)):


These parameters are normally used to measure the distance parameter (Δp,i) defined as the distance between two components, polymer and penetrant. A smaller value of D implies a greater affinity between the polymer and the penetrant, which would also increase the membrane swelling during pervaporation.


Zhou et al. studied the affinity of butanol and water for PDMS membranes. Their result revealed that hydrogen bonding and polar interaction parameters are the dominating parameters, which control the affinity of the components. Moreover, Water/PDMS has a greater distance parameter (41.4 MPa1/2) compared to butanol/PDMS (13.0 MPa1/2). The smaller distance for butanol indicates that it has a higher affinity for PDMS compared to water and will be preferentially sorbed into the membrane [88,89].

In another work, the distance parameter for acetone-butanol-ethanol (ABE) solution was reported by Niemisto et al. [90]. A pair of components having a distance value (∆PDMS,i) closer to zero are more likely to have a higher affinity to each other. According to their study, PDMS has the highest affinity towards acetone (PDMS,A=10.6), followed by butanol (PDMS,B=12.4), ethanol (PDMS,E=17.1) and water (PDMS,w=40.9).

To increase the estimation accuracy of the three-dimensional solubility parameter model presented by Hansen, weight factors could be added to Eq. (10) as shown in Eq. (11). The weight factors (W) correspond to the chemical nature of the components [91] and need to be obtained experimentally.


Froehling et al. [92] proposed a modified model to estimate the solubility parameters for a ternary system composed of a binary mixture and the membrane (Eq. (12)). In this equation, fi accounts for the volume fraction of component i in the binary mixture.




As mentioned before, the solubility parameter theory has been mostly used for qualitative and not quantitative purposes. In other words, this theory contributes more to material selection whereas it cannot be employed for process and module design or for the simulation of the mass transfer through polymeric membranes. Flory-Huggins Theory

The Flory-Huggins theory, as a semi-empirical model, has been widely used to estimate the sorption properties of pure species and binary solutions in polymeric materials. Regarding the butanol sorption properties on polymeric membranes, the Flory-Huggins model has been successfully applied for the prediction of the sorption capacity of this component in blended poly (dimethylsiloxane)–benzyl-3-butylimidazolium tetrafluoroborate membranes [90].

According to this theory, the solubility of components in membranes will be a function of Gibbs free energy of interacting substances and could be represented by a set of dimensionless parameters, which are called Flory-Huggins interaction parameters (cij). For a ternary system including a binary mixture of components i and j and the polymeric membrane (non-crosslinked high molecular weight), the Flory-Huggins model is expressed using Eq. (13).


In Eq. (13), is the interaction parameter between components i and j, while and represent the interaction between the components and the polymer. For a binary system consisting of a pure component and a polymeric membrane, Eq. (13) reduces to Eq. (14) by considering that fj is equal to zero.


Moreover, for a binary system (ai = 1 and Vi<<Vp), Eq. (15) could be substituted into Eq. (14) to obtain the Flory-Huggins binary interaction parameter between component i and the polymer (Eq. (16)). Furthermore, the activity coefficient (ai) in the Flory-Huggins equation can be determined using different thermodynamic fluid packages such as NRTL, UNIQUAC, or UNIFAC models.



In Eq. (16),  is assumed constant at a fixed temperature and could be obtained experimentally by performing swelling measurement of the polymer in contact with pure components, by sorption measurement of species in polymer or by inverse gas chromatography. However, since the latter method depends on temperature and concentration, it would be difficult to have an accurate estimation of the Flory-Huggins interaction parameter. Therefore, swelling and sorption experiments are the preferred approaches to estimate the Flory-Huggins interaction parameter.

As far as it was mentioned before, selecting the proper thermodynamic fluid package is required to estimate the activity of species and consequently predicting the sorption properties of the components in polymeric membranes via Flory-Huggins theory. Extended Raoult’s law could be used to predict the activity of components in equilibrium with the liquid feed composition at different operating conditions (Eq. (17)). In this equation, g i and xi are the activity coefficient and molar fraction in the feed, respectively.


The experimental activity coefficient of butanol has been summarized in Table 1 along with the predicted values by different models such as NRTL [93], MOSCED [94], UNIQUAC models [93] for certain specific temperatures. Although the variation between reported experimental activity coefficients is considerable, in general the reported values seem to be different from the ones predicted by theoretical models especially for estimated amounts by the MOSCED method.

2.1.2. Predictive Models for Diffusivity

Diffusion is the second step in pervaporation mass transfer, which is the movement of the molecules based on the chemical potential difference across the membrane. It is also stated that in a highly swollen polymer, the diffusion coefficient of the permeant is related to the degree of swelling, the structure of the polymer, and the permeant size. An increase in the degree of swelling is accompanied by an increase in the free volume inside the polymer which leads to a more permeable pathway for mass transfer of the penetrants and an increase in diffusivity [95–97].

All of the above reports demonstrate that the estimation of the diffusion coefficient is of paramount importance for studying the behaviour of components and the mass transfer through membranes in pervaporation separation processes. Therefore, different models used to calculate the diffusion coefficient are reviewed in the following sections. Empirical Models

Typically, in empirical models, the diffusion coefficient is represented by an exponential function of the solute concentration as shown in Eq. (18).


where , 𝜀i, and  are, respectively, the diffusion coefficient at infinite dilution of the solute (reference value), the softener [98] (empirical constant), and the volume fraction of the solute, for a single component diffusion in a rubbery polymer [45,99].

In a multicomponent separation, the model would be extended since the interactions of the polymer-component and component-component need to be taken into consideration. For instance, in a ternary system, including the membrane, Eq. (19) is applicable [99,100].



In some studies, a linear relationship between the diffusion coefficient and the solute concentration has been suggested [101,102] for glassy polymers [45].

As it is the case with all the empirical models, the complexity of the suggested model would depend highly on the studied system. In other words, as the number of components increases, the empirical parameters that need to be found experimentally increases, which implies time-consuming laboratory work and calculations. Therefore, these types of models can only be used to design the process and membrane module, while it cannot be used for membrane development. Theory of Free Volume

According to the free volume theory, the molecular movement, a physical phenomenon within the polymer matrix, is the basis of diffusion. The model was fundamentally developed from Fick’s first law of diffusion, which is applicable for ideal binary solutions. According to this theory, a penetrant molecule can only be transferred through the polymer if adequate free spaces exist. The free volumes are generated due to the random movement of the polymer segments. In other words, the free volume can be considered as a series of temporary micro-voids that are created in the polymer structure. In addition, in a pervaporation process, the sorption and desorption rates are much higher than that of the diffusion and could be assumed to happen instantly in comparison to the diffusion of species through membrane.

The diffusion coefficient for a binary system (component i and polymer) can be estimated by Eq. (20) [103,104].


where T, Adi and Bdi are the temperature (should be higher than glass transition temperature of the polymer), a parameter which is a measure of shape and size of the penetrant, and a parameter which is a measure of the size of the sorption area in the polymer, respectively. is called the free volume and it depends on the volume fraction of the permeant in the polymer(), the volume fraction of crystalline parts of the polymer (i.e. fC, p) and temperature. In several studies, different equations have been proposed for the determination of  in glassy and semi-crystalline polymers [105,106].

For a ternary system, Eq. (20) can be extended to Eq. (21) [107].


In eq. (21), the parameter  is a generalised type of Bdi to make the equation applicable for both semi-crystalline and glassy polymers, and is evaluated using Eq. (22).


The free volume theory has been reported to give a fairly good approximation of the diffusion coefficient and is applicable for membrane development in different pervaporation systems. It can be used to analize the diffusion through different types of membranes [108]. Dual-Mode Sorption Model

The dual-mode sorption model, similar to the free volume theory, considers the presence of microvoids within the polymer structure as a transportation and accumulation medium. This model is mainly used for glassy polymers even though it could also be applied for rubbery polymers. It is assumed that the microvoids are uniformly distributed throughout the membrane. The solute molecules may be captured by a neighbouring microvoid if the size of the molecule matches the size of the free space. The solute molecule is able to move through the microvoids if it has a sufficient amount of energy. If another solute molecule occupies the place of the transported molecule, diffusion will take place. Two mechanisms which have been considered for mass transfer in this model: (1) the diffusion of molecules through the microvoids (Henry’s sorption sites-CD,i) and (2) the immobilised molecules in the free space (CH,i). The total concentration of species i in the membrane is given by the summation of the two populations as expressed in Eq. (23) [45,83].


The first term of Eq. (23) corresponds to the normal dissolution of the permeant and is represented by the Henry’s law linear isotherm. The second term, corresponding to the immobilised molecules in the microvoids, is nonlinear and expressed by the Langmuir isotherm. In Eq. (23), bi is the microvoid affinity constant and represents the ratio of sorption and desorption rate constants of the penetrant through the free space and is the Langmuir maximum sorption capacity of component i in the polymeric membrane.

The permeability of component i, according to the dual-mode sorption model can be calculated using Eq. (24).




where Fi is the ratio of the mobile phase to the immobilized phase diffusion coefficients, i.e. partial immobilization and Ki is the ratio of the nonlinear to the linear sorption parameters.The extended version of the dual-mode sorption model is applicable for multicomponent systems. For a binary mixture, Eq. (25) can be used [84].



Considering Eq. (24) and Eq. (25), the permeability of component i in a binary mixture can be obtained using Eq. (26) [109].


The dual-mode sorption model has been typically used for gas separation. However, this model could be a suitable model to predict the permeability in a pervaporation process [83]. Morover, this model could be applied for the purpose of membrane development.

In addition to the models which were discussed in this section, dynamic molecular simulation was also performed for the prediction of the diffusion parameters for the mass transfer through a membrane [110,111]. Although there is no study on butanol diffusion through organic membranes using molecular simulation, it has been mentioned that the molecular simulation would be the model that could be used in the future for analizing the diffusion behaviour [45].

2.1.3. Solution-Diffusion Model for Butanol Pervaporation Separation

Instead of considering the sorption, diffusion and desorption steps separately, the overall mass transfer models take all three steps globally to calculate the permeation flux. The models do not offer a comprehensive understanding of the underlying phenomena but are nevertheless very useful for pervaporation process design and cost estimation. Different empirical models were proposed in the literature and consider the entire mass transfer mechanisms as a whole without considering directly the physico-chemical interactions. These models relies on experimental results to estimate the model parameters and the resulting models are then used for predicting the permeation flux and selectivity. The predictions are obviously valid to specific cases that are similar to the experimental system that was used to generate the data for fitting the model. Therefore, they have very limted use for membrane development.

El-Zanati et al. used the resistance-in-series model based on the solution-diffusion concept to validate the pervaporation process of a binary solution of butanol and water via Pervap 2200 (PVA crosslinked/PAN support) membranes [112]. In this work, to predict the permeation flux, the overall mass transfer resistance was considered including three resistance steps: bulk feed convection, diffusion through the membrane, and the convective removal of permeating species on the permeate side. The flux associated with these three steps can be calculated using Eq. (27).




where kL is the mass transfer coefficient at the liquid side, lM is the thickness of the membrane (m) and kv is the mass transfer coefficient at the vapour side. Under steady state conditions, the molar flux is constant at all positions within the membrane, which can be represented with a simple equation in terms of an overall mass transfer coefficient or overall mass transfer resistance (1/Kov) as presented in Eq. (28).



In the case that the partial pressure on the permeate side approaches zero, the flux in Eq. (28) can be further simplified to Eq. (29).


El-Zanati et al. used Eq. (29) to determine the change of butanol concentration as a function of time during the pervaporation process [112]. Performing a mass balance on the feed tank and using Eq. (29), Eq. (30) was obtained to estimate the butanol concentration in the feed tank as a function of time.


where Jtot is the total flux (kmol.m-2.s-1), Am is the membrane area (m2), V0T is the volume of feed tank (m3), and Cb0 is the initial butanol concentration. Based on the results presented by El-Zanati et al., the prediction model provided a good approximation of the flux and butanol concentration as a function of pervaporation time and represented the experimental results very well.

In a different work, Plaza et al. [33] used experimental sweep gas pervaporation data for the separation of butanol from ABE mixtures using a supported PTFE-gelled ionic liquid membrane and a resistance-in-series model to predict the mass transfer and diffusion coefficients. A semi-empirical procedure was used to estimate the diffusion coefficient, where the theoretical flux was assumed equal to the experimental flux and the diffusion coefficient was obtained by minimizing the mean-squared differences between the calculated and experimental values. Results showed a good fit with experimental data. However, as it was mentioned in their paper, a large difference between the calculated diffusion coefficients and values reported in similar works was observed. The values of the diffusion coefficient for butanol in ionic liquid (IL) membranes calculated in their work were approximately twice the values reported by Vopicka et al. [33,113,114].

Li et al. [18] used the solution-diffusion model based on the Fick’s first law of diffusion for butanol pervaporation separation from a binary solution of water/butanol, ABE model solutions containing different concentrations of glucose and other main metabolites, and ABE fermentation broths. The authors used a PDMS membrane and were able to quantify the coupling effect through the calculation of the overall mass transfer coefficient for each case. The overall mass transfer coefficient was estimated via the slope of the flux versus the concentration plot. Results showed that the estimated overall mass transfer coefficient for the binary solutions was larger than the estimated coefficient of the ABE model and real fermentation broths by approximately 16 and 28%, respectively. In addition, the coupling effects in ABE model solutions and fermentation broths, the complex interactions among the metabolites, media broth, and bio-macromolecules (such as blocking sorption sites and free volume) could also partly explain the difference between the calculated mass transfer coefficients. As the estimation of the overall mass transfer coefficients was strictly based on experimental data, it is not possible to generalize this information to other pervaporation separation systems. However, the insight regarding the effect of coupling and other complex interactions is useful.

Non-equilibrium solution-diffusion was applied by Yang et al. [115] for the pervaporation separation of organic compounds (acetone-butanol-ethanol) from binary aqueous solutions using PDMS membranes. In order to estimate the permeation flux, The concentration polarization and coupling effects have been considered negligible.

The diffusion coefficient of the components and the sorption rate were estimated semi-empirically based on the data from the reciprocal of the flux versus the change in membrane thickness. The model demonstrated a good agreement with the experimental data at low feed concentrations. However, some deviations were observed at high concentrations of the organic compounds in the feed due to the intensified swelling of organics with PDMS polymer chain.

Li et al. [116] used a resistance-in-series model to study the effects of concentration polarization on the pervaporation separation of butanol from an aqueous binary solution by a tri-layer PDMS membrane, where a layer of hydrophobic PE was placed in between a PDMS layer and the metal support. In their work, the overall mass transfer coefficient was estimated using solution-diffusion model. However, the mass transfer coefficient of butanol in the feed channel was calculated semi-empirically using a Sherwood number correlation for rectangular membrane module (Eq. (31)).


where dH and lm are the characteristic length of the rectangular membrane channel (hydraulic diameter) and the membrane module length, respectively. For this semi-empirical model, the mass transfer coefficient of butanol in the bulk layer is more than three times higher than that of the overall mass transfer coefficient. This clearly indicates that mass transfer is controlled by membrane diffusion, and the concentration polarization has a negligible effect on the transport process even when the system is operating in the laminar flow regime at a Reynolds number of 140.

In another work, Valentínyi et al. [117] improved the basic solution-diffusion model of Rautenbach et al. by considering that the transport coefficient as an exponential function of the permeating compound. In the basic solution-diffusion model, it is assumed that the transport coefficient has negligible concentration dependency and can be considered constant [118]. However, a large deviation between the experimental data and the model predictions was found at higher concentrations. The transport coefficient or permeance of component i (), defined in Eq. (32), was used by Rautenbach et al. because the concentration dependence of the diffusion coefficient was considered to be negligible.


In the work of Valentínyi et al., the authors assumed that the transport coefficient varies exponentially with concentration as shown in Eq. (33) since the model suggested by Rautenbach et al. could not be used at high concentrations of the permeating component.


where B* is a constant. Results obtained by Valentínyi et al. for the dehydration of butanol aqueous solutions by pervaporation using commercial hydrophilic polyvinyl acetate membranes (PVA)/PAN) showed that the estimation of the transport coefficient using Eq. (33) provides better predictions of butanol and water fluxes than the estimations using Eq. (32). The improved solution-diffusion model suggested by Valentínyi et al. could be used for azeotropic solutions with high water content as well as for modelling, design, and optimization of pervaporation processes.

Ashraf et al. used the model proposed by Valentínyi et al. to calculate the size and arrangement of the pervaporation modules [119]. The experimental results for the dehydration of an aqueous solution of butanol using a commercial pervaporation membrane were in good agreement with the proposed model. However, the model is limited to the type of membrane, components and operating conditions used in their experiments [119].

Petrychkovych et al. used the solution-diffusion model for calculating the pervaporation of butanol flux for a binary butanol/water mixture using a PE membrane. The calculated permeation flux was obtained assuming a constant diffusion coefficient, with a good agreement with the experimental permeation flux [120].

To gain a better understanding of the permeation of ABE species through polymeric membranes, it is paramount to investigate the solubility coefficient and diffusion coefficient of these components in various membranes. However, there are not many values reported in the literature for the solubility of ABE components in membranes. Table 2 presents the diffusion coefficients of ABE components reported in several investigations. It is shown that the diffusion coefficients for butanol in PDMS membranes are of the same order of magnitude in the different studies except for the values reported by Yang et al. [115] which are significantly smaller. The same trend was observed for the diffusion coefficients of acetone. For ethanol, the variation of the diffusion coefficients in PDMS membranes was more pronounced.

2.2. Maxwell-Stefan Theory

The generalized Maxwell-Stefan equation is based on the assumption that the migration of species is the result of a driving force that is counteracted by the friction of the migrating species from the other species and the immediate environment. It was formalized by Mason and Viehland for multicomponent systems and later implemented for membrane separation processes such as gas separation and pervaporation [121,121,122]. According to this theory, for multicomponent systems, the sum of the friction forces for a component, balances the driving force as expressed in Eq. (34) [45]. In this equation, represents the inverse draging force between species i and j, and vi and vj are the diffusion velocities for components i and j.


The chemical potential in Eq. (34) is related to the activity of components within the membrane assuming to have an ideal gas mixture. Eq. (35) represents a simplified form of Maxwell-Stefan equation for a ternary system including binary solution and the membrane material [45,123].




The diffusion coefficients in Eq. (37), i.e. and  are the effective diffusion coefficients for component i and j.

Eq. (35) was rearranged to calculate the permeation flux of each component through the membrane (this model needs to be solved numerically) and it was assumed that the effective diffusion coefficients are equal () due to symmetrical conditions, which led to Eq. (36) [45].



When coupling effects are negligible, (and consequently ) approaches to infinity and Eq. (36) would be simplified to the Fick’s first law of diffusion, where  and  are the Fick’s diffusion coefficients as shown in Eq. (37) [45,123].



Other researchers have used the same type of equations where the volume fraction was simply replaced by the weight or mole fraction in Eq. (36) [124,125]. Moreover, some modifications have been applied in different studies to obtain analytical solutions from the differential equations (Eq. (36)) using a few assumptions to simplify the integration procedure. The main advantage of Maxwell-Stefan theory could be its ability to predict the flux and selectivity for non-ideal multicomponent systems based on the results of single components, which significantly decreases the number of pervaporation experiments to be performed. Moreover, it has the potential to be used for module and process design together with membrane development.

Bettens et al. [34] investigated the deviations observed between the experimental results and a theoretical model that combined the extended Langmuir model for sorption and the Maxwell-Stefan diffusion for pervaporation of butanol/water and butanol/methanol. However, in their study, the desired components to be separated were water and methanol. Different types of methylated microporous silica membranes were used for this purpose. In addition, to predict the Maxwell–Stefan diffusion coefficients, three Maxwell–Stefan diffusivities were determined:,  and . Moreover, for the counter-exchange Maxwell–Stefan diffusivity, the Vignes equation (Eq. 40) was used.


The single-species Maxwell–Stefan diffusivities were used for the conditions where there were no interactions between different species and the corresponding counter-exchange coefficient was infinite (). In this case, species were not able to pass each other in the narrow pores. The single-species Maxwell–Stefan diffusivities could either be independent of the amount adsorbed (Eqs. 39a and 39b) or dependent on the amount adsorbed (Eqs. 39c and 39d).





The comparison was made by assuming different case scenarios where the diffusion coefficients between the components and the membrane in the Maxwell-Stefan model were estimated. The pure alcohol Maxwell–Stefan diffusivity at zero coverage  was derived from the experimental pure alcohol flux and the pure water Maxwell–Stefan micropore diffusivity at zero coverage  was achieved by fitting the experimental water flux data. The four scenarios that were proposed are as follows:

1) There are no interactions between feed species (). Moreover,  and  are independent from the amount adsorbed  (Eqs.39a and 39b);

2) There are no interactions between feed species (). Moreover  and  are dependent on the amount adsorbed (Eqs.39c and 39d);

3) There are interactions between feed species and  can be obtained by Vignes equation (Eq. 40).


Moreover,  and  are independent from the amount adsorbed (Eqs.39a and 39b);

4) There are interactions between feed species and  can be obtained by Eq.40. Moreover,  and  are dependent on the amount adsorbed (Eqs.39c and 39d).

According to the results, Maxwell-Stefan theory was consistent with experimental pervaporation data for the first case scenario when no interactions took place between the components, and the diffusion coefficients of the components through the membrane were independent from the amount adsorbed. In addition, prediction performed with the second scenario for the calculation of the flux of components demonstrated to be closer to the experimental data compared to the first scenario. However, the last two case scenarios did not have a good fit with experimental data.

In another work, the mass transfer through a PDMS commercial membrane for pervaporation separation of the butanol from aqueous solution has been studied by Ebneyamini et al. [126]. A semi-empirical approach on the Maxwell-Stefan model was extended to consider the effect of membrane swelling and the operating temperature on the diffusion coefficient and sorption properties. In order to estimate the effect of the operating temperature on the diffusion coefficients and the sorption of each species, Arrhenius-type equations were incorporated into the Maxwell-Stefan model. Moreover, similar to the free volume theory, an exponential relationship was used to express the diffusion coefficient of each component as a function of the degree of swelling of the membrane at a constant temperature (Eq. 41).



The result of the extended Maxwell-Stefan model showed a better agreement with the experimental data in comparison to the Maxwell-Stefan model. Moreover, the model was able to predict the membrane properties such as solubility and diffusivity at different operating temperatures and feed concentrations.

2.3. Pore-Flow Model

The other proposed model for the overall mass transfer in a pervaporation separation process is the pore-flow model. The difference between this model and the solution-diffusion model is that the solution-diffusion model considers no phase change in the membrane while, in the pore flow model, there is a phase change inside the membrane at a certain distance from the membrane surface. Considering the phase change inside the membrane, there should be a phase boundary between the liquid and the vapor for the mass transfer equations. In 1991, Okada and Matsuura [34] considered a pore flow model for describing the mass transfer through a cellulose membrane in pervaporation separation of ethyl alcohol/heptane mixtures. In this model, it is assumed that the selective layer of the membrane is formed by a series of straight cylindrical pores distributed within the membrane surface and isothermal conditions are assumed. Furthermore, the length of the pores is equal to the thickness of the active layer.

Three steps have been proposed for this model: (1) liquid transport from the pore inlet to the liquid-vapor phase boundary; (2) evaporation at the phase boundary; (3) vapor transport from the phase boundary to the pore outlet [35,36]. Moreover, it was assumed that at the liquid-vapor phase boundary, the components reach to their saturated pressure in equilibrium with the feed solution.

According to the above-mentioned assumptions, at steady state for single component permeation, the flux in the liquid and vapor segments can be calculated according to Eq. (42).




where subscripts L and V represent the liquid and vapor phases, while l is the thickness of each phase inside the pore. A and B constants are obtained using the Darcy’s equation and Henry’s law (Eq. (43)).



where lad is the thickness of the adsorption monolayer.

Since the constants A and B are determined empirically, the pore-flow model would be limited to assist only in process and module design and it is not applicable for membrane development.

Moreover, there are additional models such as pseudo phase-change solution-diffusion, which combines the specifications of both the solution-diffusion and the pore-flow models in one single framework. The pseudo phase-change solution-diffusion takes the coupling effect from the solution-diffusion model and the pseudo phase-change inside the membrane from the pore-flow model into consideration [127,128]. The pore-flow model has not been used yet for the separation of butanol by pervaporation.

3. Conclusion

Pervaporation is an appropriate process for the separation of butanol from different aqueous solutions especially ABE fermentation broths. However, to be employed at an industrial level using lab-scale experimental data, mathematical models are paramount to predict pervaporation process performances accurately. Such models can assist in the design phase of membrane modules for the purpose of optimization, flow patterns prediction, and vessel design. In addition, the membrane development using modelling could provide valuable information before membrane fabrication.

In this study, mass transfer modelling in a pervaporation separation of butanol from ABE fermentation broth, ABE model solutions and aqueous solutions were reviewed in this study.

Considering the previous works for the pervaporation of butanol, although the number of the diffusion models is very limited, the overall mass transfer was mainly used by researchers to estimate the mass transfer coefficient of the components using sorption models. Among the sorption models, which have been used for pervaporation separation process, Langmuir and Henry’s law isotherms are deemed accurate models for the sorption prediction of butanol into membranes. In addition, in comparison to the theory of free volume, the dual-mode sorption model was found to be a simpler and more applicable model for estimating the butanol diffusion coefficient. As shown in Table 3, the solution-diffusion model has been the most frequently used model for the description of the mass transfer in pervaporation separation of butanol. However, the application of the Maxwell-Stefan theory for this purpose has been very limited and, to our knowledge, the pore-flow model has not been reported in the literature for the pervaporation of butanol.

Most of the models used for the pervaporation separation of butanol are semi-empirical models, which fall in between theoretical and empirical models in terms of complexity. Although, in the semi-empirical models, the driving force is well established, the permeability would be experimentally estimated. The complexity of the model increases with the number of components involved in the mass transfer and it further increases by considering the coupling effects. In this case, the Maxwell-Stefan theory seems to be an appropriate option for considering the coupling phenomenon. In addition, it is an accurate model for membrane separation in pervaporation due to its ability to predict the flux and selectivity of the multi-component systems based on the results of single components, which significantly decreases the number of pervaporation experiments to be performed. In addition, it has the potential to be used for module and process design together with membrane development.

Furthermore, while the information obtained via the overall mass transfer models would be used for fundamental understanding of the process and module design, other information such as concentration polarization, pressure drop, flow pattern, and heat transfer also need to be considered for the purpose of process design. Moreover, molecular simulation could become in the future a very valuable tool for membrane development in the pervaporation of butanol. This study can be helpful in combining different sorption and diffusion models to achieve more rigorous models for the prediction of butanol separation by pervaporation.


ABE Acetone, Butanol, Ethanol
CIRC Circulation equilibrium still
COSMO-RS Conductor like screening model for real solvents
DDST Differential distillation
EBUL Ebuliometry
EPDM Ethylene propylene diene rubber
GEM-RS Generalized regular solution model
GLC Gas-liquid chromatography
HSA Headspace analysis
IGS Inert gas stripping
LLE Liquid-liquid extraction
LSG Local surface Guggenheim equation
MBEA Molecular beams
MMM Mixed matrix membrane
MOSCED Modified separation of cohesive energy density model
NRTL Non-random two-liquid model
NSGLC Non-steady state gas-liquid chromatography
PAI Polyamide-imide
PAN Polyacrylonitrile
PAV Polyarylene vinylene
PDMS Polydimethylsiloxane
PE Polyethylene
PET Polyethylene terephthalate
PEBA Polyether block-amide
PRV Phase ratio variation method
PV Pervaporation
PVA Poly(vinyl acetate)
PI Polyimide
POMS Poly(octhylmethyl siloxane)
PMS Poly (methoxy siloxane)
PP Polypropylene
PTFE Polytetrafluoroethylene
PTMSP Poly (1-(trimethylsilyl)-1-propyne)
RDIST Rayleigh distillation
RO Reverses osmosis
TENS Tensimetry
UNIQUAC Universal quasichemical
VPC Vapour phase calibration
WWC Wetted wall column
a activity [-]
Adi free volume parameter of component i for the glassy region in polymer [-]
Ai dimensionless constant in Eq. (40)
AM membrane area [m2]
Apore constant defined for a pure component system in the pore flow model [-]
b adsorption or hole affinity constant [bar-1] coupling coefficients in semi-empirical
B model after Meyer-Blumenroth [kmol.m-1.s-1.bar2]
B* Constant parameter in Eq. (35) [-]
Bpore constant for a pure component system in the pore flow model [-]
Bdi free volume parameter of component i for the glassy region in polymer [-]
generalised free volume parameter [-]
C concentration [kmol.m-3, kg.m-3]
CD concentration of diffusing in membrane [kmol.m-3, kg.m-3]
CH concentration in microvoids [kmol.m-3, kg.m-3]
Langmuir maximum sorption capacity in the polymeric membrane [kmol.m-3, kg.m-3]
D diffusion coefficient [m2.s-1]
D0 diffusion coefficient in infinite dilution [m2.s-1]
DD diffusion coefficient in the Henry’s Law mode [m2.s-1]
DH diffusion coefficient in the Langmuir mode [m2.s-1]
thermodynamic diffusion coefficient [m2.s-1]
modified thermodynamic diffusion coefficient [m2.s-1]
effective concentration-dependent diffusion coefficient [m2.s-1]
Maxwell-Stefan interaction parameter [m2.s-1]
permeance [kg.m-2.s-1]
ΔEdi, energy required to overcome dispersion polar interactions [J.mol-1]
ΔEhb energy required to overcome hydrogen [J.mol-1
ΔEvap energy of vaporisation [J.mol-1]
f fugacity [bar]
f FV free volume [-]
F as defined in Eq. (24) [-]
ΔH heat of adsorption [J.mol-1]
J flux [kmol.m-2.s-1]
Henry’s law constant referring to component i [bar-1]
dimensionless Henry’s law constant [g.m-3/g.m-3]
k boundary layer mass transfer coefficient [-]
(unit weight of polymer per volume of adsorbed gas molecule i) x kD,i[mol.m-3.bar-1]
Ki as defined in Eq. (24)[-]
Kov overall mass transfer coefficient across the membrane[-]
lad thickness of the adsorption monolayer [m]
lM thickness of membrane M [m]
lL length of the liquid-filled proportion of the pore in the pore flow model [m]
lV length of the vapour-filled proportion of the pore in the pore flow model [m]
M molar weight [kg.kmol-1]
Nt total number of pores per effective membrane area [-]
p partial pressure [bar, Pa]
p* Pressure in boundary layer in Eq. (27) [bar]
P permeability with reference to an activity driving force [kmol.m-1.s-1.bar-1]
rpore pore radius [m]
R gas constant [J.mol-1.K-1]
Re Reynolds number [-]
S Solubility coefficient  [kmol.m-3]
Sc Schmidt number [-]
Sh Sherwood number [-]
t time [s]
T temperature [K]
v diffusion velocity [m2.s]
V molar volume [m3.mol-1]
V0T the volume in feeding tank [m3]
W1, W2, W3 weight factors for three dimensional solubility parameters [-]
x mole fraction [-]
xi,b the bulk mole fraction [-]
x* the interface mole fraction
z z-co-ordinate
β proportional constant in free volume theory [-]
activity coefficient [-]
average activity coefficient of component i in membrane  [-]
solubility parameter [J1/2.m-3/2]
solubility parameter due to hydrogen bonds with reference to component i [J1/2.m-3/2]
solubility parameter due to dispersion forces with reference to component i [J1/2.m-3/2]
solubility parameter due to dispersion polar interactions with reference to component i [J1/2.m-3/2
𝜀 empirical constant or ‘softener’ [-]
η liquid viscosity [Pa.s]
𝜇 chemical potential [J.mol-1]
𝜇* Surface viscosity of adsorptive layer of vapour [Pa.s]
ρ density [kg.m-3]
fc crystallinity of polymer P [-]
f volume fraction [-]
as defined in Equation (40), (degree of swelling) [g.g-1]
Flory-Huggins binary interaction parameter [-]
Δp,i distance between polymer P and component i  in _-space [J1/2.m-3/2]
b butanol concentration
b0 initial butanol concentration
exp experimental
f feed
i,j component i and j
L liquid phase
m module
M membrane
p polymer
P permeate
sat saturated
tot Total
V vapour phase
0 reference
FV free volume
L liquid phase
mix Mixture


[1]      N. Abdehagh, J. Thibault, F. H. Tezel, "Separation techniques in butanol production: Challenges and developments," Biomass Bioenergy, vol. 60, pp. 222-246, 2014. View Article

[2]      Y. Zhu, F. Xin, Y. Chang, Y. Zhao, W. Weichong, "Feasibility of reed for biobutanol production hydrolyzed by crude cellulase," Biomass Bioenergy, vol. 76, pp. 24-30, 2015. View Article

[3]      V. Ujor, A. K. Bharathidasan, K. Cornish, T. C. Ezeji, "Feasibility of producing butanol from industrial starchy food wastes," Applied Energy, vol. 136 pp. 590–598, 2014. View Article

[4]      V. Erick, A. Tom, D. Sofie, Industrial biomass: Source of chemicals, materials, and energy!, BACAS, 0455-02, 2011.

[5]      M. G. Buonomenna, J. Bae, "Membrane processes and renewable energies," Renewable and Sustainable Energy Reviews, vol. 43, pp. 1343-1398, 2015. View Article

[6]      I. L. Borisov, G. S. Golubev, V. P. Vasilevsky, A. V. Volkov, V. V. Volkov, "Novel hybrid process for bio-butanol recovery: Thermopervaporation with porous condenser assisted by phase separation," Journal of Membrane Science, vol. 523, 291-300, 2017. View Article

[7]      J. Du, A. McGraw, J. A. Hestekin, "Modeling of Clostridium tyrobutyricum for Butyric Acid Selectivity in Continuous Fermentation," Energies, vol. 7, pp. 2421-2435, 2014. View Article

[8]      C. Park, C. Geng, "Mathematical modeling of fed-batch butanol fermentation with simultaneous pervaporation," Korean Journal of Chemical Engineering, vol. 13, pp. 612-619, 1996. View Article

[9]      Z. Li, Z. Shi, X. Li, "Models construction for acetone-butanol-ethanol fermentations with acetate/butyrate consecutively feeding by graph theory," Bioresource Technology, vol.  159, pp. 320-326, 2014. View Article

[10]    S. B. Bankar, S. A. Survase, H. Ojamo, T. Granstrom, "Biobutanol: the outlook of an academic and industrialist," RSC Advances, vol. 3, pp. 24734-24757, 2013. View Article

[11]    H. Dong, W. Tao, Z. Dai, L. Yang, F. Gong, Y. Zhang, Y. Li, "Biobutanol," Advances in Biochemical Engineering Biotechnology, vol. 128, pp. 85-100, 2012. View Article

[12]    G. Festel, “Competence for Second Generation Biofuels Based on Lignocellulose,” BUTALCO, GmbH company, 2014.

[13]    C. L. Gabriel, Butanol fermentation process, Industrial and Engineering Chemistry, vol. 20, pp. 1063-1067, 1928. View Article

[14]    V. Garcia, J. Pakkila, H. Ojamo, E. Muurinen, R.L. Keiski, "Challenges in biobutanol production: How to improve the efficiency?," Renewable and Sustainable Energy Reviews, vol. 15, pp. 964-980, 2011. View Article

[15]    D. Ramey, Butanol: The Other Alternative Fuel, Agric. Biofuels Technol. Sustain. Profitab., 2007 137-147.

[16]    T. Ezeji, N. Qureshi, H. Blaschek, "Butanol fermentation research: Upstream and downstream manipulations," The Chemical Record, vol. 4, no. 5, pp. 305-314, 2004. View Article

[17]    B. Bharathiraja, J. Jayamuthunagai, T. Sudharsanaa, A. Bharghavi, R. Praveenkumar, M. Chakravarthy, D. Yuvaraj, "Biobutanol - An impending biofuel for future: A review on upstream and downstream processing tecniques," Renewable and Sustainable Energy Reviews, vol. 68, pp. 788-807, 2017. View Article

[18]    S.-Y. Li, R. Srivastava, R. S. Parnas, "Study of in situ 1-Butanol Pervaporation from A-B-E Fermentation Using a PDMS Composite Membrane: Validity of Solution-Diffusion Model for Pervaporative A-B-E Fermentation," Biotechnology Progress, vol. 27, pp. 111-120, 2011. View Article

[19]    G. Liu, D. Hou, W. Wei, F. Xiangli, W. Jin, "Pervaporation Separation of Butanol-Water Mixtures Using Polydimethylsiloxane/Ceramic Composite Membrane," Chinese Journal of Chemical Engineering, vol. 19, no. 1, pp. 40-44, 2011. View Article

[20]    N. Qureshi, M. M. Meagher, J. Huang, R. W. Hutkins, "Acetone butanol ethanol (ABE) recovery by pervaporation using silicalite-silicone composite membrane from fed-batch reactor of Clostridium acetobutylicum," Journal of Membrane Science, vol. 187, pp. 93-102, 2001. View Article

[21]    N. Qureshi, M. M. Meagher, R. W. Hutkins, "Recovery of butanol from model solutions and fermentation broth using a silicalite silicone membrane," Journal of Membrane Science, vol. 158, pp. 115-125, 1999. View Article

[22]    N. Abdehagh, F. H. Tezel, J. Thibault, "Adsorbent screening for biobutanol separation by adsorption: kinetics, isotherms and competitive effect of other compounds," Adsorption, vol. 19, pp. 1263-1272, 2013. View Article

[23]    Q. Wang, Y. ding, L. Liu, J. Shi, J. Sun, Y. Xue, "Engineering Escherichia coli for autoinducible production of n-butanol," Electronic Journal of Biotechnology, vol. 18, pp. 138-142, 2015. View Article

[24]    A. P. Mariano, N. Qureshi, R. Maciel Filho, T. C. Ezeji, "Assessment of in situ butanol recovery by vacuum during acetone butanol ethanol (ABE) fermentation," Journal of Chemical Technology & Biotechnology, vol. 87, no. 3, pp. 334-340, 2012. View Article

[25]    A. Oudshoorn, L. A. M. van der Wielen, A. J. J. Straathof, "Adsorption equilibria of bio-based butanol solutions using zeolite," Biochemical Engineering Journal, vol. 48, pp. 99-103, 2009. View Article

[26]    N. Qureshi, H. P. Blaschek, "Recovery of butanol from fermentation broth by gas stripping," Renewable Energy, vol. 22, no. 4, pp. 557-564, 2001. View Article

[27]    W. J. Groot, H. S. Soedjak, P. B. Donck, R. G. J. M. van der Lans, K. C. A. M. Luyben, J. M. K. Timmer, "Butanol recovery from fermentations by liquid-liquid extraction and membrane solvent extraction," Bioprocess Engineering, vol. 5, no. 5, pp. 203-216, 1990. View Article

[28]    N. Qureshi, I. S. Maddox, "Reduction in Butanol Inhibition by Perstraction," Food and Bioproducts Processing, vol. 83, no. 1, pp. 43-52, 2005. View Article

[29]    G. Merlet, F. Uribe, C. Aravena, M. Rodríguez, R. Cabezas, E. Quijada-Maldonado, J. Romero, "Separation of fermentation products from ABE mixtures by perstraction using hydrophobic ionic liquids as extractants," Journal of Membrane Science, vol. 537, pp. 337-343, 2017. View Article

[30]    A. Garcia, E. L. Iannotti, J. L. Fischer, "Butanol fermentation liquor production and separation by reverse osmosis," Biotechnology and Bioengineering, vol. 28, pp. 785-791, 1986. View Article

[31]    N. Qureshi, H. P. Blaschek, "Butanol recovery from model solution/fermentation broth by pervaporation: evaluation of membrane performance," Biomass and Bioenergy, vol. 17, no. 2, pp. 175-184, 1999. View Article

[32]    C. Shin, X. Chelsea Chen, J. M. Prausnitz, N. P. Balsara, "Effect of block copolymer morphology controlled by casting-solvent quality on pervaporation of butanol/water mixtures," Journal of Membrane Science, vol. 523, 2017, pp. 588-595. View Article

[33]    A. Plaza, G. Merlet, A. Hasanoglu, M. Isaacs, J. Sanchez, J. Romero, "Separation of butanol from ABE mixtures by sweep gas pervaporation using a supported gelled ionic liquid membrane: Analysis of transport phenomena and selectivity," Journal of Membrane Science, vol. 444, pp. 201-212, 2013. View Article

[34]    B. Bettens, A. Verhoef, H.M. van Veen, C. Vandecasteele, J. Degrève, B. Van der Bruggen, "Pervaporation of binary water–alcohol and methanol–alcohol mixtures through microporous methylated silica membranes: Maxwell–Stefan modeling," Computers & Chemical Engineering, vol.  34, no. 11, pp. 1775-1788, 2010. View Article

[35]    T. Okada, T. Matsuura, "A New Transport Model for Pervaporation," Journal of Membrane Science, vol. 59, no. 2, pp. 133-150, 1991. View Article

[36]    T. Okada, M. Yoshikawa, T. Matsuura, "A Study on the Pervaporation of Ethanol Water Mixtures on the Basis of Pore Flow Model," Journal of Membrane Science, vol. 59, no. 2, pp. 151-168, 1991. View Article

[37]    H. Vinh-Thang, S. Kaliaguine, "Predictive Models for Mixed-Matrix Membrane Performance: A Review," Chemical Reviews, vol. 113, no. 7, 2013, pp. 4980-5028. View Article

[38]    S. Semenova, H. Ohya, K. Soontarapa, "Hydrophilic membranes for pervaporation: An analytical review," Desalination, vol. 110, pp. 251-286, 1997. View Article

[39]    P. Shao, R. Y. M. Huang, "Polymeric membrane pervaporation," Journal of Membrane Science, vol. 287, pp. 162-179, 2007. View Article

[40]    B. Smitha, D. Suhanya, S. Sridhar, M. Ramakrishna, "Separation of organic-organic mixtures by pervaporation - a review," Journal of Membrane Science, vol. 241, no. 1, pp. 1-21, 2004. View Article

[41]    L. M. Vane, "A review of pervaporation for product recovery from biomass fermentation processes," Journal of Chemical Technology & Biotechnology, vol.  80, no. 6, pp. 603-629, 2005. View Article

[42]    J. Villaluenga, A. Tabe-Mohammadi, "A review on the separation of benzene/cyclohexane mixtures by pervaporation processes," Journal of Membrane Science, vol. 169, no. 2, pp. 159-174, 2000.

[43]    J. G. Wijmans, R. W. Baker, "The solution-diffusion model: a review," Journal of Membrane Science, vol. 107, pp. 1-21, 1995. View Article

[44]    H. Karlsson, G. Tragardh, "Pervaporation of Dilute Organic-Waters Mixtures - a Literature-Review on Modeling Studies and Applications to Aroma Compound Recovery," Journal of Membrane Science, vol. 76, pp. 121-146, 1993. View Article

[45]    F. Lipnizki, G. Tragardh, "Modelling of pervaporation: Models to analyze and predict the mass transport in pervaporation," Separation and Purification Methods, vol. 30, pp. 49-125, 2001. View Article

[46]    X. S. Feng, R. Y. M. Huang, "Liquid separation by membrane pervaporation: A review," Industrial & Engineering Chemistry Research, vol. 36, no. 4, pp. 1048-1066, 1997. View Article

[47]    P. A. Kober, "Pervaporation, Perstillation and Percrystallization," Journal of the American-Chemical-Society, vol. 39, p. 944, 1917. Reprinted in Journal of Membrane Science, vol. 100, no. 1, pp. 61-64, 1995. View Article

[48]    A. Jonquières, C. Arnal-Herault, J. Babin, E. M. V. Hoek, V. V. Tarabara, "Pervaporation," Encyclopedia of Membrane Science and Technology, John Wiley & Sons, Inc., 2013. View Article

[49]    J. Néel, P. Aptel, R. Clément, "Basic aspects of pervaporation," Desalination, vol. 53, pp. 297-326, 1985. View Article

[50]    F. Lipnizki, R. W. Field, "Integration of vacuum and sweep gas pervaporation to recover organic compounds from wastewater," Separation and Purification Technology, vol. 22-23, pp. 347-360, 2001. View Article

[51]    S. M. Zhang, E. Drioli, "Pervaporation Membranes," Separation Science and Technology, vol. 30, pp. 1-31, 1995. View Article

[52]    E. J. Jeon, A. S. Kim, Y. T. Lee, "Pervaporation of butanol/water mixtures using siloxane polymer/ceramic composite membranes," Desalination and Water Treatment, vol. 48, pp. 17-26, 2012. View Article

[53]    F. Lipnizki, S. Hausmanns, G. Laufenberg, R. Field, B. Kunz, "Use of pervaporation/bioreactor hybrid processes in biotechnology," Chemie Ingenieur Technik, vol. 70, no. 12, pp. 1587-1595, 1998). View Article

[54]    J. Kujawa, S. Cerneaux, W. Kujawski, "Removal of hazardous volatile organic compounds from water by vacuum pervaporation with Hydrophobic ceramic membranes," Journal of Membrane Science, vol. 474, pp. 11-19, 2015. View Article

[55]    H. Zhou, D. Korelskiy, T. Leppäjärvi, M. Grahn, J. Tanskanen, J. Hedlund, "Ultrathin zeolite X membranes for pervaporation dehydration of ethanol," Journal of Membrane Science, vol. 399-400  pp. 106-111, 2012. View Article

[56]    Q. Liu, R. D. Noble, J. L. Falconer, H. H. Funke, "Organics/water separation by pervaporation with a zeolite membrane," Journal of Membrane Science, vol. 117, no.1-2, pp. 163-174, 1996. View Article

[57]    T. C. Bowen, R. D. Noble, J. L. Falconer, "Fundamentals and applications of pervaporation through zeolite membranes," Journal of Membrane Science, vol. 245, no. 1-2, pp. 1-33, 2004. View Article

[58]    S.-L. Wee, C.-T. Tye, S. Bhatia, "Membrane separation process—Pervaporation through zeolite membrane," Separation and Purification Technology, vol. 63, no. 3, pp. 500-516, 2008. View Article

[59]    H. Essawy, M. Tawfik, S. El-Sabbagh, A. El-Gendi, E. El-Zanati, H. Abdallah, "Novel amphiphilic conetworks based on compatibilized NBR/SBR–montmorillonite nanovulcanizates as membranes for dehydrative pervaporation of water–butanol mixtures," Polymer Engineering & Science, vol. 54, no. 7, pp. 1560-1570, 2014. View Article

[60]    K. D. Jitesh, V. G. Pangarkar, K. Niranjan, "Pervaporative stripping of acetone, butanol and ethanol to improve ABE fermentation," Bioseparation, vol. 9, no. 3, pp. 145-154, 2000. View Article

[61]    D. L. Vrana, M. M. Meagher, R. W. Hutkins, B. Duffield, "Pervaporation of Model Acetone-Butanol-Ethanol Fermentation Product Solutions using Polytetrafluoroethylene Membranes," Separation Science and Technology, vol. 28, no. 13-14, pp. 2167-2178, 1993. View Article

[62]    M. Müller, M.-N. Pons, "Coupling of gas membrane smooth pervaporation and alcoholic fermentation," Journal of Chemical Technology & Biotechnology, vol. 52, no. 3, pp. 343-358, 1991. View Article

[63]    M. Hu, L. Gao, W. Fu, X. Liu, F. Huang, Y. Luo, C. Huang, "High-performance interpenetrating polymer network polyurethane pervaporation membranes for butanol recovery," Journal of Chemical Technology & Biotechnology, vol. 90, no. 12, pp. 2195-2207, 2015. View Article

[64]    F. Liu, L. Liu, X. Feng, "Separation of acetone–butanol–ethanol (ABE) from dilute aqueous solutions by pervaporation," Separation and Purification Technology, vol. 42, no. 3, pp. 273-282, 2005. View Article

[65]    T. A. Jalal, N. M. S. Bettahalli, N. L. Le, S. P. Nunes, "Hydrophobic Hyflon AD/Poly(vinylidene fluoride) Membranes for Butanol Dehydration via Pervaporation," Industrial & Engineering Chemistry Research, vol. 54, no. 44, pp. 11180-11187, 2015. View Article

[66]    P. J. Hickey, F. P. Juricic, C. S. Slater, "The Effect of Process Parameters on the Pervaporation of Alcohols through Organophilic Membranes," Separation Science and Technology, vol. 27, no. 7, pp. 843-861, 1992. View Article

[67]    I. L. Borisov, A. O. Malakhov, V. S. Khotimsky, E. G. Litvinova, E. S. Finkelshtein, N. V. Ushakov, V. V. Volkov, "Novel PTMSP-based membranes containing elastomeric fillers: Enhanced 1-butanol/water pervaporation selectivity and permeability," Journal of Membrane Science, vol. 466 pp. 322-330, 2014. View Article

[68]    Y. Wang, T. S. Chung, H. Wang, "Polyamide–imide membranes with surface immobilized cyclodextrin for butanol isomer separation via pervaporation," AIChE Journal, vol. 57, no. 6, pp. 1470-1484, 2011. View Article

[69]    G. Liu, W. Wei, W. Jin, "Pervaporation Membranes for Biobutanol Production," ACS Sustainable Chemistry & Engineering, vol. 2, no. 4, pp. 546-560, 2014. View Article

[70]    X. Wang, J. Chen, M. Fang, T. Wang, L. Yu, J. Li, "ZIF-7/PDMS mixed matrix membranes for pervaporation recovery of butanol from aqueous solution," Separation and Purification Technology, vol. 163, pp. 39-47, 2016. View Article

[71]    H. Azimi, H. F. Tezel, J. Thibault, "Effect of embedded activated carbon nanoparticles on the performance of Polydimethylsiloxane (PDMS) membrane for pervaporation separation of butanol," Journal of Chemical Technology & Biotechnology, vol. 92, no. 12, 2017. View Article

[72]    A. Jonquières, A. Fane, "Filled and unfilled composite GFT PDMS membranes for the recovery of butanols from dilute aqueous solutions: influence of alcohol polarity," Journal of Membrane Science, vol. 125, no. 2, pp. 245-255, 1997. View Article

[73]    P. Peng, B. Shi, Y. Lan, "Preparation of PDMSSilica Nanocomposite Membranes with Silane Coupling for Recovering Ethanol by Pervaporation," Separation Science and Technology, vol. 46, no. 3, pp. 420-427, 2011. View Article

[74]    E. Drioli, S. Zhang, A. Basile, "On the coupling effect in pervaporation," Journal of Membrane Science, vol. 81, no. 1-2, pp. 43-55, 1993. View Article

[75]    S. Tan, L. Li, Z. Xiao, Y. Wu, Z. Zhang, "Pervaporation of alcoholic beverages - the coupling effects between ethanol and aroma compounds," Journal of Membrane Science, no. 264, no. 1-2, pp. 129-136, 2005. View Article

[76]    O. Kedem, "The Role of Coupling in Pervaporation," Journal of Membrane Science, vol. 47, no. 3, 277-284, 1989. View Article

[77]    M. She, S. Hwang, "Effects of concentration, temperature, and coupling on pervaporation of dilute flavor organics," Journal of Membrane Science, vol. 271, no. 1-2, pp. 16-28, 2006. View Article

[78]    M. Franke, "Auslegung und Optimierung von Pervaporationsanlagen zur Entwässerung von Lösungsmitteln und Lösungsmittelgemischen," Ph.D. Thesis, University of Aachen, 1990.

[79]    X. Feng, R. Y. M. Huang, "Concentration polarization in pervaporation separation processes," Journal of Membrane Science, vol. 92, no. 3 pp. 201-208, 1994. View Article

[80]    C.-C. Hu, T.-C. Liu, K.-R. Lee, R.-C. Ruaan, J.-Y. Lai, "Zeolite-filled PMMA composite membranes: influence of coupling agent addition on gas separation properties," Desalination, vol. 193, no. 1-3, pp. 14-24, 2006. View Article

[81]    G. M. Monsalve-Bravo, S. K. Bhatia, "Extending effective medium theory to finite size systems: Theory and simulation for permeation in mixed-matrix membranes," Journal of Membrane Science, vol. 531, pp. 148-159, 2017. View Article

[82]    D. Hofmann, C. Schepers, "Molecular simulation study on sorption and diffusion processes in polymeric pervaporation membrane materials," Molecular Simulation, vol. 32, no. 2, pp. 73-84, 2006. View Article

[83]    H. Okuno, K. Renzo, T. Uragami, "Sorption and permeation of water and ethanol vapors in poly(vinylchloride) membrane," Journal of Membrane Science, vol. 103, pp. 31-38, 1995. View Article

[84]    W. J. Koros, "Model for sorption of mixed gases in glassy polymers," Journal of Polymer Science: Polymer Physics Edition, vol. 18, no. 5, pp. 981-992, 1980. View Article

[85]    M. Buckley-Smith, "The Use of Solubility Parameters to Select Membrane Materials for Pervaporation of Organic Mixtures," Thesis, The University of Waikato, 2006. View Article

[86]    A. F. M. Barton, CRC Handbook of Solubility Parameters and Other Cohesion Parameters, 2nd ed. CRC Press, 1991.

[87]    C. M. Hansen, The three dimensional solubility parameter and solvent diffusion coefficient: Their importance in surface coating formulation. Danish Technical Press, 1967.

[88]    H. Zhou, Y. Su, X. Chen, Y. Wan, "Separation of acetone, butanol and ethanol (ABE) from dilute aqueous solutions by silicalite-1/PDMS hybrid pervaporation membranes," Separation and Purification Technology, vol. 79, no. 3, pp. 375-384, 2011. View Article

[89]    A. Hasanoglu, Y. Salt, S. Keleser, S. Ozkan, S. Dincer, "Pervaporation separation of organics from multicomponent aqueous mixtures," Chemical Engineering and Processing: Process Intensification, vol. 46, no. 4, pp. 300-306, 2007. View Article

[90]    J. Niemisto, W. Kujawski, R. L. Keiski, "Pervaporation performance of composite poly(dimethyl siloxane) membrane for butanol recovery from model solutions," Journal of Membrane Science, vol. 434, pp. 55-64, 2013. View Article

[91]    E. T. Zellers, "Three-dimensional solubility parameters and chemical protective clothing permeation. I. Modeling the solubility of organic solvents in Viton® golves," Journal of Applied Polymer Science, no. 50, no. 3, pp. 513-530, 1993. View Article

[92]    P. E. Froehling, D. M. Koenhen, A. Bantjes, C. A. Smolders, "Swelling of linear polymers in mixed swelling agents; predictability by means of solubility parameters," Polymer, vol. 17, no. 9, pp. 835-836, 1976. View Article

[93]    A. W. Islam, A. Javvadi, V. N. Kabadi, "Universal Liquid Mixture Models for Vapor−Liquid and Liquid−Liquid Equilibria in the Hexane−Butanol−Water System," Industrial & Engineering Chemistry Research, vol. 50, no. 2, pp. 1034-1045, 2011. View Article

[94]    Online Services - DDBST GmbH, (n.d.). [Online]. Available: http://www.ddbst.com/online.html [March 23, 2017].

[95]    M. Yoshikawa, H. Yokoi, K. Sanui, N. Ogata, "Pervaporation of water-ethanol mixture through poly(maleimide-co-acrylonitrile) membrane," Journal of Polymer Science: Polymer Letters Edition, vol.  22, no. 2, pp. 125-127, 1984. View Article

[96]    M. Yoshikawa, T. Yukoshi, K. Sanui, N. Ogata, "Separation of water and ethanol by pervaporation through poly(acrylic acid-co-acrylonitrile) membrane," Journal of Polymer Science: Polymer Letters Edition, vol. 22, no. 9, pp. 473-475, 1984. View Article

[97]    M. Yoshikawa, H. Yokoi, K. Sanui, N. Ogata, "Selective separation of water–alcohol binary mixture through poly(maleimide-co-acrylonitrile) membrane," Journal of Polymer Science: Polymer Chemistry Edition, vol. 22, no. 9, pp. 2159-2168, 1984. View Article

[98]    K. Boddeker, "Pervaporation durch Membranen und ihre Anwendung zur Trennung von Flüssiggemischen," Fortsch Ber VDI Reihle 3 Nr VDI-Verl. 129, 1986.

[99]    M. H. V. Mulder, C. A. Smolders, "On the mechanism of separation of ethanol/water mixtures by pervaporation I. Calculations of concentration profiles," Journal of Membrane Science, vol. 17, no. 3, pp. 289-307, 1984. View Article

[100]  R. B. Long, "Liquid Permeation through Plastic Films," Industrial & Engineering Chemistry Fundamentals, vol. 4, no. 4, pp. 445-451, 1965. View Article

[101]  F. W. Greenlaw, W. D. Prince, R. A. Shelden, E. V. Thompson, "Dependence of Diffusive Permeation Rates on Upstream and Downstream Pressures .1. Single Component Permeant," Journal of Membrane Science, vol. 2, pp. 141-151, 1977. View Article

[102]  R. A. Shelden, E. V. Thompson, "Dependence of diffusive permeation rates on upstream and downstream pressures," Journal of Membrane Science, vol. 4, pp. 115-127, 1978. View Article

[103]  A. K. Doolittle, "Studies in Newtonian Flow. II. The Dependence of the Viscosity of Liquids on Free‐Space," Journal of Applied Physics, vol. 22, pp. 1471-1475, 1951. View Article

[104]  A. K. Doolittle, "Studies in Newtonian Flow. III. The Dependence of the Viscosity of Liquids on Molecular Weight and Free Space (in Homologous Series)," Journal of Applied Physics, vol. 23, pp. 236-239, 1952. View Article

[105]  A. Kreituss, H. L. Frisch, "Free-volume estimates in heterogeneous polymer systems. I. Diffusion in crystalline ethylene–propylene copolymers," Journal of Polymer Science: Polymer Physics Edition, vol. 19, pp. 889-905, 1981. View Article

[106]  H. Fujita, "Diffusion in polymer-diluent systems," in: Fortschritte Der Hochpolymeren-Forschung, Springer, Berlin, Heidelberg, 1961, pp. 1-47. View Article

[107]  C. K. Yeom, R. Y. M. Huang, "Modelling of the pervaporation separation of ethanol-water mixtures through crosslinked poly(vinyl alcohol) membrane," Journal of Membrane Science, vol. 67, no. 1, pp. 39-55, 1992. View Article

[108]  İ. Şentarlı, S. Dinçer, Ö.T. Savaşçı, "Pervaporation Separation of Methanol-Carbontetrachloride Binary Mixtures Using Low Density Polyethylene Membranes," Chemical Engineering & Technology, vol. 21, no. 4, pp. 359-364, 1998.

[109]  W. J. Koros, R. T. Chern, V. Stannett, H. B. Hopfenberg, A model for permeation of mixed gases and vapors in glassy polymers, Journal of Polymer Science: Polymer Physics Edition, vol. 19, no. 10, pp. 1513-1530, 1981. View Article

[110]  Y. Tamai, H. Tanaka, K. Nakanishi, "Molecular Simulation of Permeation of Small Penetrants through Membranes. 1. Diffusion Coefficients," Macromolecules, vol. 27, no. 16, pp. 4498-4508, 1994. View Article

[111]  L. Fritz, D. Hofmann, "Molecular dynamics simulations of the transport of water-ethanol mixtures through polydimethylsiloxane membranes," Polymer, vol. 38, no. 5, pp. 1035-1045, 1997. View Article

[112]  E. El-Zanati, E. Abdel-Hakim, O. El-Ardi, M. Fahmy, "Modeling and simulation of butanol separation from aqueous solutions using pervaporation," Journal of Membrane Science, vol. 280, no. 1-2, pp. 278-283, 2006. View Article

[113]  H. R. Cascon, S. K. Choudhari, "1-Butanol pervaporation performance and intrinsic stability of phosphonium and ammonium ionic liquid-based supported liquid membranes," Journal of Membrane Science, vol. 429, pp. 214-224, 2013. View Article

[114]  O. Vopička, V. Hynek, K. Friess, P. Izák, "Blended silicone–ionic liquid membranes: Transport properties of butan-1-ol vapor," European Polymer Journal, vol. 46, no. 1, pp. 123-128, 2010. View Article

[115]  X. Yang, Z. Wu, F. Manquan, L. Jiding, "Nonequilibrium Dissolution-diffusion Model for PDMS Membrane Pervaporation of ABE Water Binary System," Journal of Membrane Science & Technology, vol. 6, no. 1, 2016. View Article

[116]  S.-Y. Li, R. Srivastava, R. S. Parnas, "Separation of 1-butanol by pervaporation using a novel tri-layer PDMS composite membrane," Journal of Membrane Science, vol. 363, no. 1-2, pp. 287-294, 2010. View Article

[117]  N. Valentínyi, E. Cséfalvay, P. Mizsey, "Modelling of pervaporation: Parameter estimation and model development," Chemical Engineering Research and Design, vol. 91, no. 1, pp. 174-183, 2013. View Article

[118]  R. W. Baker, "Membrane Transport Theory," in: Membrane Technology and Applications, 2nd ed, John Wiley & Sons, Ltd, 2004, pp. 15-87. View Article

[119]  M. T. Ashraf, J. E. Schmidt, J. Kujawa, W. Kujawski, H. A. Arafat, "One-dimensional modeling of pervaporation systems using a semi-empirical flux model," Separation and Purification Technology, vol. 174, pp. 502-512, 2017. View Article

[120]  R. Petrychkovych, K. Setnickova, P. Uchytil, "The influence of water on butanol isomers pervaporation transport through polyethylene membrane," Separation and Purification Technology, vol. 107, pp. 85-90, 2013. View Article

[121]  R. Krishna, L. J. P. van den Broeke, "The Maxwell-Stefan description of mass transport across zeolite membranes," The Chemical Engineering Journal and the Biochemical Engineering Journal, vol. 57, no. 2, pp. 155-162, 1995. View Article

[122]  F. Kapteijn, J. A. Moulijn, R. Krishna, "The generalized Maxwell–Stefan model for diffusion in zeolites:: sorbate molecules with different saturation loadings," Chemical Engineering Science, vol. 55  pp. 2923-2930, 2000. View Article

[123]  A. Heintz, W. Stephan, "A generalized solution—diffusion model of the pervaporation process through composite membranes Part II. Concentration polarization, coupled diffusion and the influence of the porous support layer," Journal of Membrane Science, vol. 89, no. 1-2, pp. 153-169, 1994. View Article

[124]  P. Izák, L. Bartovská, K. Friess, M. Šı́pek, P. Uchytil, "Description of binary liquid mixtures transport through non-porous membrane by modified Maxwell–Stefan equations," Journal of Membrane Science, vol. 214, no. 2, pp. 293-309, 2003. View Article

[125]  R. Krishna, J. A. Wesselingh, "The Maxwell-Stefan approach to mass transfer," Chemical Engineering Science, vol. 52, no. 6, pp. 861-911, 1997. View Article

[126]  A. Ebneyamini, H. Azimi, J. Thibault, F. H. Tezel, "Description of Butanol Aqueous Solution Transport through Commercial PDMS Membrane via Pervaporation Using Extended Maxwell–Stefan Model," J. Sep. Sci. Technol., 2017. Submitted.

[127]  J.-J. Shieh, R. Y. M. Huang, "A Pseudophase-Change Solution-Diffusion Model for Pervaporation. I. Single Component Permeation," Separation Science and Technology, vol. 33, pp. 767-785, 1998. View Article

[128]  J.-J. Shieh, R. Y. M. Huang, "A Pseudophase-Change Solution-Diffusion Model for Pervaporation. II. Binary Mixture Permeation," Separation Science and Technology, vol. 33, pp. 933-957, 1998. View Article

[129]  J. R. Snider, G. A. Dawson, "Tropospheric light alcohols, carbonyls, and acetonitrile: Concentrations in the southwestern United States and Henry’s Law data," Journal of Geophysical Research: Atmospheres, vol. 90, no. D2, pp. 3797-3805, 1985. View Article

[130]  P. Vrbka, D. Fenclová, V. Laštovka, V. Dohnal, "Measurement of infinite dilution activity coefficients of 1-alkanols (C1–C5) in water as a function of temperature (273–373 K)," Fluid Phase Equilibria, vol. 237, no. 1-2, , pp. 123-129, 2005. View Article

[131]  B. V. Ioffe, A. G. Vitenberg, "Application of the headspace analysis to systems with unknown partition coefficients," Chromatographia, vol. 11, no. 5, pp. 282-286, 1978. View Article

[132]  N. H. Sagert, D. W. P. Lau, "Limiting activity coefficients for butyl alcohols in water, n-octane, and carbon tetrachloride," Journal of Chemical & Engineering Data, vol. 31, no. 4, pp. 475-478, 1986. View Article

[133]  A. Lebert, D. Richon, "Infinite dilution activity coefficients of n-alcohols as a function of dextrin concentration in water-dextrin systems," Journal of Agricultural and Food Chemistry, vol. 32, no. 5,  pp. 1156-1161, 1984. View Article

[134]  J. A. V. Butler, C. N. Ramchandani, D. W. Thomson, "58. The solubility of non-electrolytes. Part I. The free energy of hydration of some aliphatic alcohols," Journal of the Chemical Society, pp. 280-285, 1935. View Article

[135]  G. J. Pierotti, C. H. Deal, E. L. Derr, "Activity Coefficients and Molecular Structure," Industrial & Engineering Chemistry, vol. 51, no. 1, pp. 95-102, 1959. View Article

[136]  M. G. Burnett, "Determination of Partition Coefficients at Infinite Dilution by the Gas Chromatographic Analysis of the Vapor above Dilute Solutions," Analytical Chemistry, vol. 35, no. 11, pp. 1567-1570, 1963. View Article

[137]  R. G. Buttery, L. Ling, D. G. Guadagni, "Food volatiles. Volatilities of aldehydes, ketones, and esters in dilute water solution," Journal of Agricultural and Food Chemistry, vol. 17, pp. 385-389, 1969. View Article

[138]  J. Howard Rytting, L. P. Huston, T. Higuchi, "Thermodynamic Group Contributions for Hydroxyl, Amino, and Methylene Groups," Journal of Pharmaceutical Sciences, vol. 67, no. 5, pp. 615-618 1978. View Article

[139]  I. Landau, A. J. Belfer, D. C. Locke, "Measurement of limiting activity coefficients using non-steady-state gas chromatography," Industrial & Engineering Chemistry Research, vol. 30, no. 8, pp. 1900-1906, 1991. View Article

[140]  J. Li, P. W. Carr, "Measurement of water-hexadecane partition coefficients by headspace gas chromatography and calculation of limiting activity coefficients in water," Analytical Chemistry, vol. 65, no. 10, pp. 1443-1450, 1993. View Article

[141]  S. Merk, M. Riederer, "Sorption of volatile C1 to C6 alkanols in plant cuticles," Journal of Experimental Botany, vol. 48, no. 5, pp. 1095-1104, 1997. View Article

[142]  P. G. Whitehead, S. I. Sandler, "Headspace gas chromatography for measurement of infinite dilution activity coefficients of C4 alcohols in water," Fluid Phase Equilibria, vol. 157, no. 1, 111-120, 1999. View Article

[143]  J. Altschuh, R. Brüggemann, H. Santl, G. Eichinger, O. G. Piringer, "Henry’s law constants for a diverse set of organic chemicals: Experimental determination and comparison of estimation methods," Chemosphere, vol. 39,  no. 11, pp. 1871-1887, 1999. View Article

[144]  B. R. Kim, E. M. Kalis, T. Dewulf, K. M. Andrews, "Henry’s Law Constants for Paint Solvents and Their Implications on Volatile Organic Compound Emissions from Automotive Painting," Water Environment Research, vol. 72, no. 1, pp. 65-74, 2000. View Article

[145]  K. Tochigi, M. Uchiyama, K. Kojima, "Measurement of infinite-dilution activity coefficients of alcohols in water using relative gas-liquid chromatographic method," Korean Journal of Chemical Engineering, vol. 17, pp. 502-505, 2000. View Article

[146]  L. T. Iraci, B. M. Baker, G. S. Tyndall, J. J. Orlando, "Measurements of the Henry’s Law Coefficients of 2-Methyl-3-buten-2-ol, Methacrolein, and Methylvinyl Ketone," Journal of Atmospheric Chemistry, vol. 33, no. 3, pp. 321-330, 1999. View Article

[147]  S. L. Friant, I. H. Suffet, "Interactive effects of temperature, salt concentration, and pH on head space analysis for isolating volatile trace organics in aqueous environmental samples," Analytical Chemistry, vol. 51, no. 13, pp. 2167-2172, 1979. View Article

[148]  M. Hofstee, A. Kwantes, G. W. A. Rijnders, "Determination of activity coefficients at infinite dilution by gas–liquid chromatography," Proceedings of the International Symposium on Distillation, Brighton Engl. pp. 105-109, 1960.

[149]  B. Kolb, C. Welter, C. Bichler, "Determination of partition coefficients by automatic equilibrium headspace gas chromatography by vapor phase calibration," Chromatographia, vol. 34, no. 5-8, pp. 235-240, 1992. View Article

[150]  A. K. Gupta, A. S. Teja, X. S. Chai, J. Y. Zhu, "Henry’s constants of n-alkanols (methanol through n-hexanol) in water at temperatures between 40°C and 90°C," Fluid Phase Equilibria, vol. 170, no. 2, pp. 183-192, 2000. View Article

[151]  A. Klamt, F. Eckert, "COSMO-RS: a novel and efficient method for the a priori prediction of thermophysical data of liquids," Fluid Phase Equilibria, vol. 172, no. 1,  pp. 43-72, 2000. View Article

[152]  Š. Hovorka, V. Dohnal, A.H. Roux, G. Roux-Desgranges, "Determination of temperature dependence of limiting activity coefficients for a group of moderately hydrophobic organic solutes in water," Fluid Phase Equilibria, vol. 201, no. 1, pp. 135-164, 2002. View Article

[153]  K. Fischer, J. Gmehling, "P-x and γ∞ data for the different binary butanol-water systems at 50 °C," Journal of Chemical & Engineering Data, vol. 39, no. 2, pp. 309-315, 1994. View Article

[154]  V. Dohnal, D. Ondo, "Refined non-steady-state gas–liquid chromatography for accurate determination of limiting activity coefficients of volatile organic compounds in water: Application to C1–C5 alkanols," Journal of Chromatography A, vol. 1097, no. 1-2, pp. 157-164, 2005. View Article

[155]  K. Tochigi, K. Kojima, "The Determination of Group Wilson Param Eters to Activity Coefficients by Ebulliometer," Journal of Chemical Engineering Of Japan, vol. 9, no. 4, pp. 267-273, 1976. View Article

[156]  G. M. Lobien, J. M. Prausnitz, "Infinite-dilution activity coefficients from differential ebulliometry," Industrial & Engineering Chemistry Fundamentals, vol. 21, no. 2, pp. 109-113, 1982. View Article

[157]  K. Ochi, K. Kojima, A Measurement of Vapor-Liquid Equilibria at Extreme Dilution, J. Chem. Eng. Jpn. 20 (1987) 6–10. View Article

[158]  G. Cocchi, M. G. De Angelis, F. Doghieri, "Solubility and diffusivity of liquids for food and pharmaceutical applications in crosslinked polydimethylsiloxane (PDMS) films: I. Experimental data on pure organic components and vegetable oil," Journal of Membrane Science, vol. 492, pp. 600-611, 2015. View Article

[159]  P. Radovanovic, S.W. Thiel, S.-T. Hwang, Transport of ethanol-water dimers in pervaporation through a silicone rubber membrane, Journal of Membrane Science, vol. 48, pp. 55–65, 1990. View Article

[160]  E. Favre, P. Schaetzel, Q. T. Nguygen, R. Clément, J. Néel, "Sorption, diffusion and vapor permeation of various penetrants through dense poly(dimethylsiloxane) membranes: a transport analysis," Journal of Membrane Science, vol. 92, no. 2, pp. 169-184, 1994. View Article

[161]  J. M. Watson, M. G. Baron, "The behaviour of water in poly(dimethylsiloxane)," Journal of Membrane Science, vol. 110, no. 1, pp. 47-57, 1996. View Article

[162]  J. M. Watson, P. A. Payne, "A study of organic compound pervaporation through silicone rubber," Journal of Membrane Science, vol. 49, no. 2, pp. 171-205, 1990. View Article

[163]  J. M. Watson, G. S. Zhang, P. A. Payne, "The diffusion mechanism in silicone rubber," Journal of Membrane Science, vol. 73, no. 1, 55-71, 1992. View Article

[164]  C. Zeng, J. Li, D. Wang, T. Chen, C. Zhao, C. Chen, "Infinite Dilute Activity and Diffusion Coefficients in Polymers by Inverse Gas Chromatography," Journal of Chemical & Engineering Data, vol. 51 pp. 93-98, 2006. View Article

[165]  Meintjes, M, Fermentation coupled with pervaporation: a kinetic study, (2011). [Online]. Available: http://www.academia.edu/27356851/Fermentation_coupled_with_pervaporation_a_kinetic_study (September 9, 2016). View Article

[166]  V. S. Nayak, J. B. Moffat, "Sorption and diffusion of alcohols in heteropoly oxometalates and ZSM-5 zeolite," The Journal of Physical Chemistry, vol. 92, no. 25, pp. 7097-7102, 1988. View Article


Table 1. Summary of the updated data in the literature for infinite dilution activity coefficients of butanol in water.

T (K) g Method Ref
273.15 40.21 HSA [129]
273.15 32.59 HSA [130]
283.15 40.41 HSA [130]
288.15 49.01 IGS [131]
290 20.72 MOSCED [94]
293.15 46.48 IGS [131]
293.15 41.39 HSA [132]
293.15 48.38 HSA [130]
298.15 46.48 IGS [133]
298.15 52.88 DDST [134]
298.15 48.18 TENS [135]
298.15 44.48 IGS [136]
298.15 54.98 HSA [137]
298.15 52.83 HSA [138]
298.15 47.09 IGS [131]
298.15 49.45 HSA [129]
298.15 53.68 NSGLC [139]
298.15 53.30 HSA [140]
298.15 58.85 HSA [141]
298.15 48.18 HSA [142]
298.15 53.62 WWC [143]
298.15 57.40 IGS [144]
298.15 56.90 UNIQUAC [93]
298.15 58.10 LSG [93]
298.15 56.20 GEM-RS [93]
298.15 55.80 NRTL [93]
298.15 51.37 GLC [93]
298.2 55.20 GLC [145]
298.45 51.11 HSA [146]
300 18.89 MOSCED [94]
303.15 81.61 HSA [147]
303.15 51.68 HSA [142]
303.15 53.41 TENS [130]
308.15 39.49 HSA [142]
308.15 57.23 TENS [130]
310 17.16 MOSCED (94)
313.15 49.50 GLC [148]
313.15 89.57 VPC [149]
313.15 73.19 PRV [150]
313.15 54.60 COSMO-RS [151]
318.15 60.22 TENS [130]
320 15.55 MOSCED [94]
323.15 73.11 PRV [150]
323.15 61.37 RDIST [152]
323.15 62.30 RDIST [152]
323.15 61.50 TENS [130]
323.15 58.67 HSA [130]
323.15 59.38 HSA [130]
323.2 58.67 GLC [145]
323.23 78.73 TENS [153]
328.15 61.62 NSGLC [154]
328.15 62.80 TENS [130]
330 14.09 MOSCED [94]
333.15 47.80 TENS [135]
333.15 59.32 GLC [148]
333.15 79.28 VPC [149]
333.15 73.70 PRV [150]
333.15 66.69 RDIST [152]
333.15 64.39 RDIST [152]
333.15 65.10 RDIST [152]
333.15 63.37 TENS [130]
333.2 61.13 GLC [145]
340 12.76 MOSCED [94]
343.15 59.32 EBUL [155]
343.15 67.83 EBUL [156]
343.15 79.68 VPC [149]
343.15 75.04 PRV [150]
343.15 63.82 RDIST [152]
353.15 57.23 EBUL [155]
353.15 46.48 EBUL [156]
353.15 73.19 VPC [149]
353.15 76.55 PRV [150]
353.15 62.80 RDIST [152]
363.15 55.48 EBUL [155]
363.15 79.52 PRV [150]
372.15 27.11 EBUL [156]
372.15 56.09 CIRC [130]
373.15 54.00 EBUL [155]
373.15 88.23 EBUL [157]

Table 2. Diffusion coefficients of ABE components reported in the literature for different pervaporation membranes.

Polymer/Solution Component Temperature(K) Diffusivity (m2/s) Ref
PDMS/Pure component Acetone 308 0.85x10-10 [158]
n-Butanol 308 0.3x10-10
Ethanol 308 0.60x10-10
Water 308 0.55x10-10
PDMS/Ethanol-Water mixture Ethanol 298.15 0.6x10-9 [159]
Water 298.15 4.4x10-9
Ethanol-water 298.15 0.37x10-9
PDMS/Pure component n-Butanol 303.15 3.11x10-10* [160]
t-Butanol 303.15 2.66x10-10*
s-Butanol 303.15 2.25x10-10*
Water 303.15 129.05x10-10*
PDMS (GE615)/Pure water Water 298.15 1.7x10-9 (a) [161]
PDMS (PS342.5)/Pure water Water 298.15 2x10-9 (a)
PDMS/Pure water Water 298.15 7.8x10-10 (a)
PDMS/Alcohols aqueous binary solutions n-Butanol 303.15 0.65x10-10 [162]
Ethanol 353.15 7.1x10-10*
n-Butanol 353.15 5.5x10-10
Ethanol 353.15 12x10-10
Water 353.15 12x10-10
Water 298.15 3x10-10
PDMS/Aqueous solutions n-Butanol 299.15 1.6x10 -10 (b) [163]
Crosslinked PVA/Dilute solutions n-Butanol 398.15 0.563x10-10* [164]
Ethanol 398.15 0.793x10-10*
Ethanol 393.15 0.787x10-10*
Ethanol 383.15 0.92x10-10*
Water 393.15 39.08x10-10*
Water 383.15 8.89x10-10*
PDMS/Pure component Ethanol 300 0.45x10-9 [110]
Water 300 1.45x10-9
PE/Pure Ethanol 300 0.07x10-9
Water 300 0.78x10-9
PDMS/Aqueous binary solutions Acetone 313.2 6.106x10-12 [115]
n-Butanol 313.2 2.589x10-12
Ethanol 313.2 2.05x10-12
Water 313.2 1.867x10-12
PERVAP®4060/Aqueous binary solutions Ethanol 303.15 9.55x10-9 [165]
Water 303.15 6.52x10-10
Silicalite-1 zeolite/Pure component n-Butanol 293 4.5x10-15 [166]
Ethanol 293 13.2x10-16
Silicalite-1 zeolite/Pure component Water 298 1.7x10-9 [57]
IL gel layer/ABE solutions Acetone 303.15 4.6x10-10_4.9x10-10 [33]
n-Butanol 303.15 4.2x10-11_3.9x10-10
Ethanol 303.15 5.6x10-12_5.2x10-10
Methylated silica/ Aqueous binary solutions n-Butanol 333.15 4.7x10−14( c ) [34]
Ethanol 333.15 1.1x10−13 ( c )
Water 333.15 3.7x10-12_7.1x10−12 (c )
* infinite dilute diffusion coefficient
(a) determined on the transient slope
(b) determined at half saturation of transient
(c) Maxwell–Stefan single-file diffusivities

Table 3. Mathematical models which have been used for the pervaporation separation of butanol.

Membrane Feed Modelling method Ref
Supported gelled ionic liquid ABE mixtures Solution diffusion [33]
Composite PVA membrane on PAN support butanol/water Solution diffusion [112]
Composite PDMS ABE fermentation Solution diffusion [18]
Tri-layer PDMS butanol/water Solution diffusion [116]
PVA/PAN water/butanol Solution diffusion [117]
PERVAP 2510 water/butanol Solution diffusion [119]
PE butanol/water Solution diffusion [120]
PDMS ABE water binary Non-equilibrium solution diffusion [115]
PDMS butanol/water Maxwell–Stefan [126]
Methylated microporous silica butanol/water, butanol/methanol Maxwell–Stefan [34]