非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

再人体复杂外形无粘流场数值计算研究

本文在文[1,2,3]的基础上研究了非定常Euler方程的推进迭代方法,并将无波动、无自由参数的耗散差分格式(下称NND格式)发展为隐式、迎风格式,用以计算有差曲控制翼的再入复杂飞行器的育攻角、有侧滑角的超声速无粘绕流流场。其结果是准确的。本文的方法具有以下特点: 1.可以在同一计算程序中实现亚跨超流场的计算:在超声速区域,采用推进技术;在亚跨声速区,则自动选择推进迭代求解技术。 2.NND格式不仅具有公式简单,无自由参数的特点,而且过激波时满足熵条件,不会产生非物理解,也不发生振荡。 3.隐式差分格式的无条件稳定和Gauss-Seidel迭代的快速收敛,计算省时。 本文的方法曾用于烧蚀头部凹陷外形和航天飞机简化外形的计算,能自动选择迭代和推进区域。可以预言,如计算机条件许可,本方法可以用于大攻角的无粘流场计算。     

高阶精度非线性格式WCNS-E-5在二维流动中的应用研究

本文采用具有5阶精度的加权紧致非线性显式格式(WCNS-E-5)对定常与非定常二维流动进行数值模拟,研究表明该格式对各类间断有很好的分辨捕捉能力,而且对强间断如激波的计算,即使在高马赫数与高雷诺数条件下它仍具有很好的收敛性与可靠的计算结果.此外,WCNS-E-5在粗网格条件下也体现出优越性.类如WCNS-E-5的高精度激波捕捉方法将为以后开展湍流数值模拟工作提供坚实的技术保证.     

三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

鼓包对舵的气动特性影响研究

 运用自适应八叉树结构直角网格 Euler算法对舵的流场进行了模拟。与无鼓包的舵相比较,带有鼓包的舵其气动特性,尤其是阻力有所变化。给出了直角网格生成方法以及 Euler方程解法。结果显示,本算法很准确地捕捉了激波,真实地模拟了该流场。在此基础上讨论了鼓包及其位置对舵的气动特性的影响。     

水洞收缩段流场的数值模拟及优选

介绍了一种水洞非对称收缩段的数值模拟方法和程序,计算采用有限体积预估修正格式求解非定常欧拉方程.结合具体水洞设计对几种典型收缩曲线的三维流场进行了计算和分析,给出了沿轴向的压力分布,出口截面速度分布及分布的均匀性,选出最佳曲线,为水洞的设计提供依据.结果表明:本文的计算方法和程序对于水洞或风洞的三元收缩段,是有价值的分析和辅助设计手段.     

基于刚体自由转动的连续推力消旋方法研究

针对空间失效卫星或废弃碎片转速过快,难以对其直接处理的问题,提出了一种基于刚体转动的连续小推力消旋方法。分析了轴对称消旋目标的运动特性,设计了消旋方法。将轴对称消旋目标的消旋方法推广到一般的消旋目标上,通过变换欧拉方程得到消旋目标各角速度分量的变化趋势和最终状态。通过一般消旋目标的仿真实验验证了该方法的可行性。结果表明:按给定顺序和方向施加单轴外力矩,能够有效降低消旋目标的转速。区别于传统的控制方法,最优控制方法需要精确测量消旋目标的实时状态信息,而分步消旋方法只需分辨角速度分量的正负,且控制力矩只需反向变化,推力器无需摆动,因此降低了对控制系统的要求。     

二维欧拉方程在非结构网格上的自适应有限元算法

发展了一种易于和自适应网格加密技术相结合的非结构三角形网格生成技术。这种方法的特点是域内节点可任意规定;根据欧拉方程解的误差估计,在流场某些区域增加节点后,能重新快速生成新的具有高质量三角形单元的非结构网格。欧拉方程采用Ｇａｌｅｒｋｉｎ空间离散和Ｒｉｃｈｔｍｙｅｒ两步显式时间推进相结合的有限元法求解,并以沿流线熵的变化量作为物理判据进行网格自适应加密,提高解的精度。     

基于非结构动网格的跨声速非定常流场求解

对用于非结构动网格生成的弹簧近似方法进行了研究.通过采用顶点弹簧方法,分析研究了弹簧倔强系数的取值,同时通过引入防扭转、防挤压倔强系数和边界修正,对标准弹簧近似方法进行了改进.采用改进后的方法并通过耦合求解了基于ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程,模拟了做谐和振动的刚性无后掠机翼及某翼身组合体的跨声速绕流.结果表明,改进后的方法可以大大提高网格变形能力和网格质量,且计算结果与实验结果吻合良好.     

折叠降落伞展开过程研究

为更好地从机理上研究折叠降落伞的展开过程,文章首先采用基于初始矩阵修正的直接折叠法建立了降落伞的折叠模型,并对该模型使用任意拉格朗日-欧拉法模拟了无限质量情况下不同风速、不同引导伞作用下开伞情况。计算结果证明了引导伞经验设计的有效性,但同时发现个别情况下经验设计过于粗糙,容易导致出现不良充气状态。文章的研究内容对降落伞设计、优化具有指导意义。