室内流场湍流模型对比研究

为探明不同湍流模型求解室内流场的适用性,基于计算流体力学软件Fluent和Airpak的零方程模型、室内零方程模型、标准k-ε模型、RNG k-ε模型和Realizable k-ε模型,对暖通空调领域三种典型室内流场(无障碍物时的等温强制对流、有障碍物时的等温强制对流、有内热源时的混合对流)开展数值模拟,并将数值模拟得到的风速分布与实测值进行对比分析。结果表明:基于Fluent的Realizable k-ε模型和基于Airpak的RNG k-ε模型能够准确预测三种典型室内流场的速度分布;室内零方程模型对于室内无障碍物时的等温强制对流预测结果较差,而对于有障碍物时的等温强制对流和有内热源时的混合对流,在计算精度要求不是特别高的工程计算中可以使用;零方程模型对三种典型室内流场的预测结果均较差,在室内流场特性研究中应避免使用。

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

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

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

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

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

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

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

高功率空间核电源采用棒束气冷堆具有质量轻等优势,其紧凑的栅格结构以及流动雷诺数(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)的影响更显著。

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

研究固体火箭发动机喷管中颗粒质量流率对两相流的影响规律以及不同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%。

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

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

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

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

COPRA双层熔池湍流模型优选

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

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

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

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

为探究高速受电弓气动特性仿真时不同湍流模型在不同运行环境中的适用性,基于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-l湍流模型的复杂地形流动研究

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

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

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

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

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

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

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

基于湍流模型的DC-RF混合等离子体流动及传热特性研究

直流-射频(Direct Current-Radio Frequency,DC-RF)混合等离子体具有高温、高化学活性等特点,在核用超细粉末材料制备领域有着广阔的应用前景。对这种混合等离子体的特性进行研究,可以为等离子体发生器装置的设计及稳定运行提供参考。本文采用k-ε湍流模型,对DC-RF混合等离子体发生器内的热等离子体的温度及流场空间分布进行模拟,并在此基础上分析了各工作参数对混合等离子体流动及传热特性的影响。模拟结果表明:提高DC弧电流、增大反应气及冷却气流量均能减弱混合等离子体发生器入口处的回流效应;提升RF线圈电流则能增加RF线圈附近的等离子体温度及高温弧区面积,但过大的气流量及过高的线圈电流会对发生器装置的正常工作产生不利影响。因此,需要合理调节各种工作参数来改变混合等离子体的流场分布,从而在保证装置稳定运行的前提下满足不同材料处理的需求。

TPDF-ASOM复合湍流燃烧模型及其检验

先进航空发动机燃烧室设计要求对湍流火焰精确控制,现有模拟方法需提高精度和效率。输运概率密度函数(TPDF)湍流燃烧模型精度高,代数二阶矩(ASOM)湍流燃烧模型计算成本低,类比离散涡模拟思想,基于Da数将湍流燃烧场区分“高精度”和“低成本”2个区域,在输运方程框架下采用随机场TPDF(高精度)和ASOM(低成本)方法重构TPDF-ASOM复合湍流燃烧模型,以提高模拟的整体精度和效率。在大涡模拟(LES)-TPDF程序平台创建ASOM并进一步实现TPDF-ASOM复合湍流燃烧模型,用Flame D实验数据检验所建模型和方法。结果表明:所建模型的预测结果与实验值接近,而且能够兼顾精度和计算效率。

基于时空最优低维动力系统的多尺度可压缩湍流数值模拟方法

按照周培源教授关于研究湍流数值模拟建模时必须分析和求解脉动速度场的思想,该研究基于第一性原理,系统地建立了基于时空低维最优动力系统的多尺度可压缩湍流数值模拟方法(LMS方法),并将其应用于多次冲击Richtmyer-Meshkov问题的数值模拟中,首次得到了可压缩湍流的中尺度流场和不同于DNS近似解的湍流近似解.数值结果表明,LMS方法可以用较少的网格获得更精确的湍流近似解.首先解决了研究中遇到的几个问题,为LMS方法的构建铺平了道路.这些问题是:基于湍流的物理特性,提出了湍流大、中、小尺度分解的新概念;找到了box滤波空间相关性的计算方法;指出了湍流建模理论中长期存在的逻辑错误,提出了多尺度湍流模型的概念;讨论了湍流封闭问题的本质和关键,给出了克服湍流封闭问题的数值方法.采用box滤波方法/空间网格平均方法且在大尺度网格的意义下,LMS方法的本质是一种将RANS、LES、DES和DNS等湍流数值模拟方法统一的全新湍流数值模拟方法.需要指出的是,LMS方法也可以作为湍流模型研究的辅助工具,以检验SGS尺度方程/脉动方程中各项所对应的湍流模型是否正确.