用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

M∞=1附近轴对称Euler方程的跨声速计算

本文综合应用近似因子分裂法，矢通量分裂法及ＮＮＤ格式的离散思想，对轴对称Ｅｕｌｅｒ方程构造了一种修正的隐式格式，对球柱及多截锥旋成体在Ｍ∞＝１附近跨声速绕流情况进行了计算，与其它方法比较，本方法具有收敛速度快、精度高的优点。尤其对Ｍ∞＝１时的跨声速绕流情况，上述优点就显得更加重要。     

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

轴对称跨声速绕流的凸角流动分析

本文以锥-柱-船尾形旋成俸跨声速绕流为背景,讨论了轴对称踌声速绕流的凸角流动。给出了该流动特征、规律、并着重与相应平面二维Prandtl-Meyer型膨胀流动作了对比,指出了二者异同之处。     

机载高超声速飞行器中的轴对称旋成体气动布局研究综述

气动布局设计是机载高超声速飞行器总体设计的必要内容之一，轴对称旋成体布局作为一种机载飞行器应用的主流布局形式之一，有着极为重要的研究价值。本文从总体和气动的视角概述了三款使用轴对称旋成体布局的机载高超声速飞行器。分析了采用此种布局的飞行器在研发过程中可能遇到的难题和阻碍，并站在工程应用的角度给出了解决这些难题和阻碍的手段。最后，展望了轴对称旋成体布局的机载高超声速飞行器在未来的发展趋势，为国内机载高超声速飞行器的研制工作提供借鉴与参考。     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

船尾形状对旋成体马格努斯效应的影响

尖头旋成体和船尾形状是子弹、炮弹及火箭弹等抛射体上常用的布局形式。研究表明船尾布局具有减小底部阻力、增大射程的作用,但此时旋成体的马格努斯效应增大,对运动稳定性产生不利影响。为了解释这种流动现象,对三维旋转弹流场进行了数值模拟,对从亚声速到超声速下的旋成体马格努斯力和力矩进行了分析,重点对标准形状和船尾形状两种底部进行了比较。结果表明,相对于标准形状,在所有来流下船尾形状都起到了增大马格努斯效应的作用,并且马格努斯力和力矩与船尾角成正比。为了揭示其流动机理,选择代表性计算状态对两种布局马格努斯力矩系数分布、边界层厚度分布和边界层位移厚度分布进行了对比分析,结果表明,在亚跨声速下船尾马格努斯效应由绕拐角的加速流动引起,使当地压力系数幅值增大;在超声速下船尾马格努斯效应由船尾段的气流膨胀引起,使旋成体左右两侧的边界层位移厚度畸变增大。上述两种效应都使马格努斯力矩增加,对于亚声速流动来说,该效应发生在柱段与船尾段连接位置;对于超声速流动来说,该效应发生在连接点以后的船尾段上。当来流速度在声速点附近时,上述两种效应都可能发挥作用,使船尾形状的旋成体马格努斯效应大幅增加。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

用全位势方程计算机翼的亚声速,跨声速和超声速绕流

对大后掠小展弦比细长机翼,本文对机翼纵轴垂直的每一横流截面生成O型网格,形成对机翼流场的H-O型网格,用守恒型全位势方程、差分和隐式近似因式分解迭代算法计算绕机翼的可压缩位流。自由流可从亚声速直到低超声速的全部跨声速范围。本算法要求机翼前缘有大后掠角,后缘可稍许后掠或前掠。本文算例表明,所研制的计算程序已可提供工程实用。     

垂直于物面的横截面上流态的拓扑

对于旋成体有迎角绕流，本文研究了垂直于物面的横截面上的流态结果。利用解析方法，给出了物面轮廓线上、迎风和背风对称线上以及全横截面上截面流线方程奇点的分布规律以及奇点沿流动方向的演变规律，从而可定性地确定截面流态的结构。通过求解ＮＳ方程数值模拟了钝锥有迎角高超声速绕流流场，计算获得的垂直于物面的横截面流线的拓扑结果和理论分析完全一致。     

轴对称跨音速进气道外罩绕流的高效差分算法

 本文考虑轴对称跨音速进气道外流场的有限差分计算方法。采用守恒型位流方程、贴体坐标网格和精确边界条件,根据最佳收敛准则对轴对称情形设计出新的近似因式分解迭代格式,并将所设计的格式应用到计算皮托式进气口的轴对称跨音速流场。对几种典型进气口的计算表明,本文格式收敛快,计算结果与实验符合很好。     

用三维欧拉方程计算机身与机翼绕流

本文用Jameson的三维欧拉方程有限体积法、四步Runge-Kutta时间推进格式,计算机身和大后掠细长机翼的三维可压缩绕流。对钝头机身,用C-O型网格;尖头机身与大后掠细长机翼用H-O型网格。本文介绍钝头旋成体、带座舱前机身和三角翼绕流的计算结果,显示流场等值M线分布,计算压强系数分布与实验比较,以及三角翼大迎角分离涡等。由于用了隐式残值光顺等加速收敛措施,有效地减少推进步数,节省机时。     

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

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

三维任意机身的跨音速全位势流计算

 <正> 准确、有效地估计旋成体或任意截面机身在有迎角、有侧滑情形下的跨音速绕流特性,对飞机机身、导弹等外形设计有重要意义。已趋成熟的全位势流计算方法在翼型、单独机翼及翼身组合体绕流中得到广泛应用,并产生了诸如Holst的TWINGn和TAIR,Jameson的FLO系列等实用程序;但对旋成体有迎角绕流或单独机身的全位势绕流计算在文献中还很少见到。现采用守恒型全位势方程、代数法生成的贴体坐标网格和Holst的AF2高效差分方法对旋成体或单独机身有迎角、有侧滑时的跨音速绕流进行研究;并给出分析机身跨音速特性的计算程序是很有必要的。     

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

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

用四阶Runge-Kuta法和TVD格式求解两维跨超音流场

利用坐标转换给出两维任意曲线坐标系下雅可比矩阵特征值、特征向量及Harten型TVD格式的详细表达式, 并结合四阶Runge-Kuta法发展了一种新的、有效的求解超音、跨音流场的算法。数值计算表明, 本文的方法可有效地求解跨音流场。      

绕侧滑锥体,三维机翼与机身的超音速非线性全位势流计算

本文在Grossman的非守恒全位势方程的方法中,提出以固定网格为基础的头激波安装法,节省计算机时。计算结果与实验和欧拉方程解符合很好,与线化理论结果的比较,显示了线化理论的不足。     

跨音速三维多升力面气动设计

发展了一种跨音速多升力面的气动设计方法和设计程序,它基于已成功应用于亚、跨音速机翼设计和亚音速双翼面设计的"余量修正迭代"概念.当升力面上出现超音速区和激波时自动引用迎风格式对控制方程进行修正.开发了一系列接口程序,包括目标压力设计程序.由此气动设计程羊、TAU程序以及相应的接口程序建立了跨音速多升力面气动反设计软件系统.用两个鸭翼-机翼构型验证了设计方法和设计程序,结果表明在高跨音速下设计迭代有很好的收敛性.     

叶轮机平面叶栅的气动设计及统流计算

通过全速势方程的数值积分，利用ＡＦ２格式，给出了一种用于叶轮机中叶栅型面的气动设计和绕流计算的便捷方法，对于稳态亚音速、跨音速及超音速流动均适用。与目前普及很广的、在超音区应用人工压缩性的方案相比，该算法有明显的优势。