Direct Numerical Simulation of Turbulent Rayleigh-Be´nard Convection Using a Higher Order Scheme

2007 ◽  
Author(s):  
A. K. De ◽  
V. Eswaran
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Stokes Equations ◽  
Parallel Computer ◽  
Navier Stokes ◽  
Cross Sectional ◽  
Navier Stokes Equations ◽  
Higher Order Scheme ◽  
Velocity Variances ◽  
Boussinesq Fluid

In the present paper a Direct Numerical Simulation of turbulent Rayleigh-Be´nard convection of a Boussinesq fluid with Prandtl number 0.7 in a 6 : 6 : 1 aspect ratio box at Ra = 2.106 is attempted. The incompressible Navier-Stokes equations are solved on non-uniform cartesian grids using a finite difference method with a high resolution convective scheme and 2nd order Adams-Bashforth Crank-Nicolson (ABCN) time-stepping. The numerical simulation is conducted on a parallel computer. We observed that temperature and velocity variances lead to different length and temperature scales. Negative skewnesses of vertical velocity and its derivative are observed. A non-uniform distribution of temperature near the surface proceeds to an exponential like variation in the center via an intermediate region outside the surface layer. Instantaneous snapshots of the flow field reveal intense sheet-like structures near the surface layers becoming plumes in the middle region. In the vertical cross-sectional planes, strong shearing effect near the surfaces, horizontal spreading and induction of plumes are observed.

On Direct Numerical Simulation of Turbulent Flows

2011 ◽  
Vol 64 (2) ◽  
Author(s):  
Giancarlo Alfonsi
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
High Performance Computing ◽  
Turbulent Flows ◽  
High Performance ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Turbulent Shear ◽  
Navier Stokes Equations ◽  
Performance Computing

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.

scholarly journals Numerical stability analysis of a vortex ring with swirl

2019 ◽  
Vol 878 ◽  
pp. 5-36 ◽  
Author(s):  
Yuji Hattori ◽  
Francisco J. Blanco-Rodríguez ◽  
Stéphane Le Dizès
Keyword(s):  
Numerical Simulation ◽  
Growth Rate ◽  
Direct Numerical Simulation ◽  
Vortex Ring ◽  
Stokes Equations ◽  
Base Flow ◽  
Critical Layer ◽  
Navier Stokes ◽  
Instability Mode ◽  
Navier Stokes Equations

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.

Direct Numerical Simulation of Effects of Small Angle of Incidence on Honji Instability

2011 ◽  
Author(s):  
Kun Yang ◽  
Liang Cheng ◽  
Hongwei An ◽  
Ming Zhao
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Oscillatory Flow ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Angle Of Incidence ◽  
Navier Stokes Equations ◽  
Small Incidence ◽  
Inclination Angles

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.

Direct numerical simulation of controlled transition in a flat-plate boundary layer

1995 ◽  
Vol 298 ◽  
pp. 211-248 ◽  
Author(s):  
U. Rist ◽  
H. Fasel
Keyword(s):  
Numerical Simulation ◽  
Boundary Layer ◽  
Direct Numerical Simulation ◽  
Flat Plate ◽  
Stokes Equations ◽  
Plate Boundary ◽  
Numerical Data ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Flat Plate Boundary Layer

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.

scholarly journals Embedded direct numerical simulation for aeronautical CFD

2001 ◽  
Vol 105 (1046) ◽  
pp. 193-198 ◽  
Author(s):  
N. D. Sandham ◽  
M. Alam ◽  
S. Morin
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Euler Equations ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Modified Equations ◽  
Separation Bubbles

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.

SMOOTHED PROFILE METHOD TO SIMULATE COLLOIDAL PARTICLES IN COMPLEX FLUIDS

2009 ◽  
Vol 20 (09) ◽  
pp. 1457-1465 ◽  
Author(s):  
RYOICHI YAMAMOTO ◽  
YASUYA NAKAYAMA ◽  
KANG KIM
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Colloidal Particles ◽  
Stokes Equations ◽  
Electrolyte Solutions ◽  
Rigid Body Dynamics ◽  
Spherical Particles ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Charged Colloids

A new direct numerical simulation scheme, called "Smoothed Profile (SP) method," is presented. The SP method, as a direct numerical simulation of particulate flow, provides a way to couple continuum fluid dynamics with rigid-body dynamics through smoothed profile of colloidal particle. Our formulation includes extensions to colloids in multicomponent solvents such as charged colloids in electrolyte solutions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier–Stokes equations so that the electro-hydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters.

Sign in / Sign up

Export Citation Format

Share Document