A STUDY ON CORRELATION MOMENTS OF TWO-PHASE FLUCTUATING VELOCITY USING DIRECT NUMERICAL SIMULATION

2013 ◽  
Vol 24 (10) ◽  
pp. 1350068 ◽  
Author(s):  
BING WANG ◽  
WEI WEI ◽  
HUIQIANG ZHANG

Existing models of two-phase fluctuating velocity correlation moments are unsatisfactory because of their inability to clearly identify the dependency of two-phase velocity covariance on fluid- and particle-phase velocity second moments. This is especially true of wall-bounded turbulent flows. In this paper, the statistical fluctuating velocity of both phases in particle-laden turbulent channel flows were obtained numerically by means of direct numerical simulation (DNS) coupled to the Lagrangian particle trajectory method. The effects of particle Stokes number on the scaling of two-phase fluctuating velocity correlation moments were analyzed considering effects of flow inhomogeneity. An improved two-phase correlation closure model of exponential decay with emphasis on the particle-phase kinetic energy was then proposed based on the results of an evaluation of five existing models. This new model was found to be better than previous models, which used local equilibrium assumption. The present investigations may facilitate understanding of two-phase flow physics and the construction of models capable of predicting the movements of particle-laden turbulent flows accurately using Reynolds-averaged Navier–Stokes (RANS) methods.

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.


2013 ◽  
Vol 13 (04) ◽  
pp. 1350066 ◽  
Author(s):  
O. ANWAR BÉG ◽  
M. M. RASHIDI ◽  
N. RAHIMZADEH ◽  
TASVEER A. BÉG ◽  
TIN-KAN HUNG

A two-phase thermo-hydrodynamic model is presented for transport in the vertical chamber of a porous media blood filtration device. A non-Darcy drag force formulation was employed. The Marble–Drew fluid–particle suspension model was used to simulate the plasma phase and the suspension (erythrocyte) particle phase. The non-dimensional transport equations were solved using a semi-computational procedure known as the homotopy analysis method (HAM). With the judicious use of the auxiliary parameter ℏ, HAM affords a powerful mechanism to adjust and control the convergence region of solution series. This method provides an efficient approximate analytical solution with high accuracy, minimal calculation and avoidance of physically unrealistic assumptions. Detailed computations are presented for the effects of Grashof number (Gr), momentum inverse Stokes number (Skm), Darcy number (Da), Forchheimer number (Fs), particle loading parameter (PL), buoyancy parameter (B) and temperature inverse Stokes number (SkT) on the dimensionless fluid phase velocity (U), dimensionless particle phase velocity (Up), dimensionless fluid phase temperature (Φ) and the dimensionless temperature of particle phase (Φp). A Prandtl number of 25 was used to simulate blood at room temperature. Excellent correlation was obtained between the HAM and numerical shooting quadrature solutions. The results indicated that there is a strong decrease in fluid phase velocities with increasing Darcian (first order) drag and second-order Forchheimer drag, and a weaker reduction in particle phase velocity field. Applications of the study include porous media bio-filtration devices and dialysis simulations.


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.


Author(s):  
Mosayeb Shams ◽  
Kamaljit Singh ◽  
Branko Bijeljic ◽  
Martin J. Blunt

AbstractThis study focuses on direct numerical simulation of imbibition, displacement of the non-wetting phase by the wetting phase, through water-wet carbonate rocks. We simulate multiphase flow in a limestone and compare our results with high-resolution synchrotron X-ray images of displacement previously published in the literature by Singh et al. (Sci Rep 7:5192, 2017). We use the results to interpret the observed displacement events that cannot be described using conventional metrics such as pore-to-throat aspect ratio. We show that the complex geometry of porous media can dictate a curvature balance that prevents snap-off from happening in spite of favourable large aspect ratios. We also show that pinned fluid-fluid-solid contact lines can lead to snap-off of small ganglia on pore walls; we propose that this pinning is caused by sub-resolution roughness on scales of less than a micron. Our numerical results show that even in water-wet porous media, we need to allow pinned contacts in place to reproduce experimental results.


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.


Sign in / Sign up

Export Citation Format

Share Document