用单调相关格式计算跨音速定常和非定常翼型绕流

在跨音速非定常势流计算方面,人们一直在努力寻求处理包括激波运动的非线性问题的能力和计算效率的提高。Ballhaus和Goorjian（1977年）用交替方向隐式差分法（ADI）成功地直接求解了时间域的非定常跨音速小扰动方程（此法称为“LTRAN2”）。文献[5]引用坐标变换也做了类似LTRAN2的工作。     

翼-身组合体跨音速绕流全位势积分方程数值计算

 用积分方程方法求解 Prandtl- Glauert算子表示的全位势方程 ,并计算了翼 -身组合体跨音速绕流。用 Murman- Cole差分格式计算空间场源强度 ,以捕捉激波。计算结果与相应的实验结果符合良好     

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

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

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

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

后掠机翼跨音速绕流的二级近似方法

本文采用TSDH方程计算三维后掠机翼的跨音速绕流,考虑了适用于机翼钝前缘的前缘边界条件和前缘速势方程。采用Jameson格式在不等步长格网中的推广形式,把TSDH方程离散化为差分方程组。然后,在整个计算空间内布置稀网,在机翼附近再布置密网,进行稀密网的交替迭代,以加速收敛和提高计算精度。对ONERA M6机翼的超临界无激波和有激波情况的计算表明,TSDH解与FVP解和风洞试验符合良好。     

关于“跨音速绕流的边界积分法”

文献[1]是一篇参加国际会议的文章,由第二作者在会上宣读。当时,我们只准备在会议上透露一信息:我们已可建立可压缩非线性位流方程的边界元解法。对方法的中心思想无意细说,拟在此项研究工作结束后再正式发表专文系统叙述。所以文献[1]写得比较简单显得晦涩。现仅就文献[2]和文献[3]对文献[1]所提的问题作如下补充说明,我们今后发表的正式文章中对此将作更全面的论述。     

N-S方程数值模拟弹体不可压脱体涡绕流

应用在连续方程中增加压力对时间导数项的拟压缩性方法 ,数值求解三维定常不可压N- S方程。采用隐式 Beam- Warming近似因式分解格式求解。以修正的 Baldwin- Lomax代数模型模拟湍流流动。对绕正切拱形旋成柱体的大迎角脱体涡流动进行了数值模拟 ,并分析了层流粘性和湍流粘性对脱体涡及二次分离涡的影响。在头部施加微小的非对称射流扰动可得到与实验相符的非对称脱体涡     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

三维跨音速流有限差分计算中的某些问题

一、基本方法简介 1.三维跨音速小扰动方程 柱坐标下的三维跨音速小扰动方程为     

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

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

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

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

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

振动翼型跨音速非定常Euler方程数值解方法

通过坐标变换使变化的物面边界条件及无穷远边界条件变为固定的有限范围内的边界条件。因此可在不增加计算机时的情况下,边界条件可以精确地满足。由于主控方程是Euler方程,因此,凡适合于Euler方程解的问题（跨音速、大扰动、非定常、高简化频率）均可应用此法。同时,在计算过程中采用贴体、正交网格与可控网格技术获得最佳网格分布。算例计算结果是令人满意的。     

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

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

计算绕薄翼型跨音速非定常流的积分方程法

 本文从小扰动方程出发推导出绕薄翼型的跨音速非定常流动的积分方程,为使此方程适用于具有激波的流场引入人工粘性,并对其数值求解。通过算例讨论了人工粘性和计算域大小对计算结果的影响,与实验及其它数值计算结果比较,表明本方法的准确度令人满意。计算量小,收敛性好是此法的特点。     

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

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

双三角翼前缘剖面形状对涡运动的影响

采用数值模拟和理论分析相结合的方法,研究了前缘剖面形状对双三角翼涡运动的影响,分析了前缘剖面形状对三角翼、双三角翼涡运动影响的不同机理 :对三角翼,尖前缘可以形成组织最好的涡结构,但对于双三角翼,圆前缘生成的旋涡结构较靠近翼面,涡结构紧密,诱导能力较强,可以形成有利的涡涡干扰,使内翼涡通过剪切层向外翼涡输入涡量更加容易,合并涡变得更加稳定,推迟了涡破裂,而且由于涡较靠近翼面,因而可以产生较高的非线性涡升力,这同传统的认识是不一致的。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

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

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

机翼动态气动特性实验研究

介绍了动态特性风洞实验的技术方案、测试系统和测试软件。对前缘后掠６０°的三角翼在迎角０°～９０°范围内快速上仰时的升力进行了测量。结果表明：机翼快速上仰时,随上仰速率增大最大升力系数增高,即动态升力效应更明显;失速迎角也增大。此外,给出了简缩频率、平均迎角、振荡振幅以及俯仰轴位置等对振荡机翼产生的非定常瞬时力的影响