Computation of axisymmetric confined jets in a diffuser

The paper reports on a numerical study of turbulent confined jets in a conical duct with a 5° divergence. The flow has a large ratio of jet to ambient velocities at the entrance so that it gives rise to strong recirculation. The calculations are carried out with a general finite volume method designed for calculating incompressible elliptic flows with complex boundaries. Turbulence is simulated by the standard κ–? model. The sensitivity of the solution to numerical discretization errors is examined using three convection schemes, i.e. hybrid central/upwind differencing, QUICK and SOUCUP, on two grids consisting of 68 × 50 and 102 × 82 points respectively. An examination is also made of the influence of inlet boundary conditions on the predicted flow field. The computed results are compared with experimental data for mean axial velocity, turbulent shear stress and turbulent kinetic energy profiles. It is shown that the calculations reproduce the essential features of the flow observed in the experiments.     

An experimental investigation on flow structures of confined and unconfined impinging air jets

The flow characteristics of both confined and unconfined air jets, impinging normally onto a flat plate have been experimentally investigated. The mean and turbulence velocities, and surface pressures were measured for Reynolds numbers ranging from 30,000 to 50,000 and the nozzle-to-plate spacings in range of 0.2–6. Smoke-wire technique is used to visualize the flow behavior. The effects of Reynolds number, nozzle-to-plate spacing and flow confinement on the flow structure are reported. In the case of confined jet, subatmospheric regions occur on both impingement and confinement surfaces at nozzle-to-plate spacings up to 2 for all Reynolds numbers in consideration and they lie up to nearly the same radial location at both surfaces. However, there is no evidence of the subatmospheric region in unconfined jet. It is concluded that there exists a linkage among the subatmospheric region, turbulence intensity and the peaks in heat transfer coefficients for low spacings in impinging jets.     

Computation of strongly swirling confined flows with cubic eddy‐viscosity turbulence models

An investigation on the predictive performance of four cubic eddy‐viscosity turbulence models for two strongly swirling confined flows is presented. Comparisons of the prediction with the experiments show clearly the superiority of cubic models over the linear k–εmodel. The linear k–εmodel does not contain any mechanism to describe the stabilizing effects of swirling motion and as a consequence it performs poorly. Cubic models return a lower level of Reynolds stresses and the combined forced‐free vortex profiles of tangential velocity close to the measurements in response to the interaction between swirl‐induced curvature and stresses. However, a fully developed rotating pipe flow is too simple to contain enough flow physics, so the calibration of cubic terms is still a topic of investigation. It is shown that explicit algebraic stress models require fewer calibrations and contain more flow physics. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical predictions of the bifurcation of confined swirling flows

The bifurcation of confined swirling flows was numerically investigated by employing both the k-? and algebraic stress turbulence models. Depending upon the branch solution examined, dual flow patterns were predicted at certain swirl levels. In the lower-branch solution which is obtained by gradually increasing the swirl level from a low-swirl flow, the flow changes with increasing swirl number from the low-swirl flow pattern to a high-swirl flow pattern. In the upper-branch solution which is acquired by gradually decreasing the swirl level from a high-swirl flow, on the other hand, the flow can maintain itself in the high-swirl flow pattern at the swirl levels where it exhibits the low-swirl flow pattern in the lower branch. The bifurcation of confined swirling flows was predicted with either the k-? model or the algebraic stress model being employed. Both the k-? and algebraic stress models result in comparable and sufficiently good predictions for confined swirling flows if high-order numerical schemes are used. The reported poor performance of the k-? model was clarified to be mainly attributable to the occurrence of the bifurcation and the use of low-order numerical schemes.     

Wavepacket models for subsonic twin jets using 3D parabolized stability equations

An extension of the classical parabolized stability equations to flows strongly dependent on the two cross-stream spatial directions and weakly dependent on the streamwise one is applied to model the large-scale structures present in twin-jet configurations. The existence of these unsteady flow structures, usually referred to as wavepackets, has been demonstrated in the literature for both subsonic and supersonic round jets, along with their relation to the generation of highly directional noise emitted in the aft direction. The present study considers twin-jet configurations with different separations at high Reynolds number and subsonic conditions. The existing instability modes for the twin-jet mean flow, their dependence on the separation of the two jets, and the interaction between the wavepackets originating from the two jets is investigated here. Arising from the axisymmetric mode for single round jets, two dominant modes are found for twin jets: a varicose one, relatively insensitive to jets' proximity, but likely to be efficient in radiating noise; a sinuous one, whose amplification is strongly dependent on the jets' distance, and which can be expected to produce weaker acoustic signatures.     

From single-scale turbulence models to multiple-scale and subgrid-scale models by Fourier transform

A theoretical method based on mathematical physics formalism that allows transposition of turbulence modeling methods from URANS (unsteady Reynolds averaged Navier–Stokes) models, to multiple-scale models and large eddy simulations (LES) is presented. The method is based on the spectral Fourier transform of the dynamic equation of the two-point fluctuating velocity correlations with an extension to the case of non-homogenous turbulence. The resulting equation describes the evolution of the spectral velocity correlation tensor in wave vector space. Then, we show that the full wave number integration of the spectral equation allows one to recover usual one-point statistical closure whereas the partial integration based on spectrum splitting gives rise to partial integrated transport models (PITM). This latter approach, depending on the type of spectral partitioning used, can yield either a statistical multiple-scale model or subfilter transport models used in LES or hybrid methods, providing some appropriate approximations are made. Closure hypotheses underlying these models are then discussed by reference to physical considerations with emphasis on identification of tensorial fluxes that represent turbulent energy transfer or dissipation. Some experiments such as the homogeneous axisymmetric contraction, the decay of isotropic turbulence, the pulsed turbulent channel flow and a wall injection induced flow are then considered as typical possible applications for illustrating the potentials of these models.      

Apparent transition behavior of widely-used turbulence models

The Spalart–Allmaras and the Menter k–ω SST turbulence models are shown to have the undesirable characteristic that, for fully turbulent external flow computations, a transition region can occur whose extent varies with grid density. Extremely fine two-dimensional grids over the front portion of an airfoil are used to demonstrate the effect. As the grid density is increased, the laminar region near the nose becomes larger. In the Spalart–Allmaras model this behavior is due to convergence to a laminar-behavior fixed point that occurs in practice when freestream turbulence is below some threshold. It is the result of a feature purposefully added to the original model in conjunction with a special trip function. This degenerate fixed point can also cause non-uniqueness regarding where transition initiates on a given grid. Consistent fully turbulent results can easily be achieved by either using a freestream turbulence level higher than the threshold or by making a simple change to one of the model constants. Near the area where turbulence initiates, the SST model exhibits sensitivity to numerical resolution, but its solutions are unique on a given grid. Inconsistent apparent transition behavior with grid refinement in this case does not stem from the presence of a degenerate fixed point. A nullcline analysis is used to visualize the local behavior of the model.     

内锥流量计流出系数预测方法研究

采用标准k-ε模型、RNG(Renormalization Group)k-ε模型、Realizable k-ε模型和Reynolds应力方程模型 RSM(Reynolds Stress Model) 对100 mm口径6种结构的内锥流量计内流场进行了数值模拟.在等效直径比β值为0.65的三种结构内锥流量计流出系数的仿真计算中,四种湍流模型计算结果与物理实验结果误差的平均值分别为4.19%,2.84%,2.88%和-0.822%;对β值为0.85的情况,各模型计算误差的平均值分别为11.8%,9.62%,9.30%和4.76%.研究结果表明,RSM模型在6种结构内锥流量计流出系数的预测中,计算精度较高,表现出了较好的性能,优于三种k-ε涡粘模型,更适于内锥流量计流场数值模拟与流出系数的预测.     

Particle modulations to turbulence in two-phase round jets

The particle modulations to turbulence in round jets were experimentally studied by means of two-phase velocity measurements with Phase Doppler Anemometer （PDA）. Laden with very large particles, no significant attenuations of turbulence intensities were measured in the farfields, due to small two-phase slip velocities and particle Reynolds number. The gas-phase turbulence is enhanced by particles in the near-fields, but it is significantly attenuated by the small particles in the far-fields. The smaller particles have a more profound effect on the attenuation of turbulence intensities. The enhancements or attenuations of turbulence intensities in the far-fields depends on the energy production, transport and dissipation mechanisms between the two phases, which are determined by the particle prop- erties and two-phase velocity slips. The non-dimensional parameter CTI is introduced to represent the change of turbulence intensity.     

Assessment of the performance of different classes of turbulence models in a wide range of non-equilibrium flows

The performance of thirteen benchmark turbulence models within the RANS framework has been assessed in classical non-equilibrium flows. Linear and non-linear eddy-viscosity schemes, Reynolds stress transport models and single- and two-time-scale approaches have been considered in the investigation. Among the test cases studied are homogeneous shear and normally strained flows, adverse-pressure-gradient, favourable-pressure-gradient and oscillatory boundary layer flows, fully developed oscillatory and ramp up pipe flows and steady and pulsated backward-facing-step flows. The main advantages and drawbacks of the models in each of the test cases are discussed. These discussions provide a reasonably wide understanding of the expected behaviour of the models for future applications in non-equilibrium flows, and also result in suggestions on how the effectiveness of existing models can be further improved.     

Cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density

An investigation on the predictive performance of cubic eddy‐viscosity turbulence models for strongly swirling confined flows with variable density is presented. Comparisons of the prediction with the experiments show some improvements of cubic models over the linear k–ε model. The linear k–ε model does not contain any mechanism to represent the interaction of swirl and density variation and as a consequence it performs poorly. With appropriate modelling, two‐equation cubic turbulence models can capture the subcritical nature of the flow, represent the azimuthal velocity profiles of combined forced‐free vortex motion, and predict the combined effects of swirl and density variation fairly well. However, the calibration of model coefficients is still a topic of investigation. Further amendments are also needed for the equations of k and ε to take into account the effects of swirl and density gradients correctly. Copyright © 2004 John Wiley & Sons, Ltd.     

Large eddy simulation of compressible round jets with coflow

Large eddy simulations of subsonic round jets are carried out using high order compact finite difference scheme and an explicit filtering based approximate deconvolution method. The jets have a Mach number of 0.9 and Reynolds number of $4.5\times {10}^5$ based on jet diameter and centerline velocity at inflow. Results obtained for the mean flow and turbulence intensities agree well with those in existing literature. We also study the effects of co-flow velocity ratio on the flow physics. Increase in potential core length and decrease in spreading rate of jet is observed in the presence of co-flow. The effects of co-flow velocity ratio on the axial Reynolds stress and turbulent kinetic energy budgets are also presented. It is observed that increasing co-flow velocity ratio leads to reduction in the turbulence intensities and near-field sound levels.     

Two-equation turbulence models for prediction of heat transfer on a transonic turbine blade

Two versions of the two-equation k–ω model and a shear stress transport (SST) model are used in a three-dimensional, multi-block, Navier–Stokes code to compare the detailed heat transfer measurements on a transonic turbine blade. It is found that the SST model resolves the passage vortex better on the suction side of the blade, thus yielding a better comparison with the experimental data than either of the k–ω models. However, the comparison is still deficient on the suction side of the blade. Use of the SST model does require the computation of distance from a wall, which for a multi-block grid, such as in the present case, can be complicated. However, a relatively easy fix for this problem was devised. Also addressed are issues such as (1) computation of the production term in the turbulence equations for aerodynamic applications, and (2) the relation between the computational and experimental values for the turbulence length scale, and its influence on the passage vortex on the suction side of the turbine blade.     

A structural model for the vortex tubes of isotropic turbulence

A class of exact solutions of the Navier–Stokes equations is introduced to model the fine-scale, tubular structures of isotropic turbulence. The model vortices exhibit slow algebraic fall-off of the induced velocity, and accurately reproduce the velocity signatures observed in DNS and experiments. The proposed model has interesting implications for the theoretical analysis of turbulence, supporting the view that the inertial range energy scaling may have a link with the near-singular velocity field induced by vortex tubes produced by the roll-up of vortex sheets.      

Hybrid RANS/LES computations of plane impinging jets with DES and PANS models

The qualities of a DES (Detached Eddy Simulation) and a PANS (Partially-Averaged Navier–Stokes) hybrid RANS/LES model, both based on the k–ω RANS turbulence model of Wilcox (2008, “Formulation of the k–ω turbulence model revisited” AIAA J., 46: 2823–2838), are analysed for simulation of plane impinging jets at a high nozzle-plate distance (H/B = 10, Re = 13,500; H is nozzle-plate distance, B is slot width; Reynolds number based on slot width and maximum velocity at nozzle exit) and a low nozzle-plate distance (H/B = 4, Re = 20,000). The mean velocity field, fluctuating velocity components, Reynolds stresses and skin friction at the impingement plate are compared with experimental data and LES (Large Eddy Simulation) results. The k–ω DES model is a double substitution type, following Davidson and Peng (2003, “Hybrid LES–RANS modelling: a one-equation SGS model combined with a k–ω model for predicting recirculating flows” Int. J. Numer. Meth. Fluids, 43: 1003–1018). This means that the turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the eddy viscosity formula. The k–ω PANS model is derived following Girimaji (2006, “Partially-Averaged Navier–Stokes model for turbulence: a Reynolds-Averaged Navier–Stokes to Direct Numerical Simulation bridging method” J. Appl. Mech., 73: 413–421). The turbulent length scale in the PANS model is constructed from the total turbulent kinetic energy and the sub-filter dissipation rate. Both hybrid models change between RANS (Reynolds-Averaged Navier–Stokes) and LES based on the cube root of the cell volume. The hybrid techniques, in contrast to RANS, are able to reproduce the turbulent flow dynamics in the shear layers of the impacting jet. The change from RANS to LES is much slower however for the PANS model than for the DES model on fine enough grids. This delays the break-up process of the vortices generated in the shear layers with as a consequence that the DES model produces better results than the PANS model.     

Algebraic turbulence models for the computation of two-dimensional high-speed flows using unstructured grids

The incorporation of algebraic turbulence models in a solver for the 2D compressible Navier–Stokes equations using triangular grids is described. A practical way to use the Cebeci–Smith model and to modify it in separated regions is proposed. The ability of the model to predict high-speed perfect-gas boundary layers is investigated from a numerical point of view.     

Simulation of sonic boom interaction with shear turbulence

We investigate the interaction of pressure waves with the shear turbulence in a supersonic jet flow obtained from the direct numerical integration of the time-dependent, three-dimensional Euler equations. The resulting self-sustaining turbulent jet flow matches closely the relevant time and length scales of atmospheric turbulence. N-waves, characteristic of sonic booms, are simulated by perturbing the pressure profile and following these perturbations in space and time. The results reproduce most of the wave forms reported in laboratory experiments and in flight test data.     

The low Reynolds number turbulent flow and mixing in a confined impinging jet reactor

Turbulent mixing takes an important role in chemical engineering, especially when the chemical reaction is fast compared to the mixing time. In this context a detailed knowledge of the flow field, the distribution of turbulent kinetic energy (TKE) and its dissipation rate is important, as these quantities are used for many mixing models. For this reason we conduct a direct numerical simulation (DNS) of a confined impinging jet reactor (CIJR) at Re = 500 and Sc = 1. The data is compared with particle image velocimetry (PIV) measurements and the basic flow features match between simulation and experiment. The DNS data is analysed and it is shown that the flow is dominated by a stable vortex in the main mixing duct. High intensities of turbulent kinetic energy and dissipation are found in the impingement zone which decrease rapidly towards the exit of the CIJR. In the whole CIJR the turbulence is not in equilibrium. The strong mixing in the impingement zone leads to a rapid development of a monomodal PDF. Due to the special properties of the flow field, a bimodal PDF is generated in cross-sections downstream the impingement zone, that slowly relaxes under relaminarising conditions. The time required for meso-mixing is dominating the overall mixing performance.