Study on the dynamic characteristics of flow over building cluster at high Reynolds number by large eddy simulation

In this paper,the dynamic characteristics of building clusters are simulated by large eddy simulation at high Reynolds number for both homogeneous and heterogeneous building clusters.To save the computational cost a channel-like flow model is applied to the urban canopy with free slip condition at the upper boundary.The results show that the domain height is an important parameter for correct evaluation of the dynamic characteristics.The domain height must be greater than 8h(h is the average building height)in order to obtain correct roughness height while displacement height and roughness sublayer are less sensitive to the domain height.The Reynolds number effects on the dynamic characteristics and flow patterns are investigated.The turbulence intensity is stronger inside building cluster at high Reynolds number while turbulence intensity is almost unchanged with Reynolds number above the building cluster.Roughness height increases monotonously with Reynolds number by 20%from Re*=103 to Re*=105 but displacement height is almost unchanged.Within the canopy layer of heterogeneous building clusters,flow structures vary between buildings and turbulence is more active at high Reynolds number.     

Simulating high Reynolds number flow in two-dimensional lid-driven cavity by multi-relaxation-time lattice Boltzmann method

By coupling the non-equilibrium extrapolation scheme for boundary condition with the multi-relaxation-time lattice Boltzmann method, this paper finds that the stability of the multi-relaxation-time model can be improved greatly, especially on simulating high Reynolds number (Re) flow. As a discovery, the super-stability analysed by Lallemand and Luo is verified and the complex structure of the cavity flow is also exhibited in our numerical simulation when Re is high enough. To the best knowledge of the authors, the maximum of Re which has been investigated by direct numerical simulation is only around 50,000 in the literature; however, this paper can readily extend the maximum to 1000,000 with the above combination.     

高超临界雷诺数区间内二维圆柱绕流的实测研究

实测强风工况下高度167 m的徐州彭城电厂冷却塔的表面风荷载,并归纳历史上其他研究人员给出的实测结果,以丰富高超临界雷诺数(Re)区间二维圆柱绕流的试验成果.在低湍流度均匀流场和高湍流度大气边界层流场中分别开展4种风速8类粗糙度条件下的冷却塔刚性模型测压风洞试验,通过对比低雷诺数(Re=2.1×10~5—4.19×10~5)条件下的风洞试验结果和高雷诺数(Re=5.4×10~7—1×10~8)条件下的现场实测结果研究各种静动态绕流特征随雷诺数的变化规律,重点考察雷诺数无关现象的产生条件.研究结果表明,对于物表相对粗糙度在0.01以上的圆柱绕流,雷诺数不相关现象存在于很宽的雷诺数范围(2×10~5Re1×10~8)内;增大来流湍流度亦能引起的雷诺数无关现象,但此时该现象可能仅存在于一个较窄的低雷诺数范围内.     

有限雷诺数情况下微槽道可压缩流动滑移流模型稳定性的数值模似

通过对微槽道一阶滑移流模型稳定性的理论研究，发现在不可压缩流和雷诺数趋于零的情况下，当Knudsen数Kn大于1／9时，会出现不稳定现象．这一理论结论在可压缩流和有限雷诺数情况下是否成立，尚不清楚，采用成熟的计算方法，用数值计算来研究可压缩流和有限雷诺数情况下一阶滑移流模型的稳定性，发现随着雷诺数的减小，不稳定的临界Knudsen数也不断减小，逐渐趋近于在不可压缩流和雷诺数趋于零的情况下得到的理论分析结果。     

A restricted nonlinear large eddy simulation model for high Reynolds number flows

Recent numerical studies of the restricted nonlinear (RNL) model have demonstrated its ability to reproduce important features of wall turbulence despite a severe reduction in the number of degrees of freedom. In these prior studies, the RNL model included full resolution of the viscous term. In this work, we extend the RNL model to arbitrarily high Reynolds numbers by developing a RNL large eddy simulation (LES) framework along with a method to systematically identify an appropriate streamwise wavenumber support based on spectral properties of wall turbulence. This method leads to a band-limited RNL–LES system which is successful in reproducing some of the most important statistical features captured in previous low to moderate Reynolds number simulations, e.g. the mean velocity and second-order moment profiles. The RNL–LES framework offers a new approach to understanding the connection between coherent structures and the momentum transfer mechanisms of wall turbulence at arbitrarily high Reynolds numbers, where resolution of the viscous terms can become computationally expensive even with the relatively low computational complexity afforded through the dynamical restriction of the RNL model.     

Stabilized finite element method for incompressible flows with high Reynolds number

In the following paper, we discuss the exhaustive use and implementation of stabilization finite element methods for the resolution of the 3D time-dependent incompressible Navier–Stokes equations. The proposed method starts by the use of a finite element variational multiscale (VMS) method, which consists in here of a decomposition for both the velocity and the pressure fields into coarse/resolved scales and fine/unresolved scales. This choice of decomposition is shown to be favorable for simulating flows at high Reynolds number. We explore the behaviour and accuracy of the proposed approximation on three test cases. First, the lid-driven square cavity at Reynolds number up to 50,000 is compared with the highly resolved numerical simulations and second, the lid-driven cubic cavity up to Re = 12,000 is compared with the experimental data. Finally, we study the flow over a 2D backward-facing step at Re = 42,000. Results show that the present implementation is able to exhibit good stability and accuracy properties for high Reynolds number flows with unstructured meshes.     

出口雷诺数对平面射流自保持性的影响

通过实验研究出口雷诺数对平面湍流射流自保持性的影响. 测量的射流来自相同的喷嘴但不同的雷诺数Re（≡U_jh/ν,其中U_j是出口平均速度、h是窄缝出口的厚度和ν是黏性系数）,其变化范围是Re=4582—57735.所得的数据包括沿轴线的平均速度、湍流强度、积分尺度、高阶矩和能谱. 实验发现,随着Re的增大,平面射流发展减慢,平均速度和湍流强 关键词： 平面射流 雷诺数 自保持性     

雷诺数对圆形渐缩喷嘴湍流射流的影响

本文研究雷诺数(Re)对圆形渐缩喷嘴湍流射流的影响.实验在射流出口雷诺数为 Re = 4050—20100 的范围内进行,分别测量了射流出口、中心线的平均及湍流流场以及部分径向剖面速度分布.所有测量均采用单热线恒温热线风速仪进行高频采样,所测流场范围在轴向上为 0—30d(这里d为射流出口直径).虽然出口速度分布均为"平顶帽"形,但测量结果依然反映出Re对射流出口以及下游流场有强烈的影响.当Re小于临界值(～10000)时 关键词： 雷诺数 圆形射流 热线风速仪     

Effects of upstream roughness and Reynolds number on separated and reattached turbulent flow

This paper presents experimental investigation of upstream roughness and Reynolds number effects on the recirculation region over a smooth forward facing step. The upstream rough wall was produced from 1.5 mm sand grains and the Reynolds number based on step height, Re_h, was varied from 2040 to 9130 for both the upstream smooth and rough walls. For the smooth wall, the reattachment length increased monotonically with Re_h to an asymptotic value of 2.2 step heights for Re_h ≥ 6380. Upstream roughness reduced the reattachment length by 44% because of larger momentum deficit and higher turbulence level in the rough wall boundary layer. The mean velocities and Reynolds stresses were also reduced by roughness. The Reynolds shear stress and production of turbulent kinetic energy showed high negative values at the leading edge of the step indicating counter-gradient diffusion. The implications of these results for standard eddy viscosity models are discussed.     

大雷诺数Navier-Stokes方程的两水平亚格子模型稳定化方法

基于两重网格离散方法,提出三种求解大雷诺数定常Navier-Stokes方程的两水平亚格子模型稳定化有限元算法.其基本思想是首先在一粗网格上求解带有亚格子模型稳定项的Navier-Stokes方程,然后在细网格上分别用三种不同的校正格式求解一个亚格子模型稳定化的线性问题,以校正粗网格解.通过适当的稳定化参数和粗细网格尺寸的选取,这些算法能取得最优渐近收敛阶的有限元解.最后,用数值模拟验证三种算法的有效性.     

高马赫数高压力的氧碘化学激光器阵列喷管的数值模拟

采用计算流体力学方法,研究了以氮气为载气的新型高总压氧碘化学激光器(COIL)阵列喷管。模拟结果表明:采用高马赫数的氮气流引射低马赫数的氧气流,可以提高光腔出口的驻点压力;高超声速的氮气与声速的氧气混合较慢,在喷管出口安装翼片有利于增强气流混合;喷管出口安装大翼片,翼片诱导的横向涡可以到达氮喷管的中心,光腔内混合比较充分。通过采用10组分21反应的化学反应模型,模拟了阵列喷管内多组分气体的混合和化学反应过程。模拟结果表明：光腔内生成了激发态碘原子和基态碘原子,光腔中获得了正增益,而且光腔出口的总压也由2.6 kPa提升至28.9 kPa。     

高马赫数高压力的氧碘化学激光器阵列喷管的数值模拟

 采用计算流体力学方法,研究了以氮气为载气的新型高总压氧碘化学激光器(COIL)阵列喷管。模拟结果表明:采用高马赫数的氮气流引射低马赫数的氧气流,可以提高光腔出口的驻点压力;高超声速的氮气与声速的氧气混合较慢,在喷管出口安装翼片有利于增强气流混合;喷管出口安装大翼片,翼片诱导的横向涡可以到达氮喷管的中心,光腔内混合比较充分。通过采用10组分21反应的化学反应模型,模拟了阵列喷管内多组分气体的混合和化学反应过程。模拟结果表明：光腔内生成了激发态碘原子和基态碘原子,光腔中获得了正增益,而且光腔出口的总压也由2.6 kPa提升至28.9 kPa。     

Influences of initial velocity, diameter and Reynolds number on a circular turbulent air/air jet

This paper assesses the suitability of the inflow Reynolds number defined by Re_o ≡ U_oD/ν (here U_o and D are respectively the initial jet velocity and diameter while ν is kinematic viscosity) for a round air/air jet. Specifically an experimental investigation is performed for the influences of U_o, D and Re_o on the mean-velocity decay and spread coefficients (K_u, K_r) in the far field of a circular air jet into air from a smoothly contracting nozzle. Present measurements agree well with those previously obtained under similar inflow conditions. The relations K_u ∝ U_o and K_r ∝ 1/U_o for U_o < 5 m/s appear to work, while each coefficient approaches asymptotically to a constant for U_o > 6 m/s, regardless of the magnitudes of Re_o and D. It is revealed that Re_o may not be an appropriate dimensionless parameter to characterize the entire flow of a free air/air jet. This paper is the first paper that has challenged the suitability of Re_o for turbulent free jets.     

颗粒Reynolds正应力轨道模型及后台阶两相流动数值模拟

用随机过程理论建立气固两相耦合脉动量Lagrange方程,并建立了气固两相流随机颗粒轨道模型中颗粒Reynolds正应力的Lagrange方程.将新的模型运用于各向同性的湍流衰减的流场中,模拟颗粒的湍流扩散特性,与文[1]中的模型和实验结果作比较,并使用新模型模拟了后台阶两相流动.     

Numerical analysis of fluid flow through a cylinder array using a lattice Boltzmann model

In this paper we present a detailed computational study of an incompressible Newtonian fluid flow across a periodic array of two-dimensional cylinders which is a simplest non-trivial representation of a porous media. A two-dimensional Lattice Boltzmann Method is used to solve the governing Navier－Stokes equation taking into account of viscous dissipation effects and influence of nonlinear fluid drag. Both the flow fields and the Darcy－Forchheimer drag coefficient as a function of the solid volume fraction are calculated for a wide range of flow Reynolds numbers. The predictions were compared with the results from conventional numerical and empirical models for verification. Apart from confirming that inertial effects can cause a significant deviation from Darcy's law for large velocities the results also show that the characteristics of the vorticity field vary considerably as the Reynolds number increases, which will have major implications to the transport of passive particulate substances within the pores and their removal rate.     

A multi-scale simulation method for high Reynolds number wall-bounded turbulent flows

We present an assessment and enhancement of the hybrid two-level large-eddy simulation method (A.G. Gungor and S. Menon, A new two-scale model for large eddy simulation of wall-bounded flows, Prog. Aerosp. Sci. 46 (2010), pp. 28–45), a multi-scale formulation for simulation of high Reynolds number wall-bounded turbulent flows. The assessment of the method is performed by examining role of static and dynamic blending functions used to perform hybridisation of two-level simulation (K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows, J. Comput. Phys. 220 (2006), pp. 290–311; K. Kemenov and S. Menon, Explicit small-scale velocity simulation for high-Re turbulent flows. Part 2: Non-homogeneous flows, J. Comput. Phys. 222 (2007), pp. 673–701) and large-eddy simulation methods. The sensitivity of first- and second-order turbulence statistics to the type of blending functions is investigated by simulating a fully developed turbulent flow in a channel at a friction Reynolds number Re_τ = 395 and comparing the results with those obtained using a direct numerical simulation. The first-order statistics do not show any significant differences for different blending functions, but the second-order statistics show some minor differences. The dynamic evaluation of the hybrid region and the blending function is necessary for non-equilibrium and complex flows where use of a static blending function can lead to inaccurate results. We propose two criteria for the dynamic evaluation; first evaluates extent of the hybrid region based on the subgrid turbulent kinetic energy and the second estimates the blending function based on a characteristic length scale. The computational efficiency of the method is enhanced by incorporating a hybrid programming paradigm where a standard domain decomposition by the message-passing-interface library is combined with the open multi-processing based parallelisation. A further enhancement of the method is achieved by incorporating a closure model for the unclosed hybrid terms in the governing equations, which appear due to hybridisation of two-level- and large-eddy-simulation methods. The model is based on an order of magnitude approximation and a preliminary assessment of the model shows improvement of turbulence statistics when used to simulate turbulent flow in a periodic channel. The assessment and improvements to the multi-scale method make it more suitable for simulation of practical wall-bounded turbulent flows at higher Reynolds number than a conventional large-eddy simulation. This is demonstrated by simulating two representative cases; turbulent flow at high Reynolds number in a periodic channel and flow over a bump placed on the lower surface of a channel, where a relatively coarser computational grid is found to be sufficient for reasonably accurate results.     

Direct numerical simulation of transitional flow in a finite length curved pipe

Transition to turbulent flow in a curved pipe has been well studied through experiments and numerical simulations. Numerical simulations often use a helical pipe with an infinite length such that the inlet and outlet boundary conditions can be modelled as periodic which greatly reduces computational time. In this study, we examined a finite length curved pipe with Poiseuille flow imposed at the inlet and a stress-free boundary condition at the outlet. Direct numerical simulation of the Navier-Stokes equations for rigid walls and a Newtonian fluid was performed using nek5000. Straight extensions were added to the inlet and outlet such to diminish the impact of boundary conditions on the flow field in the region with curvature. The examined model has a pipe radius of curvature that is three times the pipe radius. The model has ～355 million nodes and required an order of magnitude greater computational time when compared with an infinite length curved pipe. Results show that the critical Reynolds number, the lowest value with instabilities present in the flow, is much greater than that of a straight pipe and occurs near Re=5000–5200. This is larger than the critical Reynolds number typically reported for an infinite length curved pipe (Re=4200–4300).     

Numerical simulation of the flow within and over an intersection model with Reynolds-averaged Navier--Stokes method

In this study, the Reynolds-averaged Navier--Stokes (RANS) method is employed to simulate the flow within and over an intersection model with three kinds of $k$--$\varepsilon$ turbulence closure schemes, namely, standard model, renormalization group (RNG) model and realizable $k$--$\varepsilon$model. The comparison between the simulated and observed flow fields shows that the RANS simulation with all the three turbulence models cannot completely and accurately reproduce the observed flow field in all details. A detailed comparison between the predicted profiles of wind velocities and the measured data shows that the realizable $k$--$\varepsilon$ model is the best one among the three turbulence closure models in general. However, the extent to which the improvement is achieved by the realizable $k$--$\varepsilon$ model is still not enough to completely and accurately describe the turbulent flow in a relatively complex environment.     

A lattice Bhatnagar--Gross--Krook model for a class of the generalized Burgers equations

A new lattice Bhatnagar--Gross--Krook (LBGK) model for a class of the generalized Burgers equations is proposed. It is a general LBGK model for nonlinear Burgers equations with source term in arbitrary dimensional space. The linear stability of the model is also studied. The model is numerically tested for three problems in different dimensional space, and the numerical results are compared with either analytic solutions or numerical results obtained by other methods. Satisfactory results are obtained by the numerical simulations.     

����������IJ���������ģ����ֵģ��

采用FLUENT软件分别对外加均匀横向磁场的等截面三维充分发展液态金属管流的层流模型和低雷诺数湍流Lam/Bremhost(LB)模型进行了数值模拟,分析了外加磁场对普通方管LB模型速度分布和压降的影响。比较在相同哈特曼数下,层流和湍流模型方管截面上速度分布和管道中MHD压降。其中,对电流的计算采用磁感应方程来求得。数值模拟结果证明了用低雷诺数LB湍流模型解决方管磁流体流动的可行性。通过层流模型和湍流模型的对比可知,层流模型有较短的入口长度,但管内流体的压降却很大;而湍流模型管内速度更加平均化,管内压降较小,但管内入口长度较长。