基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.      

翼型几何模型对边界层转捩预测精度的影响

对飞行器表面的边界层转捩现象进行准确地预测,是设计高性能层流翼型的关键。基于边界层数值模拟和现代CFD技术的特点,选取与SSTk-ω湍流模型相耦合的γ-Reθt转捩模型,对不同类型的翼型边界层转捩进行预测。结果表明:在转捩计算过程中出现的摩阻曲线震荡,是由翼型表面计算网格不光顺所引起的;提高翼型几何模型光顺程度能够在不同程度上提高翼型气动特性的计算精度,并能得到更加准确的转捩预测结果。     

隐式紧耦合SST和TNT湍流模型的高速流动数值模拟

将SST(shear stress transport)和TNT(turbulent/non-turbulent)湍流模型输运方程与平均流场控制方程进行隐式紧耦合求解,结合当地时间步长方法和湍流源项隐式处理确保求解过程的快速和稳定.采用AUSMPW+(AUSM by pressure-based weight functions)格式和LU-SGS(lower-upper symmetric Gauss-Seidel)隐式紧耦合方法对高超声速压缩拐角流动、锥柱裙流动和超声速非对称激波/边界层干扰问题进行了数值模拟.计算结果与实验值的对比表明:SST模型和TNT湍流模型可以很好地预测15°压缩拐角流动的壁面压力和热流密度;随着压缩拐角的增大,计算结果与实验值偏差增大;可压缩性修正对压缩拐角流动的压力和热流密度分布有很大影响,对超声速非对称激波/边界层干扰流动影响很小;隐式紧耦合方法比显式紧耦合方法具有更好的收敛特性.      

结构因素对离心通风器性能影响的数值研究

在高转速下,航空发动机离心通风器内腔存在强烈的湍流流动,而稳定工作后通风器内腔中漩涡呈周期性产生和淹灭。小尺寸颗粒受流体微团湍流随机脉动的影响,产生相对于平均流的随机脉动运动。为了研究通风孔偏心距和辐板顶圆半径对通风器性能的影响,采用离散相模型(DPM)模拟通风器内颗粒的运动轨迹。并应用随机游走(DRW)模型模拟连续相湍流瞬时速度脉动对颗粒轨迹的影响,采用随机涡寿命模型确定随机追踪模型的积分时间。结果表明:辐板顶圆半径的增大有助于提高通风器的分离效率,同时也增大了腔内流通阻力;通风孔偏心距对减小通风阻力的作用明显,但降低了离心通风器的分离效率。     

高升力翼型复杂流动模拟中湍流模型性能评估

采用软件Fluent中工程常用的7个涡黏湍流模型对某种高升力3段翼的降落阶段绕流进行了数值模拟.通过和试验及经过试验校验过的延迟脱体模拟（delayed detached eddy simulation,DDES）结果进行详细的对比分析,包括翼型压力系数、马赫数、涡量场和湍动能分布等,系统地研究了常用涡黏湍流模型对该高升力翼型的模拟性能.结果表明：对于平均流场,standard k-ω模型的模拟性能最好,能较好预测翼型压力系数、襟翼分离区位置和大小等;SST （shear stress transport）k-ω模型性能也较好,一方程SA （Spalart-Allmaras）模型和四方程v²-f模型具有相近的性能,而k-ε系列模型性能相对较差.对于湍流场的模拟性能,和上述平均流场对比所得结论相一致,但所有湍流模型均未预测出襟翼分离区附近湍动能“最大”的分布特征.     

考虑模型抖振力跨向不完全相关性效应的气动导纳识别

首先通过考虑节段模型上抖振力跨向不完全相关性效应推导了节段模型风洞试验中作用在模型断面上的分布抖振力谱和由底支式天平测到的模型总抖振力谱之间的关系。然后,以准平板断面为例,进行了格栅湍流场节段模型的测力和同步测压试验,获得了模型总抖振力谱以及模型抖振力跨向相关性函数。接着,采用等效导纳法以及抖振力自谱和抖振力脉动风速交叉谱综合残量最小二乘法分别识别了准平板节段模型等效气动导纳和六分量气动导纳,讨论了模型抖振力跨向不完全相关性效应对气动导纳识别结果的影响,并把识别得到的准平板断面气动导纳与平板断面气动导纳的理论结果——Sears函数进行了比较。结果表明：忽略抖振力跨向不完全相关性效应（即假设节段模型分布抖振力沿跨向完全相关）而直接采用平均抖振力作为断面上分布抖振力的传统方法会导致气动导纳识别结果偏小,并且,其偏小程度会随着频率的增加而增加;此外,相对于抖振升力和扭矩相关的气动导纳分量而言,由于抖振阻力的跨向相关性比抖振升力和扭矩的跨向相关性显得更弱,因此抖振力跨向不完全相关性效应对阻力相关气动导纳分量识别结果的影响更大;利用识别得到的六分量气动导纳反算的作用在模型上的分布抖振力谱与试验中实测结果非常接近,经抖振力跨向不完全相关性效应修正后的竖向脉动风速对应的升力和扭矩气动导纳分量的识别结果与Sears函数也比较接近,从而验证了用于六分量气动导纳识别的自谱-交叉谱综合最小二乘法的可靠性。     

固体火箭内流场的时均和大涡模拟技术对比研究

分别采用时均RNG湍流模拟和大涡模拟(LES)方法对某型固体火箭发动机内流场进行了瞬态数值模拟研究,建模过程中针对推进剂燃烧以及燃面边界采用了修正壁面通量的方法,并使用了动网格技术。结果表明:大涡模拟技术可以得到发动机内湍流流场中详细拟序结构,时均湍流模拟对于火箭发动机内流场时均量的模拟仍有较好的工程精度。     

直升机附件舱流场和温度场的数值仿真研究

建立了某型直升机发动机附件舱内流动和传热的数学模型,利用 k-ε二方程紊流模型模拟舱内流动,应用非结构化网格与网格自适应技术进行区域离散化,采用有限容积法对附件舱内的空气流场和温度场进行了三维数值仿真和分析,为附件舱冷却系统的改进和优化设计提供了依据。结果表明,附件舱冷却系统可以满足舱内冷却要求。     

湍流模型对发动机排气系统流场和红外辐射特征计算的影响

计算分析了采用k-ω SST(shear-stress transport)模型、RNG(renormalization group)k-ε模型配合标准壁面函数(RNG k-ε+标准壁面函数)、RNG k-ε模型配合增强型壁面函数法(RNG k-ε+增强壁面函数)3种湍流模拟方法,对典型涡扇发动机轴对称收敛喷管排气系统流场、温度场、组分质量浓度场及红外辐射特征计算结果的影响.研究表明:对于这3种方法,排气系统腔体内壁面压力分布基本一致;但在温度场、组分质量浓度分布以及红外辐射特征方面的差别较大.RNG k-ε+增强壁面函数计算得到的内涵壁面与燃气流的温度、CO₂与H₂O的质量浓度和红外辐射强度值最大,RNG k-ε+标准壁面函数的计算结果其次,而k-ω SST 模型的计算结果最小.通过与试验结果对比,发现k-ω SST模型对排气系统红外辐射特征的模拟更加合理一些.      

基于数值模拟的复杂地形风场风资源评估方法

与采用线性模型的WAsP软件相比,CFD方法具有可以逼真模拟复杂地形三维风场大气流动的优势,因此成为了复杂地形风场风资源评估的发展方向。为了把CFD软件包FINE/TURBO和带有壁面函数的k-ε湍流模型用于复杂地形地貌风电场的风资源评估和微观选址中去,本文以南澳岛风场作为研究对象,在中性大气条件假设下,对风场区域以30°风向为间隔进行数值模拟。研究采用不同入口边界条件对数值计算结果的影响。比对模拟结果发现,风场中的风加速与湍流强度基本上不随来流条件发生变化。基于这一规律,结合全年测风数据评估风场的风资源分布状况,建立了一套风资源评估的数值模拟及后处理方法,为工程实际提供了有力参考。