scholarly journals Credibility analysis of computational fluid dynamic simulations for compound channel flow

2013 ◽  
Vol 15 (3) ◽  
pp. 926-938 ◽  
Author(s):  
M. S. Filonovich ◽  
R. Azevedo ◽  
L. R. Rojas-Solórzano ◽  
J. B. Leal

In this paper, verification and validation of a turbulence closure model is performed for an experimental compound channel flow, where the velocity and turbulent fields were measured by a Laser Doppler Velocimeter (LDV). Detailed Explicit Algebraic Reynolds Stress Model (EARSM) simulations are reported. There are numerous methods and techniques available to evaluate the numerical uncertainty associated with grid resolution. The authors have adopted the Grid Convergence Index (GCI) approach. The velocity components, the turbulence kinetic energy (TKE), the dissipation rate and the Reynolds stresses were used as variables of interest. The GCI results present low values for the u velocity component, but higher values in what concerns the v velocity component and w velocity component (representing secondary flows) and for Reynolds stresses RSxy and RSyz. This indicates that the mean flow has converged but the turbulent field and secondary flows still depend on grid resolution. Based on GCI values distribution, the medium and fine meshes were further refined. In addition to GCI analysis, the authors have performed linear regression analysis for estimating the mesh quality in what concerns small value variables. Comparison of numerical and experimental results shows good agreement.

Author(s):  
M. Kanniche ◽  
R. Boudjemadi ◽  
F. Déjean ◽  
F. Archambeau

The flow in a linear turbine cascade (Gregory-Smith et al. (1990)) is numerically investigated using a Reynolds Stress Turbulence closure. A particular attention is given to secondary flows where the normal Reynolds stresses are expected to play an important role. The most classical turbulence closure, the k-epsilon model uses the Boussinesq Eddy Viscosity concept which assumes an isotropic turbulent viscosity. The Reynolds stresses are then related to local velocity gradients by this isotropic eddy viscosity. Corollary, the principal axes of the Reynolds stress tensor are colinear with those of the mean strain tensor. The advantage of Reynolds Stress Turbulence closure is the calculation of Reynolds stresses by their own individual transport equations. This leads to a more realistic description of the turbulence and of its dependance on the mean flow. The most classical Second Order turbulence model (Launder et al. (1975)) is applied to a linear turbine cascade, and the results are compared to secondary velocity and turbulence measurements at cross-passage planes.


2022 ◽  
Vol 933 ◽  
Author(s):  
T. Bon ◽  
J. Meyers

Recent studies have demonstrated that large secondary motions are excited by surface roughness with dominant spanwise length scales of the order of the flow's outer length scale. Inspired by this, we explore the effect of spanwise heterogeneous surface temperature in weakly to strongly stratified closed channel flow (at $Ri_\tau =120$ , 960; $Re_\tau = 180$ , 550) with direct numerical simulations. The configuration consists of equally sized strips of high and low temperature at the lower and upper boundaries, while an overall stable stratification is induced by imposing an average temperature difference between the top and bottom. We consider the influence of the width of the strips ( ${\rm \pi} /8 \leq \lambda /h \leq 4{\rm \pi} $ ), Reynolds number, stability and upper boundary condition on the mean flow structure, skin friction and heat transfer. Results indicate that secondary flows are excited, with alternating high- and low-momentum pathways and vortices, similar to the patterns induced by spanwise heterogeneous surface roughness. We find that the impact of the surface heterogeneity on the outer layer depends strongly on the spanwise heterogeneity length scale of the surface temperature. Comparison to stable channel flow with uniform temperature reveals that the heterogeneous surface temperature increases the global friction coefficient and reduces the global Nusselt number in most cases. However, for the high-Reynolds cases with $\lambda /h \geq {\rm \pi} /2$ , we find a reduction of the friction coefficient. At stronger stability, the vertical extent of the vortices is reduced and the impact of the heterogeneous temperature on momentum and heat transfer is smaller.


1993 ◽  
Vol 115 (1) ◽  
pp. 109-114 ◽  
Author(s):  
T.-M. Liou ◽  
Y.-Y. Wu ◽  
Y. Chang

Laser-Doppler velocimeter measurements of mean velocities, turbulence intensities, and Reynolds stresses are presented for periodic fully developed flows in a channel with square rib-disturbed walls on two opposite sides. Quantities such as the vorticity thickness and turbulent kinetic energy are used to characterize the flow. The investigated flow was periodic in space. The Reynolds number based on the channel hydraulic diameter was 3.3×104. The ratios of pitch to rib-height and rib-height to chamber-height were 10 and 0.133, respectively. Regions where maximum and minimum Reynolds stress and turbulent kinetic energy occurred were identified from the results. The growth rate of the shear layers of the present study was compared with that of a backward-facing step. The measured turbulence anisotropy and structure parameter distribution were used to examine the basic assumptions embedded in the k–ε and k–ε–A models. For a given axial station, the peak axial mean-velocity was found not to occur at the center point. The secondary flow was determined to be Prandtl’s secondary flow of the second kind according to the measured streamwise mean vorticity and its production term.


Author(s):  
Peter Ajersch ◽  
Jian-Ming Zhou ◽  
Steven Ketler ◽  
Martha Salcudean ◽  
Ian S. Gartshore

The fluid mechanics and heat transfer characteristics of film cooling is three-dimensional and highly complex. To better understand this problem, an experimental study was conducted in a low speed wind tunnel on a row of six rectangular jets injected at 90° to the crossflow (main stream flow). The jet-to-crossflow velocity ratios (blowing ratios) examined were 0.5, 1.0, and 1.5, and the jet spacing-to-jet width ratio was 3.0. No significant temperature difference between jet and crossflow air was introduced. Mean velocities, and six flow stresses were measured using a three-component laser Doppler velocimeter operating in coincidence-mode. Seeding of both jet and cross-stream air was achieved with a commercially available smoke generator. Flow statistics are reported in the form of vector plots, contours, and x-y graphs, showing velocity, turbulence intensity, and Reynolds stresses. To complement the detailed measurements, flow visualization was accomplished by transmitting the laser beam through a cylindrical lens, thereby generating a narrow, intense sheet of light. Jet air only was seeded with smoke, which was illuminated in the plane of the light sheet. Therefore, it was possible to record on video tape the trajectory and penetration of the jets in the crossflow. Selected still images from the recordings are presented. Numerical simulations of the observed flow field were made by using a multi-grid, segmented, k-ε CFD code. Special near wall treatment included a non-isotropic formulation for the effective viscosity, a low Re model for k, and an algebraic model for the length scale. Comparisons between the measured and computed velocities show good agreement for the non-uniform mean flow at the jet exit plane. Velocities and stresses on the jet centerline downstream of the orifice are less well predicted, probably because of inadequate turbulence modeling, while values off the centerline match those of the experiments much more closely.


1997 ◽  
Vol 119 (2) ◽  
pp. 330-342 ◽  
Author(s):  
P. Ajersch ◽  
J.-M. Zhou ◽  
S. Ketler ◽  
M. Salcudean ◽  
I. S. Gartshore

The fluid mechanics and heat transfer characteristics of film cooling are three-dimensional and highly complex. To understand this problem better, an experimental study was conducted in a low-speed wind tunnel on a row of six rectangular jets injected at 90 deg to the crossflow (mainstream flow). The jet-to-crossflow velocity ratios (blowing ratios) examined were 0.5, 1.0, and 1.5, and the jet spacing-to-jet width ratio was 3.0. No significant temperature difference between jet and crossflow air was introduced. Mean velocities and six flow stresses were measured using a three-component laser-Doppler velocimeter operating in coincidence mode. Seeding of both jet and cross-stream air was achieved with a commercially available smoke generator. Flow statistics are reported in the form of vector plots, contours, and x-y graphs, showing velocity, turbulence intensity, and Reynolds stresses. To complement the detailed measurements, flow visualization was accomplished by transmitting the laser beam through a cylindrical lens, thereby generating a narrow, intense sheet of light. Jet air only was seeded with smoke, which was illuminated in the plane of the light sheet. Therefore, it was possible to record on video tape the trajectory and penetration of the jets in the crossflow. Selected still images from the recordings are presented. Numerical simulations of the observed flow field were made by using a multigrid, segmented, k–ε CFD code. Special near-wall treatment included a nonisotropic formulation for the effective viscosity, a low-Re model for k, and an algebraic model for the length scale. Comparisons between the measured and computed velocities show good agreement for the nonuniform mean flow at the jet exit plane. Velocities and stresses on the jet centerline downstream of the orifice are less well predicted, probably because of inadequate turbulence modeling, while values off the centerline match those of the experiments much more closely.


Author(s):  
Thomas Ludwig Kaiser ◽  
Kilian Oberleithner

Abstract In this paper a new method is introduced to model the transport of entropy waves and equivalence ratio fluctuations in turbulent flows. The model is based on the Navier-Stokes equations and includes a transport equation for a passive scalar, which may stand for entropy or equivalence ratio fluctuations. The equations are linearized around the mean turbulent fields, which serve as the input to the model in addition to a turbulent eddy viscosity, which accounts for turbulent diffusion of the perturbations. Based on these inputs, the framework is able to predict the linear response of the flow velocity and passive scalar to harmonic perturbations that are imposed at the boundaries of the computational domain. These in this study are fluctuations in the passive scalar and/or velocities at the inlet of a channel flow. The code is first validated against analytic results, showing very good agreement. Then the method is applied to predict the convection, mean flow dispersion and turbulent mixing of passive scalar fluctuations in a turbulent channel flow, which has been studied in previous work with Direct Numerical Simulations (DNS). Results show that our code reproduces the dynamics of coherent passive scalar transport in the DNS with very high accuracy and low numerical costs, when the DNS mean flow and Reynolds stresses are provided. Furthermore, we demonstrate that turbulent mixing has a significant effect on the transport of the passive scalar fluctuations. Finally, we apply the method to explain experimental observations of transport of equivalence ratio fluctuations in the mixing duct of a model burner.


2008 ◽  
Vol 130 (6) ◽  
Author(s):  
Martin Agelinchaab ◽  
Mark F. Tachie

A particle image velocimetry is used to study the mean and turbulent fields of separated and redeveloping flow over square, rectangular, and semicircular blocks fixed to the bottom wall of an open channel. The open channel flow is characterized by high background turbulence level, and the ratio of the upstream boundary layer thickness to block height is considerably higher than in prior experiments. The variation of the Reynolds stresses along the dividing streamlines is discussed within the context of vortex stretching, longitudinal strain rate, and wall damping. It appears that wall damping is a more dominant mechanism in the vicinity of reattachment. In the recirculation and reattachment regions, profiles of the mean velocity, turbulent quantities, and transport terms are used to document the salient features of block geometry on the flow. The flow characteristics in these regions strongly depend on block geometry. Downstream of reattachment, a new shear layer is formed, and the redevelopment of the shear layer toward the upstream open channel boundary layer is studied using the boundary layer parameters and Reynolds stresses. The results show that the mean flow rapidly redeveloped so that the Clauser parameter recovered to its upstream value at 90 step heights downstream of reattachment. However, the rate of development close to reattachment strongly depends on block geometry.


2012 ◽  
Vol 711 ◽  
pp. 161-170 ◽  
Author(s):  
Siddharth Talapatra ◽  
Joseph Katz

AbstractMicroscopic holographic PIV performed in an optically index-matched facility resolves the three-dimensional flow in the inner part of a turbulent channel flow over a rough wall at Reynolds number ${\mathit{Re}}_{\tau } = 3520$. The roughness consists of uniformly distributed pyramids with normalized height of ${ k}_{s}^{+ } = 1. 5{k}^{+ } = 97$. Distributions of mean flow and Reynolds stresses agree with two-dimensional PIV data except very close to the wall (${\lt }0. 7k$) owing to the higher resolution of holography. Instantaneous realizations reveal that the roughness sublayer is flooded by low-lying spanwise and groove-parallel vortical structures, as well as quasi-streamwise vortices, some quite powerful, that rise at sharp angles. Conditional sampling and linear stochastic estimation (LSE) reveal that the prevalent flow phenomenon in the roughness sublayer consists of interacting U-shaped vortices, conjectured in Hong et al. (J. Fluid Mech., 2012, doi:10.1017/jfm.2012.403). Their low-lying base with primarily spanwise vorticity is located above the pyramid ridgeline, and their inclined quasi-streamwise legs extend between ridgelines. These structures form as spanwise vorticity rolls up in a low-speed region above the pyramid’s forward face, and is stretched axially by the higher-speed flow between ridgelines. Ejection induced by interactions among legs of vortices generated by neighbouring pyramids appears to be the mechanism that lifts the quasi-streamwise vortex legs and aligns them preferentially at angles of $54\textdegree \text{{\ndash}} 63\textdegree $ to the streamwise direction.


2011 ◽  
Vol 667 ◽  
pp. 1-37 ◽  
Author(s):  
JIARONG HONG ◽  
JOSEPH KATZ ◽  
MICHAEL P. SCHULTZ

Utilizing an optically index-matched facility and high-resolution particle image velocimetry measurements, this paper examines flow structure and turbulence in a rough-wall channel flow for Reτ in the 3520–5360 range. The scales of pyramidal roughness elements satisfy the ‘well-characterized’ flow conditions, with h/k ≈ 50 and k+ = 60 ~ 100, where h is half height of the channel and k is the roughness height. The near-wall turbulence measurements are sensitive to spatial resolution, and vary with Reynolds number. Spatial variations in the mean flow, Reynolds stresses, as well as the turbulent kinetic energy (TKE) production and dissipation rates are confined to y < 2k. All the Reynolds stress components have local maxima at slightly higher elevations, but the streamwise-normal component increases rapidly at y < k, peaking at the top of the pyramids. The TKE production and dissipation rates along with turbulence transport also peak near the wall. The spatial energy and shear spectra show an increasing contribution of large-scale motions and a diminishing role of small motions with increasing distance from the wall. As the spectra steepen at low wavenumbers, they flatten and develop bumps in wavenumbers corresponding to k − 3k, which fall in the dissipation range. Instantaneous realizations show that roughness-scale eddies are generated near the wall, and lifted up rapidly by large-scale structures that populate the outer layer. A linear stochastic estimation-based analysis shows that the latter share common features with hairpin packets. This process floods the outer layer with roughness-scale eddies, in addition to those generated by the energy-cascading process. Consequently, although the imprints of roughness diminish in the outer-layer Reynolds stresses, consistent with the wall similarity hypothesis, the small-scale turbulence contains a clear roughness signature across the entire channel.


2014 ◽  
Vol 748 ◽  
Author(s):  
Julio M. Barros ◽  
Kenneth T. Christensen

AbstractThe characteristics of a turbulent boundary layer overlying a complex roughness topography were explored with stereo particle-image velocimetry measurements in the wall-normal–spanwise $(y\mbox{--}z)$ plane. The roughness under consideration was replicated from a turbine blade damaged by deposition of foreign materials containing a broad range of topographical scales arranged in a highly irregular manner. The single-point turbulence statistics displayed strong spanwise heterogeneity, in particular spanwise-alternating low- and high-momentum flow pathways in the mean flow marked by enhanced Reynolds stresses and turbulent kinetic energy. The spanwise regions between high- and low-momentum flow pathways were occupied by swirling motions, suggesting the generation and sustainment of turbulent secondary flows due to the spanwise heterogeneity of the complex roughness under consideration. Similar observations were recently reported for more ordered spanwise roughness transitions by Nugroho, Hutchins & Monty (Intl J. Heat Fluid Flow vol. 41, 2013, pp. 90–102) and Willingham et al. (Phys. Fluids vol. 26, 2014, 025111).


Sign in / Sign up

Export Citation Format

Share Document