欧拉方程组的隐式近似因子分解差分方法的研究

本文对求解二维欧拉方程组的隐式近似因子分解差分方法进行了简要的分析,并利用它计算了简单直扩张通道内的亚音速、超音速和跨音速管流、绕无限长拐角和对称棱形翼型的超音速绕流以及绕NACA 0012翼型的亚音速、跨音速和超音速绕流等典型定常流动。计算结果表明,它具有良好的适应性和精度,可以正确地捕获强激波和自动满足Kutta条件,Courant数可以取大于1的数值,从而加速收敛。分析和计算都是在贴体坐标下进行的。方法本身可以推广到三维流的计算。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

CE／SE并行计算研究及在Euler方程中的应用

时空守恒元和解元方法(CE/SE)是一种高精度的数值求解方法.采用改进的CE/SE方法和并行分区技术,针对非结构网格,发展了一套求解二维Euler方程的并行程序.对NACA0012翼型的亚跨声速流动和多段翼型复杂绕流进行了分区并行求解,多区并行计算的结果与实验的结果吻合较好,表明了本文研究方法的正确性和有效性.同时文中的分区方法能实现各处理器之间的负载均衡,有效地节省计算时间.     

飞机／进气道系统内外流场计算

为了对飞机/进气道系统内外流场进行数值模拟,本文提出并研制了计算定常(或非定常)外压式飞机/进气道系统内外流场的计算方法和程序。它适用于任意几何外形的飞机/进气道在亚音速、跨音速或超音速来流条件下的内外流场和气动参数计算。求解的主管方程是在一般曲线坐标下的 Euler 方程、简化/雷诺平均的 Navier-Stokes方程。用隐式时间差分和近似因式分解方法离散和求解.针对飞机/进气道的复杂流场,本文采用了贴体网格生成技术和分区域方法.对于湍流流动,把 Baldwin-Lomax 湍流模型推广到飞机/进气道内外流场计算并使用了壁面律计算公式。应用上述方法对绕翼型的亚、跨音速、绕可压缩拐角的超音速湍流以及飞机进气道系统内外流场进行了计算,计算结果与实验数据和其它计算结果符合较好。     

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。     

M=1音速流翼型设计计算

来流马赫数为1.0的翼型绕流型态不同于其它跨音速时的流动型态。在数学上,它的远场边值呈奇异性。本文利用Chaplygin转换推导出精确的速度面流动方程,又根据小扰动理论,采用Tricomi方程的奇点近似解析解建立了对应在速度面上来流奇点的边界条件。在用有限差分法进行数值计算时,创建一种新的差分格式,并在不同速度区采用不同座标系统进行线松驰迭代计算,计算表明,其结果是稳定的和收敛的。本文用Spline方法确定出五个仟意的速度面上的边界曲线,经计算,在获得五个音速流流场的同时,也得到五个封闭的音速流翼型。     

正交流线坐标系中计算二维跨音速流的流线迭代法

本文提出了流线迭代法,用以实现文献[7]所提出的正交流线坐标系中的二维定常无粘气体动力学方程组的数值求解。 此方法的基本做法是:从物面起到流场远边界逐条流线地求出流动速度值,同时求出下一条流线的形状;按此方式对全流场进行迭代。迭代的收敛性很好。 用本文方法在零攻角和跨音速绕流的条件下计算了双圆弧翼型、NACA 0012翼型和R.A.E.101翼型。双圆弧翼型结果与实验数据符合较好。钝头翼型结果与其它方法结果相比是令人满意的。     

跨音速翼型多参数多约束优化设计方法

本文通过采用EXTREM数值优化方法并基于求解全速势方程的跨音速粘流翼型计算方法, 研究发展了一种应用形状函数修形的跨音速翼型多参数多约束数值优化设计方法, 应用该方法可以从普通低速翼型和超临界翼型出发通过数值优化后得到在跨音速区阻力系数最小化的翼型几何外形。设计表明, 该方法具有收敛快、调用目标函数次数少等优点。      

跨音速翼型气动优化设计方法

 提出了一种有效的跨音速翼型气动优化设计方法。翼型的流场解由欧拉方程的数值解提供。带约束条件的优化计算分别采用了直接法和间接法两种优化算法。算例结果表明,本方法提供了一种跨音速翼型改型设计及新翼型设计的有效工具。     

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

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

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

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

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

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

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

高负荷跨音压气机转子ATS-3气动设计

ATS-3的设计工作是在已有的设计较为成功的风扇压气机 ATS-2的基础上 ,进行细致的优化 ,主要是想提高失速裕度 ,也为高负荷、高压比跨声速压气机的设计提供更为直接、有用的设计数据。优化设计的重点是在根部应用可控扩散叶型的思想进行二维反问题优化、并对整体的安装角进行调整。其特性的数值模拟结果表明 ATS-3较 ATS-2有更大的失速域度 ,效率有所提高 ,但设计点的压比和流量略微有所降低     

计算跨音速三维机翼非定常载荷的核函数法

 本文给出了在跨音速流动中,计算三维机翼非定常载荷分布的核函数法和计算结果。此方法是基于线化升力面理论,利用在亚、超音速情况下,联系载荷分布与下洗的积分方程的核函数,采用局部线化法来处理跨音速的非定常绕流问题。在计算中引进了线分布的冲波偶极子,并嵌入适当的正冲波边界条件。 我们从统一的超、亚音速核函数表达式出发,对亚音速区数值积分作了新的简便的处理,对于亚、超音速区相互的诱导下洗给出了自己的数值积分过程。我们算了三个例题,和文献上的数值符合较好。     

气体运动论BGK格式的翼型绕流数值模拟

以气体运动论BGK(Bhatnagaar-Gross-Krook)格式为基础,通过改进非均匀网格下的通量计算方法,把气体运动论BGK方法应用到二维贴体网格上模拟了可压缩的无黏和黏性流动。在重构阶段,以交界面周围4个网格单元为一组,先把每个网格单元的宏观量(密度、动量和总能量)转换在相对于交界面的坐标系下,然后使用van Leer限制器插值求出交界面的宏观量,在求解阶段再使用BGK理论求出流过交界面的通量。通过上述方法使贴体网格和均匀网格,在计算格式上一致,与其他求解二维贴体网格的气体运动论BGK方法相比,该方法相对简单,计算量减小。结合贴体网格通量计算和气体运动论BGK格式方法理论,分别对亚、跨、超声速下的NACA0012翼型绕流进行了数值模拟和对比研究,所得可压缩无黏流动的数值结果与北大西洋公约组织航天研究与发展咨询组(AGARD)的结果吻合良好。本文还模拟了低雷诺数状态下的可压缩黏性流动,所得结果与文献中的数值计算结果吻合良好,证明了该方法在数值模拟及激波捕捉方面具有鲁棒性,为在复杂流动中的应用提供了一种简单的途径。     

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

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

计算流体力学中的网格生成方法及其应用

讨论了多种网格生成方法及其流场的数值模拟。给出了一种生成现代双立尾战斗机外形的分区对接网格和Ｅｕｌｅｒ方程分区解。给出了一种带副翼偏转的三维机翼分区网格及绕流的Ｅｕｌｅｒ方程解。以二维多段翼型为对象分别讨论了分区重叠结构网格方法、非结构网格方法和自适应的笛卡尔网格法。