翼型跨音速绕流的数值计算

本文选用精确速势方程作为翼型跨音速无粘绕流的数学模型。在变换过的直角座标中,用有限差分(亚音速区用中心差分,超音速区用旋转差分)对精确速势方程离散化。差分方程形成的代数方程组用列松弛迭代法求解。 为了计入粘性效应,利用精确速势方程和附面层动量积分方程进行联合迭代求解。这一点对超临界翼型的计算显得特别重要。 算例与其它数值解及实验作了比较,对于一般翼型和超临界翼型,本文结果是良好的。     

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

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

关于跨音速零升力翼型绕流的松弛和两种自修正风洞的收敛性

本文采用横向小扰动而纵向大扰动速势方程,计算了跨音速零升力翼型的绕流。在线松弛的数值实验中,φ_γ的差分式用简单迭代和φ_(xx)的差分式用改进迭代时,稳定性较好。此结论与文献的线化理论分析相符。 本文用混合差分法数值模拟,证明了基于两个控制面上的静压和基于一个控制面及翼面上的静压的跨音速零升力翼型自修正风洞的收敛性。对前一种方案,NACA0012翼型,M∞=0.9,RAE104翼型,M∞=0.8,对后一种方案,NACA0012翼型,M∞=0.72,0.8,在迎角为零和风洞高度与翼弦之比为3时,均能收敛到无洞壁干扰的自由流。     

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

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

绕钝前缘翼型定常跨音速小扰动势流的一种计算方法

提高钝前缘翼型的跨音速压力分布计算方法的精度与效率对翼型设计十分重要。国内对跨音速小扰动势流方法进行了研究,但对M_∞>0.8的超临界情况尚未计算。国外文献[3]指出当M_∞>0.9时完全速势方程的一种方法失败了。 本文试图改进经典的小扰动势流方程,探索比较稳定的迭代方法和超松弛方法,以克服超临界流计算常常不易收敛的困难,使小扰动势流计算的应用范围扩大到更高的M     

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

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

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

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

翼型跨音速小扰动势流-边界层计算

翼型绕流的外部无粘流,采用跨音速小扰动速势方程和改进小扰动速势方程对三种差分格式进行对比计算:(1)Murman-Cole非守恒式;(2)守恒式;(3)Engquist-Osher式。在相同的初值条件下,它们的计算稳定性、收敛性和收敛解基本相同,但激波位置(1)在前,(2)在后,(3)在中间。对人工时间阻尼项εφ_(xt)的作用进行了线化分析,所得结论与数值实验结果相符。二级精度格式计算表明:采用这种方案能在几乎不增加计算机时的前提下,大大提高差分计算精度。本文最后采用边界层积分方程法,计算翼型边界层,并与势流进行迭代计算。在弱激波条件下,迭代收敛解与实验值接近。     

基于N-S方程的翼型实验侧壁边界层控制

使用三维可压、雷诺平均N-S方程数值计算了NACA0012翼型在二元翼型风洞中不同迎角和马赫数状态下的压力分布,改进了划分计算网格的Hilgenstock方法,发展了一种利用N S方程计算结果,来确定翼型实验中为消除侧壁边界层干扰的抽气压力的方法。通过比较利用本文介绍的方法计算的NACA0012翼型实验侧壁抽气压力和相同状态下,在西北工业大学0.3m×0.1m二维跨音速风洞中经油流显示证明为合适的抽气压力,结果表明该方法是可行的。     

适用于任意壁的二维亚跨音速洞壁干扰修正法

本文使用风洞壁附近控制面上的实测静压分布及模型所受的气动力来进行二维亚跨音速风洞试验的洞壁干扰修正。计算时不涉及洞壁的透气特性,可用于透气壁或实壁。根据试验马赫数范围的不同,采用不同的速势方程进行计算。提供了两种适应不同需要的方法:一种是适用于洞壁附近为亚音速流动情况的快速算法。另一种是可用于低、亚、跨音速风洞试验,并能判断试验数据“可修正”程度的有限差分解法。给出了计算实例,并比较了不同的速势方程对计算结果的影响。     

绕钝前缘机翼跨音速势流的松弛计算及其稳定性与收敛性讨论

本文在机翼钝前缘处用精确速势方程和精确的边界条件,其他地方用纵向大扰动而横向小扰动的速势方程和相应的边界条件,联立求解。数值算例1为矩形机翼,展弦比λ=12,翼剖面为NACA0012,自由流的马赫数M_∞=0.63,迎角α=2°,翼根剖面压力分布的计算结果与二元亚音速精确数值解（Sells,1968）接近。算例2为NACA RM A51G31实验的机翼,垂直于1/4弦线的翼剖面为NACA64A010,其后退角χ_1/4=45°,λ=3,根梢比η=2,M_∞=0.4,0.8,0.9,α=2°。计算与实验接近。 本文建立跨音速定常小扰动速势差分方程的线松弛改进迭代在局部线化假设下的稳定性条件和松弛解收敛到原来的微分方程解的条件。这些条件大多数与数值实验相符。     

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

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

任意翼型有小分离气泡跨音速绕流的迭代算法

 本文给出一种带小分离气泡的任意翼型粘性跨音速绕流的计算方法,采用有粘-无粘干扰迭代的概念。无粘流的全速势方程用AF差分格式在保角变换法生成的计算网格中求解,粘流附面层方程用C-S盒式法求解,用逆算法消除分离点处的奇性。本文对Ma_∞=0.8,Re_∞=2×10~6,迎角α=3.5°和4°的NACA64A010翼型粘性绕流进行了计算,结果与实验相比较,吻合良好。     

三角翼过失速非定常洞壁干扰修正

分别以固壁条件和洞壁附近实测的压力分布,模拟各类实壁和透气壁试验段的洞壁边界条件,利用N-S方程数值求解模型在风洞中的绕流场,得出洞壁干扰对跨音速模型绕流的影响和气动力修正结果。初步研究结果表明,该方法能较有效地模拟模型在跨音速风洞中的绕流场,GBM04A模型在0.6m风洞中的试验结果经洞壁干扰修正后与无干扰参考结果吻合较好     

跨音速翼型分离绕流的粘性/无粘干扰迭代计算

本文用有粘/无粘干扰迭代的概念计算了跨音速任意翼型的绕流问题。位流的速位方程用AF2格式求解,而边界层微分方程用C-S盒式法求解,逆算法的引用可以克服边界层方程在分离点处的奇性问题,对分离区湍流代数模型的修正可以得到与实验更吻合的结果。计算结果表明有粘/无粘干扰迭代概念在小分离泡的情况下也是适用的。     

跨音速飞机模型试验洞壁干扰数值模拟的初步研究

分别以固壁条件和洞壁附近实测的压力分布,模拟各类实壁和透气壁试验段的洞壁边界条件,利用N-S方程数值求解模型在风洞中的绕流场,得出洞壁干扰对跨音速模型绕流的影响和气动力修正结果。初步研究结果表明,该方法能较有效地模拟模型在跨音速风洞中的绕流场,GBM04A模型在0.6m风洞中的试验结果经洞壁干扰修正后与无干扰参考结果吻合较好。     

Euler方程与边界层积分方程耦合求解跨音速翼型绕流

 应用Euler方程求解跨音速翼型特性时考虑了粘性影响,粘性影响是通过边界层动量和能量积分方程求解的,即粘流/无粘流迭代方法。其中Euler方程采用LU-ADI方法求解;边界层方程均由正解法过渡到反解法,以解决强激波干扰区出现小分离泡的计算问题。计算中使用了贴体C网格,通过一定变换使其保持基本正交。计算结果表明,压力分布、摩阻系数分布与实验结果符合较好。     

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

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

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

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

一种改进的计算压力分布的格林函数方法

对计算压力分布的格林函数方法作了一定的改进,虽然格林函数方法能给出比较准确的速势分布计算,但在压力分布的计算方面却比较欠缺,只给出了一个根据离散点上的速势值进行有限差分计算的简单的公式,采用这个公式无法算出物面准确的压力分布。鉴于此,通过对速势积分方程离散化后所得到的扰动速势的解析表达式进行微分运算,得到一组计算扰动速势导数的解析公式,并由此计算了物面上各离散点上的压力系数,其结果相当准确,与Navier-Stokes方程算法及实验结果相比,均吻合较好。另外,还对库塔条件在数值求解速势积分方程中的应用作了一定的改进,使计算结果更加准确。