大涡模拟模型环形燃烧室污染特性

在任意曲线坐标系下,采用大涡模拟模型环形燃烧室两相燃烧流场中污染物的生成.用亚网格涡破碎(EBU)燃烧模型估算化学反应速率;分别用NO亚网格动力学模型与CO亚网格模型来模拟NO和CO的生成;采用随机离散模型模拟气液两相湍流流动,两相之间耦合采用PSIC算法.计算结果与瞬态速度场、出口温度和污染物质量分数分布的实验数据比较吻合,表明大涡模拟方法和亚网格模型可以用来预估实际燃烧室两相喷雾燃烧流场和污染物生成过程.      

超声速凹槽大涡数值模拟研究

在超声速燃烧室内常采用凹槽增强燃料与空气的混合,实现火焰稳定,本文采用大涡模拟方法对开式凹槽流场进行了数值模拟研究,计算采用Smagorinsky 亚网格模型模拟小尺度涡的作用.数值模拟结果给出了凹槽流动的自激振荡的发展过程,以及凹槽流动中拟序结构的演变过程,如涡的卷起、增长,涡与涡之间的合并、破碎等,计算得到的凹槽波动的峰值频率和理论结果一致.     

大涡模拟模型燃烧室燃烧性能计算

对带双级扩压器的模型燃烧室气液两相瞬态喷雾燃烧过程,在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法进行数值研究,同时采用多维经验分析法预估燃烧性能.采用 k 方程亚网格尺度模型模拟亚网格湍流黏性;亚网格EBU(eddy-break-up)燃烧模型预估化学反应速率;多维经验分析法计算燃烧性能;并在非交错网格体系下气相采用SIMPLE(semi-implicit method for pressure-linked equations)算法对控制方程进行求解,液相采用随机离散模型,两相之间的耦合采用PSIC(particle-source-in-cell)算法.通过大涡模拟瞬态及时均计算结果表明:与粒子图像测速仪(PIV)测量的瞬态速度场、出口温度分布试验数据吻合,表明在三维贴体坐标系下采用欧拉-拉格朗日两相大涡模拟方法,数值模拟模型燃烧室两相喷雾燃烧流场,所采用的亚网格模型可以用于燃烧室气液两相喷雾燃烧流场的大涡模拟;燃烧性能计算结果与试验测量结果基本一致,说明所采用多维经验分析法可以用来数值模拟航空发动机燃烧室燃烧性能的计算,特别是污染物的预估,为设计低污染高性能航空发动机燃烧室提供有用的设计依据.      

DDES方法在复杂旋翼流场计算中的应用

应用延迟脱体涡模拟（DDES）方法,进行了复杂分离流动旋翼流场的数值模拟研究。为控制计算规模,更好地模拟旋翼桨叶近壁面附近区域流动分离及远场尾迹发展,采用了贴体与背景自适应直角网格嵌套的网格系统;同时为提高计算速度,采用了基于分布式存储的数据并行模式,实现了贡献单元搜索、洞切割和流场计算的并行化。计算工作包含一个共轴双旋翼悬停和单旋翼下降状态,目的是分别讨论DDES方法在旋翼表面附近空间区域是否存在大范围分离时计算结果与雷诺平均Navier-Stokes结果之间的差异。计算表明,对于下降状态,DDES方法预测出了更强的尾迹区分离流动;对于桨叶表面基本为附着流动的双旋翼中等桨距角情况,两种方法预测的流场结构相似,仅在下旋翼桨叶气动力和旋翼桨根下方空间区域的诱导速度计算结果上存在微细差异。     

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

脊状表面流场数值计算方法研究

根据现有湍流数值模拟方法的特点和脊状表面减阻理论研究的实际需要,研究了基于雷诺平均N-S方程的脊状表面流场数值计算方法,并基于纵向及横向脊状表面的形状特点,提出了不同脊状表面计算域选取及网格剖分的方法。一系列典型横向及纵向脊状表面流场数值计算结果表明,所采用的这套脊状表面数值计算方法不仅能够精确获得脊状表面的减阻特性,还可以获得脊状结构内部及附近流速、压强、涡量等细微流场结构的分布规律。     

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

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

小型涡喷发动机整机流场数值模拟技术研究

为计算某小型弹用涡喷发动机流场,研究了针对涡喷发动机整机的建模、网格处理以及动态工作点确定等关键技术。通过复杂结构部件建模和拼接,建立了涡喷发动机整机流场计算网格,理论模型中考虑了叶轮流场、湍流和燃烧过程的模拟,并以作用转子上气动力矩为基础计算转子的转速。在通用计算流体力学（CFD）软件中,通过编写UDF程序构建了发动机流场的计算框架。结果表明：全机流场仿真技术为分析涡喷发动机工作过程提供了有力工具,同时揭示了几何建模的精度等因素对计算结果准确性有较大影响。     

三维Euler方程自适应八叉树结构各向异性直角网格算法

采用各向异性直角网格对相应的三维复杂流场进行了模拟 ,提出了物面三角化方法及网格与物面的相交判断过程。用以中心差分为基础的 Jameson的有限体积法对 Euler方程进行求解。计算结果表明 ,各向异性直角网格算法对相应流场进行模拟 ,计算时间较短而效果较好     

各向异性网格自适应计算在超燃模拟中的应用

采用各向异性网格自适应求解技术,将其应用于DLR超燃冲压发动机燃烧室中的超燃模拟.开展了3个算例,包括采用滑移壁面条件的燃烧室冷流场模拟、采用无滑移壁面条件的冷流场模拟及采用无滑移壁面条件的反应流模拟.模拟中,各向异性网格自适应计算捕捉到了如激波、射流、边界层、火焰面等具有各向异性特征的大梯度区域,并利用各向异性网格进行了很好的加密.对比利用各向同性网格的初始流场计算,各向异性网格自适应计算使基于滑移条件、无滑移条件的冷流计算及反应流计算的网格单元数量分别下降了36.2%,36.4%和36.8%,有效降低了计算规模,而且流场大梯度区域的计算结果更准确,辨析度更好.结果表明:对于像超燃这类具有各向异性特征的问题,各向异性网格系统比各向同性网格系统有更好的计算效率及准确性,同时也表明基于Mach数场构造的各向异性网格系统可以有效应用于超燃计算.      

翼尖涡流场特性及其控制

大型运输飞机的尾涡系是诱发后继小型飞机空难的重要原因,需要有效的涡控制装置来削弱其强度.通过风洞实验,研究了翼型为NACA23016的矩形半机翼模型翼尖尾涡流动结构和控制方法.应用七孔探针空间流场定量测试技术研究了翼尖涡的流动结构,给出了翼尖尾涡在下游两倍弦长距离内的速度和压力场分布随迎角变化的规律.在机翼翼梢布置不同组合方式的翼梢涡扩散器,来控制翼尖涡.研究结果表明,正负90°和60°安装角的双翼梢涡扩散器可将翼尖涡涡核的静压增加60%以上.其旋涡强度削弱机理为:翼梢涡扩散器将集中的翼尖涡破碎分成两个或多个强度更弱的旋涡.在流体粘性的作用下,旋涡能量耗散更快,可有效地削弱翼尖尾涡的强度.     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

双柱振动诱导定常整流旋涡流动的实验研究

本文用流动显示实验研究圆形静止不可压缩粘性流体域内以圆心对称布置的双柱作同频同相小振幅振动时诱导的Ｓｔｏｋｅｓ层外的定常整流旋涡流动。实验表明,串列双柱当柱间距足够大时整流流谱呈八涡结构。当间距很小或为零时上涡结构。双柱的斜置有利于整流旋的合并,并柱间距接近至一定范围时,整流流动呈三维特征。振幅的减小和频率的降低使流动易出现三维现象。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

三角翼前缘脱体涡和尾涡的卷起及其相互作用

本文应用二维非定常比拟和恪子涡(Vortex-in-Cell)方法,数值模拟了三角翼前缘涡层的卷起以及与尾涡的相互作用。由于使用了数千个点涡和较小的空间网格,获得了前缘涡层小涡配对出现的较大的离散涡,模拟了紧卷涡层的不稳定现象。这是以前类似的工作未曾实现的。     

栅格翼尾流中的发卡涡

介绍了在水漕中用氢气泡法显示的斜置蜂窝式栅格翼尾流流态 ,发现在中等迎角下 ,相交翼元的每个交叉处邻近壁面的尾流中 ,有伴随着壁面流动分离而周期性脱落的发卡涡 ,它们相互连接 ,形成“涡链”。还介绍了栅格翼尾流中发卡涡的主要流动特征 ,并对其产生的流动机理作了初步探讨。     

有限长圆管内刚体旋转流的直接数值模拟

通过对三维不可压缩粘性流体的直接数值模拟,研究了有限长圆管内刚体旋转流在初始扰动下,扰动的增长和动力学演化过程.首先,通过改变流动的Reynolds数和旋转角速度ω,分析了轴对称流动中Reynolds数和ω对流动失稳的影响.对于高Reynolds数轴对称流动,利用数值模拟得到的不稳定结果与无粘线性稳定性理论预测的结果一致.最后,为进一步研究流动不稳定发展的三维效应,又进行了完全三维流动的模拟.数值结果表明,三维流动中不稳定的主导模态为螺旋模态,模态在线性段指数增长,导致流动产生螺旋型的漩涡破裂,之后经过非线性段的增长后达到饱和,流动最终发展为轴对称泡型旋涡破裂的稳定状态.     

透平叶片顶部间隙流动特性的实验和数值研究

利用粒子图像测速技术（PIV）捕捉透平叶片顶部泄漏流特征,并以此数据验证湍流模型和用商业软件CFX12．0进行的数值模拟方法。所研究的叶片为典型的GE—E3叶片,为了展示泄漏涡的生成和发展过程,用实验数据展示了3个不同截面的速度分布。数值计算中使用了混合网格生成技术及5种湍流模型。通过与实验数据的比较发现：RNG肛∈模型计算所得的泄漏涡与实验所拍摄的真实流动能较好地吻合。此模型和计算方法同样适用于研究叶顶射流对泄漏流的影响。计算结果显示：通过叶片顶部气膜孔射流产生的阻挡效应,最多能降低6．12％的主流泄漏。     

不同深宽比陷窝诱导涡结构分析

陷窝所诱导的旋涡强化动量或热的交换,使得陷窝在分离流动控制和强化对流传热研究中很受关注。深宽比是影响陷窝诱导涡结构的重要参数,开展了深宽比对陷窝诱导涡结构影响的研究。研究发现,深宽比为0.06时,涡结构为马蹄涡;深宽比为0.1至0.14时,涡结构为闭式分离泡;深宽比为0.16至0.20时,涡结构为对称类龙卷风涡;深宽比为0.22至0.30时,涡结构为非对称类龙卷风涡。通过涡核涡强度对比发现,过深的陷窝和过浅的陷窝其内涡强度均比中等深度陷窝弱。