k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理, 完善等离子体激励器数学模型, 提升等离子体激励器扰动能力, 采用粒子图像测速技术, 在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究. 实验时, 将非对称布局激励器布置在平板模型上, 随后将带有激励器的模型放置在有机玻璃箱内, 从而避免环境气流对测试结果的影响. 基于激励器诱导流场, 分析了激励电压对诱导射流特性的影响, 揭示了较高电压下诱导射流近壁区的拟序结构, 获得了卷起涡、二次涡等拟序结构的演化发展过程, 计算了卷起涡脱落频率, 阐述了卷起涡与启动涡的区别, 初步探索了卷起涡的耗散机制. 结果表明: (1)层流射流不能完全概括等离子体诱导射流特性, 激励电压是影响射流特性的重要参数. 当电压较低时, 诱导射流为层流射流; 当电压较高时, 诱导射流的雷诺数提高, 射流剪切层不稳定, 层流射流逐渐发展为湍流射流. (2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构; 在固定电压下, 这些涡结构存在恒定的卷起频率. (3)当激励电压较高时, 流动不稳定使得卷起涡发生了拉伸、变形, 引起了流场湍动能增大, 从而加速了卷起涡的耗散. 研究结果为全面认识激励器射流特性, 进一步挖掘激励器卷吸掺混能力, 提升激励器控制能力积累基础.      

等离子体激励器诱导射流的湍流特性研究

为了进一步掌握等离子体流动控制机理,完善等离子体激励器数学模型,提升等离子体激励器扰动能力,采用粒子图像测速技术,在静止空气下开展了介质阻挡放电等离子体激励器诱导射流特性研究.实验时,将非对称布局激励器布置在平板模型上,随后将带有激励器的模型放置在有机玻璃箱内,从而避免环境气流对测试结果的影响.基于激励器诱导流场,分析了激励电压对诱导射流特性的影响,揭示了较高电压下诱导射流近壁区的拟序结构,获得了卷起涡、二次涡等拟序结构的演化发展过程,计算了卷起涡脱落频率,阐述了卷起涡与启动涡的区别,初步探索了卷起涡的耗散机制.结果表明:(1)层流射流不能完全概括等离子体诱导射流特性,激励电压是影响射流特性的重要参数.当电压较低时,诱导射流为层流射流;当电压较高时,诱导射流的雷诺数提高,射流剪切层不稳定,层流射流逐渐发展为湍流射流.(2)等离子体诱导湍流射流包含着卷起涡、二次涡等拟序结构;在固定电压下,这些涡结构存在恒定的卷起频率.(3)当激励电压较高时,流动不稳定使得卷起涡发生了拉伸、变形,引起了流场湍动能增大,从而加速了卷起涡的耗散.研究结果为全面认识激励器射流特性,进一步挖掘激励器卷吸掺混能力,提升激励器控制能力积累基础.     

充分发展自由湍流的直接数值模拟

讨论了充分发展自由剪切湍流在极限流动状态下数值模拟的方法。根据关于有结构的湍流的概念,对有序运动和大尺度涡的模拟提出了由非定常欧拉方程按基本网格和时间间隔作平均得到的离散耗散模型。借助于合适的耗散机理的概念反映了小尺度(小网格)脉动对不同尺度分辨率的贡献。湍流的随机分量根据动理学方程用蒙特-卡洛统计方法进行了模拟。所指出的方法大大地降低了对电子计算机能力的要求。所建议的算法的有效性,在有分离流物体的绕流问题中以及二维和空间流动的湍流尾流的演化问题中得到了证明。      

风口区空气流动的有限元模拟与实验研究

通风和空调出风口的射流特性主要取决于出风口截面之前的一段管道内的流动状况,研究送风口的射流偏转特性有重要实用价值.应用湍流大涡模拟技术结合Taylor-Galerkin有限单元法对工程中常见工况下的风口区管道内的三维流场进行了数值模拟和实验研究,数值计算结果与实验结果符合较好.表明湍流大涡模拟方法适合于边界形状复杂,存在各向异性的大尺度涡的内流情况,能可靠地预报风口区空气流动的射流偏转特性.     

具有壁面射流激励的圆管内气相流动大涡模拟

采用大涡模拟方法和单方程亚格子模式对小尺度量进行模拟。研究了不同强度壁面射流激励对圆管内气相流动的影响,模拟结果给出了射流对瞬态拟序结构发展、时平均流向速度分布的影响。随着射流强度的增加,射流入口附近流体的回流现象增强。射流强度足够大时可以减小管壁处的切应力值,同时会减小壁面附近流动速度,这种速度分布会导致气体夹带颗粒的能力下降,从而在实际两相流动中容易造成壁面附近的气粒返混现象。     

数值耗散对压气机流动分离涡模拟的影响研究

分离涡模拟DES是压气机流动模拟中常用的高保真湍流模式。为了使DES准确解析湍流,数值耗散必须限制在合理范围内。然而,当前的压气机流动DES类研究工作中仍然普遍采用高耗散的迎风格式。本文首先基于DES类方法计算的各向同性衰减湍流结果,定量比较了多种不同数值格式的耗散,证实了高耗散迎风格式严重低估中高波数湍流能量。高阶重构格式可以一定程度上改善该问题,但能量耗散仍然过高。本文在高阶重构的基础上,进一步引入自适应耗散函数修改Riemann求解器,构造了自适应耗散格式。该格式在全波数范围都能准确地预测湍流能谱。将该格式配合DES类方法模拟跨声速离心压气机流动,其预测的压比相比于三阶迎风格式,更加接近实验结果。此外,自适应耗散格式显著提高了中小尺度流动结构的分辨率。分析表明,在使用DES类方法模拟压气机流动时,有必要采用数值耗散较低的离散格式,以准确预测压气机总体性能和流动结构。本文构造的自适应耗散格式是一种良好选择。     

入口速度比对射流拟序结构的影响

为了深入了解湍流流动机理以及湍流拟序结构发现过程的影响因素，本文采用大涡模拟方法对不同入口射流伴流速度比的平面湍射流流动进行了数值模拟。采用分步投影法求解动量方程，亚格子项采用标准Smagorinsky亚格子模式模拟，压力泊松方程采用修正的循环消去法快速求解，空间方程采用二阶精度的差分格式，在时间方向上采用二阶精度的显式差分格式。模拟结果给出了平面射流中湍流拟序结构的瞬态发展演变过程，分析了入口速度比对射流拟序结构发展演化过程及宏观流场形态的影响。为进一步研究射流拟序结构及其在湍流流动中的作用提供了基础。     

关于植被中湍流的研究

本文从实验和模式理论两个方面介绍了近年来在植被湍流研究方面的进展.实验研究表明,植被内湍流是高度间歇性的,大尺度涡在湍流输运中起主导作用.植物枝、干、叶打碎了大尺度涡,产生的迹湍流中的小尺度涡更容易耗散成热。因此,植被内能量不按一般的能量级串理论预示的方式进行,而要发生所谓“短路”现象,使得功率谱曲线的斜率在惯性区比-2/3律更负.现有的高阶封闭模式虽然取得了很大成功,但在模拟湍流强度上高估了其大小,还需根据实验加以改进.      

大涡模拟及其在湍流燃烧中的应用

大涡模拟作为一种研究湍流流动和湍流燃烧的有效手段,在国际上已经得到广泛应用。本文在回顾了大涡模拟(LES)的基本思想及其实施方法的基础上着重介绍了前人在大涡模拟的亚格子湍流模式和亚格子燃烧模式中的研究成果,同时给出了采用不同亚格子模式的大涡模拟在湍流燃烧中的应用实例,指出了大涡模拟在湍流燃烧中的重要作用,为大涡模拟的进一步发展和应用提供参考。      

Mean Flow and Turbulence Measurements in a Turbulent Free Cruciform Jet

The results of measurements of all three components of the mean velocity vector, the Reynolds normal and primary shear stresses and the mean static pressure in a turbulent free jet, issuing from a sharp-edged cruciform orifice, are presented in this paper. The measurements were made with an x-array hot-wire probe and a pitot-static tube in the near flow field of the jet. The Reynolds number, based upon the equivalent diameter of the orifice, was 1.70 × 10⁵. In addition to the quantities measured directly, the mean streamwise centreline velocity decay, the jet half-velocity widths, the jet spreading rate, the mean streamwise vorticity, the mass entrainment rate, the integral momentum flux and the one-dimensional energy spectra have been derived from the measured data. The results show that the mean streamwise centreline velocity decay rate of the cruciform jet is higher than that of a round jet issuing from an orifice with the same exit area as that of the cruciform orifice. The mean streamwise velocity field changed shape continuously from a cruciform close to the orifice exit plane to circular at 12 and half equivalent diameters downstream. The mean streamwise vorticity field, up to about three equivalent diameters downstream of the orifice exit plane, consists of four pairs of counter-rotating cells, which are aligned with the four edges in the centre of the cruciform orifice.     

Measurements in the near flow field of an isosceles triangular turbulent free jet

The near field mean flow and turbulence characteristics of a turbulent jet of air issuing from a sharp-edged isosceles triangular orifice into still air surroundings have been examined experimentally using hot-wire anemometry and a pitot-static tube. For comparison, some measurements were made in an equilateral triangular free jet and in a round free air jet, both of which also issued from sharp-edged orifices. The Reynolds number, based on the orifice equivalent diameter, was 1.84×10⁵ in each jet. The three components of the mean velocity vector, the Reynolds normal and primary shear stresses, the one-dimensional energy spectra of the streamwise fluctuating velocity signals and the mean static pressure were measured. The mean streamwise vorticity, the half-velocity widths, the turbulence kinetic energy and the local shear in the mean streamwise velocity were obtained from the measured data. It was found that near field mixing in the equilateral triangular jet is faster than in the isosceles triangular and round jets. The mean streamwise vorticity field was found to be dominated by counter-rotating pairs of vortices, which influenced mixing and entrainment in the isosceles triangular jet. The one-dimensional energy spectra results indicated the presence of coherent structures in the near field of all three jets and that the equilateral triangular jet was more energetic than the isosceles triangular and round jets.     

Turbulent structure during transition to self-similarity in a round jet

 The developing turbulent region of a round jet was investigated using an improved implementation of digital particle image velocimetry (DPIV). The two-dimensional flow field in planes normal and parallel to the axial velocity was measured at locations between 15 and 30 diameters downstream, for two Reynolds numbers of 5500 and 16,000. The study consisted of instantaneous snapshots of the velocity and vorticity fields as well as measurements of velocity correlations up to third order. In this regime, the Reynolds number had a significant effect on both the instantaneous flow structure and the profiles of mean velocity across the jet. Coherent streamwise structures were present in the jet core for the lower Reynolds number. Additional structures whose evolution was governed by time scales two orders of magnitude larger than the convective scale inside the jet were observed in the entrainment field. The velocity correlations provided further support for the validity of DPIV turbulence measurements. The data was consistent with the equations of motion and momentum was conserved. DPIV measurements of turbulent kinetic energy components agreed with the hot-wire measurements of previous studies. Received: 27 November 1996/Accepted: 14 July 1997     

The development of an axisymmetric curved turbulent wall jet

An experimental study has been carried out of the low speed Coanda wall jet with both streamwise and axisymmetric curvature. A single component laser Doppler technique was used, and by taking several orientations at a given point, values of the three mean velocities and five of the six Reynolds stresses were obtained. The lateral divergence and convex streamwise curvature both enhanced the turbulence in the outer part of the jet compared with a plane two-dimensional wall jet. The inner layer exhibited a large separation of the positions of maximum velocity and zero shear stress. It was found that the streamwise mean velocity profile became established very rapidly downstream of the slot exit. The profile appeared fairly similar at later downstream positions, but the mean radial velocity and turbulence parameters showed the expected nonself preservation of the flow. Removal of the streamwise curvature resulted in a general return of the jet conditions toward those expected of a plane wall jet. The range and accuracy of the data may be used for developing turbulence models and computational techniques for this type of flow.     

Large and small scales in a turbulent orifice round jet: Reynolds number effects and departures from isotropy

An experimental investigation of the near field of a turbulent orifice jet is performed using high resolution Particle Image Velocimetry, aiming to highlight effects on the flow field due to changes in Reynolds number. The attention is focused onto departures from isotropy for large and small scales, by considering statistics of mean square velocity and velocity derivatives and specifically the non-dimensional ratios of such quantities. The results compare well with available literature data and pointed out that the effects of Reynolds number on large scales are usually small and limited to a region ranging less than seven-ten diameters from the jet outlet. For small scales, such Reynolds number dependence is extended up to ten-fifteen diameters. Farther from the jet exit, Reynolds number dependence almost disappears and all data approach similar asymptotic behaviors. On the other hand, velocity and some velocity derivative statistics clearly show that neither large nor small scale statistics strictly follow the isotropy condition; nonetheless, differences from that condition are limited to a factor which is almost constant in the whole measured field. In order to provide a link between such large and small scale departures from isotropy, a relation among mean square velocity ratios and mean square derivative ratios is proposed and proved to be well verified in the measured region and interval of Reynolds numbers. This relation allows deriving small scale derivative ratios, which are difficult to measure experimentally or to obtain numerically, due to high resolution requirements, from large scale velocity ratios, which are achieved much easier.     

Evaluation of stochastic particle dispersion modeling in turbulent round jets

ODT (one-dimensional turbulence) simulations of particle-carrier gas interactions are performed in the jet flow configuration. Particles with different diameters are injected onto the centerline of a turbulent air jet. The particles are passive and do not impact the fluid phase. Their radial dispersion and axial velocities are obtained as functions of axial position. The time and length scales of the jet are varied through control of the jet exit velocity and nozzle diameter. Dispersion data at long times of flight for the nozzle diameter (7 mm), particle diameters (60 and 90 µm), and Reynolds numbers (10, 000–30, 000) are analyzed to obtain the Lagrangian particle dispersivity. Flow statistics of the ODT particle model are compared to experimental measurements. It is shown that the particle tracking method is capable of yielding Lagrangian prediction of the dispersive transport of particles in a round jet. In this paper, three particle-eddy interaction models (Type-I, -C, and -IC) are presented to examine the details of particle dispersion and particle-eddy interaction in jet flow.     

PIV measurements of a turbulent jet issuing from round sharp-edged plate

An experimental investigation is presented of a turbulent jet issuing from a round sharp-edged orifice plate (OP) into effectively unbounded surroundings. Planar measurements of velocity were conducted using Particle Image Velocimetry (PIV) in the near and transition regions. The Reynolds number, based on the jet initial diameter and velocity, is approximately 72,000. The instantaneous and mean velocities, Reynolds normal and shear stresses were obtained. The centerline velocity decay and the half-velocity radius were derived from the mean velocity. It is revealed that primary coherent structures occur in the near field of the OP jet and that they are typically distributed asymmetrically with respect to the nozzle axis. Comparison of the present PIV and previous hot-wire measurements for the OP jet suggests that high initial turbulence intensity leads to reduced rates of decay and spread of the mean flow field and moreover a lower rate of variation of the turbulence intensity. Results also show that self-similarity of the mean flow is well established from the transition region while the turbulent statistics are far from self-similar within the measured range to 16 diameters.     

Investigation on time scales in a low Reynolds number jet flow using particle-tracking velocimetry

The near flow field of an axisymmetric water jet at Reynolds numbers between 2000 and 5000 is investigated using Particle-Tracking Velocimetry. Measurements are taken in the longitudinal section (along the mean flow) and in cross-sections (orthogonal to the mean flow). From the former, correlation coefficients of the two in-plane velocity components in a Lagrangian framework are obtained: thus Lagrangian integral scales can be computed. Those of the streamwise velocity (axial) component increase on moving away from the centreline, whereas the opposite happens for the vertical velocity (radial) component: integral time scales of the two components are almost equal at the interface between jet and ambient fluids. On the other hand, integral scales are almost constant or increase slightly with the axial direction. In cross-sections, fluid ejection and injection from the jet centreline are observed to be connected to counter-rotating vortices (mushroom): their number and size change with Reynolds number in agreement with results from other authors. The maximum ejection velocity (orthogonal to the mean jet flow), at 3 nozzle diameters downstream of the outlet, is found to be one half of the mean outlet velocity.     

Velocity statistics along the stagnation line of an axi- symmetric stagnating turbulent flow

 Velocity statistics along the stagnation line of an axi-symmetric wall stagnating turbulent flow are studied experimentally. A low turbulence, uniform air flow from a nozzle type air supply with an exit diameter of 50 mm stagnates at a wall located 50 mm downstream. A flow velocity is set to 3 m/s, 10 mm downstream from the exit of the air supply. Instantaneous values of streamwise and radial velocities are measured by laser-Doppler velocimetry. The turbulence level in the air flow is changed by use of turbulence generator. When the turbulence generator is not installed in the air supply, the mean velocity profile in the streamwise direction fits well with that of a laminar viscous flow with the rms value of velocity fluctuations low near the wall. With the turbulence generator installed, a significant turbulence structure appears near the wall. When the wall is approached, the rms value of velocity fluctuations in the streamwise direction decreases monotonically while the profile of the rms value in the radial direction reaches a maximum near the wall. The increase in the rms value of velocity fluctuations in the radial direction near the wall is attributed to the bi-modal histogram of the fluctuating velocity in the radial direction. Near the wall, the instantaneous stagnation streamline fluctuates and the probability of the mean location of the stagnation point reaches a maximum not at the stagnation line but on a circle around the stagnation line, resulting in the bi-modal histogram. Turbulence statistics, the rms value of velocity fluctuation and the turbulent kinetic energy, can be normalized successfully by similarity parameters based on the strain rate and the reference turbulent kinetic energy introduced by Champion and Libby. Received: 7 April 1995/Accepted: 27 September 1996     

Direct Numerical Simulation of Subsonic Round Turbulent Jet

Direct numerical simulation(DNS) of spatially developing round turbulent jet flow with Reynolds number 4,700 was carried out. Over 20 million grid points were used in this simulation. Fully compressible three-dimensional Navier–Stokes equations were solved. High order explicit spatial difference schemes and Runge–Kutta time integration scheme were used to calculate derivatives and time marching, respectively. Non-reflecting boundary conditions and exit zone techniques were adopted. Some refined computational grids were used in order to capture the smallest turbulent structures near the centerline of the jet. Low level disturbance were imposed on the jet inflow velocity to trigger the developing of turbulence. Turbulent statistics such as mean velocity, Reynolds stresses, third order velocity moments were obtained and compared with experimental data. One-dimensional velocity autospectra was also calculated. The inertial region where the spectra decays according to the k^− 5/3 was observed. The quantitative profiles of mean velocity and all of the third order velocity moments which were difficult to measure via experimental techniques were presented here in detail. The jet flow was proven to be close to fully self-similar around 19 jet diameters downstream of jet exit. The statistic data and revealed flow feature obtained in this paper can provide valuable reference for round turbulent jet research.