Physical mechanisms governing drag reduction in turbulent Taylor-Couette flow with finite-size deformable bubbles

The phenomenon of drag reduction induced by injection of bubbles into a turbulent carrier fluid has been known for a long time; the governing control parameters and underlying physics is however not well understood. In this paper, we use three dimensional numerical simulations to uncover the effect of deformability of bubbles injected in a turbulent Taylor-Couette flow on the overall drag experienced by the system. We consider two different Reynolds numbers for the carrier flow, i.e. $Re_i=5\times 10^3$ and $Re_i=2\times 10^4$; the deformability of the bubbles is controlled through the Weber number which is varied in the range $We=0.01 - 2.0$. Our numerical simulations show that increasing the deformability of bubbles i.e., $We$ leads to an increase in drag reduction. We look at the different physical effects contributing to drag reduction and analyse their individual contributions with increasing bubble deformability. Profiles of local angular velocity flux show that in the presence of bubbles, turbulence is enhanced near the inner cylinder while attenuated in the bulk and near the outer cylinder. We connect the increase in drag reduction to the decrease in dissipation in the wake of highly deformed bubbles near the inner cylinder

Deformable ellipsoidal bubbles in Taylor-Couette flow with enhanced Euler-Lagrange tracking

In this work we present numerical simulations of $10^5$ sub-Kolmogorov deformable bubbles dispersed in Taylor-Couette flow (a wall-bounded shear system) with rotating inner cylinder and outer cylinder at rest. We study the effect of deformability of the bubbles on the overall drag induced by the carrier fluid in the two-phase system. We find that an increase in deformability of the bubbles results in enhanced drag reduction due to a more pronounced accumulation of the deformed bubbles near the driving inner wall. This preferential accumulation is induced by an increase in the resistance on the motion of the bubbles in the wall-normal direction. The increased resistance is linked to the strong deformation of the bubbles near the wall which makes them prolate (stretched along one axes) and orient along the stream-wise direction. A larger concentration of the bubbles near the driving wall implies that they are more effective in weakening the plume ejections which results in stronger drag reduction effects. These simulations which are practically impossible with fully resolved techniques are made possible by coupling a sub-grid deformation model with two-way coupled Euler-Lagrangian tracking of sub-Kolmogorov bubbles dispersed in a turbulent flow field which is solved through direct numerical simulations. The bubbles are considered to be ellipsoidal in shape and their deformation is governed by an evolution equation which depends on the local flow conditions and their surface tension

Flow induced dissolution of femtoliter surface droplet arrays

The dissolution of liquid nanodroplets is a crucial step in many applied processes, such as separation and dispersion in food industry, crystal formation of pharmaceutical products, concentrating and analysis in medical diagnosis, and drug delivery in aerosols. In this work, using both experiments and numerical simulations, we \textit{quantitatively} study the dissolution dynamics of femtoliter surface droplets in a highly ordered array under a uniform flow. Our results show that the dissolution of femoliter droplets strongly depends on their spatial positions relative to the flow direction, drop-to-drop spacing in the array, and the imposed flow rate. In some particular case, the droplet at the edge of the array can dissolve about 30% faster than the ones located near the centre. The dissolution rate of the droplet increases by 60% as the inter-droplet spacing is increased from 2.5 $\mu$m to 20 $\mu$m. Moreover, the droplets close to the front of flow commence to shrink earlier than those droplets in the center of the array. The average dissolution rate is faster for faster flow. As a result, the dissolution time $T_{i}$ decreases with the Reynolds number Re of the flow as $T_{i}\propto Re^{-3/4}$. The experimental results are in good agreement with numerical simulations where the advection-diffusion equation for the concentration field is solved and the concentration gradient on the surface of the drop is computed. The findings suggest potential approaches to manipulate nanodroplet sizes in droplet arrays simply by dissolution controlled by an external flow. The obtained droplets with varying curvatures may serve as templates for generating multifocal microlens in one array

Fluctuation-induced force in homogeneous isotropic turbulence.

Understanding force generation in nonequilibrium systems is a notable challenge in statistical physics. We uncover a fluctuation-induced force between two plates immersed in homogeneous isotropic turbulence using direct numerical simulations. The force is a nonmonotonic function of plate separation. The mechanism of force generation reveals an intriguing analogy with fluctuation-induced forces: In a fluid, energy and vorticity are localized in regions of defined length scales. When varying the distance between the plates, we exclude energy structures modifying the overall pressure on the plates. At intermediate plate distances, the intense vorticity structures (worms) are forced to interact in close vicinity between the plates. This interaction affects the pressure in the slit and the force between the plates. The combination of these two effects causes a nonmonotonic attractive force with a complex Reynolds number dependence. Our study sheds light on how length scale-dependent distributions of energy and high-intensity vortex structures determine Casimir forces

A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes

In this paper we show and discuss the use of a versatile interaction potential approach coupled with an immersed boundary method to simulate a variety of flows involving deformable bodies. In particular, we focus on two kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii) flow in the left ventricle of the heart with either a mechanical or a natural valve. Both examples have in common the two-way interaction of the flow with a deformable interface or a membrane. The interaction potential approach (de Tullio & Pascazio, Jou. Comp. Phys., 2016; Tanaka, Wada and Nakamura, Computational Biomechanics, 2016) with minor modifications can be used to capture the deformation dynamics in both classes of problems. We show that the approach can be used to replicate the deformation dynamics of liquid-liquid interfaces through the use of ad-hoc elastic constants. The results from our simulations agree very well with previous studies on the deformation of drops in standard flow configurations such as deforming drop in a shear flow or a cross flow. We show that the same potential approach can also be used to study the flow in the left ventricle of the heart. The flow imposed into the ventricle interacts dynamically with the mitral valve (mechanical or natural) and the ventricle which are simulated using the same model. Results from these simulations are compared with ad- hoc in-house experimental measurements. Finally, a parallelisation scheme is presented, as parallelisation is unavoidable when studying large scale problems involving several thousands of simultaneously deforming bodies on hundreds of distributed memory computing processors

The Queen's awards for export, technological and environmental achievement General background note

Available from British Library Document Supply Centre-DSC:99/36493 / BLDSC - British Library Document Supply CentreSIGLEGBUnited Kingdo

Drag reduction in numerical two-phase Taylor–Couette turbulence using an Euler–Lagrange approach

Two-phase turbulent Taylor–Couette (TC) flow is simulated using an Euler–Lagrange approach to study the effects of a secondary phase dispersed into a turbulent carrier phase (here bubbles dispersed into water). The dynamics of the carrier phase is computed using direct numerical simulations (DNS) in an Eulerian framework, while the bubbles are tracked in a Lagrangian manner by modelling the effective drag, lift, added mass and buoyancy force acting on them. Two-way coupling is implemented between the dispersed phase and the carrier phase which allows for momentum exchange among both phases and to study the effect of the dispersed phase on the carrier phase dynamics. The radius ratio of the TC setup is fixed to η=0.833, and a maximum inner cylinder Reynolds number of Rei=8000 is reached. We vary the Froude number (Fr), which is the ratio of the centripetal to the gravitational acceleration of the dispersed phase and study its effect on the net torque required to drive the TC system. For the two-phase TC system, we observe drag reduction, i.e. the torque required to drive the inner cylinder is lower compared with that of the single-phase system. The net drag reduction decreases with increasing Reynolds number Rei, which is consistent with previous experimental findings (Murai et al., J. Phys.: Conf. Ser., vol. 14, 2005, pp. 143–156; Phys. Fluids, vol. 20(3), 2008, 034101). The drag reduction is strongly related to the Froude number: for fixed Reynolds number we observe higher drag reduction when Fr1. This buoyancy effect is more prominent in low Rei systems and decreases with increasing Reynolds number Rei. We trace the drag reduction back to the weakening of the angular momentum carrying Taylor rolls by the rising bubbles. We also investigate how the motion of the dispersed phase depends on Rei and Fr, by studying the individual trajectories and mean dispersion of bubbles in the radial and axial directions. Indeed, the less buoyant bubbles (large Fr) tend to get trapped by the Taylor rolls, while the more buoyant bubbles (small Fr) rise through and weaken them