Numerical Simulation of Turbulent Pipe Flow Through an Abrupt Axisymmetric Constriction

Turbulent flow through an abrupt axisymmetric contraction has been investigated in a computational study. The full three-dimensional Navier-Stokes equations have been solved in cylindrical coordinates. The ring-like obstruction was accounted for by an immersed boundary method implemented in an open source code. The results were compared with data from an experimental study and a surprisingly good agreement was found in the separated-flow region downstream of the obstruction. A distinctly higher turbulence level was predicted upstream of the obstruction as compared with the measured data. Due to the imposed streamwise periodicity the simulated configuration mimicked a sequence of equally spaced obstructions and this periodicity gives rise to substantially higher turbulence levels than in a single-ring configuration. Particular attention was paid to the possible occurrence of mean-flow asymmetries, as frequently reported in planar-symmetric configuration. No sign of the Coanda effect was observed in the present study, for which a physical explanation was provided.     

Mean Flow and Turbulence Measurements in a Turbulent Free Cruciform Jet

The results of measurements of all three components of the mean velocity vector, the Reynolds normal and primary shear stresses and the mean static pressure in a turbulent free jet, issuing from a sharp-edged cruciform orifice, are presented in this paper. The measurements were made with an x-array hot-wire probe and a pitot-static tube in the near flow field of the jet. The Reynolds number, based upon the equivalent diameter of the orifice, was 1.70 × 10⁵. In addition to the quantities measured directly, the mean streamwise centreline velocity decay, the jet half-velocity widths, the jet spreading rate, the mean streamwise vorticity, the mass entrainment rate, the integral momentum flux and the one-dimensional energy spectra have been derived from the measured data. The results show that the mean streamwise centreline velocity decay rate of the cruciform jet is higher than that of a round jet issuing from an orifice with the same exit area as that of the cruciform orifice. The mean streamwise velocity field changed shape continuously from a cruciform close to the orifice exit plane to circular at 12 and half equivalent diameters downstream. The mean streamwise vorticity field, up to about three equivalent diameters downstream of the orifice exit plane, consists of four pairs of counter-rotating cells, which are aligned with the four edges in the centre of the cruciform orifice.     

A New Divergence Free Synthetic Eddy Method for the Reproduction of Inlet Flow Conditions for LES

This paper describes a recent development of the Synthetic Eddy Method (SEM) proposed by Jarrin et al. (Int J Heat Fluid Flow 30(3):435–442, 2009) for generation of synthetic turbulence. The present scheme is designed to produce a divergence-free turbulence field that can reproduce almost all possible states of Reynolds stress anisotropy. This improved representation, when used to provide inlet conditions for an LES, leads to reduced near-inlet pressure fluctuations in the LES and to a reduced development length, both of which lead to lower computer resource requirements. An advantage of this method with respect to forcing approaches (which require an iterative approach) is the suitability for direct usage with embedded LES. Results for a turbulent channel flow are reported here and compared to those from the original SEM, and other direct approaches such as the VORTEX method of Sergent (2002) and the Synthesized Turbulence approach of Davidson and Billson (Int J Heat Fluid Flow 27(6):1028–1042, 2006), showing overall improved performance and a more accurate representation of turbulence structures immediately downstream of the inlet.     

一种紊流润滑理论分析新方法—复合型紊流模式理论

在对主要用于高速轻载工况的常用紊理论进行简要分析和充分考察润滑流场的边界条件及内部结构的基础上，采用理论上比现有紊流润滑理论更为合理的复合型紊流模式理论，即在近壁区采用低紊流雷诺数的ｋ－ε模式，而在紊流核心区采用代数雷诺应力模式，对复杂流场的紊流润滑进行了分析，同时在计入惯性效应的情况下，推导出了一种适用于高压密封和高速重载轴承等计算，计算结果与实验数据十分吻合，验证了模型的有效性，可以应用于高压     

The Scale-Adaptive Simulation Method for Unsteady Turbulent Flow Predictions. Part 2: Application to Complex Flows

The paper summarises the validation activity performed with the Scale-Adaptive Simulation turbulence model (SAS model) using the two commercial CFD solvers, ANSYS-FLUENT and ANSYS-CFX. Both the KSKL-SAS and the SST-SAS model variants have been tested, although most cases have been computed with the second. The turbulence-resolving capability of the SAS method has been validated with a representative set of test cases, covering both underlying generic flows as well as practical engineering applications. In addition to the purely aerodynamic flows with massive separation and heat transfer they include also such physical phenomena as turbulent combustion and aeroacoustics. The illustrating results show the potentials of the SAS approach for industrial flow simulations. Most of the test case simulations were conducted during the recent EU project “DESider”.     

A Very Large Eddy Simulation Model Using a Reductionist Inlet Turbulence Generator and Wall Modeling for Stable Atmospheric Boundary Layers

Fluid Dynamics - Despite many advances in numerical simulation of stable boundary layers (SBL), most of the models developed are complex and computationally expensive. A computational fluid...     

Simulation and Modelling of a Skewed Turbulent Channel Flow

A time-dependent three-dimensionally skewed flow is investigated using direct numerical simulations of the incompressible Navier-Stokes equations. The effect on the instantaneous and mean turbulent field is investigated. Instantaneous flowfields reveal that the skewing has the effect of initially reducing the strength and height of quasi-streamwise vortices of both signs of rotation with respect to the skewing. A mechanism for this process is put forward. The mean flowfields show drops in turbulence quantities such as turbulence kinetic energy. In addition to this, two-equation turbulence modelling of the flow is carried out. This highlights a deficiency, in that the standard turbulence models are unable to capture the drop in turbulence intensity due to the skewing. A modification based on the exact dissipation equation is found to significantly improve the model behaviour for this flow. This revised version was published online in July 2006 with corrections to the Cover Date.     

黏弹性流体扩展及收缩流动的格子Boltzmann模拟分析

本文应用格子Boltzmann 方法(LBM)并结合Oldroyd-B 模型,讨论了不可压缩的 Navier-Stokes 方程和平流扩散本构方程的解耦及各自求解方法,以及两类问题的边界处理格式,实现了黏弹性流体在二维1:3 扩展流道以及3:1 收缩流道中的流动的数值模拟．获得了不同雷诺数Re 和维森伯格数Wi 以及黏度vs 下流动的流线分布,计算给出了漩涡的涡心位置和大小,并分析了参数Re、Wi 和vs 对流动特点的影响．模拟结果表明本文所采用模型和边界处理方法具有良好的精度和稳定性．     

Numerical Simulation of Turbulent Flow in Concentric Annuli

In this paper we consider a fully developed turbulent flow in a round pipe with a small inner annulus. The diameter of the inner annulus is less than 10% of the diameter of the outer pipe. As a consequence, the surface area of the inner pipe compared to the outer pipe is small. The friction exerted by the wall on the flow is proportional to the surface area and the wall shear stress. Due to the small surface area of the inner annulus the additional stress on the flow due to the presence of the annulus may expected to be negligible. However, it will be shown that the inner annulus drastically changes the flow patterns and gives rise to unexpected scaling properties. In previous studies (Chung et al., Int J Heat Fluid Flow 23:426–440, 2002; Churchill and Chan, AIChE J 41:2513–2521, 1995) it was argued that radial position of the point of zero shear stress does not coincide with the radial location of the point of maximum axial velocity. In our direct numerical simulations we observe a coincidence of these points within the numerical accuracy of our model. It is shown that the velocity profile close to the inner annulus is logarithmic.     

Turbulent Transition in the Inlet Region of a Circular Pipe

Calculated and experimental data on turbulent transition in a circular pipe are analyzed. The calculations were performed using the three-parameter turbulence model. The dependence of the distance from the inlet to the point of minimum friction during transition on the Reynolds number for fixed inlet conditions and the distribution of the turbulence parameters over the pipe length and radius are obtained. The dependence of the maximum (critical) Reynolds number, Re^*, for which there is no transition in the pipe, on the inlet intensity and scale of turbulence is found. It is suggested that Re^* depends on the inlet perturbation parameters up to Re^* = 1000, where the friction coefficients for laminar and turbulent flows coincide.     

Turbulence Modelling of Unsteady Turbulent Flows Using the Stress Strain Lag Model

This paper reports the application of a recently developed turbulence modelling scheme known as the C _as model. This model was specifically developed to capture the effects of stress-strain misalignment observed in turbulent flows with mean unsteadiness. Earlier work has reported the approach applied within a linear k-ε modelling framework, and some initial testing of it within the k-ω SST model of Menter (AIAA J 32:1598–1605, 1994). The resulting k-ε-C _as and SST-C _as models have been shown to result in some of the advantages of a full Reynolds Stress transport Model (RSM), whilst retaining the computational efficiency and stability benefits of a eddy viscosity model (EVM). Here, the development of the the high-Reynolds-number version of the C _as model is outlined, with some example applications to steady and unsteady homogeneous shear flows. The SST-C _as form of the model is then applied to further, more challenging cases of 2-D flow around a NACA0012 aerofoil beyond stall and the 3-D flow around a circular cylinder in a square duct, both being flows which exhibit large, unsteady, separated flow regions. The predictions returned by a range of other common turbulence modelling schemes are included for comparison and the SST-C _as scheme is shown to return generally good results, comparable in some respects to those obtainable from far more complex schemes, for only moderate computing resource requirements.     

Large Eddy Simulation of Compressible Subsonic Turbulent Jet Starting From a Smooth Contraction Nozzle

Compressible subsonic turbulent starting jet with a relatively large Reynolds number of significant practical importance is investigated using large eddy simulation (LES), starting from a smooth contraction nozzle. The computational domain of truncated conical shape is determined through the comparison of the time-averaged numerical solution with the particle imaging velocimetry measurements for the steady jet. It is shown that the starting jet consists of a leading vortex ring followed by a quasi-steady jet, and the instantaneous velocity field exhibits contraction and expansion zones, corresponding to the high pressure (HP) and low pressure (LP) regions formed by the convecting vortex rings, and are related to the Kelvin-Helmholtz instability. The thin boundary layer inside the smooth contraction nozzle evolves into a shear layer at the nozzle exit and develops with the downstream penetration of the jet. Using λ ₂ criterion, the formation and evolution of the vortical structures are temporally visualized, illustrating distortion of vortex rings into lobed shapes prior to break-down. Rib-shape streamwise vortex filaments exist in the braid region between a pair of consecutive vortex rings due to secondary instabilities. Finally, formation and dynamics of hairpin vortices in the shear layer is identified.     

Evaluation of the Filtered Noise Turbulent Inflow Generation Method

Flows around buildings and in urban areas have the ability to exchange mass and momentum through mixing layers. The complex dynamical phenomena arising in mixing layers can be studied using Large Eddy Simulation (LES). As mixing layers depend on the turbulence conditions upstream of the buildings or urban areas, appropriate turbulent inlet conditions have to be provided to a simulation. Due to the high efficiency and level of control, the filtered noise inflow method was selected. The control over the Reynolds stresses as well as nine length scales make this method suitable to replicate conditions measured in experiments. In this paper, a formal method to obtain the filter coefficients is presented. This is achieved by relating the spatial filtering to a Finite Impulse Response (FIR) filter and the temporal filtering to an Autoregressive (AR) model. Three closed-form solutions for the spatial filter coefficients are presented having a Gaussian, double-exponential and exponential correlation function. By means of an LES simulation of a turbulent wall-bounded flow, the input-output behaviour is investigated. It was found that a combination of a Gaussian filter with length scales that increase with increasing wall distance result in the fastest downstream development of the artificial turbulence and the smallest loss of turbulent kinetic energy.     

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.     

Large-Eddy Simulation of Dispersion from Line Sources in a Turbulent Channel Flow

Large-eddy simulations of the dispersion from scalar line sources at various locations within a fully developed turbulent channel flow at Re = uh/ν = 10,400 are presented. Both mean and fluctuating scalar quantities are compared with those from the single available set of experimental data (Lavertu and Mydlarski, J Fluid Mech 528:135–172, 2005) and differences are highlighted and discussed. The results are also discussed in the context of scalar dispersion in other kinds of turbulent flows, e.g. homogeneous shear-flow. Initial computations at a much lower Reynolds number are also reported and compared with the two available direct numerical simulation data sets.     

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.     

Direct Numerical Simulation of Turbulent Flow over a Rectangular Trailing Edge

This paper describes a direct numerical simulation (DNS) study of turbulent flow over a rectangular trailing edge at a Reynolds number of 1000, based on the freestream quantities and the trailing edge thickness h; the incoming boundary layer displacement thickness δ^* is approximately equal to h. The time-dependent inflow boundary condition is provided by a separate turbulent boundary layer simulation which is in good agreement with existing computational and experimental data. The turbulent trailing edge flow simulation is carried out using a parallel multi-block code based on finite difference methods and using a multi-grid Poisson solver. The turbulent flow in the near-wake region of the trailing edge has been studied first for the effects of domain size and grid resolution. Then two simulations with a total of 256 × 512 × 64 (∼ 8.4×10⁶) and 512 × 1024 × 128 (∼ 6.7×10⁷) grid points in the computational domain are carried out to investigate the key flow features. Visualization of the instantaneous flow field is used to investigate the complex fluid dynamics taking place in the near-wake region; of particular importance is the interaction between the large-scale spanwise, or Kármán, vortices and the small-scale quasi-streamwise vortices contained within the inflow boundary layer. Comparisons of turbulence statistics including the mean flow quantities are presented, as well as the pressure distributions over the trailing edge. A spectral analysis applied to the force coefficient in the wall normal direction shows that the main shedding frequency is characterized by a Strouhal number based on h of approximately 0.118. Finally, the turbulence kinetic energy budget is analysed. Received 4 March 1999 and accepted 27 October 2000     

Eddy-Viscosity and Stress-Transport Turbulence Models in Application to a Plane Synthetic Jet

A two-dimensional unsteady plane jet is computed using four turbulence models. The strain field generated by the ‘synthetic’ (reciprocating) jet under consideration contains regions of strong normal and shear straining, while vortex pairs formed by the issuing jet are highly rotational. Turbulence closures at the eddy-viscosity level (standard k???ε and two-timescale models) and two stress-transport closures (the ‘basic’ model and a two-component-limit scheme) are examined against experimental data. It is found that the extreme nature of the strain environment presents a severe challenge to the models tested; in particular, fundamental weaknesses in the constitutive equation of linear eddy viscosity models are highlighted. An additional term in the dissipation rate equation of the two-timescale model is found to hinder the formation of clearly defined starting vortex pairs which play a critical role in determining the near-field pattern of synthetic jet development. The stress-transport models are able to reproduce measured turbulence quantities with a greater degree of accuracy, and the two-component-limit formulation shows an improvement over the basic model.     

On the Model of Generation of Vortex Structures in an Isotropic Turbulent Flow

It is known that turbulence is characterized by intermittence which is closely related to the development of unsteady nonisotropic intense small-scale vortex structures. In this study, small fluid particles from the inertial range of isotropic turbulence are considered. It is shown that the phenomenon of rotation intensification and stretching of the particles can be analyzed theoretically. In recent experimental and numerical studies, where this phenomenon was called “the pirouette effect”, its significance in the mechanism of the intense small-scale structures generation was discussed. In this study, a linear stochastic Lagrangian model for the effect is developed. In this model, the kinetic equation for the distribution function of the squared cosine of the angle between the vorticity and the eigenvector of the strain rate tensor of a fluid particle is derived and time history asymptotics of this quantity are analytically calculated at large and small times. The results are in good agreement with the recent experiments and numerical calculations. An analysis made in this study shows that the linear processes probably play the crucial role in certain processes in the isotropic turbulence, which is known to be a principally nonlinear phenomenon. The model developed makes it possible to analyze the statistics of the Lagrangian dynamics of small fluid particles in the inertial range which can be useful in some computational approaches to turbulence.