A parallel direct-forcing fictitious domain method is employed to perform fully resolved numerical simulations of turbulent channel flow laden with finite-size particles. The effects of the particle-fluid density ratio on the turbulence modulation in the channel flow are investigated at the friction Reynolds number of 180, the particle volume fraction of 0.84%, and the particle-fluid density ratio ranging from 1 to 104.2. The results show that the variation of the flow drag with the particle-fluid density ratio is not monotonic, with a larger flow drag for the density ratio of 10.42, compared to those of unity and 104.2. A significant drag reduction by the particles is observed for large particle-fluid density ratios during the transient stage, but not at the statistically stationary stage. The intensity of particle velocity fluctuations generally decreases with increasing particle inertia, except that the particle streamwise root-mean-square velocity and streamwise-transverse velocity correlation in the near-wall region are largest at the density ratio of the order of 10. The averaged momentum equations are derived with the spatial averaging theorem and are used to analyze the mechanisms for the effects of the particles on the flow drag. The results indicate that the drag-reduction effect due to the decrease in the fluid Reynolds shear stress is counteracted by the drag-enhancement effect due to the increase in the total particle stress or the interphase drag force for the large particle-inertia case. The sum of the total Reynolds stress and particle inner stress contributions to the flow drag is largest at the density ratio of the order of 10, which is the reason for the largest flow drag at this density ratio. The interphase drag force obtained from the averaged momentum equation (the balance theory) is significantly smaller than (but agrees qualitatively with) that from the empirical drag formula based on the phase-averaged slip velocity for large density ratios. For the neutrally buoyant case, the balance theory predicts a positive interphase force on the particles arising from the negative gradient of the particle inner stress, which cannot be predicted by the drag formula based on the phase-averaged slip velocity. In addition, our results show that both particle collision and particle-turbulence interaction play roles in the formation of the inhomogeneous distribution of the particles at the density ratio of the order of 10.

No abstract

Ping Wang (2016) A parallel fictitious domain method for the interface-resolved simulation of particle-laden flows and its application to the turbulent channel flow, Engineering

Interface-resolved direct numerical simulations of downward particle-laden turbulent channel flows are performed by using a direct-forcing fictitious domain method. The effects of the particle settling coefficient, the density ratio (2, 10, and 100), and the particle size on fluid-turbulence interactions are investigated at a bulk Reynolds number of 5746 and a particle volume fraction of 2.36%. Our results indicate that the significant particle-induced reduction in the turbulence intensity does not take place for the downflow at a low density ratio of 2, and the turbulence intensity generally increases with an increasing particle Reynolds number at the same other control parameters, unlike the upflow case. The total turbulent kinetic energy (TKE) in the channel is larger for the downflow than for the upflow at the same particle Reynolds number, whereas the TKE at the channel center is roughly independent of the flow direction when the particle inertia is very large. For a density ratio of 2, the particles aggregate and are preferentially located in the low-speed streaks in the near-wall region, whereas for a density ratio of 10, the particles migrate toward the channel center, similar to the zero-gravity case. The flow friction increases with an increasing settling coefficient for the same density ratio and particle size, and the friction at the density ratio of order (10) is smallest. The pair distribution function shows the transition from the turbulence-dominated feature to the sedimentation-dominated feature, as the settling coefficient increases.

Abstract

The influence of ambient fluid stratification on buoyant miscible jets and plumes is studied theoretically and experimentally. Given a fixed set of jet/plume parameters, and an ambient fluid stratification sandwiched between top and bottom homogeneous densities, a theoretical criterion is identified to show how step-like density profiles constitute the most effective mixers within a broad class of stable density transitions. This is assessed both analytically and experimentally, respectively by establishing rigorous a priori estimates on generalized Morton-Taylor-Turner (MTT) models (Morton et al., Proc. R. Soc. Lond. A, vol. 234, 1956, pp. 1-23; Fischer et al., Mixing in Inland and Coastal Waters. Academic, 1979), and by studying a critical phenomenon determined by the distance between the jet/plume release height with respect to the depth of the ambient density transition. For fluid released sufficiently close to the background density transition, the buoyant jet fluid escapes and rises indefinitely. For fluid released at locations lower than a critical depth, the buoyant fluid stops rising and is trapped indefinitely. A mathematical formulation providing rigorous estimates on MTT models is developed along with nonlinear jump conditions and an exact critical-depth formula that is in good quantitative agreement with the experiments. Our mathematical analysis provides rigorous justification for the critical trapping/escaping criteria, first presented in Caulfield & Woods (J. Fluid Mech., vol. 360, 1998, pp. 229-248), within a class of algebraic density decay rates. Further, the step-like background stratification is shown to be the most efficient mixing profile amongst a broad family of stably stratified profiles sharing the same density transition within a fixed distance. Finally, the analysis uncovers surprising differences between the Gaussian and top-hat profile closures concerning initial mixing of the jet and ambient fluid.

We study the hydrodynamics of a spherical and dumbbell-shaped microswimmer in a tube. Combined with a squirmer model generating tangential surface waves for self-propulsion, a direct-forcing fictitious domain method is employed to simulate the swimming of the microswimmers. We perform the simulations by considering the variations of the swimming Reynolds numbers (Re), the blockage ratios (κ) and the relative distances (ds) between the squirmers of the dumbbell. The results show that the squirmer dumbbell weakens the inertia effects of the fluid more than an individual squirmer. The constrained tube can speed up an inertial pusher (propelled from the rear) and an inertia pusher dumbbell; a greater distance ds results in a slower speed of an inertial pusher dumbbell but a faster speed of an inertial puller (propelled from the front) dumbbell. We also illustrate the swimming stability of a puller (stable) and pusher (unstable) swimming in the tube at Re = 0. At a finite Re, we find that the inertia and the tube constraint competitively affect the swimming stability of the squirmers and squirmer dumbbells. The puller and puller dumbbells swimming in the tube become unstable with increasing Re, whereas an unstable–stable–unstable evolution is found for the pusher and pusher dumbbells. With increasing κ, the puller and puller dumbbells become stable while the pusher and pusher dumbbells become unstable. In addition, we find that a greater ds yields a higher hydrodynamic efficiency η of the inertial squirmer dumbbell.

No abstract