On the role of return to isotropy in wall-bounded turbulent flows with buoyancy

2018 ◽  
Vol 856 ◽  
pp. 61-78 ◽  
Author(s):  
Elie Bou-Zeid ◽  
Xiang Gao ◽  
Cedrick Ansorge ◽  
Gabriel G. Katul

High Reynolds number wall-bounded turbulent flows subject to buoyancy forces are fraught with complex dynamics originating from the interplay between shear generation of turbulence ($S$) and its production or destruction by density gradients ($B$). For horizontal walls, $S$ augments the energy budget of the streamwise fluctuations, while $B$ influences the energy contained in the vertical fluctuations. Yet, return to isotropy remains a tendency of such flows where pressure–strain interaction redistributes turbulent energy among all three velocity components and thus limits, but cannot fully eliminate, the anisotropy of the velocity fluctuations. A reduced model of this energy redistribution in the inertial (logarithmic) sublayer, with no tuneable constants, is introduced and tested against large eddy and direct numerical simulations under both stable ($B<0$) and unstable ($B>0$) conditions. The model links key transitions in turbulence statistics with flux Richardson number (at $Ri_{f}=-B/S\approx$$-2$, $-1$ and $-0.5$) to shifts in the direction of energy redistribution. Furthermore, when coupled to a linear Rotta-type closure, an extended version of the model can predict individual variance components, as well as the degree of turbulence anisotropy. The extended model indicates a regime transition under stable conditions when $Ri_{f}$ approaches $Ri_{f,max}\approx +0.21$. Buoyant destruction $B$ increases with increasing stabilizing density gradients when $Ri_{f}<Ri_{f,max}$, while at $Ri_{f}\geqslant Ri_{f,max}$ limitations on the redistribution into the vertical component throttle the highest attainable rate of buoyant destruction, explaining the ‘self-preservation’ of turbulence at large positive gradient Richardson numbers. Despite adopting a ‘framework of maximum simplicity’, the model results in novel and insightful findings on how the interacting roles of energy redistribution and buoyancy modulate the variance budgets and the energy exchange among the components.

2007 ◽  
Vol 589 ◽  
pp. 103-123 ◽  
Author(s):  
G. GULITSKI ◽  
M. KHOLMYANSKY ◽  
W. KINZELBACH ◽  
B. LÜTHI ◽  
A. TSINOBER ◽  
...  

This is part 3 of our work describing experiments in which explicit information was obtained on all the derivatives, i.e. spatial derivatives, ∂/∂xj, and temporal derivatives, ∂/∂t, of velocity and temperature fields (and all the components of velocity fluctuations and temperature) at the Reynolds numberReλ~104.This part is devoted to the issues concerning temperature with the emphasis on joint statistics of temperature and velocity derivatives, based on preliminary results from a jet facility and the main results from a field experiment. Apart from a number of conventional results, these contain a variety of results concerning production of temperature gradients, such as role of vorticity and strain, eigen-contributions, geometrical statistics such as alignments of the temperature gradient and the eigenframe of the rate-of-strain tensor, tilting of the temperature gradient, comparison of the true production of the temperature gradient with its surrogate. Among the specific results of importance is the essential difference in the behaviour of the production of temperature gradients in regions dominated by vorticity and strain. Namely, the production of temperature gradients is much more intensive in regions dominated by strain, whereas production of temperature gradients is practically independent of the magnitude of vorticity. In contrast, vorticity and strain are contributing equally to the tilting of the vector of temperature gradients.The production of temperature gradients is mainly due to the fluctuative strain, the terms associated with mean fields are unimportant. It was checked directly (by looking at corresponding eigen-contributions and alignments), that the production of the temperature gradients is due to predominant compressing of fluid elements rather than stretching, which is true of other processes in turbulent flows, e.g. turbulent energy production in shear flows. Though the production of the temperature gradient and its surrogate possess similar univariate PDFs (which indicates the tendency to isotropy in small scales by this particular criterion), their joint PDF is not close to a bisector. This means that the true production of the temperature gradient is far from being fully represented by its surrogate. The main technical achievement is demonstrating the possibility of obtaining experimentally joint statistics of velocity and temperature gradients.


Author(s):  
Michael Leschziner ◽  
Ning Li ◽  
Fabrizio Tessicini

This paper provides a discussion of several aspects of the construction of approaches that combine statistical (Reynolds-averaged Navier–Stokes, RANS) models with large eddy simulation (LES), with the objective of making LES an economically viable method for predicting complex, high Reynolds number turbulent flows. The first part provides a review of alternative approaches, highlighting their rationale and major elements. Next, two particular methods are introduced in greater detail: one based on coupling near-wall RANS models to the outer LES domain on a single contiguous mesh, and the other involving the application of the RANS and LES procedures on separate zones, the former confined to a thin near-wall layer. Examples for their performance are included for channel flow and, in the case of the zonal strategy, for three separated flows. Finally, a discussion of prospects is given, as viewed from the writer's perspective.


Author(s):  
Yan Jin

Abstract The turbulent flow in a compressor cascade is calculated by using a new simulation method, i.e., parameter extension simulation (PES). It is defined as the calculation of a turbulent flow with the help of a reference solution. A special large-eddy simulation (LES) method is developed to calculate the reference solution for PES. Then, the reference solution is extended to approximate the exact solution for the Navier-Stokes equations. The Richardson extrapolation is used to estimate the model error. The compressor cascade is made of NACA0065-009 airfoils. The Reynolds number 3.82 × 105 and the attack angles −2° to 7° are accounted for in the study. The effects of the end-walls, attack angle, and tripping bands on the flow are analyzed. The PES results are compared with the experimental data as well as the LES results using the Smagorinsky, k-equation and WALE subgrid models. The numerical results show that the PES requires a lower mesh resolution than the other LES methods. The details of the flow field including the laminar-turbulence transition can be directly captured from the PES results without introducing any additional model. These characteristics make the PES a potential method for simulating flows in turbomachinery with high Reynolds numbers.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


Author(s):  
Joseph W. Hall ◽  
Charles E. Tinney ◽  
Julie M. Ausseur ◽  
Jeremy T. Pinier ◽  
Andre M. Hall ◽  
...  

2018 ◽  
Vol 28 (5) ◽  
pp. 1096-1116 ◽  
Author(s):  
Emmanuel Leveque ◽  
Hatem Touil ◽  
Satish Malik ◽  
Denis Ricot ◽  
Alois Sengissen

Purpose The Lattice Boltzmann (LB) method offers an alternative to conventional computational fluid dynamics (CFD) methods. However, its practical use for complex turbulent flows of engineering interest is still at an early stage. This paper aims to outline an LB wall-modeled large-eddy simulation (WMLES) solver. Design/methodology/approach The solver is dedicated to complex high-Reynolds flows in the context of WMLES. It relies on an improved LB scheme and can handle complex geometries on multi-resolution block structured grids. Findings Dynamic and acoustic characteristics of a turbulent airflow past a rod-airfoil tandem are examined to test the capabilities of this solver. Detailed direct comparisons are made with both experimental and numerical reference data. Originality/value This study allows assessing the potential of an LB approach for industrial CFD applications.


Author(s):  
Ghazaleh Ahmadi ◽  
Hassan Kassem ◽  
Bernhard Stoevesandt ◽  
Joachim Peinke ◽  
Stefan Heinz

1973 ◽  
Vol 59 (2) ◽  
pp. 281-335 ◽  
Author(s):  
I. J. Wygnanski ◽  
F. H. Champagne

Conditionally sampled hot-wire measurements were taken in a pipe at Reynolds numbers corresponding to the onset of turbulence. The pipe was smooth and carefully aligned so that turbulent slugs appeared naturally atRe> 5 × 104. Transition could be initiated at lowerReby introducing disturbances into the inlet. For smooth or only slightly disturbed inlets, transition occurs as a result of instabilities in the boundary layer long before the flow becomes fully developed in the pipe. This type of transition gives rise to turbulent slugs which occupy the entire cross-section of the pipe, and they grow in length as they proceed downstream. The leading and trailing ‘fronts’ of a turbulent slug are clearly defined. A unique relation seems to exist between the velocity of the interface and the velocity of the fluid by which relaminarization of turbulent fluid is prevented. The length of slugs is of the same order of magnitude as the length of the pipe, although the lengths of individual slugs differ at the same flow conditions. The structure of the flow in the interior of a slug is identical to that in a fully developed turbulent pipe flow. Near the interfaces, where the mean motion changes from a laminar to a turbulent state, the velocity profiles develop inflexions. The total turbulent intensity near the interfaces is very high and it may reach 15% of the velocity at the centre of the pipe. A turbulent energy balance was made for the flow near the interfaces. All of the terms contributing to the energy balance must vanish identically somewhere on the interface if that portion of the interface does not entrain non-turbulent fluid. It appears that diffusion which also includes pressure transport is the most likely mechanism by which turbulent energy can be transferred to non-turbulent fluid. The dissipation term at the interface is negligible and increases with increasing turbulent energy towards the interior of the slug.Mixed laminar and turbulent flows were observed far downstream for\[ 2000 < Re < 2700 \]when a large disturbance was introduced into the inlet. The flow in the vicinity of the inlet, however, was turbulent at much lowerRe. The turbulent regions which are convected downstream at a velocity which is slightly smaller than the average velocity in the pipe we shall henceforth call puffs. The leading front of a puff does not have a clearly defined interface and the trailing front is clearly defined only in the vicinity of the centre-line. The length and structure of the puff is independent of the character of the obstruction which created it, provided that the latter is big enough to produce turbulent flow at the inlet. The puff will be discussed in more detail later.


2013 ◽  
Vol 718 ◽  
pp. 1-4 ◽  
Author(s):  
B. J. McKeon

AbstractMarusic et al. (J. Fluid Mech., vol. 716, 2013, R3) show the first clear evidence of universal logarithmic scaling emerging naturally (and simultaneously) in the mean velocity and the intensity of the streamwise velocity fluctuations about that mean in canonical turbulent flows near walls. These observations represent a significant advance in understanding of the behaviour of wall turbulence at high Reynolds number, but perhaps the most exciting implication of the experimental results lies in the agreement with the predictions of such scaling from a model introduced by Townsend (J. Fluid Mech., vol. 11, 1961, pp. 97–120), commonly termed the attached eddy hypothesis. The elegantly simple, yet powerful, study by Marusic et al. should spark further investigation of the behaviour of all fluctuating velocity components at high Reynolds numbers and the outstanding predictions of the attached eddy hypothesis.


Sign in / Sign up

Export Citation Format

Share Document