The motion of a deforming capsule through a corner

A three-dimensional deformable capsule convected through a square duct with a corner is studied via numerical simulations. We develop an accelerated boundary integral implementation adapted to general geometries and boundary conditions. A global spectral method is adopted to resolve the dynamics of the capsule membrane developing elastic tension according to the neo-Hookean constitutive law and bending moments in an inertialess flow. The simulations show that the trajectory of the capsule closely follows the underlying streamlines independently of the capillary number. The membrane deformability, on the other hand, significantly influences the relative area variations, the advection velocity and the principal tensions observed during the capsule motion. The evolution of the capsule velocity displays a loss of the time-reversal symmetry of Stokes flow due to the elasticity of the membrane. The velocity decreases while the capsule is approaching the corner as the background flow does, reaches a minimum at the corner and displays an overshoot past the corner due to the streamwise elongation induced by the flow acceleration in the downstream branch. This velocity overshoot increases with confinement while the maxima of the major principal tension increase linearly with the inverse of the duct width. Finally, the deformation and tension of the capsule are shown to decrease in a curved corner

The effect of polydispersity in a turbulent channel flow laden with finite-size particles

We study turbulent channel flows of monodisperse and polydisperse suspensions of finite-size spheres by means of Direct Numerical Simulations using an immersed boundary method to account for the dispersed phase. Suspensions with 3 different Gaussian distributions of particle radii are considered (i.e. 3 different standard deviations). The distributions are centered on the reference particle radius of the monodisperse suspension. In the most extreme case, the radius of the largest particles is 4 times that of the smaller particles. We consider two different solid volume fractions, 2% and 10%. We find that for all polydisperse cases, both fluid and particles statistics are not substantially altered with respect to those of the monodisperse case. Mean streamwise fluid and particle velocity profiles are almost perfectly overlapping. Slightly larger differences are found for particle velocity fluctuations. These increase close to the wall and decrease towards the centerline as the standard deviation of the distribution is increased. Hence, the behavior of the suspension is mostly governed by excluded volume effects regardless of particle size distribution (at least for the radii here studied). Due to turbulent mixing, particles are uniformly distributed across the channel. However, smaller particles can penetrate more into the viscous and buffer layer and velocity fluctuations are therein altered. Non trivial results are presented for particle-pair statistics.Comment: Under review in the European Journal of Mechanics/B - Fluid

Sedimentation of finite-size spheres in quiescent and turbulent environments

Sedimentation of a dispersed solid phase is widely encountered in applications and environmental flows, yet little is known about the behavior of finite-size particles in homogeneous isotropic turbulence. To fill this gap, we perform Direct Numerical Simulations of sedimentation in quiescent and turbulent environments using an Immersed Boundary Method to account for the dispersed rigid spherical particles. The solid volume fractions considered are 0.5-1%, while the solid to fluid density ratio 1.02. The particle radius is chosen to be approximately 6 Komlogorov lengthscales. The results show that the mean settling velocity is lower in an already turbulent flow than in a quiescent fluid. The reduction with respect to a single particle in quiescent fluid is about 12\% and 14\% for the two volume fractions investigated. The probability density function of the particle velocity is almost Gaussian in a turbulent flow, whereas it displays large positive tails in quiescent fluid. These tails are associated to the intermittent fast sedimentation of particle pairs in drafting-kissing-tumbling motions. The particle lateral dispersion is higher in a turbulent flow, whereas the vertical one is, surprisingly, of comparable magnitude as a consequence of the highly intermittent behavior observed in the quiescent fluid. Using the concept of mean relative velocity we estimate the mean drag coefficient from empirical formulas and show that non stationary effects, related to vortex shedding, explain the increased reduction in mean settling velocity in a turbulent environment.Comment: In press on Journal of Fluid Mechanic

Transition and self-sustained turbulence in dilute suspensions of finite-size particles

We study the transition to turbulence of channel flow of finite-size particle suspensions at low volume fraction, i.e. $\Phi \approx 0.001$. The critical Reynolds number above which turbulence is sustained reduces to $Re \approx 1675$, in the presence of few particles, independently of the initial condition, a value lower than that of the corresponding single-phase flow, i.e. $Re\approx1775$. In the dilute suspension, the initial arrangement of the particles is important to trigger the transition at a fixed Reynolds number and particle volume fraction. As in single phase flows, streamwise elongated disturbances are initially induced in the flow. If particles can induce oblique disturbances with high enough energy within a certain time, the streaks breakdown, flow experiences the transition to turbulence and the particle trajectories become chaotic. Otherwise, the streaks decay in time and the particles immigrate towards the channel core in a laminar flow.Comment: 12 pages, 3 figure

Clustering and increased settling speed of oblate particles at finite Reynolds number

We study the settling of rigid oblates in quiescent fluid using interface-resolved Direct Numerical Simulations. In particular, an immersed boundary method is used to account for the dispersed solid phase together with lubrication correction and collision models to account for short-range particle-particle interactions. We consider semi-dilute suspensions of oblate particles with aspect ratio AR=1/3 and solid volume fractions $\phi=0.5\%-10\%$. The solid-to-fluid density ratio $R=1.5$ and the Galileo number (i.e. the ratio between buoyancy and viscous forces) based on the diameter of a sphere with equivalent volume $Ga=60$. With this choice of parameters, an isolated oblate falls vertically with a steady wake with its broad side perpendicular to the gravity direction. At this $Ga$, the mean settling speed of spheres is a decreasing function of the volume $\phi$ and is always smaller than the terminal velocity of the isolated particle, $V_t$. On the contrary, we show here that the mean settling speed of oblate particles increases with $\phi$ in dilute conditions and is $33\%$ larger than $V_t$. At higher concentrations, the mean settling speed decreases becoming smaller than the terminal velocity $V_t$ between $\phi=5\%$ and $10\%$. The increase of the mean settling speed is due to the formation of particle clusters that for $\phi=0.5\%-1\%$ appear as columnar-like structures. From the pair-distribution function we observe that it is most probable to find particle-pairs almost vertically aligned. However, the pair-distribution function is non-negligible all around the reference particle indicating that there is a substantial amount of clustering at radial distances between 2 and $6c$ (with $c$ the polar radius of the oblate).Comment: Submitted to Journal of Fluid Mechanic