气粒流动两相亚网格尺度动能模型的大涡模拟

目前多数两相流动的大涡模拟采用单相流动的亚网格尺度应力模型,没有考虑或者没有完整地考虑两相亚网格尺度应力间的相互作用。本文提出一种两相亚网格尺度动能模型(LES-kkp),可以完整地考虑其相互作用。将该模型用于突扩气粒两相流动的大涡模拟,同时进行统一二阶矩(USM)两相湍流模型的雷诺平均模拟(RANS)。瞬态结果显示出两相流动中的气体形成了典型的拟序结构,但颗粒流动结构完全不同于气体的,没有形成涡结构。LES-kkp和RANS-USM模拟结果实验检验表明,对时间平均速度两种模拟给出结果与实验结果吻合良好,即RANS-USM得到了LES-kkp的验证。对两相均方根脉动速度和两相脉动速度关联量,LES-k-kp的模拟结果优于RANS-USM的模拟结果。     

旋流和无旋突扩流动的LES和RANS模拟

本文用smagorinsky-Lilly亚网格尺度湍流模型对旋流突扩流动(s=0.53)和无旋突扩流动(s=0)进行了大涡模拟(LES模拟),同时分别用压力应变项为IPCM和IPCM+Wall模型的雷诺应力方程模型进行了RANS模拟,和LES的统计结果对比。LES的统计结果与雷诺应力模型的模拟结果及实验对照表明,LES结果与实验结果的吻合比雷诺应力模型的好,说明所用的亚网格尺度湍流模型对旋流流动是适用的,LES结果是可信的。LES的瞬态结果揭示出在旋流作用下,流场中存在复杂的旋涡脱落现象。大涡结构极易破碎成小涡,而在无旋突扩流动的情况下,由于剪切的作用更强,大涡结构的尺寸和范围比旋流流动的要大。     

湍流边界层拟序结构的大涡模拟研究

采用动力亚格子模型,利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。从瞬时速度和脉动 速度场、脉动速度相关、均方根脉动涡量分布、以及瞬时涡量场等多个方面,对湍流边界层流动的拟序结构进行了分析, 包括近壁区小尺度湍流结构和瞬态过程,如条纹结构、喷射和扫掠过程、以及近壁旋涡结构等。     

一种改进后的亚格子特征尺度表达式

针对大涡模拟涡粘性亚格子模式中使用网格尺度为特征尺度存在的不足,提出一种改进的亚格子特征尺度表达式,利用结合Sagaut混合尺度模式获得的改进模式与Smagorinsky模式、Germano模式和Sagaut模式等进行对比研究.时间发展混合层和中性大气边界层的大涡模拟结果表明,改进的有效亚格子特征尺度能合理反映亚格子脉动的时空分布特性,同时扩展了适用范围.在对耗散程度、流动演化过程、可解湍流强度和可解雷诺应力等方面的模拟,改进模式优于基准模式,表明亚格子脉动有效特征尺度在大涡模拟中具有一定的理论基础和广泛的应用价值.     

基于转捩模型及声比拟方法的高精度圆柱分离涡/涡致噪声模拟

基于七阶加权紧致非线性格式(WCNS-E8T7),结合延迟分离涡模拟(DDES)和Ffowcs WilliamsHawkings声比拟方法,对亚临界雷诺数下单圆柱、圆柱-翼型的分离涡/涡致噪声问题进行了数值模拟.针对亚临界雷诺数下圆柱尾迹中的转捩问题,发展了基于γ-Reθ模型高精度转捩-延迟分离涡模拟(Tran-DDES)方法,并与传统基于全湍流剪切应力输运(SST)模型的高精度DDES方法进行了对比.单圆柱模拟结果表明:传统SST-DDES方法会造成平均流场的回流区增大,压差阻力偏小等问题;而添加转捩模型的Tran-DDES方法与实验符合得很好.圆柱尾迹中添加翼型后,翼型对圆柱附近流场产生影响,使SST-DDES方法造成的圆柱后回流区偏大的问题减弱,并与Tran-DDES模拟结果差异变小.但在脉动量预测以及脉动产生的噪声预测方面, Tran-DDES方法仍与实验符合得更好.     

圆形射流湍流场的大涡模拟研究

使用发展的大涡模拟并行程序,数值模拟了圆形射流湍流场,与实验值进行详细比较,对大涡模拟程序进行检验,并分析网格尺度对计算结果的影响。结果表明,大涡模拟计算的平均速度场和雷诺应力场与实验值符合得非常好,验证了该大涡模拟计算程序的可靠性和精确度。同时发现一阶统计矩受网格尺度影响不大,二阶及高阶矩与网格尺度有较大的关系,因此使用大涡模拟方法计算湍流及湍流燃烧问题应对网格依赖性进行必要分析,至少保证二阶统计矩具有较好的精度。     

非结构网格上大涡模拟数值方法研究

随着计算能力的飞速发展,在非结构网格上进行大涡模拟是研究叶轮机械真实几何结构下复杂流动问题的有效途径。本文针对非结构网格上数值格式人工黏性过大的缺点,利用湍流中动能级联特性,发展了适用于大涡模拟的低耗散数值格式,并结合亚格子模型,考察了其在基础湍流算例中的有效性。通过与传统迎风格式的对比,说明本文发展的低耗散格式可以显著提高非结构网格上大涡模拟的精度和可靠性。     

非结构网格上大涡模拟数值方法研究EI北大核心CSCD

随着计算能力的飞速发展,在非结构网格上进行大涡模拟是研究叶轮机械真实几何结构下复杂流动问题的有效途径。本文针对非结构网格上数值格式人工黏性过大的缺点,利用湍流中动能级联特性,发展了适用于大涡模拟的低耗散数值格式,并结合亚格子模型,考察了其在基础湍流算例中的有效性。通过与传统迎风格式的对比,说明本文发展的低耗散格式可以显著提高非结构网格上大涡模拟的精度和可靠性。     

颗粒增强湍流的新模型和气粒流动的湍流变动

用雷诺应力方程模型和极细的网格系对单个颗粒受湍流气体绕流进行了数值模拟,研究了改变颗粒直径和气体相对速度时颗粒增强气体湍流的规律.据此构造了颗粒尾涡增强气体湍流的新模型.将此子模型加入到两相流动模型中,对竖直和水平通道内气粒两相流动进行了数值模拟,和实验结果的对照表明,考虑颗粒尾涡增强气体湍流效应得到的气体湍流脉动速度的模拟结果比不考虑此效应的模拟结果好得多.     

均匀入口鼓泡流化床的TFM及DEM模拟

分别运用双流体模型（TFM）及离散元模型（DEM）的软球模型，模拟了二维及三维实验尺度的均匀入口鼓泡流化床内的气固两相流动特性，其中TFM模型结合了气相κ-ε湍流模型，而DEM模型则结合气相κ-ε湍流模型或Smagorinsky亚网格涡黏模型（SGS）。通过对比本文模拟结果与他人实验及计算结果发现，TFM与DEM软球模...     

非预混燃烧大涡模拟的代数二阶矩模型

一种改进的代数二阶矩燃烧模型被应用在Flame-D的大涡模拟中。代数二阶矩模型能够将化学反应速率的脉动项进行模拟,更准确地计算化学反应速率。在实际燃烧过程中,化学反应速率不仅与当地的温度以及组分浓度等相关,更与反应物的混合速率有密切的关系。研究考虑了亚网格尺度内混合速率对化学反应速率的影响,提出了改进的代数二阶矩模型。通过与相同模拟条件下的涡耗散概念模型的模拟结果以及实验结果进行对比,改进的代数二阶矩模型的准确性得到了验证,同时亚网格化学反应速率的模拟对总化学反应速率的影响也得到了研究。     

不可压反应流动直接模拟和亚网格模型的检验

本文用谱方法对三维槽道不可压湍流反应流动进行了直接模拟,用直接模拟数据对大涡模拟亚网格质量流和燃烧模型进行了检验,结果发现,引入壁面阻尼修正的模型与精确值的符合比较好.     

两种亚网格燃烧模型的对比研究

为详细比较动态增厚火焰模型和火焰锋面密度模型的性能,本文分别采用这两种亚网格燃烧模型对工业燃气轮机PRECCINSTA的模型燃烧室内湍流预混火焰进行了大涡模拟研究。计算结果与实验数据吻合良好,且这两种模型预测得到的速度、温度和主要组分的统计数据也非常接近。然而对于CO的分布两者的预测结果差异明显,文中对此进行了简要分析。研究还表明本文的计算结果与实验数据的吻合度与文献中报道的前期结果相比有所提高。     

基于大涡模拟的小型转子发动机缸内湍流研究

本文应用大涡模拟的方法对小型转子发动机的缸内流场进行了三维数值分析,分别从拟序结构、湍流脉动和湍动能三个方面研究小型转子发动机缸内的瞬态湍流特征。计算结果表明,大涡模拟结合Q准则可以有效的识别小型转子发动机缸内流场的大尺度拟序结构。进气过程中高强度的涡团主要分布于燃烧室的前部和后部,压缩过程中高强度的涡团分布从燃烧室中部向后部转移。湍流在垂直于转子截面上的法向脉动最为剧烈,其中中部区域的剧烈程度远高于两侧。湍流在平行于前后端面的截面上的法向脉动最弱,没有明显的波动。从进气口打开到压缩上止点期间,缸内亚网格湍动能和平均气流速度出现两次明显的峰值,这分别与进气能量和大尺度涡团的破碎有直接的关系。     

贯流风扇气动声场数值模拟及分析

针对贯流风扇气动噪声传播特性,采用基于非结构网格的CE/SE算法对其气动流场进行数值模拟,并模拟了气动噪声的传播.湍流模拟采用大涡模型,搭接式滑移网格模型处理动静干涉.噪声传播采用完全欧拉方程作为控制方程,远场采用无反射边界条件.将使用该方法得到的远场噪声频谱与实验数据进行对比分析,结果表明与实验测量得到宽带频谱不同,采用大涡模拟湍流模型进行声源模拟时,会加强部分离散频率上的声级,从而产生误差.     

空间发展湍流气粒两相平面混合层的大涡模拟与统计

本文对空间发展的湍流气固两相平面混合层流动进行了大涡模拟研究,其中气相亚网格尺度（SGS）使用结构函数模型,气相控制方程组采用SIMPLE方法求解,固体颗粒运动用拉格朗日方法计算。计算结果正确重现了流体涡结构的卷起、合并和破碎过程,以及小尺寸颗粒在涡边缘（低涡度区）的局部富集现象。对直径分别为42μm、72μm和135μm分别进行了模拟,并将统计结果和实验测量结果（Hishida　et　al[1]）比较,表明两者的平均速度吻合很好,但颗粒数密度和脉动速度存在较明显的差异,因此有必要对亚网格应力和颗粒之间的耦合作用以及拟序结构的三维性对颗粒运动的影响开展深入研究。     

气动光学效应湍流脉动模型系数修正

以光学窗口的气动光学效应研究为背景,研究影响光学传输预测的流场湍流脉动量预测.选取平板、压缩拐角、凹腔流动模型,采用大涡模拟（LES）方法研究超声速湍流脉动工程预测模型的系数修正.结果表明：LES方法能够获得无激波干扰、强激波干扰及底部大分离条件下,湍流密度脉动的定量分布,并据此给出脉动工程模型的系数修正,结果已经应用于型号飞行器的光学窗口气动光学效应预测.     

加力燃烧室热态流场的大涡模拟

本文采用k方程亚网格尺度模型对带V形槽稳定器模型加力燃烧室紊流化学反应流进行大涡模拟的研究,采用亚网格EBU燃烧模型估算化学反应速率,用热通量辐射模型估算辐射通量。数值计算表明,在稳定器后面开始在短时间内出现涡的交替脱落和逐渐消失过程,然后形成稳定的回流区,数值计算结果与实验比较吻合,表明可用大涡模拟研究实际燃烧室瞬态燃烧流场。     

尾流的双尺度简化模型及其大涡模拟建模

叶轮机中存在着复杂的湍流流动结构,如何准确地模拟湍流流动一直是工程界关注的核心问题。由于湍流在大部分区域都处于强非均衡态,根据均衡能量传输假设建立起来的RANS湍流模型和大涡模拟亚格子模型都不再适用。非均衡湍流的一个经典的模型流动是尾流,但通常意义上的尾流往往是无数不同尺度相互作用的结果,难以进行简化分析。本文构造了一个极度简化的双尺度尾流模型进行了直接数值模拟(DNS),结果显示在转捩前期就存在反向的非均衡能量传输。分别采用各种传统的RANS湍流模型和大涡模拟亚格子模型进行了数值模拟,发现各种模型都无法准确预测尾流转捩。在考虑非均衡能量传输机制后,通过流向速度梯度扭率构造的理性亚格子应力模型可以使大涡模拟模拟的转捩位置相对于标准Smagorinsky模型有所改善;不采用涡粘假设的改进速度增量(IVI)模型可以得到更加接近DNS的结果。这些结果为湍流模型进一步的改进提供了思路。     

湍流燃烧二阶矩模型的大涡模拟验证

对美国Sandia国家实验室测量的甲烷-空气湍流射流火焰进行了三维大涡模拟(LES),其统计的时均量以及脉动量符合实验测量值,且和RANS代数二阶矩湍流燃烧模型模拟结果大致符合。LES统计得到的化学反应率系数脉动和甲烷浓度脉动的互关联与对应的时均量梯度的乘积分布趋势相同,大小幅值的位置也相对应,与二阶矩湍流燃烧模型的封闭假设相符,验证了二阶矩代数燃烧模型的合理性。