Particulate flows have mainly been studied under the simplifying assumption of a one-way coupling regime where the disperse phase does not modify the carrier fluid. A more complete view of multiphase ...flows can be gained calling into play two-way coupling effects, i.e. by accounting for the inter-phase momentum exchange, which is certainly relevant at increasing mass loading. In this paper we present a new methodology rigorously designed to capture the inter-phase momentum exchange for particles smaller than the smallest hydrodynamical scale, e.g. the Kolmogorov scale in a turbulent flow. The momentum coupling mechanism exploits the unsteady Stokes flow around a small rigid sphere, where the transient disturbance produced by each particle is evaluated in a closed form. The particles are described as lumped point masses, which would lead to the appearance of singularities. A rigorous regularization procedure is conceived to extract the physically relevant interactions between the particles and the fluid which avoids any ‘ad hoc’ assumption. The approach is suited for high-efficiency implementation on massively parallel machines since the transient disturbance produced by the particles is strongly localized in space. We will show that hundreds of thousands of particles can be handled at an affordable computational cost, as demonstrated by a preliminary application to a particle-laden turbulent shear flow.
The hydrodynamics of a flagellated micro-organism is investigated when swimming close to a planar free-slip surface by means of numerical solutions of the Stokes equations obtained via a boundary ...element method. Depending on the initial conditions, the swimmer can either escape from the free-slip surface or collide with the boundary. Interestingly, the micro-organism does not exhibit a stable orbit. Independently of escape or attraction to the interface, close to a free-slip surface, the swimmer follows a counter-clockwise trajectory, in agreement with experimental findings (Di Leonardo et al., Phys. Rev. Lett., vol. 106 (3), 2011, 038101). The hydrodynamics is indeed modified by the free surface. In fact, when the same swimmer moves close to a no-slip wall, a set of initial conditions exists which result in stable orbits. Moreover, when moving close to a free-slip or a no-slip boundary, the swimmer assumes a different orientation with respect to its trajectory. Taken together, these results contribute to shed light on the hydrodynamical behaviour of micro-organisms close to liquid–air interfaces which are relevant for the formation of interfacial biofilms of aerobic bacteria.
Turbulent flow separation induced by a protuberance on one of the walls of an otherwise planar channel is investigated using direct numerical simulations. Different bulge geometries and Reynolds ...numbers – with the highest friction Reynolds number simulation reaching a peak of
$Re_{\unicodeSTIX{x1D70F}}=900$
– are addressed to understand the effect of the wall curvature and of the Reynolds number on the dynamics of the recirculating bubble behind the bump. Global quantities reveal that most of the drag is due to the form contribution, whilst the friction contribution does not change appreciably with respect to an equivalent planar channel flow. The size and position of the separation bubble strongly depends on the bump shape and the Reynolds number. The most bluff geometry has a larger recirculation region, whilst the Reynolds number increase results in a smaller recirculation bubble and a shear layer more attached to the bump. The position of the reattachment point only depends on the Reynolds number, in agreement with experimental data available in the literature. Both the mean and the turbulent kinetic energy equations are addressed in such non-homogeneous conditions revealing a non-trivial behaviour of the energy fluxes. The energy introduced by the pressure drop follows two routes: part of it is transferred towards the walls to be dissipated and part feeds the turbulent production hence the velocity fluctuations in the separating shear layer. Spatial energy fluxes transfer the kinetic energy into the recirculation bubble and downstream near the wall where it is ultimately dissipated. Consistently, anisotropy concentrates at small scales near the walls irrespective of the value of the Reynolds number. In the bulk flow and in the recirculation bubble, isotropy is restored at small scales and the isotropy recovery rate is controlled by the Reynolds number. Anisotropy invariant maps are presented, showing the difficulty in developing suitable turbulence models to predict separated turbulent flow dynamics. Results shed light on the processes of production, transfer and dissipation of energy in this relatively complex turbulent flow where non-homogeneous effects overwhelm the classical picture of wall-bounded turbulent flows which typically exploits streamwise homogeneity.
Polymer chains in turbulent flows are generally modelled as dumbbells, i.e. two beads joined by a nonlinear spring. The dumbbell only maps a single spatial configuration, described by the polymer ...end-to-end vector, thus a multi-bead FENE (finitely extensible nonlinear elastic) chain seems a natural improvement for a more accurate characterisation of the polymer spatial conformation. At a large Weissenberg number, a comparison with the more accurate Kuhn chain reveals that the multi-bead FENE chain drastically overestimates the probability of folded configurations. Surprisingly, the dumbbell turns out to be the only meaningful bead-spring model to coarse-grain a polymer macromolecule in turbulent pipe flows.
Two-way momentum-coupled direct numerical simulations of a particle-laden turbulent channel flow are addressed to investigate the effect of the particle Stokes number and of the particle-to-fluid ...density ratio on the turbulence modification. The exact regularised point-particle method is used to model the interphase momentum exchange in presence of solid boundaries, allowing the exploration of an extensive region of the parameter space. Results show that the particles increase the friction drag in the parameter space region considered, namely the Stokes number $St_+ \in 2,80$, and the particle-to-fluid density ratio $\rho _p/\rho _f \in 90,5760$ at a fixed mass loading $\phi =0.4$. It is noteworthy that the highest drag occurs for small Stokes number particles. A measurable drag increase occurs for all particle-to-fluid density ratios, the effect being reduced significantly only at the highest value of $\rho _p/\rho _f$. The modified stress budget and turbulent kinetic energy equation provide the rationale behind the observed behaviour. The particles’ extra stress causes an additional momentum flux towards the wall that modifies the structure of the buffer and of the viscous sublayer where the streamwise and wall-normal velocity fluctuations are increased. Indeed, in the viscous sublayer, additional turbulent kinetic energy is produced by the particles’ back-reaction, resulting in a strong augmentation of the spatial energy flux towards the wall where the energy is ultimately dissipated. This behaviour explains the increase of friction drag in particle-laden wall-bounded flows.
The generalised Kolmogorov equation is used to describe the scale-by-scale turbulence dynamics in the shear layer and in the separation bubble generated by a bulge at one of the walls in a turbulent ...channel flow. The second-order structure function, which is the basis of such an equation, is used as a proxy to define a scale-energy content, that is an interpretation of the energy associated with a given scale. Production and dissipation regions and the flux interchange between them, in both physical and separation space, are identified. Results show how the generalised Kolmogorov equation, a five-dimensional equation in our anisotropic and strongly inhomogeneous flow, can describe the turbulent flow behaviour and related energy mechanisms. Such complex statistical observables are linked to a visual inspection of instantaneous turbulent structures detected by means of the Q-criterion. Part of these turbulent structures are trapped in the recirculation where they undergo a pseudo-cyclic process of disruption and reformation. The rest are convected downstream, grow and tend to larger streamwise scales in an inverse cascade. The classical picture of homogeneous isotropic turbulence in which energy is fed at large scales and transferred to dissipate at small scales does not simply apply to this flow where the energy dynamics strongly depends on position, orientation and length scale.
The Exact Regularized Point Particle (ERPP) method is extended to treat the interphase momentum coupling between particles and fluid in the presence of walls by accounting for the vorticity ...generation due to the particles close to solid boundaries. The ERPP method overcomes the limitations of other methods by allowing the simulation of an extensive parameter space (Stokes number, mass loading, particle-to-fluid density ratio and Reynolds number) and of particle spatial distributions that are uneven (few particles per computational cell). The enhanced ERPP method is explained in detail and validated by considering the global impulse balance. In conditions when particles are located close to the wall, a common scenario in wall-bounded turbulent flows, the main contribution to the total impulse arises from the particle-induced vorticity at the solid boundary. The method is applied to direct numerical simulations of particle-laden turbulent pipe flow in the two-way coupling regime to address the turbulence modulation. The effects of the mass loading, the Stokes number and the particle-to-fluid density ratio are investigated. The drag is either unaltered or increased by the particles with respect to the uncoupled case. No drag reduction is found in the parameter space considered. The momentum stress budget, which includes an extra stress contribution by the particles, provides the rationale behind the drag behaviour. The extra stress produces a momentum flux towards the wall that strongly modifies the viscous stress, the culprit of drag at solid boundaries.
The drag force acting on a body moving in a fluid has two components, friction drag due to fluid viscosity and form drag due to flow separation behind the body. When present, form drag is usually the ...most significant between the two, and in many applications, streamlining efficiently reduces or prevents flow separation. As studied here, when the operating fluid is water, a promising technique for form drag reduction is to modify the walls of the body with superhydrophobic surfaces. These surfaces entrap gas bubbles in their asperities, avoiding the direct contact of the liquid with the wall. Superhydrophobic surfaces have been vastly studied for reducing friction drag. We show they are also effective in reducing flow separation in turbulent flow and therefore in reducing the form drag. Their conceptual effectiveness is demonstrated by performing direct numerical simulations of turbulent flow over a bluff body, represented by a bump inside a channel, which is modified with different superhydrophobic surfaces. The approach shown here contributes to new and powerful techniques for drag reduction on bluff bodies.