No Arabic abstract
Numerical simulations of turbulent channel flows, with or without additives, are limited in the extent of the Reynolds number Re and Deborah number De. The comparison of such simulations to theories of drag reduction, which are usually derived for asymptotically high Re and De, calls for some care. In this paper we present a study of drag reduction by rodlike polymers in a turbulent channel flow using direct numerical simulation and illustrate how these numerical results should be related to the recently developed theory.
We address the Additive Equivalence discovered by Virk and coworkers: drag reduction affected by flexible and rigid rodlike polymers added to turbulent wall-bounded flows is limited from above by a very similar Maximum Drag Reduction (MDR) asymptote. Considering the equations of motion of rodlike polymers in wall-bounded turbulent ensembles, we show that although the microscopic mechanism of attaining the MDR is very different, the macroscopic theory is isomorphic, rationalizing the interesting experimental observations.
We address the phenomenon of drag reduction by dilute polymeric additive to turbulent flows, using Direct Numerical Simulations (DNS) of the FENE-P model of viscoelastic flows. It had been amply demonstrated that these model equations reproduce the phenomenon, but the results of DNS were not analyzed so far with the goal of interpreting the phenomenon. In order to construct a useful framework for the understanding of drag reduction we initiate in this paper an investigation of the most important modes that are sustained in the viscoelastic and Newtonian turbulent flows respectively. The modes are obtained empirically using the Karhunen-Loeve decomposition, allowing us to compare the most energetic modes in the viscoelastic and Newtonian flows. The main finding of the present study is that the spatial profile of the most energetic modes is hardly changed between the two flows. What changes is the energy associated with these modes, and their relative ordering in the decreasing order from the most energetic to the least. Modes that are highly excited in one flow can be strongly suppressed in the other, and vice versa. This dramatic energy redistribution is an important clue to the mechanism of drag reduction as is proposed in this paper. In particular there is an enhancement of the energy containing modes in the viscoelastic flow compared to the Newtonian one; drag reduction is seen in the energy containing modes rather than the dissipative modes as proposed in some previous theories.
The flow of fluids in channels, pipes or ducts, as in any other wall-bounded flow (like water along the hulls of ships or air on airplanes) is hindered by a drag, which increases many-folds when the fluid flow turns from laminar to turbulent. A major technological problem is how to reduce this drag in order to minimize the expense of transporting fluids like oil in pipelines, or to move ships in the ocean. It was discovered in the mid-twentieth century that minute concentrations of polymers can reduce the drag in turbulent flows by up to 80%. While experimental knowledge had accumulated over the years, the fundamental theory of drag reduction by polymers remained elusive for a long time, with arguments raging whether this is a skin or a bulk effect. In this colloquium review we first summarize the phenomenology of drag reduction by polymers, stressing both its universal and non-universal aspects, and then proceed to review a recent theory that provides a quantitative explanation of all the known phenomenology. We treat both flexible and rod-like polymers, explaining the existence of universal properties like the Maximum Drag Reduction (MDR) asymptote, as well as non-universal cross-over phenomena that depend on the Reynolds number, on the nature of the polymer and on its concentration. Finally we also discuss other agents for drag reduction with a stress on the important example of bubbles.
Turbulent flows under transcritical conditions are present in regenerative cooling systems of rocker engines and extraction processes in chemical engineering. The turbulent flows and the corresponding heat transfer phenomena in these complex processes are still not well understood experimentally and numerically. The objective of this work is to investigate the turbulent flows under transcritical conditions using DNS of turbulent channel flows. A fully compressible solver is used in conjunction with a Peng-Robinson real-fluid equation of state to describe the transcritical flows. A channel flow with two isothermal walls is simulated with one heated and one cooled boundary layers. The grid resolution adopted in this study is slightly finer than that required for DNS of incompressible channel flows. The simulations are conducted using both fully (FC) and quasi-conservative (QC) schemes to assess their performance for transcritical wall-bounded flows. The instantaneous flows and the statistics are analyzed and compared with the canonical theories. It is found that results from both FC and QC schemes qualitatively agree well with noticeable difference near the top heated wall, where spurious oscillations in velocity can be observed. Using the DNS data, we then examine the usefulness of Townsend attached eddy hypothesis in the context of flows at transcritical conditions. It is shown that the streamwise energy spectrum exhibits the inverse wavenumber scaling and that the streamwise velocity structure function follows a logarithmic scaling, thus providing support to the attached eddy model at transcritical conditions.
We compare experimental data and numerical simulations for the dynamics of inertial particles with finite density in turbulence. In the experiment, bubbles and solid particles are optically tracked in a turbulent flow of water using an Extended Laser Doppler Velocimetry technique. The probability density functions (PDF) of particle accelerations and their auto-correlation in time are computed. Numerical results are obtained from a direct numerical simulation in which a suspension of passive pointwise particles is tracked, with the same finite density and the same response time as in the experiment. We observe a good agreement for both the variance of acceleration and the autocorrelation timescale of the dynamics; small discrepancies on the shape of the acceleration PDF are observed. We discuss the effects induced by the finite size of the particles, not taken into account in the present numerical simulations.