Recent Progress on MILES for High Reynolds Number Flows

2002 ◽  
Vol 124 (4) ◽  
pp. 848-861 ◽  
Author(s):  
F. F. Grinstein ◽  
C. Fureby

A promising large-eddy simulation (LES) approach is monotonically integrated LES (MILES) which involves solving the Navier-Stokes equations using high-resolution monotone algorithms. In MILES, the subgrid scale (SGS) flow physics is provided by intrinsic, nonlinear, high-frequency filters built into the discretization and implicit SGS models. Mathematical and physical aspects of implicit SGS modeling using nonlinear flux-limiters are addressed using a formalism based on the modified LES equations approach. Detailed properties of the implicit subgrid model are related to the flux limiter, which in turn depends on the specifics of the numerical scheme; we illustrate how the latter properties can directly affect their potential in the MILES framework. Major unresolved issues relevant to LES of complex practical turbulent flows are discussed in this context, including some aspects of boundary condition modeling and overall computational model validation.

2013 ◽  
Vol 721 ◽  
pp. 58-85 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall ◽  
Andrew Walton

AbstractThe recently understood relationship between high-Reynolds-number vortex–wave interaction theory and computationally generated self-sustaining processes provides a possible route to an understanding of some of the underlying structures of fully turbulent flows. Here vortex–wave interaction (VWI) theory is used in the long streamwise wavelength limit to continue the development found at order-one wavelengths by Hall & Sherwin (J. Fluid Mech., vol. 661, 2010, pp. 178–205). The asymptotic description given reduces the Navier–Stokes equations to the so-called boundary-region equations, for which we find equilibrium states describing the change in the VWI as the wavelength of the wave increases from $O(h)$ to $O(Rh)$, where $R$ is the Reynolds number and $2h$ is the depth of the channel. The reduced equations do not include the streamwise pressure gradient of the perturbation or the effect of streamwise diffusion of the wave–vortex states. The solutions we calculate have an asymptotic error proportional to ${R}^{- 2} $ when compared to the full Navier–Stokes equations. The results found correspond to the minimum drag configuration for VWI states and might therefore be of relevance to the control of turbulent flows. The key feature of the new states discussed here is the thickening of the critical layer structure associated with the wave part of the flow to completely fill the channel, so that the roll part of the flow is driven throughout the flow rather than as in Hall & Sherwin as a stress discontinuity across the critical layer. We identify a critical streamwise wavenumber scaling, which, when approached, causes the flow to localize and take on similarities with computationally generated or experimentally observed turbulent spots. In effect, the identification of this critical wavenumber for a given value of the assumed high Reynolds number fixes a minimum box length necessary for the emergence of localized structures. Whereas nonlinear equilibrium states of the Navier–Stokes equations are thought to form a backbone on which turbulent flows hang, our results suggest that the localized states found here might play a related role for turbulent spots.


1997 ◽  
Vol 119 (4) ◽  
pp. 900-905 ◽  
Author(s):  
X. Zheng ◽  
C. Liao ◽  
C. Liu ◽  
C. H. Sung ◽  
T. T. Huang

In this paper, computational results are presented for three-dimensional high-Reynolds number turbulent flows over a simplified submarine model. The simulation is based on the solution of Reynolds-Averaged Navier-Stokes equations and two-equation turbulence models by using a preconditioned time-stepping approach. A multiblock method, in which the block loop is placed in the inner cycle of a multi-grid algorithm, is used to obtain versatility and efficiency. It was found that the calculated body drag, lift, side force coefficients and moments at various angles of attack or angles of drift are in excellent agreement with experimental data. Fast convergence has been achieved for all the cases with large angles of attack and with modest drift angles.


Author(s):  
Sébastien Bourdet ◽  
Marianna Braza ◽  
Yannick Hoarau ◽  
Rajah El Akoury ◽  
Arif Ashraf ◽  
...  

The fluid structure interaction due to the pitching motion of a NACA0012 aerofoil has been studied numerically at moderate and high Reynolds numbers. The dynamic mesh method has been employed in the code ICARE/IMFT solving the Navier-Stokes equations in compressible flows. At high Reynolds number, the phase-averaged Navier-Stokes equations have been solved, coupled with advanced URANS modelling in the NSMB code. The vortex dynamics and especially the stall are physically captured by the dynamic mesh method and by the URANS/Organised Eddy Simulation approach.


2017 ◽  
Vol 814 ◽  
pp. 570-591 ◽  
Author(s):  
Georgios Rigas ◽  
Aimee S. Morgans ◽  
Jonathan F. Morrison

A theory is presented where the weakly nonlinear analysis of laminar globally unstable flows in the presence of external forcing is extended to the turbulent regime. The analysis is demonstrated and validated using experimental results of an axisymmetric bluff-body wake at high Reynolds numbers, $Re_{D}\sim 1.88\times 10^{5}$, where forcing is applied using a zero-net-mass-flux actuator located at the base of the blunt body. In this study we focus on the response of antisymmetric coherent structures with azimuthal wavenumbers $m=\pm 1$ at a frequency $St_{D}=0.2$, responsible for global vortex shedding. We found experimentally that axisymmetric forcing ($m=0$) couples nonlinearly with the global shedding mode when the flow is forced at twice the shedding frequency, resulting in parametric subharmonic resonance through a triadic interaction between forcing and shedding. We derive simple weakly nonlinear models from the phase-averaged Navier–Stokes equations and show that they capture accurately the observed behaviour for this type of forcing. The unknown model coefficients are obtained experimentally by producing harmonic transients. This approach should be applicable in a variety of turbulent flows to describe the response of global modes to forcing.


1993 ◽  
Vol 115 (4) ◽  
pp. 638-645 ◽  
Author(s):  
Hsiao C. Kao

The problem of turbulent flows in two-inlet channels has been studied numerically by solving the Reynolds-averaged Navier-Stokes equations with the k–ε model in a mapped domain. Both the high Reynolds number and the low Reynolds number form were used for this purpose. In general, the former predicts a weaker and smaller recirculation zone than the latter. Comparisons with experimental data, when applicable, were also made. The bulk of the present computations used, however, the high Reynolds number form to correlate different geometries and inflow conditions with the flow properties after turning.


2003 ◽  
Vol 478 ◽  
pp. 227-235 ◽  
Author(s):  
J. D. GIBBON ◽  
Charles R. DOERING

Dissipation-range intermittency was first observed by Batchelor & Townsend (1949) in high Reynolds number turbulent flows. It typically manifests itself in spatio-temporal binary behaviour which is characterized by long, quiescent periods in the signal which are interrupted by short, active ‘events’ during which there are large excursions away from the average. It is shown that Leray's weak solutions of the three-dimensional incompressible Navier–Stokes equations can have this binary character in time. An estimate is given for the widths of the short, active time intervals, which decreases with the Reynolds number. In these ‘bad’ intervals singularities are still possible. However, the average width of a ‘good’ interval, where no singularities are possible, increases with the Reynolds number relative to the average width of a bad interval.


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

This is a report on a field experiment in an atmospheric surface layer at heights between 0.8 and 10m with the Taylor micro-scale Reynolds number in the range Reλ = 1.6−6.6 ×103. Explicit information is obtained on the full set of velocity and temperature derivatives both spatial and temporal, i.e. no use of Taylor hypothesis is made. The report consists of three parts. Part 1 is devoted to the description of facilities, methods and some general results. Certain results are similar to those reported before and give us confidence in both old and new data, since this is the first repetition of this kind of experiment at better data quality. Other results were not obtained before, the typical example being the so-called tear-drop R-Q plot and several others. Part 2 concerns accelerations and related matters. Part 3 is devoted to issues concerning temperature, with the emphasis on joint statistics of temperature and velocity derivatives. The results obtained in this work are similar to those obtained in experiments in laboratory turbulent grid flow and in direct numerical simulations of Navier–Stokes equations at much smaller Reynolds numbers Reλ ~ 102, and this similarity is not only qualitative, but to a large extent quantitative. This is true of such basic processes as enstrophy and strain production, geometrical statistics, the role of concentrated vorticity and strain, reduction of nonlinearity and non-local effects. The present experiments went far beyond the previous ones in two main respects. (i) All the data were obtained without invoking the Taylor hypothesis, and therefore a variety of results on fluid particle accelerations became possible. (ii) Simultaneous measurements of temperature and its gradients with the emphasis on joint statistics of temperature and velocity derivatives. These are reported in Parts 2 and 3.


2009 ◽  
Vol 2009 ◽  
pp. 1-20 ◽  
Author(s):  
Yan Zhang ◽  
Yinnian He

We propose a stabilized subgrid finite-element method for the two-dimensional (2D) nonstationary incompressible Naver-Stokes equation (NSE). This method yields a subgrid eddy viscosity which does not act on the large flow structures. The proposed eddy viscous term is constructed by a fluctuation operator based on an L2-projection. The fluctuation operator can be implemented by the L2-projection from high-order interpolation finite-element spaces to the low-order interpolation finite-element spaces. In this paper, P2/P1 mixed finite-element spaces are adopted to implement the calculation and the analysis. The error analysis is given based on some regular assumptions. Finally, in the part of numerical tests, the numerical computations show that the numerical results agree with theoretical analysis very well. Meanwhile, the numerical investigations demonstrate that the proposed subgrid is very effective for high Reynolds number fluid flows and superior to other referred subgrid models.


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.


Sign in / Sign up

Export Citation Format

Share Document