二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。     

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

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

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

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

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

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

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

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

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

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

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

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

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

脊状结构对翼型边界层分离特性影响大涡模拟研究

流动分离、湍流再附现象对翼型空气动力性能影响较大,为了改善翼型气动性能,研究了脊状结构对翼型边界层分离及尾迹速度的影响。采用数值计算方法,分别将脊状结构布置在NACA0018翼型的顺压梯度区(前段)和逆压梯度区(后段)内,分析了脊状结构对翼型边界层速度分布和尾迹速度分布的影响。研究结果表明:在迎角6°,来流速度为24m/s和12m/s下时,脊状结构前段布置时,翼型边界层分离点略有提前并且分离区域提前结束,其分离程度相对微弱。相比之下,脊状结构后段布置时,在推迟了边界层分离点的同时提前结束了边界层的分离区域,其尾迹速度亏损更小,亏损区域的面积也更小,边界层控制效果更为明显。两种脊状结构均可以有效的控制边界层的分离,缩小边界层分离区域的范围,减小尾迹速度损失。在翼型表面合理的布置脊状结构为翼型流动控制提供了新的思路。     

中介机匣附面层数值计算与试验分析

对位于低压风扇出口和高压压气机进口间的中介机匣壁面附面层的分布作了数值计算和试验研究。采用反积分和差分方法分别求解有分离的紊流附面层和层流附面层方程，对层流分离到紊流再附着，层流到紊流的转捩及紊流分离进行了判别，并测量了不同进口马赫数条件下的中介机匣壁面的静压和附面层厚度的分布，对比情况说明，计算结果和试验数据基本吻合。     

翼型尾缘噪声源空间分布与辐射特性关系

详细研究了翼型湍流边界层尾缘宽频噪声源空间分布与辐射特性的关系.采用基于雷诺平均流场的翼型尾缘宽频预测方法研究了NACA0012翼型湍流边界层尾缘宽频噪声在4种不同工况下的噪声源空间分布与辐射特性.首先计算了NACA0012翼型湍流边界层尾缘噪声源在不同频率下的空间分布.计算结果发现:边界层中湍流是翼型湍流边界层尾缘噪声声源.随着频率的增加,噪声源强度和噪声源空间尺寸都是先增加后减小,噪声源位置不断靠近翼型尾缘.同时也计算了边界层内不同位置处的噪声源对远场噪声的辐射特性,结果表明:边界层内层区域,其噪声频谱能量集中在高频;边界层外部区域,其噪声能量集中在中低频;攻角增大或者来流速度减小,噪声能量向低频转移.      

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

跨声速翼型上激波／边界层干扰的自适应控制计算

在激波区使用自适应壁对跨声速翼型上的激波/边界层干扰进行控制,可改变机翼的气动性能。这种被动控制可通过在翼型的激波区开一凹腔,其上覆盖一弹性橡胶膜柔壁来实现。本文给出用N-S方程数值模拟这一自适应控制翼型的跨声速粘性绕流,提出了一个适用于本特殊情况(物面边界局部地区在求解过程中有变化)的处理办法。并探讨了自适应柔壁对当代跨声速翼型绕流的影响。     

三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

后向台阶湍流边界层分离、再附及发展

应用激光测速仪对后向台阶湍流边界层分离、再附及发展进行了测量,得到时均流速、流向湍流度、平坦系数和偏斜系数等沿程分布。依据时均流速分布特性,得到分离区边界及再附点位置ｘ＿ｒ＝５．２ｈ＿ｓ（ｈ＿ｓ为台阶高度）,再附后湍流边界层恢复速度型相似的位置ｘ＿Ｔ＝１３．５ｈ＿ｓ。在详细分析湍流边界层区域时均及紊动特性的基础上,与光滑平板湍流边界层进行了比较。结果表明,除流速分布对数公式中积分常数较大外,无其它差别。此外,湍流边界层区存在尾流律,尾流参数π＝０．２９。     

现代自然层流翼型的设计方法

本文介绍一种设计跨声速自然层流翼型的计算流体力学(CFD)方法。本方法采用“正反迭代、余量修正”设计原理,通过将跨声速翼型设计软件NPU-TD2D中的反设计程序进行改进,并与含有层、湍流混合边界层修正的跨声速层流翼型计算程序DLRBGKWALZ耦合,实现了在跨声速粘性流动条件下直接设计层流翼型。亚、超临界的设计实例和风洞验证表明,本方法可以在几个设计迭代内设计出压力分布、转捩位置及气动参数均准确收敛于设计目标的新翼型,是一种设计现代自然层流翼型的有效而实用的CFD方法。     

扩散界面浸入边界法结合壁面模型在湍流模拟中的应用

提出了一种模拟高雷诺数湍流的扩散界面浸入边界法（IBM）。该方法采用壁面模型来减少壁面附近的网格量。为了实施壁面模型和边界条件，在物面外部设置了2个辅助层（一系列拉格朗日点）：外侧的参考层用于实施壁面模型来得到壁面剪应力，内侧的强制层用于实施离壁的边界条件。在实施壁面模型时，将动量方程沿物面法向积分，从而把物面切向的速度修正量与由壁面模型得到的壁面剪应力联系起来，而速度的法向分量则按二次曲线分布来近似重构。在实施边界条件时，为了严格满足无穿透条件，应用了基于隐式速度修正的IBM。最后，利用绕平板湍流和绕NACA0012翼型湍流来对方法的可行性进行了验证。     

计算三维湍流边界层的积分反方法

本文通过分析、计算比较了现有计算三维湍流边界层反方法几种主要版本。内容包括小横向流和大横向流两部分。在此基础上提出一较好的反方法计算版本。所得结论对于如何选取一合适方法有参考意义。