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

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

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

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

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

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

三维跨声速流粘性—无粘相互作用问题的分析和计算

本文对跨声速粘性-无粘相互作用问题进行了分析,论证了如何选择湍流边界层反模式的输入量。在FL027多重网格基础上作了局部加密。计算例子表明,本文建议的半逆方法适用于弱激波情形,可用于实际设计之中。     

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

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

基于欧拉方程的高超声速飞行器的壁面流线生成计算

提出了一种新的利用表面流函数法计算高超声速飞行器表面流线分布的方法.首先提出了表面流函数的概念,并通过理论推导,得到了表面流函数与表面流线的关系;然后运用结构化网格求解三维Eu-ler方程,计算得到高超声速钝头体的边界层外缘外部无粘流场气流参数;最后利用无粘流场气流参数和表面流函数的方法计算了高超声速飞行器的精确表面流线分布.结果表明,在无攻角和攻角小于20°的情况下均可以得到较好的壁面流线分布.高精度的表面流线的得到,为利用边界层内粘性主导区域的积分方程等方法进一步精确预测高超声速飞行器表面的气动加热奠定了基础.      

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

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

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

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

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

计算带有层流分离的粘性/无粘流相互作用的一个方法

基于粘性/无粘流相互作用模型和边界层积分方程,推导出计算二维和三维不可压层流粘性/无粘流相互作用的方程组。该方程组可以用于带有中等程度的分离和再附情况下。文中给出了计算方法和算例。所得的二维计算结果与其他计算和实验结果相符。本法具有简便迅速的优点。     

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

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

一种跨音速叶栅流场分析系统

本文建立了一种工程实用的跨音速叶栅流场分析系统,可进行跨音速叶栅无粘绕流和无粘流-边界层迭代计算。叶栅无粘绕流的计算以Denton拟流向“相对差分”方法为基础,y向采用松驰线性插值算法,所形成的解法能保证以较少的网点获得较准确的结果。边界层的计算采用积分方法,并考虑了紊流边界层的分离问题,文中还介绍了无粘流-边界层的耦合方式。应用本文方法对国内外多种压气机和涡轮叶栅进行了计算分析,所得结果与实验的吻合程度较好,因此,所提出的叶栅流场分析系统,可望为工程设计上实用的计算手段。     

跨音速平面叶栅全位势无粘-粘性迭代计算

本文发展了一种全位势与附面层、尾迹流正逆模式积分法迭代计算程序,给出了计算结果与国内外四个不同的叶栅风洞试验结果的对比分析。结果表明,对比是较为满意的,本文在网格生成,无粘跨音流求解和无粘流与附面层、尾迹流迭代等方面所作的发展和改进适用,使整个程序系统具有工程实用价值。      

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

跨音速弯曲壁面正激波与湍流边界层干涉的计算结果

本文介绍弯曲壁面上湍流边界层与正激波干涉的某些计算结果。壁面曲率对激波前后的压力比有重要影响。在干涉区中,边界层位移、动量和能量厚度都有显著的增加。这种增加主要决定于激波前后的压力比。(TBTF)     

动态网格快速生成及组合体跨音速非定常气动力计算

采用代数方法、局部变形技术及快速动态弹性变形技术生成三维翼/体组合的贴体动态网格。用时间精确近似因式分解差分方法求解守恒型非定常全位势方程。考虑附面层影响，通过粘/位迭代得到跨音速翼/体的非定常气动力。所得结果与实验数据相比相当吻合。     

计算三维粘性／无粘性干扰流动的准联立分区算法

 应用摄动理论导出了粘性／无粘性流动的匹配条件，与边界层方程联立形成能计及粘性／无粘性干扰效应的模型方程。同时引用等效吮吸理论对该模型方程与全位势方程进行交替迭代求解，建立了计算高雷诺数流动的一种快速收敛的粘性／无粘性分区算法。计算结果表明该算法适用于包括小尺度分离的流动；与传统的半逆解法相比，计算效率明显提高     

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。