3.1. Non-inertial focusing in curved channels of square cross-sections

Non-inertial migration of capsules in curved channels with square cross-sections is considered first. Figure 2(a) presents the numerically predicted streamwise fluid velocity component, denoted by ${U_S}$ in dimensionless form, without a capsule over the cross-section of a channel for a selected geometry. As seen, the velocity distribution in the cross-sectional plane is asymmetric with respect to the channel centreline, and shifted toward the interior face, unlike in a rectilinear channel where the distribution would be symmetric. The secondary flow is negligible in this limit, and hence the Dean vortices are not observed. The skewed velocity profile is a result of the constant streamwise pressure gradient $(1/{\tilde{R}_C})\textrm{d}\tilde{P}/\textrm{d}\varTheta $ used to drive the flow which causes a higher velocity of the fluid in the vicinity of the interior face as it moves over a shorter length compared with that near the exterior face. Note that the velocity distribution is, however, symmetric about the middle plane $(Z = 0)$ of the channel.

Figure 2.

Flow field in a curved channel in the non-inertial regime. (a) Streamwise fluid velocity contours in the cross-sectional plane for ${R_C} = 5$, $W = H = 4$. (b) Streamwise fluid velocity profile ${U_S}$ on the symmetry (middle) plane $(Z = 0)$ plotted against $R^{\prime}$. Solid red curve is the numerical result, and black dash-dot curve is the analytical result given by (3.1). Vessel centreline $(R^{\prime} = 0)$ and the location of the maximum streamwise velocity $({R^{\prime}_U})$ are indicated. Interior and exterior faces are indicated. (c) Location of the streamwise velocity maximum $({R^{\prime}_U})$ is shown as a function of channel curvature ratio $\kappa $. Here, black squares denote simulation data and red triangles denote the analytical result obtained from numerical differentiation of (3.1). Straight lines are linear fit through respective data. We vary ${R_C} = 5,7$ and 10, and $W = H = 1- 9$.

An analytical expression was obtained for the velocity field in curved microchannels in the absence of inertia by Norouzi, Kayhani & Biglari (Reference Norouzi, Kayhani and Biglari2010). From their work, the leading-order solution for the dimensionless streamwise velocity ${U_S}$ can be expressed using the current variables as

(3.1)\begin{equation}{U_S} = \sum\limits_{m = 1}^\infty {{m^{ - 3}}} [{( - 1)^m} - 1]\{ {\alpha _m}{I_1}({\zeta _m}R) + {\beta _m}{K_1}({\zeta _m}R) - {R^{ - 1}}\}\, \textrm{sin}({\zeta _m}Z + m{\rm \pi} /2),\end{equation}

where

(3.2)\begin{gather}{\alpha _m} = {{\left[ {\frac{1}{{{R_i}}}{K_1}({\zeta_m}{R_o}) - \; \frac{1}{{{R_o}}}{K_1}({\zeta_m}{R_i})} \right]} / {[{I_1}({\zeta _m}{R_i}){K_1}({\zeta _m}{R_o}) - {I_1}({\zeta _m}{R_o}){K_1}({\zeta _m}{R_i})]}},\end{gather}

(3.3)\begin{gather}{\beta _m} = {{\left[ {\frac{1}{{{R_o}}}{I_1}({\zeta_m}{R_i}) - \; \frac{1}{{{R_i}}}{I_1}({\zeta_m}{R_o})} \right]} / {[{I_1}({\zeta _m}{R_i}){K_1}({\zeta _m}{R_o}) - {I_1}({\zeta _m}{R_o}){K_1}({\zeta _m}{R_i})]}}.\end{gather}

Here, ${I_1}$ and ${K_1}$ are the modified Bessel functions of the first and second kinds, respectively, ${\zeta _m} = m{\rm \pi} /2W$, ${R_i} = {R_C} - W$ and ${R_o} = {R_C} + W$.

Figure 2(b) compares the numerically predicted velocity profile and the analytical expression given by (3.1) plotted as functions of $R^{\prime}$ on the $Z = 0$ plane. As seen, the predicted profile agrees well with the analytical expression.

Of interest is the location of the maximum streamwise fluid velocity $({U_{max}})$ which is also the location of zero shear rate. Figure 2(b) shows that this location lies between the centreline $(R^{\prime} = 0)$ and the interior face $(R^{\prime} ={-} 1)$, unlike in a rectilinear channel where it would be at the centre. Hereafter, this location is denoted as ${R^{\prime}_U}$. Of interest is the dependence of ${R^{\prime}_U}$ on the curved channel geometry. We simulate flow by varying the radius of curvature ${R_C}$ and cross-sectional dimensions $(W = H)$. Figure 2(c) shows ${R^{\prime}_U}$ as a function of the vessel curvature ratio $\kappa = \; W/{R_C}$ for different geometry. As seen, ${R^{\prime}_U}$ is negative in all simulations, indicating that it lies between the channel centreline and the interior face. We further obtain the analytical result for ${R^{\prime}_U}$ by numerical differentiation of (3.1). The analytical results show excellent agreement against the simulated data. As the curvature is increased, this location progressively shifts further toward the interior face. Furthermore, data for different curvature ratios collapse onto a straight line. This linear dependence is in qualitative agreement with the flow in a curved tube of circular cross-section at negligible inertia for which ${R^{\prime}_U}$ was analytically found as (Siggers & Waters Reference Siggers and Waters2005)

(3.4)\begin{equation}{R^{\prime}_U} ={-} \frac{{3\kappa }}{8} + {O}({\kappa ^3}),\end{equation}

where ${R^{\prime}_U}$ and $\kappa $ for the tube are similarly defined using the tube radius ${R_t}$ as $(R - {R_c})/{R_t}$ and ${R_t}/{R_C}$, respectively. In Ebrahimi et al. (Reference Ebrahimi, Balogh and Bagchi2021) where we considered curved tubes, we showed that our numerical results matched well against (3.4). Interestingly, a comparison of figure 2(c) and (3.4) suggests that, for the same curvature ratio, the velocity profile in a channel is more skewed toward the interior face than that in a tube.

The non-inertial cross-stream migration of a deformable capsule is now presented. Unless stated otherwise, the capillary number is fixed at $Ca = 1$, and the capsule is released with its centre-of-mass located in the middle plane $(Z = 0)$. Figure 3(a) shows a representative trajectory of a capsule that is initially located near the exterior face. It shows that the capsule moves in a spiralling trajectory and slowly migrates toward the interior face. Figure 3(b) shows the radial coordinate of the capsule centre in terms of $R^{\prime}$ and as a function of time for different initial locations. As seen here, for all cases, the capsule migrates toward the interior face, and eventually settles at the same radial location that is independent of the initial positions. This equilibrium location is not at the location of zero shear, rather it lies between the location of ${U_{max}}(R^{\prime} = {R^{\prime}_U})$ and the interior face $(R^{\prime} ={-} 1)$. Figure 3(b) also compares the capsule trajectory in a rectilinear channel of the same cross-section. For this, the simulation predicts that the capsule settles at the channel centreline $(R^{\prime} = 0)$ which is also the location of ${U_{max}}$, in agreement with previous studies. Therefore, the equilibrium location in the curved channel is quite different from that in a rectilinear channel, and is neither at the channel centreline nor at the location of zero shear rate. This result is significant since there is no secondary flow to drive the capsule toward the interior face in the absence of inertia.

Figure 3.

Non-inertial focusing of capsules in curved channels of square cross-sections. (a) A representative capsule trajectory for ${R_C} = 5,\; W = H = 4$ viewed along Z axis. Interior and exterior faces are shown as black circles. The capsule is released closer to the exterior face and on the symmetry plane. It migrates cross-stream toward the interior face while remaining on the symmetry plane, and eventually settles at a location close to the interior face. (b) Evolution of the radial location of the capsule centre-of-mass for different initial locations (red, black, green continuous lines). Final equilibrium location is indicated by ${R^{\prime}_{eq}}$. Also shown is the migration of a capsule in a rectilinear channel of same height and width (blue, dash-dot-dot). The inset shows $R^{\prime}$ as an angular plot for the same case shown in (b) using black line. (c) Velocity profiles in the middle plane as a function of $R^{\prime}$ for curved and rectilinear channels. Locations are marked at which capsule shapes are compared. (d–h) Capsule shapes are shown looking along the Z axis at different locations as indicated in (c). For the curved channels, orientations are adjusted so that the radial direction is uniform as indicated in (d). Three-dimensional shapes (with grey shades) correspond to that in the curved channel, and 2-D shapes showing cell perimeter correspond to that in the rectilinear channel. Note that $Ca = 1$ in all figures unless stated otherwise.

The inset of figure 3(b) shows $R^{\prime}$ as a function of the angle $\theta $ traversed by the capsule in the X–Y plane for the case of $R^{\prime}(t = 0) = 0$. When compared with $R^{\prime}$ vs. t in figure 3(b), the radial trajectory appears to be stretched as the capsule approaches the equilibrium. This is because of the higher axial velocity near the equilibrium position which makes the $R^{\prime}$ vs. t curve appears compressed in this region.

It may also be noted that since the equilibrium position is not at the interior wall, a capsule placed very close to the interior wall can migrate outward and settle at the same equilibrium location. For the specific geometry considered in figure 3(b), the closest location that the capsule can be placed to the interior wall is $R^{\prime} ={-} 0.75$, which is actually very close to the equilibrium position. We have done simulations with other geometries, for example, ${R_C} = 7,W = H = 4$, and observed that a capsule with initial position close to the interior wall indeed moves outward to the same equilibrium position.

Shape evolution of a migrating capsule $(Ca = 1)$ at different radial locations is shown figures 3(d–h). Specific radial locations selected are identified using the velocity profile shown in figure 3(c). Also compared is the shape evolution of a migrating capsule in a rectilinear channel. In both geometries, highly asymmetric shapes are observed away from the equilibrium locations. For the rectilinear channel, the capsule shape becomes symmetric about the major axes in the cross-sectional plane upon reaching the equilibrium location at the channel centreline. When projected on the $Z = 0$ plane, it shows the typical bullet shape. In contrast, for the curved channel, the shape remains highly asymmetric even after the capsule arrives at the equilibrium location (figure 3h). This increased shape asymmetry even after reaching the equilibrium is due to a threefold effect. The equilibrium position is closer to the interior wall. The streamline in this region has higher curvatures. Also the velocity profile and velocity gradient (figure 3c) both change rapidly with respect to $R^{\prime}$ around this location compared with other regions. Thus, the increased asymmetry is due to the proximity of interior wall, increased streamline curvature and rapid change in velocity and its gradient near the equilibrium location. The shape is, however, symmetric about $Z = 0$ due to the flow symmetry.

The radial migration of deformable capsules and their equilibrium in a curved channel in the absence of inertia are determined by three underlying hydrodynamic mechanisms: (i) the deformation-induced migration in the presence of a shear gradient that generally drives the capsule from a region of higher shear rate to that of a lower shear rate, (ii) the effect of the nearby wall that generally causes the capsule to move away from the wall and (iii) the streamline curvature-induced migration that drives the capsule from regions of lower to higher curvature (Shafer et al. Reference Shafer, Laiken and Zimm1974; Chan & Leal Reference Chan and Leal1979; Goh et al. Reference Goh, Phan Thien and Atkinson1984; Ghigliotti et al. Reference Ghigliotti, Rahimian, Biros and Misbah2011). For a rectilinear channel, only the first two mechanisms are present and they result in the capsule migration toward and settling at the channel centreline. For a curved channel, each of the three mechanisms causes a migration away from the exterior face in the region $R^{\prime} > 0$. In contrast, in the region $R^{\prime} < \; {R^{\prime}_U}$ the deformation and wall effects both cause a migration toward the exterior face, while the curvature effect causes a migration toward the interior face. Because of these opposing effects, the capsule settles at a location between ${R^{\prime}_U}$ and the interior face.

Effects of the radius of curvature ${R_C}$ and cross-section dimension $(W = H)$ are shown in figure 4. These parameters affect both the equilibrium position and the rate of migration. Figure 4(a) shows the effect of ${R_C}$ while $W = H$ is held fixed. As seen, as ${R_C}$ decreases, the capsule migrates faster, and the equilibrium position moves further closer toward the interior face. Figure 4(b) shows the effect of $W = H$ while ${R_C}$ is fixed. It shows that, as $W = H$ increases, the migration rate becomes slower, but the equilibrium position moves closer to the interior face. These trends are due to the occurrence of a higher streamline curvature near the interior face with increasing ${R_C}$ and W, resulting in the dominance of the streamline curvature effect which causes a migration toward and settling further closer to the interior face. Note that the capsules take different times to reach the equilibrium positions depending on ${R_C}$ and W. As ${R_C}$ increases, the curvature-induced migration weakens, and capsules take longer to reach equilibrium. In the limit that ${R_C} \to \infty $, which corresponds to a straight channel, the curvature-induced migration is absent.

Figure 4.

Influence of radius of curvature ${R_C}$ and half-width W on the non-inertial focusing in curved channels of square cross-sections. (a) Radial trajectories are shown in terms of $R^{\prime}$ versus time for varying ${R_C} = 4,5,7,10$, while $W = H = 3$ is fixed. (b) Radial trajectories for varying $W = H = 3,4,6$ are shown while ${R_C} = 7$ is kept fixed. (c) Equilibrium locations are shown in terms of $({R_{eq}} - {R_C})$ (solid lines) as functions of W for ${R_C} = 5$ (red, square), 7 (blue, triangle), 10 (green, circle). Also shown is the location of the streamwise velocity maximum as $({R_U} - {R_C})$ (dash lines). (d) Values of ${R^{\prime}_{eq}}$ (solid lines) and ${R^{\prime}_U}$ (dash lines) are shown versus curvature ratio $\kappa $.

Multiple simulations are performed by varying ${R_C}$ and $W = H$, and the resulting equilibrium locations are presented in figure 4(c,d). Hereafter, the equilibrium radial location is denoted as ${R_{eq}}$. Figure 4(c) shows ${R_{eq}} - {R_C}$, which indicates how far a capsule settles from the channel centreline, as a function of W. Also shown for comparison is ${R_U} - {R_C}$, the distance of the location of ${U_{max}}$ from the channel centreline. As seen here, for all cases considered, ${R_{eq}} - {R_C} < 0$, and $({R_{eq}} - {R_C}) < ({R_U} - {R_C})$. Therefore, irrespective of channel dimensions, the capsule equilibrium location is further closer to the interior face than the location of maximum streamwise velocity. The results further show that the difference between $({R_{eq}} - {R_C})$ and $({R_{eq}} - {R_U})$ increases with increasing W and decreasing ${R_C}$, meaning that the capsule settles increasingly further away from the location of zero shear. Figure 4(d) shows ${R^{\prime}_{eq}} = ({R_{eq}} - {R_C})/W$ as a function of the curvature ratio $\kappa $. Similar to the linear dependence of ${R^{\prime}_U}$ on $\kappa $ as observed earlier in figure 2, a linear variation of ${R^{\prime}_{eq}}$ with respect to $\kappa $ is also predicted. Note that ${R^{\prime}_{eq}} \to 0$ as $\kappa \to 0$, which is the case of a rectilinear channel, and can be seen in the trend of the numerical results. Also, the slope of ${R^{\prime}_{eq}}$ versus $\kappa $ is higher than that of ${R^{\prime}_U}$ versus $\kappa $ implying that the capsule settles further away from the location of zero shear and toward the interior face with increasing channel curvature and width.

The results presented so far are for a constant capillary number of $Ca = 1$. The influence of $Ca$ is considered next. Figure 5(a) presents the radial location in terms of $R^{\prime}$ against time for different values of $Ca$ in the range 0.2 to 10. As seen, the equilibrium location is nearly independent of $Ca$. The rate of migration, however, increases with increasing $Ca$ as it is generally the case in a unidirectional shear flow. The figure also shows the capsule shapes near the equilibrium for different Ca. Although the equilibrium position is independent of Ca, the shapes are not. This is because the shapes are directly affected by $Ca$ – the higher $Ca$, the greater the deformation. In contrast, the equilibrium position is dictated by several competing mechanisms, i.e. shear gradient, and the curvature-induced migration. Since these two mechanisms oppose each other in the interior region of the channel, the result is an equilibrium position that is independent of $Ca$.

Figure 5.

(a) Effect of $Ca$ in the non-inertial regime. Shown here is the radial position as a function of time for different $Ca$. Capsules migrate to the same radial location that is nearly independent of $Ca$. Here, ${R_C} = 5,\; W = 4$. Also shown are the capsule shapes near the equilibrium position for different $Ca$. (b) Effect of various initial $Z$-locations. Capsules are released at different initial Z locations that are marked as ‘O’. Trajectories are shown in the cross-sectional plane. Irrespective of initial Z locations, capsules migrate to the$\; Z = 0$ plane and settle at the same radial location indicated as ${R_{eq}}$. Here ${R_C} = 7,\; W = 6$. (c) Capsule trajectories in the cross-section plane in a rectilinear channel $(W = H = 4)$.

We now provide some scaling arguments based on previous and current results. Previous studies have considered mathematical theory of cross-stream migration of deformable particles in unidirectional shear flow in the non-inertial regime (e.g. Chan & Leal Reference Chan and Leal1979; Helmy & Barthes-Biesel Reference Helmy and Barthes-Biesel1982; Shapira & Haber Reference Shapira and Haber1990). We focus on the cross-stream velocity of the capsule, also termed the migration or lateral velocity, denoted by ${V_m}$; the equilibrium location corresponds to ${V_m} = 0$. Note that ${V_m}$ is scaled by ${U_0}$ as before. In the limit of small deformation and in an unbounded 2-D Poiseuille flow, the migration velocity of a liquid drop (Chan & Leal Reference Chan and Leal1979) and a capsule (Helmy & Barthes-Biesel Reference Helmy and Barthes-Biesel1982) follows ${V_m}\,\sim \mathcal{E}(\tilde{\eta }/\tilde{W}){({a_0}/\tilde{W})^3}$, where $\mathcal{E}$ is a dimensionless parameter that involves properties of the particle that determine its deformability, e.g. surface tension or membrane elasticity, and internal to external fluid viscosity ratio, $\tilde{W}$ is a characteristic length of the flow and $\tilde{\eta }$ is the distance of the particle centre from the location of ${U_0}$. Generally, $\mathcal{E}\sim D$, where D is a measure of particle deformation, e.g. the Taylor deformation parameter. In the limit of small deformation, $D\sim Ca$. Note that these studies did not include the direct wall interaction. Later, a numerical study by Doddi & Bagchi (Reference Doddi and Bagchi2008) showed that for capsules in a 2-D channel of half-height $\tilde{W}$, ${V_m}\sim {F_1}(Ca){F_2}({\tilde{\xi }_0}/\tilde{W}){({a_0}/\tilde{W})^3}$, where ${F_1}$ and ${F_2}$ are functions of $Ca$ and ${\tilde{\xi }_0}/\tilde{W}$, and ${\tilde{\xi }_0}$ is the distance of the capsule centre from the nearest wall. Additional studies (e.g. Coupier et al. Reference Coupier, Kaoui, Podgorski and Misbah2008; Losserand et al. Reference Losserand, Coupier and Podgorski2019) showed that in a general shear flow (linear or parabolic) ${V_m}\sim \mathcal{E}{({a_0}/{\tilde{\xi }_0})^\delta }{a_0}\tilde{G}(\tilde{\xi })/{U_0}$, where $\tilde{G}(\tilde{\xi })$ is the shear rate, $\tilde{\xi }$ is the coordinate in the cross-flow direction and the exponent $\delta $ is a constant. In the case of linear shear flow, $\delta = 2$. This form shows that the wall repulsion has a stronger effect near the wall. From this general form, one can write for a 2-D Poiseuille flow in a channel ${V_m}\sim \mathcal{E}(a_0^{\delta + 1}/\tilde{\xi }_0^\delta \tilde{W})G$ where $G(\tilde{\xi }) = \tilde{G}/{\tilde{G}_{max}}$ is the dimensionless shear rate and ${\tilde{G}_{max}}$ is the maximum (wall) shear rate. Following Chan & Leal Reference Chan and Leal1979, Helmy & Barthes-Biesel Reference Helmy and Barthes-Biesel1982, Shapira & Haber Reference Shapira and Haber1990 and Doddi & bagchi Reference Doddi and Bagchi2008 one may let $\delta = 2$, which gives ${V_m}\sim \mathcal{E}(1/{W^3}){(W/\xi )^2}G$; here, all lengths are made dimensionless using ${a_0}$. This expression gives the relative strengths of different mechanisms causing migration of a deformable particle in a wall-bounded parabolic flow: $\mathcal{E}$ for the deformation effect, G for the velocity gradient effect, ${(W/\xi )^2}$ for the wall repulsion.

As evident from above, there has been a significant mathematical investigation on the migration velocity of deformable particles in unidirectional flow. But unfortunately there exists no such study for curved vessels. Here, we present a scaling argument for the migration velocity in curved flows borrowing from the above mathematical developments and using our present numerical data. In a curved channel, there are migration due to (i) particle deformation, (ii) velocity gradient, (iii) wall repulsion and (iv) curvature effect. The first three can be taken care of by the above formulation, and we denote it as $V_m^{\textrm{(i)}-\textrm{(iii)}}\sim \mathcal{E}(1/{W^3}){(W/\xi )^2}G$. As the present numerical data suggest, the streamline curvature can induce a migration of a deformable particle even if the flow is unbounded. We propose this migration velocity arising purely from the curvature effect can be modelled by replacing the terms representing velocity gradient and wall repulsion effect by a term that is dependent on the curvature as $V_m^{(\textrm{iv})}\sim \mathcal{E}(1/{W^3}){\kappa ^m}$, where the exponent m is a constant. We now verify this expression using our numerical simulation data. Figure 6(a) shows ${V_m}$ as a function of $\kappa $ for fixed values of $Ca$ and W. Our data suggest a nonlinear dependence, which can be approximated by taking $m = 2$. Then, in figure 6(b) we plot ${V_m}$ as a function of $(1/{W^3}){\kappa ^2}$, which gives an apparent linear scaling. Finally, figure 6(c) shows ${V_m}$ as a function of $Ca$. A nonlinear variation is observed since the deformation (and, $Ca$) is not small. Note that we have assumed a similar dependence of $V_m^{\textrm{(i)}-\textrm{(iii)}}$ and $V_m^{(\textrm{iv})}$ on $\mathcal{E}$. This is based on the present finding of figure 5(a) that the equilibrium locations are independent of $Ca$. In the limit of small deformation in a unidirectional shear flow, the migration velocity varies linearly with respect to $Ca$ (Chan & Leal Reference Chan and Leal1979; Helmy & Barthes-Biesel Reference Helmy and Barthes-Biesel1982). Since the curvature-induced migration must balance the shear-induced migration, the former must also be linear with respect to $Ca$ in this limit. Furthermore, since the equilibrium location is independent of Ca, both the curvature-induced migration and shear-induced migration must have the same scaling with respect to Ca. The net migration velocity is then

(3.5)\begin{equation}{V_m}\sim{\pm} \mathcal{E}\left( {\frac{1}{{{W^3}}}} \right){\left( {\frac{W}{\xi }} \right)^2}G + \mathcal{E}\left( {\frac{1}{{{W^3}}}} \right){\kappa ^2},\end{equation}

where ${\pm} $ is used for the first term to indicate that the migration caused by the velocity gradient and wall repulsion can be either inward or outward, while that due to the curvature is always inward. For $R^{\prime} > \; {R^{\prime}_U}$, both terms have same sign, so that the net effect is additive and toward the interior wall. For $R^{\prime} < \; {R^{\prime}_U}$, they are opposite, resulting in the particle equilibrium.

Figure 6.

Migration velocity near the curved channel centre as a function of (a) $\kappa $, (b) ${\kappa ^2}/{W^3}$ and (c) $Ca$.

It may also be noted that it is not possible to isolate the effect of the four mechanisms (i)–(iv) as noted above through our simulations. This is because a curved flow (bounded or unbounded) will always have a velocity gradient. In contrast, a curved flow representing a rigid body-like rotational fluid motion will not induce any deformation, and hence there is no migration. Furthermore, the cross-stream migration is absent as $Ca \to 0$.

The results presented so far are for capsules released at the middle plane $(Z = 0)$ of the channel. Since the velocity is symmetric about this plane, there cannot be any out-of-plane motion, and the capsule remains in this plane while it radially migrates toward the interior face. The influence of different initial Z locations is shown in figure 5(b) where capsule trajectories in the cross-sectional plane are shown. For all cases, capsules migrate to the $Z = 0$ plane and settle on this plane at the same equilibrium radial location ${R_{eq}}$. Since there is no curvature effect in the Z direction, the capsule motion in Z is caused by the wall effect and the shear gradient-induced deformation effect. Both of them drive the capsule toward the $Z = 0$ plane (which is the location of the lowest shear rate with respect to $Z$), while the curvature effect along R drives it to ${R_{eq}}$. Furthermore, because of the curvature effect, the capsule migrates faster in the radial direction and arrives at the equilibrium radial location earlier than it arrives at the $Z = 0$ plane. Also shown in figure 5(c) are capsule trajectories in the cross-section plane in a rectilinear channel. For this geometry, all trajectories converge to the channel centre, implying that the equilibrium position is the channel centreline in a rectilinear channel irrespective of initial positions, but it moves to a radially off-centred location toward the interior face with the introduction of the channel curvature.

The present results for capsule migration and equilibrium locations in curved channels of square cross-sections are qualitatively similar to those of curved tubes of circular cross-sections. As predicted by Ebrahimi et al. (Reference Ebrahimi, Balogh and Bagchi2021), for a curved tube also, capsules migrate radially toward the inner side of the tube and settle at an equilibrium location that is located increasingly away from the location of the maximum streamwise velocity and closer to the inner side with increasing curvature ratio. As seen here for curved channels, the capsule equilibrium position for curved tubes was also predicted to be independent of $Ca$, and located on the $Z = 0$ plane, which is the plane of symmetry for the tube. Figure 7 shows a quantitative comparison of the equilibrium locations for curved tubes and channels. Despite the difference in geometry, the predicted equilibrium locations are surprisingly very close in the two cases.

Figure 7.

Comparison of equilibrium positions for curved channels of square cross-sections (red, solid line and squares) and curved tubes of circular cross-section (black, dash line and circles) in terms of (a) ${R^{\prime}_{eq}}$ versus $\kappa $, and (b) ${R_{eq}} - {R_C}$ versus ${R_U} - {R_C}$. The data for the curved tubes are taken from Ebrahimi et al. (Reference Ebrahimi, Balogh and Bagchi2021).

3.2. Non-inertial focusing in curved channels of rectangular cross-sections

In this section, curved channels of rectangular cross-sections are considered in the non-inertial regime. For this purpose, we introduce the aspect ratio $\gamma = W/H$. Other relevant geometric parameter is the curvature ratio $\kappa = W/{R_C}$ as defined before. The influence of varying $\gamma $ while keeping $\kappa $ constant is considered first. Figure 8 shows the velocity distribution and the location of the maximum velocity. The velocity profile becomes more skewed toward the interior face with increasing $\gamma $. The trend of ${R^{\prime}_U}$, however, is more complex. In general, ${R^{\prime}_U}$ moves closer to the interior wall as $\gamma $ increases (or H decreases) due to the confinement effect. In contrast, ${R^{\prime}_U}$ becomes independent of H (and $\gamma $) as channel height increases (or aspect ratio decreases). Variation of ${R^{\prime}_U}$ with respect to H (and $\gamma $) becomes greater as the channel width decreases.

Figure 8.

Effect of aspect ratio $\gamma = W/H$ while keeping the curvature ratio $\kappa = W/{R_C}$ constant on the velocity distribution in the non-inertial regime. Velocity contours are shown in the cross-sectional plane for (a) $\gamma = 1.5$ and (b) $1/3$, while $\kappa = 3/5$. (c) Streamwise velocity profile in $Z = 0$ plane as a function of $R^{\prime}$ for $\gamma = 1.5$ (green, dash-dot curve), 0.5 (black, continuous) and 1/3 (red, dash); $= 3/5,$ ${R_C} = 5$. (d) Location of the streamwise velocity maximum ${R^{\prime}_U}$ as a function of channel height. Here ${R_C} = 7,\; W = 3$ (black, continuous curve), ${R_C} = 7,\; W = 6$ (red, dash-dot), ${R_C} = 5,\; W = 3$ (green, dash). Inset shows the same data plotted against $\gamma $.

Radial trajectories of capsules for different aspect ratio are shown in figure 9(a). For all cases considered in the figure, capsules migrate toward the interior face. Noticeably, differences exist between capsule trajectories at larger aspect ratio (i.e. smaller channel height); but as $\gamma $ decreases, trajectories become indistinguishable. Equilibrium radial location ${R^{\prime}_{eq}}$ is shown in figure 9(b) as a function of H, and $\gamma $ (inset). Also presented in the figure the location of the streamwise velocity maximum that can be compared with the capsule equilibrium location. The equilibrium is between the location of velocity maximum and the interior face of the channel. The figure further shows that for smaller channel heights (or higher $\gamma $) the equilibrium location progressively moves closer to the interior face with decreasing channel height (or increasing $\gamma $). However, at larger values of H (i.e. as $\gamma \to 0$), the equilibrium locations asymptote to fixed values that are independent of H or $\gamma $. This trend of capsule trajectory and the equilibrium location follows the trends of ${R^{\prime}_U}$ and streamwise velocity profiles which become insensitive to H at larger values, as seen previously in figure 8(d).

Figure 9.

Effect of varying the aspect ratio $\gamma = W/H$ while keeping the curvature ratio $\kappa = W/{R_C}$ constant on the non-inertial migration. Capsule motion in $Z = 0$ plane is considered. (a) Radial trajectory in terms of $R^{\prime}$ as a function of time is shown for different $\gamma $. (b) Capsule equilibrium location ${R^{\prime}_{eq}}$ (filled symbols) as a function of channel height. Also shown for comparison is the location of streamwise velocity maximum ${R^{\prime}_U}$ (open symbols). Here ${R_C} = 7,W = 3$ (black, continuous curve), ${R_C} = 7,\; W = 6$ (red, dash-dot),${R_C} = 5,\; W = 3$ (green, dash). The inset shows the same data plotted against. Panels (c) and (e) show capsule shapes for ${R_C} = 5,\; W = 3,\; H = 2$ ($\gamma = 1.5,\; \kappa = 3/5$) near the channel centreline and equilibrium, respectively, and panels (d) and (f) show capsule shapes for ${R_C} = 5,\; W = 3,\; H = 9$ ($\gamma = 1/3,\; \kappa = 3/5$) near the channel centreline and equilibrium, respectively. (g) Time for capsules released at the channel centreline to reach 95 % of the equilibrium distance as a function of H. ${R_C} = 7,\; W = 3$ (black, continuous curve), ${R_C} = 7,\; W = 6$ (red, dash-dot), ${R_C} = 5,\; W = 3$ (green, dash). The inset shows the same data as a function of $\gamma $.

It may be noted in figure 9(b), that the curve for ${R^{\prime}_U}$ for ${R_C} = 7;\;W = 6$ does not asymptote with increasing H like other curves. This is because a higher value of W requires a higher value of H to reach an asymptote. Other cases considered in the figure are for $W = 3$, for which the asymptote is reached for $H\; < \; 8$. For $W = 6$; a much higher value of H would be needed to see an asymptotic trend. Furthermore, for this specific case, ${R^{\prime}_{eq}}$ reaches an asymptote before ${R^{\prime}_U}$ does, unlike other cases. This is because the increased curvature ratio $\kappa $ at this W causes a stronger curvature-induced inward migration of the capsule. In this case, the curvature-induced migration is primarily balanced by the wall repulsion, while the shear gradient has a weaker effect. Since confinement in H affects weakly and indirectly via the alteration of the velocity profile, the balance of these two mechanisms (which remain fixed due to the fixing of W and $\kappa $) results in the noted behaviour.

Capsule shapes are presented in figure 9(c–f). Away from the equilibrium locations, capsules attain crescent shapes irrespective of channel heights. Close to the equilibrium locations, the shape resembles an asymmetric parachute for smaller channel heights, but a bullet shape for larger heights. This difference in shape is due to the equilibrium location moving closer to the interior face (and, hence, toward the region of a higher streamline curvature) and increased velocity profile asymmetry with increasing aspect ratio.

Figure 9(a) further suggests that the rate of radial migration is dependent on $\gamma $, even when migration along $Z = 0$ plane is considered and the curvature ratio is held fixed. Although $Z = 0$ plane is the plane of symmetry, and the channel height is not expected to make a significant effect on the radial migration along this plane, it is evident from figure 9(a) that capsules migrate faster toward the equilibrium position with increasing $\gamma $. This trend may also seem counter-intuitive as increasing confinement is expected to retard the capsule motion. However, the higher migration rate is due to the greater asymmetry in capsule shape, velocity profile and the equilibrium location being further closer to the interior face and in a region of higher streamline curvature, all of which result from the increasing aspect ratio as discussed above. Of interest, then, is the time ${T_{0.95}}$ taken by a capsule to travel 95 % of the radial distance from the channel centreline to the equilibrium location, which is presented in figure 9(g). As seen, ${T_{0.95}}$ becomes shorter as the channel height decreases, but reaches an asymptote at larger heights. The variation of ${T_{0.95}}\; $with respect to H (or $\gamma $) is most prominent for channels of larger widths,because of the higher streamline curvature in such cases. This result is intriguing because although the equilibrium location moves further away from the channel centreline and toward the interior face with decreasing H, the capsule takes a shorter time to reach the equilibrium. For practical microfluidic applications in the non-inertial regime, this result is important since it implies that a rapid focusing of deformable particles, even when the motion is restricted in the midplane, can be achieved with smaller channel heights (or larger aspect ratio) for a given width.

Additional study for rectangular cross-section is done by varying W while keeping H constant. These two dimensionless parameters also represent the blockage by the capsule in the radial and normal directions, respectively, which differently affect capsule migration. The curvature-induced migration is active only in the radial direction, while migration in the normal direction is primarily dictated by the shear gradient. Where possible, we present the data using $\gamma $ and $\kappa $. The velocity field in the absence of the capsule is first described. Figures 10(a)–10(c) show the results for the streamwise velocity distribution as $\kappa $ is varied. As seen, increasing $\kappa $ causes the velocity profile to be more skewed toward the interior face. This is because increasing the width at a fixed radius of curvature causes a higher streamline curvature near the interior face. Also, the maximum axial velocity increases with increasing width. The numerical velocity profile is compared against the analytical expression ((3.1)) for a rectangular cross-section in figure 10(c) which shows that they agree very well. The location ${R^{\prime}_U}$ of the maximum axial velocity ${U_{max}}$ is shown in figure 10(d) as a function of $\kappa $ for different H. As seen, ${R^{\prime}_U}$ moves closer to the interior face with increasing $\kappa $. In contrast, ${R^{\prime}_U}$ moves away from the interior face with increasing H. Unlike for curved channels of square cross-sections, ${R^{\prime}_U}$ for rectangular cross-sections does not show a linear dependence with respect to $\kappa $. A linear dependence is not observed even when the data are plotted by redefining the curvature ratio using the hydraulic radius as ${R_h}/{R_C}$. Figure 10(e) shows ${R^{\prime}_U}$ as a function of $\kappa $ for a given $\gamma $. Data for different $\gamma $ do not collapse onto a single line.

Figure 10.

Non-inertial flow in curved channels of rectangular cross-sections. (a,b) Streamwise velocity contours for varying channel width at constant height: (a) ${R_C} = 7,\; H = 3,\; W = 6$, ($\kappa = 6/7$, $\gamma = 2$); (b) ${R_C} = 7,\; H = 3,\; W = 2$. (c) Numerical velocity profiles at $Z = 0$ plane as a function of $R^{\prime}$ for ${R_C} = 7,\; H = 3,\; W = 6$ ($\kappa = 6/7$, $\gamma = 2$, solid black curve), ${R_C} = 7,\; H = 3,\; W = 4$ ($\kappa = 4/7$, $\gamma = 4/3$, dash-dot-dot), ${R_C} = 7,\; H = 3,\; W = 2$ ($\kappa = 3/7$, $\gamma = 2/3$, dash-dot). Also shown is the analytical velocity given by (3.1) for ${R_C} = 7,\; H = 3,\; W = 6$ ($\kappa = 6/7$, $\gamma = 2$, dash, red curve). (d) Location of the maximum streamwise velocity as a function of $\kappa $ for different $H$: ${R_C} = 7,\; H = 3$ (solid black), ${R_C} = 7,\; H = 4$ (green, dash), ${R_C} = 7,\; H = 6$ (red, dash-dot-dot), ${R_C} = 5,\; H = 4$ (blue, dash-dot). (e) Value of ${R^{\prime}_U}$ as a function of $\kappa $ for different $\gamma $.

Capsule trajectory in terms of $R^{\prime}$ versus time is shown in figure 11(a) for varying width. Capsules are initially released at $Z = 0$ plane, since any other Z location would cause migration toward $Z = 0$ as observed before for square cross-sections. For all cases seen in the figure, radial migration toward the interior face and settling at an equilibrium location are predicted. For fixed ${R_C}$ and H, capsules settle increasingly closer to the interior face with increasing W. The equilibrium locations are presented in figure 11(b) in terms of ${R_{eq}} - {R_C}$. As seen, for fixed ${R_C}$ and H, values of ${R_{eq}}$ decrease with increasing W implying that capsules settle increasingly away from the channel centreline. Also indicated in the figure are the equilibrium locations in channels of square cross-sections. The settling distance from the channel centreline for rectangular cross-sections could be larger or smaller than that for square cross-sections depending on the channel width. Figure 11(c) presents ${R^{\prime}_{eq}}$ as a function of $\kappa $. The inset of figure 11(c) shows ${R^{\prime}_{eq}}$ and ${R^{\prime}_U}$ against $\kappa $ for a given aspect ratio. As seen here, ${R^{\prime}_{eq}}$ lies below ${R^{\prime}_U}$. Therefore, similar to the results of channels with square cross-sections, here also capsules settle at locations between the maximum streamwise velocity and the interior face, where the shear rate is non-zero. ${R^{\prime}_{eq}}$ versus $\kappa $ shows a linear dependence.

Figure 11.

Non-inertial migration of capsules in curved channels of rectangular cross-sections. Migration in the $Z = 0$ plane is considered. (a) Radial coordinate of a migrating capsule is shown in terms of $R^{\prime}$ as a function of time for varying channel width as $W = 4$ (dash-dot), 5 (dash) and 6 (continuous curve) for a fixed ${R_C} = 7,\; H = 6$. (b) Equilibrium location from the channel centreline shown as ${R_{eq}} - {R_C}$ as a function of W for various ${R_C}$ and H. Open circles are used to indicate data from square cross-sections. (c) Equilibrium location ${R^{\prime}_{eq}}$ (filled symbols) presented as a function of $\kappa $ and for various ${R_C}$ and H. For comparison, locations of the streamwise velocity maximum ${R^{\prime}_U}$ are also shown using open symbols. ${R_C} = 7,\; H = 3$ (black, continuous line), ${R_C} = 7,\; H = 4$ (green, dash), ${R_C} = 7,\; H = 6$ (red, dash-dot-dot), ${R_C} = 5,\; H = 4$ (blue, dash-dot). The inset shows the data plotted against $\kappa $ for a given aspect ratio. (d–f) Shapes of a migrating capsule when located at the centreline of the channel, and (g–i) at the equilibrium position; $W = 2$ (d,g), 4 (e,h), 6 (f,i); ${R_C} = 7,\; H = 3$ are fixed.

Also, figure 11(c) shows that when ${R^{\prime}_{eq}}$ and ${R^{\prime}_U}$ are plotted for a given H, ${R^{\prime}_U}$ shows a greater difference between different channel heights for a given$\; \kappa $ than ${R^{\prime}_{eq}}$. This difference becomes of the similar magnitude when the data are plotted for a given aspect ratio instead of H as shown in the inset. The difference arises because reducing the channel height makes the velocity profile more skewed toward the interior wall; but the capsule migration in the radial direction is governed by the relative strengths of the curvature, shear gradient and wall effects acting in the radial direction only; confinement in H has only indirect effect via the alteration of the radial velocity profile.

The mechanisms that lead to the settling are similar to those discussed before. Since increasing the channel width results in a higher streamline curvature as the interior face is approached, a stronger curvature-induced inward migration ensues. In contrast, the radial extent over which the deformation-induced outward migration happens is shortened because of the increasing proximity of ${R^{\prime}_U}$ to the interior face. Their combined effect leads to the capsule settling increasingly closer to the interior face as W increases.

The influence of channel width on the capsule shape evolution is considered next. Figures 11(d)–11(f) shows the shapes when the capsules are located near the centreline of the channel. An asymmetric parachute shape is seen for the smallest width, but a crescent shape is seen for the largest width considered. As the equilibrium location is approached (figure 11g–i), the shape transitions to an asymmetric bullet shape for channels with a smaller width, but to an asymmetric parachute shape for a larger width. A greater degree of asymmetry at a larger channel width is due to the increased streamline curvature neat ${R_{eq}}$, a higher asymmetry in the fluid velocity profile, and the closer proximity of ${R_{eq}}$ to the interior face.

A few comments can be made now about the confinement effect. For a highly confined channel, the presence of the capsule is expected affect the flow. The range of W considered here is $2$ to $6$. The fluid velocity distribution in the channel in presence of the capsule is shown in in the supplementary material are available at https://doi.org/10.1017/jfm.2021.868 for $W = 2$ and 6. Even for the most confined case of $W = 2$, the velocity is altered in a small region around the capsule, and remains unchanged everywhere else. We do expect this, since even for this most confined case, the capsule to channel volume fraction is approximately 0.004 (${R_C} = 7,\; H = 3$), implying a dilute flow. We have also performed simulations at higher $W = 10{a_0}$, and the same trend is observed, that is, the capsules attain an equilibrium near the interior wall. The difference that arises for different W is the rate of radial migration and the actual radial location of the equilibria as discussed above. Even if the flow is unconfined (i.e. $W/{a_0} \to \infty $), the 2-D simulation of Ghigliotti et al. (Reference Ghigliotti, Rahimian, Biros and Misbah2011) showed that capsules migrate toward the interior region of the flow, in agreement with the present results. The confinement indeed plays a role, since the capsule attains the equilibrium position because of the presence of the interior wall, which imparts a repulsion from the wall as well as the velocity gradient that provides a shear-induced migration away from the inner wall. On the contrary, in an unconfined flow with curved streamlines, the capsule would continue to migrate toward the interior region without attaining an equilibrium, because there is no wall repulsion and opposing velocity gradient.

3.3. Inertial focusing

We now consider capsule migration in the presence of inertia by taking the capsule Reynolds number $Re_{a}\gtrsim 0.1$; accordingly, the channel Reynolds number $Re_{C}\gtrsim 1$. The results will be discussed in terms of $Re_{C}$; the corresponding Dean number can be obtained as $De = Re_{C}\sqrt \kappa $ . The velocity field for square cross-sections is discussed first, followed by capsule migration in such geometry. Then, migration in channels of rectangular cross-sections will be considered. Unless stated otherwise, $Ca = 1$ is considered.

In the presence of inertia, a secondary flow develops that is directed from the interior to the exterior face along the middle plane, but in the opposite direction along the top and bottom walls, generating two counter-rotating Dean vortices, as shown in figure 12(a). The strength of the secondary flow increases with increasing inertia, as shown in figure 12(b) where the maximum of the secondary velocity ${U_D}$ obtained at the middle plane is presented as a function of $Re_{C}$ and $De$. The predicted secondary velocity in the range $De\lesssim 10$ agrees with the theoretical result that, in the leading order, ${U_D}\sim {O}(De)$ (Siggers & Waters Reference Siggers and Waters2005). Figure 12(c) shows the streamwise velocity profiles for different $Re_{C}$. For $Re_{C}\lesssim 10$ velocity profiles are nearly indistinguishable from that at negligible inertia and are skewed toward the interior face of the channel. For $Re_{C} \gtrsim 10$, the profiles become blunt, and the location of the maximum streamwise velocity moves away from the interior wall.

Figure 12.

(a) Secondary flow at finite inertia. Colour contours show the magnitude of secondary velocity ${U_D}$ in a cross-sectional plane for $Re_{C} = 40$, ${R_C} = 7,\; W = H = 4$. (b) Maximum of the secondary velocity, scaled by the maximum streamwise velocity, as a function of $Re_{C}$, and $De$, for ${R_C} = 5,\; W = H = 4$. (c) Streamwise velocity profiles scaled by the maximum for each case on $Z = 0$ plane for different $Re_{C}$. Here ${R_C} = 5,\; W = H = 4$. Corresponding Dean number $De = \sqrt {0.8} \; Re_{C}$.

Inertial migration of capsules released at the channel centreline and on the midplane $(R = {R_C},Z = 0)$ is shown in figure 13(a). Here $Ca = 1$, and ${R_C} = 7,\; W = H = 4$ are kept fixed, but $Re_{C}$ is varied over two orders. It is seen that the direction of capsule migration can be either towards the interior or exterior face. Accordingly, the equilibrium position at finite inertia is dependent on $Re_{C}$. In general, capsules settle closer to the interior face at smaller $Re_{C}$, while they settle further away at higher values. For the specific geometry considered in figure 13(a), capsules are observed to migrate toward the interior face for $Re_{C} < 2$. Within this range, they settle at different radial locations that lie increasingly away from the interior face as $Re_{C}$ increases. For $Re_{C}\gtrsim 2$, capsules migrate toward the exterior face. For $Re_{C}\gtrsim 8$, capsules first migrate closer to the exterior face, but do not settle there. Figure 13(b) shows that for this range of $Re_{C}$, capsules leave the middle plane $(Z = 0)$ after initially staying on this plane for an extended time. Figure 13(c) shows the capsule trajectory in the cross-sectional plane superimposed with the streamlines of the secondary flow. As seen, the capsule remains in the midplane for some time, and then leaves this plane to follow a spiralling trajectory in the same direction as the Dean vortex. Eventually, it settles near the centre of one of the vortices. Figure 13(d) shows the cross-sectional trajectories for capsules released at different locations with respect to the middle plane. In all cases, the final equilibrium location is the same and near the centre of the vortex. Evidently, such trends of capsule migration at finite inertia are significantly different from those seen before in the absence of inertia.

Figure 13.

Inertial focusing in curved channels of square cross-sections. (a) Radial trajectories are shown for varying channel Reynolds number as $Re_{C} = 0.08,\; 0.8,2,8,80$, for ${R_C} = 7,W = H = 4$. (b) The Z coordinate of capsule's centre-of-mass is shown for $Re_{C} = 8$ and 80. (c) Trajectory on the cross-sectional plane is shown for $Re_{C} = 80$ superimposed on streamlines of the secondary flow. (d) Trajectories on the cross-sectional plane for different initial release locations are shown, all of which merge to the same final equilibrium location near the centre of the Dean vortex.

Figure 14(a) shows capsule trajectories at different $Re_{C}$ for $Re_{C} \ge 8$, for which migration to the vortex centre is observed. The location of the vortex centre is also shown in the figure, which moves closer toward the exterior face and further away from the midplane with increasing $Re_{C}$. Accordingly, capsules settle further toward the exterior face and away from the middle plane with increasing $Re_{C}$.

Figure 14.

Migration of a capsule toward the vortex centre. (a) Influence of $Re_{C}$ on capsule trajectory and equilibrium location in comparison with the location of the vortex centre is shown for ${R_C} = 5,W = H = 4$. Filled circles indicate the vortex centre for $Re_{C} = 8$ (black), 40 (green), 80 (blue), 320 (red). (b–e) Capsule shapes near respective equilibrium locations for different $Re_{C}$; (b) $Re_{C} = 0.8$, (c) $Re_{C} = 2$, (d) $Re_{C} = 8$ and (e) $Re_{C} = 80$.

What causes such trends in the capsule migration at finite inertia is the presence of the secondary flow, in addition to the effects from capsule deformation, altered shear gradient and the presence of the streamline curvature. As noted before, the deformation effect causes a migration toward a region of low shear, the streamline curvature causes a migration toward the interior wall, and the wall effect causes migration away from the wall. It can be inferred from the last section (§ 3.1) that these three mechanisms result in a capsule migration velocity ${V_m}$ that is of the order of $V_{m}\lesssim \ {O}({10^{ - 2}})$ and toward the interior face. At finite inertia, the shear gradient effect causes a migration away from the zero shear location, just as it does for a rigid particle. In addition, the secondary flow results in a ‘drag’, termed Dean's drag, and drives the capsule toward the exterior wall and to follow the spiralling trajectory toward the vortex centre. The final equilibrium location then depends on the strengths of these competing mechanisms, which can be understood in terms of their scaling as follows. Once again, however, we are faced with the problem of how to separate these different mechanisms in a confined curved flow, since a curved flow (except the pure rotational flow) will always have a velocity gradient. To isolate the inertial lift, we perform simulations of capsule migration in the unidirectional flow in rectilinear channels at varying $Re_{C}$ and $Ca$ so that no curvature effect is present. These results are shown in figure 15 and can be viewed against the theoretical work. Ho & Leal (Reference Ho and Leal1974) derived the cross-stream migration velocity of a rigid sphere in 2-D Poiseuille flow in a channel of width $2\tilde{W}$ as

(3.6)\begin{equation}{V_m}\sim Re_{a}{\left( {\frac{a}{{\tilde{W}}}} \right)^2}{(1 - \hat{\xi })^2}{G_1}(\hat{\xi }) - (1 - \hat{\xi })\; {G_2}(\hat{\xi }),\end{equation}

where $\hat{\xi } = \tilde{\xi }/\tilde{W}$, $\tilde{\xi }$ is the distance of the sphere centre from the nearest wall, and ${G_1}$ and ${G_2}$ are functions of $\hat{\xi }$ only. The first term represents the effect of the velocity gradient and wall, and causes a migration toward the channel centre. The second term represents the effect of the shear gradient and causes a migration toward the wall. A stale equilibrium occurs at $\hat{\xi }/2 = 0.2$ and 0.8, which are locations between the channel centre $(\hat{\xi } = 0.5)$ and the walls ($\hat{\xi } = 0$ or 1). The present simulation results for deformable capsules at finite inertia as shown in figure 15 agree with this general trend of the equilibrium position. However, the results further show that the equilibrium position is dependent on $Re_{a}$ (and, hence $Re_{C}$) and $Ca$, unlike (3.6). It moves toward the channel centre with increasing $Ca$, but toward the wall with increasing $Re_{C}$. Thus, (3.6) can be modified for a deformable capsule as ${V_m}\sim Re_{a}{(a/\tilde{W})^2}G(\hat{\xi } - {\hat{\xi }_{eq}})$, with ${\hat{\xi }_{eq}}(Re_{C},\; Ca)$ being the equilibrium position; ${V_m} \to 0$ as $\hat{\xi } \to {\hat{\xi }_{eq}}$. Although the numerical data presented in figure 15 are insufficient to conclude the dependence of ${V_m}$ on $Re_{C}$ and $Ca$, the figure shows that ${V_m}$ due to the inertial lift effect is ${\sim} {O}({10^{ - 2}}- {10^{ - 3}})$ over a wide range of parameters, which is of the similar order as the curvature-induced migration discussed above.

Figure 15.

Cross-stream migration of a capsule in a rectilinear channel of square cross-section $({R_C} = 5,W = 4)$ for different (a) $Ca$, and (b) $Re_{C}$.

This scaling allows us to explain the observed direction of capsule migration in curved channels. For weak secondary flows, e.g. ${U_D}\sim {O}({10^{ - 3}})$, the curvature-induced migration velocity ${V_m} > {U_D}$; hence, this results in the migration toward and settling closer to the interior wall. For moderate secondary flows, ${U_D}\sim {O}({10^{ - 2}})$; this yields the net curvature-induced and inertial lift-induced ${V_m}\sim {U_D}$. For this, the equilibrium locations increasingly move toward the exterior wall as ${U_D}$ increases with increasing $Re_{C}$ (figure 12b). For these cases, the equilibrium locations still lie on the $Z = 0$ plane as the deformation/shear gradient and curvature effects can balance against the Dean drag. For stronger secondary flows for which $U_{D}\gtrsim \ {O}({10^{ - 1}})$, the Dean drag dominates, and capsules released on the $Z = 0$ plane quickly move to the exterior wall, then leave the plane to follow the secondary flow, and eventually settle near the vortex centre.

In general, capsules are observed to leave the middle plane for $Re_{C}\gtrsim 5-10$, as noted in figure 13(a). However, a more precise determination of $Re_{C}$ at which this occurs is difficult. In the range $Re_{C} \approx 1- 5$, for which the secondary flow is weaker, a capsule tends to reside in the $Z = 0$ plane for a long time. Any numerical error accrued during this time can cause small variations in Z coordinate of the capsule centre. If the deformation-induced migration in Z which tends to drive the capsule back toward the $Z = 0$ plane is not strong enough, the capsule may eventually leave the middle plane.

Capsule shapes are presented in figure 14(b–e) for different $Re_{C}$. Specific orientations are not considered and views from one side are shown. For $Re_{C}\mathrm{\ \mathbin{\lower.3ex\hbox{$\buildrel< \over {\smash{\scriptstyle\sim}\vphantom{_x}}$}}\ }1$, for which the capsules settle closer to the interior face, the shapes are less asymmetric than that observed in the non-inertial regime. This is because of the increasing distance between the settling location and the interior face, so that the wall effect, streamline curvature effect and shear gradient effect are all reduced. The shape asymmetry is further reduced at even higher $Re_{C}$. Furthermore, for $Re_{C} > 1$, the difference in shapes at various $Re_{C}$ is small, which further supports the notion that capsule migration in this regime is determined primarily by the secondary flow.

Effect of capsule deformability at finite inertia is considered in figure 16, where trajectories are shown for $Re_{C} = 0.8,\; 8$ and 80 but $Ca$ varying from 0.02 to 10. For $Re_{C} = 0.8$ (figure 16a) it is seen that the direction of capsule migration can be either towards the interior or exterior face of the channel. Accordingly, the equilibrium position at this $Re_{C}$ is dependent on $Ca$. This is in stark contrast to what was observed for the non-inertial regime in § 3.1, where the equilibrium position was found to be independent of $Ca$. In general, capsules settle closer to the interior face at larger values of $Ca$, but away from the interior face and closer to the exterior face at smaller values of $Ca$. This result shows the competing effects of capsule deformation, inertial lift and channel curvature. At this small $Re_{C} = 0.8$, the secondary flow, although present, is relatively weak. For large values of $Ca$, the deformation and the curvature effects are dominant over the inertial lift, resulting in the migration toward the interior face. For smaller values of $Ca$, the deformation and the curvature effects are weaker, and the inertial lift dominates, resulting in the migration toward the exterior face. Figure 16(a) further suggests that the equilibrium locations are very sensitive to $Ca$ in the range $Ca < 0.2$ and at the small $Re_{C}$, since all underlying mechanisms of migration, namely, the deformation, inertial lift and curvature effects are comparable. At a higher $Re_{C}\gtrsim 5-10$, as shown in figure 16(b–e) for $Re_{C} = 8$ and 80 for which the secondary flow is stronger, the results become nearly independent of $Ca$. For the large range of $Ca$ considered here, it is observed that the capsules first migrate toward the exterior face while remaining on the midplane. They remain at the fixed radial locations on this plane for some time, and afterward leave the midplane and follow the secondary flow direction to finally settle near the vortex centre. Thus, under the presence of a strong secondary flow, the migration behaviour in general is independent of $Ca$. A few subtle differences based on $Ca$ can be noted for $Re_{C} = 8$ in figure 16(b,c). First, the capsule moves increasingly closer to the exterior face with decreasing $Ca$, before leaving the midplane. This is because of the diminishing effects of the deformation and curvature against the inertial lift and Dean's drag. Second, the equilibrium locations, although close to the vortex centre for all $Ca$, are weakly dependent on $Ca$. Specifically, capsules settle increasingly away from the midplane, and closer to the vortex centre, with decreasing $Ca$. While there is no curvature effect along the Z direction, the deformation effect along this direction causes the capsule to settle closer to the midplane. This difference in the settling locations, as well as how far a capsule moves toward the exterior face before leaving the midplane, for different $Ca$ diminishes with increasing $Re_{C}$, as seen in figure 16(d–e) for $Re_{C} = 80$, due to Dean's drag dominating over the deformation effect.

Figure 16.

Effect of capsule deformability in the inertial regime. Trajectories are shown at different $Ca$ by considering a weaker secondary flow in (a) for $Re_{C} = 0.8$, and stronger secondary flows in (b,c) for $Re_{C} = 8$, and in (d,e) for $Re_{C} = 80$. Line symbols for different $Ca$ are same for all plots and indicated in (a). Radial trajectories versus time are shown in (a,b,d), while cross-sectional trajectories are shown in (c,e).

The focusing of capsules by the Dean vortex as predicted for curved channels is qualitatively similar to that observed earlier for curved tubes (Ebrahimi et al. Reference Ebrahimi, Balogh and Bagchi2021). The capsule equilibrium locations are also quantitatively very close in the two geometries for the range of parameters considered. For example, for $Re_{C} = 40,\; {R_C} = 5,{R_h} = 4$, they differ by less than 3 %. This is because the locations of the vortex centres are also observed to be very close in the two geometries.

Influence of channel radius of curvature ${R_C}$ is also examined. Simulations are performed over a range ${R_C} = 5$ to 15, by keeping $W = H$ and $Re_{C}$ fixed. We observe that the location of the vortex centre with respect to the channel centreline $(R = {R_C},Z = 0)$ does not significantly change in this range; accordingly, the equilibrium locations of the capsule measured relative to the channel centreline are very close even though the channel curvature is altered. For example, for $Re_{C} = 80,\; W = 4$, the equilibrium location relative to the channel centreline varies by less than 2 % of the channel width as ${R_C}$ varies from 5 to 15.

The channel width, however, has a greater influence. To study this, we simulate capsule migration in curved channels of fixed radius of curvature ${R_C} = 7$, and fixed capsule Reynolds number $Re_{a} = 10$, but vary $W = H$ in the range 3 to 6. The corresponding vortex centres and the capsule equilibrium locations are shown in figure 17(a). As seen, with increasing channel width, both of them move further toward the exterior face and further away from the midplane, due to a stronger secondary flow caused by the simultaneous increase in the channel Reynolds number ($Re_{C}$ varies from 60 to 120) as well as the curvature ratio $\kappa $. The equilibrium locations in a curved channel are compared with those in a rectilinear channel of an identical cross-section in figure 17(b). Simulations are performed for capsule migration in rectilinear channels in the inertial regime by releasing capsules from different initial locations. We observe that capsules first migrate toward a specific location with respect to the channel centreline in the cross-sectional plane. The collection of all these locations forms a squircle as shown in figure 17(b). This result is similar to the ring-like focusing position of rigid particles observed by Choi et al. (Reference Choi, Seo and Lee2011) at $Re_{C} = 12$ in their experiment. We further observe that after arriving at the squircle, a capsule can move along the perimeter. We refrain from determining whether there are higher probabilities for capsule focusing at the centres of the four sides of the squircle, as it is the case for rigid particles at higher $Re_{C}$, or along the diagonals as observed by Raffiee et al. (Reference Raffiee, Dabiri and Ardekani2017) and Schaaf & Stark (Reference Schaaf and Stark2017) for capsules. In contrast to such infinite/multiple equilibrium locations in rectilinear channels, two unique focusing locations occur in curved channels as indicated in the figure. Furthermore, for the range of $Re_{C}\ \sim\ 40-320$, we find that the equilibrium locations in the rectilinear channel are not dependent on $Re_{C}$ for a fixed cross-sectional geometry, in agreement with studies using rigid particles (e.g. Schonberg & Hinch Reference Schonberg and Hinch1989; Kim & Yoo Reference Kim and Yoo2008), as well as deformable capsules (Kilimnik et al. Reference Kilimnik, Mao and Alexeev2011; Raffiee et al. Reference Raffiee, Dabiri and Ardekani2017; Schaaf & Stark Reference Schaaf and Stark2017). In contrast, for the curved channel, we predict that the equilibrium locations are dependent on $Re_{C}$ as shown in the figure.

Figure 17.

(a) Influence of channel width (square cross-section) on vortex centres (open symbols) and capsule equilibrium locations (filled symbols). Here ${R_C} = 7$ is fixed. (b) Comparison of equilibrium locations in a curved channel $({R_C} = 5)$ and in a rectilinear channel of the same $W = H = 4$. For the rectilinear channel, the equilibrium locations lie on a squircle as indicated, whereas they are unique locations for the curved channel.

The small difference between the centres of the capsule and the Dean vortex can be explained as follows. Generally, capsules settle further toward the exterior face in the radial direction and closer to the middle plane in the Z direction compared with the vortex centre. The offset in the radial direction is due to the centrifugal force which tends to act radially outward. In contrast, the offset in the Z direction is due to the presence of deformation which tends to act toward the middle plane.

Focusing in curved channels of rectangular cross-sections in the inertial regime is considered next. For this, simulations are performed by varying channel width and height, but keeping the hydraulic diameter, and hence, the channel Reynolds number constant at $Re_{C} = 80$. The radius of curvature is also kept fixed. As shown in figure 18(a) for different cross-sections, capsules follow spiralling trajectories toward the centres of the Dean vortices, as it is the case for the square cross-sections. Figure 18(b) shows the vortex centres and capsule equilibrium locations in the cross-sectional plane for varying width and height. The location of the vortex centre moves further toward the exterior face with increasing width, and further away from the midplane with increasing height. The shift in the radial direction becomes significant at higher $Re_{C}$, for example, $Re_{C} = 320$. As such, the capsule equilibrium location also shifts further toward the exterior face with increasing channel width, and further away from the midplane with increasing height (figure 18b). A key result to note here is that, for channels of smaller widths but larger heights, the radial location of the equilibrium is close that of the vortex centre. But with increasing width and decreasing height, the equilibrium moves further away from the vortex centre and more toward the exterior face. The distance between the vortex centre and the capsule equilibrium increases with increasing $Re_{C}$. This increasing distance is due to the manner in which the streamwise velocity profile in the radial direction changes with increasing channel width. As shown in figure 18(c), the streamwise velocity becomes more skewed toward the exterior face with increasing width and decreasing height, causing the location of zero shear rate/maximum streamwise velocity to increasingly move to that side. The radial location of the streamwise velocity maximum lies further toward the exterior face than the location of the vortex centre. As a result, the capsule deformation effect results in the equilibrium position in the radial direction to lie between the vortex centre and the exterior face. As for the Z location of the capsule equilibrium, it is observed to be slightly closer to the midplane than the vortex centre, because of the deformation effect in the Z direction which tends to drive the capsule toward the midplane.

Figure 18.

Focusing in curved channels of rectangular cross-sections in the inertial regime. Here, channel radius of curvature and hydraulic diameter $({R_C} = 7,\; {D_h} = 8)$ are fixed. (a) Capsule trajectories in the cross-sectional plane are shown for two different heights and widths. The black and red boxes show the channel cross-section $(W,\; H/2)$ for respective cases. (b) Vortex centres (open symbols) and capsule equilibrium locations (filled symbols) are shown for different heights and widths, and for $Re_{C} = 80$ (squares) and $320$ (circles). The inset shows the data plotted using $Z/H$ and $R^{\prime}$. (c) Streamwise fluid velocity, scaled by the corresponding maximum, extracted at the symmetry plane and plotted as a function of the radial distance.