A direct numerical simulation study of interface propagation in homogeneous turbulence

2015 ◽  
Vol 772 ◽  
pp. 127-164 ◽  
Author(s):  
R. Yu ◽  
X.-S. Bai ◽  
A. N. Lipatnikov

A 3D direct numerical simulation (DNS) study of the evolution of a self-propagating interface in forced constant-density statistically stationary homogeneous isotropic turbulence was performed by solving Navier–Stokes and level-set equations under a wide range of conditions that cover various (from 0.1 to 2.0) ratios of the interface speed $S_{L}$ to the r.m.s. turbulent velocity  $U^{\prime }$ and various (50, 100 and 200) turbulent Reynolds numbers $\mathit{Re}$. By analysing computed data, the following issues were addressed: (i) dependence of the speed and thickness of the fully developed statistically planar mean front that envelops the interface on $U^{\prime }/S_{L}$ and $\mathit{Re}$, (ii) dependence of the fully developed mean turbulent flux of a scalar $c$ that characterizes the state of the fluid ($c=0$ and 1 ahead and behind the interface respectively) on $U^{\prime }/S_{L}$ and $\mathit{Re}$, (iii) evolution of the mean front speed, its thickness, and the mean scalar flux during the front development after embedding a planar interface into the forced turbulence and (iv) relation between canonical and conditioned moments of the velocity, velocity gradient and pressure gradient fields.

Author(s):  
Arnab Chakraborty ◽  
HV Warrior

The present paper reports numerical simulation of turbulent flow over a square cylinder using a novel scale resolving computational fluid dynamics technique named Partially-Averaged Navier–Stokes (PANS), which bridges Reynolds-Averaged Navier–Stokes (RANS) with Direct Numerical Simulation (DNS) in a seamless manner. All stream-wise and wall normal mean velocity components, turbulent stresses behavior have been computed along the flow (streamwise) as well as in transverse (wall normal) direction. The measurement locations are chosen based on the previous studies so that results could be compared. However, the Reynolds number ( Re) of the flow is maintained at 21,400 and K– ω turbulence model is considered for the present case. All the computations are performed in OpenFOAM framework using a finite volume solver. Additionally, turbulent kinetic energy variations are presented over a wide range of measurement planes in order to explain the energy transfer process in highly unsteady turbulent flow field. The fluctuating root mean square velocities in the streamwise as well as in the wall normal direction have been discussed in the present work. It has been found that Partially-Averaged Navier–Stokes (PANS) model is capable of capturing the properties of highly unsteady turbulent flows and gives better results than Reynolds-Averaged Navier–Stokes (RANS). The results obtained using Partially-Averaged Navier–Stokes (PANS) are quite comparable with Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) data available in literature. The partially-averaged Navier–Stokes results are compared with our simulated Reynolds-Averaged Navier–Stokes (RANS) results, available experimental as well as numerical results in literature and it is found to be good in agreement.


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.


2001 ◽  
Author(s):  
X. Ai ◽  
B. Q. Li

Abstract Turbulent magnetically flows occur in a wide range of material processing systems involving electrically conducting melts. This paper presents a parallel higher order scheme for the direct numerical simulation of turbulent magnetically driven flows in induction channels. The numerical method is based on the higher order finite difference algorithm, which enjoys the spectral accuracy while minimizing the computational intensity. This, coupled with the parallel computing strategy, provides a very useful means to simulate turbulent flows. The higher order finite difference formulation of magnetically driven flow problems is described in this paper. The details of the parallel algorithm and its implementation for the simulations on parallel machines are discussed. The accuracy and numerical performance of the higher order finite difference scheme are assessed in comparison with the spectral method. The examples of turbulent magnetically driven flows in induction channels and pressure gradient driven flows in regular channels are given, and the computed results are compared with experimental measurements wherever possible.


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):  
Zhenzhong Li ◽  
Jinjia Wei ◽  
Bo Yu

Multiphase flow with particles covers a wide spectrum of flow conditions in natural world and industrial applications. The experiments and the direct numerical simulation have become the most popular means to study the dilute particle-laden flow in the last two decades. In the experimental study, the mean Reynolds number is often adjusted to the value of single-phase flow for each set of particle conditions. However, the friction Reynolds number usually keeps invariable in the direct numerical simulation of the particle-laden flows for convenience. In this study the effect of the difference between given mean Reynolds number and friction Reynolds number was investigated. Two simulations were performed for each set of particle parameters, and the mean Reynolds number and friction Reynolds number were kept invariant respectively. From the results it can be found that the turbulence intensity and the dimensionless velocities are larger when keeping the friction Reynolds constant. And the results calculated from the cases of keeping the mean Reynolds number invariable agree with the experiment results better. In addition, the particle distribution along the wall-normal coordinate was found to be unchanged between two simulation conditions. As a suggestion, keeping the same mean Reynolds number in the direct numerical simulation of particle-laden flow is more appropriate.


2018 ◽  
Vol 74 (6) ◽  
pp. 840-860 ◽  
Author(s):  
Sagar Dave ◽  
Chetankumar Anghan ◽  
Shaswat Saincher ◽  
Jyotirmay Banerjee

2013 ◽  
Vol 5 (3) ◽  
pp. 435-445
Author(s):  
M. S. I. Mallik ◽  
M. A. Uddin ◽  
M. A. Rahman

Direct numerical simulation (DNS) in two-dimensional homogeneous isotropic turbulence is performed by using the Spectral method at a Reynolds number Re = 1000 on a uniformly distributed grid points. The Reynolds number is low enough that the computational grid is capable of resolving all the possible turbulent scales. The statistical properties in the computed flow field show a good agreement with the qualitative behavior of decaying turbulence. The behavior of the flow structures in the computed flow field also follow the classical idea of the fluid flow in turbulence. Keywords: Direct numerical simulation, Isotropic turbulence, Spectral method. © 2013 JSR Publications. ISSN: 2070-0237 (Print); 2070-0245 (Online). All rights reserved. doi:http://dx.doi.org/10.3329/jsr.v5i3.12665 J. Sci. Res. 5 (3), 435-445 (2013)  


Sign in / Sign up

Export Citation Format

Share Document