1. Introduction
Membrane processes consume less energy compared to traditional separation methods, require less operational space, and generate fewer by-products, making them efficient and environmentally friendly[1-3]. Additionally, membrane technology has a wide range of applications, including gas separation, organic solvent separation, seawater desalination, packaging, and pharmaceuticals. The choice of membrane material depends on the type of substance to be separated. While several factors, such as chemical resistance and cost-effectiveness, must be considered when selecting membrane materials, the pore characteristics of the membrane are generally the primary determinant of which species permeate or are retained in a mixture.
Membranes are broadly classified into porous and non-porous types based on pore size. Non-porous membranes, used in applications such as gas separation, reverse osmosis, and pervaporation, exhibit pore structures that dynamically change according to the diffusion behavior of the permeating species. This can be explained by the solution-diffusion model. On the other hand, while porous membranes maintain their physical pore structure regardless of the diffusion behavior of the permeating species, there is no single comprehensive model that describes them. The pore-flow model, which explains the behavior of substances based on pressure gradients through fine pores, may or may not apply depending on the characteristics of the porous membrane.
Given these factors, the analysis techniques used to interpret and cross-verify the separation performance of porous membranes become crucial. In general, the separation mechanism of porous membranes involves molecular sieving, where only species small enough to physically pass through the pores can permeate. The pore characteristics as molecular sieves include not only pore size but also pore size distribution (PSD), porosity, tortuosity, and connectivity. Each of these characteristics requires different analytical methods. This paper focuses on methods for analyzing pore size distribution, one of the key characteristics of porous membranes.
2. Main Discussion
2.1. Types of porous membranes
As shown in Fig. 1, porous membranes can be classified into two types based on their structure: screen filters and depth filters. Screen filters have finely structured surface pores that filter particles at the membrane surface. These membranes exhibit an asymmetric structure, composed of a surface layer with small pores and a support layer with relatively larger pores, and are typically used in ultrafiltration membranes (pore size: 2~100 nm). This indicates that screen filters are suitable for filtering fluids with a relatively narrow particle size distribution. The filtration mechanism of screen filters can be explained using the Ferry-Renkin equation (Equation 1), which accounts for the pore size (r), the particle size (a), and the parabolic velocity profile of the fluid as it passes through the pores[4,5].
Equation 1 is used to estimate the pore size of the membrane based on the radius of the filtrate particles and the membrane rejection rate (R). In Fig. 2(a), an example is provided where the equation is applied to actual spherical proteins, plotting a rejection rate curve for the membrane. Spherical proteins were chosen as marker proteins because their simple shape allows for the calculation of particle size based solely on molecular weight. In Fig. 2(b) curve indicates the rejection as a function of the radius of the protein at a fixed pore size.
Depth filters, on the other hand, allow particles to pass through the membrane thickness, where filtration occurs within the membrane structure. Unlike screen filters, depth filters have a symmetric structure and are typically used in microfiltration membranes (pore size: 100~500 nm). The filtration mechanism of depth filters includes not only size sieving, where particles are filtered based on their size as they pass through the membrane, but also adsorption of particles within the membrane structure. The filtration mechanisms can be further categorized as follows: 1) size sieving, where particles are filtered based on size; 2) inertial capture, where larger particles are trapped against the membrane wall due to their inability to follow the fluid flow; 3) adsorption due to Brownian diffusion, where smaller particles adhere to the membrane surface due to irregular motion; and 4) electrostatic adsorption, where charged particles are captured based on the electrostatic properties of the membrane. Depth filters are suitable for filtering fluids with a wide range of particle size distributions, but in cases where electrostatic adsorption is used, the rejection rate may decrease if the membrane's charge performance deteriorates over time.
2.2. Types of pore characteristics to consider
In porous membranes, pore characteristics are directly linked to separation performance, making their determination and analysis crucial. However, the values associated with these characteristics may not always accurately represent membrane performance. As shown in Fig. 3(a), even if the average porosity of a screen filter membrane is 0.7~0.8, the porosity of the actual filtering layer (skin layer) may be as low as 0.05, significantly lower than the average. Additionally, porous membranes with a range of pore sizes have different standards for defining representative pore size, depending on the membrane type (See Fig. 3(b)). For ultrafiltration membranes, the mean pore size is typically used as the representative value, while for microfiltration membranes, the size of the largest particles that can permeate is used as the representative pore size. Consequently, the actual measured pore size for microfiltration membranes is often significantly smaller than the representative value.
A common indicator used as a substitute for pore characteristics is the molecular weight cutoff (MWCO). By gradually increasing the molecular weight of a substance and passing it through the membrane, an S-shaped curve can be obtained, showing the relationship between molecular weight and rejection rate (Fig. 4). The molecular weight at which the rejection rate reaches 90% is referred to as the MWCO of the membrane, and it is often used to assess the pore characteristics of ultrafiltration membranes.
However, MWCO cannot be considered an absolute indicator for several reasons[6]. First, the MWCO value can vary depending on the choice of marker protein, as the shape and size of proteins, even with the same molecular weight, can differ under test conditions. Second, the wide variance in pore size distribution of ultrafiltration membranes can make it difficult to precisely define the MWCO. Third, even for the same membrane, MWCO values can vary depending on how test conditions such as pressure and measurement time are set by the tester. Fourth, since ultrafiltration membranes are generally hydrophobic, proteins may adsorb onto the membrane surface or interior, and the extent of this adsorption can vary depending on factors like protein concentration, solvent pH, and ionic strength. If these variables are not controlled, the membrane's rejection rate or MWCO may not reflect its true performance. Therefore, it is crucial for testers to control the testing conditions and environment to accurately determine the MWCO, allowing for meaningful comparisons of membrane characteristics.
2.3. Analysis methods and theories for pore characteristics
Methods for analyzing the pore characteristics of porous membranes can be categorized into direct and indirect approaches. Direct methods, such as scanning electron microscopy (SEM), transmission electron microscopy (TEM), and atomic force microscopy (AFM), allow for the observation of the membrane’s surface and cross-sectional morphology[8-10]. These techniques provide real-time direct visualization of surface porosity, pore shape, and pore size. However, their limitation lies in their localized analysis, which requires a large number of samples to obtain representative data for the entire membrane.
In contrast, indirect methods infer pore size distribution based on theoretical models and measurement data. Indirect approaches can be further divided into three categories according to the underlying theory:
-
1) Methods based on the Young-Laplace equation, such as the bubble point method and liquid displacement method.
-
2) Techniques based on the Kelvin equation, including permoporometry, gas adsorption-desorption, and evapoporometry.
-
3) Methods based on the Gibbs-Thomson equation, such as thermoporometry.
This mini review focuses on how these theoretical models are applied to various analytical techniques and how each method is utilized for characterizing porous membranes.
2.3.1. Young-laplace equation
The Young-Laplace equation describes the pressure gradient across the interface between two immiscible fluids and is commonly used to interpret surface tension in capillary phenomena. Assuming the pores of the membrane are cylindrical in shape and the interface of two immiscible fluids (e.g., gas-liquid or liquid-liquid) is distributed within these pores, the pressure difference (Δp) across the interface can be expressed as follows:
Here, γ represents the interfacial tension, H is the curvature of the interface, θ is the contact angle between the interface and the pore wall, and d is the pore diameter. By analyzing the distribution of these pressure differences, the pore size distribution of the membrane can be determined.
2.3.1.1. Liquid displacement method
The liquid displacement method, based on the Young-Laplace equation, includes the bubble point method, the gas-liquid displacement method, and the liquid-liquid displacement method.
The bubble point method involves gradually increasing the pressure in all pores filled with liquid until the liquid is displaced by gas. The first pore completely displaced by gas (the first bubble) is assumed to be the largest pore in the membrane. Thus, the bubble point method is commonly used as a commercial technique to assess membrane integrity rather than to determine pore size distribution (i.e., to verify whether the pore sizes fall within the acceptable range)[11-14]. However, a drawback of this method is that high pressure is required if the gas-liquid interfacial tension is high. For instance, when using water as the wetting fluid, the water-gas interfacial tension is 73 mN/m, and approximately 15 bar of pressure is needed to analyze membranes with pore sizes up to 0.2 μm[13]. If the pressure is too high, it may deform the membrane structure, reducing measurement accuracy, so lower- tension liquids (e.g., halogenated fluids) are often used for analysis[15,16].
The gas-liquid displacement method involves gradually increasing the gas pressure after the bubble point is reached, allowing smaller pores filled with liquid to be replaced by gas. By measuring the gas pressure and flow rate, the pore size distribution of the membrane can be determined. To convert the pressure and flow rate data into pore sizes, the Hagen-Poiseuille equation is applied to describe convective flow through the pores. However, if the pore size is smaller than the mean free path of the gas, the measurement may overestimate the pore density[15]. Therefore, a combination of the Hagen-Poiseuille and Knudsen models is used to correct the pore size distribution[17]. For microfiltration membranes with pore sizes ranging from 0.1 to 0.5 μm, studies have shown that simply using the Hagen-Poiseuille behavior with a liquid having a higher interfacial tension than water (73 mN/m) is sufficient[18].
The typical procedure for measuring pore size distribution using the gas-liquid displacement method involves first immersing the membrane in liquid for a set period to ensure that the liquid fills the pores completely. Then, the membrane is placed in a dead-end filtration cell, and pressure is gradually increased. Once the first bubble appears, gas begins to permeate the membrane, and the pore size is determined based on the measured gas flow rate and applied pressure. The Young-Laplace equation is used here, and the choice of wetting fluid should match the hydrophilic or hydrophobic nature of the membrane to ensure accurate pore size calculation. However, since there is no clear standard for selecting a wetting fluid, halogenated fluids—characterized by low interfacial tension, low reactivity, and low vapor pressure—are commonly used. The gas-liquid displacement method, also known as Liquid Extrusion Porometry, Flow Permporometry, or Capillary Flow Porometry, divides the pressure application into approximately 200~300 increments, producing a pore size distribution resembling a normal distribution[19-21].
The advantages of the gas-liquid displacement method include its simplicity, non-destructive nature, and ability to analyze only interconnected pores. However, limitations include: 1) The assumption that pores are cylindrical, 2) Variation in pore size distribution based on the choice of wetting fluid, and 3) Difficulty in analyzing pores smaller than 50 nm[22]. Especially when the pore shape is irregular, the cylindrical assumption (i.e., assuming no tortuosity and using membrane thickness as pore length) can lead to inaccuracies[23].
The liquid-liquid displacement method is suitable for analyzing smaller pore sizes (ca. < 50 nm) than the gas-liquid method. Similar in principle, it uses the interface of two immiscible liquids to calculate pore size distribution. Unlike the gas-liquid method, it relies solely on the Hagen-Poiseuille model to describe convective flow, allowing for calculation of both pore numbers and sizes without the interference seen in gas-liquid behavior[24-26]. Additionally, because the liquid-liquid interfacial tension is lower than that of gas-liquid, it allows for analysis at relatively lower pressures[27].
In the liquid-liquid displacement method, two immiscible liquids are used: one acts as the wetting fluid, while the other (the displacing liquid) has low affinity with the membrane, displacing the wetting fluid from the pores. Similar to the bubble point method, when the displacing liquid begins to displace the wetting liquid at the threshold pressure, pore size is determined based on the flow rate and applied pressure. Using the Hagen-Poiseuille model, the flow rate (Jn) at each pressure step (Pn) is expressed as follows[23]:
where rn is the pore diameter at step n, Nn is the number of pores at step n, η is the dynamic viscosity of the displacing fluid, and l is the pore length.
The choice and combination of wetting and displacing liquids determine the effectiveness of the liquid-liquid displacement method. Initially, water-alcohol mixtures (e.g., butanol or propanol) with low interfacial tension (ca. 1.73 mN/m) were used[27,28]. However, alcohols tend to adsorb within polymeric membranes, causing swelling and reducing measurement accuracy [29]. Recently, ammonium sulfate solution (displacing fluid) and polyethylene glycol solution (wetting fluid) have been used to prevent membrane swelling[31]. Using water-soluble polymers and salts in aqueous systems offers the advantage of easy removal after analysis, but this method is limited to hydrophilic membranes and may lack reproducibility for smaller pore sizes due to the higher viscosity of both fluids compared to water[30-32].
To improve reproducibility, the displacing and wetting liquids in the liquid-liquid displacement method should possess the following characteristics: 1) Low liquid-liquid interfacial tension, 2) Low viscosity and vapor pressure, and 3) Low adsorption within the membrane material.
The liquid-liquid displacement method's advantages include being a non-destructive analysis, targeting only interconnected pores, and its applicability to ultrafiltration and nanofiltration membranes (pore size 1~10 nm)[33]. However, the method's limitations are: 1) The assumption of cylindrical pore shapes, and 2) The lack of a clear standard for selecting wetting and displacing liquids.
2.3.2. Kelvin equation
The Kelvin equation, similar to the Young-Laplace equation, describes the interface between two immiscible fluids but focuses on the thermodynamic equilibrium when the two fluids are in contact. As shown in Fig. 5, when the meniscus of the liquid phase is convex toward the gas phase, the thermodynamic equilibrium between the two phases can be expressed as follows:
Rearranging Equation 4, it can be expressed as:
By neglecting VL (P0 - PV) and substituting the Young-Laplace equation into Equation 5, the Kelvin equation is derived as follows:
2.3.2.1. Permporometry
Permporometry, based on the Kelvin equation, involves reducing the gas pressure below the saturation vapor pressure so that a condensable gas fills the membrane’s pores. The pressure is then gradually increased, causing the condensed gas in the pores to evaporate, and the flow rate of a non-condensable gas like helium, nitrogen, or oxygen is used to map the pore size distribution[34,35]. In typical setups, nitrogen and the condensable gas permeate from both the feed and permeate sides of the membrane, while oxygen is sent to the permeate side to measure flow rates using gas chromatography[36].
Condensable gas, such as methanol, ethanol, propanol, or hexane, must meet two criteria: 1) It should be a wetting and inert gas, and 2) It must have suitable vapor pressure and evaporation rates for precise analysis[37,38].
Permporometry is commonly used for ceramic membranes, measuring pore sizes from 4 to 50 nm, while its derivative, Nanopermporometry, can measure pores as small as ~0.6 nm[35,37-39]. The limitation of this method lies in the assumption that the pores are non-intersecting and cylindrical, which may differ from the actual pore size distribution. Attempts have been made to improve accuracy by introducing an invasion percolation model to refine the analysis[40,41].
2.3.2.2. Gas adsorption-desorption
Gas adsorption-desorption is a widely used technique for analyzing the surface area and pore size distribution of not only porous membranes but also general porous materials. It has been applied to various membrane materials, including ceramics and polymers [42-45]. The experimental steps using a surface area analyzer (BET method) are as follows: 1) The membrane sample is thoroughly dried, 2) Helium gas is introduced at cryogenic temperatures to measure the dead volume of each test and reference cell, after which the helium is evacuated, 3) Nitrogen gas is introduced at a constant flow rate, and the pressure is monitored over time to determine the adsorption-desorption isotherm of the sample[46].
Pores are classified into three categories: macropores (pore size > 50 nm), mesopores (pore size: 2~50 nm), and micropores (pore size < 2 nm). The behavior and interpretation of the adsorption-desorption isotherm depend on the pore type. Analyzing macropores or micropores in membranes using this method can lead to measurement errors. Specifically, micropores cannot be explained by the Kelvin equation because the meniscus inside these pores does not form continuously[47]. Therefore, the gas adsorption-desorption method is best suited for mesoporous membranes, such as those used in ultrafiltration. The hysteresis observed in the adsorption- desorption isotherm may cause discrepancies between the actual and measured pore size distributions, often due to intersecting pores within the membrane [48].
2.3.2.3. Evapoporometry
Evapoporometry involves filling the membrane pores with a highly volatile wetting fluid and then gradually evaporating it while measuring the membrane weight in real-time to determine the pore size distribution[49-51]. To ensure measurement precision, the external environment must be temperature and humidity controlled [49]. The wetting fluid evaporates from larger to smaller pores, and common fluids include acetone, ethanol, isopropyl alcohol, and butanol. The effectiveness and precision of this method depend on the chemical properties of the wetting fluid. The preferred wetting fluids are those that: 1) have high volatility (or high vapor pressure), and 2) possess high interfacial tension and molar volume[51].
Evapoporometry has been reported mainly for analyzing polymeric membranes such as PES, PC, and PVDF with pore size distributions of approximately 20~100 nm. This technique can analyze both open and closed pores simultaneously, as it varies the behavior of evaporation by filling closed pores with a non-volatile, non-wetting fluid (e.g., silicone oil). The difference in the rate of evaporation over time then correlates with changes in pore size distribution[49-52]. Therefore, the results of evapoporometry depend on whether non-volatile and non-wetting fluids are used, and the method should be adapted according to the membrane material and pore characteristics.
The advantages of evapoporometry include its non-destructive nature, the ability to analyze both open and closed pores, and its applicability to ultrafiltration and nanofiltration membranes (pore sizes 1~10 nm). However, limitations include: 1) The assumption of cylindrical pores, and 2) The lack of clear criteria for selecting the appropriate wetting and displacing fluids.
2.3.3. Gibbs-Thomson equation
While the Kelvin equation explains thermodynamic equilibrium under isothermal conditions through pressure variations, the Gibbs-Thomson equation describes thermodynamic equilibrium under isobaric conditions, focusing on temperature changes. When a liquid confined in pores solidifies at low temperatures, the free energy at the liquid-solid interface changes. This can be expressed as a function of phase transition temperature (T) as follows[53,54]:
Here, T0 is the initial temperature, VL is the molar volume of the liquid, ΔSF is the heat of fusion, γls is the liquid-solid interfacial tension, and Cls is the curvature of the liquid-solid interface. The solid formed from the liquid, whether as a meniscus or a spherical nucleus due to nucleation, has the following curvature:
where r is the pore diameter, and t is the thickness of the condensate layer that adheres to the pore wall and does not participate in the phase transition. Combining Equations 7 and 8 allows for determining pore size as a function of temperature variation. For instance, if benzene is used as the liquid, the equation becomes [54]:
2.3.3.1. Thermoporometry
Thermoporometry is a technique capable of analyzing pore sizes in the range of 3~300 nm. While it has not been extensively reported for membrane pore analysis, it has primarily been applied to polymer materials such as polysulfone (PSF) and polycarbonate (PC) [55-57]. The pore size distribution is plotted based on phase transition curves obtained using calorimetry. The sample is soaked in a liquid, frozen, and then gradually thawed while monitoring the change in heat. Benzene or water is typically used as the fluid because these solvents facilitate observation of freezing-melting conditions in calorimetry setups.
3. Conclusions
Membranes can be considered a type of black box, where the membrane medium is continuously influenced by the physicochemical properties of the permeating species, and, in turn, the diffusion behavior of the permeating species is dynamically affected by the membrane. In porous membranes, the diffusion pathways (pores) are formed irregularly within the membrane matrix, making it challenging to define boundary conditions. This paper reviews analytical techniques that can standardize and infer these irregular pore characteristics. The most commonly reported and preferred methods for determining pore size distribution are based on the Young-Laplace equation, Kelvin equation, and Gibbs-Thomson equation, with liquid displacement methods and evapoporometry being particularly prevalent. The features and limitations of each method are summarized in Table 1.
A common limitation of these techniques is that they assume the pores are cylindrical, which may not accurately represent the actual pore structure and characteristics. Additionally, none of these methods provide a definitive criterion for selecting the fluid used to fill the pores. Since the pore characteristics of membranes depend on the properties of the wetting fluid (e.g., wettability, calorific value, interfacial tension), it is important to choose the appropriate fluid and analytical method based on the membrane’s material, structure, and properties. Furthermore, the results of pore characteristic analysis should be cross-validated with other techniques such as SEM and AFM to ensure accuracy.