高超音速化学非平衡钝劈绕流数值计算

 自1970年Davis提出粘性激波层方法以来,用数值方法求解高超音速轴对称钝体绕流问题国内外已做了大量工作,但至今仍未见到关于平面问题的计算结果。对于平面问题,虽然方程在形式上比轴对称简单,但由于二维效应,激波层较厚,用文献[2]的方法向下游区推进有困难。另外,在驻点线上采用极限关系式虽能克服方程的奇性,但驻点解对流向步长Δξ有依赖。     

轴流压缩机粘性尾迹噪声研究

本文研究了微弯机翼在正弦阵风作用下的脉动升力,并导出计算公式,分析了动、静叶栅几何及气动参数对转子脉动升力及噪声级的影响。     

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

跨声速开槽壁超临界翼型洞壁干扰的N—S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一     

虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

再入弹抛壳气动特性研究

本文用有限差分方法求解再入母弹头与其所抛壳体之间的无粘干扰流场,获得了体—壳的干扰气动力,并在准定常假设的基础上近似模拟了弹体与壳体的分离过程。数值研究的结果表明,体—壳间存在的干扰会导致壳的压心位置急剧变化,危及母弹安全;抛壳速度是体—壳安全分离的关键参数。     

基于FW-H方程的旋翼气动声学计算研究

由流体力学N S方程导出的非齐次波动方程———FfowcsWilliams Hawkings方程(简称FW H方程),可以精确地描述在静止流体中运动的物体与流体相互作用的发声问题。以FW H方程为理论模型,将旋翼桨叶运动发声问题等效为包含桨叶的任意运动控制面(声源面)的声辐射问题,并在旋翼绕流Euler方程数值模拟的基础上,在时域内计算了悬停旋翼和前飞旋翼的声场。应用于UH 1H和AH 1/OLS两种旋翼模型的气动声学计算表明:计算结果与噪声实验值符合良好;所研制的程序不仅能够较准确地计算单极子噪声和偶极子噪声,而且具有较强的跨音速四极子噪声预测能力。     

跨声速粘性流绕振荡翼型的非定常计算

本文采用时间和穴是二阶精度的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模式及ｑ－ω二方程微分模式，结合网格自适应技术，数值模拟Ｎ－Ｓ方程，计算了跨声速下的翼型非定常运动，包括俯仰，浮沉和前后平移振荡。结果表明，压力分布和气动系数与实验基本符合，微分模式和自适应网格能够显著提高激波和边界层计算精度。     

超声速绕凹角流动的压力分裂形式简化N—S方程数值解

本文从压力分裂形式的简化Ｎ－Ｓ方程出发，用Ｒｕｂｉｎ格式对超声速粘性流动进行了数值求解，给出了超声速绕凹角附着流及带有小分离气泡流动的结果，与其它数值方法所得结果相比，符合较好。