均匀湍流场下矩形断面流场及气动参数数值模拟

采用大涡模拟方法数值模拟了两种矩形断面(B/H=1,B/H=2)在不同湍流来流下的流场及气动参数,采用谐波合成方法生成与目标谱一致的脉动入口风速.数值计算结果表明来流湍流对断面的流场及气动参数有较明显的影响,通过与均匀来流的计算结果作对比可以看出,湍流导致断面的阻力系数减小明显,其主要原因为背压区负压系数的增加,与实验的结果较一致;从流场分布看湍流主要影响断面尾流涡的形成,这种影响会影响尾流涡的形成.     

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

颗粒两相流自由剪切层大涡模拟

采用大涡模拟(LES)方法数值仿真了可压缩流自由剪切层内颗粒两相流流场结构.气相控制方程采用二维盒式滤波法滤除小于网格尺度的小涡以得到LES控制方程,而流场中小涡对大涡的影响采用K-ωSST湍流模型模拟.在Lagrange坐标系下用颗粒轨道模型仿真颗粒运动,用双线性插值方法实现Euler坐标系到肌Lagrange坐标系的转换,并考虑了颗粒速度和温度的改变对气相流场的反作用.分析了仿真结果中气相自由剪切层的涡结构特性,不同St数颗粒在流场中的运动特性,颗粒对流场的反作用效果,得到了合理的结论.     

超声速燃烧的湍流流场数值模拟

选择了6个有代表性的湍流模型,对超声速燃烧湍流流场开展了数值模拟.采用的模型包括5个线性涡粘性模型和1个非线性模型,在非线性模型中同时引入了可压缩修正.数值结果表明,本文发展的数值方法可以成功地应用于求解可压缩、含化学反应的多组分N-S方程.同时引入可压缩性修正的二阶涡粘性模型给出的结果较好.     

平板湍流转捩过程中Λ涡和环状涡的周围流场研究

采用高精度湍流直接数值模拟技术,对平板流动转捩过程中的涡系结构,尤其是Λ涡和环状涡附近的流场进行了详细的研究.高分辨率的数值模拟给出了实验手段无法得到的全流场信息,以及雷诺平均纳维斯托克斯方程方法(RANS)无法给出的精细涡结构信息.数值模拟在精确捕捉Λ涡和环状涡等转捩过程中典型涡系结构、上喷和下扫等典型流动现象的同时,发现了一些与传统认识所不同的现象:1)Λ涡并不是当地涡量集中的地方;2)环状涡的附近存在一个速度显著大于远场来流速度的区域.这些现象对传统湍流转捩理论提出了挑战,并引发人们对湍流生成与维持机理的进一步研究.     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

考虑湍流影响的前缘分离涡运动与破碎特性的研究

本文以Ｎ－Ｓ方程和Ｈａｌｌ涡核模型假设为基础，导出了描述湍流的涡核运动方程。利用差分计算方法，对三角翼前缘分离涡运动及其破碎特性进行了数值求解，分析了涡核流场的结构和各物理量的变化特性，反映了三角翼前缘分离涡运动及其破碎特性的实质，说明了湍流涡核方程能更有地模拟三角翼前缘分离涡运动和准确地确定涡核破碎位置，并得出了攻钐雷诺数对破碎位置的影响，计算结果与实验数据十分吻合。     

激波诱导轴对称气动矢量喷管流场数值模拟

在落压比3～10,次流相对流量比2.5%～20%工况下,采用RNGk-ε湍流模型对扩张段开缝的激波诱导轴对称气动矢量喷管试验件进行了数值模拟.结果表明:壁面静压分布计算值和试验数据相对误差不大于10.1%.次流的注入使得气动矢量喷管内流流动非常复杂,流场结构的主要特征是在扩张段有一对旋向相反的主分离涡与射流角涡和一个位于次流与出口截面之间较大的回流区.流场结构随着落压比和次流相对流量比的变化而改变.      

细长尖头旋成体大迎角非对称涡系结构

本文采用数值计算方法,对细长尖拱旋成体大迎角背风侧非对称涡系结构以及与沿轴向交变的侧向力分布的关系进行了研究.求解的是N-S方程,采用Jameson中心格式,湍流采用修正的B-L模型.通过数值模拟揭示了在头部开始产生非对称涡,形成二次涡的机理以及在背风侧形成集中涡的过程.从而说明了采用集中涡来模拟这类流场的正确性,表明了二次涡对多个脱体集中涡形成的贡献.     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.