不可压缩湍流的多尺度模型

基于可变时间间隔平均方法,提出了一种不可压缩湍流多尺度模型及平均流动方程.与传统的单尺度湍流模型不同,该模型在建立的过程中,保留了湍流的多尺度特性,结合平均流动方程,可以更好地预测湍流流场特征.通过模拟平面后台阶流动和不对称平板扩压器流动,并将预测结果与标准k-ε模型的预测结果对比,初步验证了其可信性及优越性.结果表明:计算所得的平面后台阶流动的流向再附长度与台阶边压力系数比比标准k-ε模型提高精度约20%;平板扩压器流动的回流区位置误差约为7%、倾斜壁面摩擦因数误差约为5%,而标准k-ε模型未能预测出分离现象.可以看出该模型适用于典型的分离流动,在湍流流场的预测中表现优异,具有一定的工程应用价值      

基于人工神经网络的可压缩湍流大涡模拟模型

在国家数值风洞（NNW）工程项目的指导下,空间人工神经网络（SANN）模型被用于强可压缩湍流大涡模拟（LES）研究,其中流场的湍流马赫数分别为0.6、0.8、1.0。基于湍流的多尺度空间结构特性和人工神经网络方法发展的高精度空间神经网络（SANN）模型适用于不可压缩湍流和弱可压缩湍流。对于强可压缩湍流,流场中会出现激波结构,给大涡模拟带来了挑战。本文的研究结果表明：SANN模型适用于强可压缩湍流的大涡模拟。在先验分析中,SANN模型预测的亚格子应力和亚格子热流的相关系数超过0.995,远远高于梯度模型和近似反卷积模型等传统模型;传统模型的相对误差大于30%,而SANN模型在这方面有很大的改进,相对误差低于11%。在后验分析中,与隐式大涡模拟（ILES）、动态Smagorinsky模型（DSM）、动态混合模型（DMM）相比,SANN模型能更精确地预测能谱、各类湍流统计特性以及瞬态流动结构。因此,基于湍流多尺度空间结构特性的人工神经网络模型加深了人们对强可压缩湍流亚格子建模的认识,同时可以服务于NNW工程的流体力学模型构造。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。     

可压缩湍流的多尺度分析

作者的研究团队近几年在可压缩湍流的多尺度性质方面开展了系统的研究工作.通过多过程分解,研究了可压缩湍流中的速度和热力学量的剪切部分、胀压部分、伪声模态、声模态、熵模态的多尺度性质,并总结了各类可压缩条件下的速度和热力学量的谱的标度律.通过滤波方法,研究了动能和热力学量的多尺度传输现象,并重点分析了可压缩性对动能胀压部分...     

不可压缩翼型绕流数值研究

使用标准k-ω模型及与色散模型相耦合的k-ω模型分别计算了NACA0012翼型和NACA 4412翼型的低速绕流问题.NACA 0012翼型计算了其来流雷诺数为2.88×106,攻角从0°到15°范围内的流动结构、翼型表面压力分布和升力、阻力特性;NACA 4412翼型计算了临界雷诺数为1.52×106,攻角为13.87°时的流动分离和翼型表面压力系数 ,并与实验数据进行对比.结果表明:在同等条件下,使用与色散模型相耦合的k-ω模型计算得到的NACA 0012翼型的升力和阻力系数比标准k-ω模型提高精度约5%,NACA 4412翼型的表面压力系数精度提高了约3%,进一步验证了其可信性,可将其进一步应用到低速飞行器的气动计算中.      

可压缩湍流边界层近壁区马蹄涡的演化

采用大涡模拟方法对可压缩槽道湍流进行了数值模拟.通过给定扰动波方法确定进口边界条件,很快得到了完全发展的湍流.利用动态亚格子模型、预处理技术以及高阶格式和方法,得到了精确的数值结果.计算得到的速度剖面和脉动分量与直接数值模拟的解进行了对比,验证了计算结果的可靠性.结果清楚地展示了湍流场中马蹄涡的生成和演化过程,包括对称涡腿、单侧涡腿的马蹄涡结构,以及流场中存在由亚谐波引起的拟序结构的交错现象.分析了湍流边界层近壁区可能出现的复杂的马蹄涡结构,及其涡结构形状也不严格对称的现象.     

SST湍流模型的激波与可压缩效应改进

由于缺乏对某些重要流动特征的考虑,针对不可压流发展的标准SST湍流模型在描述超声速流场时存在明显的局限性。为改善SST模型在吸气式高超声速推进系统内流等复杂超声速流场中的预测精度,基于流动特征结构定向开展了激波和可压缩效应改进。通过激波/湍流边界层判别函数和可压缩效应判别函数定位标准SST模型参数或建模假设失效的区域,针对性地改进湍流模型。采用超声速平板边界层流动、超声速压缩拐角分离流动、超声速隔离段复杂激波串流动以及HIFiRE-2超声速内流等算例进行了测试,结果表明改进模型具有与标准SST模型一致的边界层预测能力,但显著提高了对激波干扰流动及逆压分离流的预测精度。     

用湍流模型研究Richtmyer-Meshkov不稳定性诱导的湍流混合

将k-ε模型应用到存在激波的Richtmyer-Meshkov不稳定性诱导混合的研究中,采用代数关系式补充模型的封闭假设,给出了合理描述RMI特征的湍动能产生项的表达式。采用该模型对两个激波管实验进行了模拟。通过与实验和文献中的直接模拟结果比较分析,验证了所采用的模型封闭、模型常数、数值算法、程序实现是合适的。     

湍流涡黏性模型方程中相位差的测量

基于理论上湍流相干结构复涡黏性模型对涡黏性系数的分析,应用热线测速技术测量了平板湍流边界层多尺度相干结构动力学方程中非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差。分析了湍流边界层相干结构猝发过程中,非相干结构成分对相干结构贡献的雷诺应力分量与相干结构流向速度流向变形率之间的相位差沿湍流边界层法向的变化规律,为建立更加符合实际的湍流模型提供了实验依据。     

湍流输运的张量系数模型

实验表明,除均匀各向同性湍流之外,由湍流脉动引起的标量的输运都是各向异性的,且平均量梯度与湍流输运通量之间有一夹角。因此,在用梯度准则模拟湍流输运时,其输运系数应是一个二阶张量。本文推出一个张量形式的输运系数计算公式,并用它分析均匀各向同性湍流、有应变率的均匀剪切流、管内流动和二维平面射流中的热通量问题,均获得了满意的结果。     

湍流大涡数值模拟进展

本文简要陈述湍流大涡数值模拟的原理、优点,着重讨论湍流大涡数值模拟方法的关键问题及其可能解决的途径,包括脉动的过滤、亚格子模型、近壁模型和标量湍流的大涡数值模拟中的特殊问题.文章强调大涡数值模拟中亚格子应力的本质是可解尺度湍流和不可解尺度湍流动量间的输运,并以作者最近提出的新型亚格子模型说明发展亚格子模型的正确途径.文章最后提出湍流大涡数值模拟近期需要迫切解决的问题和其他具有挑战性的方向.     

壁面湍流模型对湍流分离流动数值模拟的影响

本文用三种近壁湍流模型计算了二个二维的湍流分离流动：单包流动和多包流动。结果表明：二层模型和低Reynolds数模型具有类似的特性。它们基本上给出合理的分离流动结构。壁函数的方法由于对数律的假定基本上不适合于计算湍流的分离流动。     

Gamma-Theta经验转捩模型在DES中的应用

与RANS湍流模型相比,DES湍流模型的边界层外自由流部分存在LES区,该区域的涡粘性会有所降低。为了研究了这一特征对γ-Reθt转捩模型的影响,以平板实验为算例,选取基于k-w SST湍流模型的DES模型和DDES模型,检验γ-Reθt转捩模型在DES模型和DDES模型中的表现。结果表明：DES模型和DDES模型使转捩位置推迟,SST—DES模型由于对网格过度敏感而没有触发转捩,SST—DDES模型使转捩位置后移较少;DDES模型比DES模型更适合γ-Reθt转捩模型;转捩经验关系需要进行修正后才能将γ-Reθt应用于DDES模型。     

湍流边界层各向异性张量模式的实验研究

使用高时间分辨率粒子图像测速技术在水槽内精细测量平板湍流边界层内不同流向、法向位置的瞬时速度矢量空间场的时间序列信号。利用实验数据对湍流各向异性的涡黏张量模式进行了实验研究,给出了湍流边界层不可压缩槽道湍流中各向异性的涡黏张量各分量的空间分布规律。各向异性的涡黏张量分量的空间分布揭示了流向平均速度的法向梯度对雷诺应力三个分量的不同贡献。涡黏张量三个分量空间上的不同规律说明对涡黏张量的各向异性的考虑是合理的,而各向同性的涡黏模型不适宜描述剪切湍流雷诺应力的物理本质。     

不可压湍流反应流直接模拟和RANS反应模型检验

用谱方法对三维不可压槽道湍流反应流动进行了直接数值模拟,得到了温度和质量分数耦合的瞬态数据库.结果显示温度与质量分数脉动在近壁区都有条带结构.推导了雷诺平均方法中湍流二阶矩反应模型中所求关联量的精确输运方程,发现耗散项是化学反应率系数-质量分数关联量封闭的关键,耗散项需要考虑化学反应影响.基于数据库的统计结果,对关联量模型方程中各项进行了先验研究,发现产生项和耗散项的贡献最大,扩散项和反应项的贡献较小,化学反应对各项大小和分布形状有明显影响.在算例中,直接模拟统计得到的化学反应率系数-质量分数关联和用代数二阶矩模型的模拟值很接近,说明ASOM模型具有一定的合理性.      

用非线性涡粘性模式计算三维湍流边界层

针对非线性涡粘性模式在求解三维湍流边界层流动时的不足，本文从压力－变形率关联项中的快速项出发，考虑流动非均匀影响，在显式代数应力模式的雷诺应力表达式中引入了反映雷诺应力“松驰”效应的速度二阶导数项，构造了一个新的非线性涡粘性模式。通过对典型算例－翼体角隔流动的计算结果表明，新模式能较好地再现出三维湍流边界层内雷诺切应力方向的发展滞后于速度梯度方向发展这一流动特性。     

一种改进的类DES湍流模拟方法

构造了一种基于一方程S-A(Spalart-Allmaras)模型和一方程Yoshizawa亚格子模型的混合RANS/LES(Reynolds-averaged Navier-Stokes/large eddy simulation)湍流模拟方法.在涡黏假设的基础上,将Yoshizawa亚格子湍动能方程转化为等效的亚格子湍流涡黏性输运方程,并采用混合函数将其与S-A模型方程进行混合,从而改进了DES(detached eddy simulation)模型的亚格子行为,同时克服了其依靠网格控制模型转换的缺点.模拟了超声速的带斜坡凹腔流动,并与相同网格下的LES及DES结果进行了比较,结果表明该混合RANS/LES方法在远离壁面的自由剪切流区域与LES行为一致,而在附着边界层区域表现优于LES和DES方法.      

用二阶矩两相湍流模型模拟鼓泡床内气泡—液体湍流两相流动

建立了二阶矩气液两相湍流模型，模拟了不同工况下二维矩形断面气液鼓泡床中气泡－液体湍流两相流动，给出了气泡和液体的速度场、气泡体积分数和两相雷诺应力分布，基本工况的模拟结果与美国俄亥俄州立大学PIV测量结果符合较好。文中研究了不同气体表观速度对两相流场的影响。模拟结果显示了鼓泡床内液体的回流流动和气泡的上升运动、各向异性的上湍流，气泡湍流脉动比液体的强，以及气泡体积分数和两相湍流强度随着气体表观速度的增大而增大等规律。