A Parallel HOFD Scheme for Direct Numerical Simulation of Turbulent Magnetically Driven Flows

2001 ◽  
Author(s):  
X. Ai ◽  
B. Q. Li
Keyword(s):  
Numerical Simulation ◽  
Finite Difference ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Parallel Machines ◽  
Higher Order ◽  
Electrically Conducting ◽  
Wide Range ◽  
Order Scheme ◽  
Higher Order Scheme

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.

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

2020 ◽  
Vol 234 (14) ◽  
pp. 2821-2832 ◽  
Author(s):  
Arnab Chakraborty ◽  
HV Warrior
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Turbulent Flow ◽  
Turbulent Flows ◽  
Normal Direction ◽  
Navier Stokes ◽  
Square Cylinder ◽  
Mean Velocity ◽  
Wide Range ◽  
Reynolds Averaged Navier Stokes

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.

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 Flow topologies in bubble-induced turbulence: a direct numerical simulation analysis

2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty
Keyword(s):  
Numerical Simulation ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Volume Fraction ◽  
Fluid Velocity ◽  
Small Scale ◽  
Flow Topology ◽  
Two Phase ◽  
Local Flow ◽  
Interface Curvature

This paper presents a detailed investigation of flow topologies in bubble-induced two-phase turbulence. Two freely moving and deforming air bubbles that have been suspended in liquid water under counterflow conditions have been considered for this analysis. The direct numerical simulation data considered here are based on the one-fluid formulation of the two-phase flow governing equations. To study the development of coherent structures, a local flow topology analysis is performed. Using the invariants of the velocity gradient tensor, all possible small-scale flow structures can be categorized into two nodal and two focal topologies for incompressible turbulent flows. The volume fraction of focal topologies in the gaseous phase is consistently higher than in the surrounding liquid phase. This observation has been argued to be linked to a strong vorticity production at the regions of simultaneous high fluid velocity and high interface curvature. Depending on the regime (steady/laminar or unsteady/turbulent), additional effects related to the density and viscosity jump at the interface influence the behaviour. The analysis also points to a specific term of the vorticity transport equation as being responsible for the induction of vortical motion at the interface. Besides the known mechanisms, this term, related to surface tension and gradients of interface curvature, represents another potential source of turbulence production that lends itself to further investigation.

Direct numerical simulation of turbulent channel flow up to

2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser
Keyword(s):  
Numerical Simulation ◽  
Reynolds Number ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Streamwise Velocity ◽  
Mean Velocity ◽  
Layer Structures ◽  
High Reynolds

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.

Sign in / Sign up

Export Citation Format

Share Document