Study of turbulent flow past a square cylinder using partially-averaged Navier–Stokes method in OpenFOAM

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.


2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.


2012 ◽  
Vol 16 (5) ◽  
pp. 1510-1514
Author(s):  
Tian Li ◽  
Li-Hao Zhao ◽  
Xiao-Ke Ku ◽  
Helge Andersson ◽  
Terese Lovas

This paper investigates the performance of Reynolds-averaged Navier-Stokes model on dispersion of particles in wall turbulence. A direct numerical simulation of wall-bounded channel flow with particles suspensions was set as a benchmark. The standard k-? model coupled with two different eddy interaction models was used in Reynolds-averaged Navier-Stokes model and compared to the direct numerical simulation. Detailed comparisons between direct numerical simulation and Reynolds-averaged Navier-Stokes model on particle distribution evolving over time were carried out.


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.


2005 ◽  
Vol 73 (3) ◽  
pp. 413-421 ◽  
Author(s):  
Sharath S. Girimaji

A turbulence bridging method purported for any filter-width or scale resolution—fully averaged to completely resolved—is developed. The method is given the name partially averaged Navier-Stokes (PANS) method. In PANS, the model filter width (extent of partial averaging) is controlled through two parameters: the unresolved-to-total ratios of kinetic energy (fk) and dissipation (fε). The PANS closure model is derived formally from the Reynolds-averaged Navier-Stokes (RANS) model equations by addressing the following question: if RANS represents the closure for fully averaged statistics, what is the corresponding closure for partially averaged statistics? The PANS equations vary smoothly from RANS equations to Navier-Stokes (direct numerical simulation) equations, depending on the values of the filter-width control parameters. Preliminary results are very encouraging.


2011 ◽  
Vol 2011 ◽  
pp. 1-13 ◽  
Author(s):  
Nilanjan Chakraborty ◽  
Andrei N. Lipatnikov

The statistics of mean fluid velocity components conditional in unburned reactants and fully burned products in the context of Reynolds Averaged Navier Stokes (RANS) simulations have been studied using a Direct Numerical Simulation database of statistically planar turbulent premixed flame representing the corrugated flamelets regime combustion. Expressions for conditional mean velocity and conditional velocity correlations which are derived based on a presumed bimodal probability density function of reaction progress variable for unity Lewis number flames are assessed in this study with respect to the corresponding quantities extracted from DNS data. In particular, conditional surface averaged velocities(ui)¯Rsand the velocity correlations(uiu)j¯Rsin the unburned reactants are demonstrated to be effectively modelled by the unconditional velocities(ui)¯Rand velocity correlations(uiuj)¯R, respectively, for the major part of turbulent flame brush with the exception of the leading edge. By contrast, conditional surface averaged velocities(ui)¯Psand the velocity correlations(uiu)j¯Psin fully burned products are shown to be markedly different from the unconditional velocities(ui)¯Pand velocity correlations(uiuj)¯P, respectively.


Sign in / Sign up

Export Citation Format

Share Document