基于通流方法的某涡喷发动机整机数值仿真

采用轴对称方法对带有黏性力的三维Euler方程组进行降维，利用时间推进方法求解，得到适用于航空发动机整机计算的准三维数值仿真软件，并对某涡喷发动机整机进行设计点和非设计点数值模拟。首先，对地面静止状态节流特性进行研究，将计算结果与试验数据对比可知：推力的最大相对误差为-5.1%，单位燃油消耗率的最大相对误差为+4.8%，相对转速为95%时，单位燃油消耗率最低；其次，获取了飞行马赫数为0.7工况下的高度特性以及飞行高度为3 km工况下的速度特性，将计算结果与设计参数对比可知：对于高度特性，推力的最大相对误差为-4.61%，单位燃油消耗率的最大相对误差为+5%，对于速度特性，推力的最大相对误差为-5.83%，单位燃油消耗率的最大相对误差为+5.92%；再次，分别对压气机/涡轮进行部件模拟，预测了发动机的共同工作线；最后，对发动机设计工况下的流场以及气动参数的展向分布进行了分析。     

超声速来流与侧向射流作用下的三维流场结构

研究超声速流动与侧向射流作用下的三维流场结构对飞行器表面流体分离控制以及超燃发动机燃料混合等具有重要作用。基于大涡模拟(LES)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=4,平板上侧向射流的流场结构特性进行了数值模拟。模拟结果表明,流场内的主要波系由弓形激波、分离诱导的激波以及桶形激波组成,这与相关实验结果相符。同时,计算结果还清晰显示了射流与主流相互作用下,射流上游的分离区流动呈现高度瞬态特性,分离区内旋涡结构不断变化,数目由2个到4个再到6个不断发展,平板表面流线复杂。射流尾涡摆动,涡管上扬,呈现非定常特性。     

DNS/LES方法在剪切湍流模拟中的应用

提出采用DNS/LES混合方法进行自由剪切类混合流动湍流研究,该方法兼顾了直接数值模拟(DNS)方法在转捩线性范围内的高精确度模拟和大涡模拟(LES)方法在转捩非线性及全湍流区的高效模拟两方面优点,通过对低对流Mach数流场计算表明该方法可以较好地模拟流场的转捩过程.同时,文中还对比实验分析了湍流场对光学传输效应的影响,给出了与实验一致的结论.     

基于大涡模拟的圆转矩喷管尾喷流强化掺混机理研究

为揭示2元圆转矩喷管尾喷流强化掺混的内在机制,应用大涡模拟(LES)方程对2种相同进、出口直径的喷管模型(轴对称、2元圆转矩)在Ma=0.8、高雷诺数(2×10~5)条件下进行了数值模拟计算。结果表明:与轴对称喷管相比,圆转矩喷管射流掺混效应增强,速度衰减快,核心区长度和高温区域面积减小。同时尾喷流拟序结构变化说明:2种喷管主要拟序结构均包含涡环、涡辫、发卡涡、螺旋涡等相似结构;但圆转矩喷管在射流近场诱导出的涡旋更丰富,边角剪切涡发展更快,形成明显的CVP结构,导致其射流柱失稳时刻提前、距离缩短;同时,喷管形式的改变使得射流剪切层内雷诺应力增大,速度脉动增强。拟序结构发展及雷诺剪切应力变化说明在射流流场中涡旋发展耗散速度增大、速度边界层脉动增强、射流柱易失稳是导致射流掺混增强的本质因素,为异形喷管的强化掺混机理提供了依据。     

不同等离子体激励强度下气膜冷却特性的大涡模拟研究

为进一步改善燃气轮机叶片气膜冷却效果,采用大涡模拟(LES)方法对不同等离子体激励强度情况下的平板气膜冷却流场进行了数值模拟研究。结果表明:与无等离子体激励时相比,等离子体激励强度逐渐增至10时射流出口最大流向与法向速度分别增大了16%和7%左右,并移向气膜孔的尾缘,而气膜孔前缘附近的法向速度约减小了4%,从而减少了射流迎风面上冷、热气流的掺混;等离子体对气膜孔下游回流区的动量注入效应使得回流区内的流向速度增大,抑制了横流绕流分离旋涡的发展;等离子体气动激励削弱了肾形涡对的强度及其抬升冷却射流的能力,从而提高了气膜冷却效率,中心线气膜冷却效率随激励强度的增大而升高,当激励强度为10时中心线气膜冷却效率最大提高了55%。     

固体火箭发动机转角涡脱落分类

应用三维大涡模拟(LES)方法对C1x模拟试验器进行了数值模拟,以研究转角涡脱落引起的声涡耦合现象.结果表明:与二维和轴对称LES相比,三维LES方法可以准确捕获转角涡脱落运动的三维特性及锁频现象.根据转角涡脱落是否在腔内形成声反馈循环过程可将其分为两类问题,第1类形成声反馈循环涡;第2类形成自耗散涡.燃烧室的台阶扩张比、台阶前后长径比是影响两类转角涡脱落现象的重要结构参数.     

亚声速等膨胀平面射流的初始流场结构

基于大涡模拟方法与5阶精度混合TCD/WENO格式,对等膨胀射流在马赫数为Ma=0.6条件下的初始流场结构进行了数值模拟。数值结果描述了射流等膨胀过程中,射流初始流场的形成、变化及发展过程。并发现了射流初始流场中涡核为射流失稳的起始点,研究了射流混合层与涡核剪切层上涡的产生、发展与破碎以及它们卷吸环境气体并与射流气体混合的整个过程,揭示了射流初始过程中涡核与涡环直径的变化趋势,得到了射流轴向速度与压力的分布。数值结果可为相关发动机的设计提供重要参考。     

亚声速等膨胀方管射流轴置换现象的数值研究

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度WENO/TCD混合格式,对Ma=0.6的亚声速等膨胀方管射流的初始流场进行了数值研究。数值结果清晰地描述了初始主涡环的形成与三维演变过程。方管射流具有平面管与圆管射流的一般流场特征,但方管周向曲率的不一致导致轴向流动速度不均匀,使方形主涡环出现复杂的Biot-Savart自诱导变形,并在方形涡环截面上诱导形成4对反向流向涡对。在主涡环的Biot-Savart自诱导变形与流向涡对的诱导速度共同作用下致使主涡环截面形状相对于初始时旋转45°,完成轴置换。另外,在亚声速等膨胀射流中,流向涡对的诱导速度在轴置换中占主要作用。     

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

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

二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

突扩燃烧室流场RANS与LES模拟

以中心突扩燃烧室为模拟对象,分别从旋涡运动、空间速度分布和时间速度脉动等方面对RANS模拟结果、LES模拟结果和试验结果进行了比较分析。表明RANS方法没有模拟出燃烧室中的旋涡运动,与试验结果相差甚远;而LES方法则比较好地模拟出了燃烧室中旋涡运动和速度脉动。得出LES方法是对突扩燃烧室内流场数值研究的有效方法。     

MVG控制斜激波/湍流边界层干涉的大涡模拟

采用浸没边界法(IBM)对带有微型涡发生器(MVG)控制器的激波/湍流边界层干涉流动进行了大涡模拟(LES)。以来流马赫数为2.3的斜激波(由平板上方8°楔产生)为基本流动入射平板湍流边界层,通过在干涉区前布置MVG阵列来控制激波诱导的边界层分离。采用浸没边界法处理MVG的复杂几何,分析了MVG尾迹区平均流速度剖面,雷诺应力,瞬态旋涡结构。结果表明:时均流场显示MVG尾迹区存在一对对转的主流向涡,流向涡加剧了边界内的动量交换从而增加了边界层抗分离能力,而瞬态流场则反映出MVG尾迹区的剪切层由于Kelvin-Helmholtz(K-H)不稳定性会卷起为一列展向旋涡。     

LESS燃烧室非定常旋流流动

结合实验和CFD数值方法,研究了LESS燃烧室旋流器出口附近的非定常流动特征.一维热线测试结果表明:预燃级出口存在着2249Hz周期性速度振荡,而主燃级出口流动未见明显脉动.通过瞬态LES数值结果检验了在旋流器出口附近采用一维热线测量三维流动速度偏差为8%.另外,数值模拟计算得出预燃级出口速度脉动频率为2331Hz,与热线测试结果偏差为3.6%.最后通过频谱分析和相关性分析,得出周期性速度脉动由预燃级内进动涡核主导.      

低旋流预混燃烧稳燃机理的大涡模拟

采用大涡模拟(LES)方法对比了高低旋流的预混燃烧特性,并分析了低旋流的稳燃机理。结果表明:LES耦合部分搅拌反应燃烧模型(PaSR)可以准确的捕捉旋流预混燃烧的细小火焰结构和流场分布。不同旋流强度引发不同的流场特性,低旋流诱发较弱的回流区,而高旋流会诱发很强的回流区,这种强烈的流场特征往往会引起振动频率极低的进动涡核拟序结构(PVC)。通过涡脱落和燃烧温度分布的瞬时耦合关系,发现火焰剪切层引起的内外涡脱落频率差是低旋流火焰稳定的最主要因素,低旋流燃烧受涡旋结构的影响交替出现W和V型火焰。     

Kármán涡街中旋涡三维变形的初步研究

本文用涡动力学及LIA方法作为基本理论模型,数值研究了尾迹中孤立涡和细涡丝的三维演化规律。结果表明,圆柱分离尾迹中的旋涡存在三维不稳定性。卡门涡在尾迹平均流场中演化产生三维的类似于马蹄形-勺子形的流场结构。细涡丝在涡辫区的三维演化形成趋向流场最大拉伸变形方向的流向涡结构。     

间断Galerkin方法计算叶型转捩流动

使用间断Galerkin方法研究叶型转捩流动,进行了大涡模拟(LES)和转捩模型的求解,对T106低压透平和Zierke-Deutsch压气机叶栅内的流动进行了计算.通过T106低压透平的计算对LES和转捩模型进行了比较,结果表明两种方法得到的压力分布和分离泡位置均与实验值吻合良好.LES得到的分离泡的轴向位置为0.145~0.165m,转捩模型得到的分离泡的轴向位置为0.150~0.156m.LES可以再现复杂的瞬时流动细节,对于深入研究流动机理很有意义,而转捩模型尽管无法获得足够的流动细节,但是也能预测边界层的分离和转捩现象,并且结果与LES时均结果相差不大,对于工程应用很有价值.使用转捩模型计算Zierke-Deutsch压气机叶栅内的流动也获得了与实验值符合的结果.      

对转桨扇三维流场特性数值研究

对转桨扇（Contra-rotating propfan,CRP）是下一代民用航空推进备选方案开式转子发动机最重要的气动部件,其气动性能对整机性能影响显著。本文对不同进距比下的对转桨扇三维流场进行数值模拟,结合压气机及螺旋桨相关理论分析了对转桨扇内部流场及其滑流区涡结构和滑流特征。结果表明,对转桨扇后排流动特征及性能参数变化幅度均超过前排。对转桨扇实际进口气流角受到轮毂附面层、诱导速度、抽吸效应的共同影响,可根据不同叶高位置轴向速度的分布规律判断三种影响分别起主导作用的位置。在对转桨扇滑流区中,桨尖涡是导致损失的主要原因,径向涡量衰减比周向和轴向涡量衰减更快。对转桨扇滑流在径向上影响至3.5倍叶高位置。气流出后排桨扇后会持续加速直至静压达到环境压力,加速区域长度约为5倍桨扇半径。     

k-ε（RNG）、LES模拟建筑物周围气流特征及湍流扩散的风洞试验验证

采用k-ε（RNG）与LES湍流模型在来流与建筑物迎风侧呈不同角度的情形下,模拟了位于立方体建筑物顶部污染源所排放污染物的流动和扩散规律,并与相应的风洞试验结果进行了比较。流场分析结果表明：数值模拟能够较好地模拟建筑物顶部回流、背风侧空腔区以及再附着点等。浓度场分析结果表明：来流与建筑物成45°时,建筑物顶部回流区与背风侧空腔区的数值模拟结果略低于风洞试验结果;来流与建筑物成90°时,建筑物顶部回流区数值模拟结果略高于风洞试验结果,而背风侧空腔区的数值模拟结果与风洞试验结果基本一致。综合分析表明：建筑物周围的流场影响浓度场的分布,LES、k-ε（RNG）模型都能够较好地模拟建筑物周围的流动和扩散规律,两种模型相比,LES模型与风洞试验吻合得更好。总之,风洞试验和数值模拟相结合能较好地研究建筑物对流动和扩散的影响。     

空调器室外机轴流风机内部流动的PIV研究

采用粒子图像测速仪PIV,对具有半管道式结构特点的空调器室外机轴流风机内部流场进行了实验研究,并结合实验结果分析了叶片顶部的叶尖涡和叶片出口尾缘涡的流动特性。实验结果显示在轴流风机流道内部叶顶区域存在与叶轮旋转方向相反的叶尖涡结构。叶尖涡产生于叶片前缘叶顶近吸力面侧,在流道内部与主流发生干涉后朝向周向和出口传播并逐渐耗散。叶尖涡涡心轨迹与叶顶弦长方向的夹角为10°,在叶高方向上叶尖涡的径向位置并不固定。与普通管道内部流动不同,叶片顶部与导风罩间的间隙中未捕捉到明显的叶顶泄漏涡现象。叶片出口近尾缘处30％以上叶高明显捕捉到尾缘涡结构,叶片压力面和吸力面侧的径向速度存在明显的方向变化,切向速度在尾迹区增加。