Numerical Simulations of Particle-Laden Flow Based on Given Friction Reynolds Number and Mean Reynolds Number Respectively

Author(s):  
Zhenzhong Li ◽  
Jinjia Wei ◽  
Bo Yu

Multiphase flow with particles covers a wide spectrum of flow conditions in natural world and industrial applications. The experiments and the direct numerical simulation have become the most popular means to study the dilute particle-laden flow in the last two decades. In the experimental study, the mean Reynolds number is often adjusted to the value of single-phase flow for each set of particle conditions. However, the friction Reynolds number usually keeps invariable in the direct numerical simulation of the particle-laden flows for convenience. In this study the effect of the difference between given mean Reynolds number and friction Reynolds number was investigated. Two simulations were performed for each set of particle parameters, and the mean Reynolds number and friction Reynolds number were kept invariant respectively. From the results it can be found that the turbulence intensity and the dimensionless velocities are larger when keeping the friction Reynolds constant. And the results calculated from the cases of keeping the mean Reynolds number invariable agree with the experiment results better. In addition, the particle distribution along the wall-normal coordinate was found to be unchanged between two simulation conditions. As a suggestion, keeping the same mean Reynolds number in the direct numerical simulation of particle-laden flow is more appropriate.

2014 ◽  
Vol 136 (4) ◽  
Author(s):  
Zhi-Gang Feng

The laminar forced convection of a heated rotating sphere in air has been studied using a three-dimensional immersed boundary based direct numerical simulation method. A regular Eulerian grid is used to solve the modified momentum and energy equations for the entire flow region simultaneously. In the region that is occupied by the rotating sphere, a moving Lagrangian grid is used, which tracks the rotational motion of the particle. A force density function or an energy density function is introduced to represent the momentum interaction or thermal interaction between the sphere and fluid. This numerical method is validated by comparing simulation results with analytical solutions of heat diffusion problem and other published experimental data. The flow structures and the mean Nusselt numbers for flow Reynolds number ranging from 0 to 1000 are obtained. We compared our simulation results of the mean Nusselt numbers with the correlations from the literature and found a good agreement for flow Reynolds number greater than 500; however, a significant discrepancy arises at flow Reynolds number below 500. This leads us to develop a new equation that correlates the mean Nusselt number of a heated rotating sphere for flows of 0≤Re≤500.


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.


2017 ◽  
Vol 829 ◽  
pp. 392-419 ◽  
Author(s):  
V. Kitsios ◽  
A. Sekimoto ◽  
C. Atkinson ◽  
J. A. Sillero ◽  
G. Borrell ◽  
...  

The statistical properties are presented for the direct numerical simulation of a self-similar adverse pressure gradient (APG) turbulent boundary layer (TBL) at the verge of separation. The APG TBL has a momentum thickness-based Reynolds number range from $Re_{\unicode[STIX]{x1D6FF}_{2}}=570$ to 13 800, with a self-similar region from $Re_{\unicode[STIX]{x1D6FF}_{2}}=10\,000$ to 12 300. Within this domain the average non-dimensional pressure gradient parameter $\unicode[STIX]{x1D6FD}=39$, where for a unit density $\unicode[STIX]{x1D6FD}=\unicode[STIX]{x1D6FF}_{1}P_{\!e}^{\prime }/\unicode[STIX]{x1D70F}_{w}$, with $\unicode[STIX]{x1D6FF}_{1}$ the displacement thickness, $\unicode[STIX]{x1D70F}_{w}$ the mean shear stress at the wall and $P_{\!e}^{\prime }$ the far-field pressure gradient. This flow is compared with previous zero pressure gradient and mild APG TBL ($\unicode[STIX]{x1D6FD}=1$) results of similar Reynolds number. All flows are generated via the direct numerical simulation of a TBL on a flat surface with far-field boundary conditions tailored to apply the desired pressure gradient. The conditions for self-similarity, and the appropriate length and velocity scales, are derived. The mean and Reynolds stress profiles are shown to collapse when non-dimensionalised on the basis of these length and velocity scales. As the pressure gradient increases, the extent of the wake region in the mean streamwise velocity profiles increases, whilst the extent of the log-layer and viscous sublayer decreases. The Reynolds stress, production and dissipation profiles of the APG TBL cases exhibit a second outer peak, which becomes more pronounced and more spatially localised with increasing pressure gradient. This outer peak is located at the point of inflection of the mean velocity profiles, and is suggestive of the presence of a shear flow instability. The maximum streamwise velocity variance is located at a wall normal position of $\unicode[STIX]{x1D6FF}_{1}$ of spanwise wavelength of $2\unicode[STIX]{x1D6FF}_{1}$. In summary as the pressure gradient increases the flow has properties less like a zero pressure gradient TBL and more akin to a free shear layer.


1983 ◽  
Vol 48 (8) ◽  
pp. 2232-2248 ◽  
Author(s):  
Ivo Roušar ◽  
Michal Provazník ◽  
Pavel Stuhl

In electrolysers with recirculation, where a gas is evolved, the pumping of electrolyte from a lower to a higher level can be effected by natural convection due to the difference between the densities of the inlet electrolyte and the gaseous emulsion at the outlet. An accurate balance equation for calculation of the rate of flow of the pumped liquid is derived. An equation for the calculation of the mean volume fraction of bubbles in the space between the electrodes is proposed and verified experimentally on a pilot electrolyser. Two examples of industrial applications are presented.


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.


2021 ◽  
Vol 1877 (1) ◽  
pp. 012035
Author(s):  
Shengxiang Lin ◽  
Huanxiong Xia ◽  
Zhenyu Zhang ◽  
Jianhua Liu ◽  
Honglei Wang

2012 ◽  
Vol 2012 ◽  
pp. 1-13 ◽  
Author(s):  
Mohit Katragadda ◽  
Nilanjan Chakraborty ◽  
R. S. Cant

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Rethas been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Retand Ka before reaching an asymptotic value for large values of Retand Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel’dovich flame thicknessηi/δzdoes not exhibit any significant dependence on Retfor the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width.


2016 ◽  
Vol 30 (8) ◽  
pp. 6727-6737 ◽  
Author(s):  
Cong Xu ◽  
Zhihua Wang ◽  
Wubin Weng ◽  
Kaidi Wan ◽  
Ronald Whiddon ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document