No Arabic abstract
We discuss an inertial migration of oblate spheroids in a plane channel, where steady laminar flow is generated by a pressure gradient. Our lattice Boltzmann simulations show that spheroids orient in the flow, so that their minor axis coincides with the vorticity direction (a log-rolling motion). Interestingly, for spheroids of moderate aspect ratios, the equilibrium positions relative to the channel walls depend only on their equatorial radius $a$. By analysing the inertial lift force we argue that this force is proportional to $a^3b$, where $b$ is the polar radius, and conclude that the dimensionless lift coefficient of the oblate spheroid does not depend on $b$, and is equal to that of the sphere of radius $a$.
In a shear flow particles migrate to their equilibrium positions in the microchannel. Here we demonstrate theoretically that if particles are inertial, this equilibrium can become unstable due to the Saffman lift force. We derive an expression for the critical Stokes number that determines the onset of instable equilibrium. We also present results of lattice Boltzmann simulations for spherical particles and prolate spheroids to validate the analysis. Our work provides a simple explanation of several unusual phenomena observed in earlier experiments and computer simulations, but never interpreted before in terms of the unstable equilibrium.
The mechanical deformability of single cells is an important indicator for various diseases such as cancer, blood diseases and inflammation. Lab-on-a-chip devices allow to separate such cells from healthy cells using hydrodynamic forces. We perform hydrodynamic simulations based on the lattice-Boltzmann method and study the behavior of an elastic capsule in a microfluidic channel flow in the inertial regime. While inertial lift forces drive the capsule away from the channel center, its deformability favors migration in the opposite direction. Balancing both migration mechanisms, a deformable capsule assembles at a specific equilibrium distance depending on its size and deformability. We find that this equilibrium distance is nearly independent of the channel Reynolds number and falls on a single master curve when plotted versus the Laplace number. We identify a similar master curve for varying particle radius. In contrast, the actual deformation of a capsule strongly depends on the Reynolds number. The lift-force profiles behave in a similar manner as those for rigid particles. Using the Saffman effect, the capsules equilibrium position can be controlled by an external force along the channel axis. While rigid particles move to the center when slowed down, very soft capsules show the opposite behavior. Interestingly, for a specific control force particles are focused on the same equilibrium position independent of their deformability.
We analyze theoretically and experimentally the triadic resonance instability (TRI) of a plane inertial wave in a rotating fluid. Building on the classical triadic interaction equations between helical modes, we show by numerical integration that the maximum growth rate of the TRI is found for secondary waves that do not propagate in the same vertical plane as the primary wave (the rotation axis is parallel to the vertical). In the inviscid limit, we prove this result analytically, in which case the change in the horizontal propagation direction induced by the TRI evolves from $60^circ$ to $90^circ$ depending on the frequency of the primary wave. Thanks to a wave generator with a large spatial extension in the horizontal direction of invariance of the forced wave, we are able to report experimental evidence that the TRI of a plane inertial wave is three-dimensional. The wavevectors of the secondary waves produced by the TRI are shown to match the theoretical predictions based on the maximum growth rate criterion. These results reveal that the triadic resonant interactions between inertial waves are very efficient at redistributing energy in the horizontal plane, normal to the rotation axis.
Previous studies on nonspherical particle-fluid interaction were mostly confined to elongated fiber-like particles, which were observed to induce turbulence drag reduction. However, with the presence of tiny disk-like particles how wall turbulence is modulated and whether drag reduction occurs are still unknown. Motivated by those open questions, we performed two-way coupled direct numerical simulations of inertialess spheroids in turbulent channel flow by an Eulerian-Lagrangian approach. The additional stress accounts for the feedback from inertialess spheroids on the fluid phase. The results demonstrate that both rigid elongated fibers (prolate spheroids) and thin disks (oblate spheroids) can lead to significant turbulence modulations and drag reduction. However, the disk-induced drag reduction is less pronounced than that of rigid fibers with the same volume fraction. Typical features of drag-reduced flows by additives are observed in both flow statistics and turbulence coherent structures. Moreover, in contrast to one-way simulations, the two-way coupled results of spheroidal particles exhibit stronger preferential alignments and lower rotation rates. At the end we propose a drag reduction mechanism by inertialess spheroids and explain the different performance for drag reduction by fibers and disks. We find that the spheroidal particles weaken the quasistreamwise vortices through negative work and, therefore, the Reynolds shear stress is reduced. However, the mean shear stress generated by particles, which is shape-dependent, partly compensates for the reduction of Reynolds shear stress and thus affects the efficiency of drag reduction. The present study implies that tiny disk-like particles can be an alternative drag reduction agent in wall turbulence.
Deformation-induced lateral migration of a bubble slowly rising near a vertical plane wall in a stagnant liquid is numerically and theoretically investigated. In particular, our focus is set on a situation with a short clearance $c$ between the bubble interface and the wall. Motivated by the fact that numerically and experimentally measured migration velocities are considerably higher than the velocity estimated by the available analytical solution using the Fax{e}n mirror image technique for $a/(a+c)ll 1$ (here $a$ is the bubble radius), when the clearance parameter $varepsilon(= c/a)$ is comparable to or smaller than unity, the numerical analysis based on the boundary-fitted finite-difference approach solving the Stokes equation is performed to complement the experiment. The migration velocity is found to be more affected by the high-order deformation modes with decreasing $varepsilon$. The numerical simulations are compared with a theoretical migration velocity obtained from a lubrication study of a nearly spherical drop, which describes the role of the squeezing flow within the bubble-wall gap. The numerical and lubrication analyses consistently demonstrate that when $varepsilonleq 1$, the lubrication effect makes the migration velocity asymptotically $mu V_{B1}^2/(25varepsilon gamma)$ (here, $V_{B1}$, $mu$, and $gamma$ denote the rising velocity, the dynamic viscosity of liquid, and the surface tension, respectively).