The Reynolds number (Re) is defined as Re =
ρ Umax Dh/
μ (
ρ is the fluid density,
Umax is the characteristic velocity,
Dh is the hydraulic diameter,
μ is the mean viscosity), which is used to quantify the relative importance of the inertial and viscous effect for a given flow condition. A dimensionless Weissenberg number (Wi) is defined as Wi = λ
= 2 λ
Umax/
Dh (λ is the relaxation time of the polymer solution,
is the shear rate), which is used to quantify the viscoelastic effects of a fluid and describes the ratio of elastic and viscous forces. To compare the relative importance of the elastic effect to the inertial effect, the fluid elasticity (El) is defined as the ratio of the Weissenberg number to the Reynolds number, i.e., El = Wi/Re = 2 λ
μ/(
ρ Dh2) [
19], which is a dimensionless number and can determine the number of particle focusing points under a non-Newtonian viscoelastic fluid.
3.1. Viscoelastic Particle Focusing under Non-Newtonian Fluid
First, we investigated and observed the focusing and migration of particles under a viscoelastic fluid in three kinds of microchannels with different T-shaped cross-sections. Furthermore, experiments were carried out for two sizes of particles (5 and 13 μm diameter), which is one of the important parameters for determining elasto-inertial particle focusing behaviors [
20].
Figure 4 shows the particle focusing at a flow rate of 100 μL/min according to the location from the inlet of the Type 3 microchannel in non-Newtonian viscoelastic fluid (PEO500). At the location near the inlet (1 cm), no focusing was observed for either particle size due to the insufficient elastic forces developed. As the channel length was increased to 7 cm, particles were focused in a single stream located at the center, with a flow rate of 100 μL/min. The lateral migration speed was strongly dependent upon particle size [
21,
22], so tighter particle focusing was obtained for the 13 μm particle 7 cm downstream from the inlet.
Figure 5 presents the optical images of particle focusing (top and side view) in the microchannel with different T-shaped cross-sections for different flow rates and particle sizes in non-Newtonian viscoelastic fluid. All optical images were captured 7 cm downstream. The dotted lines are the boundaries of the microchannel from the top and side view, which were obtained from the bright-field optical images.
Figure 6 presents the normalized fluorescent intensity profiles presenting the particle distribution at various flow rates in the microchannel with different T-shaped cross-sections including a microchannel with a square cross-section which was prepared for comparison.
As shown in
Figure 5a,c,e, when viewed from the side, 13 μm particles were observed to focus to a single line near the bottom side for all types of channels. Furthermore, the effect of elasto-inertial force was dominant for low flow rates ((below 50 μL/min, Re ≈ 0 and Wi > 0) so the particles moved to the centerline for all types of channels (top view). The experimental results show that there was only one stream from both front and side views, meaning that true 3D elasto-inertial particle focusing was achieved for the tested condition. On the other hand, when the flow rate increased and reached 50 μL/min (Re > 0 and Wi > 0), the single-line particle focusing changed to double-line focusing in the Type 1 channel, but there was no change in the Type 2 and Type 3 channels. For the higher flow rate of 100 μL/min (inertial force dominant regime), due to the relatively strong inertial force at the high flow rate [
6], the focused particles were more or less dispersed around the centerline in the Type 2 channel and the triple-line focusing was observed in the Type 3 channel. This implies that the lateral migration of particles in a non-Newtonian fluid is strongly affected by the inertial force rather than the elastic force [
20].
As shown in
Figure 5b,d,f, particles with a relatively small particle size (5 μm) were randomly dispersed and not focused at specific points in the Type 1 and Type 2 channels for low flow rates (below 50 μL/min, Re ≈ 0 and Wi > 0). In the case of the Type 3 channel, however, small particles were slightly focused in the centerline for low flow rates. As the flow rate increased up to 50 μL/min (Re > 0 and Wi > 0), the particles were focused in the center of the microchannels under the influence of elasto-inertial force. However, unlike the large particles (13 μm), the small particles were not completely focused in the center, showing wider focused bands. This was due to the competition between the non-negligible inertial lift force and the relevant elastic force under viscoelastic fluid flowing [
6,
20].
Consequently, the particle migration and focusing were more distinct for large particles because the force on the particle induced by the normal stress was proportional to the cube of the particle diameter [
20]. In the cases of the Type 1 and Type 2 channels, tighter focusing was generally observed for the large particles such as the square channel, but the small particles were rather tightly focused in the center of the channel at high flow rates (over 100 μL/min) in the Type 3 channel. Moreover, the Type 3 channel showed a peculiar phenomenon in which large particles did not gather but moved to special focusing points. In total, there were three focusing points.
Figure 7 shows the positions of 5 μm and 13 μm particles distributed within the channel at the flow rate of 100 μL/min for three types of channels. As the aspect ratio increased (from Type 1 to Type 3), small particles were more tightly focused in the center. However, there was not enough force acting on the large particles in the side part of the channel to move them to the center, and they were focused in the side part of the channel. This will be explained in the section reporting the simulation.
3.2. Simulation Results
In order to first assess the flow pattern without the inertial particle on the channel cross-plane,
Figure 8a,b show contour maps of the axial velocity (
U) at the flow rate of 10 μL/min as a representative case for Type 2 and Type 3, respectively. The axial velocity was normalized by its maximum (
Umax), while we used an actual physical unit for both the horizontal and vertical axis.
Figure 8 shows that the aspect ratio of the channel significantly affects the mean flow structure on the cross-plane. For example, the horizontal location of
Umax for both Type 2 and Type 3 is the same, residing on the symmetric plane (i.e.,
x = 0). However, the vertical location of
Umax is different depending on the aspect ratio of the channel cross-section. For Type 2, the local maximum in
U occurs below the junction (i.e.,
y = 10 μm), whereas Type 3 has a local maximum above the junction on the vertical axis (i.e.,
y = 25 μm). This is because the shallow bottom channel for Type 3 induces a higher-pressure zone that moves the location of
Umax further away from the bottom wall. Furthermore, the presence of the higher-pressure zone for Type 3 affects the mean flow topology, as shown in
Figure 8a,b. Since H/W = 1 for Type 2, most of the fluid travels through the bottom channel, while the fluid in Type 3 mainly goes through the upper channel.
The contour maps of the shear rate for both channels are presented in
Figure 8c,d. As reported in the previous study [
23], the shear rate as a differential of the fluid velocity plays a crucial role in particle migration in the viscoelastic regime. This is because the shear rate is closely correlated with the competing mechanism of the inertial lift force between the shear gradient and wall interaction. Furthermore, previous studies [
7,
23] have shown a preferential migrating pattern of the inertial particles toward the region where the shear rate is zero, particularly when the second normal stress difference in a viscoelastic fluid is negligible. With these notions,
Figure 8c,d give a qualitative picture of the migration pattern of the elasto-inertial particles for the cross-section shapes considered in this study. It is clear that each wall created a relatively higher shear rate, whereas the center and each corner of the channel cross-section had the lowest shear rate for both Type 2 and Type 3, as observed in the previous studies [
3,
11,
20]. Most importantly, the protruding sharp edges at the junction caused the highest shear rate for both channels. The effect of these sharp edges became more significant when they were coupled with the bottom wall effect as the aspect ratio increased (see
Figure 8d). For example, the topology of the low shear rate (i.e.,
< 1000 s
−1) seems bottleneck-like for Type 3. Based on these observations, it can be expected that the migrating pattern of the elasto-inertial particles passing through the Type 3 channel would be different to those in the Type 2 channel.
Figure 9 presents the variation in x- and y-directional net elasto-inertial force (F
net) for 5 and 13 μm particles in the Type 2 channel along the aforementioned lines, as illustrated in
Figure 2. The normalization was made for F
net by ρ
f U
max2 d
4/D
h2, where ρ
f and U
max denote the fluid density and the maximum velocity of the channel. As mentioned earlier, the horizontal axis in
Figure 9 used an actual physical unit, which is better for describing the particle location with respect to the origin. Symbols for each case are denoted in the legend in
Figure 9. Solid lines connecting through the data points were made by piecewise cubic interpolation to show a shape-preserving trend in F
net along the x and
y axis. Note that the current simulation could not resolve a valid F
net after a certain point when the spherical particle resided in the proximity of the channel wall. It should also be noted that the negative value in F
net when x (or y) > 0 represents a directional force toward the origin, while its positive value when x (or y) < 0 indicates a directional force toward the wall.
For both the 5 and 13 μm particles shown in
Figure 9a,b,
Fx,net is always negative and has a single zero crossing at the origin at a flow rate of 100 μL/min for Type 2 (i.e., T
2D5Q100, T
2D13Q100, red symbols). This indicates that both 5 and 13 μm particles suspended in a viscoelastic fluid are likely to drift toward the origin in a lateral direction when the fluid inertia is sufficient. This numerical observation agrees well with our experiments, showing a single train of particles for both sizes at the same flow rate, as shown in
Figure 5c,d. However, at the lower flow rate (i.e., 1 and 10 μL/min) in
Figure 9a,b, a different pattern in
Fx,net between two inertial particles can be observed. For example,
Fx,net for the larger particle (i.e., 13 μm) still shows a negative with a single equilibrium point at the origin despite some fluctuations in
Fx,net at the lowest flow rate (i.e., 1 μL/min). This means the lateral drift consistently occurred for large particles (13 μm) independent of the flow rate within our measured range. For the case of small particles (5 μm), on the other hand,
Fx,net has multiples of zero crossings at the lower flow rates, as shown in
Figure 9a, indicating wide dispersion of inertial particles along the lateral direction. These numerical results confirm a unique migrating behavior of the inertial particles depending on the particle size and the flow rate reported in our experiments (see
Figure 5d).
Similar drifting patterns of elasto-inertial particles can also be observed for the vertical direction in
Figure 9c,d. For the large particle size (i.e., 13 μm), the equilibrium point in
Fy,net is made independent of the current flow rates and consistently resides at 10 μm away from the vertical origin for all cases (see
Figure 9d). The case of the small particle also shows a similar trend in
Fy,net and an equilibrium point at the highest flow rate, as shown in
Figure 9c. This observation in our numerical study explains the particle focusing when the effect of the elasto-inertial force is large enough due to either a high degree of the flow rate or a large particle surface area. Furthermore, the vertical equilibrium point observed herein is where the local minimum
and the local maximum
U are positioned.
In this study, our experimental result shows a different migrating pattern of particles depending on their size, particularly in the Type 3 channel, and thus has a potential implication for particle separation by leveraging the aspect ratio of the current T-shaped channel. Therefore, in this numerical study, it is of particular importance to confirm such migrating behavior can be observed in the Type 3 channel and explore their background mechanism compared to the one in the Type 2 channel. Note that the x-directional net forces are only considered in this case since y-directional net forces for Type 3 have the same trend and equilibrium point as those of Type 2. It should also be noted that the maximum flow rate increases at 150 μL/min as the characteristic drift for the larger particles started to occur at the flow rate of 150 μL/min in the present experiment.
Figure 10 shows the variation in the
x-directional net force exerted on two different particles varying with the flow rates. For the case of the 5 μm particle, the first two lowest flow rates (i.e., T
3D5Q1 and T
3D5Q10) generate multiples of the zero crossing in the x-directional force, indicating an absence of the particle focusing in the lateral direction, as shown in
Figure 10a. With an increase in flow rate, however, a consistent directional force occurs and pushes the 5 μm particle toward the origin of Type 3, indicating a single train of the particle focusing along the lateral centerline from 100 μL/min. These migrating patterns for the smaller particle agreed well with the results of the experiment, as seen in
Figure 5f. The case of 13 μm particle, however, shows unique migrating behavior as the flow rate increases. When increasing the flow rate up to 100 μL/min, as shown in
Figure 10b,
Fx,net becomes more consistent, moving the particle toward the lateral origin with a single equilibrium point. This is supported by our experimental observations, shown in
Figure 5e, that the intensity of the fluorescent particles focused at the lateral centerline becomes higher, up to a flow rate of 100 μL/min. Interestingly, at the flow rate of 150 μL/min, the migrating pattern in the lateral direction for 13 μm particle case remarkably changes in the Type 3 channel. As shown in
Figure 10b,
Fx,net at the maximum flow rate is characterized by two negative convexes and a positive plateau between them. This indicates two possible points of the particle focusing each at
x = 0 and
x = 55 μm, which is not observed in the case of 5 μm particles at the same flow rate. The interesting pattern of the particle focusing in Type 3 depending on the particle size could be attributed to the coupled effect between the sharp edge at the junction and the bottom wall (see
Figure 8d). As was discussed earlier, the shallow channel depth with a protruded sharp edge as in Type 3 resonates with the effect of the shear rate on its cross-plane. Thus, the Type 3 channel preferentially segregates the large particle (i.e., 13 μm) with two lateral locations for the particle focusing. In contrast, the small particle (i.e., 5 μm) consistently drifts toward a single focusing point at the lateral origin at a flow rate of 150 μL/min. This numerical result has excellent agreement with our experimental result at the same flow rate in the same microchannel. Furthermore, the current observation suggests a potential implication of the particle separation depending on its size when the aspect ratio of the T-shaped microchannel and the flow rate are appropriately leveraged.