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

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

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

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

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

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

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

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

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

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

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

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

固体燃料冲压发动机导弹推进性能的数字模拟

固体燃料冲压发动机导弹推进性能的数字模拟加拿大Ｃａｒｌｅｔｏｎ大学，为模拟固体燃料冲压发动机推动的导弹在零攻角下的超音速飞行研制了一种计算机程序，采用经验阻力数据计算射程、飞行速度及时间。首先计算通过锥形脱体激波的变化，在其下游，超音速气流进入进气道...     

微通道流动转捩数值模拟

应用CFD计算软件ANSYS CFX对水力直径为0.4mm的矩形截面微尺度通道中的气体流动进行了数值模拟.采用了3种流动模型,即γ-Re_θt转捩模型、层流模型和剪切应力输运（SST）模型,并将模拟结果与实验结果进行了对比.模拟结果表明:层流模型及SST模型不具有预测转捩的能力;γ-Re_θt转捩模型捕捉到了临界雷诺数且摩擦因数相对误差在20%以下.在未达到临界雷诺数之前,通道后部已经出现局部湍流区,且局部湍流区随着雷诺数的增大而增大.壁面摩擦因数的最低点对应的流向位置与通道中前后湍流区的交结位置几乎一致.      

拉瓦尔型微喷管性能的DSMC模拟

运用DSMC方法,对轴对称拉瓦尔型微喷管流动现象进行了仿真模拟,分析了不同壁面边界条件对喷管流动性能的影响。研究结果表明,不同的壁面反射模型对喷管内的流场结构及推力有较大影响;在同等壁面热条件下,随着动量调节系数的增加,喷管粘性损失增大,喷管的推力和比冲都将下降;对不同光滑度与清洁度的喷管表面,一概采用完全漫反射模型,将低估喷管推力性能。