Higher-order and length-scale statistics from DNS of a decelerated planar wall-bounded turbulent flow

A DNS database is presented to document third- and fourth-order moments and their budgets for fully developed plane-channel flow and for strained plane-channel flow. The effect of straining has a similar effect on statistics, such as the skewness and flatness of velocities, as in an adverse-pressure-gradient (APG) boundary layer. In addition to higher-order statistics, some modeling implications are also described, including issues related to decomposition of the velocity–pressure gradient correlations, the assumption of dissipation isotropy and a fourth-order turbulence time scale. An analysis of two-point correlations along the inhomogeneous direction is made to include the effects of APG straining. This reveals the advantage of a one-point length-scale based on wall-normal velocity fluctuations as an alternative to the traditionally used length-scale based on turbulence kinetic energy. The present study should prove to be useful for turbulence modelers in need of data to develop and assess higher-than-second-order Reynolds-averaged Navier–Stokes closure models.     

DNS and LES of turbulent flow in a closed channel featuring a pattern of hemispherical roughness elements

Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) were performed for fully-developed turbulent flow in channels with smooth walls and walls featuring hemispherical roughness elements at shear Reynolds numbers Re_τ = 180 and 400, with the goal of studying the effect of these roughness elements on the wall-layer structure and on the friction factor. The LES and DNS approaches were verified first by comparison with existing DNS databases for smooth walls. Then, a parametric study for the hemispherical roughness elements was conducted, including the effects of shear Reynolds number, normalized roughness height (k⁺ = 10–20) and relative roughness spacing (s⁺/k⁺ = 2–6). The sensitivity study also included the effect of distribution pattern (regular square lattice vs. random pattern) of the roughness elements on the walls. The hemispherical roughness elements generate turbulence, thus increasing the friction factor with respect to the smooth-wall case, and causing a downward shift in the mean velocity profiles. The simulations revealed that the friction factor decreases with increasing Reynolds number and roughness spacing, and increases strongly with increasing roughness height. The effect of random element distribution on friction factor and mean velocities is however weak. In all cases, there is a clear cut between the inner layer near the wall, which is affected by the presence of the roughness elements, and the outer layer, which remains relatively unaffected. The study reveals that the presence of roughness elements of this shape promotes locally the instantaneous flow motion in the lateral direction in the wall layer, causing a transfer of energy from the streamwise Reynolds stress to the lateral component. The study indicates also that the coherent structures developing in the wall layer are rather similar to the smooth case but are lifted up by almost a constant wall-unit shift y⁺ (∼10–15), which, interestingly, corresponds to the relative roughness k⁺ = 10.     

Influence of gravity and lift on particle velocity statistics and transfer rates in turbulent vertical channel flow

The present study reports detailed statistics for velocity and transfer rates of heavy particles dispersed in turbulent boundary layers. Statistics have been extracted from a homogeneous source of data covering a large target parameter space and are used here to analyze the effects of gravity and lift on particle dispersion and deposition in a systematic way. Datasets were obtained performing Direct Numerical Simulation (DNS) of particle-laden turbulent upward/downward flow in a vertical channel. Six values for the particle timescale (the particle Stokes number, St) ranging three orders of magnitude were considered to analyze the deposition process as the controlling mechanism was shifting from turbulent diffusion to inertia-moderated crossing trajectories. For the particle timescales examined, gravity and lift do not influence the qualitative behavior of particles even though velocity profiles and deposition coefficients are modified in a non-monotonic fashion, reaching an optimum for St ? 15. Physical mechanisms for the different behavior are discussed. Raw data and statistics obtained from the present DNS are made available at http://cfd.cineca.it (mirror site: http://158.110.32.35/download/database) and can be used to test simple models and closure equations for multiphase RANS and Large Eddy simulations.     

Near-wall measurements of turbulence statistics in a fully developed channel flow with a novel laser Doppler velocity profile sensor

This paper reports on the measurements of the near-wall turbulence statistics in a fully developed channel flow. The flow measurements were carried out with a novel laser Doppler velocity profile sensor with a high spatial resolution. The sensor provides both the information of velocity and position of individual tracer particles inside the measurement volume. Hence, it yields the velocity profile inside the measurement volume, in principle, without the sensor being mechanically traversed. Two sensor systems were realized with different techniques. Typically the sensor has a relative accuracy of velocity measurement of 10⁻³ and the spatial resolution of a few micrometers inside the measurement volume of about 500 μm long. The streamwise velocity was measured with two independent sensor systems at three different Reynolds number conditions. The resulting turbulence statistics show a good agreement with available data of direct numerical simulations up to fourth order moment. This demonstrates the velocity profile sensor to be one of the promising techniques for turbulent flow research with the advantage of a spatial resolution more than one magnitude higher than a conventional laser Doppler technique.     

Numerical modeling of fluid flow and heat transfer in a narrow Taylor-Couette-Poiseuille system

We consider turbulent flows in a differentially heated Taylor-Couette system with an axial Poiseuille flow. The numerical approach is based on the Reynolds Stress Modeling (RSM) of [Elena and Schiestel, 1996] and [Schiestel and Elena, 1997] widely validated in various rotor-stator cavities with throughflow ( [Poncet, 2005], [Poncet et al., 2005] and [Haddadi and Poncet, 2008]) and heat transfer (Poncet and Schiestel, 2007). To show the capability of the present code, our numerical predictions are compared very favorably to the velocity measurements of Escudier and Gouldson (1995) in the isothermal case, for both the mean and turbulent fields. The RSM model improves, in particular, the predictions of the k-ε model of Naser (1997). Then, the second order model is applied for a large range of rotational Reynolds (3744 ? Re_i ? 37,443) and Prandtl numbers (0.01 ? Pr ? 12), flow rate coefficient (0 ? C_w ? 30,000) in a very narrow cavity of radius ratio s = R_i/R_o = 0.961 and aspect ratio L = (R_o − R_i)/h = 0.013, where R_i and R_o are the radii of the inner and outer cylinders respectively and h is the cavity height. Temperature gradients are imposed between the incoming fluid and the inner and outer cylinders. The mean hydrodynamic and thermal fields reveal three distinct regions across the radial gap with a central region of almost constant axial and tangential mean velocities and constant mean temperature. Turbulence, which is weakly anisotropic, is mainly concentrated in that region and vanishes towards the cylinders. The mean velocity distributions are not clearly affected by the rotational Reynolds number and the flow rate coefficient. The effects of the flow parameters on the thermal field are more noticeable and considered in details. Correlations for the averaged Nusselt numbers along both cylinders are finally provided according to the flow control parameters Re_i, C_w, and Pr.     

Experimental and numerical investigation of turbulent cross-flow in a staggered tube bundle

This paper presents the results of measurements and numerical predictions of turbulent cross-flow in a staggered tube bundle. The bundle consists of transverse and longitudinal pitch-to-diameter ratios of 3.8 and 2.1, respectively. The experiments were conducted using a particle image velocimetry technique, in a flow of water in a channel at a Reynolds number of 9300 based on the inlet velocity and the tube diameter. A commercial CFD code, ANSYS CFX V10.0, is used to predict the turbulent flow in the bundle. The steady and isothermal Reynolds–Averaged Navier–Stokes (RANS) equations were used to predict the turbulent flow using each of the following four turbulence models: a k-epsilon, a standard k-omega, a k-omega-based shear stress transport, and an epsilon-based second moment closure. The epsilon-based models used a scalable wall function and the omega-based models used a wall treatment that switches automatically between low-Reynolds and standard wall function formulations.

The experimental results revealed extremely high levels of turbulence production by the normal stresses, as well as regions of negative turbulence production. The convective transport by mean flow and turbulent diffusion were observed to be significantly higher than in classical turbulent boundary layers. As a result, turbulence production is generally not in equilibrium with its dissipation rate. In spite of these characteristics, it was observed that the Reynolds normal stresses approximated from the k-based two-equation models were in a closer agreement with experiments than values obtained from the second moment closure. The results show that none of the turbulence models was able to consistently reproduce the mean and turbulent quantities reasonably well. The omega-based models predicted the mean velocities better in the developing region while the epsilon-based models gave better results in the region where the flow is becoming spatially periodic.     

Assessment of the performance of different classes of turbulence models in a wide range of non-equilibrium flows

The performance of thirteen benchmark turbulence models within the RANS framework has been assessed in classical non-equilibrium flows. Linear and non-linear eddy-viscosity schemes, Reynolds stress transport models and single- and two-time-scale approaches have been considered in the investigation. Among the test cases studied are homogeneous shear and normally strained flows, adverse-pressure-gradient, favourable-pressure-gradient and oscillatory boundary layer flows, fully developed oscillatory and ramp up pipe flows and steady and pulsated backward-facing-step flows. The main advantages and drawbacks of the models in each of the test cases are discussed. These discussions provide a reasonably wide understanding of the expected behaviour of the models for future applications in non-equilibrium flows, and also result in suggestions on how the effectiveness of existing models can be further improved.     

Characteristic turbulent structure of a modified drag-reduced surfactant solution flow via dosing water from channel wall

The modification of the near-wall structure is very important for the control of wall turbulence. To ascertain the effect of near-wall modulation on the viscoelastic drag-reduced flow, the modified characteristics of a surfactant solution channel flow were investigated experimentally. The modulation was conducted on the boundary of the channel flow by injecting water from the whole surface of one side of the channel wall. The diffusion process of the injected water was observed by using the planar laser-induced fluorescence technique. The velocity statistics and characteristic structure including the spatial distributions of instantaneous streamwise velocity, swirling strength, and Reynolds shear stress were analyzed based on the velocity vectors acquired in the streamwise wall-normal plane by using the particle imaging velocimetry technique. The results indicated that the disturbance of the injected water was constricted within a finite range very near the dosing wall, and the Reynolds shear stress was increased in this region. However, the eventual drag reduction rate was found to be increased due to a relatively large decrement of viscoelastic shear stress in this near-wall region. Moreover, the flow structure under this modulation presented obvious regional characteristics. In the unstable disturbed region, the mixing of high-speed and low-speed fluids and the motions of ejection and sweep occurred actively. Many clockwise vortex cores were also found to be generated. This characteristic structure was similar to that in the ordinary turbulence of Newtonian fluid. Nevertheless, outside this disturbed region, the structure still maintained the characteristics of the drag-reduced flow with non-Newtonian viscoelastic additives. These results proved that the injected Newtonian fluid associated with the modified stress distribution creates a diverse characteristic structure and subsequent enhanced drag reduction. This investigation can provide the experimental basis for further study of turbulence control.     

The nature of near-wall convection velocity in turbulent channel flow

A novel notion of turbulent structure the local cascade structure-is introduced to study the convection phenomenon in a turbulent channel flow. A space-time cross-correlation method is used to calculate the convection velocity. It is found that there are two characteristic convection speeds near the wall, one associated with small-scale streaks of a lower speed and another with streamwise vortices and hairpin vortices of a higher speed. The new concept of turbulent structure is powerful to illustrate the dominant role of coherent structures in the near-wall convection, and to reveal also the nature of the convection-the propagation of patterns of velocity fluctuations-which is scale-dependent.     

Orthogonal wavelet analysis of counter gradient transport phenomena in turbulent asymmetric channel flow

In this paper four families of orthogonal wavelets are applied to analyze the turbulent counter gradient transport phenomena in fully developed asymmetric channel flows. The results show that: (1) In the instance of counter gradient transport, the principal scale of the coherent structure is responsible for the strong local counter gradient transport; (2) Counter gradient transport phenomena have a strong effect on the intermittency of turbulence; (3) Non-Gaussian part of the principal coherent structure is essential for counter gradient transport phenomena.The project supported by the National Natural Science Foundation of China (10272071, 10472063)     

Direct numerical simulation of free falling sphere in creeping flow

In the present study, direct numerical simulations (DNS) are performed on single and a swarm of particles settling under the action of gravity. The simulations have been carried out in the creeping flow range of Reynolds number from 0.01 to 1 for understanding the hindrance effect, of the other particles, on the settling velocity and drag coefficient. The DNS code is a non-Lagrange multiplier-based fictitious-domain method, which has been developed and validated by Jin et al. (2008 Jin, S. A parallel algorithm for the direct numerical simulation of 3D inertial particle sedimentation. Conference proceedings of the 16th annual conference of the CFD Society of Canada. Edited by: Bergstrom, D. J. and Spiteri, R. 9–11, June. Saskatoon, Saskatchewan, , Canada [Google Scholar]; A parallel algorithm for the direct numerical simulation of 3D inertial particle sedimentation. In: Conference proceedings of the 16th annual conference of the CFD Society of Canada). It has been observed that the time averaged settling velocity of the particle in the presence of other particles, decreases with an increase in the number of particles surrounding it (from 9 particles to 245 particles). The effect of the particle volume fraction on the drag coefficient has also been studied and it has been observed that the computed values of drag coefficients are in good agreement with the correlations proposed by Richardson and Zaki (1954 Richardson, J. F. and Zaki, W. N. 1954. Sedimentation and fluidization: part I. Transactions of the Institution of Chemical Engineers, 32: 35–53.  [Google Scholar]; Sedimentation and fluidization: part I. Transactions of the Institution of Chemical Engineers, 32, 35–53) and Pandit and Joshi (1998 Pandit, A. B. and Joshi, J. B. 1998. Pressure drop in packed, expanded and fluidized beds, packed columns and static mixers – a unified approach. Reviews in Chemical Engineering, 14: 321–371. [Crossref], [Web of Science ®] , [Google Scholar]; Pressure drop in packed, expanded and fluidised beds, packed columns and static mixers – a unified approach. Reviews in Chemical Engineering, 14, 321–371). The suspension viscosity-based model of Frankel and Acrivos (1967 Frankel, N. A. and Acrivos, A. 1967. On the viscosity of a concentrated suspension of solid spheres. Chemical Engineering Science, 22: 847–853. [Crossref], [Web of Science ®] , [Google Scholar]; On the viscosity of a concentrated suspension of solid spheres. Chemical Engineering Science, 22, 847–853) shows good agreement with the DNS results.     

DNS of passive scalar transport in turbulent channel flow at high Schmidt numbers

We perform DNS of passive scalar transport in low Reynolds number turbulent channel flow at Schmidt numbers up to Sc = 49. The high resolutions required to resolve the scalar concentration fields at such Schmidt numbers are achieved by a hierarchical algorithm in which only the scalar fields are solved on the grid dictated by the Batchelor scale. The velocity fields are solved on coarser grids and prolonged by a conservative interpolation to the fine-grid.

The trends observed so far at lower Schmidt numbers Sc  10 are confirmed, i.e. the mean scalar gradient steepens at the wall with increasing Schmidt number, the peaks of turbulent quantities increase and move towards the wall. The instantaneous scalar fields show a dramatic change. Observable structures get longer and thinner which is connected with the occurrence of steeper gradients, but the wall concentrations penetrate less deeply into the plateau in the core of the channel.

Our data shows that the thickness of the conductive sublayer, as defined by the intersection point of the linear with the logarithmic asymptote scales with Sc^−0.29. With this information it is possible to derive an expression for the dimensionless transfer coefficient K⁺ which is only dependent on Sc and Re_τ. This expression is in full accordance to previous results which demonstrates that the thickness of the conductive sublayer is the dominating quantity for the mean scalar profile.     

Numerical simulation of an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder

This article presents a numerical investigation of turbulent flow in an axisymmetric separated and reattached flow over a longitudinal blunt circular cylinder. The governing equations were discretized by the finite-volume method and SIMPLER method was applied to solve the equations on a staggered grid. The turbulent flow was numerically simulated using the standard k–ε, Abe–Kondoh–Nagano (AKN) and Shear Stress Transport (SST) turbulence models. The comparisons made between numerical results and experimental measurements showed that the SST model is superior to other models in the present calculation.Computations were performed for three different Reynolds numbers of 6000, 10 000 and 20 000 based on the cylinder diameter. To our knowledge, this study represents the first numerical investigation of the present flow configuration. The computational results were validated with the available experimental data of reattachment length, mean velocity distribution and wall static pressure coefficient in the turbulent blunt circular cylinder flows. Further, other characteristics of the flow, such as turbulent kinetic energy, pressure, streamlines, and the velocity vectors are discussed.The results show that the main characteristics of the turbulence flow in the separation region, such as reattachment length or velocity profiles, are nearly independent of the Reynolds number. The obtained results showed that a secondary separation bubble may appear in the main separation bubble near the leading edge. Furthermore, it was found that the turbulent kinetic energy has a large effect on the formation of the secondary bubble.     

Analysis of velocity equation of steady flow of a viscous incompressible fluid in channel with porous walls

Steady flow of a viscous incompressible fluid in a channel, driven by suction or injection of the fluid through the channel walls, is investigated. The velocity equation of this problem is reduced to nonlinear ordinary differential equation with two boundary conditions by appropriate transformation and convert the two‐point boundary‐value problem for the similarity function into an initial‐value problem in which the position of the upper channel. Then obtained differential equation is solved analytically using differential transformation method and compare with He's variational iteration method and numerical solution. These methods can be easily extended to other linear and nonlinear equations and so can be found widely applicable in engineering and sciences. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental and high-order LES analysis of the flow in near-wall region of a square cylinder

A coupled experimental/numerical analysis of turbulent flow past a square cylinder is performed at the ERCOFTAC Reynolds number Re = U_∞D/ν = 21,400, where U_∞ is the inflow velocity and D the cylinder height. Complementary Laser Doppler Velocimetry (LDV) and high-order large-eddy simulations (LES) approaches, based on a spectral vanishing technique (SVV-LES), provide a comprehensive data base including both instantaneous data and post-processed statistics. Beyond these results, an achievement of the paper is to investigate the coherent structures developing on the sides and in the wake of the cylinder with a special focus on the flow features in the near-wall region. The flow is found to separate at the leading edge of the cylinder with the occurence of three-dimensional Kelvin-Helmholtz (KH) pairings localized in the separating shear layer. The interaction between these KH vortical structures and Von Kármán vortex shedding (VK) in the near wake is discussed based on both visualisations and frequency analysis. In particular, signatures of VK and KH vortical structures are found on velocity time samples.     

Cobra probe measurements of mean velocities, Reynolds stresses and higher-order velocity correlations in pipe flow

This paper is concerned with the validation of a four-hole pressure probe, known as a cobra probe, for turbulence measurement. For the first time in the literature, third- and fourth-order velocity correlations measured using a pressure probe are presented. The probe measurements are compared with established data for fully developed pipe flow, and good agreement is found. A new probe calibration methodology and improvements to the data acquisition and processing system are also presented.     

Wake response of a stationary finite-sized particle in a turbulent channel flow

Here we consider the effect of a finite-sized stationary particle in a channel flow of modest turbulence at _Reτ=178.12${\mathit{Re}}_{\tau }=178.12$. The size of particle is varied such that the particle Reynolds number ranges from about 40 to 450. The location of the particle is chosen to be either in the buffer layer (_yp+=17.81)$(y_{p+}=17.81)$ or at the channel center. Fully resolved direct numerical simulations of the turbulent channel flow around the particles is performed. Here the ambient turbulence intensity relative to the mean velocity seen by the particle is large (I=23.16%)$(I=23.16\%)$ in the buffer region, while it is substantially lower (I=4.09%)$(I=4.09\%)$ at the channel center. We present results on turbulence modulation due to the particle in terms of wake dynamics and vortex shedding.     

Dynamic instabilities of boiling two-phase flow in a single horizontal channel

Dynamic instabilities of two-phase flow associated with refrigerant R-11 in a uniformly heated horizontal in-tube boiling system were experimentally investigated. An experimental setup was designed and built to work in a wide range of mass fluxes G [75–1050 kg/(m² s)], heat fluxes q (0–100 kW/m²), and fluid inlet temperatures T_inlet (2–24°C). Dynamic instability data were obtained under various working conditions. The dependence of oscillation amplitude and period on system parameters is discussed, and the boundaries of various oscillations are located on the steady-state characteristic curves.     

Statistical structure of the velocity field in cavitating flow around a 2D hydrofoil

Transformation of flow turbulence structure with cavitation occurrence, determination of the flow conditions favorable for nucleation of cavitation bubbles, influence of the statistical structure of turbulence on this process and the inverse effect of cavitation on the flow dynamics are challenging problems in modern fluid mechanics. The paper reports on the results of statistical processing of the velocity fields measured by a PIV technique in cavitating flow over a 2D symmetric hydrofoil for four flow conditions, starting from a cavitation-free regime and finishing by unsteady cloud cavitation. We analyze basic information on the statistical structure of velocity fluctuations in the form of histograms and Q-Q diagrams along with profiles of the mean velocity and turbulent kinetic energy. The research reveals that the flow turbulence pattern and distributions of turbulent fluctuations change significantly with the cavitation development. Under unsteady cloud cavitation conditions, the probability density function of the fluctuating velocity has a two-mode distribution, which indicates switching of two alternating flow conditions in a region above the hydrofoil aft part due to periodic passing of cavitation clouds. Behaviors of the mean and most probable velocities unexpectedly appear to be different with a monotonous increase of the incoming flow velocity. This finding must be caused by modification of the skewness coefficient of the fluctuating velocity.