DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

旋转圆柱绕流流场特性分析

对雷诺数Re = 20000 ~ 90000、相对转速ɑ = 0 ~ 0.72的旋转圆柱后方流场进行了实验测量, 分析了旋转圆柱后方不同剖面处的速度分布规律和湍流度分布规律. 采用LES方法对旋转圆柱绕流问题进行了数值模拟, 分析旋转圆柱周围流场特性和自由剪切层变化规律, 最后通过理论模型对流场变化进行分析, 得出如下结论: 当圆柱逆时针旋转时, 同一雷诺数下随着相对转速的增加, 旋转圆柱尾迹区域下方速度突变处的位置随着相对转速的增加而上移, 而上方速度突变处的位置不变, 雷诺数的增加使旋转圆柱尾迹区域下方速度突变处位置有小幅度的下移. 通过数值模拟发现, 圆柱旋转之后, 圆柱后方下侧涡的位置明显上移, 且幅度较大. 下方的自由剪切层有明显的上移, 上方的自由剪切层位置变化较小. 最后通过理论分析发现, 圆柱后侧下方涡位置的上移对圆柱升力影响十分显著, 在高雷诺数、低相对转速的条件下, 旋转圆柱后侧下方涡位置的改变对旋转圆柱的升力、尾流区自由剪切层的变化起到了重要的影响.      

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.     

单一水平轴风电机组尾迹的模拟方法与流动机理研究综述

为实现碳达峰、碳中和“3060”目标, 风能将在我国能源体系发挥重要作用. 风力机尾迹是影响风电性能和度电成本的关键因素, 需在风力机布置和控制设计中充分考虑. 本文首先介绍风力机尾迹的数值模拟方法, 包括解析模型、低阶模型、大涡模拟和来流湍流生成方法. 解析模型和低阶模型可快速计算风力机尾迹, 但依赖于模型参数, 且不能或不能准确预测尾迹湍流特性. 结合风力机参数化模型的大涡模拟可准确预测尾迹蜿蜒等湍流特征, 是流动机理研究的有力工具, 可为发展快速预测模型提供数据和理论支撑. 接着, 本文介绍了叶尖涡、中心涡和尾迹蜿蜒并讨论其产生机理. 对于湍流来流, 叶尖涡主要存在于近尾迹. 蜿蜒是远尾迹的主要特征, 影响下游风力机的来流特征. 尾迹蜿蜒的产生有两种机制: 来流大尺度涡和剪切层失稳. 数值和观测结果显示两种机制共同存在. 机舱和中心涡对尾迹蜿蜒有重要影响. 采用叶片和机舱的致动面模型可准确预测尾迹蜿蜒. 研究显示不同风力机尾迹间的湍流特征存在相似性, 为发展尾迹湍流的快速预测模型提供了理论依据. 当前研究多关注平坦地形上的风力机尾迹, 复杂地形和海洋环境下的大气湍流和风力机尾迹的机理复杂, 现有工程模型无法准确预测, 有待深入研究.      

均匀、振荡及组合流绕平板流动的尾迹研究

本文利用离散涡模型及改进的新生涡产生机制对三种不同来流绕平板的近尾迹进行数值研究。计算讨论了定常流中平板绕流流动的总体特性和近尾迹流场;对于简谐振荡来流,相应于K_c=2.0、4.0 和10.0 分别得到两种不同的尾迹形态。给出了小 K_c 数平板尾迹涡配对、运动的新模式而相应的阻力、惯性力系数计算比以前涡模拟结果更接近于 U 型管实验结果。对于流向组合来流本文模拟了涡锁定及其动力特性并于实验相符,给出了流向扰动对平板绕流流动的影响。     

高亚临界雷诺数圆柱绕流的尺度自适应模拟

借助γ-Re_θ转捩模型,实现了高亚临界雷诺数（Re=1.4×105）下圆柱层流分离流动的尺度自适应模拟.统计平均结果表明数值计算和实验测量较为接近,尤其在圆柱后半段的分离区中,压力系数和实验符合得很好,误差主要源于分离点预测的不准确. 瞬态流动则显示,层流分离的剪切层中出现了展向不稳定,且在向下游的输运过程中不断增强,最后转捩为完全湍流. 在湍流分离模拟中,由于缺乏剪切层失稳的非定常性,SST-SAS 模型的尺度分辨能力变弱,因此在分离区以及下游尾迹中求解出的湍流尺度要明显较层流分离时大.      

SCALE ADAPTIVE SIMULATION OF FLOWS AROUND A CIRCULAR CYLINDER AT HIGH SUB-CRITICAL REYNOLDS NUMBER

借助γ-Re_θ转捩模型,实现了高亚临界雷诺数（Re=1.4×10⁵）下圆柱层流分离流动的尺度自适应模拟.统计平均结果表明数值计算和实验测量较为接近,尤其在圆柱后半段的分离区中,压力系数和实验符合得很好,误差主要源于分离点预测的不准确. 瞬态流动则显示,层流分离的剪切层中出现了展向不稳定,且在向下游的输运过程中不断增强,最后转捩为完全湍流. 在湍流分离模拟中,由于缺乏剪切层失稳的非定常性,SST-SAS 模型的尺度分辨能力变弱,因此在分离区以及下游尾迹中求解出的湍流尺度要明显较层流分离时大.     

湍流涡致振动频率特性

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

均匀旋转对圆柱水动力及流动结构的影响

基于格子Boltzmann方法 (LBM)对均匀旋转控制下的低雷诺数(Re=100)圆柱绕流问题进行了数值模拟,得到了转速比从0～10变化下,旋转控制对圆柱水动力及流动结构的影响规律.使用动态模态分解(DMD)对流场特征进行提取,并分析了施加旋转控制之后转速比对流场不同模态和增长率的影响.结果表明,随着转速比增大,圆柱下游流动结构依次呈现出卡门涡街、剪切层、反向剪切层、单侧涡和附着涡5种结构;阻力系数时均值先减小,随后在转速进入单侧涡区间后增大,升力系数与力矩系数的时均值均单调增加,同时,在出现涡脱落的两个转速区间内,水动力出现了明显的波动,且二次失稳时波动幅度更大. DMD的结果表明,圆柱下游的流动结构主要受圆柱壁面的旋转影响而发生改变并产生全新流动模态;旋转会对流动稳定性产生影响:在未充分发展阶段,旋转对流动稳定性的影响不显著,而在充分发展后,各转速下的流场不稳定模态数均远少于未充分发展阶段,随着转速比的增大,流动稳定性会产生不同程度的增强或减弱,且无涡脱落时的稳定性高于有涡脱落时,因此,通过旋转控制抑制尾涡脱落可以有效增强流动的稳定性.     

Population,characteristics and kinematics of vortices in a confined rectangular jet with a co-flow

Vortex behavior and characteristics in a confined rectangular jet with a co-flow were examined using vortex swirling strength as a defining characteristic. On the left side of the jet, the positively (counterclockwise) rotating vortices are dominant, while negatively rotating vortices are dominant on the right side of the jet. The characteristics of vortices, such as population density, average size and strength, and deviation velocity, were calculated and analyzed in both the cross-stream direction and the streamwise direction. In the near-field of the jet, the population density, average size and strength of the dominant direction vortices show high values on both sides of the center stream with a small number of counter-rotating vortices produced in the small wake regions close to jet outlet. As the flow develops, the wake regions disappear, these count-rotating vortices also disappear, and the population of the dominant direction vortices increase and spread in the jet. The mean size and strength of the vortices decrease monotonically with streamwise coordinate. The signs of vortex deviation velocity indicate the vortices transfer low momentum to high-velocity region and high momentum to the low velocity region. The developing trends of these characteristics were also identified by tracing vortices using time-resolved particle image velocimetry data. Both the mean tracked vortex strength and size decrease with increasing downstream distance overall. At the locations of the left peak of turbulent kinetic energy, the two-point spatial cross-correlation of swirling strength with velocity fluctuation and concentration fluctuation were calculated. All the correlation fields contain one positively correlated region and one negatively correlated region although the orientations of the correlation fields varied, due to the flow transitioning from wake, to jet, to channel flow. Finally, linear stochastic estimation was used to calculate conditional structures. The large-scale structures in the velocity field revealed by linear stochastic estimation are spindle-shaped with a titling stream-wise major axis.     

Multisensor Hot Wire Vorticity Probe Measurements of the Formation Field of Two Corotating Vortices

Experimental evidence is reported, regarding the formation of a pair of co-rotating tip vortices by a split wing configuration, consisting of two half wings at equal and opposite angles of attack. Simultaneous measurements of the three-dimensional vector fields of velocity and vorticity were conducted on a cross plane at a downstream distance corresponding to 0.3 cord lengths (near wake), using an in-house constructed 12-sensor hot wire anemometry vorticity probe. The probe consists of three closely separated orthogonal 4-wire velocity sensor arrays, measuring simultaneously the three-dimensional velocity vector at three closely spaced locations on a cross plane of the flow filed. This configuration makes possible the estimation of spatial velocity derivatives by means of a forward difference scheme of first order accuracy. Velocity measurements obtained with an X-wire are also presented for comparison. In this near wake location, the flow field is dictated by the pressure distribution established by the flow around the wings, mobilizing large masses of air and leading to the roll up of fluid sheets. Fluid streams penetrating between the wings collide, creating on the cross plane flow a stagnation point and an “impermeable” line joining the two vortex centres. Along this line fluid is directed towards the two vortices, expanding their cores and increasing their separation distance. This feeding process generates a dipole of opposite sign streamwise mean vorticity within each vortex. The rotational flow within the vortices obligates an adverse streamwise pressure gradient leading to a significant streamwise velocity deficit characterizing the vortices. The turbulent flow field is the result of temporal changes in the intensity of the vortex formation and changes in the position of the cores (wandering).     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

Cross-Stream Vorticity Field Induced by Streamwise Vortices in Unbounded and Wall-Bounded Shear Flows

This computational study examines the unsteady cross-stream vorticity structures that form when one or more streamwise vortices are immersed in homogeneous and boundary-layer shear flows. A quasi-two-dimensional limit is considered in which the velocity and vorticity fields, while still possessing three nonzero components, have vanishing gradient in the streamwise direction. This idealization is suitable to applications such as streamwise vortices that occur along a ship hull or airplane fuselage and it can be used as an idealized representation of the quasi-streamwise vortices in the near-wall region of a turbulent boundary layer. In this quasi-two-dimensional idealization, the streamwise velocity has no effect on the cross-stream velocity associated with the vortex. However, the vortex acts to modify the cross-stream vorticity component, resulting in regions of the flow with strong deviations in streamwise velocity. This paper examines the complex structures that form as the cross-stream vorticity field is wrapped up by the vortex and the effect of these structures on the streamwise velocity field, first for vortices immersed in homogeneous shear flow and then for vortices immersed in a boundary layer along a flat wall. Received 2 January 2002 and accepted 13 August 2002 Published online 3 December 2002 RID="*" ID="*" This project was supported by the Office of Naval Research under Grant Number N00014-01-1-0015. Dr. Thomas Swain is the program manager. Communicated by T.B. Gatski     

圆柱尾迹涡的三维演化及结构特征

应用无粘涡丝运动学理论和局部诱导近似(LIA)方法,以Lagrange观点跟踪涡丝在背景流场中运动,用数值方法研究了中等Re数(≈10~3)下圆柱分离尾迹中Kármán涡和涡辫区涡丝的三维演化的机制和动力学过程,及其结构特征。背景流场考虑为尾迹时间平均速度流场和Kármán涡街流场。初始展向小扰动为指数形式和谐波形式。结果指出:Kármán涡和涡辫区中的涡丝具有展向不稳定性,形成流向涡量。在尾迹的初期输运过程中,表现出有序的大、小尺度涡结构。并进一步分析了其产生的机理。     

Proper orthogonal decomposition analysis of coherent motions in a turbulent annular jet

A three-dimensional incompressible annular jet is simulated by the large eddy simulation(LES) method at a Reynolds number Re = 8 500. The time-averaged velocity field shows an asymmetric wake behind the central bluff-body although the flow geometry is symmetric. The proper orthogonal decomposition(POD) analysis of the velocity fluctuation vectors is conducted to study the flow dynamics of the wake flow.The distribution of turbulent kinetic energy across the three-dimensional POD modes shows that the first four eigenmodes each capture more than 1% of the turbulent kinetic energy, and hence their impact on the wake dynamics is studied. The results demonstrate that the asymmetric mean flow in the near-field of the annular jet is related to the first two POD modes which correspond to a radial shift of the stagnation point. The modes 3 and 4 involve the stretching or squeezing effects of the recirculation region in the radial direction. In addition, the spatial structure of these four POD eigenmodes also shows the counter-rotating vortices in the streamwise direction downstream of the flow reversal region.     

TR-PIV measurement of the wake behind a grooved cylinder at low Reynolds number

A comparative study of the wakes behind cylinders with grooved and smooth surfaces was performed with a view to understand the wake characteristics associated with the adult Saguaro cacti. A low-speed recirculation water channel was established for the experiment; the Reynolds number, based on the free-stream velocity and cylinder diameter (D), was kept at Re_D=1500. State-of-the-art time-resolved particle image velocimetry (TR-PIV) was employed to measure a total of 20 480 realizations of the wake field at a frame rate of 250 Hz, enabling a comprehensive view of the time- and phase-averaged wake pattern. In comparison to the wake behind the smooth cylinder, the length of the recirculation zone behind the grooved cylinder was extended by nearly 18.2%, yet the longitudinal velocity fluctuation intensity was considerably weakened. A global view of the peaked spectrum of the longitudinal velocity component revealed that the intermediate region for the grooved cylinder, which approximately corresponds to the transition region where the shear layer vortices interact, merge and shed before the formation of the Karman-like vortex street, was much wider than that for the smooth one. The unsteady events near St=0.3-0.4 were detected in the intermediate region behind the grooved cylinder, but no such events were found in the smooth cylinder system. Although the formation of the Karman-like vortex street was delayed by about 0.6D downstream for the grooved cylinder, no prominent difference in the vortex street region was found in the far wake for both cylinders. The Proper Orthogonal Decomposition (POD) method was used extensively to decompose the vector and swirling strength fields, which gave a close-up view of the vortices in the near wake. The first two POD modes of the swirling strength clarified the spatio-temporal characteristics of the shear layer vortices behind the grooved cylinder. The small-scale vortices superimposed on the shear layers behind the grooved cylinder were found to be generated and convected downstream in the same phase, which would significantly reduce the fluctuating force on the cylinder surface.     

Topology of the vorticity field in three-dimensional shear layers and wakes

An experimental and numerical study of the three-dimensional transition of plane wakes and shear layers behind a flat plate is presented. Flow visualization techniques are used to monitor the response of laminar flows at moderate Reynolds numbers (≈100) to perturbations periodically distributed along the span. In this way, the formation and evolution of streamwise vortex tubes and their interaction with the spanwise vortices are analyzed. The flow was studied numerically by means of three-dimensional inviscid vortex dynamics. Assuming periodicity in the spanwise and the streamwise direction, we discretize the vorticity field into two layers of vortex filaments with finite core diameter. Comparison between experiment and visualization indicates that important features of the three-dimensional evolution can be reproduced by inviscid vortex dynamics. Vortex stretching in the strain field of the spanwise rollers appears to be the primary mechanism for the three-dimensional transition in this type of flows.     

Time resolved flow-field measurements of a turbulent mixing layer over a rectangular cavity

High Reynolds number, low Mach number, turbulent shear flow past a rectangular, shallow cavity has been experimentally investigated with the use of dual-camera cinematographic particle image velocimetry (CPIV). The CPIV had a 3 kHz sampling rate, which was sufficient to monitor the time evolution of large-scale vortices as they formed, evolved downstream and impinged on the downstream cavity wall. The time-averaged flow properties (velocity and vorticity fields, streamwise velocity profiles and momentum and vorticity thickness) were in agreement with previous cavity flow studies under similar operating conditions. The time-resolved results show that the separated shear layer quickly rolled-up and formed eddies immediately downstream of the separation point. The vortices convect downstream at approximately half the free-stream speed. Vorticity strength intermittency as the structures approach the downstream edge suggests an increase in the three-dimensionality of the flow. Time-resolved correlations reveal that the in-plane coherence of the vortices decays within 2–3 structure diameters, and quasi-periodic flow features are present with a vortex passage frequency of ~1 kHz. The power spectra of the vertical velocity fluctuations within the shear layer revealed a peak at a non-dimensional frequency corresponding to that predicted using linear, inviscid instability theory.