Numerical simulation for rotating internal weakly viscoelastic flows in rectangular ducts

The present work develops a numerical method for the solution of rotating internal weakly viscoelastic flows in rectangular ducts for dimensionless parameters such as the Reynolds, Rossby and Weissenberg numbers, taken respectively in the intervals between 171 and 12000, 0.047 and 1/12 and up to 1/10000. It is shown that the usual counter‐rotating double‐vortex configuration of secondary flow breaks down with the increase of the Reynolds number (over the threshold of 171). For higher Reynolds numbers such as 7500 and 12000 the secondary flow diffuses to the interior of the duct where it assumes a fully developed configuration and the transition to the turbulence structure is observed. The Sobolev norms increase almost proportionally to the increase of the Reynolds number, and play an essential role for more complex problems involving transition to turbulence modelling. Copyright © 2002 John Wiley & Sons, Ltd.     

Moment of Resistance of a Disk Rotating in a Closed Axisymmetric Cavity

A turbulent flow in a closed axisymmetric cavity with a rotating disk is considered. The moment of resistance to disk rotation is calculated as a function of the dimensionless gap between the motionless casing and rotating disk and of the Reynolds number. Results calculated on the basis of different turbulence models are compared with data of a physical experiment and with available correlations. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 153–160, January–February, 2006.     

气体动理学格式与多尺度流动模拟

介绍了气体动理学格式(GKS)的基本构造原理及其在两种典型多尺度流动模拟中的应用。GKS利用介观BGK方程的跨尺度演化解来构造网格界面上的数值通量,从而发展出能随计算网格尺度变化自动切换物理模型的多尺度方法。对湍流这种宏观多尺度流动,发展了高精度GKS方法并成功用于低雷诺数湍流的直接数值模拟;为实现对高雷诺数湍流的高效精细模拟,基于拓展BGK方程和已有的RANS,LES模型建立了新型多尺度模拟框架。对跨流域稀薄流动,发展了适合大规模并行的三维统一气体动理学格式(UGKS),并建立了适合轴对称稀薄流动的UGKS。研究表明,GKS在多尺度流动高效模拟中的优异性能,具有很好的发展前景。     

Some Swirling-flow Challenges for Turbulent CFD

The paper examines some of the continuing challenges, within a RANS framework, of computing turbulent swirling flows such as are encountered in industry and the environment. The principal focus is on modelling turbulent transport processes but serious problems also arise in handling numerical issues, too. Recent researches of two of these types of flow by the authors and their colleagues in the Turbulence Mechanics Group at Manchester are examined; namely, the confined flow within a rotor–stator disc cavity and the trailing wing-tip vortex. The former flow, while geometrically axisymmetric, has been found to create multiple rotating vortices necessitating a three-dimensional time-dependent analysis. The wing-tip vortex is extremely sensitive to the choice of turbulence model and only a second-moment closure that complies with the constraints of two-component turbulence has been found capable of handling both the flow over the wing and the wake vortex. Moreover, because of the large distances downstream of the aircraft to which, for practical cases, computations need to be carried, the numerical strategy is brought into question. Finally, arising from these two test cases, outline remarks are made about a swirling flow that poses one of the major computational challenges of the twenty-first century.     

Mixing in horizontally heterogeneous flows

An experimental study of mixing across density interfaces produced by laterally heterogenous turbulence is presented in this paper. The turbulence is generated by a flow or air bubles rising through a density interface produced by brine and fresh water. The mixing efficiency, , of the process is measured comparing the increase in potential energy with the available kinetic energy. We find that there is a decrease in the global mixing efficiency of the process with the length of the tank, the shape of (Ri) depends also on the air flow producing the turbulence, showing a geometrical limit to the ammount of kinetic energy which may be used for mixing.     

统一二阶矩模型用于模拟旋流湍流两相流动

用统一二阶矩模型（ＵＳＭ）模拟了旋流数为０４７和１５的气粒两相流动，并和实验结果以及ｋ ε ｋｐ模型的模拟结果进行了对比．研究结果表明，提高旋流数减小了轴向速度反流区，增大了切向速度似固核区．ＵＳＭ和ｋ ε ｋｐ模型预报旋流数为０４７时的两相速度场差别不大，并都和实验结果接近，但前者预报的旋流数为１５的两相速度场比后者有改进，在两种情况下，前者都能揭示出后者无法预报的两相湍流各向异性规律．     

基于颗粒湍流和颗粒碰撞相互作用规律的湍流模型的研究

颗粒湍流和颗粒碰撞的相互作用规律是两相流动中的核心问题。用颗粒湍流模型和颗粒碰撞的动力论模型叠加的方法在研究两相湍流流动方面取得了一定的成效,但是还有待改进。本文基于颗粒湍流形成大尺度脉动和颗粒间碰撞引起小尺度脉动的概念,从双流体模型出发,建立了两相流动的双尺度kp-pε两相湍流模型。利用该模型对下行床和突扩室内的气固...     

Advances in studies on two-phase turbulence in dispersed multiphase flows

Particle/droplet/bubble fluctuation and dispersion are important to mixing, heat and mass transfer, combustion and pollutant formation in dispersed multiphase flows, but are insufficiently studied before the 90 years of the last century. In this paper, the present author reports his systematic studies within nearly 20 years on two-phase turbulence in dispersed multiphase flows, including particle fluctuation in dilute gas-particle and bubble-liquid flows, particle-wall collision effect, coexistence of particle turbulence and inter-particle collisions, fluid turbulence modulation due to the particle wake effect and validation of the two-fluid RANS modeling using large-eddy simulation.     

Rotation effects on a fully-developed turbulent pipe flow

A fully-developed turbulent pipe flow is allowed to pass through a rotating pipe section, whose axis of rotation coincides with the pipe axis. At the exit end of the rotating section, the flow passes into a stationary pipe. As a result of the relaxation of surface rotation, the turbulent flow near the pipe wall is affected by extra turbulence production created by the large circumferential shear strain set up by the rapid decrease of the rotational velocity to zero at the wall. However, the flow in the most part of the pipe is absent of this extra turbulence production because the circumferential strain is zero as a result of the solid-body rotation imparted to the flow by the rotating pipe section. The combined effect of these two phenomena on the flow is investigated in detail using hot-wire anemometry techniques. Both mean and turbulence fields are measured, together with the wall shear and the turbulent burst behavior at the wall. A number of experiments at different rotational speeds are carried out. Therefore, the effects of rotation on the behavior of wall shear, turbulent burst at the wall, turbulence production and the near-wall flow can be documented and analysed in detail.     

A turbulence model for inclined, bubbly flows

A new turbulence model for the flow of a two phase (liquid-liquid) flow in an inclined pipe is presented. An eddy viscosity is used to model the effects of shear induced turbulence and bubble induced turbulence. The cross-pipe momentum transport arising from the buoyant rise of bubbles across the axial flow is also modelled. Numerical simulations have been carried out in both one and two dimensions. One and two dimensional numerical simulations are presented.On leave from the University of Leeds, Leeds LS2 9JT, U.K.     

工程湍流模式理论综述及展望

本文讨论了国内外湍流模拟的现状和发展趋势．指出湍流模式的建立除了应遵循理性力学原则外,还必须密切结合工程流体的复杂流动现象,如对具有浮力的回流、分离流及强旋流的模拟,对逆梯度的输运模拟,对单相流、多相流,单流体、多流体的湍流牛顿流体及湍流的非牛顿流体的模拟．应加深对湍流机理的认识,改进湍流模拟手段,结合工程实际,提出较为通用的工程湍流模式．      

Kinematic simulation for stably stratified and rotating turbulence

The properties of one-particle and particle-pair diffusion in rotating and stratified turbulence are studied by applying the rapid distortion theory (RDT) to a kinematic simulation (KS) of the Boussinesq equation with a Coriolis term.Scalings for one- and two-particle horizontal and vertical diffusions in purely rotating turbulence are proposed for small Rossby numbers.Particular attention is given to the locality-in-scale hypothesis for two-particle diffusion in purely rotating turbulence both in the horizontal and the vertical directions. It is observed that both rotation and stratification decrease the pair diffusivity and improve the validity of the locality-in-scale hypothesis. In the case of stratification the range of scales over which the locality-in-scale hypothesis is observed is increased.It is found that rotation decreases the diffusion in the horizontal direction as well as, though to a much lesser extent, in the vertical direction.     

Studies on gas turbulence and particle fluctuation in dense gas-particle flows

Particle fluctuation and gas turbulence in dense gas-particle flows are less studied due to complexity of the phenomena. In the present study, simulations of gas turbulent flows passing over a single particle are carried out first by using RANS modeling with a Reynolds stress equation turbulence model and sufficiently fine grids, and then by using LES. The turbulence enhancement by the particle wake effect is studied under various particle sizes and relative gas velocities, and the turbulence enhancement is found proportional to the particle diameter and the square of velocity. Based on the above results, a turbulence enhancement model for the particle-wake effect is proposed and is incorporated as a sub-model into a comprehensive two-phase flow model, which is then used to simulate dilute gas-particle flows in a horizontal channel. The simulation results show that the predicted gas turbulence by using the present model accounting for the particle wake effect is obviously in better agreement with the experimental results than the prediction given by the model not accounting for the wake effect. Finally, the proposed model is incorporated into another two-phase flow model to simulate dense gasparticle flows in a downer. The results show that the particle wake effect not only enhances the gas turbulence, but also amplifies the particle fluctuation.     

Toward a Turbulence Constitutive Relation for Geophysical Flows

Rapidly rotating turbulent flows are frequently in approximate geostrophic balance. Single-point turbulence closures, in general, are not consistent with a geostrophic balance. This article addresses and resolves the possibility of a constitutive relation for single-point second-order closures for classes of rotating and stratified flows relevant to geophysics. Physical situations in which a geostrophic balance is attained are described. Closely related issues of frame-indifference, horizontal divergence, and the Taylor–Proudman theorem are discussed. It is shown that, in the absence of vortex stretching along the axis of rotation, turbulence is frame-indifferent. Unfortunately, no turbulence closures are consistent with this frame-indifference that is frequently an important feature of rotating or quasi-geostrophic flows. A derivation and discussion of the geostrophic constraint which ensures that the modeled second-moment equations are frame-invariant, in the appropriate limit, is given. It is shown that rotating, stratified, and shallow water flows are situations in which such a constitutive relation procedure is useful. A nonlinear nonconstant coefficient representation for the rapid-pressure strain covariance appearing in the Reynolds stress and heat flux equations, consistent with the geostrophic balance, is described. The rapid-pressure strain closure features coefficients that are not constants determined by numerical optimization but are functions of the state of turbulence as parametrized by the Reynolds stresses and the turbulent heat fluxes as is required by tensor representation theory. These issues are relevant to baroclinic and barotropic atmospheric and oceanic flows. The planetary boundary layers in which there is a transition, with height or depth, from a thermally or shear driven turbulence to a geostrophic turbulence is a classic geophysical example to which the considerations in this article are relevant. Received 14 October 1996 and accepted 9 June 1997     

Report on Turbulence and Mixing in Geophysical Flows II

We describe the presentations at the international meeting/workshop on Mixing in Geophysical Flows that took place at Vilanova i la Geltru, near Barcelona during the 20, 21 and 22nd of March 1997. There were more than 100 participants from 20 countries with 66 oral and poster presentations covering experimental and theoretical aspects of rotating and stratified fluids as well as field observations. The main topics discussed at the workshop were stratified flows, rotating stratified flows, gravity waves, instabilities and mixing, convection, experiments and numerical simulations of Geophysical flows and turbulent mixing. The papers are summarised in this report giving a state-of-the-art overview of present research in geophysical turbulent mixing.     

Particle and droplet deposition in turbulent swirled pipe flow

In this work we study deposition of particles and droplets in non-rotating swirled turbulent pipe flow. We aim at verifying whether the capability of swirl to enhance particle separation from the core flow and the capability of turbulence to efficiently trap particles at the wall can co-exist to optimize collection efficiency in axial separators. We perform an Eulerian–Lagrangian study based on Direct Numerical Simulation (DNS) of turbulence, considering the effect of different swirl intensities on turbulence structures and on particle transfer at varying particle inertia. We show that, for suitably-chosen flow parameters, swirl may be superimposed to the base flow without disrupting near-wall turbulent structures and their regeneration mechanisms. We also quantify collection efficiency demonstrating for the first time that an optimal synergy between swirl and wall turbulence can be identified to promote separation of particles and droplets.     

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Second-moment closure modelling of turbulence in a non-inertial frame

This article investigates the performance of second-moment closures in a rotating reference frame and presents a new closure for the rapid pressure-strain rate correlation based on the recently developed materially frame indifference principle. It is observed here that the existing second-moment closures with appropriate near-wall treatment can adequately predict flows in a rotating channel and in an axially rotating pipe for moderate rotation rate. Analysis of the newly proposed model indicates that it is capable of reflecting flow features under strong rotation.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

Recent progress in the development of the Elliptic Blending Reynolds-stress model

The Elliptic Blending Reynolds Stress Model (EB-RSM), originally proposed by Manceau and Hanjalić (2002) to extend standard, weakly inhomogeneous Reynolds stress models to the near-wall region, has been subject to various modifications by several authors during the last decade, mainly for numerical robustness reasons. The present work revisits all these modifications from the theoretical standpoint and investigates in detail their influence on the reproduction of the physical mechanisms at the origin of the influence of the wall on turbulence. The analysis exploits recent DNS databases for high-Reynolds number channel flows, spanwise rotating channel flows with strong rotation rates, up to complete laminarization, and the separated flow after a sudden expansion without and with system rotation. Theoretical arguments and comparison with DNS results lead to the selection of a recommended formulation for the EB-RSM model. This formulation shows satisfactory predictions for the configurations described above, in particular as regards the modification of the mean flow and turbulent anisotropy on the anticyclonic or pressure side.