Detached eddy simulation of flow around two wall-mounted cubes in tandem

We investigate the performance of unsteady Reynolds-averaged Navier–Stokes (URANS) computation and various versions of detached eddy simulation (DES) in resolving coherent structures in turbulent flow around two cubes mounted in tandem on a flat plate at Reynolds number (Re) of 22,000 and for a thin incoming boundary layer. Calculations are carried out using four different coherent structure resolving turbulence models: (1) URANS with the Spalart–Allmaras model; (2) the standard DES [Spalart, P.R., Jou, W.H., Strelets, M., Allmaras, S.R., 1997. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In: Liu, C., Liu, Z., (Eds.), Advances in DNS/LES. Greyden Press, Columbus, OH]; (3) the Delayed DES (DDES); and (4) the DES with a low-Re modification (DES-LR) [Spalart, P., Deck, S., Shur, M., Squires, K., Strelets, M., Travin, A., 2006. A new version of detached eddy simulation, resistant to ambiguous grid densities. Theor. Comput. Fluid Dyn. 20 (3), 181–195]. The grid sensitivity of the computed solutions is examined by carrying out simulations on two successively refined grids. The computed results for all cases are compared with the experimental measurements of Martinuzzi and Havel [Martinuzzi, R., Havel, B., 2000. Turbulent flow around two interfering surface-mounted cubic obstacles in tandem arrangement. ASME J. Fluids Eng. 122, 24–31] for two different cube spacings. All turbulence models reproduce essentially identical separation of the approach thin boundary layer and yield an unsteady horseshoe vortex system consisting of multiple vortices in the leading edge region of the upstream cube. Significant discrepancies between the URANS and all DES solutions are observed, however, in other regions of interest such as the shear layers emanating from the cubes, the inter-cube gap and the downstream wake. Regardless of the grid refinement, URANS fails to capture key features of the mean flow, including the second horseshoe vortex in the upstream junction and recirculating flow on the top surface of the downstream cube for the large cube spacing, and underestimates significantly turbulence statistics in most regions of the flow for both cases. On the coarse mesh, all three DES approaches appear to yield very similar results and fail to reproduce the second horseshoe vortex. The standard DES and DDES solutions obtained on the fine meshes are essentially identical and both suffer from premature switching to unresolved DNS, due to the mis-interpretation of grid refinement as wall proximity, which leads to spurious vortices in the inter-cube region. Numerical solutions show that the low-Re modification (DES-LR) is critical prerequisite in DES on the ambiguously fine – not fine enough for full LES – mesh to prevent excessive nonlinear drop of the subgrid eddy viscosity in low cell-Re regions like in the inter-obstacle gap. Mean flow quantities and turbulence statistics obtained with DES-LR on the fine mesh are in good overall agreement with the measurements in most regions of interest for both cases.     

Separated flow structures around a cylindrical obstacle in a narrow channel

A detailed experimental study is performed on the separated flow structures around a low aspect-ratio circular cylinder (pin-fin) in a practical configuration of liquid cooling channel. Distinctive features of the present arrangement are the confinement of the cylinder at both ends, water flow at low Reynolds numbers (Re = 800, 1800, 2800), very high core flow turbulence and undeveloped boundary layers at the position of the obstacle. The horseshoe vortex system at the junctions between the cylinder and the confining walls and the near wake region behind the obstacle are deeply investigated by means of Particle Image Velocimetry (PIV). Upstream of the cylinder, the horseshoe vortex system turns out to be perturbed by vorticity bursts from the incoming boundary layers, leading to aperiodical vortex oscillations at Re = 800 or to break-away and secondary vorticity eruptions at the higher Reynolds numbers. The flow structures in the near wake show a complex three-dimensional behaviour associated with a peculiar mechanism of spanwise mass transport. High levels of free-stream turbulence trigger an early instabilization of the shear layers and strong Bloor–Gerrard vortices are observed even at Re = 800. Coalescence of these vortices and intense spanwise flow inhibit the alternate primary vortex shedding for time periods whose length and frequency increase as the Reynolds number is reduced. The inhibition of alternate vortex shedding for long time periods is finally related to the very large wake characteristic lengths and to the low velocity fluctuations observed especially at the lowest Reynolds number.     

Flow regimes in a single dimple on the channel surface

The boundaries of the domains of existence of flow regimes past single dimples made as spherical segments on a flat plate are determined with the use of available experimental results. Regimes of a diffuser-confuser flow, a horseshoe vortex, and a tornado-like vortex in the dimple are considered. Neither a horseshoe vortex nor a tornado-like vortex is observed in dimples with the relative depth smaller than 0.1. Transformations from the diffuser–confuser flow regime to the horseshoe vortex regime and from the horseshoe vortex flow to the tornado-like vortex flow are found to depend not only on the Reynolds number, but also on the relative depth of the spherical segment. Dependences for determining the boundaries of the regime existence domains are proposed, and parameters at which the experimental results can be generalized are given.     

DNS of heat transfer increase in a cylinder stagnation region due to wake-induced turbulence

Heat transfer in the stagnation point area of a heated cylinder is investigated using Direct Numerical Simulation (DNS). The heated cylinder is subjected to the turbulent wake of a smaller cylinder placed upstream. Two Reynolds numbers based on the diameter of the heated cylinder of 13,200 and 48,000 are chosen. In accordance with correlations in the literature, an increase in heat transfer compared to fully laminar flow is found for all angles along the front circumferential area of the heated cylinder. However, due to the presence of the wake, the maximum increase is shifted away from the centerline. The characteristic turbulence level and Nusselt number in the present study are an order of magnitude higher than those reported in previous simulations. The DNS results obtained, are in good agreement with an existing experimental correlation. Finally, relevant flow structures and instantaneous temperature fields are visualized.     

Assessment of stretched vortex subgrid‐scale models for LES of incompressible inhomogeneous turbulent flow

The physical space version of the stretched vortex subgrid scale model is tested in LES of the turbulent lid‐driven cubic cavity flow. LES is carried out by using a higher order finite‐difference method. The effects of different vortex orientation models and subgrid turbulence spectrums are assessed through comparisons of the LES predictions against DNS. Three Reynolds numbers 12000, 18000, and 22000 are studied. Good agreement with the DNS data for the mean and fluctuating quantities is observed. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessment of dual plane PIV measurements in wall turbulence using DNS data

Experimental dual plane particle image velocimetry (PIV) data are assessed using direct numerical simulation (DNS) data of a similar flow with the aim of studying the effect of averaging within the interrogation window. The primary reason for the use of dual plane PIV is that the entire velocity gradient tensor and hence the full vorticity vector can be obtained. One limitation of PIV is the limit on dynamic range, while DNS is typically limited by the Reynolds number of the flow. In this study, the DNS data are resolved more finely than the PIV data, and an averaging scheme is implemented on the DNS data of similar Reynolds number to compare the effects of averaging inherent to the present PIV technique. The effects of averaging on the RMS values of the velocity and vorticity are analyzed in order to estimate the percentage of turbulence intensity and enstrophy captured for a given PIV resolution in turbulent boundary layers. The focus is also to identify vortex core angle distributions, for which the two-dimensional and three-dimensional swirl strengths are used. The studies are performed in the logarithmic region of a turbulent boundary layer at z ⁺ = 110 from the wall. The dual plane PIV data are measured in a zero pressure gradient flow over a flat plate at Re _τ = 1,160, while the DNS data are extracted from a channel flow at Re _τ = 934. Representative plots at various wall-normal locations for the RMS values of velocity and vorticity indicate the attenuation of the variance with increasing filter size. Further, the effect of averaging on the vortex core angle statistics is negligible when compared with the raw DNS data. These results indicate that the present PIV technique is an accurate and reliable method for the purposes of statistical analysis and identification of vortex structures.     

LES-based filter-matrix lattice Boltzmann model for simulating fully developed turbulent channel flow

In this paper, a three-dimensional filter-matrix lattice Boltzmann (FMLB) model based on large eddy simulation (LES) was verified for simulating wall-bounded turbulent flows. The Vreman subgrid-scale model was employed in the present FMLB–LES framework, which had been proved to be capable of predicting turbulent near-wall region accurately. The fully developed turbulent channel flows were performed at a friction Reynolds number Re_τ of 180. The turbulence statistics computed from the present FMLB–LES simulations, including mean stream velocity profile, Reynolds stress profile and root-mean-square velocity fluctuations greed well with the LES results of multiple-relaxation-time (MRT) LB model, and some discrepancies in comparison with those direct numerical simulation (DNS) data of Kim et al. was also observed due to the relatively low grid resolution. Moreover, to investigate the influence of grid resolution on the present LES simulation, a DNS simulation on a finer gird was also implemented by present FMLB–D3Q19 model. Comparisons of detailed computed various turbulence statistics with available benchmark data of DNS showed quite well agreement.     

Surface pressure on a cubic building exerted by conical vortices

Experiments were performed to study surface pressure on a cubic building underlying conical vortices, which are known to cause severe structural damage and failure. The focus is on the effects of turbulence in the incident flow. Three turbulent boundary layers were created in a boundary layer wind tunnel. A wall-mounted cube, i.e. a cube situated on the horizontal ground floor surface of the wind-tunnel test section, was used as an experimental model. The cube was subjected to the incidence flow at 40°. Steady and unsteady pressure measurements were performed on the cube surface. The analysis suggests that conical vortices developed above the top surface of the wall-mounted cube. A larger mean suction was observed on the top cube surface in the less turbulent boundary layer. With an increase in turbulence in the incoming flow, the strong suction zones decreased in size. The fluctuating pressure coefficient profiles retained their shape when the turbulence in the upstream flow of the cube increased. The fluctuating pressure coefficient was observed to be larger in more turbulent flows. The pressure fluctuations were larger on the cube surface underlying outer boundaries of the conical vortex. The fluctuating pressure coefficient under the conical vortex was three to four times larger than in the weak suction zone on the central area of the top cube surface. Close to the leading cube corner, the pressure spectra were dominated by a single low frequency peak. As the conical vortex developed, this primary peak weakened and a secondary peak emerged at a higher reduced frequency. There is a general trend of shifting the pressure spectra towards higher reduced frequencies when the turbulence in the undisturbed incident flow increases.     

Direct numerical simulation of turbulence over resolved and modeled rough walls with irregularly distributed roughness

To simulate turbulent flow over a rough wall without resolving complicated rough geometries, a macroscopic rough wall model is developed based on spatial (plane) averaging theory. The plane-averaged drag force term, which arises through averaging the Navier–Stokes equations in a plane parallel to a rough wall, can be modeled using a plane porosity and a plane hydraulic diameter. To evaluate the developed model, direct and macroscopic model simulations for turbulence over irregularly distributed semi-spheres at Reynolds number of 300 are carried out using the D3Q27 multiple-relaxation time lattice Boltzmann method. The results show that the developed model can be used to predict rough wall skin friction. The results agree quantitatively with standard turbulence statistics such as mean velocity and Reynolds stress profiles with the fully resolved DNS data. Since velocity dispersion occurs inside the rough wall and is found to contribute to turbulence energy dissipation, which the developed model cannot account for, the developed model fails to reproduce dispersion-related turbulence energy dissipation. However, it is found that the plane-averaged drag force term can successfully recover the deficiency of dispersion-related turbulence energy dissipation.     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

应用小干扰柱体控制角区马蹄涡结构的实验研究

本文提出一种简单的抑制和控制角区马蹄涡的被动控制方法.即在角区平板上游放置一个远小于主柱体的小干扰柱体,用其产生的弱马蹄涡来抑制和控制角区的马蹄涡结构.目的是使角区原来的马蹄涡结构由强变为弱、由大变为小、由多变为少、由非定常变为定常,以获得减小冲刷、抑制湍流、降低噪声、避免振动的工程效果.作者在风洞中采用烟线法和激光片光流动显示的方法开展研究,实验表明,在平板上游适当位置放置小干扰柱体的确可以有效抑制和控制角区马蹄涡结构.实验发现,当小干扰柱体放置在原角区马蹄涡生成区时,其抑制和控制效果最佳;当小干扰柱体放置在上游区或下游区时,控制效果不好.本文讨论了小干扰柱体控制角区马蹄涡的机理.此外,实验还研究了小干扰柱体相对尺度和截面形状对角区马蹄涡结构抑制和控制的影响.     

Recent progress in the development of the Elliptic Blending Reynolds-stress model

The Elliptic Blending Reynolds Stress Model (EB-RSM), originally proposed by Manceau and Hanjalić (2002) to extend standard, weakly inhomogeneous Reynolds stress models to the near-wall region, has been subject to various modifications by several authors during the last decade, mainly for numerical robustness reasons. The present work revisits all these modifications from the theoretical standpoint and investigates in detail their influence on the reproduction of the physical mechanisms at the origin of the influence of the wall on turbulence. The analysis exploits recent DNS databases for high-Reynolds number channel flows, spanwise rotating channel flows with strong rotation rates, up to complete laminarization, and the separated flow after a sudden expansion without and with system rotation. Theoretical arguments and comparison with DNS results lead to the selection of a recommended formulation for the EB-RSM model. This formulation shows satisfactory predictions for the configurations described above, in particular as regards the modification of the mean flow and turbulent anisotropy on the anticyclonic or pressure side.     

An experimental investigation of the flow structure over a corrugated waveform in a transpired air collector

An experimental investigation of the flow dynamics in a channel with a corrugated surface is presented. Particle image velocimetry was used to obtain two-dimensional velocity fields at three different locations along the channel length, over a range of Reynolds numbers. The results show a significant impact of the corrugation waveform on the mean and turbulent flow structure inside the channel. Strong bursting flow originating from the trough, sweeping flow from the bulk region and the vortex shedding off the crest were observed. Their interactions created a complex three-dimensional flow structure extended over almost the entire channel. The mean velocity profiles indicate a strong diffusion of shear. The profiles of various turbulent properties show the enhancement of turbulence in the vicinity of the waveform. It was found that the turbulence in the channel was almost entirely produced in this region above the corrugation trough. Significant momentum transfer from the corrugation wall by the turbulent velocity field was also observed. The mean and turbulent flow behaviour was found to be periodic with respect to the waveform over most of the channel length. The results show the presence of strong turbulence even at the Reynolds number that falls within the conventional laminar range.     

Large eddy simulation and a simple wall model for turbulent flow calculation by a particle method

A turbulent channel flow and the flow around a cubic obstacle are calculated by the moving particle semi‐implicit method with the subparticle‐scale turbulent model and a wall model, which is based on the zero equation RANS (Reynolds Averaged Navier‐Stokes). The wall model is useful in practical problems that often involve high Reynolds numbers and wall turbulence, because it is difficult to keep high resolution in the near‐wall region in particle simulation. A turbulent channel flow is calculated by the present method to validate our wall model. The mean velocity distribution agrees with the log‐law velocity profile near the wall. Statistical values are also the same order and tendency as experimental results with emulating viscous layer by the wall model. We also investigated the influence of numerical oscillations on turbulence analysis in using the moving particle semi‐implicit method. Finally, the turbulent flow around a cubic obstacle is calculated by the present method to demonstrate capability of calculating practical turbulent flows. Three characteristic eddies appear in front of, over, and in the back of the cube both in our calculation and the experimental result that was obtained by Martinuzzi and Tropea. Mean velocity and turbulent intensity profiles are predicted in the same order and have similar tendency as the experimental result. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical study of the turbulent channel flow under space-dependent electromagnetic force control at different Reynolds numbers

In this paper, the control of turbulent channel flow by space-dependent electromagnetic force and the mechanism of drag reduction are investigated with the direct numerical simulation(DNS) methods for different Reynolds numbers. A formulation is derived to express the relation between the drag and the Reynolds shear stress. With the application of optimal electromagnetic force, the in-depth relations among characteristic structures in the flow field, mean Reynolds shear stress, and the effect of drag reduction for different Reynolds numbers are discussed. The results indicate that the maximum drag reductions can be obtained with an optimal combination of parameters for each case of different Reynolds numbers. The regular quasi-streamwise vortex structures, which appear in the flow field, have the same period with that of the electromagnetic force.These structures suppress the random velocity fluctuations, which leads to the absolute value of mean Reynolds shear stress decreasing and the distribution of that moving away from the wall. Moreover, the wave number of optimal electromagnetic force increases,and the scale of the regular quasi-streamwise vortex structures decreases as the Reynolds number increases. Therefore, the rate of drag reduction decreases with the increase in the Reynolds number since the scale of the regular quasi-streamwise vortex structures decreases.     

Large eddy simulation of gas-particle turbulent channel flow with momentum exchange between the phases

This paper presents results of a large eddy simulation (LES) combined with Lagrangian particle tracking and a point-force approximation for the feedback effect of particles on the downward turbulent gaseous flow in a vertical channel. The LES predictions are compared with the results obtained by direct numerical simulation (DNS) of a finer computational mesh. A parametric study is conducted for particles with two response times in simulations with and without streamwise gravitational settling and elastic, binary interparticle collisions. It is shown that the classical and the dynamic Smagorinsky turbulence models adequately predict the particle-induced changes in the mean streamwise velocity and the Reynolds stresses of the carrier phase for the range of parameters studied. However, the largest discrepancies between the LES and DNS results are found in the cases of particle-laden flows. Conditional sampling of the instantaneous resolved flow fields indicates that the mechanisms by which particles directly oppose the production of momentum and vorticity of the organized fluid motions are also observed in the LES results. However, the geometric features of the near-wall quasistreamwise vortices are overestimated by the use of both turbulence models compared to the DNS predictions.     

Numerical simulation for rotating internal weakly viscoelastic flows in rectangular ducts

The present work develops a numerical method for the solution of rotating internal weakly viscoelastic flows in rectangular ducts for dimensionless parameters such as the Reynolds, Rossby and Weissenberg numbers, taken respectively in the intervals between 171 and 12000, 0.047 and 1/12 and up to 1/10000. It is shown that the usual counter‐rotating double‐vortex configuration of secondary flow breaks down with the increase of the Reynolds number (over the threshold of 171). For higher Reynolds numbers such as 7500 and 12000 the secondary flow diffuses to the interior of the duct where it assumes a fully developed configuration and the transition to the turbulence structure is observed. The Sobolev norms increase almost proportionally to the increase of the Reynolds number, and play an essential role for more complex problems involving transition to turbulence modelling. Copyright © 2002 John Wiley & Sons, Ltd.     

Dissipation of Turbulent Kinetic Energy in a Cylinder Wall Junction Flow

The subject of this study is the discussion of the dissipation of turbulent kinetic energy and its Reynolds number scaling in front of a wall-mounted cylinder. We employed highly resolved Large-Eddy Simulation and ensured that the computational grid was fine enough to resolve most of the scales. A perceptible fraction of the total dissipation is modeled. However, this fraction - about one third - is small enough so that the total dissipation suffers only marginally from some potential shortcomings of the turbulence model. Individual terms of the pseudo dissipation tensor and their Reynolds number scaling are discussed and compared. This tensor and thus the turbulent small scale structures are not isotropic at the Reynolds numbers investigated. Furthermore, the near-wall anisotropy under the horseshoe vortex is likely to persist to larger Reynolds numbers as it can be linked to a flapping of the near-wall layer. The turbulent length scale shows a strong spatial variability. In the region of the vortex system in the cylinder front, the distribution reveals a similar shape as the one of the turbulent kinetic energy and its amplitude is in the order of magnitude of the cylinder diameter. In contrast to the region dominated by the approach flow, the turbulent length scale is independent of the Reynolds number in the region dominated by the vortex system. Even though the flow investigated is in non-equilibrium, common a priori estimations and scalings of the Kolmogorov length scale based on macro scales give satisfying results.     

Shear‐improved Smagorinsky modeling of turbulent channel flow using generalized Lattice Boltzmann equation

Generalized Lattice Boltzmann equation (GLBE) was used for computation of turbulent channel flow for which large eddy simulation (LES) was employed as a turbulence model. The subgrid‐scale turbulence effects were simulated through a shear‐improved Smagorinsky model (SISM), which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. The results of mean velocity distribution across the channel showed good correspondence with direct numerical simulation (DNS) data. Negligible discrepancies were observed between the present computations and those reported from DNS for the computed turbulent statistics. Three‐dimensional instantaneous vorticity contours showed complex vortical structures that appeared in such flow geometries. It was concluded that such a framework is capable of predicting accurate results for turbulent channel flow without adding significant complications and the computational cost to the standard Smagorinsky model. As this modeling was entirely local in space it was therefore adapted for parallelization. Copyright © 2010 John Wiley & Sons, Ltd.