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

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

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

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

加力室液雾燃烧的模型与应用

提出了描写航空发动机加力燃烧室内气液两相湍流燃烧的控制方程,采用了双流体模型,高温辐射通量模型,双时间尺度湍流K-ε模型,应用SIMPLE方法解这组方程,对二维轴对称的涡喷加力燃烧室进行计算,证明应用本文的液雾燃烧模型可以得出合理的计算结果。      

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

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

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

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

显式代数雷诺应力模型在超声速/高超声速流场计算中的应用

为了提高显式代数雷诺应力模型在高超声速流场中的计算精度,以修正后的CG k-ε模型作为湍流时间尺度决定方程,将Wallin和Johansson显式代数雷诺应力模型(WJ2000模型)应用于超声速/高超声速流场数值计算,并与线性湍流模型的计算结果进行了对比分析。结果表明WJ2000模型对高超声速压缩拐角附近壁面压力和热流分布的预测要优于线性湍流模型,同时网格依赖性明显低于后者;对CG k-ε湍流时间尺度决定方程所做修正不影响WJ2000模型对超声速流场的预测结果,但使其对高超声速流场的预测精度明显提高。     

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

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

针对叶轮机流场分析的不同Spalart-Allmaras湍流模型评价

Spalart-Allmaras（SA）湍流模型在Reynolds-Average Navier-Stokes（RANS）方程求解中得到了广泛的应用。针对叶轮机数值模拟，国内外学者提出了很多修正的SA模型，但它们之间的对比少有研究。为此，综述了原始的SA（SA-standard），SA-neg，SA-Helicity，SA-Ning及SA-R五种不同SA湍流模型。针对NASA Rotor 67和NASA Rotor 37，考察这五种湍流模型对数值计算稳定性的影响及对叶轮机内部流动细节捕捉的能力。研究表明：SA-standard和SA-neg湍流模型计算结果几乎相同；相比于SA-standard湍流模型，SA-Ning湍流模型计算的压比、效率和堵塞流量均较大；SA-R湍流模型模拟的激波/边界层干涉更严重，分离区域更大，计算的压比小，模型修正效果不理想；SA-Helicity湍流模型不仅能极大地提高计算的压比和流量范围，而且能提高计算的稳定性和失速工况下的计算准确性，但计算的堵塞流量和效率较小。     

液氢/液氧火箭发动机喷雾燃烧过程三维数值模拟

用三维湍流Ｎ－Ｓ方程及拟流体模型描述液体火箭发动机内部喷雾两相燃烧流动过程,两相之间的质量交换采用蒸发模型计算,气相化学反应速率采用Ａｒｒｈｅｎｉｕｓ公式计算。采用ＳＩＭＰＬＥＲ与ＩＰＳＡ算法求解两相控制方程并对氢氧发动机燃烧流动过程进行了三维数值分析,得到了发动机推力室内燃气参数的详细分布情况。结果表明对氢氧发动机而言,蒸发过程非常快,推进剂混合过程是燃烧过程的决定因素。     

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

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

基于S-A湍流模型和间歇因子输运方程的转捩流数值模拟

考虑到γ-Re_θt转捩模型中间隙因子输运方程并不依赖于具体的湍流模型,因而提出了耦合S-A湍流模型和间歇因子输运方程的转捩流模拟方法.其主要思想是构造与间歇因子相关的两个耦合函数,分别作用到S-A模型的生成项和耗散项用以控制湍流的产生与发展,从而实现层流到湍流的数值转捩.结果表明:该方法中模型方程具有完全当地性的特点,易于实施,可直接用以求解三维转捩流动.4个典型算例表明该方法模拟结果与γ-Re_θt模型准确度相当,但少求解了两个输运方程,计算耗时减少了15%.      

具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

离心叶轮内稀释固液两相流动的数值模拟

 对固液两相流动进行了数值模拟。液相和固相分别采用 k- ε湍流模型和代数颗粒流( AP)湍流模型并考虑了液固两相间的相互作用。用此模型结合 SIMPLE- C算法对离心叶轮内的稀释固液两相流动进行了计算并研究了颗粒对流场的影响。计算结果同试验结果基本一致     

两种湍流模型在跨声速绕流计算的应用研究

湍流模型在对复杂流场的数值计算中起着非常重要的作用.本文采用SA一方程湍流模型和SST k-ω二方程湍流模型,通过解耦求解雷诺平均N-S方程和湍流模型方程实现对亚跨超声速湍流流场的数值模拟.对NACA0012翼型和ONERA-M6机翼跨声速绕流流场进行了计算,对压力分布和激波位置与实验结果进行了细致的比较,并分析了不同离散格式、不同网格疏密及壁面函数对计算结果的影响,在计算过程中这两个模型体现出了较好的简捷性和健壮性.     

前后翼型装药结构SRM气固两相流动特性

基于可压缩流动强守恒型Navier-Stokes方程,运用颗粒轨道模型(PTM)和shear stress transfer(SST)k-ω湍流模型,采用计算单元内颗粒源法(PSIC)进行气固两相耦合计算,建立了某型前后翼型装药结构固体火箭发动机(SRM)工作初期的三维两相内流场数值计算模型.对比分析了纯气相和气固两相...     

三维加力燃烧室两相湍流燃烧的数值模拟

本文对涡扇发动机三维加力燃烧室内的两相湍流燃烧过程进行了数值模拟。两相流模型采用双流体模型,湍流模型采用标准 k-ε模型,湍流燃烧采用涡旋破碎 ( EBU)模型,数值方法采用 LEAGAP算法。计算结果定性合理。      

基于S-A与SSG/LRR-ω两种湍流模型的CHN-T1标模计算与分析

通过求解雷诺平均Navier-Stokes方程,对第一届AeCW-1研讨会上要求的所有工况进行了计算,着重研究S-A湍流模型和雷诺应力模型在CHN-T1标模上的表现,评估两种湍流模型的网格收敛性,对抖振特性的模拟能力以及对雷诺数的敏感性。根据与实验数据的对比分析,得出主要结论如下:两种湍流模型都能趋向于网格收敛,雷诺应力模型的网格收敛更快。在较大迎角的工况,两种湍流模型对机翼上激波诱导的分离流有较大的差别。CHN-T1标模的尾撑,改变了飞机平尾的压力分布,并且在飞机机尾上产生复杂的三维分离流动,对俯仰力矩系数的准确计算产生了较大的影响。考虑静气弹效应,两种湍流模型计算的力矩系数比不考虑静气弹效应更大,雷诺应力模型的偏差更大。在高雷诺数工况,两种湍流模型的计算结果趋于一致。     

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

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

用椭圆型、抛物型两种方法计算平面叶栅湍流尾迹

本文用连续方程、稳定不可压缩椭圆型和抛物型时均N-S方程及K-ε双方程湍流模型计算了平面叶栅湍流尾迹。计算表明,上述两种方法都能得到较好的数值结果,但抛物方法可节省大量计算时间和内存空间。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.