跨声速翼型绕流计算边界层反方法的研究

本文用有限差分方法,通过有粘/无粘迭代计算了二元翼型的跨声速绕流问题。在边界层粘性区域内考虑了层流、转捩及湍流流动,当边界层内出现分离时,使用边界层反方法,采用代数湍流模型。算例表明,对激波/边界层弱干扰和强干扰情况,该方法的结果与风洞实验结果吻合良好,对于求解边界层小分离流场是一种好的近似方法。     

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

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

跨声速翼型有粘／无粘强干扰计算

本文介绍一种二元跨声速激波-边界层强干扰的计算方法。边界层计算采用湍流边界层积分反方法,它借助Whitfield和Swafford提供的既适合附着流,也适合分离流的速度剖面表达式。跨声速无粘流用全速势方程模拟。通过边界上排溢速度来考虑粘性的影响,用有粘/无粘迭代得到粘性流解。本方法计算的结果与其它方法以及实验的结果进行了比较,证明该方法可以在工程上推广使用。     

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

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

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

平面跨音叶栅中粘流-无粘流的强相互干扰计算

为了设计高效率高负荷的跨音叶轮机械,气动计算中必须计入粘性的影响。当叶片表面流动不发生分离,或分离区很小时,流道内的流动可以分为叶片表面以及尾缘下游的粘性薄层——边界层和尾迹区以及除此以外的无粘流区。因而可以分别应用欧拉方程组及边界层方程组求解这两个区内的流动。无粘流-粘流的相互影响是通过无粘流计算及边界层计算的一系列迭代计算来实现的:由给定的气动参数及叶栅几何形状进行无粘计算,得到叶片及尾迹表面的流动参数分布,作为首次边界层及尾迹计算的输入数据,求得边界层及尾迹的排挤厚度δ_1,并以此对无粘流动的几何边界进行修正,作为下一次无粘计算的流场边界。这样的迭代重复进行,直到满足全场的总体收敛条件为止,其迭代框图如图所1示。     

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

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

计及边界层汇流效应的多段翼型失速特性计算

 本文给出计及边界层汇流效应的多段翼型失速特性的解法。用高阶奇点分布面元法（Panel Method）求位流解,然后进行各翼段粘性尾迹形状迭代,并解出正常边界层和汇流边界层特性。当翼段上有后缘分离时,还要确定分离尾迹的形状,用位移厚度当量源（汇）模拟粘性效应。进行粘/位流迭代直至收敛。上述方法在超过多段翼型失速迎角时仍然有效。计算结果与实验数据比较,符合良好。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

(高)超声速流动试验技术及研究进展

近年来,与高速飞行器相关的(高)超声速流动受到了极大的关注。这类流动所具有的非定常性、强梯度和可压缩性对试验方法和风洞设计技术提出了挑战。超声速纳米示踪平面激光散射(NPLS)技术是由作者所在团队研发的非接触光学测试技术。它能够以较高的空间分辨率来揭示超声速三维流场的一个瞬态剖面的时间解析的流动结构。介绍了NPLS技术以及基于NPLS开发的密度场测量、雷诺应力测量和气动光学波前测量等方法,并回顾了这些技术在超声速边界层、超声速混合层、超声速压缩拐角、激波/边界层相互作用和光学头罩绕流等流动中的应用,清晰地再现了边界层、混合层、激波等典型流场结构及其时空演化特性。另外,为了模拟和研究高空大气条件下边界层自然转捩和超声速混合层的转捩特性,介绍了高超声速静风洞、超-超混合层风洞的设计技术以及层流化喷管的设计方法。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

二维超／高超声速进气道流场数值模拟

对超／高超声速三级压缩进气道流场进行了数值模拟,来流马赫数为4、6,进气道内流动为层流状态,根据二维Navier-Stokes方程,采用二阶精度Roe格式进行离散。按照流场特点,合理地设计网格分布及调整不同黏性范围的熵修正,防止了壁面附近过大的数值耗散,使计算结果更加合理。在进气道模型的各级压缩折转角处,获得了清晰的激波结构,在进气道内部的各种波系的相交、反射和激波诱导的边界层分离等现象都得到合理的描述。计算得到的压力分布,在各级压缩斜板上同简单波理论结果十分接近。用本文方法计算了另一个二级压缩进气道,沿上、下壁面的压力分布与试验比较符合得较好。