激励小尺度模式在湍流圆管射流中的应用

采用非涡黏性的激励小尺度（Stimulated Small Scale)模式对空间发展的轴对称湍流圆管射流进行了大涡模拟。以雷诺数为10000的流动为例,考证了激励小尺度模式在自由剪切流模拟中的可行性,描述了湍流强度、雷诺应力和湍流耗散量的变化,同时与标准的Smagorinsky涡黏性模式的计算结果进行了比较。数值结果显示,激励小尺度模式能够更为合理地描述湍流的耗散特性和能量传输特性,从而较为准确地展示出空间发展射流中由于流动不稳定而出现的旋涡产生、发展、破碎及合并等过程。     

方形管道中湍流水流某些特征的激光多普勒法测量及其与圆形管道的对比

用激光多普勒测速法测量了雷诺数从4.76×10~3到4.76×10~4的矩形管道中的湍流水流。其结果与相同雷诺数时的圆形管道的流动就某些流动特征进行了对比。引人注目的是,在矩形管道内速度分布更“平坦”,且随雷诺数的增加而强烈地趋向更加均匀,湍流强度也比圆管流动的大得多。在横截面的中线上,对四种不同雷诺数所得的速度分布数据与普适的湍流速度外部定律有较好的符合。为了保持湍流的功率谱以及减小频带加宽,在测量仪器方面使用了某些技术。由Berman和Dunning所给出的有关速度梯度加宽的效应的式子也从本测量中得到了实验数据的支持。     

基于全局动态Vreman模型的复杂非定常湍流的大涡模拟

在湍流数值模拟方法中,大涡模拟方法可以提供丰富的大涡旋信息,已逐渐成为复杂湍流问题数值研究的重要方法。而大涡模拟中,最重要的一环是尽量准确地构建能反映流场物理本质特征的亚格子应力模型。基于该思想,将一种新型的大涡模拟亚格子应力模型-Vreman亚格子应力模型用于高雷诺数三维后台阶流动的求解,计算结果与实验结果进行对比分析结果较吻合,验证了该模型的可靠性。这是对该模型用于无任何均匀流动方向的高雷诺数复杂湍流非定常流动的首次检验,计算结果优于基于传统的Smagorinsky涡粘性的动态亚格子模型。     

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

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

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

湍流涡致振动频率特性

用数值模拟方法对固定圆柱湍流涡脱落频率与弹性圆柱湍流涡致振动频率特性进行了研究,湍流计算模型采用标准κ-ε模型,压力泊松方程提法基于非交错网格系统.研究结果表明:固定圆柱湍流绕流涡脱落频率基本不随雷诺数而变,对于同一固有频率弹性圆柱,涡振频率基本不随雷诺数而变;对于某一固定雷诺数流动涡振频率在一定范围内与系统固有频率有关.     

谱消去黏性谱元法大涡模拟

引入一种新的利用谱元法进行湍流大涡模拟的方法: 谱消去黏性法.谱消去黏性法原是为了解决双曲型问题谱逼近的稳定性而引进的,最近人们发现它还可用于湍流大涡模拟. 与其它大涡模拟方法相比,这种方法几乎不必修改原代码便可在标准的谱元法中实现,而且几乎不增加计算量. 文章使用谱元法结合谱消去黏性法对雷诺数12\,000时的三维驱动方腔流进行湍流大涡模拟,并提供了模拟的初步数值结果及其统计分析,湍流统计特性表明得到的结果与已知的实验和直接数值模拟结果有较好的一致性. 另外,还考察了不同的谱消去黏性参数对稳定性和模拟结果的影响.      

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

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

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

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

管内湍流旋流的数值计算

本文用有限差分法对直管内的湍流旋流进行了数值模拟。计算中采用Ｂｏｕｓｓｉｎｅｓｑ湍流涡粘性假设的基本思想和Ｋ－ε双方程模型来求解雷诺应力各分量。为了反映旋流中湍流转输的非均匀性和各向异性特征，对雷诺应力各分量及与之相主尖的各湍流粘性系数分别进行计算。计算结果表明该模型能较好地反映直管内湍流旋流的流动结构。     

The spanwise spectra in wall-bounded turbulence

The pre-multiplied spanwise energy spectra of streamwise velocity fluctuations are investigated in this paper. Two distinct spectral peaks in the spanwise spectra are observed in low-Reynolds-number wall-bounded turbulence.The spectra are calculated from direct numerical simulation(DNS) of turbulent channel flows and zero-pressure-gradient boundary layer flows. These two peaks locate in the nearwall and outer regions and are referred to as the inner peak and the outer peak, respectively. This result implies that the streamwise velocity fluctuations can be separated into large and small scales in the spanwise direction even though the friction Reynolds number Reτ can be as low as 1000. The properties of the inner and outer peaks in the spanwise spectra are analyzed. The locations of the inner peak are invariant over a range of Reynolds numbers. However, the locations of the outer peak are associated with the Reynolds number,which are much higher than those of the outer peak of the pre-multiplied streamwise energy spectra of the streamwise velocity.     

Turbulence manipulation in air–water flows on a stepped chute: An experimental study

In a stepped channel operating with large flow rates, the flow skims over the pseudo-bottom formed by the step edges as a coherent stream. Intense three-dimensional recirculation is maintained by shear stress transmission from the mainstream to the step cavities, while significant free-surface aeration takes place. The interactions between free-surface aeration and cavity recirculation are investigated herein with seven step cavity configurations. The experiments were conducted in a large stepped channel operating at large Reynolds numbers. For some experiments, triangular vanes, or longitudinal ribs, were placed across the step cavities to manipulate the flow turbulence to enhance the interactions between the mainstream flow and the cavity recirculation region. The results showed a strong influence of the vanes on the air–water flow properties in both free-stream and cavity flows. The findings demonstrate some passive turbulence manipulation in highly turbulent air–water flows.     

Turbulent thermal diffusion in a multi-fan turbulence generator with imposed mean temperature gradient

We studied experimentally the effect of turbulent thermal diffusion in a multi-fan turbulence generator which produces a nearly homogeneous and isotropic flow with a small mean velocity. Using particle image velocimetry and image processing techniques, we showed that in a turbulent flow with an imposed mean vertical temperature gradient (stably stratified flow) particles accumulate in the regions with the mean temperature minimum. These experiments detected the effect of turbulent thermal diffusion in a multi-fan turbulence generator for relatively high Reynolds numbers. The experimental results are in compliance with the results of the previous experimental studies of turbulent thermal diffusion in oscillating grid turbulence (Buchholz et al. 2004; Eidelman et al. 2004). We demonstrated that the turbulent thermal diffusion is an universal phenomenon. It occurs independently of the method of turbulence generation, and the qualitative behavior of particle spatial distribution in these very different turbulent flows is similar. Competition between turbulent fluxes caused by turbulent thermal diffusion and turbulent diffusion determines the formation of particle inhomogeneities.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Large-eddy structures of turbulent swirling flows and methane-air swirling diffusion combustion

Turbulent swirling flows and methane-air swirling diffusion combustion are studied by large-eddy simulation (LES) using a Smagorinsky-Lilly subgrid scale turbulence model and a second-order moment (SOM) SGS combustion model, and also by RANS modeling using the Reynolds Stress equation model with the IPCM+wall and IPCM pressure-strain models and SOM combustion model. The LES statistical results for swirling flows give good agreement with the experimental results, indicating that the adopted subgrid-scale turbulence model is suitable for swirling flows. The LES instantaneous results show the complex vortex shedding pattern in swirling flows. The initially formed large vortex structures soon break up in swirling flows. The LES statistical results of combustion modeling are near the experimental results and are as good as the RANS-SOM modeling results. The LES results show that the size and range of large vortex structures in swirling combustion are different from those of isothermal swirling flows, and the chemical reaction is intensified by the large-eddy vortex structures. The project supported by the Special Funds for Major State Basic Research (G-1999-0222-07). The English text was polished by Keren Wang.     

Investigation of downstream and sideline subsonic jet noise using Large Eddy Simulation

The sound fields radiated by Mach number 0.6 and 0.9, circular jets with Reynolds numbers varying from 1.7×10³ to 4×10⁵ are investigated using Large Eddy Simulations. As the Reynolds number decreases, the properties of the sound radiation do not change significantly in the downstream direction, whereas they are modified in the sideline direction. At low Reynolds numbers, for large angles downstream from the jet axis, the acoustic levels are indeed remarkably lower and a large high-frequency part of the sound spectra vanishes. For all Reynolds numbers, the downstream and the sideline sound spectra both appear to scale in frequency with the Strouhal number. However their peak amplitudes vary following two different velocity exponents according to the radiation direction. The present observations suggest the presence of two sound sources: a Reynolds number-dependent source, predominant for large radiation angles, connected to the randomly-developing turbulence, and a deterministic source, radiating downstream, related to a mechanism intrinsic to the jet geometry, which is still to be comprehensively described. This view agrees well with the experimental results displaying two distinguishable components in turbulent mixing noise [1, 2].     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.     

An experimental study of the effect of free-stream turbulence on a trailing vortex

A study of turbulence evolution and spectra within and just outside the core of a trailing vortex is performed. The vortex is generated by a vortex generator consisting of four blades positioned orthogonally to each other with the same angle of attack and placed in a low-speed wind tunnel. A grid is placed upstream of the vortex generator to produce free-stream turbulence, which wraps around and interacts with the columnar vortex. Instantaneous measurements of the three velocity components are obtained using a miniature four-sensor hot-wire probe. The study focuses on the distribution of turbulence energy and Reynolds stress among the different spectral components of the flow at different positions across the vortex core and different axial positions along the tunnel. The effect of background grid turbulence on the spectral energy distribution of the vortex is examined in comparison to the vortex alone.     

Turbulence,dynamic similarity and scale effects in high-velocity free-surface flows above a stepped chute

In high-velocity free-surface flows, air entrainment is common through the interface, and intense interactions take place between turbulent structures and entrained bubbles. Two-phase flow properties were measured herein in high-velocity open channel flows above a stepped chute. Detailed turbulence measurements were conducted in a large-size facility, and a comparative analysis was applied to test the validity of the Froude and Reynolds similarities. The results showed consistently that the Froude similitude was not satisfied using a 2:1 geometric scaling ratio. Lesser number of entrained bubbles and comparatively greater bubble sizes were observed at the smaller Reynolds numbers, as well as lower turbulence levels and larger turbulent length and time scales. The results implied that small-size models did underestimate the rate of energy dissipation and the aeration efficiency of prototype stepped spillways for similar flow conditions. Similarly a Reynolds similitude was tested. The results showed also some significant scale effects. However a number of self-similar relationships remained invariant under changes of scale and confirmed the analysis of Chanson and Carosi (Exp Fluids 42:385-401, 2007). The finding is significant because self-similarity may provide a picture general enough to be used to characterise the air–water flow field in large prototype channels.     

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.