NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

γ-Reθ转捩模型在二维低速问题中的应用

介绍了基于当地关联的γ-Reθ转捩模型.通过Spalart建议的环境源项法对自由流湍流度衰减的有效控制,实现了该转捩模型的实用化.利用中国空气动力研究与发展中心的TRIP软件平台上实现的该转捩模型,对S809翼型、NLR7301两段翼型和30P30N三段翼型3个低速翼型绕流问题进行了计算,将不同迎角下的计算结果与试验结...     

不可压缩实际粘性流体绕流孤立翼型的数值研究

本文提出了基于微分-积分方程组求解N-S方程的有限差分法求解不可压缩实际粘性流体绕孤立翼型流动。计算中采用了K-ε湍流模型结合考虑逆压梯度的壁面函数法,并导出任意非正交曲线坐标系下的压力修正方法和适用于大流通量的混合格式。分别计算了层流和湍流情况下带分离翼型绕流,与流动显示及实验数据比较表明,在预测分离点位置,翼型表面压力、平均速度分布、雷诺切应力uv方面,两者吻合,结果是满意的。     

基于非结构重叠网格的二维N—S方程求解与应用研究

本文发展了非结构重叠网格技术.通过构造任意单元类型网格重叠的网格间边界确定算法、宿主单元搜索算法以及网格间插值算法,把非结构重叠网格推广应用于求解NS方程,并结合湍流模型模拟二维湍流流场.N-S方程的求解采用有限体积法,通量计算采用Osher格式,湍流模型采用SA模型.为验证本文方法的正确性和有效性,模拟了RAE2822翼型和某三段翼型绕流湍流流场,结果表明采用非结构重叠网格模拟复杂外形湍流流场不但网格生成效率高,而且容易推广应用于具有相对运动的非定常问题,具有较高的工程应用价值.     

低速两段翼型的数值模拟

在多段翼型风洞实验中,很难观察翼型不同迎角下主翼、襟翼上的绕流及缝道流动。使用了O-H型混合的结构化网格,采用k-ε二方程湍流模型求解可压的N—S方程来模拟二元风洞中两段翼型的流场。计算结果与实验进行了比较,结果表明：计算与实验结果吻合良好,说明本方法可以较好地模拟两段翼型的绕流,并能很好地显示流场的变化情况。另外,从结果中可以看出：随着迎角的变化,主翼尾流和缝道间的流动对襟翼附面层分离有很大的影响作用。     

高超声速复杂外形热流的NND有限元数值计算方法研究

有限元法是一种加权积分的方法,这个特点和权函数的特性使得有限元法在计算诸如热流(与温度的一阶导数相关)等物理量时有其内在的优越性.本文研究了高超声速复杂外形热流的NND有限元数值计算方法,并对高超声速球钝锥外形和带翼飞行器外形的绕流进行了数值模拟,算例结果令人满意.     

非定常自由来流对二维翼型气动特性的影响研究

利用计算软件和实验测量方法,研究了非定常自由来流对静态二维翼型气动特性的影响,分析研究了来流速度脉动频率变化对气动特性产生的作用.结果表明,来流速度以短周期脉动时,升力系数随来流速度的减小而增加;来流速度以长周期脉动时,升力系数随来流速度的减小而减小.分析表明这与前缘分离涡在翼面上的传递过程有关.又利用二维翼型动态实验台,研究了非定常自由来流对做动态运动的二维翼型气动特性的影响.结果表明,来流风速的脉动使升力系数的迟滞包线进一步扩大,最大升力系数增加.     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

二维凹腔超声速流动的混合RANS/LES模拟

基于Menter的k-ω SST湍流模型构建了一种混合RANS/LES模拟方法,通过采用一个与到壁面距离相关的衔接函数将处理近壁区的SST模型过渡到处理主流区的Yoshizawa一方程亚格子模型.利用此方法对用于超燃冲压发动机的凹腔的二维超声速流动进行模拟,模拟结果再现了二维凹腔剪切层的拟序结构,计算得到的凹腔自激振荡频率、时均统计的压力分布以及压力剖面与实验结果吻合较好.     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

超声速湍流边界层中横向声速喷流的混合LES/RANS模拟

通过对超声速湍流边界层中横向声速喷流的计算,对一种混合大涡/雷诺平均Navier-Stokes(LES/RANS)模拟方法进行了测试,该方法采用一个依赖于到壁面的距离及当地湍流参数的混合函数结合两方程k-ω SST(shear stress transport)湍流模型和混合尺度亚格子模型来封闭湍流项.计算结果表明:混合模拟方法能够捕捉到喷流/湍流边界层相互干扰的非定常大尺度结构,且对分离区长度、壁面静压峰值和膨胀区静压分布的计算精度要高于RANS(Reynolds-averaged Navier-Stokes)方法.      

基于RANS/LES方法的超声速底部流场数值模拟

分别采用基于两方程k-ω剪切应力输运(SST)湍流模型的延迟DES(DDES)、更改的DDES(MDDES)和改进的DDES(IDDES)方法,并引入可压缩修正,结合三阶MUSCL-Roe和五阶WENO-Roe两种空间离散格式,针对超声速底部的复杂流动现象,开展了数值模拟研究。计算结果表明本文方法能够捕捉到超声速底部流动中丰富的湍流结构,通过分析计算结果对超声速底部的流动机理有了进一步的认识,为下一步的超声速底部流动减阻改进和雷诺平均NavierStokes/大涡模拟(RANS/LES)方法在非定常高可压缩性流动中的应用提供了参考。通过对比分析不同空间离散格式的计算结果研究了数值耗散对计算的影响,五阶WENO-Roe格式的计算结果与实验结果吻合良好;对不同RANS/LES混合方法的计算结果进行了对比分析,结果表明IDDES方法在近壁区的表现优于DDES和MDDES方法。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

多段翼型大迎角分离流动的DELAYED RANS/LES混合算法

多段翼型的大迎角绕流发生大范围附面层分离,具有明显的三维与非定常流动特性。RANS/LES混合算法继承了LES对流动分离区大尺度漩涡准确模拟的优点,避免了纯LES算法需求网格量巨大与亚格子模型壁面函数不成熟等问题,对分离流动的模拟效果优于RANS算法。以S-A湍流模型与Smagorinsky亚格子模型为基础,借鉴DDES的附面层延迟控制思想,构造了可用于对接网格、重叠网格的DELAYED RANS/LES混合算法。研究了GA（W）-1多段翼型的大迎角分离流动及其气动特性。     

战斗机垂尾脉动压力数值模拟

在亚跨超计算流体力学（CFD）软件平台（TRIP）上开发了基于RANS/LES混合思路的IDDES流动模拟技术,并通过NACA0021翼型60°大迎角分离流动与串列圆柱绕流模拟对RANS/LES混合方法的精确度进行了验证,针对某战斗机外形的垂尾脉动压力开展了数值模拟研究。战斗机来流马赫数为0.1,基于全机长度的雷诺数为2×10⁶,模型迎角为20°、30°和40°。分别通过脉动压力系数、脉动压力功率谱密度、空间流动结构以及侧向力响应曲线等对战斗机的垂尾脉动压力进行了分析。脉动压力模拟结果表明：当垂尾完全沉浸在边条翼脱体涡破碎后的宽频湍流脉动气流中时,垂尾翼梢位置的脉动压力会发生明显的增大。     

复杂湍流流动的混合RANS/LES方法研究

混合RANS/LES方法是近年来复杂湍流模拟的重要方法。简要回顾了混合RANS/LES模拟方法的发展历程,着重总结和分析了分离涡混合模型、类Menter SST混合模型以及湍流能量谱一致混合方法三类混合模型的构造方法和发展历程,对这三种方法的应用以及优缺点进行了简要评述,最后指出了混合RANS/LES方法发展应该兼顾的问题,为下一步的混合RANS/LES模拟提供了参考。     

一种滤波SST方法在翼型深失速模拟中的应用

 为了提高原始剪切应力输运(SST)湍流模型对于分离流动的求解精度,将大涡模拟(LES)中的滤波因子和SST方程相结合构造出一种滤波SST方法,利用湍流尺度对流场求解区域进行划分,近壁面附近的稳态流动由湍流模型控制,远壁面采用LES方法进行模拟。与传统混合RANS/LES方法相比,该方法的特点是:滤波因子的选取不再依赖于网格尺度,可以有效地降低网格诱导分离现象发生的概率。采用该方法对NACA0021翼型深失速特性进行了仿真研究,对比了非定常雷诺平均Navier-Stokes(URANS)方法和SST-DES方法,从仿真结果可以看出滤波SST方法有效地激活了分离区域的脉动,充分展现了分离的三维特性;同时算例求解结果证明该方法的精度高于URANS方法,与试验结果吻合较好,显示其具有一定的工程应用价值。     

超声速燃烧凹腔质量交换特性的混合RANS/LES模拟

对超声速冷流及燃烧流条件下超燃冲压发动机火焰稳定凹腔的质量交换特性进行了研究.采用混合RANS/LES(Reynolds-averaged Navier-Stokes/large-eddy simulation)方法对非定常流场进行数值模拟,系统研究了凹腔结构参数和横向喷流燃烧放热对凹腔质量交换特性的影响.首次采用"原子追踪法"对化学反应流条件下的凹腔质量交换过程进行了刻画.结果表明:在冷流条件和燃烧流条件下,凹腔驻留时间均随长深比的增大而增加;凹腔后缘倾角对驻留时间的影响在冷流条件下比较显著,而在燃烧流条件下明显减弱;相同结构的凹腔,在燃烧流条件下的驻留时间相对于冷流条件明显减小.      

RANS/LES在超声速突起物绕流中的应用研究

安装在超声速/高超声速飞行器表面的突起物如机翼、控制舵等通常会导致复杂的激波/边界层干扰,对突起物的局部气动特性甚至飞行器整体的气动特性产生较大的扰动.在采用计算流体力学(CFD)数值模拟此类问题时,传统的求解雷诺平均Navier-Stokes(RANS)方程方法由于不能准确预测湍流脉动流场并且精度有限,在应用上受到一定的限制.本文在研究B-L (Baldwin-Lomax)内层模型和Smagorinsky亚格子模型优缺点的基础上,提出了一种新型的RANS/LES(Large Eddy Simulation)混合模型,并进行了算例验证,证实了该方法的可行性.在此基础上,对火箭表面突起物的干扰流场进行了数值模拟研究,细致地刻画了突起物附近的激波/边界层干扰、剪切层失稳和底部分离涡形成的非定常过程,获得了突起物及火箭表面上的压力脉动历程并进行了频谱分析.研究发现,相对于突起物底部的非定常分离流动,突起物前缘的激波和边界层相互干扰的非定常过程是突起物周围压力脉动的主导因素,这种高频的压力脉动可能对火箭内设备的正常工作产生不利的影响.