使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.      

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

用GAO-YONG湍流模式对分离流动的数值模拟

采用GAO-YONG湍流模式对二维管道内三角形钝体绕流问题以及槽道内后台阶分离流动进行了数值模拟.求解程序基于开源数值计算平台OpenFOAM (Open Field Operation and Manipulation),数值模拟结果很好地预测了钝体绕流问题中漩涡脱落的尾流的流动趋势以及后台阶流动中分离再附的流动结构,同时分析了速度分布以及摩擦系数等参数并与实验值进行了对比,结果符合很好.这表明GAO-YONG湍流模式对大分离流动有较好的预测能力,对工程实践具有指导性作用.     

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

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

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

使用Gao-Yong湍流方程计算翼体角偶流动

Gao-Yong湍流方程基于侧偏统计平均方法保留了湍流脉动量的一阶统计信息,并引入加权漂移速度对称性及正交各向异性,导出漂移流连续方程、动量方程和机械能方程,最后依据湍流物理的唯象论,使用动量传输链概念模化封闭了整个方程组.大量算例验证了此方程对广泛范围的复杂湍流问题的适定性.为了进一步的工程应用,将其加入开源计算流体力学软件OpenFOAM(Open Field Operation and Manipulation)中,并就翼体角偶流动进行了数值模拟.翼体角偶流动具有三维分离、马蹄涡等十分复杂的流动特征,通过与实验测量结果的对比研究表明,Gao-Yong湍流方程不论在定性还是定量上,都成功捕捉到这类复杂流动现象的主要特征,得到了令人满意的结果.     

湍流燃烧反应机理的简化

采用速率敏感度分析法分别对燃料为CO/H2和CH3OH的钝体绕流湍流扩散火焰的燃烧机理进行了简化.判断出了冗余组分,并通过对敏感度矩阵F的主成分分析,找出了次要反应,使原始的复杂湍流燃烧机理得到了大幅度的简化.简化机理与详细机理的计算结果吻合得很好.结果表明,速率敏感度分析的方法能大大简化复杂的湍流反应机理,且计算效率高,精度满足要求.     

基于熵限制的Baldwin-Lomax湍流模型

针对Baldwin-Lomax(BL)模型在模拟超声速复杂流动中的不足,提出一种基于熵限制的新型修正方式.分别采用原始的BL模型(BL-origin)、修正后的BL模型(BL-entropy)以及Wilcox k-ω两方程模型对超声速平板、超声速二维压缩拐角以及膨胀-压缩拐角3个典型算例进行计算,结果表明BL-entropy 能够较好地获取长度尺度进而得到均匀合理的涡粘性分布.同时,该修正方式简单高效并且适用于其他与超声速边界层有关的模型,具有较大的发展潜力.     

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

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

多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

各向同性湍流能量级串中的旋涡分岔机制

充分发展各向同性湍流能量级串和多尺度相互作用一直是湍流理论研究的核心问题.目前,对于该物理过程的完全理解或精确的数学描述缺乏基于第一原理的理论.简要介绍了湍流能量级串的概念、起源、发展历程及面临的挑战问题,着重阐述了目前各种现有描述方法的局限性.基于三维不可压缩流体的Karman-Howarth方程,根据新得到的各向同性湍流尺度演化方程以及在这一方向上的理论进展,证明存在以湍流Taylor微尺度为动力学量的非线性动力系统.根据上述新的理论,可以认为:湍流能量级串由一系列的旋涡非线性分岔过程刻画,呈现Feigenbaum倍周期分岔的途径.     

通气位置对潜射航行体流体动力特性影响分析

基于均质平衡流理论,通过求解混合介质的RANS方程、SST湍流输运方程和各相之间的质量输运方程,开展了通气位置对潜射航行体流体动力特性影响的三维数值模拟研究,对比分析了不同通气位置条件下空泡形态特性、表面压力分布以及阻力变化特性.结果表明:当两个通气口间距增加到一定值后,在两个通气口之间区域空泡发生断裂,在断裂区域迎流面和背流面表面压力值升高,且在空泡断裂闭合位置出现压力峰值;在相同通气量条件下,随着通气口间距的增大,压差阻力系数和粘性阻力系数均呈减小趋势变化.     

谐波平衡法在低速非定常流模拟中的应用

谐波平衡法是一种有效的周期性非定常流的计算方法.采用基于可压缩流的谐波平衡方程在计算低速不可压流动时,会由于对流通量计算格式中的数值粘性污染,降低解的精度和收敛性.采用预处理技术,使得基于可压缩流的谐波平衡方程可以直接用于低速周期性非定常流的计算中.选取典型的不可压方腔驱动流和低雷诺数圆柱绕流为例进行了时间推进法和谐波平衡法的计算对比.计算结果表明预处理后的谐波平衡方程适合于低速流的计算,在谐波平衡法中采用较少阶数的谐波计算就可以还原出几乎准确的非定常流场.