Open access peer-reviewed chapter

Dissolution Mass Transfer of Trapped Phase in Porous Media

By Anindityo Patmonoaji, Yingxue Hu, Chunwei Zhang and Tetsuya Suekane

Submitted: September 11th 2020Reviewed: December 10th 2020Published: January 25th 2021

DOI: 10.5772/intechopen.95448

Downloaded: 154


Dissolution mass transfer of trapped phase (TP) to flowing phase (FP) in porous media plays significant roles in hydrogeology, e.g., groundwater contamination by non-aqueous phase liquids, groundwater in-situ bioremediation, and geological carbon sequestration. In this chapter, this phenomenon is described. First, the physical and mathematical models are given. Afterwards, various conditions affecting this process, i.e., porous media characteristics, capillary trapping characteristics, flow bypassing, TP characteristics, and FP velocity, are discussed. These various conditions are described based on three parameters affecting the dissolution mass transfer: TP interfacial area (A), TP dissolution ratio (ξ), and mass transfer coefficient (k). Eventually, models to predict the mass transfer are formulated based on non-dimensional model. All of the data in this chapter are based on the experiments obtained by using micro-tomography and a series of image processing techniques from our latest works.


  • dissolution
  • mass transfer
  • porous media
  • micro-tomography
  • interfacial area
  • trapped phase cluster

1. Introduction

1.1 Applications of dissolution mass transfer in porous media

Dissolution mass transfer of trapped phase (TP) to flowing phase (FP) in porous media plays significant roles in hydrogeology (Figure 1), e.g., groundwater contamination by non-aqueous phase liquid (NAPL) [1, 2, 3], groundwater in-situ bioremediation [4, 5], and geological carbon sequestration (GCS) [6].

Figure 1.

Dissolution mass transfer in porous media phenomena in hydrogeology (modified from Patmonoaji et al. [7] with permission).

In groundwater contamination by NAPL, NAPL could leak from industrial plant and contaminate the groundwater stream [1, 2, 3]. When the NAPL enters the groundwater, it displaces the existing groundwater, and then it was trapped due to capillary trapping. As it is trapped, mass transfer from the trapped NAPL to the groundwater flow occurs, transferring the dissolved NAPL to the groundwater. As a result, groundwater is contaminated, and the period of contamination depends on the mass transfer rate of trapped NAPL to the groundwater. Low mass transfer rate results in a prolong contamination, and vice versa.

In-situ bioremediation is one of the decontamination method of NAPL contamination [4, 5]. This method harnesses the ability of native bacteria in the groundwater to biodegrade the contaminants. To enhance the biodegradation process, O2, CH4, NH3, CO3 and N2 gases are injected at the upstream of the contamination area to supply nutrients and stimulate the bacteria. After the injection, these gases are trapped due to capillary trapping and thus mass transfer occurs. The mass transfer rate is important because it affects the efficiency of the gas transport. Higher mass transfer rate means faster delivery but faster depletion of trapped gases. Therefore, it should be controlled to match the bacteria needs.

In GCS, the captured CO2 from the point source, e.g., coal power plant, steel factory, and cement factory, is injected into deep saline aquifers for sequestration. After the injection, the CO2 is trapped in the groundwater due to capillary trapping [6, 8]. Although the CO2 can be sequestered this way, the reservoir pressure increases due to the injection. However, by the time, the trapped CO2 will dissolve, reducing the reservoir pressure. As a result, the mass transfer rate is critical to approximate the pressure reduction rate for safer CO2 sequestration process.

Given the scale of the phenomena is massive, the observation scale can be categorized into pore-scale, core-scale, and field-scale (Figure 1). Pore-scale is in the order of μm, and core-scale is in the order of cm. Field-scale, however, is in the order of m to km. Each of the scale is used to explain the scales above it. Pore-scale and core-scale are performed in a laboratory, whereas field-scale is performed in the field. Only pore-scale and core-scale will be discussed in this chapter.

Given the opaqueness of porous media system, direct observation of dissolution process is difficult to be performed. To evaluate the mass transfer process in detail, the interfacial area (A) of TP needs to be measured. With the information of interfacial area of TP, the mass transfer coefficient (k) can be calculated (described later in Eq. (5)). In addition, the distribution of the TP clusters is also need to be monitored to assess the dissolution mass transfer process.

Earlier studies mainly relied only on core-scale studies by the measurement of dissolved phase concentration from the effluent and amount of trapped phase through gravimetrical measurement without knowing the trapped phase characteristics inside [1, 2, 3, 4, 5]. However, because the trapped phase characteristics are required to measure the Afor the mass transfer, and its characteristics are strongly controlled by the porous media characteristics, the studies become phenomenological because kcannot be calculated. In addition, the overall dissolution process cannot be observed, and thus flow bypassing cannot be monitored although it also affects the mass transfer process.

In this chapter, dissolution mass transfer in porous media is discussed based on experimental observation on the dissolution mass transfer of various trapped gas into flowing water in porous media by using micro-tomography and image processing [7, 9, 10, 11, 12]. By using this methods, Aand distribution of TP clusters can be monitored throughout dissolution process, and thus kcan be calculated. In addition, by monitoring dissolution process, flow bypassing was investigated, and thus its effect to kwas found. Moreover, unique dissolution process of TP with high dissolution ratio (ξ), which is the ratio of solubility limit (Csol,i) and density (ρi) was observed. Its effect to kwas also found. Therefore, all of three parameters effecting dissolution mass transfer (A,ξ, and k), were elucidated.

The contents of this chapter are described as follow. In Section 2, the physical and mathematical model are given. In Section 3, the general description about the micro-tomography and image processing methods are described. In Section 4, various parameters affecting this process, i.e., porous media characteristics, capillary trapping characteristics, flow bypassing, TP properties, and FP velocity, are discussed. All of these condition will be investigated based on the point of view of three parameters affecting dissolution mass transfer: A,ξ, and k. Porous media together with capillary trapping characteristics affect the amount ofAfor mass transfer process. Flow bypassing affects the uniformity of mass transfer process and eventually k. TP properties affect the ξthrough its Csol,iand ρi. Lastly, FP velocity affects the advection rate of the solute, influencing k. Eventually, the mass transfer model based on non-dimensional number is also formulated at the end of the chapter.


2. Physical and mathematical model

2.1 What is dissolution mass transfer?

At first, we can start from the basic by describing mass transfer. Mass transfer is a process of a movement of mass from one medium to another medium. The examples can be found easily in the surrounding of our daily lives. For example, the mass transfer of camphor in a bathroom into the surrounding air to improve the smells of the room (Figure 2a). Another example is the mass transfer of oxygen from the air into a lake (Figure 2b). Without the supply of oxygen through mass transfer, the fishes living in the lake cannot breathe the oxygen and stay alive.

Figure 2.

(a) Mass transfer of camphor to surrounding air, (b) mass transfer of oxygen from air to water in lake, and (c) illustration of mass transfer process.

Similar with other transport phenomena process, i.e., heat transfer and momentum transfer, mass transfer occurs because of non-equilibrium condition [13, 14, 15]. In case of mass transfer, equilibrium can be achieved once the system possesses equal chemical potential (η) of the considered component (i):


with subscripts 1 and 2 represent the phases. This equation can also be described with partition coefficient or Henry’s coefficient (Figure 2c) as follow:


with Cias concentration and Hiccas dimensionless Henry coefficient of component i, which is constant with pressure but varies with temperature.

The concept of mass transfer nowadays is mainly based on the works of Thomas Graham (1805–1869) and Adolf Fick (1829–1901). Thomas Graham was the pioneer in the study of mass transfer through his works in the diffusion of H2 gas in a tube [16]. By using analogies with Fourier’s work [17] in heat conduction, Adolf Fick formulated the based equation for mass transfer that is the Fick’s law of diffusion mass transfer [18] defined as follow:


with jas mass flux and Dias diffusion mass transfer coefficient of compound iin the medium. Another form is based on Newton’s law of cooling as follow:


with kas mass transfer coefficient. In this approach, kand ΔCiis used instead of Diand Ci. As a result, the mass transfer is judged based on concentration difference without knowing the concentration gradient. Because the concentration gradient in this system is difficult to obtain in the system, this form of mass transfer equation is more preferable in mass transfer from an interface into bulk solution, such as in the hydrogeology processes mentioned in Section 1.

In case of dissolution mass transfer, an additional process present. Dissolution here refers to the decrease in volume due to mass transfer. As a result, the Eq. (4) can be modified to volume perspective as follow:


with ρias density of compound i, Vas volume, and tas time. Therefore, high mass transfer rate also results in high volume decrease rate or dissolution rate. The example of camphor in Figure 2a provides an example of dissolution mass transfer, whereas oxygen mass transfer in a pool of water in Figure 2b provides an example of mass transfer without dissolution.

2.2 Physical model of mass transfer from an interface into a bulk solution

There are several mechanisms models explaining mass transfer from an interface into a bulk solution [14]. Due to the analogy between mass transfer and heat transfer, similar models with heat transfer are also used. However, because most of heat transfer application is in fluid–solid system, whereas mass transfer applications occur in fluid–fluid system, not all model can be used for both.

The simplest physical model is the stagnant film theory (Figure 3a) [19] that assumes the existence of a stagnant film near the interface of mass transfer process. This stagnant film is also referred as an unstirred layer because no fluid motion exists in this film. However, this film is hypothetical because fluid motions still occur beyond the no-slip condition at interface. The effects of velocity and other parameters are then represented by the thickness of this stagnant film. Higher mass transfer rate results in thinner δ, and vice versa.

Figure 3.

(a) Illustration of stagnant film model, (b) penetration theory model at low and high velocity, and (c) surface renewal theory model.

Second mechanisms model is penetration theory (Figure 3b) [20]. The key assumption is that a thick moving film is continuously generated in the bulk fluid, which the mass transfer occurs by diffusion across this film. The effect of velocity to the mass transfer is represented by the concentration distribution. As the solute diffuse into the bulk liquid, it will be flushed away by the flow. Higher velocity results in faster removal of the solute. As a result, higher velocity results in steeper concentration distribution, leading to higher mass transfer rate, and vice versa.

Third mechanisms model is surface renewal theory (Figure 3c) [21], which is the further development of penetration theory. The bulk liquid is depicted into two regions, a region near the interface that is renewed quickly and a region of well-mixed bulk fluid. The renewed region is renewed by the flow and similar with penetration theory. However, it is constantly exchanged with new elements from a second bulk region, making this theory closer to the realistic physical situation.

Among these models, however, no theory is completely satisfactory. Nevertheless, these models provide valuable insight to the mechanisms in mass transfer. For more detail, readers are advised to check the articles in the references [19, 20, 21].

For more empirical modeling of the mass transfer process, an empirical equation based on dimensionless number is usually used. The experimental data of mass transfer is represented by power-law relation of dimensionless number as follow:


with ϵ1, ϵ2, and ϵ3as constants, Shas Sherwood number representing the mass transfer rate and Reas Reynolds number representing the bulk fluid velocity and geometry of the system. These non-dimensional number are defined as follow:


with ρFPas FP density, μFPas FP absolute viscosity, Uas FP velocity, and Las characteristic length.

In case of porous media, the physical mechanism of mass transfer is actually similar. In an open space, when dissolution occurs, the solute can be flushed away easily. In porous media, on the other hand, these processes are governed by the solute transport in porous media. In addition, in the application of hydrogeology, the TP clusters are distributed throughout the porous media (Figure 4), resulting in continuous mass transfer process along the porous media. As a result, a special attention to the solute transport throughout the porous media is required before assessing the mass transfer process with Eq. (5).

Figure 4.

Illustration of dissolution mass transfer in porous media and the predictedCi,x,tbased onequation 16.

2.3 Fluid flow and solute transport in porous media

The fluid flow in hydrogeology application is relatively slow and can be classified as Stokes flow or creeping flow, which the ratio between fluid velocity and viscosities is relatively low. This condition is represented by Relower than one. Steady-state Stokes equation is given as follow:


with Pas pressure and gas gravitational acceleration. This equation is then further derived into Darcy law, which is introduced by Henry Darcy [22] as a relationship to describe fluid flux (q) in sand filters, which was first proposed by Nutting [23] and Wyckoff et al. [24], described as follow:


with Kas permeability. Given its original development, Darcy equation has been used mainly in the field of fluid flow in porous media. Similarly, with Stokes flow, to satisfy Darcy regime, the Reneeds to be low. Experimental tests demonstrated that it can still be achieved when the Reis lower than ten [25].

For the solute transport in porous media, it is affected by the tortuosity and velocity profile in the pore networks. Nevertheless, under Darcy regime, it can be predicted by using advection-dispersion equation [26] as follow:


with, Dhdlas hydrodynamic dispersion coefficient and Ras external mass transfer source or sink. With this equation, solute transport along the porous media can be predicted, and it can be used to assess the mass transfer with Eq. (5).

2.4 Calculation of dissolution mass transfer in porous media

To evaluate the dissolution mass transfer in porous media, combination of modified version of Eqs. (5) and (11) are used [1, 2, 3, 4, 5, 9]. To calculate the mass transfer coefficient in porous media, Eq. (5) is slightly modified as follow:


with ϕas porosity, STPas TP saturation, and aTPas specific interfacial area, which is the amount of Aper porous media volume. However, the Ci, which is the driving force of the mass transfer, is unknown. To approximate it, Cican be calculated using Eq. (11). Because the experiment can be assumed as a one-dimensional column, one-dimensional advection-dispersion equation as follow:


with STP,csas cross-sectional average of TP saturation and xas axial direction, can be used. In addition, it can be simplified with some assumptions [1, 2, 3, 9]. First, because dissolution mass transfer is a prolong process, the problem can be approached as a pseudo-steady state, neglecting the time accumulation term at left-hand side of the equation. Second, the dispersion term, which is the first term at the right-hand side of the equation, can be neglected due to the dominance in advection. It can be assessed by calculating the Peclet number (Pe) as follow:


with dpar,50%as median particle size. Perepresents the ratio of advection and diffusion of solute transport. High Penumber represents advection dominant regime, whereas low Perepresents diffusion dominant regime. When the value of Peis higher than five, it represent advection dominant regime [27]. Thus, the dispersion term can be neglected. This was also proven with order of magnitude analysis [3, 9]. The order of magnitudes of accumulation and dispersion terms are significantly lower than the others. Hence, it can be simplified to:


By discretizing ∂xto Δx, Ciat the location of x+Δxand at the given time of t(Ci,x+Δx,t) as shown in Figure 4, can be calculated as follow:


The average of concentration in between xand x+Δxat the time of t(Ci,xavg,t) (Figure 4) then can be calculated by averaging the Ci,x,tand Ci,x+Δx,tas follow:


This Ci,xavg,tthen can be used to calculate the mass transfer coefficient as follow:


with aTP,cs,tavgas the cross-sectional average of aTPthat is averaged in between the time of tt(aTP,cs,tt)and t(aTP,cs,t)calculated as follow:


For the modeling of mass transfer, Reand Shin Eqs. (7) and (8) were used with Lequal to dpar,50%.

3. Experimental procedure

3.1 Experimental setup

As mentioned in Section 1, this chapter is based on the experimental observation by using micro-tomography and image processing [7, 9, 10, 11, 12]. For the micro-tomography, an X-ray microtomography scanner (micro-CT) was used. This micro-CT is able to generate 3D images in the form of 992 cross section slice image with the dimension of 992 × 992 voxels with the voxel size of 16.427 μm. As a result, observation area with the size of 16.34 × 16.34 × 16.34 mm3 was generated.

To generate the porous media, cylindrical container with an inner diameter of 10 mm was used as the casing of porous media. Inside this container, various particles will be packed inside to generate various porous media. All of the packing was performed under vibration to generate closed random packing classification [28]. As a result, the quality of the packing can be controlled. Additionally, sintered glass discs were placed at inlet and outlet to generate uniform flow.

To generate capillary trapping and similar condition with hydrogeology conditions in Section 1, the porous media will be sequentially saturated with degassed deionized (DDI) water and the gas. First, the gas was injected into the porous medium to saturate the porous media and remove all of other gases. Second, DDI water was injected upward to saturate the porous media, generating a condition similar with groundwater. Third, the experimental gas was injected again to remove most of the DDI water, resulting in similar condition with gas injection into groundwater. Finally, the DDI water was injected upward at the desired experimental flow rate. With this method, the generated trapped gas is uniform along the sample. Once the DDI water had breakthrough, and no gas flow was observed in the outlet tube, creating a quasi-steady state, the experiment began, and the dissolution process was monitored with micro-CT. The schematic about the experimental setup is given in Figure 5.

Figure 5.

Schematics of experimental setup with micro-tomography (modified from Patmonoaji et al. [9] with permission).

3.2 Image processing

For the data processing, image processing techniques were carried out to quantify the data from the micro-CT image. The main method is to differentiate the gas from the surroundings by performing image binarization (Figure 6). Although the contrast between water and particles is indistinguishable, the contrast of the gas among other phases is easily distinguishable. Therefore, a direct thresholding method was enough to separate gas from other phases. With this binarized images, then quantitative data, i.e., amount of trapped gas, and surface area of the trapped gas, can be simply derived with image processing to assess the dissolution mass transfer. Furthermore, by identifying separated trapped gas clusters, each of the trapped gas clusters can be characterized based on its volume and surface area. In addition, by performing similar method to porous media saturated with gas, porosity of the porous media can be obtained. Moreover, by applying watershed (WS) segmentation method [11], the throats and pores of the porous media can be identified and measured.

Figure 6.

Binarization image processing methods: (a) cross sectional image of the porous media, (b) histogram of the gray value, (c) binarized image with black as the gas, and (d) 3D binarized image of gas (modified from Patmonoaji et al. [1] with permission).

4. Factors affecting dissolution mass transfer in porous media

4.1 Porous media characteristics

A porous medium can be described as a solid containing void networks inside it. Structure-wise, the void network can be seen as a series of pores and throats (Figure 7a) [28, 29]. Pores are connected by a throat, which is a constriction in the network. Because capillary trapping mainly occupies pores, pore characteristics affect the morphology of the capillary trapping that eventually controls the amount Afor the mass transfer process (Figure 7b).

Figure 7.

(a) Illustration of pore and throat in pore network and (b) effect of porous media characteristics to TP morphology.

These pore and throat in porous media are mainly governed by the shape, size, and distribution of the granular particles forming the porous media. Shown in Table 1 are the properties of nine granular particles used to generate porous media [10, 11]. Those properties Cuas particle distribution uniformity coefficient, Ccas particle distribution gradation coefficient, and Φas porosity. Three particle shapes, which are plastic beads (PB), glass beads (GB), and Toyoura sands (TS), are used. PB represents angular shape particle like quartz sands, whereas GB represents sphere particles. TS, on the other hand, represents natural silica sands.

TypeSize range (μm)dpar,50%(μm)CuCcΦ

Table 1.

Properties of nine porous media generated from different shape, size, and distribution of granular particles.

Figure 8a shows the pore size distribution (PSD) of the porous media generated from those particles. The PSD is shown as probability density function (PDF) of pore equivalent diameter (dpore). To obtain this PSD, combinations of micro-tomography and WS method were performed [11]. GB tends to generate larger dporethan PB, whereas TS generate similar dporewith GB. For particle size, larger particles generate larger dpore. However, in case of particle size distribution as compared by PB 180–212 μm, PB 125–300 μm, and PB 125–710 μm, similar PSD was found. Although the range of particle size distribution is different, the dpar,50%shown in Table 1, is actually similar. Pores are voids generated due to contact of granular particles (Figure 7b). Contact of large particles generates large voids. However, if the distribution of particles is wide. Large voids that is generated by large particles can be filled with the smaller particles. As a result, dporetends to be controlled by the smaller particle size, particularly dpar,50%. This is also consistent with the Φ. Wider particle range generates lower Φbecause the large pores generated from the contact of large particles are filled by smaller particles.

Figure 8.

(a) PSD and (b)dpar,50%-normalized PSD.

Proportionality of the PSD is also found by normalizing PSD with the dpar,50%(Figure 8b). By normalizing the PSD with dpar,50%, the PSD tend to converge into a single curve but not exactly. The differences are due to particle shape variation. Nevertheless, it demonstrates proportionality to dpar,50%for each particle shape.

The pore shape can also be evaluated through its surface area (Apore) and volume (Vpore) as shown in Figure 9. Each point is the result obtained from averaging Aporeof all pores in a log-scale range of Vpore. Typical pores shape detected in porous media with PB, GB, and TS is also given as yellow, blue, and green objects, respectively. As shown in normal-scale, pores in smaller dpore,50%generates larger Apore. Related with particle shape, the Aporeof PB tends to be higher than GB and TS. GB generates pore with the smallest Apore, whereas TS generate pore with Aporein between PB and GB. This higher Aporeis due to the angular shapes of PB.

Figure 9.

VporeandAporein the porous media and the typical pore shapes in each particle shape (modified from Patmonoaji et al. [11]).

4.2 Capillary trapping characteristics

Capillary, buoyancy, and viscous forces play important roles in the formation of capillary trapping (Figure 10a). These forces can be represented with capillary number (Ca) and Bond number (Bo)as follow:

Figure 10.

(a) Illustration of forces exerting on TP in porous media and morphology of (b) trapped NWP and (c) WP in porous media.


Carepresents the comparison between viscous and capillary forces, whereas Borepresents the comparison between buoyancy and capillary forces. Buoyancy and viscous forces exert on the TP to escape from the trapped location, whereas capillary force resist these forces to keep the TP in the pore. If Boand Cais higher than 0.35 and 10−3, respectively, no capillary trapping occurs due to excessive buoyancy and viscous forces [30, 31]. In hydrogeology applications, the Boand Caare much lower, resulting in the possibility of capillary trapping generation.

Wettability of the porous media also affects the TP clusters morphology (Figure 10b and c). In case the TP is non-wetting phase (NWP), the TP is located in the center of the pore. On the other hand, in case the TP is wetting phase (WP), the TP tends to surround the solid. In this work, we mainly focus on the former condition.

The trapped phase size distribution (TPSD) in the porous media described in Table 1 is given in Figure 11 as PDF of TP equivalent diameter (dTP). As shown in Figure 11a, with the increase in dpar,50%, dTPincreases. Compared with the PSD, the TPSD tends to be larger due to the possibility to form multi-pores bubbles. By normalizing with dpar,50%as shown in Figure 11b, the TPSD also converge demonstrating proportionality to dpar,50%. This proportionality of TPSD with dpar,50%can be traced back to the proportionality of PSD to dpar,50%in Section 4.1. Because capillary trapping mainly occupies pore, the proportionality of PSD to dpar,50%also represents the proportionality of TPSD to dpar,50%.

Figure 11.

(a) TPSD and (b)dpar,50%-normalized TPSD (modified from Patmonoaji et al. [10]).

The TP cluster morphology can also be characterized through its volume (VTP,cltr) and surface area (ATP,cltr) as shown in Figure 12. Each point is the result obtained from averaging the ATP,cltrof all TP clusters in a log-scale range of VTP,cltr. The examples of the TP clusters shape are also given in Figure 12. At smaller size, the markers from all porous media coincides one another. However, as the VTP,cltrincrease, they branched out with TP clusters from smaller particle size exhibit larger ATP,cltr. The point of the branch out is approximately at the size of dpore,50%. When the TP cluster size is smaller than dporeof the porous media, it tends to occupy only one pore with a sphere-like shape. However, for larger cluster, it has to occupy more than one pore. As a result, it extends to the neighboring pores forming a channel of multi-pores TP cluster. These single-pore and multi-pores clusters can be modeled with a sphere and an elongated capsule model as follow:

Figure 12.

The averageATP,cltrat the givenVTP,cltrand examples TP clusters (Patmonoaji et al. [10]).


For the aTP, it is given in Figure 13a. The line represents the cross-sectional average aTP(aTP,cs), whereas the marker represents the overall average aTP(aTP,ovr). Figure 13a shows that smaller particles generate larger aTP. As described previously, smaller particles generate smaller PSD and TPSD. The single-pore TP clusters are sphere with smaller diameter, and the multi-pores TP clusters are elongated capsule with smaller diameter. As a result, the TP clusters generate higher aTP. For the effect of particle size distribution range, similarly with PSD and TPSD, it is mainly controlled by the dpar,50%. For the particle shape, PB generates the highest aTP, whereas GB and TS generate similar aTP. Given the PSD and TPSD, PB generates smaller pore and TP clusters, and thus higher aTPwas produced. GB and TS generate similar PSD but smaller than PB. Therefore, similar aTPwas produced between GB and TS but lower than PB. In addition, as shown in Figure 9, pores of PB is more irregular, and thus it generates TP clusters with more surface area. On the other hand, the pores of GB and TS are much less irregular than PB, and look relatively similar. Therefore, based on these two reasons, PB generate higher aTP, whereas GB and TS generate similar aTPbut lower than PB. In addition, by normalizing aTPwith dpar,50%, the aTPconverges into two lines (Figure 13b). All aTPof PB converge into a single line, whereas TS and GB tend to converge into another line. Hence, two models were generated:

Figure 13.

(a)aTPin all porous media and (b)dpar,50%normalized-aTPwith lines and markers representaTP,csandaTP,ovr, respectively (taken with permission from Patmonoaji et al. [10]).


For the linearity between aTPand STP, it can be explained with these relations:


with ATPas trapped phase interfacial area, VTPas trapped phase volume, and nTP,cltras the number of TP clusters. When the trapped gas is mainly multi-pores TP clusters, which the function of ATP,cltris linear with VTP,cltr(Eq. (23)), linear function between aTPand STPwill be generated. In contrast, if it consists of large number of single TP cluster, the relation will be controlled by nTP,cltr. Because nTP,cltris a function of VTP, the relation of aTPand STPwill be linear also.

As described in this section, the aTPgenerated in porous media can be predicted by using Eqs. (24) and (25), which is based on STPand dpar,50%. Given these two parameters only, the aTPcan be predicted. Higher STPgenerates higher aTP. Smaller dpar,50%also generates higher aTP. Porous media generated from more angular particles will also generates higher aTP. As described in Eq. (12), higher aTPresults in faster dissolution mass transfer. Therefore, porous media generated from angular and smaller particles leads to faster dissolution, especially at the condition of high STP.

4.3 Flow bypassing

In general, flow bypassing can be divided into two categories. First is bypassing due to heterogeneity of porous media [32, 33], and second is bypassing induced by dissolution fingering in homogeneous porous media [34, 35, 36]. Both of them direct water flow to preferential flow path, resulting in non-uniform dissolution mass transfer. As a result, it will affect the k. In this section, bypassing induced by dissolution fingering is examined through dissolution experiment performed in porous media from PB with various sizes and distributions [7] shown in Table 2. For the TP, N2 gas was used.

Size range (μm)125–150180–212355–425600–710125–300250–425125–710

Table 2.

Characteristics of particle, porous medium, flow, and mass transfer inside each porous medium sample.

In general, the mechanism of dissolution fingering development can be explained as follow. Initially, capillary trapping is relatively homogeneous, and thus the permeability tends to be homogeneous as well. However, a slight difference in permeability remains present. This slight difference generates preferential water flow, directing the water to flow to region with higher permeability. At this higher permeability region, water flow slightly faster, and thus the TP clusters dissolve slightly faster as well. As a consequence, it gradually increases the permeability of this slightly permeable region, enhancing the permeability differences gradually. This results in a positive feedback, and thus the region of this preferential path grows with the dissolution time, resulting in the development of dissolution finger.

In short, the possibility of dissolution fingering development can be evaluated by measuring the approximated width of dissolution finger (λ) and compare it with the given size of porous medium diameter sample (Ly), which in this work is 10 mm. By using similar analysis for mineral dissolution instabilities based on linear stability analysis [37, 38], we can approximate the size of the dissolution finger by measuring the dimensionless growth rate (ω¯) using the calculation from Zhao et al. [36] as follow:


with m¯as dimensionless wave number, Swias irreducible water saturation, SN2,0as initial gas saturation, τas tortuosity, αLas longitudinal dispersivity, αTas transversal dispersivity. By calculating the ω¯at the given value of m¯, the graph of ω¯at the given m¯can be generated as shown in Figure 14. At the value of ω¯lower than zero, the finger will stabilize, whereas at the value ω¯larger than zero, the finger will develop. The peak value of ω¯represents the value of critical dimensionless wave number (m¯critical), which the wave will grow fastest and most likely to be observed. Therefore, by finding the m¯criticalin each experiment cases, λin each experiment cases can be calculated as follow:

Figure 14.

The function ofω¯andm¯for the prediction of dissolution finger width calculated for PB 600–710 μm.


The dissolution process inside PB 600–710 μm and PB 180–212 μm are given in Figure 15a as 3D images. In PB 600–710 μm, dissolution process occurs uniformly, whereas in PB 180–212 μm, dissolution occurs non-uniformly. The kthroughout the dissolution process is given in Figure 15b with the function toSTP,cs. Two groups of kare found. First group is higher k, that occurs in PB 125–710 μm, PB 250–425 μm, PB 355–425 μm, and PB 600–710 μm. Second group is lower kthat occur in PB 125–300 μm, PB 125–150 μm, and PB 180–212 μm. From the predicted λ, the first group belong to the value of λlarger than Ly, whereas the second group belong to the value of λsmaller than Ly. When the size of λis smaller than Ly, dissolution fingers could fully develop, resulting in bypassing and lower k. This trend occurs except in PB 125–710 μm that was influenced by the particle size distribution. The classification of kand predicted λis given in Table 2.

Figure 15.

(a) Dissolution progress of (top) particle 600-710 μm and (bottom) particle 180-212 μm and (b)kat the given saturation in all porous media (modified from Patmonoaji et al. [7]).

As explained in this section, in the presence of bypassing induced by dissolution fingering, kis lower due to non-uniform dissolution process. Because this is also occur in homogeneous porous media, and the limiting condition is just the ratio of the container to the predicted size of dissolution fingering, it is expected to occurs in field-scale of hydrogeology processes.

4.4 Trapped phase properties

Another parameter that affect the dissolution mass transfer rate is ξ, which is the ratio of Ci,soland ρi. ρiaffects the decrease rate in volume of dissolution process, and Ci,solserves as the driving force of the mass transfer process. Higher ξresults in faster dissolution mass transfer process. As comparison, ξvaries for about 10−9–10−3 in mineral dissolution, 10−4–10−3 in NAPL dissolution, and 10−3–10−1 in gas dissolution [36]. Timewise, this implicates that the dissolution of minerals could take months and years, whereas dissolution of NAPL could take days to weeks. In contrast, dissolution of gases could only take minutes to hours.

Dissolution in porous media is not only mass transfer process, but also solute transport. When mass transfer occurs from TP, the dissolve solute need to be transported away from the surrounding TP by the FP. Therefore, high value in ξcould generate problems in the balance of mass transfer and solute transport.

To compare the effect ξ, CO2 and N2 gases were used. The ρiof N2 and CO2 gases are 1.842 and 1.165 kg/m3, whereas the Ci,solare 1.449 and 0.0175 kg/m3, resulting in ξof 0.787 and 0.015, respectively. The differences in ξis more than fifty times. The experiments were performed with the same porous media and FP velocity, resulting in the same Re.

The axial average TP saturation throughout the dissolution process is given in Figure 16. The dissolution of N2, and CO2 gases behave differently. In case of N2 gas, as the time proceeds, STP,csdecrease from the initial value to zero. For CO2, on the other hand, the STP,csdoes not dissolve completely but stop at the value about 0.03 before it dissolves completely, resulting in two dissolution stages. First dissolution stage occurs as the dissolution of STP,csto about 0.03. Second dissolution stage occurs as the complete dissolution of the remaining STP,cs. By varying Re, the same dissolution pattern of CO2 is still observed [9].

Figure 16.

STP,csof N2 and CO2 gas during dissolution process under Re of 0.04.

By observing the dissolution process from the 3D images of segmented TP clusters as shown in Figure 17, the different dissolution behavior can be observed in more detail. In case of N2, TP clusters dissolve completely from multi-pores size. However, in case of CO2, at first, the multi-pores TP clusters dissolve to single-pore TP clusters. Afterwards, second stage occurs, dissolving the remaining clusters.

Figure 17.

TP cluster of N2 and CO2 observed at different dissolution time. Different colors indicate different clusters (modified from Patmonoaji et al. [9]).

For the k, additionally, there is a significant different between CO2 and N2 gas as shown in Figure 18. Given the same Re, the kof CO2 is more than one order of magnitude lower than N2 gas. In addition, kshows a dipping point in around 0.03, which is consistent with the two stages dissolution of CO2. When the STP,csreaches 0.03, the dissolution process slowed down, and thus kdecreases. Then second dissolution stages occurs, dissolving the rest STP,cs, resulting in the increase of k.

Figure 18.

Averagekof all experiments with different gases andRe.

Compared with other gases, CO2 exhibit about fifty times larger. When TP cluster dissolves, the solute should be transported away by fresh incoming FP. If the solute is not transported, the solute concentration surrounding the TP cluster accumulates, slowing down the mass transfer. During the dissolution of CO2, this rapid dissolution is believed to generate area with high CO2 solute concentration due to much faster mass transfer rate than advection rate. In addition, as the trapped gas dissolves rapidly, it quickly retracts to the center of pores, increasing the difficulty of the solute to be removed by FP. As the dissolution occurs rapidly, the incoming FP cannot compete with the rapid dissolution of TP. As a result, two dissolution stages occur. First stage occurs before the accumulation of solute, whereas the second stage occurs after the accumulated solute is advected away from the surrounding of TP.

As discussed in this section, although higher ξresults in higher dissolution mass transfer rate, when it becomes too high, it could result in lower k. This condition can be represented by two dissolution stages that occurs due to competition between mass transfer rate and advection rate.

4.5 Flowing phase velocity

To investigate the effect of FP velocity, the experiments were also performed under different FP velocity (Rebetween 0.0016 and 0.16) with N2 and CO2 gases. They were performed with porous media from PB 250–425 μm to avoid bypassing and focus the discussion to the effect of FP velocity. By using the Shand Re, the mass transfer can be modeled with Eq. (6) as follow (Figure 19):

Figure 19.

Comparison of predicted-Shand measured-Shfrom the model.


The Eq. (40) is for N2 gas, whereas Eq. (41) is for CO2 gas. The difference is due to the effect of high ξin CO2, resulting in two dissolution stages. Because the kis also lower, the non-dimensional number model result in lower Shthan N2 gas.

In this section, the effect of FP velocity was shown. With the increase in FP velocity, kincrease. This relation can also be represented in non-dimensional number, as represented by Reand Shin Eqs. (40) and (41). The effect of high ξis also observable as lower Shdue to lower k.

5. Conclusion

Throughout this work, the physical mechanism and mathematical model of dissolution mass transfer of trapped phase in porous media was described. Various condition affecting it were also discussed from experiment data obtained from micro-tomography and a series of image processing techniques. The effects of these conditions were represented in three parameters: A,ξ, and k.

In case of porous media characteristics, the characteristics of porous media generated from various particle shape, size, and distribution were investigated. Smaller dpar,50%results in higher aTP. Porous media generated from more angular particles will also generates higher aTP. In addition, STPlinearly effect the amount of aTP. Therefore, porous media generated from angular and smaller particles leads to faster dissolution, especially at the condition of high STP.

Bypassing induced by dissolution fingering was also described. When the predicted λis smaller than Ly, dissolution finger can fully develop, resulting in bypassing and decrease in k. Because this occurs in homogeneous porous media, it is expected to occurs in field-scale of hydrogeology processes.

In addition, the effect of ξwas also investigated. If the ξof the TP is high, different dissolution process, which is represented by two stage dissolutions, occurs, leading to decrease in k. The competition between mass transfer and advection rate is believed to be the reason of this different dissolution process.

Eventually, the effect of FP velocity based on Shand Rewas formulated. With the increase of FP velocity, kincrease as represented by the increase in Shat the given Re. In case of CO2, the Shalso increase with Re, but it is still lower than N2. This is due to the effect two stage dissolution observed in CO2, resulting in lower k.



This work was supported by JSPS KAKENHI with grant numbers 17H00790 and 20 J14975.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Anindityo Patmonoaji, Yingxue Hu, Chunwei Zhang and Tetsuya Suekane (January 25th 2021). Dissolution Mass Transfer of Trapped Phase in Porous Media, Porous Fluids - Advances in Fluid Flow and Transport Phenomena in Porous Media, Vallampati Ramachandra Prasad, IntechOpen, DOI: 10.5772/intechopen.95448. Available from:

chapter statistics

154total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Heat Transfer Enhancement Using Unidirectional Porous Media under High Heat Flux Conditions

By Kazuhisa Yuki

Related Book

First chapter

3D Boundary Layer Theory

By Vladimir Shalaev

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us