后台阶流动的Hybrid RANS/LES模拟

采用雷诺平均的Navier-Stokes和大涡模拟混合的方法对后台阶流动进行了数值模拟,并与实验结果进行了对比;分区尺度采用IDDES(Improved Delayed Detached Eddy Simulation)模式:来流包含背景湍流度时,混合方法转化为WMLES(Wall Modeling in LES)模拟,除了近壁区外的大部分区域都是直接求解;不包含来流湍流度时,转化为DDES(Delayed Detached Eddy Simulation)模式,附面层采用纯RANS模拟,大分离流动区用类脱体涡模拟方法模拟;对流通量的离散采用迎风型低耗散通量分裂格式和中心型四阶偏斜对称格式组合的模式;湍流模型采用Spalart-Allmaras(S-A)模型;亚格子长度尺度取决于网格分布和到壁面的距离.      

XY-SAS模型对于分离流动的性能分析

XY-SAS(Xu-Yan-Scale Adaptive Simulation)模型是一种新的基于平均涡量的一方程SAS类模型,不但解决了分离涡模拟DES(Detached Eddy Simulation)的雷诺平均N-S方程RANS(Reynolds Averaged Navier-Stokes)-大涡模拟LES(Large Eddy Simulation)交界面问题,而且计算量远小于早期SAS.通过对XY-SAS模型的机理分析,认为该模型克服了早期SAS模型预测的湍流能谱在高波数衰减不足的缺陷,并在分析的基础上对该模型进行了改进,取消了其限制器.选用平板边界层和方柱绕流作为算例,进一步研究了该模型对于稳态流动和大分离流动的性能,数值模拟结果同实验值符合良好,证明XY-SAS模型自身完备,是预测稳态流动和大尺度分离流动的有效方法,实现了涡粘性系数根据当地流动状态动态地自适应调整和光滑过渡.     

XY-SAS模型中不同网格尺度限制器的影响分析

XY-SAS(Xu & Yan Scale Adaptive Simulation)模型通过在SA(Spalart-Allmaras)模型耗散项的湍流尺度中引入基于平均涡量的冯卡门长度尺度,在不稳定流动区域具有类似LES(Large Eddy Simulation)的计算性能,且不存在DES(Detached Eddy Simulation)方法中RANS(Reynolds Averaged Navier-Stokes)与LES的交界面问题.由于SAS类模型预测的湍流能谱存在高波数衰减不足的缺陷,需要进一步在湍流尺度中引入网格尺度作为限制器.采用不同形式的网格尺度限制器构造了3种XY-SAS模型,对Re=3 900的圆柱绕流进行数值模拟,并与DES及实验数据进行比较.结果表明网格尺度限制器对XY-SAS模型的计算性能有较大影响,以最小网格间距作为高波数衰减限制器是较为合理的选择,此时模型对于大分离流动能够达到DES的模拟能力.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-Stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA(Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

DES方法对返回舱气动热环境数值模拟

脱体涡模拟方法(DES,Detached Eddy Simulation)由于结合了湍流模型(RNAS,Reynolds Averaged Navier-stokes)和大涡模拟(LES,Large Eddy Simulation)两者各自优势,是模拟非定常大分离流动的有效方法.采用基于SA (Spalart-Allmaras)湍流模型的DES方法,对高超声速返回舱外形进行数值模拟,计算所得的返回舱表面压力及热流分布与实验结果吻合,证实了DES方法的优越性.最后,在DES方法中加入可压缩修正,结果证实在高超声速流动中,可压缩修正方法在一定程度上能够提高原始模型的预测性能.     

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

基于湍流尺度的混合RANS/LES模型

比较BL(Baldwin-Lomax)湍流模型的内层模型与Smagorinsky亚格子模型,两者的主要差别在于湍流尺度的计算方法不同,而RANS(Reynolds Average Numerical Simulation)方程组与LES(Large Eddy Simulation)方程组在形式上又是相同的.基于这两点,发展了RANS/LES混合模型.该模型是在网格滤波尺度与到壁面的距离相等处将流场分区,在近壁区直接求解RANS方程组,在远离壁面区域,求解LES方程组.为了保证湍流粘性系数在内外层分界处的光滑过渡,内外层湍流粘性系数通过混合函数连接.混合方法减轻了对壁面附近网格分辨率的要求和时间步长的限制,因此比LES所需的计算成本少.用该混合模型模拟了绕双椭球的高超音复杂流场,结果表明采用混合模型可以准确地模拟可压缩湍流流场.      

DES与DDES在湍流分离中的原理与性能研究

随着工程上流动结构的日益复杂，兼具雷诺平均Navier-Stokes(RANS)方法高效率与大涡模拟（LES）高精度的分离涡模拟（DES）类混合方法成为现阶段工程中最有效的湍流模拟方法之一。围绕DES类混合方法中的DES与延迟分离涡模拟（DDES）方法开展工作，分析二者开关函数构造上的不同，研究延迟因子作用机理，并考察DES与DDES方法的求解能力。研究表明：DES与DDES方法在模拟表现上存在一定差异，DDES方法通过引入延迟因子，保护RANS求解区域，改善模化应力不足，降低了DDES方法对交界面系数CDES敏感程度；DDES方法在计算过程中容易出现过度保护，导致求解瞬时涡结构能力不如DES方法，分析与延迟因子引入比重及开关函数构造形式有关。     

LES／RANS混合方法对超声速底部流动的应用

为克服大涡模拟(LES,Large Eddy Simulation)对近壁区域使用的局限,近年来LES/RANS(Reynolds Averaged Navier-Stokes)混合方法成为研究热点,利用LES与RANS两者自身的优点克服对方的不足,有效实现计算精度和效率的统一.基于LES/RANS混和方法,采用Smagorinsky-BL(Baldwin-Lomax)模型和新型Vreman-BL模型,对超音速底部流动这一难题进行了模拟,分析了其流场结构.计算结果表明:得到的速度型、底面压力分布与实验值吻合,优于雷诺平均方法的结果.     

尺度自适应模拟的网格尺度关联分析

采用尺度自适应模拟（SAS）对雷诺数3 900的圆柱绕流展开数值模拟，对比分离涡模拟（DES）的计算结果和已有文献中的实验数据，系统研究了SAS和网格尺度的关联性问题。详细研究了不同网格分辨率和展向计算域的影响，分析了冯卡门尺度在尾迹区的时均湍流统计特性和瞬时分布规律。结果表明:在相同的网格分辨率下，SAS预测的回流区长度小于大涡模拟（LES），表现出较早的剪切层失稳；网格加密后，SAS预测的回流速度增大、雷诺应力峰值降低，计算结果与Lourenco & Shih的实验结果相接近。此外，在相同网格分辨率下改变展向计算域大小对SAS结果的影响很小。对SAS的网格尺度关联分析可以为该方法的工业应用提供指导。      

基于释热率分析的钝体贫油熄火过程分析

针对锥形钝体稳定的甲烷-空气预混湍流火焰复杂的熄火过程，采用大涡模拟(LES)与输运方程概率密度函数(TPDF)湍流燃烧模型相结合的模拟方法，研究远离熄火、近熄火及熄火点的火焰和释热率(HRR)数值变化情况，定量分析熄火判据。结果表明:冷态速度模拟结果和实验的相对均方根误差在10%以内，热态误差在20%以内；释热率是否出现在OH和CH₂O重叠的区域，是判断熄火的一个重要参数；远离熄火时，释热率高的区域主要沿内侧剪切层出现；近熄火工况下，释热率在流向轴上闭合，回流区下游也出现较大的释热率；熄火点工况下，较大释热率的区域在回流区下游和上游均有出现；模拟预测的熄火情况和实验PLIF结果一致；平均释热率可作为判断熄火的定量依据，即当钝体后方0.2d处内侧剪切层平均释热率与回流区平均释热率的比值小于4时，发生熄火。      

战斗机武器外挂投放与内埋投放比较

为了研究弹体外挂投放与内埋投放的区别,利用基于Menter SST湍流模式的改进延迟分离涡模拟(IDDES)方法以及重叠网格技术,分别对亚声速和超声速来流条件下,同一弹体外挂投放和内埋投放进行了数值模拟,得到了亚声速和超声速条件下外挂投放与内埋投放弹体的下落规律。通过对比分析表明:亚声速和超声速来流条件下,内埋投放由于受舱体内强非定常流场以及舱体唇口剪切层的影响,弹体受很大的抬头力矩,弹体姿态角变化较大,投放特性劣于外挂投放。进一步研究表明:在弹射时给弹体一定的低头角速度,使弹体以低头姿态穿越剪切层,则可以大幅度降低剪切层带来的不利影响,提高内埋投放弹体分离品质。     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.      

高精度激波捕捉格式的性能分析

在数值模拟激波-湍流干扰问题时,为了精确模拟湍流脉动的小尺度结构,数值格式应尽量减小自身的数值耗散.但是为了稳定、无震荡地捕捉激波,数值格式又必须具有一定程度的数值耗散.基于上述这2种相互矛盾的要求,采用无黏泰勒格林涡、可压缩湍流等经典算例对目前广受好评的几种高精度格式进行了系统地对比.一系列的计算分析表明,当前工程中广泛应用的高精度格式虽然都能一定程度地满足上述2种要求,但与理想目标仍有较大的差距.     

改进SA模型对翼型分离流动的数值模拟

采用(CFD,Computational Fluid Dynamics)软件Fluent,选取(SA,Spalart-Allmaras)模型,对风力机翼型S825进行了二维数值模拟,并与实验数据进行比较.针对SA模型不能准确预测翼型尾部流动分离的情况,分析了分离区域内湍流的强非平衡输运特性.研究表明,SA模型中系数C_b1与湍流的非平衡输运特性密切相关,进而提出了修正模型系数C_b1改进SA模型的方法,结果表明该方法能够更准确地模拟风力机翼型分离流动.