关于跨音速定常小扰动平面势流线松弛的稳定性和收敛性的讨论

引言 对跨音速势流方程采用混合差分线松弛求解,常碰到稳定性和收敛性问题。差分方程的舍入误差,若在求解过程中不会增长,则格式稳定。当选取的步长趋于零时,差分方程的解趋于微分方程的解,则格式收敛。要得到所要求的解,必须使差分格式既稳定又收敛。由于跨音速势流方程的非线性,严格的稳定性与收敛性证明十分困难。大量的计算实验指出:在松弛求解中,计算是否稳定,是否收敛和当收敛时的收敛速度,     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

机翼无粘跨音速绕流数值计算

本文介绍机翼无粘定常跨音速绕流的一种数值计算方法。选用全势方程作为数学模型。用适当的坐标变换把后掠翼变成矩形翼,再把无限物理空间变成有限的计算空间。在计算空间中,用混合型有限差分格式对方程离散化。差分方程形成的代数方程组用线松弛迭代求解。逐次加密计算网格,使计算较为经济。     

机翼跨音速无粘绕流数值计算

本文介绍机翼跨音速定常无粘绕流的一种数值计算方法。选用精确速势方程作为问题的数学模型。在直角坐标系中,经过适当的坐标变换,先把后掠翼变成矩形翼,再把无限的物理空间变成有限的计算空间。在计算空间中,用有限差分格式(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用线松弛迭代求解。     

用三维N－S方程计算钝体绕流

本文采用有限体积方法求解三维低速Ｎ－Ｓ方程，研究了三维钝体和类汽车物体的绕流。用中心差分和迎风差分的混合差分格式构造代数方程组，采用了隐式和显式相结合的方法使主对角元素占优以提高迭代求解过程的稳定性。利用压力校正算法，采用面松弛迭代来实现方程的求解。湍流的计算采用了标准的双方程模型。本文给出了钝体和类小汽车外形的计算结果。钝体的计算结果和实验结果进行了对比，与其它方法的计算结果也进行了对比，它们之间具有合理的一致性。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

烧蚀外形方程差分计算方法研究

分析了烧蚀外形方程的性质,研究了人工粘性项对烧蚀外形计算结果的影响,发现人工粘性在抹平非物理波动的同时,会严重影响计算精度.借鉴于CFD中NND格式的思想,构造了求解烧蚀外形方程的无波动、无自由参数的耗散差分格式,得到了满意的计算结果.     

求解双曲型守恒律的五阶松弛格式

给出了一种求解一维双曲型守恒律的五阶松弛格式.该格式以五阶WENO重构和显隐式Runge-Kutta方法为基础.本文格式保持了松弛格式简单的优点,即不用Riemann解算器和计算非线性通量函数的雅可比矩阵.用该格式对一维Euler方程进行了数值试验,并与三阶和四阶松弛格式的计算结果进行了比较,结果表明本文的格式具有更低的数值耗散和更高的分辨率.     

环境压强对电弧喷射推力器性能的影响

为了详细研究环境压强对电弧喷射推力器性能的影响,对低功率氮氢电弧喷射推力器在不同环境压强下的工作性能进行了比较计算。采用二阶精度无波动、无自由参数的耗散差分格式(NND格式)和经典龙格库塔方法求解耦合电磁源项的N-S方程组,采用高斯-塞德尔超松弛迭代方法数值求解电磁场方程。计算结果给出了不同环境压强下推力器的推力、比冲、推进效率和出口平面冲量密度分布。结果分析定量地阐释了环境压强对推力器性能的影响,能为推力器的实验研究及空间应用提供参考。      

任意翼型跨音速绕流守恒全速势方程的快速有限差分解法

一、引言 目前求解翼型跨音速绕流全速势方程比较好的有限差分解法,主要有二种:一种是Garabedian,Korn和Jameson所发展的方法。这种方法先将物理平面上翼型外部区域保角转绘到圆内,然后在圆内交替使用快速直接解法和线松弛技术,求解全速势有限差分方程。由于圆心对应物理平面的无限远点,速势在该点为无限大,为了避免这种奇性,Jameson引进了扰动速势以除去包含这种奇性的自由流速势,使主管方程     

一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

叶轮机械中三维跨声速流动的正,反混合问题

文中提出了跨声速叶轮机械中三维流动的一种正、反混合问题。首先，将作者于１９８３年的平面叶栅跨声速流动的正、反混合问题推广到三维流动情况，其中反问题的特点是既能计及叶片的气动性能，又能在一定程度上兼顾强度、冷却、工艺等方面的要求，蕨，由于应用了作用设计的一种适合叶轮机械中三维流动特点的广义ＶｏｎＭｉｓｅｓ坐标系，克服了反问题所特有的不定求解域边界的困难，并将三维正问题分析、叶片的完全三维设计以及三维     

跨音速叶轮机械叶片的完全三维设计方法

设计方法由一种平面叶栅跨音速流动反问题推广得到,其特点是既计及叶片气动性能,又能在一定程度上兼顾强度、冷却、工艺等方面的要求。由于应用了一种适合叶轮机械中三维流动特点的广义ＶｏｎＭｉｓｅｓ坐标系,克服了反问题所特有的求解域边界不定的困难;并使正问题计算和反问题设计统一在同一算法之中。典型算例的结果表明：提出的跨音速叶片完全三维设计理论模型是合理的。     

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

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

跨音涡轮中三维流动的数值试验

 本文从积分型的基本方程组出发,以我们的关于平面叶栅和任意旋成面叶栅的工作为基础,将文献（3、4）中的计算方案推广到三维流动情况。编制了相应的计算机程序,进行了一些数值试验,得到了初步的计算结果。 数值试验表明,拟三维计算结果可以导致相当大的误差,从而说明了完全三维流动计算的必要性。另一方面,如果取2000个左右的计算点,则在国产的013机上算一个例子,约需100分钟。从实际应用的观点来看也是可行的。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

跨音速透平叶栅速度面设计方法的研究

 低能量损失系数跨音速透平叶栅的设计和实现是提高跨音速透平性能的关键问题之一。本文给出了基于有限面积法的跨音速透平叶栅速度面反问题设计方法。首先推导了基本方程的对称形式,用于上游奇点处理,并对速度面上不连续边界条件作了分析和数值验证说明。然后详细叙述了流函数场的有限面积解法包括主方程的积分变换,基本有限面积格式和边角有限面积格式处理以及解精确度的粗略分析。所编制的程序简明快速。计算结果表明解相当稳定并且不需要在音速线附近等处加密网格或作其他处理。算例也指明音速线上临界奇点位置的选取,对流场及叶型型线形状有一定影响,应在设计时注意。     

带操纵面机翼的非定常跨声速流数值解法

本文介绍一种带操纵面机翼的非定常跨声速流的有限差分计算方法。采用的方程是三元非定常跨声速修正小扰动位势方程,使用时间积分法,求解格式是近似因式分解交替方向隐式(ADI)格式。使用这种方法计算了F-5机翼在来流马赫数为0.9和0.925时的定常气动力和操纵面振荡的非定常气动力。计算结果与国外的NLR试验结果和XTRAN3S方法的计算结果进行了比较,表明计算是成功的。     

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

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