有弓形冲波影响情况下的层流附面层方程的数值解

本文给出了一种轴对称层流附面层方程的数值解法。利用Levy-Lees变换,将物理平面上的附面层方程变换到新的坐标平面内,再利用三点隐式有限差分法来解这组方程。在考虑熵吞影响的情况下,本文发展了一种与国内外惯用的流线吞嚥法不同的不需迭代即可给出附面层外缘参数的方法,从而使计算时间缩短了三分之二。本方法给出的结果同其他作者的数值解和实验结果作了比较,结果是令人满意的。本文还给出了驻点热流在等熵和变熵情况下随雷诺数和马赫数变化的计算结果,并同实验结果作了比较。     

二元全速势跨音速方程的数值解

本报告提供了二元全速势跨音速方程的一种数值解法。本方法的特点是放弃了繁杂的曲线座标系,而在简洁直观的直角座标系中,利用Jameson的旋转差分格式,并在物面边界上采用速势内嵌场,建立全流场的差分方程。二元翼型的样例计算表明本方法和其他全速势方程的数值解法具有同等的精度。对物面外形作一次附面层位移厚度修正后的算例表明计算结果和试验结果符合较好。 本方法比起采用曲线座标系的方法,更易推广到三元流。     

绕低频振动旋成体的跨音速非定常流核函数方法

本文对无粘跨音速流中作俯仰振动的细长旋成体,给出了一种计算表面非定常压力分布和动稳定性导数的近似方法,即用时间线化近似,在引入局部线化假设后,将跨音速小扰动方程化为一个积分方程,并用精确物面边界条件求数值解。本方法算例与其它理论和实验结果作了比较,表明用它可以得出较合理的结果。     

弹性飞机空气动力及动态特性的数值计算和研究

本文给出了一种弹性飞机动态特性的数值计算和分析方法。文中考虑了飞机受扰后弹性振动与飞机总体运动的耦合作用,并着重研究了前几阶低频结构振型对飞机总体运动的影响。本方法采用了格林定理所导出的速势积分方程方法计算定常和完全非定常空气动力。采用特征根迭代法求解弹性运动方程,只需经较少次数的迭代即可收敛。文中附有算例,其计算结果是合理的。     

叶轮机械全三维气动杂交问题的数值求解方法

本文介绍了一种求解叶轮机械全三维气动杂交问题新的数值方法。本方法使用了非正交曲线座标系和非正交速度分量的势函数方程, 并构造了守恒形式的近似分解因式格式, 通过不断修正叶片通道宽度, 使叶面上速度达到预先要求的速度分布。根据本方法使用 FORTRANTT编写了计算程序, 并对单转子压气机叶片进行了修正计算, 其结果十分令人满意。      

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

烧蚀传热体有限元法

本文介绍一种适合于烧蚀传热特点的简便和实用的有限元法。首先从能量守恒原理出发,建立带有导数的有限元素方程,然后介绍了采用伽辽金(Galerkin)法对时间域进行有限元离散而得到的各种格式。 为了解决烧蚀移动边界问题,采用了在动坐标系上建立有限元素方程。 用本方法对典型问题进行了计算,并与理论解、差分解进行了比较,对方法作了初步讨论。本方法应用实例是高速飞行器轴对称端头再入烧蚀体的传热计算。     

包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

推进剂供应管路内液体瞬变流一维有限元计算

将液路计算中偏微分方程模型化成有限阶数的常微分方程组,由一维液体瞬变管流方程的特征线差分格式,提出了一种计算管腔互联结构形式的液体推进剂供应管路中瞬变流的一维有限元方法,它兼有特征线方法的特点。数值计算表明对单根等截面圆管阀门关断问题的计算结果与特征线方法计算结果吻合。利用本文方法计算了一个模型发动机的脉冲工况动态过程,取得了很好的效果。     

考虑粘性相互作用的跨音速机翼分析

本文将一些已有的三维全位势方程计算方法和一种三维积分边界层方法作了弱相互作用的耦合计算。计算了三种机翼并与已有的试验数据作了比较。     

面向稀薄流非线性本构预测的机器学习方法

稀薄非平衡流域内连续介质假设已经失效,主要围绕Boltzmann方程及模型方程对稀薄非平衡流开展理论与计算研究,统一气体动理论格式（UGKS）是其中一种代表性方法。在稀薄非平衡流数值模拟中,Navier-Stokes （N-S）方程连续介质假设已经失效,不能有效描述流场非平衡特征。UGKS方法虽然计算精度高,但速度空间离散导致计算效率低下,多维高速条件下数值计算难以开展。基于数据驱动的思想,在N-S方程与UGKS方法的研究基础上发展出了一种稀薄非平衡流非线性本构关系求解方法（DNCR）。该方法以N-S与UGKS求解器获得的流场数值模拟计算结果作为训练数据集,基于流场特征参数采用极端随机树算法生成机器学习模型,对预测流场中线性黏性应力项与热流项进行非线性修正,并耦合非线性本构关系求解宏观守恒方程得到目标状态稀薄非平衡流动数值解。针对DNCR方法中所采用的机器学习方法-极端随机树模型,通过二维顶盖驱动方腔流算例对高维非线性建模涉及的特征参数选取、参数调优开展了相关验证工作,选取若干典型状态对极端随机树模型的泛化性能开展研究,并评估了相关模型与方法的计算精度与计算效率。     

扫描有限元法求解二元喷管的三维可压粘性流场

发展了一种二元喷管流场计算的三维Ｎ－Ｓ方程扫描有限元方法,给出了Ｎ－Ｓ方程中二阶导数项新的离散方法,以及更完善、简洁的扫描有限元方程的表达式。开发了相应的计算程序,用所发展的方法对二维平面收扩喷管、两种大宽高比的二元喷管及其圆变方转接段的三维粘性流场进行了计算,计算结果表明本方法具有较宽广的适用范围,计算结果良好     

简化Navier-Stokes方程的数值解法

众所周知,L.Prandt1的边界层理论对于早期的粘性流研究和应用起了重要的作用。时至今日,边界层方程的数值计算仍然是一个十分活跃的方面。但是,目前实践中提出的许多新课题(如激波与边界层的相互作用问题,分离流问题等)远非边界层方程所能描述,而必须借助更复杂的微分方程。在这方面,六十年代初期提出的简化Navier-Stokes方程(SNSE)是对古典边界层方程的重要改进,它适用于十分广泛的流动问题。我们以高超音速粘性绕流向题为主,对二维及三维简化NS方程的数值解法作了一些研究,本文不打算全面叙述所有计算问题,而集中对计算方法方面的问题作一小结,以供参考。     

多级轴流压气机特性的一种计算方法

本描述了多级轴流压气机特性计算的一种方法，流场计算采用流线曲率法解完全径向平衡方程。本在特性计算中引进了环壁附面损失模型，考虑了叶片去亏损力，二次流和径向间隙漏流等影响。采用了新的流场计算模型，流场计算稳定，收敛较快。并对计算结果与试验数据进行了对比，认为这种计算方法可以用于轴流压气机特性计算。     

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

积分方程法计算翼型的跨音速绕流

从跨音速小扰动方程出发推导积分方程的过程中,本文用任意形状的封闭曲线C_Q（其极限趋于零）挖去奇点Q,最终得到无奇性（指无穷奇性,不包括Cauchy奇性）的积分方程。 对于跨音速流中的圆头翼型的前缘问题,提出了一种解决办法。 证明了Nixon给出的反演公式对于超临界有激波的小扰动流动也成立。 关于积分方程法中的人工粘性方法,对Sachdev和Lobo提出的方法做了改进。 最后给出了NACA0012翼型在有无升力和有无激波各种状态下的计算结果。比较表明,本方法的计算结果与其它方法的计算结果符合得较好,且计算量很小。     

管道振动问题中一种辛变换的计算方法

以符合哈密顿(Hamiltonian)特性的管道振动方程的简化问题给出了一种具体的计算方法.为了保持原方程哈密顿结构,这种简化过程中使用的矩阵必须是辛矩阵(sympletic matrix).针对管道振动方程给出一种辛矩阵的计算方法,该方法利用了哈密顿方程线性部分的特征空间,用特征空间的一些性质来完成辛矩阵的计算.     

过渡态带隔热涂层的燃气涡轮叶片热状态计算方法研究

过渡态条件下热端部件热状态计算在评估发动机性能、强度寿命及内部间隙等方面占有十分重要的地位。参考了俄罗斯关于过渡态涂层叶片的1种计算方法,但对其中的数学处理作了更广泛的集成和推广,使得该方法能更充分地考虑过渡态变化过程的复杂性以及换热边界条件的适时变化等问题;以该方法编制了带涂层热端部件的工程计算程序,将计算结果与国外数据进行了对比,吻合较好。     

导弹空中挂飞分离试验仿真计算

导弹设计时必须考虑其投放过程中与载机分离时的安全性问题.利用计算机数字仿真技术进行导弹分离安全性分析是一种非常经济、有效的方法.在某型导弹挂飞试验中,我们以Euler方程和刚体六自由度运动方程为理论基础,完成了对该某型导弹在试验选定的各种速度、迎角下的投放分离轨迹进行了计算仿真,取得了非常满意结果,载机实际试飞验证表明,试验数据与理论计算结果十分吻合.该结果也证明了Euler方程和刚体六自由度运动方程在工程应用上是非常有效的,特别是对于具有一定长细比结构体的空中运动轨迹分析计算是非常适用的.     

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。