Direct Numerical Simulation of Particle-Laden Turbulent Flows in Two-Way Coupling Using Equilibrium Assumption

2006 ◽  
Author(s):  
Babak Shotorban ◽  
S. Balachandar

A two-fluid approach is proposed for direct numerical simulation of particle-laden turbulent flows in two-way coupling through equilibrium assumption that is valid for particles with sufficiently small time constants. Making this assumption, a Eulerian velocity field is calculated for the particle phase through a truncated series expansion in terms of the velocity and acceleration of the fluid phase and the gravity acceleration. The transport equation of the Eulerian concentration field of particles (particle volume fraction) is solved along with the fluid phase equations for which the effect of particles on the fluid phase is taken into account through source terms in the momentum equations. For the assessment purposes, particle-laden isotropic turbulence is simulated. The results obtained through this approach are compared against those obtained by a trajectory approach in which the particle equations are solved in the Lagrangian framework. It is shown that there is a good agreement between the results obtained by the proposed two-fluid model and the trajectory approach by comparing the mean turbulent kinetic energy and its dissipation rate of the fluid phase as well as their spectra.

2018 ◽  
Vol 857 ◽  
pp. 270-290 ◽  
Author(s):  
Josef Hasslberger ◽  
Markus Klein ◽  
Nilanjan Chakraborty

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.


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.


2018 ◽  
Vol 19 (4) ◽  
pp. 401 ◽  
Author(s):  
Ahmed Zeeshan ◽  
Nouman Ijaz ◽  
Muhammad Mubashir Bhatti

This article addresses the influence of particulate-fluid suspension on asymmetric peristaltic motion through a curved configuration with mass and heat transfer. A motivation for the current study is that such kind of theory is helpful to examine the two-phase peristaltic motion between small muscles during the propagation of different biological fluids. Moreover, it is also essential in multiple applications of pumping fluid-solid mixtures by peristalsis, i.e., Chyme in small intestine and suspension of blood in arteriole. Long wavelength, as well as small Reynolds number, have been utilized to render the governing equations for particle and fluid phase. Exact solutions are presented for velocity (uf,p), temperature (θf,p) and concentration distributions (φf,p). All the parameters such as Prandtl number (Pr), particle volume fraction (C), suspension parameter (M1), curvature parameter (k), volumetric flow rate (Q), Schmidt number (Sc), phase difference (φ), Eckert number (Ec), and Soret number (Sr) discussed graphically for peristaltic pumping (Δp), pressure gradient (dp/dx), velocity (uf,p), temperature (θf,p) and concentration distributions (φf,p). The streamlines are also plotted with the aid of contour.


Sign in / Sign up

Export Citation Format

Share Document