Abstract
We present model equations for the trapping and accumulation of particles in a cylindrical channel with nanostructured inner walls when a fluid passes through, carrying a moderate load of impurities. The basic ingredient of the model is the introduction of a phenomenological ‘effectivecharge density’ of the walls, related to the electrical charges exposed in the nanotexture. The effective charge is gradually reduced as the flow runs through the channel and the trapped impurities cover the internal walls. Based on the proposed equations, the position and time dependence of the areal density of trapped impurities, and the filtration performance, may be calculated. It is proposed that experimentally testing these results may help to understand the enhanced trapping capability observed in many diverse nanotextured channel structures.
Keywords:
Nanostructures embedded in larger systems; Fluid impurity trapping in nanostructures; Effectivecharge densityBackground
Recently, [18] researchers from both academia and industry have experimentally demonstrated that a variety of nanostructures and nanotextured media (see some examples later in this introduction) can efficaciously trap nanoimpurities carried by fluids when a flow is induced by external hydrostatic pressure. These findings are not only scientifically interesting, but also promising for the socially and economically important application of purification of drinking water and other liquids [4,79]. When compared to conventional porous filters, the new media have the important advantages of retaining impurities of sizes typically in the tens of nanometers and, at the same time, presenting a resistance to hydrodynamic flow orders of magnitude smaller than what conventional models would predict for channels of diameters as small as the particles being trapped.
Roughly, we can divide the structures presenting such enhanced impurity trapping capability into two groups: (a) The first group corresponds to those formed by nanometricdiameter channels through which the fluid flows [14]. A wellknown example is the nanotube arrays grown and experimentally tested by Srivastava and coworkers [1]. Other specially interesting examples are graphene membranes although, by now, they have been probed only through molecular dynamics simulations [2]. In any nanometricdiameter channel, simple size exclusion will play a major role in the retention of nanoimpurities. However, in addition, these structures also exhibit remarkable capability to trap some ions significantly smaller than the channels’ diameter [1,2]. The resistance to flow is observed to be well lower than what conventional models predict for these diameters, a phenomena often attributed to waternanostructure interactions (see, e.g., [1]) though not yet fully understood at the quantitative calculation level. (b) The second group corresponds to nanostructures embedded in larger structures, resulting in filters composed by channels with micrometric diameters and inner walls coated with nanoparticles. Examples are conventional microfilters coated with Y_{2}O_{3}[5], ZrO_{2}[6], or Al_{2}O_{3}[7,8] nanopowders (further examples can be found in the reviews [3,4,9]). These structures have been observed by their growers to have a surprisingly good filtration performance for nanometric impurities, as small as approximately 10 nm, in spite of the relatively large diameter of the channels (note that in a channel with a diameter of 1 μm only about 0.04% of the fluid will transit closer than 10 nm from the walls) [39]. Their hydrodynamic resistance is quite low, similar to the one of conventional micrometric filters. Their trapping capability is observed to depend on pH and zeta potential [58] and, thus, electrostatic and polar attraction may be suspected to play a significant role in the filtration mechanism and dynamics. However, attempts to modelize them have been scarce. The authors of [7,8] empirically characterized their filters using generalpurpose plugflow adsorption models, like those used for column chromatography, and fitting the Langmuir and BET isotherms. While these models may be very useful for the design of improved filters, further theoretical work modelling the trapping capability and time evolution explicitly in terms of the charge of the walls, and also more detailed measurements, seems to be needed to further understand the impurity trapping mechanism and dynamics of this type of channels with nanostructured walls.
The purpose of this ‘Nano Idea Letter’ is to propose a specific model for the nanoimpurity trapping capability of cylindricallike channels with nanostructured inner walls of the type composing filters of category ‘b’ in the previous paragraph. We explore theoretically a simplified but realistic view in which the improved filtration capability is primarily due to the fact that the nanotexturing exposes electrical charges in the walls which induce both electrostatic and van der Waals attractions over the impurities in the fluid. This nanostructuring also provides chemical anchors for the binding of those impurities once they collide with the channel walls. Correspondingly, our basic ingredients will be the introduction of an effectivecharge density, z_{e}, of the inner walls of the channels and writing down as a function of z_{e} the impurity trapping probability. As it could be expected, z_{e }will depend on the areal density n of impurities already trapped in the inner walls of the channel. We obtain within the model the evolution of n and z_{e }with position x and with time t when the liquid is flowing through the channel.
The model produces agreement with the initial trapping performances quantitatively reported by experimentalists in various systems. Also, we propose that further detailed measurements as a function of time may be crucial to test these ideas more thoroughly. We believe that some aspects of the model could also be useful to partly explain the trapping of the smaller ions in the nanodiameter channels of category ‘a’. However, its full applicability to that case is limited by our use of classical dynamics for the carrying fluid. Hence, we do not focus here on that category (also, for these nanodiameter channels, in which the number of fluid atoms is manageably small, molecular dynamics simulations as those in [2] could be a more reliable, albeit not general, approach).
Obtainment of an equation for the areal density of trapped impurities in a channel with nanostructured walls
Initial modelling and notations
Our starting point, and most of our basic notations, is illustrated in Figure 1. We consider a channel with nanostructured inner walls, its nominal shape being cylindricallike with average radius r_{0}and length L. We divide it into slices along the axial coordinate x, each with differential thickness dx. A fluid flows through the channel due to externally applied hydrostatic pressure, carrying a load of impurities. Some of those impurities will become trapped by the inner wall of the channel, then reducing its effective radius to a value r_{e}(x,t) = r_{0}−r_{1}n(x,t), where t is the time, n is the number areal density of trapped impurities in the inner wall, and r_{1}is a constant proportional to the average volume of impurities. Throughout this letter, by ‘areal density’ we refer to quantities normalized using the nominal area of the inner wall (2Πr_{0}dx for a differential slice) and not the cross section of the channel. Also, for simplicity, we consider all impurities equal among them (subsequent generalization to multiple chemical species should be easy). The average radius of the impurities is noted ρ_{0}. The impurity concentration in the fluid is considered to be moderate enough as to not significantly affect its viscosity and as for the impurities in the fluid to be noninteracting with each other (specially when colliding with the channel wall).
Figure 1. Representation of a nanostructured channel filter as modelled in the present letter. The nominal shape of the channel is supposed to be cylindrical with length L, and the figure shows only the differential slice with axial coordinate from x to x + dx. The radiuses r_{0 }and ρ_{0 }correspond to the average dimensions of the bare channel and impurities. The effective radiuses r_{e }and ρ_{e }vary as trapped impurities cover the inner wall, via their dependences on, respectively, the areal density n of trapped impurities and on the areal density z_{e }of effective charge of the inner wall. This z_{e }reflects that exposed charges in a nanostructured surface attract the impurities in the fluid and also constitute binding anchors for those impurities. It is expected to diminish as impurities cover the surface, for which we assume the simple z_{e}(n) dependence given by Equation 1 of the main text.
Effectivecharge density of the inner wall, z_{e}
We now introduce the important concept of a phenomenological ‘effective charge’ of the inner wall of the channel. We quantify this effective charge via its areal density z_{e}, and as already commented on in the introduction, it reflects the fact that the nanostructured walls expose charges that induce both electrostatic and van der Waals attractions over the components of the impurities in the fluid. Indeed, z_{e} will depend on the areal density of already trapped impurities n (which will screen out the wall) and also on the chemistry specifics of the wall and impurities. Let us focus on the mutual interplays between n and z_{e} and in obtaining an equation for their evolution with time as flow passes through the channel. In particular, the interdependence z_{e}(n) may be naturally expected to be continuously decreasing when n increases, to take a finite value z_{0} at n = 0 (clean filter), and to saturate to zero when n reaches some critical value n^{sat} at which all active centers of the wall become well covered by impurities. We thus postulate the simplest z_{e}(n) dependence fulfilling such conditions:
where the notation ∥…∥ stands for min{1,…}. Obviously, other sensible choices for z_{e}(n) are possible such as, e.g., z_{e}(n) = z_{0}(1−∥n/n^{sat}∥)z_{1}, with z_{1} as a positive coefficient that probably covers at a good approximation most actual possibilities depending on its value, as it corresponds to a z_{e}(n) functionality with downwards curvature if 0 < z_{1 }< 1, to no curvature for z_{1 }= 1 (i.e., Equation 1), and to upwards curvature for z_{1 } >1. For simplicity, we shall consider in this letter only the case z_{1 }= 1. Note also that z_{e} may depend on position and time via the n(x,t) dependence.
Impurity trapping probabilities as a function of z_{e}and n
The role played by z_{e }in our model will be in fact twofold. First, it affects how large the distance is
within which if the impurity approximates the inner wall then the latter attracts
the former so much as to consider it as a collision. This attraction distance may
be seen as an effective radius, ρ_{e}, of the impurity (see Figure
1), so if the distance from the center of the impurity to the center of the channel
is larger than r_{e}−ρ_{e}, the impurity will actually touch the wall (dressed with already trapped impurities).
Let us discuss the ρ_{e}(z_{e}) dependence. We consider first the simplest case of an unscreened electrostatic interaction,
in which the potential energy of an impurity at a distance ρ_{e} from the wall is
In the righthand side of this equation, for convenience, we have expressed the thermal kinetic energy in units of the unscreened potential in the clean channel at a distance ρ_{0} from the surface, so ρ_{1} is an nondimensional coefficient proportional to T. We have also taken into account the finite bare size of the impurities by using ρ_{e}−ρ_{0} instead of ρ_{e} in the potential energy term. From the above equation, ρ_{e }can be obtained with the help of the principal Lambert W function as follows:
Although W(x) can be easily evaluated by modern computers, it is worthwhile to mention its asymptotes
W(x)≃x for x ≪ 1 and
We may now obtain the fraction f_{e }of impurities that flow, at given t and x values, near a collision distance from the impuritydressed wall. For that, we assume
that the fluid velocity profile is given by the Poiseuille law,
[10]
In the considerations leading to Equation 4, we have implicitly taken the concentration
of impurities as constant along the radial coordinate r. However, in principle, it could be expected that near the walls the electric potential
will influence the distance between impurities. To test whether this effect may be
of relevance, a Debyelike concentration profile
The second influence played by z_{e }in our model concerns the probability that an impurity gets actually bound to the
inner wall of the channel once it actually is within a collision distance from that
wall. We express the probability that a given impurity entering a differential slice
of the channel with thickness dx gets trapped in that slice as
Equation for ∂n(x,t)/∂t
Let us now build, on the basis of the above relationships, equations for the evolution of the areal density of trapped impurities, n, as a function of time t and position x when an impure fluid flows through the channel due to hydrostatic pressure. We start by considering the differential channel slice going from x up to x + dx, and we write down the expression for ∂n/∂t evaluated at that slice, which may evidently be expressed as Φ_{imp}Ω_{trap}/(2Πr_{0}dx), where Φ_{imp} is the flow of impurities brought by the incoming fluid (in units of s^{−1}; the factor (2Πr_{0}dx)^{−1 }is due to the areal density normalization in the definition of n). In its turn, Φ_{imp} can be written as Φ_{imp }= C_{imp}Φ where Φ is the fluid flow and C_{imp }the incoming number concentration of impurities. Gathering together the previous results in this letter, we get
with the z_{e}(n) and ρ_{e}(z_{e}) dependences given by Equations 1 and 3.
Equations for Φ(t)and ∂C_{imp}(x,t)/∂x
In order to solve the filtration dynamics (i.e., to obtain n(xt) and C_{imp}(xt)), it is necessary to supplement Equation 5 with formulas for Φ(t) and C_{imp}(xt). Regarding the fluid flow, we apply the Poiseuille law for incompressible fluids of viscosity η in a cylindrical channel of length L and radius r_{e}(xt): [10]
In this equation, P is the pressure difference between both ends of the finitelength channel, which
we take constant with time. Note that Φ becomes zero when at some x, the n value becomes n^{clog }≡ r_{0}/r_{1}, i.e., r_{e} becomes zero at that location and the channel becomes fully closed by impurities.
Note also that Equation 6 reduces in the particular case r_{0 }≫ r_{1}n(x,t) (which is common in experiments) to
We construct now the supplementary equation for C_{imp}(x,t). For that, we again consider the differential channel slice going from x up to x + dx. The number of impurities that become trapped in its walls during an interval dt is (2Πr_{0}dx)(∂n/∂t)dt (the factor 2Πr_{0}dx is again due to the areal normalization in the definition of n). The numbers of impurities entering and exiting the slice in the liquid flow are Φ(t)C_{imp}(x,t)dt and Φ(t)C_{imp}(x + dx,t)dt respectively. Mass conservation balance therefore gives
Notice that Equations 5 to 7 are coupled to each other. In fact, they form now a closed set that can be numerically integrated by providing the specific values for the characteristics of the filter, for any given pressure difference P and incoming impurity concentration C_{imp}(0,t). In what follows, for simplicity, we will always consider for the latter a constant value C_{0}. The computation to numerically integrate Equations 5 to 7 is relatively lightweight (e.g., calculating our Figure 2 took about 15 min in a current personal computer that considered 2 × 10^{4 }finiteelement xslices).
Figure 2. Time dependence. (a) Results, obtained by integrating Equations 5 to 7, for the time dependence of the
areal density of trapped impurities (continuous lines) at the entrance of the channel
n(x = 0,t) and at its exit point n(x = L,t), and also the global average areal density of trapped impurities
Discussion of the results obtained by integrating the model equations
Numerical integration and comparison with some existing partial measurements
We show in Figure
2 an example of the results obtained by numerically integrating Equations 5 to 7 using
some representative values for the parameters involved (and always in the case of
constant P and C_{imp}, and starting from a clean initial state n(xt = 0)=0). In particular, we have chosen parameter values that reproduce the case of
channels coated with Y_{2}O_{3} nanopowders as measured in
[5] (they are essentially valid also for the quite similar case of channels with ZrO_{2} nanocoating reported by the same group in
[6]). In these filters, the channels have a typical value of the nominal radius r_{0 }= 500 nm and length L = 7.25 mm. They were shown
[5] to efficaciously retain MS2 viruses (of radius ρ_{0 }= 13 nm) carried by water with NaCl as background electrolyte and a conductivity of
400μS/cm (corresponding then to λ_{D}≃5.1 nm) feed at a pressure P = 3 bar. The incoming impurity number concentration was
Figure
2a presents the results corresponding to integrating Equations 5 to 7 using these parameter
values, for the areal density of trapped impurities at the entrance of the channel
n(x = 0,t) and at its exit point n(x = Lt) and also for the global average areal density of trapped impurities
As a further test, we have repeated the same numerical integration as in Figure
2a,b but considering a radial impurity concentration profile
Linear, logarithmic, and saturated approximations
In Figure
2a, it is possible to identify in our results for the areal density of trapped impurities
some tranges in which the tdependence is relatively simple: (1) The initial time behavior is an approximately
linear n(t) growth; (2) in the intermediate regime, the growth of n(t) becomes approximately logarithmic; and (3) at sufficiently large t values, the saturation limit is reached, in which n approaches a value n^{sat} at a slow pace. These regimes are easily seen in Figure
2a for n(x = 0,t), n(x = L,t), and
In the case of a very short cylindrical channel (so that all xderivatives may be neglected), it is possible to find analytical expressions for the n(t) evolution in the linear and logarithmic regions: For the linear regime, by just introducing in Equation 5 the condition t ≃ 0, we find:
with
The logarithmic regime can be found by using the condition n ≃ n^{sat}/2:
with
In obtaining the above Equations 8 to 11, we have assumed that n(0) = 0 and that ρ_{e }< r_{e} at t = 0 or t_{1/2}.
Conclusions and proposals for future work
This letter has proposed a model for the main generic features of the channels with nanostructured inner walls with respect to trapping and accumulation of impurities carried by fluids. This includes, e.g., their capability to clean the fluid from impurities of a size much smaller than the channels’ nominal radius, with comparatively small resistance to flow (much smaller than in conventional channels with a radius as small as the impurities). The model attributes the enhanced filtration capability to the longrange attraction exerted by the exposed charges in the nanostructured walls and also to their binding capability once the impurities actually collide with them. Both features were quantitatively accounted for by means of a phenomenological ‘effectivecharge density’ of the nanostructured wall. The model also predicts the time evolution of the trapped impurity concentration and of the filtering capability, including three successive regimes: a linear regime, a logarithmic regime, and the saturated limit.
We believe that our equations could make possible some valuable future work, of which two specific matters seem to us more compelling: First, it would be interesting to check at the quantitative level the agreement with experiments of the time evolutions predicted above. For that, we propose to perform timedependent measurements made in controlled flow setups. We have chosen in our equations flow constraints which seem appropriate for this purpose; however, the model can be tested with different setups by just numerically integrating Equation 5 under the corresponding experimental constraints. One of the prime purposes of this letter is in fact to suggest such measurements.
A second interesting future work (already in progress in our research group) is the design of optimal geometries for the combinations of channels forming filters. Our model opens this possibility because of the explicit use of the Poiseuille relations that allow the calculation of the resistance to flow of complex associations of those channels, in series and/or parallel. The effective diffusivity and tortuosity of the pathways’ network are also accounted for by these equivalentcircuit analyses.
Competing interests
The author declares that he has no competing interests.
Acknowledgements
This work has been supported by the MICINN project FIS201019807 and by the Xunta de Galicia 2010/XA043 and 10TMT206012PR projects. All projects are cofunded by ERDF from the European Union.
References

Srivastava A, Srivastava ON, Talapatra S, Vajtai R, Ajayan PM: Carbon nanotube filters.
Nat Mater 2004, 3:610614. PubMed Abstract  Publisher Full Text

CohenTanugi D, Grossman JC: Water desalination across nanoporous graphene.
Nano Lett 2012, 12:36023608. PubMed Abstract  Publisher Full Text

Humplik T, Lee J, O’Hern SC, Fellman BA, Baig MA, Hassan SF, Atieh MA, Rahman F, Laoui T, Karnik R, Wang EN: Nanostructured materials for water desalination.
Nanotechnology 2011, 22:292001292019. PubMed Abstract  Publisher Full Text

Theron J, Walker JA, Cloete TE: Nanotechnology and water treatment: applications and emerging opportunities.
Crit Rev Microbiol 2008, 34:4369. PubMed Abstract  Publisher Full Text

Wegmann M, Michen B, Graule T: Nanostructured surface modification of microporous ceramics for efficient virus filtration.
J Eur Ceram Soc 2008, 28:16031612. Publisher Full Text

Wegmann M, Michen B, Luxbacher T, Fritsch J, Graule T: Modification of ceramic microfilters with colloidal zirconia to promote the adsorption of viruses from water.
Water Research 2008, 42:17261734. PubMed Abstract  Publisher Full Text

Tepper F, Kaledin L: Nanostructured chembio nonwoven filter. In Nanoscience and Nanotechnology for Chemical and Biological Defense, Volume 1016. Edited by Nagarajan R, Zukas W, Hatton TA, Lee S. Washington: American Chemical Society; 2009:273288.

Tepper F, Kaledin L, Kaledin T: Nonwoven electrostatic media for chromatographic separation of biological particles.
J Liq Chromatogr Related Technol 2009, 32:607627. Publisher Full Text

Meridian Institute: Workshop on nanotechnology, water and development (Chennai, India 2006) [http://www.merid.org/nanowaterworkshop].

Landau LD, Lifschitz EM: Fluid Mechanics. Oxford: Pergamon Press; 1987.

Sparreboom W, van den Berg A, Eijkel JCT: Transport in nanofluidic systems: a review of theory and applications.
New J Phys 2010, 12:015004015027. Publisher Full Text