三维槽道湍流的大涡模拟

利用大涡模拟方法模拟了雷诺数为13000的充分发展槽道湍流流动。为了比较网格和亚格子模型对计算结果的影响,采用了三种不同的网格数目,以及Smagorinsky模型和动力模型等两种亚格子模型,并与直接数值模拟结果进行了比较。结果表明,大涡模拟可以比较准确地给出速度分布、脉动速度均方根分布等湍流统计参数。在此基础上,对近壁区小尺度湍流结构进行了研究。     

高雷诺数槽道湍流的壁面模化大涡模拟研究

选取基准壁湍流的槽道流动,研究了多种模型的壁面模化大涡模拟.模型包括经典的大涡模拟、Spalart-Allmaras、分离涡模拟和一种动态混合模型.基于摩擦速度的雷诺数范围为395~12000,采用3组粗糙网格,流向和展向维数分别同取37,49和65,法向维数保证y+(1)~1.主要研究平均速度、雷诺切应力分布、详细分析了各模型的特性差异并展示了相应的湍流结构.研究表明:在高雷诺数粗糙网格下,大涡模拟失去求解精度,分离涡模拟出现对数律不匹配,动态混合模型的计算接近直接数值模拟,其对数率区可解应力约占雷诺切应力的93%,边界层外层可解应力约占99%.这说明合适的混合模型可以在经济成本下保证计算精度,具有解决实际问题的潜力.      

大涡PIV方法测量水平矩形槽道内湍流耗散率

湍能耗散率的准确测量对工程实际和湍流的理论研究都有着重要的意义。以往的研究多是基于单点速度测量来估算单点的湍能耗散率或者有限体积内的平均值，而不能提供整个流场范围内的耗散率的分布情况。PIV能够测量瞬时流场的速度分布情况，它更适合于用来测量整个流场内的耗散率分布。本文分析了PIV测量得到的速度与单点测量方法之间的关系，提出了大涡PIV方法，运用大涡模拟中的亚网格尺度(SCS)应力模型，对PIV测量结果进行处理得到了水平槽道内的耗散率的分布。对采用不同SCS模型的结果以及量纲分析的结果进行了比较。     

大涡模拟中亚格子模型的改进及其在槽道湍流中的应用

在用大涡模拟的方法计算具有强剪切的槽道湍流时, 常用的亚格子应力模型, 包括考虑壁面修正的模型, 都会出现平均流剖面偏高的现象.这表明在壁面附近亚格子应力模型还不能描述实际情况.针对这一问题, 修正了亚格子雷诺应力模型的壁面函数, 得到了较好的计算结果.用修正后的模型计算出的平均速度分布、均方根速度的分布以及雷诺应力的分布, 均与直接数值模拟(DNS)的结果吻合较好.      

湍流大涡数值模拟的现状与展望

近年来,随着计算机技术和计算流体力学(CFD)技术的发展,大涡数值模拟(LES)越来越多地应用于解决湍流流动问题.本文介绍了当前LES研究与应用的主要领域以及取得的成果,分析了LES的局限及实施的难点,并对LES未来的研究方向和重点进行了总结.     

旋转槽道湍流的格子Boltzmann方法模拟

基于多松弛格子Boltzmann方法的大涡模拟对雷诺数为194,旋转数从0~5.0的旋转槽道湍流进行数值模拟.采用动态亚格子应力模型模化滤波后的不封闭项,修正二阶矩作用力模型计算压力梯度、哥氏力.对平均速度、均方根脉动速度、雷诺应力以及湍流结构进行了计算.结果显示哥氏力使流场平均速度呈现不对称性:在压力面,随着旋转数的增加,湍流度增强;而在吸力面湍流脉动减弱,具有层流化的趋势.将格子Boltzmann模型与直接数值模拟求解进行对比,结果验证了格子Boltzmann方法在旋转湍流模拟中的可行性.      

一种剪切修正的涡黏性亚格子模型及其在槽道湍流中的应用

提出了一种利用壁面附近强剪切对涡黏系数进行修正的亚格子模型,称为剪切修正的涡黏模型.该模型能使涡黏系数在壁面附近自动满足壁面修正,在脉动为零的地方,涡黏系数自动为零,保证了壁面附近涡黏系数趋于零的特性.用该模型对槽道湍流壁面无吹吸、有吹吸,较低和较高雷诺数的情况进行了数值模拟,与直接数值模拟的结果进行了对比,结果吻合较好.      

湍流大涡数值模拟进展

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

直方管内充分发展湍流的大涡模拟 第二部分

本文第一部分用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行了数值模拟。所得结果证实大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。在本文第二部分用所得的数据库对湍流运动的结构和形态进行了分析，并且用流体质点运动图形显示了湍流运动的二次流。     

直方管内充分发展湍流的大涡模拟 第一部分

本文用大涡模拟的方法对直方管内Ｒｅ＝４．９×１０￣４的充分发展湍流运动进行数值模拟，所得结果清楚地显示了二次流的存在。统计平均量的数值结果也与实验结果吻合甚好，这表明大涡模拟是高雷诺数复杂湍流运动数值模拟的一个有效手段。     

用于飞行仿真的二维紊流场的数字生成法

 在飞行仿真的某些特殊情况(例如编队飞行、空中加油情况)下有必要生成二维紊流场,其中的紊流速度是两个坐标(例如沿飞行方向的x和沿翼展方向的y)的随机函数。在已有的一维紊流场生成技术的基础上,本文提出了一种简易而有效的生成二维紊流场,即在矩形网络上生成紊流速度序列的方法。经抽样检验,这样生成的紊流场的相关函数基本上符合紊流模型的相关特性,因而证实了该方法的可行性。还讨论了在飞行仿真中二维紊流生成器的两种运行方式——在线式和离线式。     

用蒙特卡罗法仿真生成三维空间大气紊流场

基于三维紊流场相关函数矩阵用蒙特卡罗法仿真生成三维空间大气紊流场,并用理论证明和紊流场数值检验的方法,证明这样的空间大气紊流场较好地符合所选用的大气紊流模型 (如 Dryden模型),且具有良好的均匀性和各向同性。三维空间大气紊流场可用于飞机飞行仿真,尤其是多架飞机在大气紊流中飞行的仿真 (如空中加油、编队飞行)。     

混合流场控制的大涡模拟

采用大涡模拟方法研究了三维小扰动下,混合流场大尺度拟序结构的产生和演化过程,捕捉了展向涡的卷起、配对、合并,以及二次流向涡的出现等大尺度的三维拟序结构,分析了拟序结构与入口扰动方式之间的内在联系,再现了涡卷自身撕裂而引发转捩的现象。     

标准的与非线性的k-ε模型在二维剪切流中的应用

 采用标准 k-ε和非线性 k-ε模型对二维直通道及小曲率弯道中的剪切流进行了数值计算。与实验结果对比后得到 :由于非线性 k-ε模型考虑了流体变形与应力关系中的非线性项,因而能正确地反映湍流各向异性,能正确预测二维剪切流中湍流强度及湍流动能等参数的分布     

两种雷诺应力湍流模式的比较及其在涡流计算中的应用

本文将要讨论的两种雷诺应力湍流模式为:(1)雷诺应力代数模式;(2)雷诺应力输运模式。这两种模式分别应用于弱强涡射流的计算中,其结果表明,代数模式并不能很好地近似输运模式,特别是在一些重要的流场区域,如对称轴附近,前者给出非真实的物理结果。本文将探讨代数模式的这一缺陷,同时提出一新的输运模式。     

高湍流度畸变流场的模拟

应用“双板结构”模拟板成功地模拟了高湍流度畸变流场。模拟结果表明:板间距离增加将使板后流场平均总压恢复系数下降;第二块板有改变稳态畸变值和调整畸变图谱的作用。由此得到“双板结构”模拟板模拟高湍流度畸变流场的设计原则:第~块板应以宽板为主,以产生高湍流度畸变流场;第二块板起调整图谱,改变畸变值的作用。实验还表明“双板结构”模拟板畸变特性随Ma数的变化与“单板结构”明显不同。     

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

固体发动机喷管内两相流的有限元数值模拟

本文给出了固体发动机喷管内两相流的数值模似。从非守恒型两相流甚本方程出发,采用加权积分技术,导出了离散点的流动物理量的有限元数值计算公式。数值模拟分别与JPL喷管的两相流壁面静压分布测量结果以及与差分法计算结果作了比较,比较结果表明:二者是吻合的。本文还以中型喷管为例,计算了粒子的轨迹,粒子的速度滞后和燃气的温度滞后,计算结果表明;入口粒子冲击角对粒子轨迹有明显的影响。