平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

洁净室流场大涡模拟

应用湍流大涡模拟方法对矢流洁净室和全顶棚送风下部两侧回风式洁净室进行了数值模拟。计算中使用了隐含流线抑风耗散作用的具有三阶部分展开的Taylor-Galerkin离散格式，导出了经Gauss滤波后的基本方程组的有限元列式，应用涡粘性亚格子模式的湍流大涡模拟方法进行数值模拟。数值计算结果与实验结果符合较好，表明湍流大涡模拟能够较准确地捕捉到洁净室流场中旋涡的位置和尺度，在洁净室的数值仿真方面，与标准的K－ε湍流模式相比，湍流大涡模拟结果具有较高的可靠性。     

三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

大涡模拟中亚格子模型的改进及其在槽道湍流中的应用

在用大涡模拟的方法计算具有强剪切的槽道湍流时, 常用的亚格子应力模型, 包括考虑壁面修正的模型, 都会出现平均流剖面偏高的现象.这表明在壁面附近亚格子应力模型还不能描述实际情况.针对这一问题, 修正了亚格子雷诺应力模型的壁面函数, 得到了较好的计算结果.用修正后的模型计算出的平均速度分布、均方根速度的分布以及雷诺应力的分布, 均与直接数值模拟(DNS)的结果吻合较好.      

高对流Mach数三维混合层转捩特性分析及小激波结构模拟

采用空间大涡模拟方法对超/超混合、超/亚混合两类三维可压缩平面混合层转捩及其全湍流流场进行了研究,认为混合层动量厚度饱和点可作为流场转捩完成的标志.计算所得到的线性扰动波激励下流场转捩拟序结构与随机扰动下自然失稳结构以及文献结果进行了比对,其结果是一致的,表明了引入线性扰动激励来研究流场转捩结构是合理的.同时,本文还在较高对流Mach数流动下得到了三维流场动态小激波结构,其分布具有非对称特性,且形状与实验及直接数值模拟结果相似.不同条件混合层转捩计算表明:高对流Mach数下混合层转捩以Λ涡结构的形成和发展为主导机制,受扰动及对流条件的影响Λ涡结构不尽相同,某些情况下流场出现二维与三维涡结构共存现象.充分发展湍流区域,流场脉动速度分量量级相同,湍流压缩效应随着对流Mach数提高而明显增强.     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

基于格子Boltzmann方法的方通道湍流的大涡模拟

基于格子Boltzmann方程的大涡模拟方法,对以摩擦速度、方通道水利直径为特征尺度,雷诺数为300的直方通道内湍流流动进行数值计算.利用多松弛时间格子Boltzmann方法来模拟流场的流动,切应力改善亚格子应力模型来模化滤波后的非封闭项.将模化后的亚格子应力与格子Boltzmann方法中的松弛时间相关联,使得松弛时间当地化,从而能够准确地模拟湍流.对湍流的平均流向速度、平均二次流速度以表征湍流强度的均方根速度以及不同截面流向瞬时涡做了计算和评估.计算结果与直接数值模拟、实验数据相吻合,证明了格子Boltzmann方法在计算通道湍流中的精度.      

DBD涡流发生器及其在角区流动控制中的数值研究

为明确介质阻挡放电(DBD)涡流发生器对马蹄涡的影响,在采用唯象模型的基础上,通过数值方法研究了DBD涡流发生器诱导产生流向涡的结构特性及其对马蹄涡的控制特性。研究结果表明,流动在激励器上游边缘处形成羊角涡,自由剪切层卷入流向涡的涡核并为其提供持续涡量;在柱体根部角区流动中,当对称面两侧激励器诱导流动指向对称面,诱导涡与马蹄涡环绕方向相反时,马蹄涡可以得到有效抑制,反之,则控制效果不佳。最后得出,诱导涡对下游马蹄涡的控制机制体现在其黏性扩散作用、掺混作用以及低压效应3个方面。     

轴流风机的进口径向畸变研究

对轴流风机的进口速度畸变情况进行了分析,利用改进的双激盘理论通过进口速度分布的改变来分析进口畸变对于性能的影响,并且将模型的计算结果与实际的实验结果相比较,得到一种合理的速度分布.      

Velocity Distribution of Interstellar H Atoms in the Heliospheric Interface

In this paper we present first results of a numerical computation of the velocity distribution function of interstellar H atoms in the heliospheric interface, the region of the solar and interstellar wind interaction. The velocity distribution is a key tool to evaluate uncertainties introduced by various simplified models of the interface. We numerically solve the kinetic equation for gas of H-atoms self-consistently with the hydrodynamic equations for plasma. Neutral and plasma components are efficiently coupled by charge exchange. The interaction disturbs the atom velocity distribution, which is assumed to be Maxwellian in the circumsolar local interstellar medium. It is shown that besides ‘original’ interstellar atoms, there are three other important atom populations originating in the heliospheric interface. Velocity distribution functions of these populations at the heliopause are presented and discussed. This revised version was published online in August 2006 with corrections to the Cover Date.     

非均匀速度来流条件下液体射流喷雾特性

实验研究了非均匀速度来流条件下液体横向射流喷雾特性.实验利用多孔板实现速度的不均匀,应用片光照相得到喷雾中心截面照片,使用相位多普勒粒子分析仪(PDPA)测量液滴的直径和速度.实验结果表明:在非均匀速度来流中,射流的弯曲和破碎机理有所不同,导致穿透深度和雾化特性发生改变.与均匀速度来流相比,正速度梯度来流时射流穿透深度...     

随机恒速运动目标的搜索方程及持续探测概率

根据目标位置的初始分布和速度分布,给出了随机恒速运动目标在任意时刻位置和速度分布密度函数。定义了探测函数、目标存在和探测不成功的联合概率密度函数和目标存活概率函数,建立了基于目标速度分布的搜索方程并给出了其特征迹线解。在目标初始分布和速度分布均为圆正态分布条件下,分析了搜索方程特征迹线的基本特征和持续探测发现概率积分域的形态。给出了一个以直升机吊放声纳探测潜艇为背景的持续探测发现概率的算例,算例表明,对随机恒速运动目标持续探测的发现概率,与对静止目标进行持续探测的发现概率相似,服从时间的指数规律。     

低比转速泵叶轮流道内部流动的PIV试验

介绍了粒子图像测速仪(PIV)、泵试验装置和试验方案.采用PIV技术,对3种流量下叶轮流道内部3个不同测量平面上的流动进行了试验研究,获得了速度矢量分布.探讨了速度分布沿叶轮周向和轴向的变化规律以及流量对速度分布的影响.结果表明,随着流量的减小,吸力面速度略有减小,压力面速度变化较为明显,当流量达到额定流量的0.6倍时,在压力面附近存在回流现象;在叶轮出口处,流量对相对速度方向没有明显影响;叶轮流道内部速度沿轴向分布均匀,绝对速度沿半径方向呈增加趋势,无量纲绝对速度的周向分量v_u/v从压力面到吸力面呈下降趋势,且随流量的减小而增大.研究结果为叶轮设计以及离心泵内部流动的试验研究提供了借鉴.      

Numerical investigation of windage heating within shrouded rotor-stator cavity system with central inflow

The rotating disk surface temperature rise due to windage heating effect by numerically modeling the turbulent flow within a rotor-stator cavity which is available with a peripheral shroud and imposed through airflow was dealt with. The windage heating may be defined as viscous friction heating caused by relative velocity differences across the boundary layers between the fluid and the rotating disk surface. The kinetic energy dissipation process could transform the rotating shaft power into thermal heating. Commercial finite volume based solver, ANSYS/CFX was employed to numerically simulate this physical process by using the shear stress transport (SST) turbulence model. CFD results include the rotating disk surface temperature axial distribution and tangential velocity distribution of the fluid domain. The velocity difference between the result obtained by particle image velocimetry (PIV) experiments and CFD simulation are within 5%. The adiabatic disk temperature rise can be calculated by the tangential velocity of disk and fluid in large gap ratio and turbulent parameter. CFD temperature distribution results and those estimated via velocity differences are within 10%.     

单通道民用飞机客舱气流组织数值仿真研究

建立了单通道民用飞机客舱气流组织数值仿真模型,对不同的送风角度、送风速度和送风比例等状态的客舱气流组织进行了数值仿真研究,通过对比分析不同状态的仿真结果的流场、温度场、预期不满意百分率PPD、预测平均反应PMV、吹风感所引起的不满意率DR及其垂直温差导致的不满意率PD等热舒适性参数,发现水平角度送风、上风口较低风速和下风口较大风速、上风口风量少于下风口的设计较为合理。     

流通面积比对旋流杯油雾速度场的影响

设计了3种不同一二级流通面积比的旋流杯试验件，在气压自模区工况下，运用粒子图像测速仪（PIV）对其出口油雾速度场进行测量，研究了不同流通面积比对旋流杯出口油雾速度场的影响，分析了油雾速度场的结构及其速度分布的变化规律。结果表明:随着一二级流通面积比的增加，旋流杯出口油雾速度场呈现出回流区面积和回流区长度逐渐变小、回流速度逐渐变大、出口射流张角逐渐减小的趋势。此外，从径向速度沿径向分布可以看出，随着流通面积比的增加油雾速度场的抗偏斜能力具有增强的趋势。      

雷诺数对涡轮部件性能的影响

利用数值模拟手段研究了雷诺数对低压涡轮部件性能的影响,分析了低压涡轮内部流动的变化。结果表明低压涡轮效率随雷诺数降低而下降,且其变化趋势为非线性变化;叶片表面速度分布亦发生较明显的变化,导致叶型损失与二次流损失增大;结果还表明叶片表面速度分布的设计对涡轮性能随雷诺数变化的规律具有重大影响,在低雷诺数条件下涡轮部件的设计中需加以注意。      

跨流域高超声速绕流Boltzmann模型方程并行算法

通过对Boltzmann方程碰撞积分进行模型化处理,提出了统一描述各流域复杂高超声速流动输运现象的气体分子速度分布函数控制方程,使用离散速度坐标法对分布函数方程所依赖的速度空间离散降维,构造出直接求解分子速度分布函数的气体动理论耦合迭代数值格式,研制了复杂飞行器高超声速绕流气动热力学计算模型。基于对气体动理论数值计算方法内在并行性、变量依赖关系、数据通信与并行可扩展性的分析研究,使用区域分解并行化方法提出了新型的气体动理论数值算法并行方案;研究了数据的并行分布与并行执行特征,开展了大规模的并行化程序设计,构造了可稳定运行于成千上万CPU的高性能并行算法,用以模拟各流域复杂飞行器的高超声速绕流问题。以稀薄流到连续流环境下不同Knudsen数、不同马赫数的可重复使用类球锥卫星体及翼身组合复杂飞行器等气动力、热绕流问题为研究对象展开大规模并行计算,并进行算法验证,所得计算结果与理论分析、直接模拟蒙特卡罗方法(DSMC)的模拟值及有关实验数据吻合较好,揭示了飞行器跨流域高超声速下的复杂流动机理与变化规律,提供了一条能够可靠模拟高超声速飞行器跨流域气动力及热问题的统一的算法应用研究途径。