On Direct Numerical Simulation of Turbulent Flows

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.

2010 ◽  
Vol 663 ◽  
pp. 100-119 ◽  
Author(s):  
BINH K. LIEU ◽  
RASHAD MOARREF ◽  
MIHAILO R. JOVANOVIĆ

This study builds on and confirms the theoretical findings of Part 1 of this paper (Moarref & Jovanović, J. Fluid Mech., 2010, doi:10.1017/S0022112010003393). We use direct numerical simulation of the Navier–Stokes equations to assess the efficacy of blowing and suction in the form of streamwise travelling waves for controlling the onset of turbulence in a channel flow. We highlight the effects of the modified base flow on the dynamics of velocity fluctuations and net power balance. Our simulations verify the theoretical predictions of Part 1 that the upstream travelling waves promote turbulence even when the uncontrolled flow stays laminar. On the other hand, the downstream travelling waves with parameters selected in Part 1 are capable of reducing the fluctuations' kinetic energy, thereby maintaining the laminar flow. In flows driven by a fixed pressure gradient, a positive net efficiency as large as 25 % relative to the uncontrolled turbulent flow can be achieved with downstream waves. Furthermore, we show that these waves can also relaminarize fully developed turbulent flows at low Reynolds numbers. We conclude that the theory developed in Part 1 for the linearized flow equations with uncertainty has considerable ability to predict full-scale phenomena.


2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès

The linear instability of a vortex ring with swirl with Gaussian distributions of azimuthal vorticity and velocity in its core is studied by direct numerical simulation. The numerical study is carried out in two steps: first, an axisymmetric simulation of the Navier–Stokes equations is performed to obtain the quasi-steady state that forms a base flow; then, the equations are linearized around this base flow and integrated for a sufficiently long time to obtain the characteristics of the most unstable mode. It is shown that the vortex rings are subjected to curvature instability as predicted analytically by Blanco-Rodríguez & Le Dizès (J. Fluid Mech., vol. 814, 2017, pp. 397–415). Both the structure and the growth rate of the unstable modes obtained numerically are in good agreement with the analytical results. However, a small overestimation (e.g. 22 % for a curvature instability mode) by the theory of the numerical growth rate is found for some instability modes. This is most likely due to evaluation of the critical layer damping which is performed for the waves on axisymmetric line vortices in the analysis. The actual position of the critical layer is affected by deformation of the core due to the curvature effect; as a result, the damping rate changes since it is sensitive to the position of the critical layer. Competition between the curvature and elliptic instabilities is also investigated. Without swirl, only the elliptic instability is observed in agreement with previous numerical and experimental results. In the presence of swirl, sharp bands of both curvature and elliptic instabilities are obtained for $\unicode[STIX]{x1D700}=a/R=0.1$, where $a$ is the vortex core radius and $R$ the ring radius, while the elliptic instability dominates for $\unicode[STIX]{x1D700}=0.18$. New types of instability mode are also obtained: a special curvature mode composed of three waves is observed and spiral modes that do not seem to be related to any wave resonance. The curvature instability is also confirmed by direct numerical simulation of the full Navier–Stokes equations. Weakly nonlinear saturation and subsequent decay of the curvature instability are also observed.


Author(s):  
Kun Yang ◽  
Liang Cheng ◽  
Hongwei An ◽  
Ming Zhao

This paper concerns Honji instability generated around a circular cylinder in an oscillatory flow with a small oblique angle. In this study, direct numerical simulation has been conducted for an oscillatory flow past a stationary cylinder with small incidence angles (α) of 5° and 10° at KC number of 2 and β number of 200. The three-dimensional Navier-Stokes equations are solved using the Petrov-Galerkin finite element method. Flow structures around the cylinder are visualized through using streamlines, velocity vectors and vorticity contours. Honji instability has been captured at both chosen inclination angles. However Honji vortex pairs are asymmetric at α = 5° and 10° due to the inclination of the oscillation direction and can only be observed during the flow reversal. It is also found that the flow inclination appears to suppress the three-dimensional instability.


1995 ◽  
Vol 298 ◽  
pp. 211-248 ◽  
Author(s):  
U. Rist ◽  
H. Fasel

The three-dimensional development of controlled transition in a flat-plate boundary layer is investigated by direct numerical simulation (DNS) using the complete Navier-Stokes equations. The numerical investigations are based on the so-called spatial model, thus allowing realistic simulations of spatially developing transition phenomena as observed in laboratory experiments. For solving the Navier-Stokes equations, an efficient and accurate numerical method was developed employing fourth-order finite differences in the downstream and wall-normal directions and treating the spanwise direction pseudo-spectrally. The present paper focuses on direct simulations of the wind-tunnel experiments by Kachanov et al. (1984, 1985) of fundamental breakdown in controlled transition. The numerical results agreed very well with the experimental measurements up to the second spike stage, in spite of relatively coarse spanwise resolution. Detailed analysis of the numerical data allowed identification of the essential breakdown mechanisms. In particular, from our numerical data, we could identify the dominant shear layers and vortical structures that are associated with this breakdown process.


2001 ◽  
Vol 105 (1046) ◽  
pp. 193-198 ◽  
Author(s):  
N. D. Sandham ◽  
M. Alam ◽  
S. Morin

Abstract A method is proposed by which a direct numerical simulation of the compressible Navier-Stokes equations may be embedded within a more general aeronautical CFD code. The method may be applied to any code which solves the Euler equations or the Favre-averaged Navier-Stokes equations. A formal decomposition of the flowfield is used to derive modified equations for use with direct numerical simulation solvers. Some preliminary applications for model flows with transitional separation bubbles are given.


Sign in / Sign up

Export Citation Format

Share Document