计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

关于“跨音速流的边界积分法”一文的某些问题

 最近杨岞生教授和J.C.Wu教授提出了一种跨音速流的边界积分法。该方法以脉冲压强为基础建立非线性可压缩位流的边界积分表示式,利用边界元法给出了二维、三维机翼跨音速位流的计算方法。我们认为文献[1]在推导该方法的基本方程（即主管脉冲压强的方程）和求此方程的基本解这两个方面存在一些问题。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

机翼跨声速绕流的压强极小积分有限元算法及应用

本文从全速位方程出发。利用压强极小积分,给出了机翼跨声速绕流的有限元算法,并编制出相应的FORTRAN语言程序。文中给出了适于有限元法的网格生成技术,采用了线松弛迭代法和逐步加密法进行数值求解,提供了M6机翼和一个三角翼的计算结果。     

二元全速势跨音速方程的数值解

本报告提供了二元全速势跨音速方程的一种数值解法。本方法的特点是放弃了繁杂的曲线座标系,而在简洁直观的直角座标系中,利用Jameson的旋转差分格式,并在物面边界上采用速势内嵌场,建立全流场的差分方程。二元翼型的样例计算表明本方法和其他全速势方程的数值解法具有同等的精度。对物面外形作一次附面层位移厚度修正后的算例表明计算结果和试验结果符合较好。 本方法比起采用曲线座标系的方法,更易推广到三元流。     

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

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

跨音速翼型和机翼的反设计计算方法

给出一种跨音速翼型和机翼的反设计计算方法。对所应用的积分方程反方法引入人工粘性项;采用Riegels因子法消除前缘奇性;对强激波问题采用光滑-松弛过程;并将方程中的系数积分成解析形式;对二维翼型反设计计算还提出了一种封闭形式的正则化条件。算例结果表明,该方法对跨音速翼型和机翼设计是一种有效的工具。     

绕升力机身跨音速流动的有限元素法数值分析

 <正> 本文应用全速位方程最小压强积分的有限元素法解绕升力机身的跨音速流动。机身头部可以是尖的或带有进气道的。用头部伸出无穷长圆柱来模拟头部进气。尖头机身在元素形状处理上具有一定的复杂性,从数值实验角度考虑我们分别按尖头和头部伸出无穷长细圆柱来近似尖头头部两种方法对尖头机身进行计算。采用人工密度法捕获激波,但由于绕机身流线形状比较复杂,直接采用机翼问题中的人工粘性公式,解往往不收敛。我们推导出一种较精确的,适用于复杂流动情况的人工粘性计算公式,为了加速收敛,应用网格逐次加密技术,可使收敛速度提高2～3倍。本文还给出绕升力机身跨音速流动的远场速位解析式。     

三维薄翼非对称绕流的跨音速积分方程及其有限元解法

与有限差分法相比,积分方程的有限元解法具有计算工作量小,对计算机容量要求低而又能达到一定的计算精度等优点,但国外的研究尚限于二维翼型。最近,文[4]发展了积分方程的有限元解法并成功地应用于求解三维薄翼跨音速对称绕流问题。在此基础上,本文进一步发展该方法并用于求解三维薄翼跨音速有升力的非对称绕流问题。     

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

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

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

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

计算机翼跨音速绕流的加罚Galerkin有限元法

1.引言 全速位Galerkin有限元法是跨音速有限元法中重要方法之一。处理超临界流时,为了自动捕捉到激波,通常采用上风技术,但因此导致了质量守恒条件得不到完全满足。文献[1]从Bateman变分原理出发,应用加罚函数方法,找出了强制满足质量守恒条件在Galerkin有限元方程中所应加入的附加项。本文将文献[1]的方法应用到超临界跨音速绕机翼流动,克服了一般Galerkin有限元法中由于采用上风技术而出现的激波不够陡和位     

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

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

关于“跨音速流的边界积分法”一文的商榷

 近年来,杨岞生、吴镇远提出了一种新的解决非线性可压缩流的边界元法,他们在引进压力冲量之后将一个非线性问题转变为线性问题,然后利用边界积分表示式给出了二维、三维机翼跨音速绕流的计算方法。若干数值算例说明该方法的结果与其它数值方法和实验结果相比符合得较好。由于该方法能够将一个非线性可压缩位势流转换为对压力冲量的线性方程,因此理所当然地引起了国内外学者的关注。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

复合材料超临界机翼跨音速颤振实验与数值计算研究

颤振课题是飞机设计过程中常常遇到的一个关键技术问题。以支线客机ARJ21超临界机翼颤振特性研究为背景,在俄罗斯TsAGI的T-106风洞中完成了该复合材料机翼跨音速颤振实验,基于N-S方程和无限插值方法（TFI）生成三维贴体运动网格对超临界机翼跨音速颤振进行了并行计算。结果表明：复合材料的超临界机翼在跨音速区域具有跨音速颤振＂凹坑＂现象;与风洞实验结果相比,有较好的一致性,为使用超临界机翼的运输类飞机跨音速颤振特性预计提供了一定的参考。     

大扰动跨音速流积分方程的导出及亚临界流的数值解法

本文对一方向上大扰动的平面跨音速势流微分方程提出了新的近似简化方案,从数学上完成了相应积分方程的推导;其中关于边界及激波条件作了小扰动近似处理;最后提供了关于亚临界流的算例。     

计算三维振动机翼上跨-亚音速压力分布的局部线化面元法

本文以小扰动位势流理论为基础,采用局部线化假设,建立了三维振动机翼的跨音速下洗积分方程;给出了跨-亚音速流核函数的计算方法;用推广的偶极子栅格法计算了矩形谐振动机翼上的非定常压力分布。结果表明,本方法可以较准确地计算非线性的跨-亚音速非定常流动,关于跨-超音速流动的数值解法,尚需作进一步的研究。