粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

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

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

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

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

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

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

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

用粘流—无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进，改进的半反－局部联立耦合方法。考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好，算例表明，翼型高升力状态计算必须包括尾流的作用，也应用对势流计算压力与实际压力之间     

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

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

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

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

球柱三维分离及涡流的数值模拟

对三维可压非定常Ｎ－Ｓ方程进行有限体积中心离散,并用Ｒｕｎｇｅ－Ｋｕｔａ显式格式进行数值求解。数值模拟了球柱分离及涡流的形态结构,给出了具有激波、边界层、大分离区、涡流及其间相互干扰的复杂流型     

Box法用于三维边界层计算的数值稳定性分析

本文对 Box法用于三维边界层计算的数值稳定性做了 Neumman分析,并借助于数值计算,证实了目前公认的一些稳定性问题和三维边界层计算中所发现的现象,得到了一些有参考价值的结论。      

用弱非线性理论研究二次稳定性问题

本文把研究三波共振的弱非线性理论应用于平板边界层和平面Poiseuille流三维扰动的二次稳定性问题。其数值结果与文献[1]、[2]用Floquet形式的二次稳定性理论得出的结论是一致的。     

边界层分离流边值问题的算法研究

本文给出了边界层方法边值问题的一种快速求解方法,该方法可用于二维与三维问题,它是在正方法边界层算法的基础之上增加了一个矩阵算法,该矩阵很容易逆,总的ＣＰＵ时间仅增加约２０％,使用该算法对反方法问题求解和进行数值稳定性分析都将十分方便。     

轴向前倾叶片流场的数值分析

对一种新式的风机叶片即轴向前倾叶片流场进行数值计算, 采用时间推进法计算无粘流场, 采用非正交曲线坐标系附面层积分方程求解附面层参数, 得到了一些有意义的结果。      

超音速绕凹角湍流分离流动的计算

本文用统一的Levy-Lees变换以及正算法与逆算法相结合,求解了超音速绕凹角湍流分离流动。 对附着流区用边界层正算法,压强分布用流过尖劈统一的高超音速与超音速公式,湍流模型取代数涡粘性模型;对凹角分离区用边界层逆算法,给定位移厚度δ~*分布,湍流模型取代数松弛模型;边界层计算采用Cebeci-Keller Box方法;计算成功地算得分离流场,较好地预估了分离点与重附点位置以及壁面压强分布与表面摩擦应力分布。     

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

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

利用流线法决定绕扁长椭球流动分离的位置

本文按照流线法把所测量的扁长椭球的压力分布用于决定涡层型三维流动分离,按边界层边缘处流线会聚的位置近似地决定分离线的位置。作为算例,计算了扁长椭球体在α=30°状态下的绕流,并在水洞中作了流场显示,数值结果和实验结果的对比表明,二者十分一致。分析表明,涡层型分离不仅与壁面极限流线的收敛有关,而且还与边界层边缘处流线会聚有关,前者的出现并不是引起涡层型流动分离的充要条件。