Two-dimensional viscous vortex flow around a circular cylinder

A code based on a vortex method for simulating viscous vortex flow around a circular cylinder has been developed. The original method to obtain the no-through and no-slip condition introduced here is described in detail. The diffusive part of the vorticity transport equation is treated using a deterministic method based on the solution of the heat equation. Results have been obtained and are presented concerning the vorticity field, the velocity field, the evolution of drag and lift coefficients and the Strouhal number.     

绕流圆柱的自由表面磁流体流动实验

利用PIV(Particle image velocimetry)光学测量系统,在雷诺数分别为7 767和11 60时,实验研究了磁场对绕流圆柱的自由表面流动的影响.实验研究结果表明,磁场使边界层内的无量纲平均速度分布偏离对数分布,尾迹变窄变短,分离点向后移;当雷诺数增加到11 60时,在相同磁场强度作用下,磁场的影响减弱,从而使边界层内的无量纲平均速度分布偏离对数分布的程度减小.      

Blockage effect on the flow around a cylinder probe in calibration

Flow around a 2-D cylinder pressure probe placed in uniform flow,free jet flow,and wind tunnel flow was analyzed with potential flow theory and simulated with numerical method.Blockage effect was investigated under several typical flow Mach numbers.The result from numerical simulation shows a similar trend to the one from potential flow method while varies in quantity.Wind tunnel walls accelerate the flow near the probe and thus produce a blockage effect;Boundary of free jet flow,however,decelerates the flow and thus produces a "negative" blockage effect.A maximum incoming Mach number exists when the probe is calibrated in wind tunnel in high subsonic condition due to choking caused by shocks and shock induced separation.The critical Mach number varies with blockage ratio,which makes high Mach number impossible to achieve in large blockage ratio condition.The blockage effect itself is unavoidable for calibration or measurement although a sufficiently small blockage ratio brings minor effect.Correction can be implemented based on the numerical simulation result presented in this paper and further works.      

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

绕圆头回转体通气空化流型的实验研究

为了深入了解通气空化流动现象,利用高速全流场显示技术,对绕圆头回转体通气空化流型进行实验研究。结果表明,重力效应和通气量对通气空化的多相流流型起主要作用。定义了弗洛德数和通气率两个无量纲数,将绕圆头回转体通气空化分为5种多相流流型,即透明空泡、透明气弹、透明分层、水气混合以及半透明水气混合。流动参数对流型的影响分为2个阶段,即重力起主要作用阶段和重力效应不明显阶段。在重力起主要作用阶段,通气率一定时,随着弗洛德数的增大,附着弹体的空泡倾斜程度变小,弹体上表面的断裂空泡转变为贴着弹体壁面的稳定空泡;弗洛德数一定时,随着通气率的增大弹体上表面断裂空泡的尺度不断增大。在重力效应不明显阶段,通气率一定时,随着弗洛德数的增大,雷诺数变大,流场的湍流强度增大,空泡尾流区域水气交换的程度加剧;弗洛德数一定时,随着通气率的增大,通气空化数减小,绕弹体的云雾状空泡逐渐转变为透明空泡。最后,进一步分析了重力影响下透明空泡脱落的非定常过程,以及反向射流作用下云雾状空泡交替脱落的非定常过程。     

亚临界圆柱绕流的DES方法比较

针对地面结冰试验设施喷雾系统尾迹湍流模拟需求,考察了3种脱体涡模拟(DES)方法对三维亚临界圆柱绕流预测的准确性,比较和分析了瞬时流动特征和流场统计量。研究发现:1)从瞬时流动特征来看,k-ω雷诺平均NavierStokes(RANS)控制的区域对剪切应力输运k-ωDES(SSTk-ωDES)求解的准确性有显著影响;2)从流动参数统计量的误差范围来看,回流区长度和流向最小速度的预测质量能在一定程度上反应圆柱绕流数值模拟的准确性;3)综合比较,SSTk-ωDES预测值与试验值和大涡模拟(LES)数据最为吻合,具有应用于地面结冰试验设施喷雾系统尾迹湍流计算的潜力。     

Simulation of the flow around a rotating circular cylinder taking into account laminar-turbulent transition

We consider a problem of simulating a viscous incompressible fluid flow past a rotating cylinder in the region of critical Reynolds numbers as well as the simulation of separated bubbles in the boundary layer and the reverse Magnus effect. The Fluent 12 software is used for numerical simulation. The transition k-kl-ω model of turbulence makes it possible to consider the laminarturbulent flow regime transition.     

Simulation of the flow around a stationary circular cylinder taking into account laminar-turbulent transition

We consider a problem of a stationary incompressible viscous fluid flow around a flat circular cylinder. In the vicinity of the critical Reynolds number Re _cr a stepwise drop of the cylinder drag takes place, which is called the drag crisis.     

带射流冲击短扰流柱排内的流动和损失

用五孔探针测量了带射流冲击的短扰流柱排内的流场，并用压力扫描阀测量了端壁和柱面的压力分布，分析了在涡轮叶片尾缘区内射流冲击强化扰流柱排通道内换热的机理，研究发现，带射流冲击时流动在靠近扰流柱表面附近速度较大，在对称中心区域有大片低速区，压力系数远远小于无射流时情况，气流经过孔板后压力系数迅速下降，达到最小值，沿流动方向压力系数逐渐恢复。随喷射雷诺数增大，扰流柱表面的分离点位置后移，总结出了实验工况范围内带射流冲击的短扰流柱排内压力损失系数与雷诺数之间的经验公式，便于工程实际应用。     

低亚声速和跨声速矩形柱绕流的大涡模拟研究

采用基于非结构网格的有限体/有限元混合格式和大涡模拟的方法求解可压缩Navier-Stokes方程,研究了不同长宽比矩形柱低亚声速和跨声速绕流的流动特性。在雷诺数为22000时,对来流马赫数等于0.1和0.75,截面长宽比分别为1∶1、2∶1、3∶1和4∶1的矩形柱绕流进行了大涡模拟,以研究长宽比和压缩性对矩形柱绕流流场的影响。马赫数为0.1时,Strouhal数随着长宽比的增大先降低再增大然后再降低;长宽比为3∶1和4∶1时会有流动的再附产生;柱体上表面的三维特性在长宽比大时更明显。马赫数为0.75时,Strouhal数随着长宽比的增大逐渐减小;湍流脉动和涡脱落受到抑制;方柱的近尾迹区域,有两种形成机制不同的局部超声速区。     

Efficiency of a gas screen in the supersonic dispersed flow under conditions of inertial deposition of particles on a protected surface

The results of numerical studying the influence of inertially deposited particles of the condensed phase on the efficiency of a gas screen in a supersonic dispersed flow are presented. The possibility of realizing the Leont’ev paradox on the adiabatic surface region streamlined by a supersonic dispersed flow, involved in attaining the temperature of the surface being protected smaller than the coolant temperature on the permeable section of screen formation is specified.     

惯性颗粒在非均匀加速流中的输运问题

基于新型碳基颗粒增强或碳基纤维材料的再入体热烧蚀问题对飞行器的气动力/热性能影响机制尚不明确。初步研究不考虑化学反应和气体引射效应,把热烧蚀问题中的颗粒剥离问题模化为驻点流中的离散颗粒启动问题。由于高速气流在驻点温度达上万度量级,局部马赫数接近于0,研究近似假设靠近壁面的流体动力学为不可压缩的,并使用不可压缩颗粒解析的直接数值模拟方法（PR-DNS）研究一个惯性颗粒在平面挤压流和壁面驻点流中的动力学,从而揭示在驻点流中的颗粒输运机制。研究发现颗粒在壁面驻点流中的输运机制是水平输运,几乎没有垂直输运,与平行剪切颗粒床输运非常不同,这是驻点流特有的垂直向下挤压作用导致的。     

高速气流冲击柱状泡沫多孔体的传热特性分析

采用流体/多孔区域一体化单区域算法,数值研究了高速绕流条件下前置于圆柱体前缘表面的柱状泡沫多孔体内部的传热特性。基于蒙特卡罗法考虑多孔域内的辐射热效应,分析了变化多孔区域长度和多孔阻力特性对模型激波阻力和前缘多孔区域气动热的影响。结果表明:在圆柱体前缘安置一定长度及带有适当阻力特性的泡沫多孔材料,可同时减小整体激波阻力并降低前缘表面的气动热效应。在模拟工况下,无量纲长度1.0、黏性阻力系数0.2×107m-2及惯性阻力系数200m-1的前缘泡沫多孔可减小激波阻力13.5%,降低约75%的前缘表面的平均气动热流密度。保持无量纲长度不变,减小泡沫多孔区域惯性阻力系数会降低激波阻力,但会略微增加前缘壁面气动热流密度。      

烧蚀环境下的圆柱绕流计算模型

基于多学科理论建立了高温高压气流环境下圆柱烧蚀、剥蚀的数理模型。利用离散涡方法计算流场与圆柱表面压力分布;采用三方程烧蚀模型计算热化学控制机制下的烧蚀速率,并将其与扩散控制机制下的结果相比较以确定烧蚀量;引入颗粒轨道模型求解剥蚀颗粒的运动。基于该模型对高温高压燃烧室内圆柱形烧蚀试件的绕流场、烧蚀量以及剥蚀颗粒的运动轨迹进行了编程计算分析。研究表明:低雷诺数条件下烧蚀环境对绕流流场的影响不大,而较高雷诺数条件下烧蚀绕流流场与无烧蚀绕流流场的涡街分布存在显著差别;在低雷诺数条件下流场对烧蚀速率的影响甚微,较高雷诺数条件下烧蚀速率略大于无流场烧蚀情况;涡量的存在决定着剥蚀颗粒的运动轨迹与分布,进而影响尾流热场与烧蚀速率。该模型可为热防护材料的烧蚀实验提供参考,并为烧蚀、剥蚀过程的非线性分析提供数据支持。     

多尺度湍流模型在圆柱绕流中的应用

使用基于可变时间间隔平均方法的多尺度湍流模型,对高超临界区雷诺数下非定常圆柱绕流进行了数值模拟.结果表明:多尺度湍流模型在升阻力系数预测上与实验值一致,而标准k-ε模型和SST(shear stress transport)k-ω模型却只能兼顾其一;在旋涡脱落特性预测上与实验值相一致,其误差小于SST k-ω模型的5.8%和标准k-ε模型的37.6%;在表面平均压力系数分布预测上其误差仅为3.6%,明显优于SSTk-ω模型的13%和标准k-ε模型的53.7%;在表面摩擦因数分布及分离角度方面与实验结果相吻合,且分离角度误差仅为0.78%,结果优于标准k-ε模型的1.04%和SSTk-ω模型的1.83%,充分验证了该湍流模型应用于复杂湍流模拟预报的潜力.     

绕轴旋转圆柱结冰特性结冰风洞试验

绕轴旋转物体结冰是一种常见的结冰现象,如风力机叶片和直升机旋翼结冰等。针对该旋转模型的结冰问题,以绕轴旋转圆柱为对象,进行了结冰风洞试验研究。试验在自行设计的利用自然低温的结冰风洞系统中完成。在对结冰风洞试验能力进行验证后,对具有不同直径、转速和结冰时间的旋转圆柱结冰进行了试验,建立了旋转圆柱的结冰分析与评价方法,分析了转速、圆柱直径、结冰时间等对圆柱结冰形状主要特征量的影响规律,包括结冰面积、无因次结冰面积、驻点厚度、无因次驻点厚度、驻点偏转角、无因次结冰上下极限等,获得了旋转圆柱的结冰特性。在此基础上,通过回归正交试验设计,构建了旋转圆柱结冰的预测模型,并进行了验证与分析。研究结果可为绕轴旋转物体结冰特性研究提供参考和借鉴。     

轴向流中圆柱体的非线性动力学分析

通过引入流体摩擦力的二次阻尼项,推导出了轴向流中两端简支圆柱体的新的运动微分方程。对添加二次阻尼项的运动微分方程进行了无量纲化和Ritz-Galerkin离散。通过数值模拟,对两端铰支圆柱体的稳定性进行了系统的分析,并得出颤振发生的临界流速条件。     

A combined method for calculating the trajectories of suspended particles

A combined method for calculating the trajectories of suspended particles in the hydrodynamic velocity fields is proposed. The method utilizes capabilities of existing CFD codes and parallel computing on a computer with a graphics card that supports CUDA technology. The software package that implements the method is described.     

发动机缸内滚流运动的稳流实验与三维数值模拟

采用变压差法的涡流动量计式实验台和Ricardo评价方法,对四气门发动机缸内滚流运动以及气道流通系数进行了实验研究;并且根据气道稳流实验的状态与环境参数,采用RNG k-ε湍流模型、混合网格、有限体积法以及SIMPLE算法,建立了与之相对应的三维数值模拟分析模型,得到了气道内部与缸内流场更为详尽的数据信息.结果表明,实验气道可以在缸内产生滚流,但也存在局部的缺陷;数值模拟与稳流实验的结果吻合良好,模拟分析模型具有较高的准确性和实用价值.