变比热激波关系式及其在高超声速进气道计算中的应用

本文根据空气热力性质表建立了空气热力性质的解析表达式,从气动热力学基本方程和激波前后比热比不同出发建立并推导出完全变比热激波方程组及其解法。对分别使用变比热激波关系式与定比热激波和等熵流关系式计算结果进行的比较表明,使用变比热激波关系式进行高超声速进气道气动热力学分析计算,能够减少定比热假定产生的偏差,提高计算准确性。      

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

二维亚、跨音速数字解的流线座标系方程组

本文导出了适用于亚、跨音速二维绕流问题数字解的流线座标系气体动力学方程组。由于在跨音速流场内可能存在激波,所以也导出了流线座标系二维激波关系式。最后附有算例。     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

在激波管中粒子阻力系数的实验研究

在激波管中用消光法测量了气固两相流动激波松弛区的粒子浓度。基于一维静止激波理论和理想“等效气体”模型,实验确定了粒子阻力系数。拟合实验结果,得到了阻力系数与粒子雷诺数之间的实验关系式:C_D=4300/Re~(1.57)。     

斜激波总压损失率极小值理论解与物理意义

以法向马赫数作为激波强度表征量,对斜激波关系式进行重新推导,得到了穿过激波总压损失率极小值的理论解。控制方程表达式为激波角对物面角的线性函数。依据斜激波总压损失率极小值解析公式,首先,绘制了针对超声速流总压损失率应用的楔形角-激波角-马赫数的斜激波效率图。其次,通过生成斜激波三维总压损失率等值线图,呈现了总压损失率在楔形角-激波角-特征马赫数空间上的分布规律。此外,利用斜激波效率图,揭示了等总压损失率条件下马赫数与激波角的对称双解现象。     

具有分流叶片跨声速离心叶轮激波及能量损失分析

对具有分流叶片跨声速离心叶轮,采用激波基本关系式对数值信息进行分析和判定,理论证实了叶轮流道中激波的存在,并确定流道内高速气流易在主叶片前缘、分流叶片前缘以及主叶片的吸力面与分流叶片的压力面通道位置形成激波.进一步对比分析了叶顶间隙和运行工况对激波位置分布特征的影响及引起的能量损失.结果表明:固定叶顶间隙下,随着流量减小,激波位置向流道上游移动;叶顶间隙导致激波线位置向流道下游偏移;分流叶片前缘处有最强激波并引起较明显能量损失.      

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

飞机RCM维修方式改革

就如何进行飞机维修方式改革,介绍了 RCM 的定义和基本情况、飞机 RCM 维修体系资料编制方法和内容、RCM 在国际上的推广和应用,并结合实践经验就飞机维修方式改革及RCM 维修体系资料编制方法和内容进行了论述。     

COMPUTATION OF COMPRESSIBLE FLOW PAST SLENDER WING－BODY USING FULL－POTENTIAL EQUATION

ＣＯＭＰＵＴＡＴＩＯＮＯＦＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＰＡＳＴＳＬＥＮＤＥＲＷＩＮＧ－ＢＯＤＹＵＳＩＮＧＦＵＬＬ－ＰＯＴＥＮＴＩＡＬＥＱＵＡＴＩＯＮＨｕａｎｇＭｉｎｇｋｅ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｏｄｙｎａｍｉｃｓ，ＮａｎｊｉｎｇＵｎｉｖｅ...     

用全位势方程计算细长翼-身组合体的可压缩绕流

用Ｈ－Ｏ网格、守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕大后掠细长翼－身组合体的可压缩流动。对ＡＧＡＲＤ－Ｂ翼－身组合体的计算表明,计算的自由流Ｍａｃｈ数可从亚音速直到低超音速。算出的升力与力矩特性与实验数据符合较好     

一种新型多维通量函数在求解平面叶栅跨音速流动的应用

采用一种改进的Ｒｏｅ格式对平面叶栅的无粘流动进行了数值模拟,并与Ｒｏｅ格式的计算结果进行了比较。采用四阶Ｒｕｎｇｅ－Ｋｕｔａ时间差分格式,空间三阶有限体积法。网格生成采用贴体网格近壁面修正的方法。计算结果表明,该格式在减少格式对网格的依赖性,抑制数值粘性方面具有明显的作用。     

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

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

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

跨音速三维机翼及翼身组合体非守恒型全位势流AF3算法

 本文对跨音速非守恒型全位势流,用Baker的快速有限差分法,改变其网络拓扑,研制出分析机翼及翼身组合体跨音速绕流的实用程序。其中包括:对“C-H”型贴体坐标网格提出AF3算法的具体实施方案取得了快速收敛效果;采用张量形任意曲线坐标上的非守恒全位势方程,使得程序中跨音速求解部分既能用于分析单独机翼又能用于翼身组合体跨音速绕流计算。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

两种椭圆型方程求源项方法在喷管内流场网格生成中的应用

用椭圆型方程生成二维网格，发展两种Hilgenstock求源项方法，实现对网格与边界间距及与边界正交的控制。     

突扩后突缩的三维流场计算

根据有关固体火箭发动机装药后端到喷管间的空间特征,对三维突扩后突缩内部流场计算进行了探索和研究。计算中成功地应用了SIMPLE算法解决复杂的三维流场计算问题,并在初场构成、源项离散及松弛系数的选取等方面都进行了调试改进。通过建立计算模型、划分网格、对方程及边界条件离散化和非线性处理,以及线性方程组求解等,计算得到了符合物理规律的数值解,指出了空间涡及二次涡的存在。     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。