后向台阶层流边界层分离实验研究

本文应用偏振差动式激光测速仪对后向台阶层流边界层分离、再附及发展为湍流边界层的整个过程进行了精细的量测,得到了时均速度、湍流强度等的分布特性。我们发现,在本实验条件下,层流边界层分离后在再附点处为过渡流,进一步发展最后形成湍流边界层。在湍流边界层区域,流速分布遵循对数律,但积分常数比平板湍流边界层中所得值要大。     

来流湍流度对边界层发展的影响

应用ＬＤＶ（Ｌａｓｅｒ Ｄｏｐｐｌｅｒ Ｖｅｌｏｃｉｍｅｔｅｒ）测速技术及氢气泡流动显示技术对水流中平板边界层的发展进行了详细的观测，研究了来流湍流度对边界层转捩、湍流边界层特性等的影响，得到了层流边界层、转捩区、湍流边界层区的分区界限，在本实验条件下，来流湍流度对边界层转捩有非常明显的作用，但对湍流边界层的特性影响不大。     

后向台阶层流分离剪切层特性研究

应用激光测速仪对生向台阶层流边界层分离产生的剪切层在再附前的发展进行了测量，得到了时均速度、湍流强度及剪切层厚度等的分布特性，并研究了来流湍流度对台阶后流场湍强度分布及剪切层发展的影响，实验中发现台阶后有一个流速降低的区域，对此从涡动力学的角度作了解释。     

磁云边界层内朗缪尔波活动现象的初步研究

对35例磁云边界层内朗缪尔波活动现象进行了初步分析, 发现两类磁云边界层内特有的朗缪尔波活动现象. 一类是相对于邻近鞘区和磁云本体, 整个边界层内朗缪尔波活动增强; 另一类是短时间的朗缪尔波爆发现象, 同时伴随着宽频带的多普勒频移离子声波活动. 考察了其中一例朗缪尔波爆发事件对应的高分辨率电子分布函数数据, 发现速度约为7×103 km/s的电子束流形成尾峰分布不稳定性导致了朗缪尔波的爆发.      

基于Transition SST模型的高雷诺数圆柱绕流数值研究

为了研究高雷诺数下圆柱绕流边界层的转捩现象和圆柱尾迹近壁区的流动特征,首先通过在典型雷诺数下采用Transition SST四方程转捩模型模拟圆柱绕流得到的结果与实验结果及采用SST k-ω两方程湍流模型模拟结果的对比分析,验证了Transition SST模型在模拟高雷诺数下圆柱绕流的优越性,并较为准确地预测出了圆柱绕流边界层的转捩现象及尾迹近壁区的流动特征。然后分别对亚临界区、临界区、超临界区和过临界区的圆柱绕流问题进行了数值模拟,分析了不同雷诺数下圆柱绕流的流场结构及圆柱表面压力系数、摩擦力系数的变化规律,研究了圆柱绕流近壁区的流动特征、边界层转捩的流动机理、转捩位置及其随雷诺数的变化规律。结果表明,亚临界区二维圆柱绕流边界层发生层流分离,无分离泡和转捩现象;临界区和超临界区二维圆柱绕流边界层先产生了分离泡现象,之后流动发生了转捩并在转捩后发生湍流分离;过临界区二维圆柱绕流边界层流动在转捩之后发生湍流分离,无分离泡现象;在临界区、超临界区和过临界区,二维圆柱绕流边界层转捩位置随雷诺数增大向前驻点移动。     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

超临界层流翼型优化设计策略

针对超临界层流翼型设计问题,提出一种两轮优化策略。采用γ-Reθt转捩模型耦合剪切应力输运(SST)模式的湍流模型对翼型边界层转捩进行预测。翼型几何参数化建模采用形状分类函数转换(CST)方法,设计变量为描述翼型几何特征的参数。第1轮优化的目的是尽量提高层流区域的比例,气动分析模型为基于Kriging模型的代理模型,优化算法为遗传算法,通过优化获得满足约束要求的层流翼型。第2轮优化目的是对第1轮优化获得的翼型进行微调,进一步提高翼型的升阻比,气动分析直接采用CFD程序,优化算法采用基于梯度的优化算法。算例表明,应用本文提出的两轮优化策略,可将超临界翼型NASA SC(2)0412优化设计成超临界层流翼型,翼型的上下表面层流区比例分别达到了55.5%和47.0%,升阻比提高了38.1%。     

高超声速尖锥边界层转捩数值模拟

边界层转捩对高超声速飞行器气动力和热产生重要影响.通过联立求解Favré平均层流脉动能方程与Favré平均Navier-Stokes方程,开展了高马赫数条件下尖锥边界层转捩位置的数值预测研究.将数值计算结果与可用的实验结果进行对比,结果表明增大单位Reynolds数可使转捩提前发生,攻角效应可使迎风面转捩延迟,背风面转捩提前.     

防冰表面的对流换热计算分析

利用动量和热边界层积分分析的方法对采用电热防冰的二维机翼对流换热系数进行了数值计算,比较了等温突捩以及变温光滑两种不同的边界层模型.层流向湍流转捩的起始点位置和过渡区域的长度由经验公式判断,并通过内外热流耦合迭代求解得到了防冰表面的平衡温度,与文献试验值进行了比较.所得结果表明:加入过渡区模型可较好地预测对流换热系数及防冰表面温度.同时,结果还表明环境压力变化对换热系数的的影响较大.      

粗糙度对涡轮叶片流动转捩及传热特性的影响

为研究表面粗糙度对涡轮叶片流动转捩以及传热特性的影响,在自行开发的CFD程序平台上提出了对γ-Re_θ转捩模型的粗糙度修正方法,并参考平板绕流和涡轮叶栅的实验数据对该方法进行验证。考虑粗糙度效应的影响,对Mark Ⅱ涡轮导叶5411工况进行数值模拟,得到如下结论:表面粗糙度对层流边界层换热系数影响不大,而对湍流边界层则有较大影响,进而显著改变壁面温度分布;与光滑壁面相比,5μm的等效沙粒粗糙度使吸力面湍流区域壁面温度升高约5.7K,100μm粗糙度使壁面温度升高28.4 K,增幅达5%左右;当壁面粗糙度较低时,激波干涉对吸力面边界层的转捩起主导作用,而当粗糙度大于某临界值时,其作用会使转捩位置突然变化,本算例中该临界值近似为150μm。      

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

行星际磁云边界层中的电子分布函数特征

基于WIND 飞船观测的1995---2006 年间的磁云事件, 研究了磁云边界层中电子的流动图样, 以及电子速度分布函数的特点与电子加热和加速的关系, 得出以下结果. ①磁云边界层中存在的电子流动图样, 包括各向同性、双向流动和单向流动等. ② 相比于背景太阳风和磁云本体, 通常情况下磁云边界层中电子分布函数的核心热电子成分 (E<60eV) 增多, 超热电子成分 (E>60eV) 在沿磁场垂直方向上增多, 而在沿磁场平行或反平行方向以单方向增加为主, 此外, 还在近1/10 的磁云边界层中观测到了高能电子的明显增多. ③对比研究了磁云边界层与磁云驱动的激波对电子速度分布函数的调制作用, 经过激波, 电子分布函数的超热电子成分在各方向上都有增加, 不同于磁云边界层中在沿磁场平行或反平行方向上超热电子成分以单方向增加为主, 表明二者有不同的形成机制. ④考察了磁云边界层中的波活动增强和电子分布函数及离子流量增加的对应关系. 上述观测和对比分析进一步表明了磁云边界层是一种重要的动力学结构, 磁重联是一种可能的形成机制.      

角区三维漩涡流动的新分离结构

分别用数值计算和PIV(Particle Image Velocimetry)实验方法研究角区层流边界层三维定常分离的流动结构,证实了角区确实存在有别于传统分离现象的附着鞍点结构,对称面上游的流线并非从壁面向上抬起从壁面分离,而是经由空间某个奇点向壁面附着.角区马蹄涡的传统分离鞍点结构和附着鞍点结构之间存在着一定的演化规律,影响参数包括模型头部钝度、边界层厚度和雷诺数.一定的边界层和雷诺数条件下随着模型头部钝度减小,角区马蹄涡将从传统的分离鞍点结构逐步过渡为附着鞍点结构;对一定的模型若雷诺数越大、边界层越厚则角区流动越趋向于传统的分离鞍点结构,反之则倾向于附着鞍点结构.     

北京航空学院学报1985年总目录(总第32期—总第35期)

280181282 北京航空学院学报1985年第飞期(总第32期)镍基纤料纤焊的不锈钢接头重熔温度和组织的研究庄鸿寿守恒型有限元法及其在跨音速流中的应用李椿聋含有l型中心裂纹的矩形薄板在突然位移作用下的瞬时效应 张孝思李维伯一种结构优化的迭代线性规划方法张承煦增进微分方程数值稳定性的途径热力参数计算法评定航机陆用机组L衬一1500现踢性能验收试验刘思永数合金恒幅载荷下腐蚀疲劳裂纹扩展影响因素和机理的综述郭淑铭建立合理运输结构的Fuzz夕模式李志东陈一青对流换热边界层内嫡产率分布研究沃尔什码分多址数字信息传输系统的性能分析沈…     

磁云边界层与间断的观测研究

利用WIND飞船的磁场和等离子体观测数据,分析了1995年2月至2003年8月之间82个磁云中的边界层事件.我们认为(1)磁云边界层中方向间断DD(Directional Discontinuity)类型中RD(Rotational Discontinuity),TD(Tangential Discontinuity),ED(Either Discontinuity),ND(NeitherDiscontinuity)的分布为37:18:44:1(%),与背景太阳风中的分布51:12:35:2(%)不同;主要区别在于RD与TD的比例变化. (2)磁云边界层的内外边界切向间断的比例很高,占总数的20%,而且两侧的密度和温度平均相对变化分别为|△N/N|=0.24和|△T/T|=0.19,大于边界层内部的平均值,显示了更多的切向间断特征;此外,磁云边界层中方向间断的出现频次约为太阳风中的1.87倍. (3)磁云边界层中方向间断的法向分布在θ-φ平面中不是随机分布,而是以θ=-24.90°,φ=217.49°为几何中心,主要是指向远离太阳的方向,而不是简单的各向同性分布.初步结果表明,间断是磁云边界层中的重要结构,它有着不同于背景太阳风间断类型比,为诊断磁云边界层的形成机理提供了依据.      

“跨音速非等熵全位势方法及粘流/无粘流相互作用的计算”通过鉴定

流体所朱自强教授为首的课题组自1986年历经5年研究的“跨音速非等熵全位势方法及粘流/无粘流相互作用的计算”课题于1991年底通过校级鉴定,1992年3月初通过部级鉴定。课题组发展了一套以非等熵全位势方法解和三维可压缩湍流边界层积分方程反方法解相耦合     

带喷流的超声速光学头罩流动显示

在马赫数3.8的超声速风洞中,以高时空分辨率的基于纳米示踪的平面激光散射(NPLS,Nano-tracer based Planar Laser Scattering)技术为实验手段,研究了有无喷流的超声速光学头罩流场的精细结构,清晰地再现了流场中的激波、膨胀波、剪切层和湍流边界层等复杂结构.通过分析时间相关的流场NPLS图像,可以发现流场结构随时间的演化特性.结果表明:无喷流情况下光学窗口上方的大部分流场处于层流状态;有喷流情况下剪切层的层流区域较短,在很短的距离内转捩至湍流状态;喷流出口压力高于外界压力情况下剪切层的转捩位置比压力匹配情况下较为靠前,光学窗口上方的涡结构也较为复杂.比较而言,后者对气动光学性能的影响更大.