Explicit algebraic Reynolds stress model for shock-dominated flows

Shock waves drastically alter the nature of Reynolds stresses in a turbulent flow, and conventional turbulence models cannot reproduce this effect. In the present study, we employ explicit algebraic Reynolds stress model (EARSM) to predict the Reynolds stress anisotropy generated by a shockwave. The model by Wallin and Johansson (2000) is used as the baseline model. It is found to over-predict the post-shock Reynolds stresses in canonical shock turbulence interaction. The budget of the transport equation of Reynolds stresses computed using linear interaction analysis shows that the unsteady shock distortion mechanism and the pressure–velocity correlations are important. We propose improvement to the baseline model using linear interaction analysis results and redistribute the turbulent kinetic energy between the principle Reynolds stresses. The new model matches DNS data for the amplification of Reynolds stresses across the shock and their post-shock evolution, for a range of Mach numbers. It is applied to oblique shock/boundary-layer interaction at Mach 5. Significant improvements are observed in predicting surface pressure and skin friction coefficient, with respect to experimental measurements.     

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

雷诺应力对平均和脉动压力影响的数值模拟

采用增加雷诺应力边界条件方法考察了雷诺应力对流场的影响。当雷诺应力边界条件存在时,流场内平均和脉动压力的变化体现了雷诺应力的影响。本文首先在没有雷诺应力边界条件和有雷诺应力边界条件的情况下得到TTU低矮建筑标准模型表面的平均压力系数,结果表明:考虑雷诺应力可使建筑物表面的平均压力系数的绝对值增大。本文还利用Selvam提出的脉动压力系数数学模型,计算了TTU模型表面的脉动压力系数,并对计算结果进行了分析。     

Hydrodynamic and thermal fields analysis in gas-solid two-phase flow

The present work aims to investigate numerically the flowfield and heat transfer process in gas-solid suspension in a vertical pneumatic conveying pipe. The Eulerian-Lagrangian model is used to simulate the flow of the two-phases. The gas phase is simulated based on Reynolds Average Navier-Stokes equations (RANS) with low Reynolds number k-ε model, while particle tracking procedure is used for the solid phase. An anisotropic model is used to calculate the Reynolds stresses and the turbulent Prandtl number is calculated as a function of the turbulent viscosity. The model takes into account the lift and drag forces and the effect of particle rotation as well as the particles dispersion by turbulence effect. The effects of inter-particles collisions and turbulence modulation by the solid particles, i.e. four-way coupling, are also included in the model. Comparisons between different models for turbulence modulation with experimental data are carried out to select the best model. The model is validated against published experimental data for velocities of the two phases, turbulence intensity, solids concentration, pressure drop, heat transfer rates and Nusselt number distribution. The comparisons indicate that the present model is able to predict the complex interaction between the two phases in non-isothermal gas-solid flow in the tested range. The results indicate that the particle-particle collision, turbulence dispersion and lift force play a key role in the concentration distribution. In addition, the heat transfer rate increases as the mass loading ratio increases and Nusselt number increases as the pipe diameter increases.     

BOUNDARY TREATMENT AND AN EFFICIENT PRESSURE ALGORITHM FOR INTERNAL TURBULENT FLOWS USING THE PDF METHOD

The generalized Langevin model, which is used to model the motion of stochastic particles in the velocity–composition joint probability density function (PDF) method for reacting turbulent flows, has been extended to incorporate solid wall effects. Anisotropy of Reynolds stresses in the near-wall region has been addressed. Numerical experiments have been performed to demonstrate that the forces in the near-wall region of a turbulent flow cause the stochastic particles approachi ng a solid wall to reverse their direction of motion normal to the wall and thereby, leave the near-wall layer. This new boundary treatment has subsequently been implemented in a full-scale problem to prove its validity. The test problem considered here is that of an isothermal, non-reacting turbulent flow in a two-dimensional channel with plug inflow and a fixed back-pressure. An efficient pressure correction method, developed in the spirit of the PISO algorithm, has been implemented. The pressure correction strategy is easy to implement and is completely consistent with the time- marching scheme used for the solution of the Lagrangian momentum equations. The results show remarkable agreement with both k–ϵ and algebraic Reynolds stress model calculations for the primary velocity. The secondary flow velocity and the turbulent moments are in better agreement with the algebraic Reynolds stress model predictions than the k– ϵ predictions. © 1997 by John Wiley & Sons, Ltd.     

Study on anisotropic buoyant turbulence model

Ｂｕｏｙａｎｔｆｌｏｗｉｓｏｎｅｏｆｔｈｅｆｕｎｄａｍｅｎｔａｌｆｌｏｗｓ．Ｔｈｅｄｉｆｆｅｒｅｎｃｅｏｆｄｅｎｓｉｔｙｂｅｔｗｅｅｎｄｉｓｃｈａｒｇｅｄｆｌｕｉｄａｎｄｔｈｅａｍｂｉｅｎｔｆｌｕｉｄｃａｎｃａｕｓｅｂｕｏｙａｎｔｔｕｒｂｕｌｅｎｔｆ...     

Numerical analysis of developing turbulent flow in a 180° bend tube by an algebraic Reynolds stress model

A numerical analysis has been performed for three‐dimensional developing turbulent flow in a 180° bend tube with straight inlet and outlet section used by an algebraic Reynolds stress model. To our knowledge, numerical investigations, which show the detailed comparison between calculated results and experimental data including distributions of Reynolds stresses, are few and far between. From this point of view, an algebraic Reynolds stress model in conjunction with boundary‐fitted co‐ordinate system is applied to a 180° bend tube in order to predict the anisotropic turbulent structure precisely. Calculated results are compared with the experimental data including distributions of Reynolds stresses. As a result of this analysis, it has been found that the calculated results show a comparatively good agreement with the experimental data of the time‐averaged velocity and the secondary vectors in both the bent tube and straight outlet sections. For example, the location of the maximum streamwise velocity, which appears near the top or bottom wall in the bent tube, is predicted correctly by the present method. As for the comparison of Reynolds stresses, the present method has been found to simulate many characteristic features of streamwise normal stress and shear stresses in the bent tube qualitatively and has a tendency to under‐predict its value quantitatively. Judging from the comparison between the calculated and the experimental results, the algebraic Reynolds stress model is applicable to the developing turbulent flow in a bent tube that is known as a flow with a strong convective effect. Copyright © 2005 John Wiley & Sons, Ltd.     

Particle-turbulence interaction in a boundary layer

Particle-turbulence interaction in wall turbulent flows has been studied. A series of experiments varying particle size, particle density, particle loading and flow Re has been conducted. The results show that the larger polystyrene particles (1100 μm) cause an increase in the number of wall ejections, giving rise to an increase in the measured values of the turbulence intensities and Reynolds stresses. On the other hand, the smaller polystyrene particles (120 μm) bring about a decrease in the number of wall ejections, causing a decrease in the measured intensities and Reynolds stresses. These effects are enhanced as the particle loading is increased. It was also found that the heavier glass particles (88 μm) do not bring about any significant modulation of turbulence. In addition, measurements of the burst frequency and the mean streak-spacing show no significant change with increase in particle loading. Based on these observations, a mechanism of particle transport in wall turbulent flows has been proposed, in which the particles are transported (depending on their size, density and flow Re) by the bursting events of the wall regions.     

Dynamics of dual-particles settling under gravity

As a first step towards understanding particle–particle interaction in fluid flows, the motion of two spherical particles settling in close proximity under gravity in Newtonian fluids was investigated experimentally for particle Reynolds numbers ranging from 0.01 to 2000. It was observed that particles repel each other for Re>0.1 and that the separation distance of settling particles is Reynolds number dependent. At lower Reynolds numbers, i.e. for Re<0.1, particles settling under gravity do not separate.The orientation preference of two spherical particles was found to be Reynolds number dependent. At higher Reynolds numbers, the line connecting the centres of the two particles is always horizontal, regardless of the way the two particles are launched. At lower Reynolds numbers, however, the particle centreline tends to tilt to an arbitrary angle, even of the two particles are launched in the horizontal plane. Because of the tilt, a side migration of the two particles was found to exist. A linear theory was developed to estimate the side migration velocity. It was found that the maximum side migration velocity is approximately 6% of the vertical settling velocity, in good agreement with the experimental results.Counter-rotating spinning of the two particles was observed and measured in the range of Re=0–10. Using the linear model, it is possible to estimate the influence of the tilt angle on the rate of rotation at low Reynolds numbers. Dual particles settle faster than a single particle at small Reynolds numbers but not at higher Reynolds numbers, because of particle separation. The variation of particle settling velocity with Reynolds number is presented. An equation which can be used to estimate the influence of tilt angle on particle settling velocity at low Reynolds number is also derived.     

Development of a Two-velocities Hybrid RANS-LES Model and its Application to a Trailing Edge Flow

The flow around a trailing edge is computed with a new hybrid method designed to more clearly separate the effects of total and sub-grid turbulent stress-modelling on the time-averaged and instantaneous velocity fields, and in turn, mean momentum and kinetic energy balances. These two velocity fields independently define Reynolds averaged and sub-grid-scale viscosities, and distinct stresses, at the same location. In particular, resolved eddies can emerge, or sweep in and out of the Reynolds averaged near wall layer, without being dampened by higher levels of the viscosity in this RANS dominated layer. The two-field hybrid model, first tested on channel flows, gives accurate predictions of mean velocities and stresses for different Reynolds numbers and coarse meshes. For the trailing edge flow the results of the hybrid model are close to the reference fine LES for mean velocity and turbulent content, whereas the DES-SST on the same coarse mesh gives too early separation.     

Riblet modelling using a second-moment closure

A low Reynolds number second-moment closure has been used to calculate a turbulent boundary layer which develops over a riblet surface with zero pressure gradient. The calculated mean velocity distributions compare favourably with measurements. Calculated Reynolds stresses away from the riblet surface region are also in agreement with measurements. In the vicinity of the riblets, the model reflects the increased anisotropy of the Reynolds stress tensor inadequately. Possible reasons for this shortcoming are discussed and suggestions for improving the model are made.     

AN EXTENDED k-ε MODEL FOR NUMERICAL SIMULATION OF WIND FLOW AROUND BUILDINGS

It is assumed in this paper that for a high Reynolds number nearly homogeneouswind flow, the Reynolds stresses are uniquely related to the mean velocity gradientsand the two independent turbulent scaling parameters k and E. By applying dimensionalanalysis and owing to the Cayley-Hamilton theorem for tensors, a new turbulenceenclosure model so-called the axtended k-ε model has been developed. The coefficientsof the model expression were detemined by the wind tunnel experimental data ofhomogeneous shear turbulent flow. The model was compared with the standard k-εmodel in in composition and the prediction of the Reynold’s normal Stresses. Using thenew model the numerical simulation of wind flow around a square cross-section tallbuilding was performed. The results show that the extended k-ε model improves theprediction of wind velocities around the building the building and wind pressures on the buildingenvelope.     

Modelling turbulence around and inside porous media based on the second moment closure

To predict turbulence in porous media, a new approach is discussed. By double (both volume and Reynolds) averaging Navier–Stokes equations, there appear three unknown covariant terms in the momentum equation. They are namely the dispersive covariance, the macro-scale and the micro-scale Reynolds stresses, in the present study. For the macro-scale Reynolds stress, the TCL (two-component-limit) second moment closure is applied whereas the eddy viscosity models are applied to the other covariant terms: the Smagorinsky model and the one-equation eddy viscosity model, respectively for the dispersive covariance and the micro-scale Reynolds stress. The presently proposed model is evaluated in square rib array flows and porous wall channel flows with reasonable accuracy though further development is required.     

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.     

Break-Up of Aerosol Agglomerates in Highly Turbulent Gas Flow

Agglomerate aerosols in a turbulent flow may be subjected to very high turbulent shear rates which through the generation of lift and drag can overcome the adhesive forces binding the constituents of an agglomerate together and cause it to break-up. This paper presents an analysis of the experimental measurements of the breakup of agglomerates between 0.1?C10???m in size, in a turbulent pipe flow followed by an expansion zone with a Reynolds numbers in the range 10⁵ to 10⁷. The analysis shows that even in wall bounded turbulence, the high turbulent shear stresses associated with the small scales of turbulence in the core can be the main source of breakup preceding any break-up that may occur by impaction at the wall. More importantly from these results, a computationally fast and efficient solution is obtained for the General Dynamic Equation (GDE) for agglomerate transport and breakup in highly turbulent flow. Furthermore the solution for the evolution of the aerosol size distribution is consistent with the experimental results. In the turbulent pipe flow section, the agglomerates are exposed continuously to turbulent shear stresses and experience more longer term breakup than in the expansion zone (following the pipe flow) where the exposure time is much less and break-up occurs instantaneously under the action of very high local turbulent shear stresses. The validity of certain approximations made in the model is considered. In particular, the inertia of the agglomerates characterised by a Stokes Number from 0.001 for the smallest particles up to 10 for 10???m particles and the fluctuations of the turbulent shear stresses are important physical phenomena which are not accounted for in the model.     

CALCULATION OF WALL-BOUNDED COMPLEX FLOWS AND FREE SHEAR FLOWS

Various wall-bounded flows with complex geometries and free shear flows have been studied with a newly developed realizable Reynolds stress algebraic equation model. The model development is based on the invariant theory in continuum mechanics. This theory enables us to formulate a general constitutive relation for the Reynolds stresses. Pope (J. Fluid Mech., 72 , 331–340 (1975)) was the first to introduce this kind of constitutive relation to turbulence modelling. In our study, realizability is imposed on the truncated constitutive relation to determine the coefficients so that, unlike the standard k–ϵ eddy viscosity model, the present model will not produce negative normal stresses in any situations of rapid distortion. The calculations based on the present model have shown encouraging success in modelling complex turbulent flows.     

用电化学方法测试动脉模型壁面剪应力

应用电化学方法,对动脉模型Ｔ型分叉部位流场壁面剪应力进行测试研究。测试了对于现有理论分析和数值计算都比较困难的高雷诺数（ＲＥ＝１０００－２０００）流动流场的壁面剪应力,并且对苦干不同雷诺数及不同支管分流情况进行了系列测试。通过实验发现,此部痊同时存在高剪应力区和低剪应力区,确定了它们的位置和剪应力的大上。系列测试还显示：随着雷诺数的变化,无量纲管应力有一定的变化;而当支管分流变化时,无量纲剪应力的     

Structure of the nonisothermal swirling gas-droplet flow behind an abrupt tube expansion

Flow structure and heat and mass transfer in a swirling two-phase stream is numerically modeled using the Reynolds stress transport model. The gas phase is described by the 3DRANS system of equations with account for the inverse influence of particles on the transport processes in the gas. The gas phase turbulence is calculated using the Reynolds stress transport model with account for the presence of disperse particles. The two-phase nonswirling flow behind an abrupt tube expansion contains a secondary corner vortex which is absent from the swirling flow. The disperse phase is redistributed over the tube cross-section. Large particles are concentrated in the wall region of the channel under the action of the centrifugal forces, while the smaller particles are in the central zone of the chamber.     

Anisotropy measurements in the boundary layer over a flat plate with suction

Laser Doppler velocity measurements are carried out in a turbulent boundary layer subjected to concentrated wall suction (through a porous strip). The measurements are taken over a longitudinal distance of 9× the incoming boundary layer thickness ahead of the suction strip. The mean and rms velocity profiles are affected substantially by suction. Two-point measurements show that the streamwise and wall-normal autocorrelations of the streamwise velocity are reduced by suction. It is found that suction alters the redistribution of the turbulent kinetic energy k between its components. Relative to the no-suction case, the longitudinal Reynolds stress contributes more to k than the other two normal Reynolds stresses; in the outer region, its contribution is reduced which suggests structural changes in the boundary layer. This is observed in the anisotropy of the Reynolds stresses, which depart from the non-disturbed boundary layer. With suction, the anisotropy level in the near-wall region appears to be stronger than that of the undisturbed layer. It is argued that the mean shear induced by suction on the flow is responsible for the alteration of the anisotropy. The variation of the anisotropy of the layer will make the development of a turbulence model quite difficult for the flow behind suction. In that respect, a turbulence model will need to reproduce well the effects of suction on the boundary layer, if the model is to capture the effect of suction on the anisotropy of the Reynolds stresses.     

THE EXAMINATION OF TURBULENCE MODELING WITH LES DATABASE

ＴＨＥＥＸＡＭＩＮＡＴＩＯＮＯＦＴＵＲＢＵＬＥＮＣＥＭＯＤＥＬＩＮＧＷＩＴＨＬＥＳＤＡＴＡＢＡＳＥＳｕＭｉｎｇ－ｄｅ（苏铭德）（ＱｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ）Ｒ．Ｆｒｉｅｄｒｉｃｈ（ＬｅｈｒｓｉｒｈｌｆｕｒＦｌｕｉｄｓｍｅｃｈａ...