过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

采用欧拉两相流法对翼型表面霜冰的数值模拟

基于欧拉两相流理论提出一种翼型表面霜冰的数值模拟方法.采用同位网格上的SIMPLE (Semi-Implicit Method for Pressure Linked Equations) 算法求解空气控制方程;提出一种可穿透型壁面边界模拟水滴对翼型表面的撞击,通过求解过冷水滴的控制方程得到翼型表面的水滴收集特性;利用冰层时间推进法模拟积冰过程,使用积冰法向生长假设生成积冰外形.对NACA 0012翼型在0°和4°迎角下的积冰情况进行了研究,冰形预测结果与实验结果一致.通过对结冰翼面的压力分布进行分析,表明积冰对翼型气动性能会产生不利影响.     

高超声速流动壁面催化复合气动加热特性

针对高超声速流动壁面催化特性,计算了不同壁面催化复合系数条件下的球锥驻点热环境。引入了经验证的数值求解Navier Stokes方程的方法,在不同壁温500K～2500K的条件下分别分析了O 2和N 2气体在壁面处的催化复合气动加热特性,得到如下结论：（1） 原子复合放热将提高近壁面温度梯度,改变近壁面组分分布;原子复合放热一部分加热飞行器形成组分扩散热流,一部分加热近壁气体提高近壁温度梯度。（2） 在壁面催化复合系数较小时,原子复合放热主要转化为组分扩散加热,对于不同壁面温度,壁面催化复合系数α<0.1时,单一气体反应组分扩散热流小于总热流的20%。     

重叠网格洞面优化技术的改进与应用

针对传统割补法在壁面狭缝处网格重叠失败的问题,提出了以单元顶点状态进行洞面优化的vertex方法及以单元大小进行洞面优化的volume方法的两种改进方法。这两种方法将传统割补法中切割和填补两步搜索方法改为一步搜索方法进行重叠网格的洞面优化。通过数值试验证明:两种洞面优化方法都能很好地解决壁面狭缝网格重叠问题,网格重叠区域流场数值传递正确,计算结果与试验值吻合,优化效率高,收敛速度快。volume优化方法与vertex优化方法相比可以减少由网格大小不对等引起的数值误差,计算结果和计算精度均优于vertex优化方法。     

霍尔推力器预电离对低频振荡及壁面腐蚀影响的研究

ATON型霍尔推力器提出后,基于缓冲区的预电离现象的系列研究随之展开。但是迄今为止,预电离率的选择仍然未有定论。文章基于SPT-100型霍尔推力器的真实工作尺寸,建立了二维轴对称全粒子网格质点法（PIC）模型对推力器放电室内的放电过程进行模拟,并通过施加不同的缓冲区预电离率着重研究了预电离率与放电通道低频振荡和推力器壁面腐蚀之间的关系。研究结果显示,缓冲区预电离率的提高可以有效抑制放电通道内的低频振荡,使推力器工作性能提高。但与此同时,由于撞击通道壁面的离子能量增加导致通道壁面腐蚀加剧。     

几种流体壁面剪切力测压管的性能研究

对四种壁面型流体剪切力测量差压计Ｓｔａｎｔｏｎ管，Ｋ型管，穴型管，及微孔表面差压管在湍流边界层条件下进行试验并对其结构特点进行比较表明，这四种差压计具有相同的特征，所测量的压差均为流量动量变化所致，其流体剪切力可以表示为τ＝Ｃ△Ｐ＾ｎ，为避免对标定数据过于繁复的多段拟合，并使差压计能在较大测量范围内保持拟合公式的精度，各类差压管在与固体壁面铅垂方向上的尺度值应尽量低，微孔表面差压管结构复杂，灵敏度     

部分覆盖构型对凹腔冷态和燃烧流场的影响

为了研究凹腔部分覆盖构型对燃烧室流场的影响,采用SST k-ε湍流模型和有限化学速率/涡破碎模型,对部分覆盖型凹腔的冷态和燃烧的二维流场进行了数值仿真。比较了贫燃和富燃条件下,覆盖板对凹腔燃料分布、流场温度和燃烧效率的影响。研究发现贫燃下,覆盖板有利于凹腔内燃料的积累,提高了凹腔的火焰稳定性,同时由于点火适宜的当量比区域位于凹腔内低速区,增强了凹腔点火能力;富燃下,覆盖板使凹腔内富燃环境进一步恶化,降低了凹腔的火焰稳定性和燃烧效率,主要燃烧区部分向下游转移。     

壁面温度对斜爆震发动机进气道流动特性影响的数值研究

为研究壁面温度条件对层流、转捩、湍流状态下斜爆震发动机进气道流场结构、流场参数的影响,选取Ma10级、具有曲面压缩段的斜爆震发动机进气道为研究对象进行数值模拟,对进气道壁面附近激波诱导分离区、热边界层的变化进行了深入探讨。数值模拟结果表明,进气道肩部圆弧过渡段出现的再层流化现象,壁面冷却对其起抑制作用,绝热壁面条件下再层流化程度最为严重。壁面温度的增加有利于延缓流动转捩,同时也导致了分离区尺寸的增加以及转捩、湍流状态下分离区主体位置逐渐前移,进气道内通道的转捩为分离诱导转捩,转捩位置主要受到分离点位置的影响,整体表现为壁面温度增加转捩位置前移。进气道出口顶板侧热边界层厚度随着壁面温度的增加逐渐变厚,转捩状态下热边界层厚度变化可达5%,温度峰值也随着壁面温度的增加逐渐增加,且峰值位置逐渐靠近壁面。壁面温度条件相同时,层流状态下热流、热边界层厚度均较小。转捩、湍流状态下进气道出口顶板侧热边界层较厚,约为层流状态3倍,同时转捩、湍流状态下热边界层厚度相差可达2%。     

基于L-M算法的液体火箭推力室壁面热流密度测量方法

为获取液体火箭推力室热防护方案构建的基础数据,结合单点瞬态法,提出了基于L-M算法测量壁面热流密度的方法。通过模拟热流密度分析算法中关键参数对求解精度的影响,给出关键参数的选取建议。在此基础上,对过氧化氢煤油推力室典型试验工况进行热流密度测量。结果表明：圆筒段末端壁面热流密度为3.35MW/m²,喉部位置单组元工况下壁面热流密度为2.21MW/m²,双组元工况下壁面热流密度为10.59MW/m²,热流密度变化过程与室压变化过程完全对应。变工况过程中近壁处燃气温度存在恢复过程,燃气温度稳定后,壁面热流密度也达到稳定状态。本文提出的方法可迅速、准确地测量推力室壁面热流密度。