Assessment of Low-Re Turbulence Models and Analysis of Turbulent Flow in Porous Media Consisting of Square Cylinders With Different Diameter Ratios

2020 ◽  
Vol 143 (1) ◽  
Author(s):  
Alejandro Alonzo-Garcia ◽  
Ana T. Mendoza-Rosas ◽  
Martín A. Díaz-Viera ◽  
Sergio A. Martínez-Delgadillo ◽  
Edgar G. Martínez-Mendoza

Abstract This paper presents a comparative study of volume average predictions between low-Reynolds-number (LRN) turbulence models: Abe–Kondoh–Nagano (AKN), Lam–Bremhorst, Yang–Shih, standard k–ϵ, and k–ω. A porous medium, which represents conditions in which the flow path changes rapidly, was defined as an infinite array of square cylinders. In addition, to explore the effect of particle size on the rapid expansion and contraction of the flow paths, the diameter ratio (DR) of the square cylinders was systematically varied from 0.2 to 0.8. This generalization revealed new insights into the flow. The Reynolds number (ReD) covered a turbulent range of 500 to 500 × 103, and the porosity ϕ was varied from 0.27 to 0.8. The correlations of the turbulent kinetic energy (k), its dissipation rate (ε), and macroscopic pressure gradient as a function of ϕ, which are useful in macroscopic turbulence modeling, are presented. The results show that the AKN model yields better predictions of the volume-averaged flow parameters because it is better suited to reproduce recirculation zones. For all the DRs, at high ϕ, the distances between walls are high, and the interstitial velocities are low. Consequently, wake flows are produced, and energy losses by friction are moderate. As the flow becomes increasingly bound, the wakes are suppressed and disrupted, and k and ε increase owing to shear layer interactions and frictional forces. Distinctive low-velocity recirculation patterns appear inside pores depending on DR.

2021 ◽  
Vol 409 ◽  
pp. 158-178
Author(s):  
Abdelkader Feddal ◽  
Abbes Azzi ◽  
Ahmed Zineddine Dellil

This paper deals with studying numerically two circular turbulent jets impinging on a flat surface with a low velocity cross flow by using ANSYS CFX 16.2, with the aim of proving the effect ofReynolds number on the flow demeanor in a vertical circular free turbulent jet with cross flow. Five turbulence models of the RANS (Reynolds Averaged Navier–Stokes) approach were tested and the k -ω SST model was chosen to validate CFD results with the experimental data. Average velocity profiles, velocity and turbulent kinetic energy contours and streamlines are presented for four case configurations. In the first three cases, the following parameters have been varied: Reynolds number at the level of the two jets ( ), wind velocity at the level of the cross-flow ( ), and the distance between the two jets (S = 45mm, 90mm and 135mm). In the last case, a new configuration of the phenomenon not yet studied so far was treated, where horizontal cross-flows were introduced from both sides in order to simulate gusts of wind disrupting a VSTOL aircraft which tries to operate close to the ground. This case was carried out for Reynolds number based on the crossflow of 4 104, 10 104 and 20 104 .The numerical results obtained show that the deflection of the jets is minimal when the Reynolds number at the level of the jets is greater than that of the cross-flow. The increase of Reynolds number at the level of the cross-flow reveals a significant deviation of the two jets with an intensity which always remains less for the second jet. As for the space parameter between the two jets, it turns out that the fact of further spacing the two jets makes the first jet even more vulnerable and leads to a greater deflection. Finally, the simulation of the wind gusts from the front and the back caused a zone of turbulence which resulted from a form of "interlacing" of the two jets under the effect of the transverse current imposed by the two sides.


2011 ◽  
Vol 680 ◽  
pp. 361-397 ◽  
Author(s):  
C. M. SEWATKAR ◽  
ATUL SHARMA ◽  
AMIT AGRAWAL

A numerical study of flow across a row of transversely oscillating square cylinders (of diameter d) has been undertaken using the lattice Boltzmann method, for a better understanding of fluid–structure interaction problems. The effects of cylinder oscillation frequency ratio (fe/fo, where fe is the cylinder oscillation frequency and fo is the corresponding vortex shedding frequency for stationary row of cylinders), amplitude ratio (A/d), non-dimensional spacing between the cylinders (s/d) and Reynolds number (Re) on ensuing flow regimes and flow parameters have been studied to understand the flow physics. Six different flow regimes observed in this study are the quasi-periodic non-lock-on-I, synchronous lock-on, quasi-periodic lock-on, quasi-periodic non-lock-on-II, synchronous non-lock-on and chaotic non-lock-on. It is observed that the range of the lock-on regime depends upon the relative dominance of incoming flow and cylinder motion. Although the lock-on regime in the case of Re = 80, s/d = 4 and A/d = 0.2 is substantially larger as compared to that for a single oscillating cylinder, the range of the lock-on regime shrinks with a reduction in the cylinder spacing, increase in the Reynolds number or decrease in the oscillation amplitude. It is also observed that the wake interaction behind the cylinders weakens with an increase in fe/fo, Re, A/d or s/d, leading to the formation of independent wakes and synchronous nature of the flow. For fe/fo ≥ 1.2, independent and intact oscillating wakes are noted and an additional frequency (wake oscillation frequency) is obtained in the time series of the lift coefficient. Although it was expected that the complexity in the wake interaction would increase with cylinder oscillation or amplitude ratio, an opposite effect (that is, formation of independent wakes) is noted from the results.


1999 ◽  
Vol 121 (4) ◽  
pp. 824-833 ◽  
Author(s):  
A. Chernobrovkin ◽  
B. Lakshminarayana

Variation of the flow Reynolds number between the take off and cruise conditions significantly affects the boundary layer development on low-pressure turbine blading. A decreased Reynolds number leads to the flow separation on the suction surface of the blading and increased losses. A numerical simulation has been carried out to assess the ability of a Navier-Stokes solver to predict transitional flows in a wide range of Reynolds numbers and inlet turbulence intensities. A number of turbulence models (including the Algebraic Reynolds Stress Model) and transition models have been employed to analyze the reliability and accuracy of the numerical simulation. A comparison between the prediction and the experimental data reveals good correlation. However, the analysis shows that the artificial dissipation in the numerical solver may have a profound effect on the prediction of the transition in a separated flow.


2001 ◽  
Vol 123 (3) ◽  
pp. 640-648 ◽  
Author(s):  
Raymond M. Jones ◽  
Albert D. Harvey ◽  
Sumanta Acharya

In this study, the predictive performance of six different two-equation turbulence models on the flow in an unbaffled stirred tank has been investigated. These models include the low Reynolds number k-ε model of Rodi, W., and Mansour, N. N., “Low Reynolds Number k-ε Modeling With the Aid of Direct Simulation Data,” J. Fluid Mech., Vol. 250, pp. 509–529, the high and low Reynolds number k-ω models of Wilson, D. C., 1993, Turbulence Modeling for CFD, DCW Industries, La Canada, CA., the RNG k-ε model, and modified k-ω and k-ε models which incorporate a correction for streamline curvature and swirl. Model results are compared with experimental laser Doppler velocimetry (LDV) data for the turbulent velocity field in an unbaffled tank with a single paddle impeller. An overall qualitative agreement has been found between the experimental and numerical results with poor predictions observed in some parts of the tank. Discrepancies in model predictions are observed in the anisotropic regions of the flow such as near the impeller shaft and in the impeller discharge region where the model overpredicts the radial velocity component. These results are discussed and a strategy for improving two-equation models for application to impeller stirred tanks is proposed.


Author(s):  
Karsten Tawackolian ◽  
Martin Kriegel

AbstractThis study looks to find a suitable turbulence model for calculating pressure losses of ventilation components. In building ventilation, the most relevant Reynolds number range is between 3×104 and 6×105, depending on the duct dimensions and airflow rates. Pressure loss coefficients can increase considerably for some components at Reynolds numbers below 2×105. An initial survey of popular turbulence models was conducted for a selected test case of a bend with such a strong Reynolds number dependence. Most of the turbulence models failed in reproducing this dependence and predicted curve progressions that were too flat and only applicable for higher Reynolds numbers. Viscous effects near walls played an important role in the present simulations. In turbulence modelling, near-wall damping functions are used to account for this influence. A model that implements near-wall modelling is the lag elliptic blending k-ε model. This model gave reasonable predictions for pressure loss coefficients at lower Reynolds numbers. Another example is the low Reynolds number k-ε turbulence model of Wilcox (LRN). The modification uses damping functions and was initially developed for simulating profiles such as aircraft wings. It has not been widely used for internal flows such as air duct flows. Based on selected reference cases, the three closure coefficients of the LRN model were adapted in this work to simulate ventilation components. Improved predictions were obtained with new coefficients (LRNM model). This underlined that low Reynolds number effects are relevant in ventilation ductworks and give first insights for suitable turbulence models for this application. Both the lag elliptic blending model and the modified LRNM model predicted the pressure losses relatively well for the test case where the other tested models failed.


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.


2010 ◽  
Vol 132 (7) ◽  
Author(s):  
Henrique Stel ◽  
Rigoberto E. M. Morales ◽  
Admilson T. Franco ◽  
Silvio L. M. Junqueira ◽  
Raul H. Erthal ◽  
...  

This article describes a numerical and experimental investigation of turbulent flow in pipes with periodic “d-type” corrugations. Four geometric configurations of d-type corrugated surfaces with different groove heights and lengths are evaluated, and calculations for Reynolds numbers ranging from 5000 to 100,000 are performed. The numerical analysis is carried out using computational fluid dynamics, and two turbulence models are considered: the two-equation, low-Reynolds-number Chen–Kim k-ε turbulence model, for which several flow properties such as friction factor, Reynolds stress, and turbulence kinetic energy are computed, and the algebraic LVEL model, used only to compute the friction factors and a velocity magnitude profile for comparison. An experimental loop is designed to perform pressure-drop measurements of turbulent water flow in corrugated pipes for the different geometric configurations. Pressure-drop values are correlated with the friction factor to validate the numerical results. These show that, in general, the magnitudes of all the flow quantities analyzed increase near the corrugated wall and that this increase tends to be more significant for higher Reynolds numbers as well as for larger grooves. According to previous studies, these results may be related to enhanced momentum transfer between the groove and core flow as the Reynolds number and groove length increase. Numerical friction factors for both the Chen–Kim k-ε and LVEL turbulence models show good agreement with the experimental measurements.


2002 ◽  
Vol 12 (03) ◽  
pp. 393-406 ◽  
Author(s):  
ZI-NIU WU ◽  
SONG FU

The k-epsilon turbulence model for incompressible flow involves two advection–diffusion equations plus point-source terms. We propose a new method for positivity analysis. This method uses an iterative procedure combined with an operator splitting. With this method we recover the well-known positivity result for the standard high Reynolds number model. Most importantly, we are able to prove the positivity result for general low Reynolds number k-epsilon models.


1997 ◽  
Vol 119 (4) ◽  
pp. 794-801 ◽  
Author(s):  
J. Luo ◽  
B. Lakshminarayana

The boundary layer development and convective heat transfer on transonic turbine nozzle vanes are investigated using a compressible Navier–Stokes code with three low-Reynolds-number k–ε models. The mean-flow and turbulence transport equations are integrated by a four-stage Runge–Kutta scheme. Numerical predictions are compared with the experimental data acquired at Allison Engine Company. An assessment of the performance of various turbulence models is carried out. The two modes of transition, bypass transition and separation-induced transition, are studied comparatively. Effects of blade surface pressure gradients, free-stream turbulence level, and Reynolds number on the blade boundary layer development, particularly transition onset, are examined. Predictions from a parabolic boundary layer code are included for comparison with those from the elliptic Navier–Stokes code. The present study indicates that the turbine external heat transfer, under real engine conditions, can be predicted well by the Navier–Stokes procedure with the low-Reynolds-number k–ε models employed.


Sign in / Sign up

Export Citation Format

Share Document