Reynolds数22 000的孤立方柱绕流的大涡模拟

方柱绕流是典型的钝体绕流问题,蕴含了丰富的流体力学现象,对这类流动的准确预测面临着诸多挑战.采用自主发展的大涡模拟程序,对来流Mach数M=0.3,Reynolds数Re_D=22 000的绕孤立方柱流动进行了细致模拟,亚格子模型使用动力涡黏模型.对计算结果的分析表明,大涡模拟所得的平均流场及Reynolds应力分布与已有实验数据和直接数值模拟结果均吻合较好,验证了预测结果的可靠性;在此基础上对瞬态流场进行了研究,展示了计算条件下方柱绕流分离转捩及尾迹区旋涡交替脱落形成Karman涡街的全过程,为更细致的流动机理探索奠定了基础.      

PISO算法在求解方柱绕流非定常尾迹中的应用

应用PISO算法及非交错斜交网格技术,编制了求解二维定常、非定常流动的通用程序。计算分析了定常不可压缩Couette流动以及水平通道内方柱绕流尾迹的非定常流动,并详细分析了不同时刻方柱后尾迹的旋涡结构及发展过程,得到了合理的结果。验证了非迭代PISO算法及非交错斜交网格在分析非定常流动中的有效性,为进一步开展非定常流动的数值分析奠定了基础．     

基于UGKS的过渡区方柱绕流数值研究

本文基于统一气体动理学格式(Unified Gas Kinetic Scheme,UGKS),对微尺度过渡区气体绕流方柱开展数值模拟研究,计算分析了Kn数对气体流动传热过程的影响规律。研究发现,随着Kn数增加,方柱壁面气体速度滑移和温度跳跃增大,壁面上压力、剪切力和热流也相应增大,方柱壁面-气体之间的换热得到强化;方柱对气体流动的阻碍作用减小,方柱前滞止区影响范围相对增大,方柱温度对柱后区域气体温度影响相对减小。     

激波沿H2／空气界面绕圆、方柱强化混合的数值研究

用迎风TVD格式求解二维、层流全N－S方程,对激波沿H2和空气界面绕圆、方柱流动及其诱导的剪切混合进行了数值模拟,得到了流场的压力和组分密度分布,计算结果表明：激波在H2传播得快,剪切层中出现吸涡和调节激波,卷吸涡与柱体撞击后,反射出一道激波,H2沿柱体表面向下游扩散,H2／空气接触面与柱体分离后,形状畸变并产生新的卷吸涡。H2分布表明：在办面上加圆柱或方柱,可有效地强化燃料混合,方柱的增强效果更明显此,在圆柱表面,H2、空气中激波均发生由RR向MR的转变,两Mach杆在下游相互透射,对于方柱,H2中激波中激波沿下表面传播几乎不受影响,空气中激波沿上表面发生Mach反射,其Mach杆和H2中绕射激波相互透射,柱体左侧最终形成一脱体激波,流场存在激波、卷吸涡、接触面向的相互作用,但波系结构相似。     

激光熔覆中同轴粉末流温度场的数值模拟

在激光同轴送粉熔覆中,由于激光与粉末流相互作用,粉末流整体温度分布直接影响激光熔覆的质量.基于非预混燃烧模型,将激光相处理为连续性介质,粉末颗粒相看作离散相物质,建立了激光作用下粉末流的质量、动量和能量方程.用Fluent软件进行了不同激光功率和粉末流速度条件下粉末流整体温度场数值模拟,讨论了各种参数对温度场分布的影响.为了验证该模型的准确度,利用CCD比色测温方法测量了粉末流整体温度场分布.结果表明,数值模拟与CCD检测结果具有良好的一致性,数值模拟结果对激光熔覆具有指导意义.     

高负荷对转压气机尾迹涡对叶顶泄漏流的影响

采用三维非定常数值模拟方法,针对高负荷对转压气机低压转子尾迹脱落涡对高压转子叶顶泄漏流的影响开展研究。研究发现:低压转子尾迹涡以"对涡"的形式在高压转子通道中输运,顺时针与逆时针旋向尾迹涡间隔分布,对叶顶泄漏流与主流交界面法向动量产生不同影响,进而使叶顶泄漏流与主流交界面形状呈现波浪形。设计工况下,尾迹涡会使叶顶泄漏流与主流交界面的位置偏向吸力面以及后移结尾激波位置进而后移二次泄漏流产生的位置,减小压力面与吸力面结尾激波入射点之间弦长区域叶顶泄漏流与轴向的夹角,提高相应弦长区域叶顶泄漏流的轴向速度,减小高压转子叶尖区的堵塞,降低二次泄漏流产生的损失,进而提升高压转子叶尖区的等熵效率。     

非定常激励抑制轴流叶栅分离的数值分析

本文对大转角压气机叶栅大攻角工况在来流尾迹非定常激励作用下的时空结构进行了数值分析.合理的非定常激励明显提升了叶栅的时均性能,与激励频率对应的涡结构得到强化,其它频率的杂涡被卷吸,无序的吸力面分离流结构变得有序.文中分析了周期性来流尾迹控制分离结构的作用机理,波涡共振促进涡的卷起和配对;湍动能强化动量交换.基于此,比较了两种来流尾迹模式下分离结构对来流周期性的响应问题.另外,针对管道风扇引起的远场噪声,在来流尾迹非定常激励作用下对其降低的可能性进行了设想.     

稀薄气体二维外部柱体绕流特性研究

本文采用直接模拟蒙特卡罗方法对稀薄气体二维外部柱体绕流问题进行了数值模拟。结果表明:外部绕流问题,在特定情况下会产生激波,激波的产生,不仅与气体的稀薄程度有关,还与来流马赫数有关。而气流与壁面之间的换热,随来流马赫数增加而增加,随气体稀薄程度增加而减小。     

方柱绕流大涡模拟

采用有限体/有限元混合格式、非结构网格和大涡模拟方法求解可压缩的N-S方程,对Re=22 000的方柱绕流进行数值模拟,并对不同的边界条件进行详细的分析比较.通过对以往研究经验的总结和利用精细的边界条件,使得采用二阶精度的数值格式和较稀疏的网格仍然得到了令人满意的计算结果,甚至优于以往采用密网格的模拟结果.     

U-RANS/PDF方法的并行计算研究

本文对非结构网格中U-RANS/PDF混合方法的并行计算进行了研究,通过对方柱绕流的数值模拟,讨论了并行程序的有效性和计算效率.非结构网格上U-RANS/PDF混合方法的并行计算,提高了PDF方法的计算效率和处理复杂几何结构的能力,使得PDF方法能够方便地应用于实际燃烧问题.     

低雷诺数下钝体三维尾迹中的涡量符号律

钝体是目前各种工程中广泛应用的一种结构.钝体绕流的尾迹涡动力学也是经典的流体力学研究对象之一.本文通过直接数值模拟,针对低雷诺数下各种钝体结构的不可压缩绕流,当形成三维尾迹时,研究具有特定符号的涡量分布特征.通过分析两类钝体结构,基本的直柱体和受到几何扰动的柱体,总结并得到了更为广泛适用的涡量符号律.通过对比并分析这两类钝体结构,结合理论证明的结果,进一步厘清了对产生涡量符号律的这两类钝体结构之间的内在物理关联,即引起自然失稳的小扰动在惯性力作用下产生的表面涡量只能向下游演化发展,而几何扰动则根据扰动位置,产生的表面涡量可以向扰动上游或下游演化发展.从而可以推测所有钝体结构尾迹中的各种型式的涡脱落模态,从涡量符号律的演化角度来看,实际上是一致的,都是起源于壁面产生特定符号组合规律的∏型涡.     

A study on drag reduction of a rotationally oscillating circular cylinder at low Reynolds number

The flow field around a rotationally oscillating circular cylinder in a uniform flow is studied by using a particle image velocimetry to understand the mechanism of drag reduction and the corresponding suppression of vortex shedding in the cylinder wake at low Reynolds number. Experiments are conducted on the flow around the circular cylinder under rotational oscillation at forcing Strouhal number 1, rotational amplitude 2 and Reynolds number 2,000. It is found from the flow measurement by PIV that the width of the wake is narrowed and the velocity fluctuations are reduced by the rotational oscillation of the cylinder, which results in the drag reduction rate of 30%. The mechanism of drag reduction is studied by phase-averaged PIV measurement, which indicates the formation of periodic small-scale vortices from both sides of the cylinder. It is found from the cross-correlation measurement between the velocity fluctuations that the large-scale structure of vortex shedding is almost removed in the cylinder wake, when the small-scale vortices are generated at the unstable frequency of shear layer by the influence of rotational oscillation.     

Viscous near-wall flow in a wake of circular cylinder at moderate Reynolds numbers

Here we present the results of experimental investigation of a cross flow around a circular cylinder mounted near the wall of a channel with rectangular cross section. The experiments were carried out in the range of Reynolds numbers corresponding to the transition to turbulence in a wake of the cylinder. Flow visualization and SIV-measurements of instantaneous velocity fields were carried out. Evolution of the flow pattern behind the cylinder and formation of the regular vortex structures were analyzed. It is shown that in the case of flow around the cylinder, there is no spiral motion of fluid from the side walls of the channel towards its symmetry plane, typical of the flow around a spanwise rib located on the channel wall. The laminar-turbulent transition in the wake of the cylinder is caused by the shear layer instability.     

High-definition PIV analysis on vortex shedding in the cylinder wake

A high-definition analysis based on flow topology is made on the vortex motions under natural and lock-on conditions in the near-wake region of a circular cylinder, where two-dimensional flow fields in the wake-transition regime are measured by a time-resolved PIV system. The Reynolds stress distributions are examined in view of the mean separation streamline and the trajectory ofthe vortex center. It is shown that, by the lock-on, the Reynolds stresses become stronger and their dispositions match well with the shortened wake bubble, indicating perfect synchronization of shedding to the oscillatory forcing flow in the near field, which causes increased lift and drag forces.     

Vortex generation in the cross-flow around a cylinder attached to an end-wall

A cylinder attached to an end-wall normal to its axis is a common feature of many practical flow systems, e.g. in turbo-machinery or when a bridge is supported by a pillar from the bed of a river. In this situation, the nominally two-dimensional boundary layer flow incident upon the cylinder develops strong three-dimensional features and a very pronounced vortex structure may arise in the upstream flow close to the wall. For the appropriate Reynolds number range, the upstream vortical structure is nominally steady and is commonly referred to as the “horseshoe vortex system”. In contrast, the flow downstream is unsteady and periodic over a wide range of Reynolds numbers and vortices aligned with the cylinder axis are shed at a regular frequency into the wake. The generation of both these vortex systems requires energy to be extracted from the incident flow with the result that the drag force on the cylinder is increased.This paper concentrates on the upstream region of the cylinder and discusses an investigation in which two-component Particle Image Velocimetry (PIV) has been used to visualise the flow behaviour for a circular cylinder on a plane end-wall. The use of PIV has enabled two orthogonal velocity components to be measured in planes defined by the upstream flow direction and the axis of the cylinder. The third (out-of-plane) velocity component was then calculated by integrating the continuity equation. Subsequently, the velocity field information has been manipulated and converted into time-averaged information.Discussion of the measured results confirms that colour displays are an invaluable aid to understanding this complex fluid flow situation since they reveal substantially more information than grey-scale plots of the same data. In particular, the source of the horseshoe vortex system can be identified when colour plots of the time-averaged velocity and vorticity distributions are obtained. A limited amount of information on the unsteady vortex structures appearing in the end-wall region upstream of the cylinder is also presented. Finally, the experimental findings are discussed in relation to the results of previous workers.     

绕圆柱体自由表面磁流体流动和传热的研究

本文对在不同雷诺数下,绕圆柱体的磁流体自由表面流动及传热进行了模拟,分析了磁场对绕流圆柱尾迹和涡分离的影响,获得了两种雷诺数下的电磁力密度、流场和温度场分布。结果表明,磁场不仅影响了流动的形态,而且对湍流有抑制作用,降低了自由表面的更新机制,从而减少了传热能力;在相同的Hartmann数下,相比低雷诺数下的流动换热情况,高雷诺数下的湍流不能被完全抑制,自由表面与尾迹的相互作用也较强,因而自由表面换热也较强。     

A numerical study of transient flow around a cylinder and aerodynamic sound radiation

The fully 3D turbulent incompressible flow around a cylinder and in its wake at a Reynolds number Re = = 9×10⁴ based on the cylinder diameter and Mach number M = 0.1 is calculated using Large Eddy Simulations (LES). Encouraging results are found in comparison to experimental data for the fluctuating lift and drag forces. The acoustic pressure in far-field is commutated through the surface integral formulation of the Ffowcs Williams and Hawkings (FWH) equation in acoustic analogy. Five different sound sources, the cylinder wall and four permeable surfaces in the flow fields, are employed. The spectra of the sound pressure are generally in quantitative agreement with the measured one though the acoustic sources are pseudo-sound regarding the incompressible flow simulation. The acoustic component at the Strouhal number related to vortex shedding has been predicted accurately. For the broad band sound, the permeable surfaces in the near wake region give qualitative enough accuracy level of predictions, while the cylinder wall surface shows a noticeable under-prediction. The sound radiation of the volumetric sources based on Lighthill tensors at vortex shedding is also studied. Its far-field directivity is of lateral quadrupoles with the weak radiations in the flow and cross-flow directions.     

Surface roughness effects on the turbulence statistics in a low Reynolds number channel flow

A high-resolution particle image velocimetry was used to characterize a low Reynolds number turbulent flow in a channel. Experiments were conducted over a sand grain-coated surface of large relative roughness, and the results were compared with measurements over a smooth surface. The roughness perturbation significantly modified the outer layer. Even though the streamwise Reynolds stress shows less sensitivity in the outer layer to the boundary condition, significant enhancements were observed in the wall-normal Reynolds stress and the Reynolds shear stress. These modifications were considered as footprints of the larger-scale eddies transporting intense wall-normal motions away from the rough wall. A quadrant decomposition shows that strong and more frequent ejections are responsible for the larger values of the mean Reynolds shear stress over the rough wall. The results also indicate that spanwise vortex cores with mean vorticity of the same sign as the mean shear are the dominant smaller-scale vortical structures over the smooth and rough walls. A linear stochastic estimation-based analysis shows that the average larger-scale structure associated with these vortices is a shear layer that strongly connects the outer layer flow to the near-wall flow. A proper orthogonal decomposition of the flow suggests that the large-scale eddy is more energetic for the rough wall, and contributes more significantly to the resolved turbulent kinetic energy and the Reynolds shear stress than the smooth wall.