scholarly journals RANS Numerical Simulation of Turbulent Particulate Pipe Flow for Fixed Reynolds Number

Author(s):  
Alexander Kartushinsky ◽  
Ylo Rudi ◽  
Igor Shcheglov ◽  
Sergei Tisler ◽  
Igor Krupenski
2019 ◽  
Vol 879 ◽  
pp. 579-603 ◽  
Author(s):  
Kohei Kunii ◽  
Takahiro Ishida ◽  
Yohann Duguet ◽  
Takahiro Tsukahara

Annular Couette flow is the flow between two coaxial cylinders driven by the axial translation of the inner cylinder. It is investigated using direct numerical simulation in long domains, with an emphasis on the laminar–turbulent coexistence regime found for marginally low values of the Reynolds number. Three distinct flow regimes are demonstrated as the radius ratio $\unicode[STIX]{x1D702}$ is decreased from 0.8 to 0.5 and finally to 0.1. The high-$\unicode[STIX]{x1D702}$ regime features helically shaped turbulent patches coexisting with laminar flow, as in planar shear flows. The moderate-$\unicode[STIX]{x1D702}$ regime does not feature any marked laminar–turbulent coexistence. In an effort to discard confinement effects, proper patterning is, however, recovered by artificially extending the azimuthal span beyond $2\unicode[STIX]{x03C0}$. Eventually, the low-$\unicode[STIX]{x1D702}$ regime features localised turbulent structures different from the puffs commonly encountered in transitional pipe flow. In this new coexistence regime, turbulent fluctuations are surprisingly short-ranged. Implications are discussed in terms of phase transition and critical scaling.


1999 ◽  
Vol 387 ◽  
pp. 39-60 ◽  
Author(s):  
H. SHAN ◽  
B. MA ◽  
Z. ZHANG ◽  
F. T. M. NIEUWSTADT

A direct numerical simulation of transitional pipe flow is carried out with the help of a spectral element method and used to investigate the localized regions of ‘turbulent’ flow that are observed in experiments. Two types of such regions can be distinguished: the puff and the slug. The puff, which is generally found at low values of the Reynolds numbers, is simulated for Re = 2200 where the Reynolds number Re is based on the mean velocity UB and pipe diameter D. The slug occurs at a higher Reynolds number and it is simulated for Re = 5000. The computations start with a laminar pipe flow to which is added a prescribed velocity disturbance at a given axial position and for a finite time. The disturbance then evolves further into a puff or slug structure.The simulations confirm the experimentally observed fact that for a puff the velocity near the leading edge changes more gradually than for a slug where an almost discontinuous change is observed. The positions of the leading and trailing edges of the puff and slug are computed from the simulations as a function of time. The propagation velocity of the leading edge is found to be constant and equal to 1.56UB and 1.69UB for the puff and slug, respectively. For the trailing edge the velocity is found to be 0.73UB and 0.52UB, respectively. By rescaling the simulation results obtained at various times to a fixed length, we define an ensemble average. This method is used to compute the average characteristics of the puff and slug such as the spatial distribution of the mean velocity, the turbulent velocity fluctuations and also the wall shear stress. By computing particle trajectories we have investigated the entrainment and detrainment of fluid by a puff and slug. We find that the puff detrains through its trailing edge and entrains through its leading edge. The slug entrains fluid through its leading and through most of its trailing edge. As a consequence the fluid inside the puff is constantly exchanged with fluid outside whereas the fluid inside a slug remains there. These entrainment/detrainment properties which are in agreement with the measurements of Wygnanski & Champagne (1973) imply that the puff has the characteristics of a wave phenomenon while the slug can be characterized more as a material property which travels with the flow.Finally, we have investigated in more detail the velocity field within the puff. In a coordinate system that travels with the mean velocity we find recirculation regions both near the trailing and leading edges which agrees at least qualitatively with experimental data. We also find streamwise vortices, predominantly in the trailing-edge region which have been also observed in experiments and which are believed to play an important role in the dynamics of the transition process.


Author(s):  
Alessandro Chiarini ◽  
Maurizio Quadrio

AbstractA direct numerical simulation (DNS) of the incompressible flow around a rectangular cylinder with chord-to-thickness ratio 5:1 (also known as the BARC benchmark) is presented. The work replicates the first DNS of this kind recently presented by Cimarelli et al. (J Wind Eng Ind Aerodyn 174:39–495, 2018), and intends to contribute to a solid numerical benchmark, albeit at a relatively low value of the Reynolds number. The study differentiates from previous work by using an in-house finite-differences solver instead of the finite-volumes toolbox OpenFOAM, and by employing finer spatial discretization and longer temporal average. The main features of the flow are described, and quantitative differences with the existing results are highlighted. The complete set of terms appearing in the budget equation for the components of the Reynolds stress tensor is provided for the first time. The different regions of the flow where production, redistribution and dissipation of each component take place are identified, and the anisotropic and inhomogeneous nature of the flow is discussed. Such information is valuable for the verification and fine-tuning of turbulence models in this complex separating and reattaching flow.


2008 ◽  
Vol 602 ◽  
pp. 209-218 ◽  
Author(s):  
J. J. J. GILLISSEN ◽  
B. J. BOERSMA ◽  
P. H. MORTENSEN ◽  
H. I. ANDERSSON

We use direct numerical simulation to study turbulent drag reduction by rigid polymer additives, referred to as fibres. The simulations agree with experimental data from the literature in terms of friction factor dependence on Reynolds number and fibre concentration. An expression for drag reduction is derived by adopting the concept of the elastic layer.


Sign in / Sign up

Export Citation Format

Share Document