用非线性k-ε-kp两相湍流模型[1*9/9]模拟旋流两相流动

建立了两相湍流的代数应力模型,并由此出发, 导出非线性k-ε-kp两相湍流模型, 目的是合理地模拟各向异性较强的旋流两相流动, 保持二阶矩模型的优点, 同时比二阶矩模型简单. 文中得到了气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式.这些代数式和两相各自的湍动能k, kp, 以及两相脉动关联湍动能kpg的方程联立, 就构成非线性k-ε-kp模型. 将该模型用于模拟旋流两相流动, 给出两相时均速度场及雷诺应力各分量,并且将模拟结果和实验数据以及二阶矩模型的模拟结果对照. 研究结果表明,该模型预报旋流两相流动的能力和二阶矩模型的能力相差不多,但计算量比二阶矩模型的小.     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

两相湍流的非线性k-ε-kp模型

由两相湍流的代数应力模型，提出一种两相湍流的非线性k-ε-kp模型，可以合理地模拟各向异性较强的旋流，浮力流等两相流动，具备二阶矩模型的优点，又比二阶矩模型简单。文中推导出气相、颗粒相的雷诺应力和两相脉动速度关联的非线性应力应变关系式。这些代数式和两相各自的湍动能k，kp以及两相脉动关联湍动能kpg的方程联立，就构成非线性k-ε-kp模型。本文将该模型用于模拟突扩湍流两相流动，给出两相时均速度场及雷诺应力各分量，并且将模拟结果和PDPA实验数据以及二阶矩模型的模拟结果对照。结果表明，该模型预报各向异性两相湍流的能力和二阶矩模型的能力接近，但是计算量比二阶矩模型的小得多。     

用非线性κ—ε—κp两相湍流模型模拟旋流两相流动

建立了两相湍流的代数应力模型，并由此出发，导出非线性κ-ε-κp两相湍流模型，目的是合理地模拟各向异性较强的旋流两相流动，保持二阶矩模型的优点，同时比二阶矩模型简单，文中得到了气相、颗粒相的雷诺应力和两相脉动速度关联的性应力应变关系式。这些代数式和两相各自的湍动能κ，κp，以及两相脉动关联湍动以κpg的方程联立，就构成非线性κ-ε-κp模型。将该模型用于模拟旋流两相流动，给出两相时均速度声及雷诺应力各分量，并且将模拟结果和实验数据以及二阶矩模型的模拟结果对照，研究结果表明，该模型预报旋流两相流动的能力和二阶矩模型的能力相差不多，但计算量比二阶矩模型的小。     

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

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

进口速度和含气率对气泡-液体两相湍流影响的数值研究

气泡-液体两相湍流规律的一个重要问题是,气泡究竟是加强还是削弱液体湍流,或者何时加强何时削弱液体湍流,尚有待深入探讨.本文用作者们提出的二阶矩液体-气泡两相湍流模型即两相雷诺应力方程模型模拟了二维通道内气泡-液体闭式多股射流在不同含气率和不同进口速度下的气泡和液体湍流脉动速度.预报结果和文献中给出的实验结果趋势一致.研究结果显示出低含气率及低进口速度下气泡增强液体湍流,高含气率及高进口速度下气泡削弱液体湍流的规律,澄清了众说纷纭的看法.     

SSG模型在压缩拐角激波/湍流边界层干扰模拟中的应用

采用非线性方案模化雷诺应力再分配项的二阶矩雷诺应力输运模型——SSG模型和两方程线性涡粘模型SST模型,对24°高雷诺数二维压缩拐角激波/湍流边界层相互干扰流动进行数值模拟。通过在壁面压力分布、壁面摩擦阻力系数分布、截面速度分布及分离区大小等方面与实验数据比较,综合评估了上述两种模型模拟该问题的能力。结果表明:两种湍流模型在上述各方面都取得了与实验较为一致的结果,但相比较而言,除SST模型在流动再附后对壁面摩擦阻力系数分布的模拟更接近实验值外,SSG模型结果在上述其他方面均更接近实验值。本文工作为压缩拐角激波/湍流边界层干扰这类流动数值模拟的湍流模型选择提供了参考。     

新型单方程湍流模型构造及其应用

为了提高雷诺应力本构关系式对于非平衡湍流的预测精度并且兼顾求解效率,发展了一种基于湍动能k的单方程(KDO)湍流模型。其主要思路为:采用平板直接数值模拟(DNS)数据对原始Bradshaw假设进行重新标定,使得当地湍动能和雷诺主应力之比能够根据当地流动条件进行自适应调节;同时,对标准k-ε模型中的湍流耗散率输运方程采用代数形式进行模化,进而形成一种一方程湍流模型。算例结果表明:KDO湍流模型对于对数率能够准确反馈,而在带有激波或部件干扰等流动现象的RAE-2822、ONERA-M6和DLR-F6算例中,KDO湍流模型能够准确控制湍动能的增长和衰减,相比于Spalart-Allmaras和Menter k-ω剪切应力输运(SST)模型,KDO湍流模型的计算结果有了较为明显的改善。     

三维边界层的湍流特性及其二阶矩模拟

根据现有实验数据,本文进一步分析了三维边界层内纵向涡对雷诺应力分布的影响,从而建立了一个足以模拟三维边界层湍流特性的二阶矩封闭模式。在新模式中,着重考虑了三维平均涡场引起的雷诺应力输运的非线性作用。以文[3]给出的模拟无限后掠翼实验为例,本文分别应用基准的二阶矩模式和新建立的非线性二阶矩模式对三维边界层内的雷诺应力分布作了数值预计。结果表明:应用基准二阶矩模式所得雷诺应力诸分量的预计值均比相应的测量值大得多,且完全不能预计切应力方向相对于速度梯度方向的滞后;而非线性二阶矩模式对雷诺应力分布的预计则基本上与测量数据符合,且能定性地模拟切应力方向相对于速度梯度方向的滞后。     

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

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

用于可压缩自由剪切流动的湍流混合长度

抓住可压缩流动变密度特性,构造出基于有效涡量的三维von Karman混合长度。湍流模型采用仅依赖湍动能k的单方程KDO(Kinetic Dependent Only)模型,引入新构造的混合长度替换旧尺度得到CKDO模型。为了验证其描述可压缩自由剪切湍流的能力,选择无壁面束缚、密度梯度大和可压缩效应强的自由剪切混合层为算例,其对流马赫数Ma_c=0.8。计算结果表明,KDO模型对混合层的速度分布有着良好的控制和模拟,而经可压缩修正后的CKDO模型与原模型及其他可压缩修正模型相比,所计算的速度分布、主雷诺剪切力和混合层厚度与试验结果更加接近,说明了该混合长度对可压缩混合层这种自由剪切湍流有着良好的刻画能力。     

改变结构参数对旋流喷嘴内外两相流场的影响

针对一种超细粉体旋流分散喷嘴,在不同进气角度和出口锥角的条件下,采用雷诺应力湍流模型(RSM)并考虑气体的可压缩性,对喷嘴内部及出口附近的强旋、跨声速流场进行了模拟,同时结合离散相模型研究了1~10μm粒子的运动轨迹和质量浓度分布规律.结果显示:改变旋流喷嘴的进气角度对旋流强度和湍流强度影响不大,而出口锥角改变时,喷嘴内的旋流强度有明显变化;随着颗粒密度或出口锥角增大,或进气角度减小,旋流喷嘴内粒子的临界逃逸粒径相应减小;在喷嘴后方,颗粒相达到较好的扩散均匀性所需的掺混距离约为20倍喷嘴出口直径.      

两相圆湍射流近场动态特性的大涡模拟(英文)

采用大涡模拟方法研究了负载颗粒两相圆湍射流的近场动态特性。为了揭示颗粒对湍流的调制规律,数值模拟中采用了点力双向耦合模型。脉动速度统计结果展现了颗粒对流动结构调制导致的相关统计特性的变化。在两相射流中,颗粒主要激起了射流初始的不稳定性并引起近场涡环的提前破碎。轴向和径向流动脉动强度因此被增强。文章还给出了颗粒对近场平均速度调制的结果。该研究对于两相湍射流的应用和控制具有重要参考价值。     

高超声速尾迹雷达散射湍流场研究

为给高超声速再入尾迹亚密湍流雷达散射分析提供所需的脉动背景场参数 ,本文提出计算非平衡再入湍流尾迹脉动等离子体场的理论方法。在研究高超声速尾迹流动特征的基础上 ,推导、使用包括化学组份浓度脉动强度的k ε g湍流模型 ,用以封闭高超声速尾迹雷诺平均控制方程 ,并用全隐式有限差分求解。以M∞ =2 1 .2 6、Re∞D =1 .33× 1 0 6的小钝锥体流动为例 ,得到的结果说明 :流场的流向和径向参数分布合理 ;在转捩点后较近距离内湍流脉动影响较大 ,随着向下游流动脉动影响迅速减弱 ;本文计算尾迹湍流脉动等离子体场的方法是可行的     

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

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

一种新的非线性K-ε两方程湍流模型

本文基于理论分析和对标准Ｋ－ε两方程湍流模型数值模拟结果的分析,在标准Ｋ－ε两方程模型的雷诺应力表达式中加入６高阶非线性项;并在模型的涡粘性系数μ１中考虑流动各向异性的影响,将涡粘性系数中参数Ｃμ取为表征流动各向异性的参数Ａ的一个线性函数,从而得到了一个新的非线性Ｋ－ε两方程湍流模型。     

用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

旋流燃烧室内多组分气体湍流混合的数值模拟

为合理地描述旋流条件下组分质量的湍流输运,将考虑湍流-旋流相互作用的新代数质量流通量模型与新代数Reynolds应力模型结合,构成了新代数应力/质量流通量模型。应用该模型对旋流燃烧室内两股不同组分射流的相互作用与湍流混合进行了数值模拟,得到了与实验相符合的气体轴向与切向速度、轴向与切向脉动速度均方根值和氦气浓度分布。模拟结果与k-ε模型的结果进行了对比。     

预混燃烧大涡模拟和燃烧模型的检验

用代数二阶矩亚网格(ASOM-SGS)燃烧模型对文献中测量的钝体后方丙烷-空气预混燃烧进行了大涡模拟,模拟统计的时平均速度、速度脉动均方根值和温度分布与实验数据结果吻合很好,表明所采用的ASOM-SGS亚网格燃烧模型是合理的。模拟的瞬态结果显现了钝体后方湍流流动和火焰结构。将大涡模拟数据统计得到的反应率系数-浓度关联量的分布规律,与代数二阶矩RANS(ASOM-RANS)燃烧模型的模拟值进行对比,结果发现,大涡模拟统计值和ASOM-RANS模型的模拟值很接近,从而证明了湍流燃烧代数二阶矩RANS模型的合理性。     

湍流燃烧模型对双旋流燃烧室喷雾燃烧的影响

采用数值模拟与试验测量相结合的方法,研究扩展旋涡破碎模型、扩展二阶矩模型和涡团耗散概念模型等三种湍流燃烧模型对双旋流湍流喷雾燃烧流场的影响.在任意曲线坐标系下数值研究双级轴向旋流器环形燃烧室全流程流场,采用粒子图像测速仪测量燃烧流场气流速度分布,热电偶测量燃烧室出口温度分布.计算结果与验证试验数据比较表明:不同湍流燃烧模型对双旋流湍流喷雾燃烧影响较大,所得的回流区形状、速度、温度场以及出口温度分布等都不太相同,其中扩展二阶矩模型所得的结果与试验值符合最好,更适用于模拟双旋流环形燃烧室湍流喷雾燃烧.