二维有限长度柔性壁面上T-S波演化的数值研究

受自然界鸟类羽毛的柔性特征启发,利用数值模拟的方式开展了柔性壁面对亚声速边界层中T-S波演化的影响研究。刚性壁面上的数值结果与线性理论吻合得很好,验证了数值方法的可靠性。在此基础上,将部分刚性壁面替换为柔性壁面,结果表明,柔性壁面可以抑制T-S波在空间上的增长,从而推迟边界层流动转捩。壁面的变形不只跟随T-S波的波形,还因为柔性段与刚性段相接的前缘和后缘引起与扰动源频率相同的更大尺度的壁面波动,壁面的实际变形由这几种波叠加而成。开展的参数研究结果表明,增大表面的质量密度对于柔性壁面衰减扰动的效果几乎没有影响;增大表面张力和增加底部支撑的弹性系数可以增加壁面的刚性,减小壁面变形的幅度;增加阻尼可以抑制柔性段前后缘产生的大尺度壁面波动的传播,而对跟随T-S波的变形影响不大。总体上,柔性壁面的变形程度越大,其扰动的抑制效果越强。      

带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

基于欧拉方程的高超声速飞行器的壁面流线生成计算

提出了一种新的利用表面流函数法计算高超声速飞行器表面流线分布的方法.首先提出了表面流函数的概念,并通过理论推导,得到了表面流函数与表面流线的关系;然后运用结构化网格求解三维Eu-ler方程,计算得到高超声速钝头体的边界层外缘外部无粘流场气流参数;最后利用无粘流场气流参数和表面流函数的方法计算了高超声速飞行器的精确表面流线分布.结果表明,在无攻角和攻角小于20°的情况下均可以得到较好的壁面流线分布.高精度的表面流线的得到,为利用边界层内粘性主导区域的积分方程等方法进一步精确预测高超声速飞行器表面的气动加热奠定了基础.      

非对称来流下隔离段动态压力特性

为了探讨非对称来流下矩形隔离段内动态压力特性,用直联实验方式以及动态压力测量技术进行了试验,并用统计分析方法分析了数据。实验结果表明,入口来流的非对称性对上下壁面压力脉动大小以及传播平均速度有较大影响,而对频率无多大影响。激波串区域内壁面压力脉动向下游传播比向上游传播衰减得快。在隔离段出口超声速条件下,观测到的压力脉动频率主要在70 Hz以下,而在出口亚声速条件时,压力脉动的频率不仅有70Hz以下的部分,而且还有100~200 Hz之间的部分。     

基于相干结构的大涡模拟壁面模型

首先对湍流边界层近壁区相干结构群进行研究,根据共振三波理论使用沿流向周期性发展的多个三维波构造其初始结构,通过直接数值模拟得到该相干结构群的演化特性.在此基础上,提出了基于近壁区相干结构流动参数的新型大涡模拟壁面模型.确定了此壁面模型与主流大涡模拟衔接的原理方案.计算结果表明:该壁面模型在减少计算量的同时,能够为主流区大涡模拟提供重要的壁面大尺度结构信息,使计算取得准确的结果,实现了壁面模型的基本功能.      

湍流模型对叶片表面换热计算的影响

为了比较湍流模型对涡轮叶片外换热计算结果的影响,采用五种湍流模型及两种壁面函数计算了NASA-MarkII导向叶片表面换热并与实验数据进行了对比,结果表明:在相同边界条件下不同湍流模型的计算结果有很明显的差别,即使是同一种湍流模型,如果采用不同的壁面处理函数其计算结果也是有很大差别的;某些湍流模型的计算值只是在某个区域较为理想,还不能找到在整个叶片表面计算结果与实验数据较为接近的湍流模型.在尚无普适性较好的湍流模型的情况下,研究在不同计算域采用不同湍流模型和壁面处理函数的计算技术,是一种较好的可行的方案.      

局部/组合抽吸对压气机叶栅分离流动控制的机理研究

为了研究不同抽吸形式对压气机内部分离流动的控制效果和机理,基于数值模拟方法,对内部同时存在角区分离和附面层分离的直列叶栅进行抽吸计算,共设置8套不同的方案进行了详细的探究。结果表明:对于控制分离的效果和机理,不同抽吸形式之间存在着一定差异;在分离发展的不同阶段,实施抽吸对其控制效果也有所不同。近端壁吸力面抽吸方案SS3同时减小了角区分离和近端壁处的尾缘附面层分离,尾部近吸力面沿展向静压分布为正"C"型,使角区内的低能流体向叶中流动,角区流动得到改善,但导致叶展中部流场的恶化;端壁抽吸和全叶高组合抽吸则使尾部近吸力面静压沿展向分布更均匀,低能流体的迁移现象减弱。基于以上结果可得:组合抽吸可以更有效地控制叶栅内部流动分离,降低损失,提高叶栅流通能力。     

基于剪切敏感液晶涂层的壁面摩擦力矢量场全局测量方法

针对壁面摩擦力测量问题,建立了一种基于剪切敏感液晶(SSLC)涂层技术的平板表面摩擦力矢量场全局测量方法。该方法利用SSLC涂层在摩擦力作用下的颜色变化特性(不同方向显示不同颜色)并结合其颜色变化与摩擦力大小之间的校准关系解算摩擦力矢量的方向和大小,能够测量整个待测区域的摩擦力矢量分布。应用所述方法测量了平板表面凸起物绕流的摩擦力矢量场。试验结果表明,该方法不仅能够高分辨率测量出平板表面凸起物绕流的摩擦力矢量场,而且能够研究凸起物的尾迹区随着流动速度增加的发展过程。     

离子推力器放电室阳极壁面电流密度分布特性研究

为了准确掌握离子推力器放电室阳极壁面电流密度分布特性，并深入理解阳极壁面处等离子体运动特性，设计了近阳极壁面等离子体诊断的具体实施方案，并基于LIPS-200离子推力器开展了近阳极壁面处等离子体诊断试验研究，得到了主要磁极附近壁面等离子体参数，并得到阳极壁面吸收电流密度分布特性。试验结果表明：LIPS-200离子推力器阳极壁面处主要磁极附近的等离子体密度范围为，测试点的电子温度范围为，壁面电流密度范围为；柱段壁面电子温度相对锥段较低，但电流密度较大，尤其在中间极靴位置电流密度最大，约为阴极极靴处电流密度的3倍，约为屏栅极靴处电流密度的2倍，阳极电流主要在放电室中间极靴处发生损失。     

扩散界面浸入边界法结合壁面模型在湍流模拟中的应用

提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。