Prediction of polymer extension, drag reduction, and vortex interaction in direct numerical simulation of turbulent channel flows

Abstract

Hydrodynamic and viscoelastic interactions between the turbulent fluid within a channel at [Formula: see text] and a polymeric phase are investigated numerically using a multiscale hybrid approach. Direct numerical simulations are performed to predict the continuous phase and Brownian dynamics simulations using the finitely extensible nonlinear elastic dumbbell approach are carried out to model the trajectories of polymer extension vectors within the flow, using parallel computations to achieve reasonable computation timeframes on large-scale flows. Upon validating the polymeric configuration solver against theoretical predictions in equilibrium conditions, with excellent agreement observed, the distributions of velocity gradient tensor components are analyzed throughout the channel flow wall-normal regions. Impact on polymer stretching is discussed, with streamwise extension dominant close to the wall, and wall-normal extension driven by high streamwise gradients of wall-normal velocity. In this case, it is shown that chains already possessing high wall-normal extensions may attempt to orientate more in the streamwise direction, causing a curling effect. These effects are observed in instantaneous snapshots of polymer extension, and the effects of the bulk Weissenberg number show that increased [Formula: see text] leads to more stretched configurations and more streamwise orientated conformities close to the wall, whereas, in the bulk flow and log-law regions, the polymers tend to trace fluid turbulence structures. Chain orientation angles are also considered, with [Formula: see text] demonstrating little influence on the isotropic distributions in the log-law and bulk flow regions. Polymer–fluid coupling is implemented through a polymer contribution to the viscoelastic stress tensor. The effect of the polymer relaxation time on the turbulent drag reduction is discussed, with greater Weissenberg numbers leading to more impactful reduction. Finally, the velocity gradient tensor invariants are calculated for the drag-reduced flows, with polymers having a significant impact on the Q–R phase diagrams, with the presence of polymers narrowing the range of [Formula: see text] values in the wall regions and causing flow structures to become more two-dimensional.