Effects of roughness on particle dynamics in turbulent channel flows: a DNS analysis

2014 ◽  
Vol 739 ◽  
pp. 465-478 ◽  
Author(s):  
Barbara Milici ◽  
Mauro De Marchis ◽  
Gaetano Sardina ◽  
Enrico Napoli

AbstractDeposition and resuspension mechanisms in particle-laden turbulent flows are dominated by the coherent structures arising in the wall region. These turbulent structures, which control the turbulent regeneration cycles, are affected by the roughness of the wall. The particle-laden turbulent flow in a channel bounded by irregular two-dimensional rough surfaces is analysed. The behaviour of dilute dispersions of heavy particles is analysed using direct numerical simulations (DNS) to calculate the three-dimensional turbulent flow and Lagrangian tracking to describe the turbophoretic effect associated with two-phase turbulent flows in a complex wall-bounded domain. Turbophoresis is investigated in a quantitative way as a function of the particle inertia. The analysis of the particle statistics, in term of mean particle concentration and probability density function (p.d.f.) of wall-normal particle velocity, shows that the wall roughness produces a completely different scenario compared to the classical smooth wall. The effect of the wall roughness on the particle mass flux is shown for six particle populations having different inertia.

Author(s):  
Kadeem Dennis ◽  
Kamran Siddiqui

The boundary layers are known for their significance in several engineering systems. In particular, the inner region of the turbulent boundary layer has been shown to play a significant role in controlling the dynamics of turbulent structures that are responsible for the transport of mass, heat and momentum. While substantial work has been done in the past to characterize the structure of turbulent flow in this region, the characterization of the three-dimensional turbulent flow structure is limited. This study reports a multi-plane particle image velocimetry (PIV) approach to investigate three-dimensional dynamics of the turbulent boundary layer in the near-wall region. Planar PIV is used to capture two-dimensional fluid velocity fields in several planes with respect to the fluid flow direction. These results are used to describe three-dimensional turbulent events given by key quantities such as mean and turbulent velocities and turbulent kinetic energy.


Author(s):  
Veeraraghava R Hasti ◽  
Prithwish Kundu ◽  
Sibendu Som ◽  
Jay P Gore

The turbulent flow field in a practical gas turbine combustor is very complex because of the interactions between various flows resulting from components like multiple types of swirlers, dilution holes, and liner effusion cooling holes. Numerical simulations of flows in such complex combustor configurations are challenging. The challenges result from (a) the complexities of the interfaces between multiple three-dimensional shear layers, (b) the need for proper treatment of a large number of tiny effusion holes with multiple angles, and (c) the requirements for fast turnaround times in support of engineering design optimization. Both the Reynolds averaged Navier–Stokes simulation (RANS) and the large eddy simulation (LES) for the practical combustor geometry are considered. An autonomous meshing using the cut-cell Cartesian method and adaptive mesh refinement (AMR) is demonstrated for the first time to simulate the flow in a practical combustor geometry. The numerical studies include a set of computations of flows under a prescribed pressure drop across the passage of interest and another set of computations with all passages open with a specified total flow rate at the plenum inlet and the pressure at the exit. For both sets, the results of the RANS and the LES flow computations agree with each other and with the corresponding measurements. The results from the high-resolution LES simulations are utilized to gain fundamental insights into the complex turbulent flow field by examining the profiles of the velocity, the vorticity, and the turbulent kinetic energy. The dynamics of the turbulent structures are well captured in the results of the LES simulations.


2021 ◽  
Author(s):  
Christina Tsai ◽  
Kuang-Ting Wu

<p>It is demonstrated that turbulent boundary layers are populated by a hierarchy of recurrent structures, normally referred to as the coherent structures. Thus, it is desirable to gain a better understanding of the spatial-temporal characteristics of coherent structures and their impact on fluid particles. Furthermore, the ejection and sweep events play an important role in turbulent statistics. Therefore, this study focuses on the characterizations of flow particles under the influence of the above-mentioned two structures.</p><div><span>With regard to the geometry of turbulent structures, </span><span>Meinhart & Adrian (1995) </span>first highlighted the existence of large and irregularly shaped regions of uniform streamwise momentum zone (hereafter referred to as a uniform momentum zone, or UMZs), regions of relatively similar streamwise velocity with coherence in the streamwise and wall-normal directions.  <span>Subsequently, </span><span>de Silva et al. (2017) </span><span>provided a detection criterion that had previously been utilized to locate the uniform momentum zones (UMZ) and demonstrated the application of this criterion to estimate the spatial locations of the edges that demarcates UMZs.</span></div><div> </div><div>In this study, detection of the existence of UMZs is a pre-process of identifying the coherent structures. After the edges of UMZs are determined, the identification procedure of ejection and sweep events from turbulent flow DNS data should be defined. As such, an integrated criterion of distinguishing ejection and sweep events is proposed. Based on the integrated criterion, the statistical characterizations of coherent structures from available turbulent flow data such as event durations, event maximum heights, and wall-normal and streamwise lengths can be presented.</div>


1991 ◽  
Vol 113 (1) ◽  
pp. 34-41 ◽  
Author(s):  
G. J. Yoo ◽  
R. M. C. So ◽  
B. C. Hwang

Internal rotating boundary-layer flows are strongly influenced by large circumferential strain and the turbulence field is anisotropic. This is especially true in the entry region of a rotating pipe where the flow is three dimensional, the centrifugal force due to fluid rotation is less important, and the circumferential strain created by surface rotation has a significant effect on the turbulence field near the wall. Consequently, viscous effects cannot be neglected in the near-wall region. Several low-Reynolds-number turbulence closures are proposed for the calculation of developing rotating pipe flows. Some are two-equation closures with and without algebraic stress correction, while others are full Reynolds-stress closures. It is found that two-equation closures with and without algebraic stress correction are totally inadequate for this three-dimensional flow, while Reynolds-stress closures give results that are in good agreement with measurements over a wide range of rotation numbers.


Author(s):  
B. Song ◽  
R. S. Amano ◽  
S. Sitarama ◽  
B. Lin

Numerical study on a three-dimensional turbulent flow in a turbine-rotor passage is presented in this paper. The standard k-ε model was used for the first phase of the turbulence computations. The computations were further extended by employing the full Reynolds-stress closure model (RSM). The computational results obtained using these models were compared in order to investigate the turbulence effect in the near-wall region. The governing equations in a generalized curvilinear coordinate system are discretized by using the SIMPLEC method with non-staggered grids. The oscillations in pressure and velocity due to non-staggered grids are eliminated by using a special interpolation method. The predicted midspan pressure coefficients using the k-ε model and the RSM are compared with the experimental data. It was shown that the present results obtained by using either model are fairly reasonable. Computations were then extended to cover the entire blade-to-blade flow passage, and the three-dimensional effects on pressure and turbulence kinetic energy were evaluated. It was observed that the two turbulence models predict different results for the turbulence kinetic energy. This variation was identified as being related to some non-isotropic turbulence occurring near the blade surface due to the severe acceleration of the flow. It was thus proven that the models based on the RSM give more realistic predictions for highly turbulent cascade flow computations than a Boussinesq viscosity model.


Author(s):  
Ling Zhen ◽  
Claudia del Carmen Gutierrez-Torres

The question of “where and how the turbulent drag arises” is one of the most fundamental problems unsolved in fluid mechanics. However, the physical mechanism responsible for the friction drag reduction is still not well understood. Over decades, it is found that the turbulence production and self-containment in a boundary layer are organized phenomena and not random processes as the turbulence looks like. The further study in the boundary layer should be able to help us know more about the mechanisms of drag reduction. The wavelet-based vector multi-resolution technique was proposed and applied to the two dimensional PIV velocities for identifying the multi-scale turbulent structures. The intermediate and small scale vortices embedded within the large-scale vortices were separated and visualized. By analyzing the fluctuating velocities at different scales, coherent eddy structures were obtained and this help us obtain the important information on the multi-scale flow structures in the turbulent flow. By comparing the eddy structures in different operating conditions, the mechanism to explain the drag reduction caused by micro bubbles in turbulent flow was proposed.


2011 ◽  
Vol 2011 ◽  
pp. 1-20 ◽  
Author(s):  
Wang Wenquan ◽  
Zhang Lixiang ◽  
Yan Yan ◽  
Guo Yakun

An innovative computational model is presented for the large eddy simulation (LES) of multidimensional unsteady turbulent flow problems in complex geometries. The main objectives of this research are to know more about the structure of turbulent flows, to identify their three-dimensional characteristic, and to study physical effects due to complex fluid flow. The filtered Navier-Stokes equations are used to simulate large scales; however, they are supplemented by dynamic subgrid-scale (DSGS) models to simulate the energy transfer from large scales toward subgrid-scales, where this energy will be dissipated by molecular viscosity. Based on the Taylor-Galerkin schemes for the convection-diffusion problems, this model is implemented in a three-dimensional finite element code using a three-step finite element method (FEM). Turbulent channel flow and flow over a backward-facing step are considered as a benchmark for validating the methodology by comparing with the direct numerical simulation (DNS) results or experimental data. Also, qualitative and quantitative aspects of three-dimensional complex turbulent flow in a strong 3D blade passage of a Francis turbine are analyzed.


2014 ◽  
Vol 758 ◽  
Author(s):  
Matteo Bernardini

AbstractThe effect of the Reynolds number on the behaviour of inertial particles in wall-bounded turbulent flows is investigated through large-scale direct numerical simulations (DNS) of particle-laden canonical channel flow spanning almost a decade in the friction Reynolds number, from $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\mathit{Re}_{\tau } = 150$ to $\mathit{Re}_{\tau } = 1000$. Lagrangian particle tracking is used to study the motion of six different particle sets, described by a Stokes number in the range $\mathit{St} = 1\text {--}1000$. At all Reynolds numbers a strong segregation in the near-wall region is observed for particles characterized by intermediate Stokes number, in the range $\mathit{St} =10\text {--}100$. The wall-normal concentration profiles of such particles collapse in inner scaling, thus suggesting the independence of the turbophoretic drift from the large-scale outer motions. This observation is also supported by the spatial organization of the suspended phase in the inner layer, which is found to be universal with the Reynolds number. The deposition rate coefficient increases with $\mathit{Re}_{\tau }$ for a given $\mathit{St}$. Suitable inner and outer scalings are proposed to collapse the deposition curves across the available ranges of Reynolds and Stokes numbers for the different deposition regimes.


2016 ◽  
Author(s):  
Joris T. Eggenhuisen ◽  
Matthieu J. B. Cartigny ◽  
Jan de Leeuw

Abstract. The inability to capture the physics of solid-particle suspension in turbulent fluid flow is holding back application of multiphase computational fluid dynamics techniques to the many problems involving particle suspension in nature and society. We present a theory for particle suspension capacity near no-slip frictional boundaries of turbulent flows. The suspension capacity parameter Γ includes universal turbulent flow scales and material properties of the fluid and particles only. Comparison to measurements shows that Γ = 1 gives the upper limit of observed suspended particle concentrations in a broad range of flume experiments and field settings. The condition of Γ > 1 coincides with complete suppression of coherent turbulent structures near the boundary in Direct Numerical Simulations of sediment-laden turbulent flow. The theory outperforms previous empiric relations when compared to data. It can be applied as a concentration boundary condition in modelling studies of dispersion of particulates in environmental and man-made flows.


Author(s):  
Genglin Tang ◽  
Roger L. Simpson ◽  
Qing Tian

Experimental results are presented from a study of the tip-gap turbulent flow structure in a low-speed linear compressor cascade wind tunnel at Virginia Tech by utilizing surface oil flow visualization, endwall pressure measurements, and instantaneous velocity measurements with a custom-made 3-orthogonal-velocity-component fiber-optic laser-Doppler velocimetry (LDV) system. Tip gap flows are pressure-driven and highly skewed three-dimensional turbulent flows. The crossflow velocity normal to the blade chord is nearly uniform in the mid tip gap and changes substantially from the pressure to suction side due to the local tip pressure loading while the TKE does not vary much across the mid tip gap. The tip gap flow correlations of streamwise and wall normal velocity fluctuations decrease significantly from the leading edge to the trailing edge of the blade due to flow skewing.


Sign in / Sign up

Export Citation Format

Share Document