Motion of an elastic capsule in a square microfluidic channel.

In the present study we investigate computationally the steady-state motion of an elastic capsule along the centerline of a square microfluidic channel and compare it with that in a cylindrical tube. In particular, we consider a slightly over-inflated elastic capsule made of a strain-hardening membrane with comparable shearing and area-dilatation resistance. Under the conditions studied in this paper (i.e., small, moderate, and large capsules at low and moderate flow rates), the capsule motion in a square channel is similar to and thus governed by the same scaling laws with the capsule motion in a cylindrical tube, even though in the channel the cross section in the upstream portion of large capsules is nonaxisymmetric (i.e., square-like with rounded corners). When the hydrodynamic forces on the membrane increase, the capsule develops a pointed downstream edge and a flattened rear (possibly with a negative curvature) so that the restoring tension forces are increased as also happens with droplets. Membrane tensions increase significantly with the capsule size while the area near the downstream tip is the most probable to rupture when a capsule flows in a microchannel. Because the membrane tensions increase with the interfacial deformation, a suitable Landau-Levich-Derjaguin-Bretherton analysis reveals that the lubrication film thickness h for large capsules depends on both the capillary number Ca and the capsule size a; our computations determine the latter dependence to be (in dimensionless form) h ~ a(-2) for the large capsules studied in this work. For small and moderate capsule sizes a, the capsule velocity Ux and additional pressure drop ΔP+ are governed by the same scaling laws as for high-viscosity droplets. The velocity and additional pressure drop of large thick capsules also follow the dynamics of high-viscosity droplets, and are affected by the lubrication film thickness. The motion of our large thick capsules is characterized by a Ux-U ~ h ~ a(-2) approach to the undisturbed average duct velocity and an additional pressure drop ΔP+ ~a(3)/h ~ a(5). By combining basic physical principles and geometric properties, we develop a theoretical analysis that explains the power laws we found for large capsules.