二维跨音速流动计算的显式多步有限面积方法

本文按Jameson等的思路,采用新的人工粘性项系数,建立显式多步有限面积方法,给出绕NACA0012翼型和RAE2822翼型跨音速流动的计算结果。     

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

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

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

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

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

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

基于CST参数化的翼型优化遗传算法研究

用外形建模的CST参数化方法,构建翼型几何可直接利用其后缘角、前缘半径等几何特征,其控制参数更能反映翼型特有的气动敏感性,有助于遗传算法搜索寻优。基于这一参数化方法,结合遗传算法,构造了用于翼型优化的设计方法。算法中,CST控制参数作为设计变量,采用二进制编码,并通过引入精英策略,提高了遗传算法的收敛性能。算法适应度评估涉及的流场求解则采用了基于Jameson有限体积法的Euler方程解算程序。先以NACA0012翼型为例,以其某一已知的表面压力分布为目标,进行了遗传算法的重构运算,给出了重构的翼型几何外形,验证了方法。在此基础上,进行了带约束的跨音速翼型优化设计,给出了升力系数极大化和阻力系数极小化等设计算例,展示出翼型优化设计的效果。     

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

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

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

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

叶栅跨音速流场的Euler解

应用Jameson的有限体积法求解二维Euler方程，模拟叶栅跨音速流场，并运用当地时间步长和隐式残差平均技术加快收敛速度，计算结果与理论解和试验数据吻合良好，本文比较了采用两种不同的通量计算公式的计算结果。     

在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

WENO格式在隐态问题中的应用

高阶、高分辨率的WENO格式大大加强了对于复杂流场结构的分辨能力和计算精度，但计算时间较TVD格式有了明显的增加。近似隐式分裂的LU方法能够加快求解稳态问题的收敛速度。本文结合WENO格式和LU方法实现了求解稳态问题的有效方法，它使得求解稳态问题在收敛速度和计算精度上都有较大的提高。通过算例的数值分析可以看出，上述方法不但收敛速度有了明显的提高，而且对于复杂流场现象仍然具有良好的分辨能力。     

喷管跨声速流场计算

本文用时间相关法求解定常跨声速喷管流场。对守恒型和非守恒型方程采用全场统一时间步长、当地时间步长和修正超阻尼法进行了计算。计算表明:采用当地时间步长可大大提高收敛速度,采用全场统一时间步长时,数值解非守恒型方程收敛较快;采用修正超阻尼法也能提高收敛速度,但与阻尼系数和时间间隔有关。     

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

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

跨音速叶栅粘流计算的多级广义有限元法

将 Taylor-Galerkin广义有限元法和多级有限元的思想结合起来 ,并在人工粘性的处理上作了改进 ,构成了收敛速度和稳定性均较好的多级广义有限元法。利用该方法 ,并基于 N-S方程研究了二维跨音速叶栅流动。计算结果与实验结果符合较好。通过计算表明 ,该方法计算稳定、收敛较快 ,是叶轮机械内部跨音速流动强有力的计算手段     

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

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

跨声速叶栅粘流计算的多级Taylor－Galerkin有限元法

本文对Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限无法在收敛速度和稳定性两方面作了改进，与多级有限元法相结合，构成了多级Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法。基于Ｎａｖｉｅｒ－Ｓｔｏ－ｋｅｓ方程，研究了平面跨声速叶栅流动，并与实验结果作了比较。计算结果表明新方法具有较好的稳定性及较快的收敛速度。     

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

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

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

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

亚音速升力面气动敏感性导数计算

具有任意曲线前缘的亚音速升力面的气动敏感性导数由核函数法给出。用自适应积分法计算弦向积分,用Multhopp法结合抽去奇点,计算Mangler积分主值。将积分核展成Chebyshef多项式的渐近展开式以保证结果的收敛性。最后将广义力系数及其敏感性导数表示成简单形式,对椭圆、矩形和后掠机翼作了计算,所得结果在升力面理论精度范围内与直接由核函数法得到的结果一致;而且所得到的偏导数可在飞机设计中分析综合用于多学科优化。     

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

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