Wavelet-based adaptive unsteady Reynolds-averaged turbulence modelling of external flows

2018 ◽  
Vol 837 ◽  
pp. 765-787 ◽  
Author(s):  
Giuliano De Stefano ◽  
Oleg V. Vasilyev ◽  
Eric Brown-Dymkoski

The recent development of the adaptive-anisotropic wavelet-collocation method, which incorporates the use of coordinate transforms, opens new horizons for wavelet-based simulations of wall-bounded turbulent flows. The new wavelet-based adaptive unsteady Reynolds-averaged Navier–Stokes approach for computational modelling of turbulent flows is presented. The proposed methodology that is integrated with anisotropic wavelet-based mesh refinement is demonstrated for a two-equation eddy-viscosity turbulence model. The performance of the method is assessed by conducting numerical simulations of the turbulent flow past a circular cylinder at subcritical Reynolds number. The present study demonstrates both the feasibility and the effectiveness of the new wavelet-based adaptive unsteady Reynolds-averaged turbulence modelling procedure for external flows.

Author(s):  
Charles G. Speziale

The high-Reynolds-number turbulent flows of technological importance contain such a wide range of excited length and time scales that the application of direct or large-eddy simulations is all but impossible for the foreseeable future. Reynolds stress models remain the only viable means for the solution of these complex turbulent flows. It is widely believed that Reynolds stress models are completely ad hoc, having no formal connection with solutions of the full Navier-Stokes equations for turbulent flows. While this belief is largely warranted for the older eddy viscosity models of turbulence, it constitutes a far too pessimistic assessment of the current generation of Reynolds stress closures. It will be shown how secondorder closure models and two-equation models with an anisotropic eddy viscosity can be systematically derived from the Navier-Stokes equations when one overriding assumption is made: the turbulence is locally homogeneous and in equilibrium. A brief review of zero equation models and one equation models based on the Boussinesq eddy viscosity hypothesis will first be provided in order to gain a perspective on the earlier approaches to Reynolds stress modeling. It will, however, be argued that since turbulent flows contain length and time scales that change dramatically from one flow configuration to the next, two-equation models constitute the minimum level of closure that is physically acceptable. Typically, modeled transport equations are solved for the turbulent kinetic energy and dissipation rate from which the turbulent length and time scales are built up; this obviates the need to specify these scales in an ad hoc fashion. While two-equation models represent the minimum acceptable closure, second-order closure models constitute the most complex level of closure that is currently feasible from a computational standpoint. It will be shown how the former models follow from the latter in the equilibrium limit of homogeneous turbulence. However, the two-equation models that are formally consistent with second-order closures have an anisotropic eddy viscosity with strain-dependent coefficients - a feature that most of the commonly used models do not possess.


Author(s):  
C. Henoch ◽  
Martin Hoffert ◽  
A. Baron ◽  
D. Klaiman ◽  
Semion Sukoriansky ◽  
...  

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.


SeMA Journal ◽  
2012 ◽  
Vol 60 (1) ◽  
pp. 51-74
Author(s):  
Christine Bernardi ◽  
Tomás Chacón Rebollo ◽  
Macarena Gómez Mármol

2011 ◽  
Vol 133 (5) ◽  
Author(s):  
F. P. P. Tan ◽  
N. B. Wood ◽  
G. Tabor ◽  
X. Y. Xu

In this study, two different turbulence methodologies are investigated to predict transitional flow in a 75% stenosed axisymmetric experimental arterial model and in a slightly modified version of the model with an eccentric stenosis. Large eddy simulation (LES) and Reynolds-averaged Navier–Stokes (RANS) methods were applied; in the LES simulations eddy viscosity subgrid-scale models were employed (basic and dynamic Smagorinsky) while the RANS method involved the correlation-based transitional version of the hybrid k-ε/k-ω flow model. The RANS simulations used 410,000 and 820,000 element meshes for the axisymmetric and eccentric stenoses, respectively, with y+ less than 2 viscous wall units for the boundary elements, while the LES used 1,200,000 elements with y+ less than 1. Implicit filtering was used for LES, giving an overlap between the resolved and modeled eddies, ensuring accurate treatment of near wall turbulence structures. Flow analysis was carried out in terms of vorticity and eddy viscosity magnitudes, velocity, and turbulence intensity profiles and the results were compared both with established experimental data and with available direct numerical simulations (DNSs) from the literature. The simulation results demonstrated that the dynamic Smagorinsky LES and RANS transitional model predicted fairly comparable velocity and turbulence intensity profiles with the experimental data, although the dynamic Smagorinsky model gave the best overall agreement. The present study demonstrated the power of LES methods, although they were computationally more costly, and added further evidence of the promise of the RANS transition model used here, previously tested in pulsatile flow on a similar model. Both dynamic Smagorinsky LES and the RANS model captured the complex transition phenomena under physiological Reynolds numbers in steady flow, including separation and reattachment. In this respect, LES with dynamic Smagorinsky appeared more successful than DNS in replicating the axisymmetric experimental results, although inflow conditions, which are subject to caveats, may have differed. For the eccentric stenosis, LES with Smagorinsky coefficient of 0.13 gave the closest agreement with DNS despite the known shortcomings of fixed coefficients. The relaminarization as the flow escaped the influence of the stenosis was amply demonstrated in the simulations, graphically so in the case of LES.


2011 ◽  
Vol 21 (03) ◽  
pp. 421-457 ◽  
Author(s):  
RAPHAËL DANCHIN ◽  
MARIUS PAICU

Models with a vanishing anisotropic viscosity in the vertical direction are of relevance for the study of turbulent flows in geophysics. This motivates us to study the two-dimensional Boussinesq system with horizontal viscosity in only one equation. In this paper, we focus on the global existence issue for possibly large initial data. We first examine the case where the Navier–Stokes equation with no vertical viscosity is coupled with a transport equation. Second, we consider a coupling between the classical two-dimensional incompressible Euler equation and a transport–diffusion equation with diffusion in the horizontal direction only. For both systems, we construct global weak solutions à la Leray and strong unique solutions for more regular data. Our results rest on the fact that the diffusion acts perpendicularly to the buoyancy force.


2015 ◽  
Vol 766 ◽  
pp. 337-367 ◽  
Author(s):  
Bartosz Protas ◽  
Bernd R. Noack ◽  
Jan Östh

AbstractWe propose a variational approach to the identification of an optimal nonlinear eddy viscosity as a subscale turbulence representation for proper orthogonal decomposition (POD) models. The ansatz for the eddy viscosity is given in terms of an arbitrary function of the resolved fluctuation energy. This function is found as a minimizer of a cost functional measuring the difference between the target data coming from a resolved direct or large-eddy simulation of the flow and its reconstruction based on the POD model. The optimization is performed with a data-assimilation approach generalizing the 4D-VAR method. POD models with optimal eddy viscosities are presented for a 2D incompressible mixing layer at $\mathit{Re}=500$ (based on the initial vorticity thickness and the velocity of the high-speed stream) and a 3D Ahmed body wake at $\mathit{Re}=300\,000$ (based on the body height and the free-stream velocity). The variational optimization formulation elucidates a number of interesting physical insights concerning the eddy-viscosity ansatz used. The 20-dimensional model of the mixing-layer reveals a negative eddy-viscosity regime at low fluctuation levels which improves the transient times towards the attractor. The 100-dimensional wake model yields more accurate energy distributions as compared to the nonlinear modal eddy-viscosity benchmark proposed recently by Östh et al. (J. Fluid Mech., vol. 747, 2014, pp. 518–544). Our methodology can be applied to construct quite arbitrary closure relations and, more generally, constitutive relations optimizing statistical properties of a broad class of reduced-order models.


1994 ◽  
Vol 47 (6S) ◽  
pp. S3-S13 ◽  
Author(s):  
Parviz Moin ◽  
Thomas Bewley

A brief review of current approaches to active feedback control of the fluctuations arising in turbulent flows is presented, emphasizing the mathematical techniques involved. Active feedback control schemes are categorized and compared by examining the extent to which they are based on the governing flow equations. These schemes are broken down into the following categories: adaptive schemes, schemes based on heuristic physical arguments, schemes based on a dynamical systems approach, and schemes based on optimal control theory applied directly to the Navier-Stokes equations. Recent advances in methods of implementing small scale flow control ideas are also reviewed.


2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi

The direct numerical simulation of turbulence (DNS) has become a method of outmost importance for the investigation of turbulence physics, and its relevance is constantly growing due to the increasing popularity of high-performance-computing techniques. In the present work, the DNS approach is discussed mainly with regard to turbulent shear flows of incompressible fluids with constant properties. A body of literature is reviewed, dealing with the numerical integration of the Navier-Stokes equations, results obtained from the simulations, and appropriate use of the numerical databases for a better understanding of turbulence physics. Overall, it appears that high-performance computing is the only way to advance in turbulence research through the front of the direct numerical simulation.


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.


Sign in / Sign up

Export Citation Format

Share Document