Heat Transfer Influenced by Turbulent Airflow Inside an Axially Rotating Diffuser

The paper presents a study of heat transfer between the turbulent airflow and the inner wall surface of an axial diffuser rotating around its longitudinal axis. Heat transfer was assessed through the measurement of a time-dependent temperature field of the diffuser inner wall surface. Measurements of the instantaneous flow velocity components were performed by a laser–Doppler anemometry system, which delivered information on mean velocity components as well as on the turbulence intensity. A significant increase of all three mean velocity components was observed near the rotating diffuser wall in comparison with a non-rotating diffuser. Temperature field measurements were carried out by means of infrared thermography. The experiment showed a significant dependence of the temperature field on the turbulent flowfield induced by diffuser rotation. A strong influence of the flow separation and reattachment on the temperature distribution was observed, while rotation was found to suppress the occurrence of flow separation from the diffuser wall. Properties of the velocity field such as turbulent kinetic energy were directly coupled with the temperature distribution in order to gain the information on how to enhance or reduce heat transfer by changing the integral parameters of the diffuser (e.g. rotation frequency or amount of flow).     

Large Eddy Simulation of Turbulent Convective Cavity Flow

The large eddy simulation model with Smagorinsky subgrid-scale model was applied to two-dimensional turbulent convective cavity flow. The Reynolds number is lying from 1×10⁴ to 4×10⁵ and Archimedes number from 0 to 0.4. The simulation results were compared with the k?? model results and experimental results wherever possible. Flow results were in good agreement with experimental data across the mid-planes. Effects of Smagorinsky constant and grid resolution were investigated.     

Turbulent Flow and Dispersion of Inertial Particles in a Confined Jet Issued by a Long Cylindrical Pipe

In this work we examine first the flow field of a confined jet produced by a turbulent flow in a long cylindrical pipe issuing in an abrupt angle diffuser. Second, we examine the dispersion of inertial micro-particles entrained by the turbulent flow. Specifically, we examine how the particle dispersion field evolves in the multiscale flow generated by the interactions between the large-scale structures, which are geometry dependent, with the smaller turbulent scales issued by the pipe which are advected downstream. We use Large-Eddy-Simulation (LES) for the flow field and Lagrangian tracking for particle dispersion. The complex shape of the domain is modelled using the immersed-boundaries method. Fully developed turbulence inlet conditions are derived from an independent LES of a spatially periodic cylindrical pipe flow. The flow field is analyzed in terms of local velocity signals to determine spatial coherence and decay rate of the coherent K–H vortices and to make quantitative comparisons with experimental data on free jets. Particle dispersion is analyzed in terms of statistical quantities and also with reference to the dynamics of the coherent structures. Results show that the particle dynamics is initially dominated by the Kelvin–Helmholtz (K–H) rolls which form at the expansion and only eventually by the advected smaller turbulence scales.     

二维曲壁扩压器内湍流边界层分离的研究

本文应用二维双色四光束激光多普勒测速仪和压强探针样细地测量了二维非对称曲面扩张通道内的不可压湍流边界层分离流动，得到了时均速度和雷诺剪应力以及正反向间歇流动因子和静压分布。实验结果分析表明；湍流边界层分离时，沿边界层高度方向存在着明显的压强差。压强差的极小值对应于位移厚度曲率的极大值和瞬时间歇分离点。Ｂａｒｄｉｎａ对数尾迹律可以较好地描述瞬时间歇分离点之前的边界层速度分布，但无法描述分离的边界层速     

Large Eddy and Reynolds-Averaged Navier - Stokes Simulations of Turbulent Flow Over an Airfoil

A Large Eddy Simulation (LES) of turbulent flow over an airfoil near stall is performed. Results of the LES are compared with those of Reynolds-Averaged Navier-Stokes (RANS) simulations using two well-known turbulence models, namely the Baldwin-Lomax model and the Spalart-Allmaras model. The subgrid scale model used for the LES is the filtered structure function model. All simulations are performed using the same structured multi-block code. In order to reduce the CPU time, an implicit time stepping method is used for the LES. The purpose of this study is to show the possibilities and limitations of LES of complex flows associated with aeronautical applications using state of the art simulation techniques. Typical flow features are captured by the LES such as the adverse-pressure gradient and flow retardation. Visualization of instantaneous flow fields shows the typical streaky structures in the near-wall region.     

A hybrid dynamic Smagorinsky model for large eddy simulation

A hybrid dynamic subgrid-scale model (HDSM) pertaining to Large-eddy simulation (LES) has been developed. The coefficient obtained by German dynamic Smagorinsky model (DSM) was integrated with a new dynamic coefficient, based on the dynamic subgrid characteristic length and controlled by the subgrid-scale (SGS) motions. In HDSM, the characteristic wave number determining the characteristic length of the dynamic subgrid is calculated from a new energy weighted mean method when the subgrid scale turbulent kinetic energy and the dissipation wave number are known. The dissipation wave-number is derived from the SGS turbulent kinetic energy spectrum equation. The total dissipation rate spectrum equation is based on the Pao energy spectrum and local equilibrium assumption. The dynamic subgrid characteristic length could take into account the rapidly fluctuating small scale behaviours and the spatial variation of turbulent characteristics. HDSM was used to simulate the fully developed channel and turbulent flow past a circular cylinder, and to determine the impact of the dam-break flow on downstream structure. The HDSM is robust in respect to anisotropic mesh and is less sensitive to grid resolution, and would accurately describe the energy transfer from large-scale to SGS fluctuations and capture more fluctuations of turbulence with same meshes compared to the DSM.     

Large Eddy Simulation of Turbulent Flow Through Submerged Vegetation

Large Eddy Simulations (LES) are performed for an open channel flow through idealized submerged vegetation with a water depth (h) to plant height (h _p) ratio of h/h _p = 1.5 according to the experimental configuration of Liu et al. (J Geophys Res Earth Sci, 2008). They used a 1D laser Doppler velocimeter (LDV) to measure longitudinal and vertical velocities as well as turbulence intensities along several verticals in the flow and the data are used for the validation of the present simulations. The code MGLET is used to solve the filtered Navier–Stokes equations on a Cartesian non-uniform grid. In order to represent solid objects in the flow, the immersed boundary method is employed. The computational domain is idealized with a box containing 16 submerged circular cylinders and periodic boundary conditions are applied in both longitudinal and transverse directions. The predicted streamwise as well as vertical mean velocities are in good agreement with the LDV measurements. Furthermore, fairly good agreement is found between calculated and measured streamwise and vertical turbulence intensities. Large-scale flow structures of different shapes are present in the form of vortex rolls above the vegetation tops as well as locally generated trailing and von- Karman-type vortices due to flow separation at the free end and the sides of the cylinders. In this paper, the flow field is analyzed statistically and evidence is provided for the existence of these structures based on the LES.     

Turbulent transport of particles in a straight square duct

Turbulent flow through a duct of square cross-section gives rise to off-axis secondary flows, which are known to transfer momentum between fluid layers thereby flattening the velocity profile. The aim of this study is to investigate the role of the secondary flows in the transport and dispersion of particles suspended in a turbulent square duct flow. We have numerically simulated a flow through a square duct having a Reynolds number of Re_τ = 300 through discretization of the Navier–Stokes equations, and followed the trajectories of a large number of passive tracers and finite-inertia particles under a one-way coupling assumption. Snapshots of particle locations and statistics of single-particle and particle pair dispersion were analyzed. It was found that lateral mixing is enhanced for passive tracers and low-inertia particles due to the lateral advective transport that is absent in straight pipe and channels flows. Higher inertia particles accumulate close to the wall, and thus tend to mix more efficiently in the streamwise direction since a large number of the particles spend more time in a region where the mean fluid velocity is small compared to the bulk. Passive tracers tend to remain within the secondary swirling flows, circulating between the core and boundary of the duct.     

Prediction of Turbulent Separated Flow in a Turnaround Duct Using Wall Functions and Adaptivity

This paper presents an application of adaptive remeshing to the prediction of turbulent separated flows. The paper shows that the κ - ε model with wall functions can predict separated flows along smooth curved surfaces. Success is achieved if the wall functions exhibit values of y⁺ close to 30, and if meshes are fine enough to guarantee that wall function boundary conditions are grid converged. Adaptive remeshing proves to be a very cost effective tool in this context. The methodology is demonstrated on a problem possessing a closed form solution to establish the performance and reliability of the proposed approach. The method is then applied to prediction of turbulent flow in an annular, axisymmetric turnaround duct (TAD). Predictions from two computational models of the TAD are compared with experimental measurements. The importance of appropriate meshes to achieve grid independent solutions is demonstrated in both cases. Better agreement with measurements is obtained when partially developed profiles of u, κ, and ε are specified at the TAD inlet.     

Second-Order Conditional Moment Closure Simulations of Autoignition of an n-heptane Plume in a Turbulent Coflow of Heated Air

Autoignition of an n-heptane plume in a turbulent coflow of heated air has been studied using the conditional moment closure (CMC) method with a second-order closure for the conditional chemical source term. Two different methodologies have been considered: (i) the Taylor expansion method, in which the second order correction was based on the solution of the full covariance matrix for the 31 reactive species in the chemical mechanism and hence was not limited to a few selected reactions, and (ii) the conditional PDF method, in which only the temperature conditional variance equation has been solved and its PDF assumed to be a β-function. The results compare favorably with experiment in terms of autoignition location. The structure of the reaction zone in mixture fraction space has been explored. The relative performance of the two methodologies is discussed.     

Turbulent flow in a tube containing an offset rod

A numerical finite volume prediction method for arbitrary-shaped passages has been applied to the case of fully developed axial turbulent flow past a rod eccentrically placed in a circular tube. The numerical method was based on an orthogonal curvilinear mesh and employed an algebraic stress transport model to calculate the full three-dimensional velocity field directly from the governing partial differential equations. This study is one of a series of applications of this prediction method to a range of different non-circular passages that have been made in order to establish the capabilities and usefulness of this type of procedure. The present eccentric rod case was the subject of a comprehensive experimental investigation by Kacker¹ which has enabled a detailed comparison to be made between the present predictions and the measurements. This comparison included local distributions of axial velocity, wall shear stress and secondary velocities; and although found to be satisfactory overall, some differences in detail revealed possible shortcomings in the measurement of secondary flow. This, together with other previously reported cases, indicates, that, although the present method cannot be expected to replace experiment in providing turbulent passage flow data, it has an important role to play in interpreting and supplementing experiments.     

Velocity Measurements of a Free-Surface Turbulent Flow Penetrating a Porous Medium Composed of Uniform-Size Spheres

We present a laboratory experimental investigation of the interaction between the turbulent flow in an open channel and the turbulent flow within its very permeable bed. The bed was composed of uniform-size spheres packed in a cubic pattern. Fluid velocities were measured by Particle Image Velocimetry (PIV), which allowed us to investigate the spatial distribution of the time-averaged flow properties in the zone where they are strongly affected by the geometry of the rough bed. We investigate the effect of bed porosity on these flow properties by comparing the results of two experimental configurations: one with an impermeable bed composed of a single layer of spheres and another with a permeable bed composed of five layers. For the latter case, PIV measurements of velocities were also carried out inside two pores adjacent to the bed surface. This data provides an insight into the mechanisms of momentum transfer between the turbulent open channel flow and the turbulent flow within its very permeable bed.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Simulations of dispersed turbulent multiphase flow

Direct numerical simulations of homogeneous isotropic turbulence are used to investigate the effects of turbulence on the transport of particles in gas flows or bubbles in liquid flows. The inertia associated with the bubbles or the particles leads to locally strong concentrations of these in regions of instantaneously strong vorticity for bubbles or strain-rate for particles. This alters the average settling rates and other processes. If the mass-loading of the dispersed phase is significant a random “turbulent” flow is generated by the particle settling. A simple demonstration of this is given, showing the statistically axisymmetric character of this flow and how it can modify an ambient turbulent flow.     

Direct Numerical Simulation of Transitional Turbulent Flow in a Closed Rotor-Stator Cavity

The transitional turbulent regime in confined flow between a rotating and a stationary disc is studied using direct numerical simulation. Besides its fundamental importance as a three-dimensional prototype flow, such flows frequently arise in many industrial devices, especially in turbomachinary applications. The present contribution extends the DNS simulation into the turbulent flow regime, to a rotational Reynolds number Re =3 × 10⁵. An annular rotor-stator cavity of radial extension ΔR and height H, is considered with L = 4.72(L = ΔR/H) and Rm = 2.33 (Rm = (R ₁+ R ₀)/ΔR). The direct numerical simulation is performed by integrating the time-dependent Navier–Stokes equations until a statistically steady state is reached. A three-dimensional spectral method is used with the aim of providing both very accurate instantaneous fields and reliable statistical data. The instantaneous quantities are analysed in order to enhance our knowledge of the physics of turbulent rotating flows. Also, the results have been averaged so as to provide target turbulence data for any subsequent modelling attempts at reproducing the flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical Simulation of a Turbulent Flow in a Channel with Surface Mounted Cubes

In this paper we report on a fourth-order, spectro-consistent simulation of a complex turbulent flow. A spatial discretization of a convection-diffusion equation is termed spectro-consistent if the spectral properties of the convective and diffusive operators are preserved, i.e. convection skew-symmetric; diffusion symmetric positive definite. We consider a fully developed flow in a channel, where a matrix of cubes is placed at a wall of the channel. The Reynolds number (based on the channel width and the mean bulk velocity) is equal to Re = 13,000. The three-dimensional flow around the surface mounted cubes has served at a test case at the 6th ERCOFTAC/IAHR/COST workshop on refined flow modeling (Delft, June 1997). Here, mean velocity profiles as well as Reynolds stresses at various locations in the channel have been computed without using any turbulence models. The results agree well with the available experimental data.     

Spray Penetration in a Turbulent Flow

Analytical expressions for mass concentration of liquid fuel in a spray are derived taking into account the effects of gas turbulence, and assuming that the influence of droplets on gas is small (intitial stage of spray development). Beyond a certain distance the spray is expected to be fully dispersed. This distance is identified with the maximum spray penetration. Then the influence of turbulence on the spray stopping distance is discussed and the rms spray penetration is computed from a trajectory (Lagrangian) approach. Finally, the problem of spray penetration is investigated in a homogeneous two-phase flow regime taking into account the dispersion of spray away from its axis. It is predicted that for realistic values of spray parameters the spray penetration at large distances from the nozzle is expected to be proportional to t ^2/3 (in the case when this dispersion is not taken into account this distance is proportional to t ^1/2). The t ^2/3 law is supported by experimental observations for a high pressure injector. This revised version was published online in July 2006 with corrections to the Cover Date.     

Advanced Turbulence Modelling of Separated Flow in a Diffuser

The paper describes an investigation into the predictive performance of linear and non-linear eddy-viscosity models and differential stress-transport closures for separated flow in a nominally two-dimensional, asymmetric diffuser. The test case forms part of a broader collaborative exercise between academic and industrial partners. It is demonstrated that advanced turbulence models using strain-dependent coefficients and anisotropy-resolving closure offer tangible advantages in predictive capability, although the quality of their performance can vary significantly, depending on the details of closure approximations adopted. Certain features of the flow defy resolution by any of the closures investigated. In particular, no model resolves correctly the flow near the diffuser's inclined wall immediately downstream of the inlet corner, which may reflect the presence of a “flapping” motion associated with a highly-localised process of unsteady separation and reattachment. This revised version was published online in July 2006 with corrections to the Cover Date.