Next Article in Journal
Laboratory Observations of Repeated Interactions between Ruptures and the Fault Bend Prior to the Overall Stick-Slip Instability Based on a Digital Image Correlation Method
Previous Article in Journal
A Comprehensive Feature Comparison Study of Open-Source Container Orchestration Frameworks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Particle Size Distribution of Colloidal Particles on Contaminant Transport in Porous Media

1
Department of Civil and Environmental Engineering, University of Ulsan, Daehak-ro 93, Nam-gu, Ulsan 680-749, Korea
2
School of Civil, Environmental and Architectural Engineering, Anam-dong 5-ga, Korea University, Seoul 136-713, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(5), 932; https://doi.org/10.3390/app9050932
Submission received: 25 December 2018 / Revised: 24 February 2019 / Accepted: 27 February 2019 / Published: 5 March 2019
(This article belongs to the Section Environmental Sciences)

Abstract

:
The presence of retained colloidal particles causes the retardation of contaminant transport when the contaminant is favorably adsorbed to colloidal particles. Although the particle size distribution affects the retention behavior of colloidal particles, the impact of particle size distribution on contaminant transport has not been reported to date. This study investigates the impact of the particle size distribution of the colloidal particles on contaminant transport through numerical simulation by representing the particle size distribution as a lognormal distribution function. In addition, the bed efficiency and contaminant saturation of simulated breakthrough curves were calculated, and a contaminant transport model with the Langmuir isotherm for the reaction between the contaminant–sand and contaminant–colloidal particle was introduced and validated with experimental data. The simulated breakthrough curves, bed efficiency, and contaminant saturation indicated that an increase in the mean and standard deviation of the particle size distribution causes the retardation of contaminant transport.

1. Introduction

It is now well-known that the transport and retention behavior of organic and inorganic colloidal particles (<10 µm) in porous media facilitates or retards contaminant transport when the contaminant is favorably adsorbed to colloidal particles [1,2,3]. Recent field and laboratory studies demonstrated that even a small amount of colloidal particles lead to the facilitation and retardation of contaminants such as radionuclides (e.g., 137Cs and 239Pu) or heavy metals (Pb(II) and Ni(II)) [4,5,6]. Therefore, the presence of colloidal particles as well as the transport behavior of colloidal particles should be taken into account in order to properly predict the contaminant transport. This would minimize underestimation or overestimation of the travel distance of contaminants [7,8].
Relatively small sizes of colloidal particles and large hydrodynamic forces exerted upon them cause the particles to be transported rather than retained, leading to the facilitation of contaminant transport. In contrast, if a relatively large number of colloidal particles are retained and not mobile, contaminants are more likely to be retarded and transport more slowly compared to the pore fluid. In this case, the larger amount of retained colloidal particles leads to the slower transport of contaminant because of more available adsorption sites in colloidal particles–sand media [9,10]. Therefore, an accurate estimation of the quantity of immobile colloidal particles cannot be overemphasized for the prediction of contaminant transport.
In order to predict the amount of immobile colloidal particles, many preceding studies of colloid transport have documented the impact of size ratio between colloidal particles and sand, solution chemistry (ionic strength or pH), fluid drag force, and degree of saturation [11,12,13,14,15]. However, the monodisperse latex colloids used in most of these experimental studies cannot represent the polydisperse colloids in subsurface such as clay particles, commonly observed in nature, that contaminant (e.g., heavy metals and radionuclides) are favorably adsorbed. In addition, sand also can be polydisperse, showing a coefficient of uniformity of about 10 [16]. This implies that the particle size distributions of sand and colloidal particles should be considered for the accurate prediction of contaminant transport.
The objective of this study was to numerically investigate the impact of the particle size distributions of sand and colloidal particles on contaminant transport. The breakthrough behavior of contaminant at the given mean and standard deviation of particle size distributions of sand and colloid particles was evaluated and discussed. In addition, the amount of immobile colloidal particles at given particles size distributions was evaluated using the model proposed by Won and Burns (2018) [17]. Consequently, the contaminant saturation and the bed efficiency of particle-retained sand medium were evaluated in order to represent simulated breakthrough curves quantitatively.

2. Mathematical Model

2.1. Colloidal Particle Transport in Saturated Porous Media

Three main mechanisms can be taken into account for colloidal particle–sand interaction during the transport of particles in saturated porous media: straining, attachment, and detachment [12,18]. Straining is the retention mechanism of particles caused by a relatively large particle size compared to the size of pore space [19,20]. Attachment is another retention mechanism of particles induced by the net attraction energy between the colloidal particles and sand, typically explained by the Derjaguin, Landau, Verwey, and Overbeek (DLVO) theory [21,22]. The detachment mechanism represents the separation of retained colloidal particles from sand induced by fluid drag forces [23,24]. The one-dimensional continuum scale advection–dispersion-type governing equation with consideration of the above three main mechanisms can be expressed as:
( n C c ) t = q ( n C c ) z + z ( n D c C c z ) n k a t t ψ a t t C c + ρ b k d e t S a t t n k s t r ψ s t r C c ,
ρ b S c t = n k a t t ψ a t t C c ρ b k d e t S a t t + n k s t r ψ s t r C c ,
where n (-) is the porosity, Cc (M L−3) is the particle concentration in aqueous phases, t (T) is the time, q (L T−1) is the Darcy’s velocity, z (L) is the depth, Dc (L2 T) is the dispersion coefficient of particles, ρb is the bulk density of sand (M L−3), and katt (T−1), kdet (T−1) and kstr (T−1) are the first-order coefficients accounting for attachment, detachment, and straining, respectively. ψatt and ψstr are the dimensionless functions for colloidal particle attachment and straining, respectively, and Sc (M M−1) is the solid-phase concentration of retained colloidal particles (Sc = Satt + Sstr, where Satt (M M−1) and Sstr (M M−1) are the solid-phase concentrations of attached and strained colloidal particles). ψatt [25,26], ψstr [12], katt [27], and kstr [12] in Equation (1) can be expressed as functions of the radius of colloidal particle (rc) and the radius of sand (rs) (see details in [17]):
ψ a t t ( r c , r s ) = 1 S a t t S max ( r c , r s ) ,
ψ s t r = ( d s 50 + z d s 50 ) β ,
k s t r ( r c , r s ) = 4.495 ( r c r s ) 1.42 ,
k a t t ( r c , r s ) = 3 ( 1 θ w ) 2 d s 50 α ( r c ) η 0 ( r c , r s ) v s ,
where Smax (M M−1) is the maximum attached mass of colloidal particles per unit mass of sand, ds50 (L) is the median size of sand, β (-) is a fitting parameter (β = 0.43, see [12]), α (-) is the attachment efficiency, η0 (-) is the single-collector collision efficiency [28,29], and vs (L T−1) is the average seepage velocity. Note that the Smax adopted in this paper was determined by the equation presented in [18], derived from the geometry of spherical colloidal particles and sand.

2.2. Colloidal-Particle-Associated Contaminant Transport in Saturated Porous Media

Assuming the all pore space is accessible to colloidal particles (i.e., no dead pores for colloidal particle transport) and colloidal particles are irreversibly retained during contaminant transport (i.e., no detachment), the mass balance equation for the colloidal-particle-associated contaminant transport in saturated porous media can be expressed as [9,30]:
( n C ) t + ρ S e t + ρ S k t + ρ S c S i c t   = q ( n C ) z + z ( n D C z ) ,
where C (M L−3) is the dissolved contaminant concentration in aqueous phase, ρ (M L−3) is the bulk density of sand medium, Se (M M−1) is the concentration of contaminant sorbed instantaneously to the sand, and Sk (M M−1) is the concentration of contaminant sorbed kinetically to the sand. In addition, Sic (M M−1) is the contaminant sorbed by retained colloidal particles. For three of the time derivative terms on the left-hand side associated with the adsorption of contaminant, the application of the Langmuir isotherm model on the contaminant–sand adsorption and the contaminant-retained colloidal particles adsorption leads to [31]:
ρ S e t   = ρ f K 1 q max ( 1 ) C t ,
ρ S k t   = ρ ω 1 [ ( 1 f ) q max ( 1 ) K 2 C 1 + K 2 C S k ] ,
ρ S c S i c t   = ρ ω 2 [ q max ( 2 ) K 2 C 1 + K 2 C S c S i c S c ] + k c ρ S c S i c ,
where f (-) is the fraction of instantaneously sorbed contaminant in the two-site adsorption model (f was assumed to be 0.01 in this study), Ki (L−3 M) is the distribution coefficient of dissolved contaminant to the adsorbent (i = 1 for sand, i = 2 for colloidal particle), qmax(i) (M M−1) is the adsorption capacity, ω (T−1) is the first-order rate constant, and kc (T−1) is the first-order reaction rate constant. The application of the Langmuir isotherm in Equations (8)–(10) keeps the amount of adsorbed contaminant less than qmax.

2.3. Sampling of Colloidal Particles and Sand

The Latin hypercube sampling (LHCS) method ([32]) was applied for sampling colloidal particles and sand at given particle size distributions. This method provides an even distribution of sample colloidal particles throughout the particle size distribution curve. The particle size distributions are represented by lognormal distribution functions and the mathematical formulation of mth probability values in the LHCS method (Pm) and mth sampled rc (rc(m)) at the given lognormal distribution function F is given as:
P m = ( 1 N ) u m + ( m 1 N ) ,
r c ( m ) = F 1 ( P m ) ,
where m is an integer between 1 and N, N is the sampling size, and um is the random number generated from 0 to 1 (0.5 was assumed in this paper for the reproducibility of sampling result). Using Equations (11) and (12) ensures the reflection of particle size distribution to the sampled rc at the given mean (µc) and standard deviation (σc), which will eventually affect the amount of retained colloidal particles by the calculated katt and kstr as presented in Equations (5) and (6), respectively.

2.4. Numerical Modeling

The equation for colloidal particle transport presented in Equation (1) was solved first at the inlet concentration of 5 g L−1 for 30 pore volumes (PVs, 1 PV represents the volume of colloidal particle suspension equal to the pore volume of sand medium) of injection. The obtained retention profile of colloidal particles (Sc in Equation (2)) was applied to the initial retained colloidal particles for solving the contaminant transport equation (Equation (7)). An inlet contaminant concentration of 5 mg L−1 was applied to the boundary condition for 30 PVs of injection. All parameters used in simulations were identical, except the mean and the standard deviation of the particle size distribution of colloidal particles. The backward Euler method was used for discretizing the governing equations, and the Picard iteration [33,34] was implemented in the model for the convergence in each time step. The dimension of the sand medium domain in all the numerical simulations was set to 30.48 cm (=1 foot) long, except the model validation (the dimension of domain = 4.61 cm). All parameters adopted in the numerical simulation are summarized in Table 1.

2.5. Bed Efficiency and Contaminant Saturation

The evaluation of breakthrough curves of contaminant allows the qualitative description of contaminant transport at varied µc and σc. However, obtaining the bed efficiency (β) and the contaminant saturation (θ) from the evaluated breakthrough curves provides the quantitative description of breakthrough curves in each condition. The β value represents the fraction of adsorbed contaminant during its transport:
β = 1 ( j = 1 P V f λ 0.5 λ ( ( C C 0 ) j + ( C C 0 ) j + 1 ) ) / P V f ,
where PVf (-) is the final elapsed pore volume (PVf = 30 in this work), λ (-) is the pore volume interval, and j is an order of normalized concentration starting from 0 PV (j = PVf/λ + 1), and C0 (M L−3) is the inlet concentration of contaminant. On the other hand, θ quantifies the saturation of bed at PVf (i.e., the fraction of occupied adsorption sites in the bed), expressed as:
θ = ( C i V P V β P V f ) / ( k = 1 N q max ( k ) W k ) ,
where VPV (L3) is the volume of 1 pore volume calculated by the porosity of the bed, N is the total number of materials (N = 2 in this study), qmax(k) (M M−1) is the maximum adsorption capacity of contaminant, and Wk (M) is the weight of the kth material in the bed.

3. Results and Discussion

3.1. Model Validation

While the equations for contaminant adsorption to sand adopted in this paper (Equations (8) and (9)) have been reported and used in the literature [31,37], the equation describing the adsorption mechanism of contaminant to retained colloidal particles presented in Equation (10) has not been reported. Therefore, the model presented in the previous sections was validated by comparing with laboratory experimental data as follows.
Two cases of the experimental data presented in [9] were selected to validate the model documented in the previous sections and to obtain ω2 and kc, which describe the adsorption rate of contaminant to retained colloidal particles. A column experiment was performed to investigate the transport of Pb(II) under the presence of retained kaolinite colloidal particles. The column with 5.08 cm diameter and 4.61 cm height was designed and uniformly packed with 1%–5% kaolinite content. After equilibrating the system using the background solution without Pb(II), 5 mg L−1 of Pb(II) solution was injected for 100 PVs. Conditions of the selected experiment cases for the model validation in this paper were 1% and 5% of kaolinite content at 6 mL min−1 of flow rate and pH ~6. Detailed information on the experimental procedure is documented in [9].
The trust-region algorithm was used to obtain the best curve-fitted breakthrough curves and ω2 and kc at given observed breakthrough curves (the step and function tolerance for the stopping criteria was set as 10−6). As seen in Figure 1, the proposed model in this paper simulated the observed breakthrough curves accurately. Even though the back-calculated coefficients in Figure 1 may be able to exclusively describe the Pb(II) transport under the presence of kaolinite, the values of ω2 = 0.016 s−1 and kc = 0.0152 s−1 were used in further simulation in order to investigate the impact of µc and σc of particle size distribution on contaminant transport at given ω2 and kc.

3.2. Impact of Median Colloidal Particle Size (rc)

Figure 2 and Figure 3 illustrate the contaminant breakthrough curves and the retention profile of colloidal particles at rc = 0.1, 0.5, 1, and 2 µm. While Figure 2 shows the concentration of contaminant at the outlet of the sand medium during the injection, Figure 3 presents the amount of retained particles with depth after injecting 30 pore volumes of particle suspension. As seen in Figure 2, the transport of contaminant was retarded when rc increased because of a more significant amount of retained colloidal particles at higher rc as illustrated in Figure 3. The higher value of rc led to greater katt and kstr in Equation (1) as presented in Figure 4 (except rc < 0.4 µm because of large η0 value at relatively low rc, see details in [29]), which resulted in a larger amount of retained colloidal particles. It can be anticipated that increases in the number of retained colloidal particles would correspond to slower the contaminant transport at the given adsorption rate. The results presented in Figure 2 imply that the size of the colloidal particles should be taken into account as one of the critical parameters in contaminant transport through porous media. That is, a slight increase in rc can retard the contaminant transport significantly. The retardation of transport would be more substantial if ω2 and kc were larger than the values considered in this paper.

3.3. Impact of the Standard Deviation of the Particle Size Distribution of Colloidal Particles

The median particle radii of 0.18 and 0.55 µm were selected, which correspond to the median particle diameter of kaolinite SA1 and kaolinite RP2 in [38], in order to investigate the impact of the standard deviation of particle size distribution. Figure 5 illustrates simulated breakthrough curves of contaminants and the retention profiles of particles at varied σc and µc = −1.72 (corresponding to the median particle size = 0.18 µm). It is clearly observed that the increase in the standard deviation of particle size distribution (i.e., an increase in the coefficient of uniformity) retarded the transport of contaminant (Figure 5a). This is attributed to the larger amount of retained particles at higher σc (Figure 5b), which is caused by larger values of sample rc at higher σc (Figure 6). As presented in [17], in the identical median size, the amount of retained particles increased as σc increased. In other words, no consideration of σc in colloidal particle transport in porous media would lead to a significant underestimation of the amount of retained particles. The results obtained in Figure 5 indicate that σc can be a critical factor not only for the colloidal particle transport but also for the contaminant transport when the contaminants are favorably adsorbed to retained colloidal particles.
The impact of σc became more significant when µc increased, as demonstrated in Figure 7. At µc = −0.6 (corresponding to the median particle size of 0.55 µm), the retardation of contaminant transport was more significant compared to the case of µc = −1.72, as illustrated in Figure 7a. This result may be explained mainly by the fact that an increase in the amount of strained particles due to increasing σc is greater at larger µc (showing a more exponential retention profile according to σc at µc = −0.6 (Figure 7b) than at µc = −1.72 (Figure 5b). The significant increase in retained particles at the top layer of the sand medium corresponding to the increase of σc provided a large number of adsorption sites at the top layer, which resulted in the significant retardation of contaminant transport.

3.4. Bed Efficiency and Contaminant Saturation

Figure 8 and Figure 9 illustrate the calculated β and θ using Equations (13) and (14) for Figure 2, Figure 5a, and Figure 7b. The value of β almost linearly increased as rc increased because of the retardation of contaminant transport, while θ decreased as rc increased, caused by the increase in retained colloidal particles. Forty percent of available adsorption sites (θ ~ 0.6) in the bed at rc = 2 µm after 30 PV of contaminant injection implies that the increase in adsorption sites affects the long-term transport behavior of contaminants as well as the retardation of transport. In addition, the evaluation of β provides the threshold rc at which the contaminant is not detected at the outlet after the given PVs. For instance, the extrapolation of β values presented in Figure 8 indicates that the contaminant would not be detected during 30 PV of injection at rc > 2.3 µm (i.e., β = 1 at rc ~ 2.3 µm).
In contrast, the increase in σc led to an increase in β and θ at µc = −1.72 and −0.6. However, while β increased significantly as σc increased, a slight increase in θ was observed at µc = −0.6. This implies the exponential increase in the adsorption capacity of the bed as σc increases when µc is relatively high. In other words, substantial adsorption sites were still available at high σc after the significant retardation of contaminant transport. As stated above, the calculation of β and θ provides the quantitative description of observed/simulated breakthrough curves as well as the prediction of the breakthrough point.

4. Conclusions

This work investigated the impact of particle size distribution for colloidal particles on contaminant transport by performing a numerical simulation. The colloidal-particle-associated contaminant transport model along with the Langmuir adsorption model for the reaction between contaminant–sand and contaminant–colloidal particles was introduced, and a parametric study was performed at varied mean and standard deviation of particle size distribution. Based on the observations from simulated breakthrough curves, the key findings can be summarized as follows:
  • The experimental breakthrough curves of Pb(II) under the presence of retained kaolinite particles were described well by the colloidal-particle-associated contaminant transport model introduced in this paper. This implies that the model can be used for the prediction of transport of heavy metals under the presence of immobile kaolinite particles or other particles that favorably adsorb the heavy metals.
  • As the median size of the colloidal particles and the standard deviation increased, the transport of contaminants was retarded. The impact of a standard deviation on the retardation of transport became more significant as the median size of the colloidal particles increased. This is attributed to the variation in the amount of retained colloidal particles by the median size and the standard deviation of particle size distribution. Therefore, the particle size distribution of particles should be taken into account in the prediction of contaminant transport.
  • A quantitative description of simulated breakthrough curves from the calculated β and θ provides the remaining available adsorption sites of the bed and the prediction of the breakthrough point at a given condition.
  • Further investigation is needed for the impact of particle size distribution on contaminant transport when particle and contaminant transport simultaneously. In addition, the unsaturated flow can also be considered in the modeling work to predict the contaminant transport from the surface to groundwater.

Author Contributions

Writing draft paper, J.W.; Investigation and discussion, J.W., D.L, K.P, H.L; Visualization, H.L; Supervision and editing the manuscript, H.C.

Acknowledgments

The authors appreciate the support partially by the Korea Agency for Infrastructure Technology Advancement under the Ministry of Land, Infrastructure and Transport of the Korean government (No. 18SCIP-B108153-04). The authors also wish to thank the anonymous reviewers for their careful review of the manuscript and valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McCarthy, J.F.; Zachara, J.M. Subsurface transport of contaminants. Environ. Sci. Technol. 1989, 23, 496–502. [Google Scholar] [CrossRef]
  2. Corapcioglu, M.Y.; Jiang, S. Colloid-facilitated groundwater contaminant transport. Water Resour. Res. 1993, 29, 2215–2226. [Google Scholar] [CrossRef]
  3. Sen, T.; Khilar, K.C. Review on subsurface colloids and colloid-associated contaminant transport in saturated porous media. Adv. Colloid Interface Sci. 2006, 119, 71–96. [Google Scholar] [Green Version]
  4. Kretzschmar, R.; Borkovec, M.; Grolimund, D.; Elimelech, M. Mobile Subsurface Colloids and Their Role in Contaminant Transport. In Advances in Agronomy; Elsevier Inc.: Amsterdam, The Netherlands, 1999; pp. 121–193. [Google Scholar]
  5. Torok, J.; Buckley, L.; Woods, B. The separation of radionuclide migration by solution and particle transport in soil. J. Contam. Hydrol. 1990, 6, 185–203. [Google Scholar] [CrossRef]
  6. Penrose, W.R.; Polzer, W.L.; Essington, E.H.; Nelson, D.M.; Orlandini, K.A. Mobility of plutonium and americium through a shallow aquifer in a semiarid region. Environ. Sci. Technol. 1990, 24, 228–234. [Google Scholar] [CrossRef]
  7. Grolimund, D.; Borkovec, M.; Barmettler, K.; Sticher, H. Colloid-Facilitated Transport of Strongly Sorbing Contaminants in Natural Porous Media: A Laboratory Column Study. Environ. Sci. Technol. 1996, 30, 3118–3123. [Google Scholar] [CrossRef]
  8. Denaix, L.; Semlali, R.M.; Douay, F. Dissolved and colloidal transport of Cd, Pb, and Zn in a silt loam soil affected by atmospheric industrial deposition. Environ. Pollut. 2001, 114, 29–38. [Google Scholar] [CrossRef]
  9. Won, J.; Burns, S.E. Role of Immobile Kaolinite Colloids in the Transport of Heavy Metals. Environ. Sci. Technol. 2018, 52, 2735–2741. [Google Scholar] [CrossRef] [PubMed]
  10. Sen, T.K.; Mahajan, S.P.; Khilar, K.C. Colloid-Associated contaminant transport in porous media: 1. Experimental studies. AIChE J. 2002, 48, 2366–2374. [Google Scholar] [CrossRef]
  11. Bradford, S.A.; Torkzaban, S.; Shapiro, A. A Theoretical Analysis of Colloid Attachment and Straining in Chemically Heterogeneous Porous Media. Langmuir 2013, 29, 6944–6952. [Google Scholar] [CrossRef] [PubMed]
  12. Bradford, S.A.; Simunek, J.; Bettahar, M.; van Genuchten, M.T.; Yates, S.R. Modeling Colloid Attachment, Straining, and Exclusion in Saturated Porous Media. Environ. Sci. Technol. 2003, 37, 2242–2250. [Google Scholar] [CrossRef] [PubMed]
  13. Elimelech, M.; O’Melia, C.R. Kinetics of deposition of colloidal particles in porous media. Environ. Sci. Technol. 1990, 24, 1528–1536. [Google Scholar] [CrossRef]
  14. Xu, S.; Gao, B.; Saiers, J.E. Straining of colloidal particles in saturated porous media. Water Resour. Res. 2006, 42. [Google Scholar] [CrossRef] [Green Version]
  15. Lenhart, J.J.; Saiers, J.E. Transport of Silica Colloids through Unsaturated Porous Media: Experimental Results and Model Comparisons. Environ. Sci. Technol. 2002, 36, 769–777. [Google Scholar] [CrossRef] [PubMed]
  16. Mitchell, J.K.; Soga, K. Fundamentals of Soil Behavior, 3rd ed.; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
  17. Won, J.; Burns, S.E. Stochastic Modeling of Kaolinite Transport through a Sand Filter. Can. Geotech. J. 2018. cgj-2018-0394. [Google Scholar] [CrossRef]
  18. Won, J.; Burns, S.E. Influence of Ionic Strength on Clay Particle Deposition and Hydraulic Conductivity of a Sand Medium. J. Geotech. Geoenvironmental Eng. 2017, 143, 04017081. [Google Scholar] [CrossRef]
  19. Xu, S.; Saiers, J.E. Colloid straining within water-saturated porous media: Effects of colloid size nonuniformity. Water Resour. Res. 2009, 45. [Google Scholar] [CrossRef] [Green Version]
  20. Bradford, S.A.; Simunek, J.; Bettahar, M.; van Genuchten, M.T.; Yates, S.R. Significance of straining in colloid deposition: Evidence and implications. Water Resour. Res. 2006, 42. [Google Scholar] [CrossRef] [Green Version]
  21. Derjaguin, B.V.; Landau, L.D. Theory of the stability of strongly charged lyophobic sols and the adhesion of strongly charged particles in solutions of electrolytes. Acta Physicochim. USSR 1941, 14, 633–662. [Google Scholar] [CrossRef]
  22. Verwey, E.J.W.; Overbeek, J.T.G. Theory of the Stability of Lyophobic Colloids; Elsevier: Amsterdam, The Netherlands, 1948. [Google Scholar]
  23. Bergendahl, J.; Grasso, D. Mechanistic Basis for Particle Detachment from Granular Media. Environ. Sci. Technol. 2003, 37, 2317–2322. [Google Scholar] [CrossRef] [PubMed]
  24. Bergendahl, J.; Grasso, D. Prediction of colloid detachment in a model porous media: Thermodynamics. AIChE J. 1999, 45, 475–484. [Google Scholar] [CrossRef]
  25. Treumann, S.; Torkzaban, S.; Bradford, S.A.; Visalakshan, R.M.; Page, D. An explanation for differences in the process of colloid adsorption in batch and column studies. J. Contam. Hydrol. 2014, 164, 219–229. [Google Scholar] [CrossRef] [PubMed]
  26. Sun, Y.; Gao, B.; Bradford, S.A.; Wu, L.; Chen, H.; Shi, X.; Wu, J. Transport, retention, and size perturbation of graphene oxide in saturated porous media: Effects of input concentration and grain size. Water Res. 2015, 68, 24–33. [Google Scholar] [CrossRef] [PubMed]
  27. Raychoudhury, T.; Tufenkji, N.; Ghoshal, S. Straining of polyelectrolyte-stabilized nanoscale zero valent iron particles during transport through granular porous media. Water Res. 2014, 50, 80–89. [Google Scholar] [CrossRef] [PubMed]
  28. Rajagopalan, R.; Tien, C. Trajectory analysis of deep-bed filtration with the sphere-in-cell porous media model. AIChE J. 1976, 22, 523–533. [Google Scholar] [CrossRef]
  29. Tufenkji, N.; Elimelech, M. Correlation Equation for Predicting Single-Collector Efficiency in Physicochemical Filtration in Saturated Porous Media. Environ. Sci. Technol. 2004, 38, 529–536. [Google Scholar] [CrossRef] [PubMed]
  30. Šimůnek, J.; He, C.; Pang, L.; Bradford, S.A. Colloid-Facilitated Solute Transport in Variably Saturated Porous Media. Vadose Zo. J. 2006, 5, 1035. [Google Scholar] [CrossRef]
  31. van Genuchten, M.T.; Wagenet, R.J. Two-Site/Two-Region Models for Pesticide Transport and Degradation: Theoretical Development and Analytical Solutions. Soil Sci. Soc. Am. J. 1989, 53, 1303. [Google Scholar] [CrossRef]
  32. McKay, M.D.; Beckman, R.J.; Conover, W.J. Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code. Technometrics 1979, 21, 239–245. [Google Scholar]
  33. Celia, M.A.; Bouloutas, E.T.; Zarba, R.L. A general mass-conservative numerical solution for the unsaturated flow equation. Water Resour. Res. 1990, 26, 1483–1496. [Google Scholar] [CrossRef]
  34. Huang, K.; Mohanty, B.P.; Leij, F.J.; van Genuchten, M.T. Solution of the nonlinear transport equation using modified Picard iteration. Adv. Water Resour. 1998, 21, 237–249. [Google Scholar] [CrossRef]
  35. Delgado, J.M.P.Q. Longitudinal and Transverse Dispersion in Porous Media. Chem. Eng. Res. Des. 2007, 85, 1245–1252. [Google Scholar] [CrossRef]
  36. Pang, L.; Šimůnek, J. Evaluation of bacteria-facilitated cadmium transport in gravel columns using the HYDRUS colloid-facilitated solute transport model. Water Resour. Res. 2006, 42. [Google Scholar] [CrossRef] [Green Version]
  37. Schijven, J.F.; Hassanizadeh, S.M.; de Bruin, R.H.A.M. Two-site kinetic modeling of bacteriophages transport through columns of saturated dune sand. J. Contam. Hydrol. 2002, 57, 259–279. [Google Scholar] [CrossRef]
  38. Palomino, A.M.; Burns, S.E.; Santamarina, J.C. Mixtures of fine-grained minerals kaolinite and carbonate grains. Clays Clay Miner. 2008, 56, 599–611. [Google Scholar] [CrossRef]
Figure 1. Observed and modeled breakthrough curves of Pb(II) under the presence of retained kaolinite colloidal particles at 1% and 5%. Back-calculated ω2 and kc are also presented with R2 values. PV: pore volume.
Figure 1. Observed and modeled breakthrough curves of Pb(II) under the presence of retained kaolinite colloidal particles at 1% and 5%. Back-calculated ω2 and kc are also presented with R2 values. PV: pore volume.
Applsci 09 00932 g001
Figure 2. Simulated breakthrough curves of the contaminant at inlet colloidal particle sizes ranging from 0.1 to 2 µm.
Figure 2. Simulated breakthrough curves of the contaminant at inlet colloidal particle sizes ranging from 0.1 to 2 µm.
Applsci 09 00932 g002
Figure 3. Simulated retention profiles of colloidal particles at in inlet colloidal particle sizes ranging from 0.1 to 2 µm.
Figure 3. Simulated retention profiles of colloidal particles at in inlet colloidal particle sizes ranging from 0.1 to 2 µm.
Applsci 09 00932 g003
Figure 4. Calculated katt and kstr using Equations (5) and (6) at varied rc (0.1 µm interval) and rs = 360 µm.
Figure 4. Calculated katt and kstr using Equations (5) and (6) at varied rc (0.1 µm interval) and rs = 360 µm.
Applsci 09 00932 g004
Figure 5. (a) Simulated breakthrough curves of contaminant and (b) the retention profiles of particles at µc = −1.72 and σc = 0, 0.1, 0.2, 1.0, 1.2, and 1.5.
Figure 5. (a) Simulated breakthrough curves of contaminant and (b) the retention profiles of particles at µc = −1.72 and σc = 0, 0.1, 0.2, 1.0, 1.2, and 1.5.
Applsci 09 00932 g005
Figure 6. Particle size distributions for colloidal particles (µc = −1.72) used in the simulation.
Figure 6. Particle size distributions for colloidal particles (µc = −1.72) used in the simulation.
Applsci 09 00932 g006
Figure 7. (a) Simulated breakthrough curves of contaminant and (b) the retention profiles of particles at µc = −0.60 and σc = 0, 0.1, 0.2, 1.0, 1.2, and 1.5.
Figure 7. (a) Simulated breakthrough curves of contaminant and (b) the retention profiles of particles at µc = −0.60 and σc = 0, 0.1, 0.2, 1.0, 1.2, and 1.5.
Applsci 09 00932 g007
Figure 8. Calculated β and θ for Figure 2.
Figure 8. Calculated β and θ for Figure 2.
Applsci 09 00932 g008
Figure 9. Calculated β and θ for Figure 5a and Figure 7a.
Figure 9. Calculated β and θ for Figure 5a and Figure 7a.
Applsci 09 00932 g009
Table 1. Parameters adopted in numerical simulation.
Table 1. Parameters adopted in numerical simulation.
ParametersValueNote and Related Reference
n0.365Updated every time step. Decreased over time as colloidal particles retained
N1 or 100N = 100 when σc > 0
Dc (10−3 cm2 s−1)6.1–6.4Calculated using equations proposed in [35], mean value was used if N > 1
D (10−3 cm2 s−1)2.9Experimentally obtained using non-reactive tracer (bromide)
q (cm s−1)0.0123-
ρb (g cm−3)1.68The specific gravity of sand was assumed 2.65
ds50 (cm)0.072Typical median size of coarse sand
µc0.1–2Typical median size of clay colloidal particles (rc50 = eµ)
σc−2.3 to 0.69Equivalent to the coefficient of uniformity ranged from 0 to 10
f (-)0.01Assumed value ([36])
K1, K2 (10−3 L mg−1)2.202, 6.102Obtained by batch experiment ([9])
qmax(1), qmax(2) (mg g−1)1.081 × 10−3, 2.854Obtained by batch experiment ([9])
ω11.859 × 10−3Obtained by batch experiment ([9])

Share and Cite

MDPI and ACS Style

Won, J.; Lee, D.; Pham, K.; Lee, H.; Choi, H. Impact of Particle Size Distribution of Colloidal Particles on Contaminant Transport in Porous Media. Appl. Sci. 2019, 9, 932. https://doi.org/10.3390/app9050932

AMA Style

Won J, Lee D, Pham K, Lee H, Choi H. Impact of Particle Size Distribution of Colloidal Particles on Contaminant Transport in Porous Media. Applied Sciences. 2019; 9(5):932. https://doi.org/10.3390/app9050932

Chicago/Turabian Style

Won, Jongmuk, Dongseop Lee, Khanh Pham, Hyobum Lee, and Hangseok Choi. 2019. "Impact of Particle Size Distribution of Colloidal Particles on Contaminant Transport in Porous Media" Applied Sciences 9, no. 5: 932. https://doi.org/10.3390/app9050932

APA Style

Won, J., Lee, D., Pham, K., Lee, H., & Choi, H. (2019). Impact of Particle Size Distribution of Colloidal Particles on Contaminant Transport in Porous Media. Applied Sciences, 9(5), 932. https://doi.org/10.3390/app9050932

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

Article Metrics

Back to TopTop