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.     

Modelling of Advection-Dominated Transport in a Porous Column with a Time-Decaying Flow Field

This paper examines the problem of the advective-dispersive movement of a non-decaying, inert chemical dye solution through the pore space of a fluid saturated porous column. The objective of the paper is to present a complete study of the one-dimensional advective-dispersive transport problem by considering certain analytical solutions, experimental results and their comparisons with specific computational simulations. Dye concentrations obtained by means of an image processing method are used in conjunction with an analytical solution to identify the hydrodynamic dispersion coefficient that governs the advective-dispersive transport problem. The experimental results and identified parameters are also used to assess the computational estimates derived from several stabilized computational schemes available in the literature, for examining advection-dominated transport processes in porous media.     

DNS and Modelling of Passive Scalar Transport in Turbulent Channel Flow with a Focus on Scalar Dissipation Rate Modelling

A direct numerical simulation of turbulent channel flow with an imposed mean scalar gradient is analyzed with a focus on passive scalar flux modelling and in particular the treatment of the passive scalar dissipation equation. The Prandtl number is 0.71 and the Reynolds number based on the wall friction velocity and the channel half width is 265. Budgets are presented for the passive scalar variance and its dissipation rate, as well as for the individual scalar flux components. These form a basis for a discussion of modelling issues related to explicit algebraic scalar flux modelling. This revised version was published online in July 2006 with corrections to the Cover Date.     

An Experimental Study of the Rotational Effects on Separated Turbulent Flow During Stall Delay

Three-dimensional velocity fields were measured using tomographic particle image velocimetry (Tomo-PIV) on a model of the blade of a small-scale horizontal axis wind turbine (HAWT) to study the effects of rotation on separated turbulent flows during stall delay at a global tip speed ratio (TSR) of 3 and a Reynolds number of 4800. The flow fields on a static airfoil were also measured at a similar angle-of-attack (AOA) and Reynolds number for comparison. It was observed that the blade’s rotation in the streamwise direction significantly affected both the mean flow and the turbulence statistics over the suction surface. The mean velocity fields revealed that, different from the airfoil flow at large AOA, the recirculation region with reversed flow did not exist on the suction surface of the blade and the flow was rather attached. Mean spanwise flow from blade’s root to its tip was also generated by the rotation. The mean vorticity vector of the blade flow was found to be tilted in the rotational direction of the blade, as well as in the wall-normal direction. Of particular effects of the rotation on Reynolds stresses were the enhancement of 〈w ²〉 and the creation of strong 〈v w〉. The production of Reynolds stresses was also affected by blade’s rotation directly through the rotational production terms and indirectly by dramatically changing the fluctuating velocity fields. The distribution of enstrophy was observed to be modified by rotation, too.     

A Direct Numerical Study of Transport and Anisotropy in a Stably Stratified Turbulent Flow with Uniform Horizontal Shear

Direct numerical simulations of homogeneous turbulence in stably stratified shear flow have been performed to aid the understanding of turbulence and turbulent mixing in geophysical flow. Two cases are compared. In the first case, which has been studied in the past, the mean velocity has vertical shear and the mean density is vertically stably stratified. In the second case, which has not been studied systematically before, the mean velocity has horizontal shear and the mean density is again vertically stably stratified. The critical value of the gradient Richardson number, for which a constant turbulence level is obtained, is found to be an order of magnitude larger in the horizontal shear case. The turbulent transport coefficients of momentum and vertical mass transfer are also an order of magnitude larger in the horizontal shear case. The anisotropy of the turbulence intensities are found to be in the range expected of flows with mean shear with no major qualitative change in the range of Richardson numbers studied here. However, the anisotropy of the turbulent dissipation rate is strongly affected by stratification with the vertical component dominating the others. This revised version was published online in July 2006 with corrections to the Cover Date.     

Network Model of Flow,Transport and Biofilm Effects in Porous Media

In this paper, we develop a network model to determine porosity and permeability changes in a porous medium as a result of changes in the amount of biomass. The biomass is in the form of biofilms. Biofilms form when certain types of bacteria reproduce, bond to surfaces, and produce extracellular polymer (EPS) filaments that link together the bacteria. The pore spaces are modeled as a system of interconnected pipes in two and three dimensions. The radii of the pipes are given by a lognormal probability distribution. Volumetric flow rates through each of the pipes, and through the medium, are determined by solving a linear system of equations, with a symmetric and positive definite matrix. Transport through the medium is modeled by upwind, explicit finite difference approximations in the individual pipes. Methods for handling the boundary conditions between pipes and for visualizing the results of numerical simulations are developed. Increases in biomass, as a result of transport and reaction, decrease the pipe radii, which decreases the permeability of the medium. Relationships between biomass accumulation and permeability and porosity reduction are presented.     

Computational Modelling of Non-Newtonian Effects on Flow in Channels with Moving Wall Indentations

Non-Newtonian effects in a channel with moving wall indentations are assessed numerically by a finite volume method for solving the unsteady incompressible Navier-Stokes equations and using a power-law model exhibiting shear thinning viscosity and Casson's model as the constitutive equations for the non-Newtonian fluid. The computations show that for a non-Newtonian fluid, there are differences in the velocity profiles and in the structure and size of the reversed flow regions as compared with the corresponding Newtonian fluid. The comparison of non-Newtonian and Newtonian wall shear stress reveals a slight decrease in the magnitude on the average for the non-Newtonian case, eventually resulting in the strength of the “wave train” being slightly weaker than those corresponding to a Newtonian fluid.     

Modelling the Flow of Yield-Stress Fluids in Porous Media

Yield-stress is a problematic and controversial non-Newtonian flow phenomenon. In this article, we investigate the flow of yield-stress substances through porous media within the framework of pore-scale network modelling. We also investigate the validity of the Minimum Threshold Path (MTP) algorithms to predict the pressure yield point of a network depicting random or regular porous media. Percolation theory as a basis for predicting the yield point of a network is briefly presented and assessed. In the course of this study, a yield-stress flow simulation model alongside several numerical algorithms related to yield-stress in porous media were developed, implemented and assessed. The general conclusion is that modelling the flow of yield-stress fluids in porous media is too difficult and problematic. More fundamental modelling strategies are required to tackle this problem in the future.     

Anisotropy Invariant Reynolds Stress Model of Turbulence (AIRSM) and its Application to Attached and Separated Wall-Bounded Flows

Numerical predictions with a differential Reynolds stress closure, which in its original formulation explicitly takes into account possible states of turbulence on the anisotropy-invariant map, are presented. Thus the influence of anisotropy of turbulence on the modeled terms in the governing equations for the Reynolds stresses is accounted for directly. The anisotropy invariant Reynolds stress model (AIRSM) is implemented and validated in different finite-volume codes. The standard wall-function approach is employed as initial step in order to predict simple and complex wall-bounded flows undergoing large separation. Despite the use of simple wall functions, the model performed satisfactory in predicting these flows. The predictions of the AIRSM were also compared with existing Reynolds stress models and it was found that the present model results in improved convergence compared with other models. Numerical issues involved in the implementation and application of the model are also addressed.     

Investigation of Anisotropy-Resolving Turbulence Models by Reference to Highly-Resolved LES Data for Separated Flow

The predictive properties of several non-linear eddy-viscosity models are investigated by reference to highly-resolved LES data obtained by the authors for an internal flow featuring massive separation from a curved surface. The test geometry is a periodic segment of a channel constricted by two-dimensional (2D) `hills' on the lower wall. The mean-flow Reynolds number is 21560. Periodic boundary conditions are applied in the streamwise and spanwise directions. This makes the statistical properties of the simulated flow genuinely 2D and independent from boundary conditions, except at the walls. The simulation was performed on a high-quality, 5M-node grid. The focus of the study is on the exploitation of the LES data for the mean-flow, Reynolds stresses and macro-length-scale. Model solutions are first compared with the LES data, and selected models are then subjected to a-priori studies designed to elucidate the role of specific model fragments in the non-linear stress-strain/vorticity relation and their contribution to observed defects in the mean-flow and turbulence fields. The role of the equation governing the length-scale, via different surrogate variables, is also investigated. It is shown that, while most non-linear models overestimate the separation region, due mainly to model defects that result in insufficient shear stress in the separated shear layer, model forms can be derived which provide a satisfactory representation of the flow. One such model is identified. This combines a particular quadratic constitutive relation with a wall-anisotropy term, a high-normal-strain correction and a new form of the equation for the specific dissipation ω = ∈/k. This revised version was published online in July 2006 with corrections to the Cover Date.     

Assessing the Predictive Capabilities of Combustion LES as Applied to the Sydney Flame Series

Large Eddy Simulation (LES) and flamelet-based combustion models were applied to four bluff-body stabilized nonpremixed and partially premixed flames selected from the Sydney flame series, based on Masri’s bluff-body test rig (University of Sydney). Three related non-reacting flow cases were also investigated to assess the performance of the LES solver. Both un-swirled and swirled cases were studied exhibiting different flow features, such as recirculation, jet precessing and vortex breakdown. Due to various fuel compositions, flow rates and swirl numbers, the combustion characteristics of the flames varied greatly. On six meshes with different blocking structure and mesh sizes, good prediction of flow and scalar fields using LES/flamelet approaches and known fuel and oxidizer mass fluxes was achieved. The accuracy of predictions was strongly influenced by the combustion model used. All flames were calculated using at least two modeling strategies. Starting with calculations of isothermal flow cases, simple single flamelet based calculations were carried out for the corresponding reacting cases. The combustion models were then adjusted to fit the requirements of each flame. For all flame calculations good agreement of the main flow features with the measured data was achieved. For purely nonpremixed flames burning attached to the bluff-body’s outer edge, flamelet modeling including strain rate effects provided good results for the flow field and for most scalars. The prediction of a partially premixed swirl flame could only be achieved by applying a flamelet-based progress variable approach.     

Advective Transport in Heterogeneous Formations: The Impact of Spatial Anisotropy on the Breakthrough Curve

Water flow and solute transport take place in formations of spatially variable conductivity K. The logconductivity Y?= ln K is modeled as a random stationary space function, of normal univariate pdf (of mean In K _G and variance ${\sigma_{Y}^{2}}$ ) and of axisymmetric autocorrelation of integral scales I _h,I _v (anisotropy ratio f?=?I _v/I _h?<?1). The head gradient and the velocity are uniform in the mean, parallel to bedding, and of constant and given as J and U, respectively. Transport is ruled by advection, which typically overwhelms pore scale dispersion in the breakthrough curve (BTC) determination. In the present study we analyze the impact of anisotropy f on the BTC of a passive solute, which is related to the mass flux??? (t, x) at a control plane at x. While a considerable body of literature dealt with weakly heterogeneous formations ( ${\sigma _{Y}^{2} <1 }$ ), the present study addresses the case of ${\sigma _{Y}^{2} >1 }$ , which is of interest for many aquifers and is more difficult to solve either numerically or by approximations. We approach the three dimensional problem by modeling the structure as an ensemble of densely packed oblate spheroids of semi-major and semi-minor axis R and f R, respectively, and independent lognormal K, submerged in a matrix of uniform conductivity K _ef, the effective conductivity of the ensemble. The detailed numerical simulations of transport show that the BTC is insensitive to the value of the anisotropy ratio f, i.e.,??? (t, x) I _h/U depends only on ${\sigma _{Y}^{2}}$ (except for small differences in the tail). This important result implies that transport, as quantified by BTCs or spatial longitudinal mass distributions, can be modeled accurately by the much simpler solutions developed in the past for isotropic media, like e.g., the semi-analytical self-consistent approximation.     

Fourier Representation of Intermittent Flow in a Porous Medium

The field measurements and numerical results for intermittent flow regime in a sandy soil show that the time distributions of the soil water flux q(z, t), and the soil water content θ(z, t)at various depths are periodic in nature, where t is time and z is the depth (i.e., at the surface z = 0 and at depths z = − 5, − 10, − 15 cm, etc). The period of q(z, t) and θ(z, t) variations are generally determined by the sum of the duration of pulse and the duration between the initiation of two consecutive pulses of water at the soil surface. Fourier series models have been given for q(z, t) and θ(z, t) variations. The predicted Fourier results for these variations have been compared with the experimentally verified numerical results—designated as observed values. The results show that the amplitudes of these variations were damped exponentially with depth, and the phase shift increased linearly with depth.     

On the Modelling of the Electro-Hydrodynamic Flow Field in Electrostatic Precipitators

Electro-hydrodynamic (EHD) flows are investigated theoretically and numerically in this paper and results are presented for the flow field in model electrostatic precipitators (EPs). The resulting flow fields are shown in various representations and explained qualitatively. Numerical calculations with different flow models (non-turbulent and RANS) were conducted to investigate the influence of the flow model on the resulting secondary flows. Furthermore, a perturbation analysis is presented, leading to a simple differential equation of the Helmholtz type. This allows a more detailed view of the important mechanisms forming the secondary flows as well as being able to obtain a very fast estimation of the resulting flow field. The calculations reveal a strong influence of a vortex formation at the beginning of the precipitation zone on the whole flow field. Furthermore, a strong effect of the boundary conditions of the electric field and the operating parameters is shown. This revised version was published online in July 2006 with corrections to the Cover Date.     

Analysis of Algebraic Closures of the Mean Scalar Dissipation Rate of the Progress Variable Applied to Stagnating Turbulent Flames

In the turbulent premixed reactive flows considered in this study, i.e. large Damköhler and Reynolds numbers, the flamelet regime of turbulent combustion applies and the scalar dissipation rate and mean reaction rate are inter related. In this situation various algebraic models for the mean chemical rate that are obtained from an equilibrium of the dominant terms of the transport equation for the scalar dissipation rate, are evaluated through their application to flames stabilized in a turbulent stagnating flow. An asymptotic analysis is first performed and results obtained through the resulting one-dimensional calculation are compared with the experimental data of Li et al. (Proc Combust Inst 25:1207–1214, 1994). Eventually, three-dimensional CFD calculations including suited algebraic closures to represent the turbulent transport terms are carried out. Results are satisfactorily compared to the experimental data of Cho et al. (Proc Combust Inst 22:739–745, 1988). As a first outcome, the analysis confirms the interest and the relevance of the corresponding algebraic closures to deal with turbulent premixed combustion in such conditions. In the search of a satisfactory representation of such premixed impinging flames, the computational results also clearly emphasize the strong intertwinment that exits between the mean reaction rate, i.e. scalar dissipation rate or micro-mixing taking place at the smallest scale of the reactive flowfield, and the Reynolds fluxes modelling, i.e. turbulent macro-mixing.     

The Role of Probabilistic Approaches to Transport Theory in Heterogeneous Media

A physical picture of contaminant transport in highly heterogeneous porous media is presented. In any specific formation the associated governing transport equation is valid at any time and space scale. Furthermore, the advective and dispersive contributions are inextricably combined. The ensemble average of the basic transport equation is equivalent to a continuous time random walk (CTRW). The connection between the CTRW transport equation, in a limiting case and the familiar advection–dispersion equation (ADE) is derived. The CTRW theory is applied to the results of laboratory experiments, field observations, and simulations of random fracture networks. All of these results manifest dominant non-Gaussian features in the transport, over different scales, which are accounted for quantitatively by the theory. The key parameter controlling the entire shape of the contaminant plume evolution and breakthrough curves is advanced as a more useful characterization of the transport than the dispersion tensor, which is based on moments of the plume. The role of probabilistic approaches, such as CTRW, is appraised in the context of the interplay of spatial scales and levels of uncertainty. We then discuss a hybrid approach, which uses knowledge of non-stationary aspects of a field site on a larger spatial scale (trends) with a probabilistic treatment of unresolved structure on a smaller scale (residues).     

The Effect of Sweep-Angle Variation on the Turbulence Structure in a Separated,Three-Dimensional Flow

A three-dimensional separated flow behind a swept, backward-facing step is investigated by means of DNS for Re _H = C _∞ H/ν = 3000 with the purpose to identify changes in the statistical turbulence structure due to a variation of the sweep angle α from 0° up to 60°. With increasing sweep angle, the near-wall turbulence structure inside the separation bubble and downstream of reattachment changes due to the presence of an edge-parallel mean flow component W. Turbulence production due to the spanwise shear ∂W/∂y at the wall becomes significant and competes with the processes caused by impingement of the separated shear-layer. Changes due to a sweep angle variation can be interpreted in terms of two competing velocity scales which control the global budget of turbulent kinetic energy: the step-normal component U _∞ = C _∞cosα throughout the separated flow region and the velocity difference C _∞ across the entire shear-layer downstream of reattachment. As a consequence, the significance of history effects for the development into a two-dimensional boundary layer decreases with increasing sweep angle. For α ≥50°, near-wall streaks tend to form inside the separated flow region. Received 7 November 2000 and accepted 9 July 2002 Published online 3 December 2002 RID="*" ID="*" Part of this work was funded by the Deutsche Forschungsgemeinschaft within Sfb 557. Computer time was provided by the Konrad-Zuse Zentrum (ZIB), Berlin. Communicated by R.D. Moser