面向液态铅铋数值传热的湍流热通量模型构建综述

液态铅铋(LBE)数值传热受制于极低普朗特数特性,传统雷诺比拟方法不能准确封闭和描述平均能量方程中的湍流热通量(THF)及温度输运过程,因此需构建针对低普朗特数流体THF的封闭模型。模型分为4类:湍流普朗特数模型、代数热通量模型、二阶矩微分热通量模型及利用机器学习的数据驱动模型。本文通过梳理其研究进展、建模思想,对模型所展现的复杂程度、在工程应用中的适用性等特点展开评述,为后续THF封闭建模研究、优化及创新提供见解。在燃料棒束间LBE的强制对流换热模拟中,显式代数热通量模型展现出较高的应用前景,机器学习则为THF封闭提供了全新的视角。

湍流燃烧二阶矩模型DNS检验研究进展的回顾

现有的湍流燃烧模型中,有的缺乏普适性,也有的用于大涡模拟的计算量太大.一种基于湍流流动模拟思路的统计矩方法被提出来,称为二阶矩(second-order moment,SOM)燃烧模型,用于非预混射流火焰和旋流火焰的雷诺平均模拟(Reynolds-averaged Navier-Stokes modeling,RANS)和大涡模拟(large-eddy simulation,LES),以及钝体后方预混火焰的大涡模拟,已经得到实验结果的良好检验.进一步的研究是该模型的直接数值模拟(direct numerical simulation,DNS)检验.本文对作者以及文献中报导的SOM模型DNS检验的研究进展进行了回顾.

多孔介质传热模型在多孔壁湍流中的适用性

为了考查不同的多孔介质传热模型在不同工况下的适用性,对带有高孔隙率多孔介质壁面槽道湍流及其传热进行了直接数值模拟研究。在多孔介质层外流体区域,通过有限差分方法求解不可压缩Navier-Stokes方程和温度对流扩散方程;在多孔介质层内,使用修正的Darcy-Brinkman-Forchheimer模型描述高孔隙率多孔介质阻力,以及分别采用局部热平衡(local thermal equilibrium,LTE)模型、局部非热平衡(local thermal non-equilibrium,LTNE)模型、理想金属(ideal metal foams,IMF)模型计算温度分布。通过对所得热场的统计特性的分析比较,探究了不同Biot数下水和空气两类流体介质的多孔介质传热模型的有效性。研究表明:LTE模型不足以准确预测金属泡沫多孔介质内传热问题,其等效导热系数因仅考虑孔隙率因素而低估了多孔介质层的传热能力;IMF模型在小比热容流体介质的算例中表现良好,可以代替LTNE描述多孔介质层内的传热,而在大比热容流体介质的算例中表现不佳,需要考虑比热容以及流固两相间的传热能力对预估的固体相温度分布进行修正。

基于物理信息神经网络的气膜冷却湍流模型反演学习

由于气膜冷却问题中湍流的复杂特性,传统雷诺平均(RANS)方法会低估湍流的热扩散强度,导致冷却效果计算不准确。对此提出了一套基于物理信息神经网络(PINN)的湍流建模框架,基于RANS流场和大涡模拟(LES)温度场,建立了数据驱动的湍流普朗特数神经网络模型,在RANS求解器中嵌入该模型,可以动态调整湍流的热扩散强度,获得了与LES高度一致的温度场。结果表明:PINN是构建数据驱动湍流模型的良好方法,对于湍流普朗特数的建模可以有效提升RANS方法对温度预测的准确性。

面向内流的激波/边界层湍流模型数据同化及应用

为研究压缩拐角激波/边界层干扰问题和机器学习方法在湍流模型参数辨识中的有效性,提出一种面向内流的激波/边界层湍流模型数据同化方法,以Kriging代理模型传播参数不确定量化过程,基于贝叶斯框架构建似然函数作为评判标准,最后利用粒子群优化算法近似获取参数的最大似然估计并进行参数验证。结果表明,通过校准大角度(24°)压缩拐角获取的湍流模型参数,可以应用到相同条件下相对较小的小角度(20°,16°和8°)压缩拐角,获取的壁面压力、摩阻系数和速度剖面均与试验值基本吻合。在Ma=2.85下校准的壁面压力,均方根误差由60.29%下降到16.56%。将大角度下获取的参数应用到Ma=2.9和不同的入射边界层厚度的条件下,获取的壁面压力和速度剖面仍与试验值基本吻合,验证了小范围马赫数内湍流模型参数的适用性。

湍流模型对风琴管喷嘴空化射流数值模拟的影响分析

针对不同湍流模型对空化射流流场数值模拟影响问题,基于ZGB(Zwart-Gerber-Belamri)空化模型和WALE(wavelet analysis of Lagrangian errors)亚格子模型,采用大涡模拟(large eddy simulation,LES)和RANS(Reynolds-averaged Navier-Stokes)模型对风琴管喷嘴空化射流流场进行数值模拟研究,分析内外流场的压力分布、速度分布以及空化云演变规律,并将空化云演变规律与实验结果进行对比分析。结果表明,LES模型速度、压力分布更符合实际情况,空化云演变规律与高速摄像拍摄结果基本吻合,LES模型可以更准确地模拟风琴管喷嘴空化射流流场,对风琴管喷嘴结构优化、打击力度与目标靶距的预测更准确。

基于直接数值模拟的压力容器入口射流湍流模型评价

采用基于谱元法的开源软件Nek5000对压力容器入口射流进行直接数值模拟(DNS)。压力容器模拟在环向简化为1/3的尺度,主管段雷诺数为63700,下降段雷诺数为14014,工质为碘化钠溶液。通过对常见物理量进行统计,构建高精度湍流数据库。并采用FLUENT对常见的壁面函数和湍流模型进行适应性评价,评价标准包括速度、湍动能和湍流黏度。结果表明:在上升段Stress-ω模型适用性最好;在下降段,Standard k-ω模型适用性最优。

棒束结构气冷换热的湍流模型适用性评价

高功率空间核电源采用棒束气冷堆具有质量轻等优势,其紧凑的栅格结构以及流动雷诺数(Re)较低的特点会影响堆芯的流动换热规律,采用CFD开展数值分析时需评价湍流模型的适用性。在气体工质棒束结构流动换热实验的基础上,利用ANSYS Fluent建立了试验段的数值模型,针对入口Re在688~2 986之间的实验工况,选择4种湍流模型开展了数值模拟,对比了加热棒包壳温度实验测量值与计算值。结果表明,4种湍流模型计算的棒温整体上均低于实验值,其中转捩SST模型结果与实验值最接近,整体平均偏差为-2.0%,较好地反映了横流等特征,可用于Re在2 000左右的棒束气冷堆芯热工水力数值计算。

笛卡尔网格下不同湍流模型的壁面函数方法研究

自适应笛卡尔网格湍流模拟时,网格数目呈雷诺数的指数方量级增长,带来网格规模“灾难”问题。壁面函数方法可以将壁面第一层网格尺度放宽1~2个量级,是有效的解决办法之一,然而,由于笛卡尔网格的非贴体特性,壁面函数的实现及其与不同湍流模型的匹配是关键。针对这一问题,从经典的Spalding壁面函数出发,采用显式赋初值的方式,解决了壁面摩擦速度迭代求解初值敏感及分层壁面函数速度不光滑问题;基于切应力恒定的壁面条件,构造了多层虚拟网格的切向速度和湍流粘性系数,实现了笛卡尔网格下的壁面函数方法;采用结构系综理论,通过DNS数据标定,构造了湍动能随湍流粘性比率变化的关系式,实现了壁面函数与KDO湍流模型的匹配,并以平板T3b、湍流凸包和跨音速NACA0012翼型绕流为算例,与常用的SST、SA湍流模型壁面函数方法进行了对比分析。结果表明:该文提出的KDO湍流模型壁面函数法准确、可靠,且因其形式简单、计算量相对较小,有较大的推广应用价值。

Menter SST湍流模型系数在高速流动模拟中的不确定度量化研究

Menter SST湍流模型系数的不确定性会对高速流动数值模拟的结果产生显著影响。本文针对二维压缩拐角、二维圆柱和三维双椭球体3个典型高速模型,定量评估湍流模型封闭系数不确定性对数值模拟的影响。首先,采用非嵌入多项式混沌方法量化由Menter SST湍流模型封闭系数不确定性导致的气动力/热数值模拟的不确定度;然后,依据Sobol敏感性指标对封闭系数的敏感性进行排序,辨识导致不确定度的关键因素。研究结果表明气动热的不确定度普遍大于气动力。在考察的9个封闭系数中,σ_(ω1)、β^(*)、κ、a_(1)是对数值模拟结果不确定度贡献较大的关键模型参数,其敏感性大小与当地流动特征相关:在流动的大部分区域内,生成项系数σ_(ω1)和κ是导致数值模拟结果不确定度的主要因素;而在流动发生分离或出现激波后的区域,数值模拟结果受耗散项系数β^(*)与剪切应力限制器参数a_(1)的影响更显著。

结合机器学习的SA湍流模型闭合系数修正

将修正Morris分类筛选法与极端梯度提升(extreme gradient boosting,XGBoost)相结合,在计算流体动力学(computational fluid dynamics,CFD)数据驱动下,用于SA(Spalart-Allmaras)湍流模型闭合系数的修正.利用分类筛选法有效缩小闭合系数研究范围,同时依据XGBoost方法在小规模数据集下取得精度较高的拟合模型,有效提升系数修正效率.在三维DLR-F6-WB构型下进行了数值实验,实验结果显示利用本方法能够在三维复杂模型上基于小样本数据进行系数修正,修正后的升阻力系数计算精度得到了显著提升.

COSINE两流体八方程子通道程序湍流交混与空泡漂移模型有效性验证

COSINE是我国自主研发的核电厂堆芯与系统安全分析的一体化软件包,其中子通道程序cosSubc-2F采用两流体,三相场(连续液相、液滴相和气相)的建模方法,除了连续液相和液滴相共用一个能量守恒方程外,其余各相均包含质量、动量和能量三个守恒方程,共计八个独立的守恒方程。在子通道程序中,燃料组件被划分为一定数目的平行子通道,这些通道并不是封闭的,它们之间由于压力梯度等因素的影响发生着流体的质量、动量和能量交换。作为通道间相互作用的关键模型,湍流交混与空泡漂移模型对于子通道程序模拟结果的准确性十分重要。文章详细介绍了子通道程序cosSubc-2F中的湍流交混与空泡漂移模型机理,选取一系列实验进行测试,程序中湍流交混与空泡漂移模型有效性得到了充分验证,并通过结果的对比分析,为模型的开发优化提出建议。

喷管两相湍流数值模拟及湍流模型性能评估

研究固体火箭发动机喷管中颗粒质量流率对两相流的影响规律以及不同RANS模型对该问题的预测性能,可为喷管设计等工程应用提供重要参考。在欧拉-拉格朗日(Euler-Lagrange)框架下,通过3D大涡模拟(3D LES)研究了颗粒质量流率对喷管两相流场的影响规律,以3D LES结果作为参考依据,分析了不同RANS模型针对喷管两相流在2D轴对称和3D模拟中的性能。结果表明:由于扩张段内颗粒集中在中轴线附近区域,颗粒质量流率越大,该区域马赫数越低,温度越高;壁面附近存在无颗粒区,故边界层内流动几乎不受颗粒影响。另外,采用不同湍流模型计算得到的颗粒分布与LES基本相同。在2D轴对称RANS模拟中,发现RNG k-ε模型所预测的气相质量流率和喷管比冲与LES结果吻合最好;在高颗粒质量分数下(约30%),不同湍流模型预测的相对比冲损失差别可超过3%。在3D RANS模拟中,发现Realizable k-ε模型表现出最佳的综合性能,所预测的物理量沿中轴线分布的准确性均优于其他模型。当颗粒质量分数约为31.2%时,Realizable k-ε模型所预测的比冲为与LES结果相比误差仅1.56%。

COPRA双层熔池湍流模型优选

为了预测下封头内双层熔池的流动和传热过程,基于不同湍流模型,同时采用凝固熔化模型对堆芯熔池研究装置(corium pool research apparatus,COPRA)双层熔池实验进行计算流体力学(CFD)数值模拟,通过数值计算获得准稳态下熔池的温度、沿壁面的热流密度与内壁面壳层的分布,将模拟结果与实验值进行比较,评价不同湍流模型的适用性和准确性,并进行湍流模型优选。结果表明,壁面模化大涡模拟(WMLES)湍流模型对下封头内双层熔融池流动与传热模拟的准确性和适用性最好;基于WMLES湍流模型,氧化层温度随着熔池高度增大而增大,氧化层上部存在强烈的湍流,在熔池底部的壳层最厚。

SAS湍流模型在瓦斯气云爆炸数值模拟问题上的应用可行性研究

SAS模型是一种新兴的CFD湍流模型。为了探究SAS模型在瓦斯气云爆炸问题上的可行性,对预混气云进行点火,与标准k-ε模型的模拟结果进行对比。结果发现SAS模型所得的模拟结果在气云爆炸发展过程、正负超压值和正负超压到达时间均与标准k-ε模型相近,验证了SAS模型的可靠性。本研究为SAS模型在工程上的应用提供一定的参考。

圆柱绕流的三大类型湍流模型及稳定性数值模拟研究

海港桩基码头下方桩柱周围的绕流特性对港口水动力与泥沙输运具有重大影响。为探究湍流数值模拟方法对三维桩柱绕流精度的影响,基于CFD方法建立OpenFOAM湍流模型,在雷诺数Re=3900条件下,研究不同近壁取值范围y+=0.95~5.00内的大涡模拟LES、原始雷诺时均法(RANS)湍流模型和增加涡流黏度限制公式的RANS湍流模型的阻力系数、斯特劳哈尔数、平均压力系数、平均速度和涡流结构。研究结果表明,与前人RANS各湍流模型及试验计算结果比较,RNGκ-ε湍流模型对阻力系数和斯特劳哈尔数模拟较准确,对涡捕捉精细程度最好;κ-ω、κ-ωSST湍流模型对斯特劳哈尔数模拟较准确,在y+=0.95~2.00时对平均压力系数模拟较准确;各RANS湍流模型对平均速度的模拟均存在回流区内速度偏差较大的问题;大涡模拟LES湍流模型在近壁网格精度较高时,所有水动力模拟结果相对准确。研究还发现增加涡流黏度限制后的RANS湍流模型,对涡结构捕捉的精细程度更高,更符合真实流动。通过对不同湍流模型下单桩水动力特性的分析研究,可为海港码头下圆形桩群的水动力计算提供参考。

基于修正参数k-l湍流模型的复杂地形流动研究

针对风资源精细化评估需求,提出采用实际测风数据修正k-l湍流模型的关键参数,以提高数值仿真对大尺度复杂地形流动的准确性和适用性。在此过程中,首先介绍k-l模型的求解过程并开展网格分辨率研究;其次利用实际测风数据修正模型参数并使用典型案例进行验证;最后以中国华南地区某实际风电场为对象,将其与k-ε模型进行对比分析以验证其有效性和优越性。结果表明,相比于k-ε模型,原始k-l模型由于高估了尾流影响,导致预测误差偏大;而修正参数后的k-l模型能有效提升精度,风加速因子和发电量等效小时数的平均相对误差分别降低5.1%和7%。因此,在k-l模型的实际应用过程中,需利用实际测风数据修正参数,或将参数B1的量级调整为102。

针对高速列车受电弓气动特性仿真的不同湍流模型适用性研究

为探究高速受电弓气动特性仿真时不同湍流模型在不同运行环境中的适用性,基于k-ε,SST k-ω和Spalart-Allmaras 3个湍流模型的理论基础和计算流体力学仿真方法,通过线路试验对仿真结果进行验证后,对比不同模型在开阔环境和隧道环境中的精确度和计算效率。结果表明:在300 km·h^(-1)速度下,SST k-ω与k-ε模型在开阔环境中的闭口和开口方向气动抬升力计算结果差异分别为0.7%和1.5%,与Spalart-Allmaras模型计算结果差异分别为7.1%和5.5%;SST k-ω和k-ε模型在边界层预测均匀,Spalart-Allmaras模型在边界层分离区的预测准确性较低,但在隧道环境中它的精确度高于k-ε模型;在计算效率方面,Spalart-Allmaras模型的迭代次数和收敛速度优于SST k-ω模型;在适用性方面,开阔环境中k-ε模型的精确度和计算效率表现优异,更为适用,而在隧道环境中Spalart-Allmaras模型具有较高的精确度和计算效率,更为适用,但在对精确度要求高且可忽略计算效率因素的情况下SST k-ω模型具有较好适用性和高精确度,更为适用。

工业离心压缩机数值模拟方法及湍流模型研究

为建立一套能够快速、准确预测工业离心压缩机基本级性能的数值模拟方法,针对某离心压缩机基本级试验装置开展了数值模拟方法的研究。探讨了数值模拟所采用的模型简化、网格型式、网格密度、边界条件、离散方法、湍流模型等,重点讨论了目前广泛使用的k-ε、RNG k-ε、k-ω、SST、SST-RM 5种湍流模型的特点和适用性。采用数值模拟方法对基本级试验装置在不同机器马赫数下的气动性能进行了预测,并与试验结果进行对比,发现数值模拟方法可以准确预测大部分工况,仅在近失速和近阻塞区与试验结果存在较大偏差。其中,采用SST-RM湍流模型的预测结果与试验结果最为吻合,这得益于其对流动分离的起始位置和分离区域的预测较其他湍流模型更为准确。研究结果对提高工业离心压缩机数值模拟精度有借鉴意义。

大气湍流模型与大气信道模型的研究与展望

为了克服大气湍流效应对无线光通信系统性能的影响,必须建立能准确描述大气折射率起伏的湍流模型。总结了国内外学者在大气湍流模型方面的研究进展,对不同大气湍流模型的理论建立进行深入分析并对比了其适用范围,同时介绍了在实际大气湍流中实验测量所建立的大气信道模型,最后对大气湍流模型的发展趋势进行了展望。