对流边界层湍流特性的数值研究

采用大涡模拟方法研究了存在逆温层的情况下大气对流边界层的湍流特性。实际大气边界层中出现逆温层是较常见的，逆温层会导致大气边界层湍流结构的变化，从而影响大气的湍流扩散和输运特性。本文比较了不同逆温梯度的工况，着重分析了逆温层对边界层中热量逆梯度输运（counter gradient heat transportation，CGHT）的影响。计算结果表明：逆温梯度越大，对流边界层的发展越受到抑制；逆温层高度降低会影响整个对流边界层的温度抬升；逆温梯度越大，垂直速度方差越小；在逆温梯度较大的情况下，其逆梯度输运区域要略微低一些，初步分析认为是由于逆温层对热对流的抑制造成的；对于逆温层高度不同的情况，高度越低的逆温层对逆梯度输运的抑制作用更明显。     

大气边界层模拟的湍流相似

大气边界层的风洞模拟,是目前研究建筑物风载,污染质的大气扩散,以及解决某些环境科学和军事科学中大型课题的重要手段之一。在许多国家,风洞模拟已成为环境部门设计大型工程、林区、厂区和城区规划的主要依据。最初的大气边界层模拟装置,仅考虑形成一定的平均风速廓线,并未考虑到气流的湍流特性。微气象学的研究表明,大气的湍流特性是影响扩散的主要因素。特别是Monin和Obukhov通过大量实测资料提出大气表面层中的相似性理论,成为研究大气扩散的重要依据。从50年代开始,一些研究大气扩散和建筑物风载问题的长实验段风洞相继筹划和建立。其中较著名的有,美国科罗拉多州立大学的1.8米×1.8米×30米的回流式长实验段气      

气相和微颗粒污染物在城市环境中的传输扩散

城市大气边界层中污染物传输扩散, 包括气相和微颗粒(PM), 是一个日益关注的问题. 本文结合笔者及合著者近些年来的工作, 试图综述在这方面研究的最新进展. 主要包括: 天气尺度(meso-scale) 和街道尺度耦合的大涡模拟, 温度稳定性层结效应, 高层建筑的影响, 在极端天气下的污染物浓度极值的预报, 化学反应和湍流的耦合计算和气溶胶动力学过程. 本文侧重于介绍英国和欧洲的数值模型研究的进展情况.      

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

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

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应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

湍流研究国家重点实验室研究动态

湍流研究国家重点实验室研究动态梁彬北京大学湍流研究国家重点实验室,北京１００８７１湍流研究国家重点实验室一直是以基础研究为主,着眼于湍流发生发展的机制,从理论、计算和实验方面进行研究．实验室也兼顾应用基础研究,在环境流体力学、大气边界层模拟和风工程研...      

稳定分层流动中湍流特性的研究综述

本文总结了近60 年来分层流动中湍流特性研究的成果. 主要从两个方面进行了综述:(1) 分层流动中湍流场的演变和混合. 在这方面主要分析稳定分层对湍流混合和湍流结构的影响, 以及混合层内湍流结构的特性和混合层的演化规律. (2) 分层流动中湍流的扩散和输运. 动量和标量的逆梯度输运特性是分层湍流研究的一个重要方向;分析分层对湍流扩散的影响. 并指出了一些值得今后进一步研究的方向.     

一类平衡大气边界层边界条件研究

基于标准k-ε湍流模型,首先利用湍流粘度方程和剪切应力在整个边界层内恒定的假设,推导出一类耗散率表达式,并根据常用的湍动能入口剖面方程以及平均风速剖面方程,计算获得相应的耗散率方程;然后在输运方程中添加自定义源项,通过已经确定的平均速度方程、湍动能方程、耗散率方程计算得到相应输运方程的自定义源项表达式,并进行空风洞数值模拟,从而得到了一类满足平衡大气边界层的来流边界条件．通过将这种边界条件与由湍流平衡条件得到的边界条件进行比较,表明本方法获得的边界条件更适用．并且,本方法无需考虑修正壁面函数和修正湍流模型常数,因而计算更为简单,可为平衡大气边界层的研究提供一种新的思路．     

激波诱导的燃烧粉尘云边界层的结构

足够强的激波扫过铺有可燃粉尘的界面时，波后形成燃烧的粉尘云边界层。为揭示其内部结构，本文对该现象进行了实验和理论两个方面的研究。理论计算表明燃烧粉云边界层可分为三个区域：诱导区、反应区和扩散区，诱导区的粉尘浓度最高，计算获得的粉尘云轮廓和点火延迟与实验结果基本吻合．     

我国风工程研究现状与展望

 介绍了我国风工程的现状，提出了在大气边界层内风特性，风对 建筑物的作用，风对桥梁的作用，车辆空气动力特性，风力机空气动力 特性，风引起的污染扩散和质量迁移等研究中尚待解决的主要问题.     

Turbulent diffusion in a medium subjected to uniform deformation

Turbulent diffusion in media subjected to uniform deformation caused by the presence of average-velocity gradients which are constant throughout the space is an idealization of real processes, in particular, of processes such as diffusion in channels of variable cross-sectional area [1], in the lowest layers of the atmosphere [2], etc. In this article we formulate the problem of the connection between the statistical characteristics of the transfer of a passive substance in turbulent diffusion in deformed media with the statistical characteristics of the turbulence. The statistical transfer characteristics generally used are the first two moments of the vector of random displacement of a liquid particle under the action of turbulent pulsations in velocity: the average displacement and the components of the dispersion tensor of the displacement of a liquid particle. We obtain connecting relations for the dispersion tensor of a liquid particle in turbulent diffusion of a passive substance in a uniform turbulent medium subjected to uniform deformation caused by average-velocity gradients which are constant throughout the space. These relations are a generalization of known expressions for undeformed media [2, 3]. We investigate the case of rapid deformation when the turbulent characteristics of the medium vary in accordance with the linear theory [4].     

Experimental and numerical analysis of the jet dispersion from a bent chimney around an obstacle

An experimental study and a numerical modelling analysis were carried out simultaneously to study the flow field structure issuing from a chimney around an obstacle. The main purpose of this study is to evaluate the impact of the jet emitted from a chimney (bent or straight) on the dynamics and the turbulent features of the surrounding flow. The consideration of these features is particularly pertinent to the understanding of mixing between the interacting flows which may be very important in controlling pollutant dispersion in the atmosphere. The experimental data are depicted by means of a PIV technique; whereas the numerical three-dimensional model is simulated through the resolution of the different governing Navier–Stokes equations. The volume finite method, together with the second order turbulent closure model (RSM), was adopted. Variations in obstacle form (cylindrical or parallelepiped) and chimney configuration (bent or straight) were tested and features studied were: the global jet plume, the windward and leeward jet spread; the size, location and magnitude of the reverse flow region; the penetration and the deflection of the jet trajectory around the obstacle. All these considerations allowed us to characterize well the impact of the injection of the jet emitted from the chimney within the crossflow, and its spreading around the obstacle and within the whole domain. Such characterization is very important with regard to pollutant dispersion and consequently to the environmental impact. Indeed, the different species contained within the emitted fumes are mainly directed by the velocity components and their mixing and progression within the domain and around the obstacle are closely related.     

Large-eddy simulation of buoyancy-driven natural ventilation in an enclosure with a point heat source

Rising buoyant plumes from a point heat source in a naturally ventilated enclosure have been investigated using large-eddy simulation (LES). The aim of the work is to assess the performance and the accuracy of LES for modelling buoyancy-driven displacement ventilation of an enclosure and to shed more light on the transitional behaviour of the plume and the coherent structures involved. The Smagorinsky sub-grid scale model is used for the unresolved small-scale turbulence. The Rayleigh number, Ra is chosen to be in the range where spatial transition from laminar to turbulent flow takes place (Ra = 1.5 × 10⁹). The plume properties (source strength and rate of spread) as well as the ventilation properties (stratification height and temperature of stratified layer) estimated using the theory of Linden et al. are found to agree reasonably well with the LES results. The variation of the plume width with height indicates a linear variation of the entrainment coefficient rather than a constant value used by Linden et al. for a fully turbulent thermal plume. Flow visualisation revealed the nature of the large-scale coherent structures involved in the transition to turbulence in the plume. The most excited modes observed in the velocity, pressure and temperature fields spectra correspond to Strouhal number in the range 0.3 ≤ St ≤ 0.55 which is in agreement with those observed by Zhou et al. for a turbulent forced plume. Excited modes less than thisvalue (St = 0.2) were observed and may be due to low-frequency motions felt throughout the flow.     

DNS study of passive plume interference emitting from two parallel line sources in a turbulent channel flow

Turbulent mixing of dual plumes emitting simultaneously from line sources in a turbulent channel flow has been studied using direct numerical simulation (DNS). Three test cases have been compared to investigate the effects of the source separation on turbulent mixing of the two instantaneous plumes. The dispersion and interference of dual plumes are investigated in both physical and spectral spaces, which include an analysis of statistical moments of the concentration field, cross-correlation between the two instantaneous plumes, pre-multiplied spectra of the velocity and concentration fields, and co-spectrum and coherency spectrum of the dual plumes. As the downstream distance from the line source increases, the plume development associated with a single source emission transitions from a turbulent convective stage to a turbulent diffusive stage. It is observed that a plume released from a ground-level source reaches the turbulent diffusive stage faster than that released from an elevated source. It is also observed that a smaller separation between the two line sources tends to facilitate a more rapid growth in the cross-correlation coefficient of two instantaneous plumes. In the near-source region, the maximum coherency spectrum is produced at lower frequencies indicating that dual-plume mixing is dominated by the external flapping effects of large-scale eddy motions. However, in the far downstream region of the sources, the coherency spectrum in the higher frequency range increases significantly, indicating that the spread of the total plume is larger than all scales of turbulent eddies, such that they all contribute to the in-plume mixing of the dual plumes.     

Taylorian diffusion in mildly inhomogeneous turbulence

Taylor’s theory of dispersion was extended to produce estimates of the far–field growth rate of the plume of a passive scalar in grid turbulence (GT) and in uniformly sheared flow (USF), both of which evolve in the streamwise direction. Expressions for the evolution of the plume width relative to the integral length scale of turbulence were also derived. The predictions of plume growth rate were tested against available data in both of these types of flows and were found to be accurate in an extensive region of USF and compatible with an extrapolated trend in GT, in which the available data did not extend very far from the scalar source. Although in both cases the measured half-width of the plume was comparable in magnitude to the streamwise integral length scale of the turbulence, the far–field approximation seemed to hold.     

Coherent motion of turbulent thermal plume in stably stratified fluid

The purpose of this study is to clarify the existence of an ordered and large scale coherent motion in a turbulent plane thermal plume in a thermally-stable stratified fluid inside a comparatively large enclosure. First, the upper part of the thermal plume was carefully observed by a flow visualization. Secondly, a wave form of plume temperature variation was measured. Thirdly, a spectrum analysis was carried out on time series data of the thermal plume. Finally, physical characteristics were investigated on vortices in the thermal plume based on results of the wave form and the spectrum analysis of the plume temperature. As a result, the main conclusions are obtained as follows. (1) An existence of vortices near the upper part of the thermal plume was firstly found by careful flow visualization. (2) From the wave form of temperature variation and the spectrum analysis of the thermal plume, it was clarified that the vortices are generated in the transition state and are transported to the turbulent state. (3) The vortices are ordered and they behave as a large scale coherent motion in the turbulent thermal plume.     

大涡模拟建筑物对近源大气污染物扩散的影响

采用大涡模拟(LES)方法对一个由上风低源排放的烟流在孤立建筑物周围扩散的现象进行了数值模拟。模拟中,离散网格的尺度控制在100倍Kolmogorov湍流尺度的量级以内,亚格子应力由Smagorinsky模型加以模拟,并采用了随机数方法生成符合指定湍流特征的入流边界条件。数值模拟结果表明:建筑物尾流区平均浓度场的分布和大小不仅取决于局地流场的影响,在很大程度上也受到因建筑物阻挡出现的烟流抬升和马蹄涡对烟流的捕获的影响。通过对比风洞试验数据显示:大涡模拟方法能较好地再现气流绕建筑物的流动形态和建筑物周围的浓度分布,因而具有广泛的应用前景。     

Large-eddy simulation of a turbulent forced plume

This paper reports on an application of large-eddy simulation (LES) to a spatially-developing round turbulent buoyant jet. The numerical method used is based on a low-Mach-number version of the governing equations for compressible flow which can account for density variations. The second-order centre-difference scheme is used for spatial discretization and an Adams–Bashforth scheme for temporal discretization. Comparisons are made between LES results, experimental measurements and plume theory for the forced plume under moderate Reynolds number and good agreement has been achieved. It is found that the plume spreading and the centerline maximum mean velocity strongly depend on the forcing conditions imposed on the inflow plane. The helical mode of instability leads to a larger spreading rate as compared to an axisymmetric mode. The enhanced entrainment is directly related to the strong turbulent momentum and energy transports between the plume and surrounding fluid induced by vortex dynamics. The entrainment ratio is about 0.09 and falls into the range of experimentally determined values. Budgets of the mean momentum and energy equations are analyzed. It is found that the radial turbulent transport nearly balances the streamwise convection and the buoyancy force in the axial momentum equation. Also, the radial turbulent stress is balanced by the streamwise convection in the energy equation. The energy-spectrum for the axial velocity fluctuations shows a −5/3 power law of the Kolmogorov decay, while the power spectrum for the temperature fluctuations shows both −5/3 and −3 power laws in the inertial-convective and inertial-diffusive ranges, respectively.     

地球物理流体动力学的发展和室内实验研究的影响

系统介绍了大气、海洋与新兴的地质系统中的流体力学,尤其是通过室内实验的研究进展．详细描述了诸如羽流、湍流、卷挟、异重流、海洋环流、混合层、双扩散、盐指、锋面、岩浆库、地幔对流、板块运动、热斑等有趣的地球物理现象,它们对于了解因浮力引起的对流的机理是十分重要的．同时在气候、环境、灾害、成矿等领域有广泛的应用．论文反映了在 Ｇ．Ｉ． Ｔａｙｌｏｒ, Ｇ．Ｋ． Ｂａｔｃｈｅｌｏｒ领导下剑桥研究组的风格──通过小型室内实验了解机理,再用应用数学方法求解问题,他们对流体力学的发展作出了重大贡献．     

城市大气环境的大涡模拟研究进展

本文回顾城市大气污染扩散的研究实践与进展, 介绍城市大气环境流动的特点和用现代计算流体力学手段开展城市大气环境的大涡模拟研究方法, 包括: 数学模型, 控制方程、亚网格湍流模式、定解条件及其数值方法. 其中, 重点介绍中尺度到微尺度的耦合模型与适用于复杂城市下垫面(满足大涡模拟分辨率要求且计算量较小)的组合模型;并给出实际算例和结果分析, 包括准确度的统计估算和湍流特性等. 最后, 讨论城市大气环境数值模拟方法进一步改进的方向和城市大气环境流动与污染物扩散数值研究应重点关注的几个问题.