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

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

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

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

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

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

用非线性k-ε-k_p两相湍流模型模拟旋流两相流动

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

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

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

航空发动机燃烧室流动数值计算中湍流模型的比较

分别采用标准k-ε湍流模型、RNGk-ε湍流模型、Realizablek-ε湍流模型以及雷诺应力模型,对某型航空发动机燃烧室流动进行了数值计算.近壁处理采用标准壁面函数法,计算得到速度矢量分布以及质量流量、湍流粘度比和湍流强度等参数.四种湍流模型计算的总体流动差别较小,但流场的细节有较明显的不同.标准k-ε模型、Realizable k-ε模型和雷诺应力模型的湍流粘性计算结果较为接近,而RNGk-ε模型计算的湍流粘性较小.     

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

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

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

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

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

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

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

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

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

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

雷诺应力模型在三维湍流流场计算中的应用

从雷诺应力模型出发,通过求解雷诺平均N-S方程组获得三维湍流流场的数值解。计算中比较了多种湍流模式,并进行了相应的流场计算。本文完成了两个典型算例。从算例1的计算结果与实验值比较中发现:采用雷诺应力模型(RSM)计算的三维流场比采用k-ε模型更贴近实验值;算例2采用了RSM模型及三维非结构网格,对一典型内流问题进行了三维流场计算。算例的数值实践表明:采用雷诺应力模型可以有效的计算各向异性的湍流流场;另外,发展非结构网格有助于模拟壁面附近的流动,并节省计算机内存。      

运用低雷诺数非线性k-ε模型研究湍流分离流动

给出了一个低雷诺数非线性ｋ－ε湍流模型，其中采用了一个非线性雷诺应力表达式并引入了两个考虑湍流近壁效应的修正函数。应用该模型和标准ｋ－ε模型计算了两种不同膨胀比的后台阶分离流。与实测结果比较表明，低雷诺数非线性模型结果较标准ｋ－ε模型结果有明显改进，它对流场的各特征量包括分离—再附长度、时均速度及各雷诺应力分量均给出了较好的预测。     

计算风工程中k-ε模型的一类边界条件

给定适当的湍流来流边界条件是计算风工程中一个重要的问题。然而通常情况下，由于问题的复杂性，很难给定一个既满足边界层自保持要求又与试验相符的边界条件。本文首先从模型方程本身出发，推导出一类近似满足k-ε模型自保持边界条件要求的湍动能表达式。此表达式与有关文献建议的湍流边界条件比较，可以考虑湍动能随高度的变化，因而更符合实测规律；然后基于风洞试验数据，定义了一组k-ε模型中的模型常数；最后通过模拟风洞试验的算例，验证并分析了此类湍动能边界条件配合试验常数，在标准k-ε模型中模拟均匀大气边界层的适用性。     

应用修正的k-ε模型研究超声速H2/Air燃烧

为了研究超声速燃烧中流体可压缩性的影响,对标准k-ε湍流模型进行可压缩性修正(包括结构可压缩性修正和膨胀可压缩性修正两部分)。分别应用标准k-ε模型、修正的k-ε模型和雷诺应力模型(RSM),考虑氢气/空气详细化学反应机理(GR I-M ech 2.11机理,10组分,28基元反应),数值模拟有壁面限制的超声速混合层冷态及热态流场。结果表明:壁面和燃烧对湍流影响都很大;修正模型对冷态以及燃烧场的预测结果优于其它两个;修正模型预测的混合层厚度更薄,燃烧区域更窄,与实验结果吻合地更好。     

液体煤油超声速燃烧数值分析的湍流模型

利用计算流体力学软件对煤油在所设计的双模态燃烧室内的喷雾燃烧进行了数值模拟。采用k-ε三种形式、k-ω两种形式和雷诺应力湍流模型,计算了在飞行马赫数为5,煤油与空气的当量比为0.551时的情况。通过将数值模拟结果与实验数据进行对比,证明这些湍流模型都可用于煤油超声速燃烧的数值计算。并且,k-εRNG模型和k-ωSST模型的计算结果与实验数据最为接近,可以认为这两个模型在模拟煤油超声速燃烧流场中是非常有效的;而雷诺应力模型效果较差。     

两种k-ω型湍流模型在无网格方法中的应用研究

基于在高度各向异性点云结构中取得良好结果的点云重构技术,采用移动最小二乘无网格法解耦计算了k-ω型湍流模型方程和雷诺平均N-S方程(RANS),实现了对湍流流动的数值模拟。通过对NACA0012翼型、RAE2822翼型的绕流以及后向台阶大分离流动的数值模拟,比较分析了k-ωSST(Shear Stress Transport)和k-ωTNT(Turbulent/Non-Turbulent)两种模型的湍流预测能力。结果表明:对弱逆压梯度情况下的湍流,两者具有相当的预测能力;而k-ωSST较k-ωTNT具有更好的激波捕捉能力。对于大分离流动,两种模型较实验数据均存在一定差异,但总体趋势吻合良好。     

Sajben扩压器复杂流动中湍流模型的性能评估

用8个常用的湍流模型对Sajben扩压器中跨声速流动进行了数值模拟，评估了Spalart-Allmaras， 标准k-ε， RNG (re-normalization group) k-ε，realizable k-ε，标准k-ω，SST(shear stress transport) k-ω，v2-f，Reynolds stress共8个湍流模型对激波/湍流边界层相互作用的模拟预测能力.通过与实验数据比较发现:SST k-ω模型和v2-f模型比其他模型模拟的更准确，其中SST k-ω模型比v2-f更能准确地预测壁面压力，然而对于分离点、再附点以及分离区长度v2-f比SST k-ω预测得更准确.      

气固两相流k-ε-kp-εp双流体模型及应用

以流体ｋ－ε双方程模型为基础，提出了一种包括颗粒湍动能输运方程和耗散率方程的ｋ－ε－ｋｐ－εｐ两相湍流模型。并对两个垂直上升管道内气固两相湍流进行了数值计算，计算结果和已有实验结果吻合得很好。     

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

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