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

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

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

高超音速飞行器前体优化设计方法研究

基于响应面法对高超音速飞行器前体进行气动优化设计,流场求解采用欧拉方程,并结合非结构动网格技术,便于适应复杂外形而且在每次外形改变后不必重新生成网格,有利于保证计算效率,响应面模型采用二次多项式来构造,数值结果表明由本文设计的前体/进气道布局结构简单,能获得更高的压缩率和总压恢复系数,以及更高的升阻比,该方法在多约束多目标条件下进行气动优化设计,设计质量较高,有一定的工程应用价值.     

雨流计数法及其在程序中的具体实现

根据雨流计数法的规则和在实际中应用的体会，介绍了雨流计数法在程序中实现的具体方法。它适合用各种语言编写，在计取循环数时采用的“四点法”使程序的实现比用其他方法更加准确可靠，简单明了。     

顺序函数法求解二维非稳态热传导逆问题

基于有限控制体积法对二维非稳态热传导问题的数值模拟，导出了处理二维非稳态热传导逆问题的顺序函数法，并用该方法来对一典型的圆环域边界条件反演问题进行了反演计算，结果表明，顺序函数法是求解二维非稳态热传导逆问题的有效方法，当测量噪声比较小时，顺序函数法能得出较高精度的反演结果；当测量噪声增大时，反演结果仍能较好地再现出精确解在空间和时间方向上的变化趋势，具有一定精度。     

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

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

吸气条件对圆柱非定常分离流影响的数值研究

本文根据二维可压非定常Ｎ－Ｓ方程，基于Ｂｅａｍ ａｎｄ Ｗａｒｍｉｇｎ及Ｓｔｅｇｅｒ的方法，在圆柱绕流出现周期分离后，施加吸气条件，研究气对于圆柱非定常绕流的影响，计算结果表明，在恰当地区域施加恰当的吸气条件，可以希望通过减小压差阻力来减小总阻力，吸气还有抑制尾迹中非对称流动的作用，从而使圆柱所受的垂直于来流方向的升力的变化幅度减小，变化频率减小，但是吸气会增加磨擦阻力，过份的吸气会导致摩阻增大太     

滑流区中三角翼大攻角气动力特性的实验研究

一种以后掠75.7°薄三角翼为主要特征的典型航空航天飞行器模型,在激波管风洞马赫数为11.9和15.4两种条件下,攻角范围20°～50°,用模型自由飞方法测量了它们的轴向力系数、法向力系数和俯仰力矩特性。相应的实验雷诺数分别为3.19×10~4和1.64×10~4,这两种流动条件均属于稀薄气流的滑流区。 实验结果表明在M_∞=11.9和15.4两种条件下,两种剖面外形模型的升力系数和阻力系数均随攻角加大而递增,其变化规律有很好的一致性,且对马赫数并不敏感;但从体轴系来看,不仅两种模型的轴向力系数不同,而且因粘性干扰的缘故,同一模型A在M_∞=15.4时比M_∞=11.9时有相对较大的轴向力系数,但两者随攻角变化的规律一致,且当α>45°时接近牛顿值。此外,实验表明两种模型的压心系数随攻角均没有明显变化。     

钝锥有攻角分离流动的数值模拟及其分析

本文定性分析了开式分离的性状,并对钝锥有攻角超声速绕流的开式分离作了数值模拟。分析指出,开式分离可能存在两种形态,第一种分离线的起点为正常点,第二种分离线的起始为鞍、结点(包括螺旋点)的组合。对于第一种形态,分离线的起点是横向分离的起始点,除分离线外,分离面上的流线不是从分离线的起点发出的。对文中计算的情况,流动属第一种开式分离。计算证实了定性分析的结论。计算和分析均指出,对第一种开式分离,在分离的起始区域,分离流面尚未卷曲,但在下游,则变成卷曲面。文中还研究了围绕物体的流管在分离诱导下的变形情况。     

简谐激励下弯曲薄板的非线性运动方程

推导了两端固定和两端简支的弯曲薄板在其支承简谐激励下的运动微分方程,并利用Galerkin方法对其进行了离散化处理。