Turbulent drag reduction in plane Couette flow with polymer additives: a direct numerical simulation study

2018 ◽  
Vol 846 ◽  
pp. 482-507 ◽  
Author(s):  
Hao Teng ◽  
Nansheng Liu ◽  
Xiyun Lu ◽  
Bamin Khomami

Drag reduction (DR) in plane Couette flow (PCF) induced by the addition of flexible polymers has been studied via direct numerical simulation (DNS). The similarities and differences in the drag reduction features of PCF and plane Poiseuille flow (PPF) have been examined in detail, particularly in regard to the polymer-induced modification of large-scale structures (LSSs) in the near-wall turbulence. Specifically, it has been demonstrated that in the near-wall region, drag-reduced PCF has features similar to those of drag-reduced PPF; however, in the core region, intriguing differences are found between these two drag-reduced shear flows. Chief among these differences is the significant polymer stretch that arises from the enhanced exchanges between elastic potential energy and turbulent kinetic energy and the commensurate observation of peak values of the conformation tensor components $\unicode[STIX]{x1D60A}_{yy}$ and $\unicode[STIX]{x1D60A}_{zz}$ in this region. This finding is in stark contrast to that of drag-reduced PPF where the polymer stretch and the exchanges between elastic potential energy and turbulent kinetic energy in the core region are insignificant; to this end, in drag-reduced PPF, peak values of the conformation tensor components appear in the near-wall region. Therefore, this study paves the way for understanding the underlying flow physics in drag-reduced PCF, particularly in the context of elastic theory. Moreover, the longitudinal large-scale streaks at the channel centre of drag-reduced PCF are greatly strengthened due to the increased production/dissipation ratio; the LSS imprint effects on the near-wall flow of drag-reduced PCF monotonically increase as the Weissenberg number is enhanced.

2011 ◽  
Vol 687 ◽  
pp. 376-403 ◽  
Author(s):  
Seong Jae Jang ◽  
Hyung Jin Sung ◽  
Per-Åge Krogstad

AbstractThe flow in an axisymmetric contraction fitted to a fully developed pipe flow is experimentally and numerically studied. The reduction in turbulence intensity in the core region of the flow is discussed on the basis of the budgets for the various turbulent stresses as they develop downstream. The contraction generates a corresponding increase in energy in the near-wall region, where the sources for energy production are quite different and of opposite sign compared to the core region, where these effects are caused primarily by vortex stretching. The vortices in the pipe become aligned with the flow as the stretching develops through the contraction. Vortices which originally have a spanwise component in the pipe are stretched into pairs of counter-rotating vortices which become disconnected and aligned with the mean flow. The structures originating in the pipe which are inclined at an angle with respect to the wall are rotated towards the local mean streamlines. In the very near-wall region and the central part of the contraction the flow tends towards two-component turbulence, but these structures are different. The streamwise and azimuthal stresses are dominant in the near-wall region, while the lateral components dominate in the central part of the flow. The two regions are separated by a rather thin region where the flow is almost isotropic.


2014 ◽  
Vol 751 ◽  
Author(s):  
V. Avsarkisov ◽  
S. Hoyas ◽  
M. Oberlack ◽  
J. P. García-Galache

AbstractA new set of numerical simulations of turbulent plane Couette flow in a large box of dimension ($\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}20\pi h,\, 2h,\, 6\pi h$) at Reynolds number $(\mathit{Re}_{\tau }) =125$, 180, 250 and 550 is described and compared with simulations at lower Reynolds numbers, Poiseuille flows and experiments. The simulations present a logarithmic near-wall layer and are used to verify and revise previously known results. It is confirmed that the fluctuation intensities in the streamwise and spanwise directions do not scale well in wall units. The scaling failure occurs both near to and away from the wall. On the contrary, the wall-normal intensity scales in inner units in the near-wall region and in outer units in the core region. The spectral ridge found by Hoyas & Jiménez (Phys. Fluids, vol. 18, 2003, 011702) for the turbulent Poiseuille flow can also be seen in the present flow. Away from the wall, very large-scale motions are found spanning through all the length of the channel. The statistics of these simulations can be downloaded from the webpage of the Chair of Fluid Dynamics.


Author(s):  
Al Hassan Afailal ◽  
Jérémy Galpin ◽  
Anthony Velghe ◽  
Rémi Manceau

CFD simulation tools are increasingly used nowadays to design more fuel-efficient and clean Internal Combustion Engines (ICE). Within this framework, there is a need to benefit from a turbulence model which offers the best compromise between prediction capabilities and computational cost. The Hybrid Temporal LES (HTLES) approach is here retained within the perspective of an application to ICE configurations. HTLES is a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) model based on a solid theoretical framework using temporal filtering. The concept is to model the near-wall region in RANS and to solve the turbulent structures in the core region if the temporal and spatial resolutions are fine enough. In this study, a dedicated sub-model called Elliptic Shielding (ES) is added to HTLES in order to ensure RANS in the near-wall region, regardless of the mesh resolution. A modification of the computation of the total kinetic energy and the dissipation rate was introduced as first adaptions of HTLES towards non-stationary ICE configurations. HTLES is a recent approach, which has not been validated in a wide range of applications. The present study intends to further validate HTLES implemented in CONVERGE code by examining three stationary test cases. The first validation consists of the periodic hill case, which is a standard benchmark case to assess hybrid turbulence models. Then, in order to come closer to real ICE simulations, i.e., with larger Reynolds numbers and coarser near-wall resolutions, the method is validated in the case of a channel flow using wall functions and in the steady flow rig case consisting in an open valve at a fixed lift. HTLES results are compared to RANS k-ω SST and wall-modeled LES σ simulations performed with the same grid and the same temporal resolution. Unlike RANS, satisfactory reproduction of the flow recirculation has been observed with HTLES in the case of periodic hills. The channel flow configuration has underlined the capability of HTLES to predict the wall friction properly. The steady flow rig shows that HTLES combines advantages of RANS and LES in one simulation. On the one hand, HTLES yields mean and rms velocities as accurate as LES since the scale-resolving simulation is triggered in the core region. On the other hand, hybrid RANS/LES at the wall provides accurate pressure drop in contrast with LES performed on the same mesh. Future work will be dedicated to the extension of HTLES to non-stationary flows with moving walls in order to be able to tackle realistic ICE flow configurations.


This study addresses the possible existence of significant counter-gradient diffusion of turbulent energy from the centre towards the walls in plane Couette flow. Data from a direct numerical simulation disprove the counter-gradient diffusion hypothesis and reveal that the turbulent energy production is nearly balanced by the energy dissipation throughout the core region. The results of the simulation are more plausible than the most comprehensive set of measurements available today.


2000 ◽  
Author(s):  
Li Wenzhong ◽  
B. C. Khoo ◽  
Xu Diao

Abstract The present paper is to determine the correction of hot-wire measurements when it is used to measure the shear flows region very close to the non-conducting wall. By numerical simulation of the Navier-Stokes and energy equations using the control volume method, we found that reasonably deployed grid distribution could largely reduce the computational domain size (for a typical Reynolds number for hot-wire near-wall measurements 4.0×10−3∼1.2, the domain boundary placing 650 diameters from the cylinder in front, rear and top is fair enough for accurate simulation, other than the domain boundary which places the 2000 diameters from the cylinder in front and top, and 3000 diameters from the cylinder in rear), and obtain the similar accuracy results for the correction of hot-wire measurements in the near-wall region. Numerical simulation results also show that, only taking the εf,εw (the maximum difference between the respective values of stream function and vorticity on successive iterations) as the criterion for convergence without judge to convergence of the temperature field seems not enough to obtain a convergent simulation result. This may be the possible reason which caused the discrepancy between the simulation results for hot-wire correction when using hot wire to measure the shear flows very close to the non-conducting wall.


Author(s):  
Dongmei Zhou ◽  
Kenneth S. Ball

This paper has two objectives, (1) to examine the effects of spatial resolution, (2) to examine the effects of computational box size, upon turbulence statistics and the amount of drag reduction with and without the control scheme of wall oscillation. Direct numerical simulation (DNS) of the fully developed turbulent channel flow was performed at Reynolds number of 200 based on the wall-shear velocity and the channel half-width by using spectral methods. For the first objective, four different grids were applied to the same computational domain and the biggest impact was observed on the logarithmic law of mean velocity profiles and on the amount of drag reduction with 28.3% for the coarsest mesh and 35.4% for the finest mesh. Other turbulence features such as RMS velocity fluctuations, RMS vorticity fluctuations, and bursting events were either overpredicted or underpredicted through coarse grids. For the second objective, two different minimal channels and one natural full channel were studied and 3% drag reduction difference was observed between the smallest minimal channel of 39.1% and the natural full channel of 36.2%. In the near-wall region, however, the minimal channel flow did not exhibit significant difference in the mean velocity profiles and other lower-order statistics. Finally, from this systematical study, it showed that the accuracy of DNS depends more on the spanwise resolution, and it also confirmed that a minimal channel model is able to catch key structures of turbulence in the near-wall region but is much less expensive.


The time-dependent structure of the wall region of separating, separated, and reattaching flows is considerably different than that of attached turbulent boundary layers. Large-scale structures, whose frequency of passage scales on the freestream velocity and shear layer thickness, produce large Reynolds shearing stresses and most of the turbulence kinetic energy in the outer region of the shear layer and transport it into the low velocity reversed flow next to the wall. This outer flow impresses a near wall streamwise streaky structure of spanwise spacing λ z simultaneously across the wall over a distance of the order of several λ z . The near wall structures produce negligible Reynolds shear stresses and turbulence kinetic energy.


Author(s):  
Zhen Qian ◽  
Minghui Zhang ◽  
Hao Yu ◽  
Fei Wei

Radial profiles of particle velocity in a large scale (418 mm I.D.) downward Circulating Fluidized Bed (CFB downer) were obtained via a Laser Doppler Velocimetry (LDV) system. Results show that particle velocity is gradually increasing along the radial direction and there exists a peak value in the near wall region. Such unique radial profile shape can be explained by the solids accumulating trend in the near wall region of the downer. Experiment results in this large scale downer are also compared with those obtained by other researchers in small scale units so as to investigate the scale-up effect on the radial particle velocity distribution in the downer.


2017 ◽  
Vol 829 ◽  
pp. 751-779 ◽  
Author(s):  
Jinyul Hwang ◽  
Hyung Jin Sung

Direct numerical simulation data of a turbulent boundary layer ($Re_{\unicode[STIX]{x1D70F}}=1000$) were used to investigate the large-scale influences on the vortical structures that contribute to the local skin friction. The amplitudes of the streamwise and wall-normal swirling strengths ($\unicode[STIX]{x1D706}_{x}$and$\unicode[STIX]{x1D706}_{y}$) were conditionally sampled by measuring the large-scale streamwise velocity fluctuations ($u_{l}$). In the near-wall region, the amplitudes of$\unicode[STIX]{x1D706}_{x}$and$\unicode[STIX]{x1D706}_{y}$decreased under negative$u_{l}$rather than under positive$u_{l}$. This behaviour arose from the spanwise motions within the footprints of the large-scale low-speed ($u_{l}<0$) and high-speed structures ($u_{l}>0$). The intense spanwise motions under the footprint of positive$u_{l}$noticeably strengthened the small-scale spanwise velocity fluctuations ($w_{s}$) below the centre of the near-wall vortical structures as compared to$w_{s}$within the footprint of negative$u_{l}$. The streamwise and wall-normal components were attenuated or amplified around the modulated vortical motions, which in turn led to the dependence of the swirling strength on the$u_{l}$event. We quantified the contribution of the modulated vortical motions$\langle -w\unicode[STIX]{x1D714}_{y}\rangle$, which were related to a change-of-scale effect due to the vortex-stretching force, to the local skin friction. In the near-wall region, intense values of$\langle -w\unicode[STIX]{x1D714}_{y}\rangle$were observed for positive$u_{l}$. By contrast, these values were low for negative$u_{l}$, in connection with the amplification of$w_{s}$and$\unicode[STIX]{x1D706}_{y}$by the strong spanwise motions of the positive$u_{l}$. The resultant skin friction induced by the amplified vortical motions within$u_{l}^{+}>2$was responsible for 15 % of the total skin friction generated by the change-of-scale effect. Finally, we applied this analysis to a drag-reduced flow and found that the amplified vortical motions within the footprint of positive$u_{l}$were markedly diminished, which ultimately contributed to the total drag reduction.


Sign in / Sign up

Export Citation Format

Share Document