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

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

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

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

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

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

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

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

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

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

用双流体模型模拟湍流两相流场

首先针对颗粒浓度在两相流动中的重要性，建立了颗粒体积浓度方程，阐述了颗粒浓度在求解过程中的特点及颗粒体积浓度方程在计算方法中的贯穿。其次计算了沙层上的稀相可压气固两相流、直管内的不同压液固两相流及叶片扩压器中的密相可压力固两相流的流动，发现了颗粒浓度两种截然不同的分布，并验证了用不可压或稀相的模型来计算可压或密相的模型会带来误差，说明了在计算可压密相两相流动时，可压及密相的因素必须考虑到。     

用非线性K—ε两方程湍流模型计算后向台阶分离流动

本文利用ＢＧａｔｓｋｉ和ＣＧＳｐｅｚｉａｌｅ根据张量不变性理论推出的二维流动站情况下的非线性Ｋ－ε两方程模型来对这两种不同扩比下的后向台阶分离湍流进行了数值模拟,从计算结果可以看出,用非线性Ｋ－ε两方程模型计算得到的回流区长度,壁面压力分布和摩阻分布等物理量的精度比用标准Ｋ－ε两方程模型计算得到的结果用了明显改进,而且,该模型还具有标准Ｋ－ε两方程模型的结构简单,使用方便的优点,易于工程应用。     

用一种非线性k-ε模型计算突扩台阶后的流动

对一种新的非线性k-ε模型的无量纲系数进行了重新拟合.用新的系数计算了槽道流动、方截面管流和等应变率,平均流场无旋流动的正应力分布.计算结果与实验资料符合很好.将该模型并入常规的k-ε模型的计算程序中,并对突扩台阶后的湍流流动进行了数值计算,同线性k-ε模型的计算结果和实验结果进行了比较.结果表明,这种非线性k-ε模型同线性模型相比有较大的改进,同已有的其它非线性模型相比,也有一定的优点.     

大功率氧煤燃烧器内湍流气固两相流动数值模拟

本文对新型同轴射流缩口式氧煤枪所组成的大功率氧煤燃烧器内冷态轴对称二维湍流气固两相流动用ｋ－ｅ－ＡＰ两相湍流模型进行了数值预报，揭示了燃烧器内两相流场特征，预报结果和实验结果相符。进一步对相应高炉实际工况高速条件下的冷态两相流场进行了预报，其结果为优选结构尺寸和流动参数提供了重要依据。     

基于非交错有限体积格式和 k-ω 湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

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

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

用湍流传热模型模拟室内空气流动和换热

利用低雷诺数NH模型、NK传热模型及模型系数改进后的NK传热模型,分别对室内空气流动和换热的速度场、温度场、湍动能场进行了模拟,并将计算结果与实验结果进行了对比分析.研究表明,低雷诺数NH模型、NK传热模型及模型系数改进后的NK传热模型均能对室内空气流动及换热的速度场、温度场、湍动能场进行较好的数值模拟,其中模型系数改进后的NK传热模型模拟的准确性最高.     

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

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

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.      

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

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

基于非交错有限体积式和k—w湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

大涡模拟超声速两相流混合层的非线性分析

为研究超声速两相流混合机理,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构.用二维盒式滤波法处理可压缩气相N-S方程,得到亚格子形式的LES控制方程;用轨道模型仿真液滴运动,并用单液滴蒸发模型模拟液滴蒸发过程.利用伪邻近法和互信息量法确定了相空间重构的嵌入维数和时间延迟,对混合层压力时间序列进行了相空间重构.分析了仿真结果中混合层的涡结构特性和液滴对流场的反作用效果,并讨论了混合层失稳过程的非线性特征.发现混合层发展的不同阶段与其系统的非线性特征值密切相关.      

带有剪切应力输运性质的kGε两方程湍流模型构造与应用

由于传统的壁面衰减函数峰值并不考虑雷诺数的变化,因此仅能在特定的雷诺数下还原边界层的近壁衰减效应。为了改善这个缺陷,一个结合不同湍流雷诺数的因子被引入壁面衰减函数,从而使其峰值对雷诺数的变化能够合理感知。另外,基于DNS数据,将Bradshaw假设中的常数进行了重新标定,引入平滑的湍涡黏性系数切换函数,最终形成了一个带有剪切应力输运性质的k-ε两方程湍流模型。通过平板、翼型、二维鼓包以及翼身组合体进行了校验,结果显示新模型的精度较高,具有一定的工程实用价值。     

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

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