气固两相流k-E-kp-Ep双流体模型及应用

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

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

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

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

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

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

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

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

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

突扩回流与大速差射流回流湍流气－固两相流动的数值模拟

本文建立了气－固两相流的多连续介质模型，并分别对突扩回流与大速差射流回流这两种与实际工程应用密切相关的复杂流动条件下的受限湍流气一固两相流动进行了数值模拟。单相湍流回流流场的模拟结果与实验符合良好。两相流动的模拟结果表明，与突扩回流相比，大速差射流回流更有利于实现较高的颗粒浓度与气相回流区及低速尾流区的匹配，增加两相滑移速度，延长颗粒停留时间及强化两相间的混合。     

突扩回流与大速差射流回流湍流气—固两相流动的数值模拟

本文建立了气－固上流的多连续介质模型，并分别对突扩回流与大速差射流回流这两种与实际工程应用密切相关的复杂流动条件下的受限湍流气－固两相流动进行了数值模拟。单相湍流回流流场的模拟结果与实验符合良好。两相流动的模拟结果表明，与突扩回流相比，大速差射流回流更有利于实现较高的颗粒浓度与气相回流区及低速尾流区的匹配，增加两相滑移速度。延长颗粒停留时间及强化两相间的混合。     

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

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

固液火箭冲压发动机内两相反应流场数值计算

为了计算固液混合式火箭冲压发动机补燃室内的三维反应流场，用块隐式法求解气相Navier-Stokes方程组，用连续介质模型和k-ε-Ap模型计算颗粒相的湍流流动与蒸发过程，用修正的k-ε-g模型描述燃料的燃烧。为了分析发动机设计参数对反应流场的影响，用不同的条件进行计算，并由此分析了补燃室几何结构和液体燃料切始颗粒直径对燃烧效率的影响。算例表明，计算方法有效可行。数值结果能够反映流场结构、液体燃料的蒸发和两种燃料的燃烧过程。     

CFD软件-PHOENICS及在燃气涡轮流场计算中的应用

了解不同CFD软件,是选好CFD软件、并经济地解决问题的关键.本文从离散方法、求解过程、湍流模型、两相流处理、结果输出等方面对Phoenics软件作一解读,并给出某燃气轮机内气固两相流动的流场计算和颗粒轨迹计算实例,以Eulerian-Lagrangian方法的成功应用,验证了Phoenics软件较强的适用性.     

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

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

90°强曲率弯管内三维湍流流场的数值分析

本文根据作者发展的一种能计算任意截面形状 90°弯管内三维不可压缩湍流流场的分析方法,利用作者实验的圆截面弯管数据 ( R_c/D=0.87)检验 k- ε双方程湍流模型预示强曲率三维湍流流场的精度。计算了 R_c/D=2.3的矩形截面和 R_c/D=0.87的圆截面 90°弯管,并将计算结果与相应实验数据进行了比较分析。结果表明,即使对曲率很强的弯管,吸力面 (内侧壁 )不出现二次流迁移所形成的低速区前,标准 k- ε双方程模型仍能较好地预示出平均速度场。      

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

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

密相、湍流、可压、气固两相流理论及在管道中的应用

就密相、湍流、可压、气固两相流的基本方程用双流体模型进行了较说尽地归纳。与单相流体相比，在两相流的基本方程中，多了一个变量（体积浓度）及两相间的耦合作用，并用CDIPSA方法对可压密相上升管道湍流流动进行了计算，结果表明：随颗粒直升的增加，使得单个颗粒脉动速度减小；但是气体及颗粒的速度差却随之增大。     

Numerical study of compression corner flowfield using Gao-Yong turbulence model

A numerical simulation of shock wave turbulent boundary layer interaction induced by a 24° compression corner based on Gao-Yong compressible turbulence model was presented.The convection terms and the diffusion terms were calculated using the second-order AUSM (advection upstream splitting method) scheme and the second-order central difference scheme,respectively.The Runge-Kutta time marching method was employed to solve the governing equations for steady state solutions.Significant flow separation-region which indicates highly non-isotropic turbulence structure has been found in the present work due to intensity interaction under the 24° compression corner.Comparisons between the calculated results and experimental data have been carried out,including surface pressure distribution,boundary-layer static pressure profiles and mean velocity profiles.The numerical results agree well with the experimental values,which indicate Gao-Yong compressible turbulence model is suitable for the prediction of shock wave turbulent boundary layer interaction in two-dimensional compression corner flows.      

旋转离心叶轮内湍流的非结构化网格数值计算

给出了同位非结构化网格数值求解旋转离心叶轮内湍流流场的SIMPLE算法。网格划分采用改进的阵面推进法。控制方程采用通用的对流扩散方程并利用有限控制容积法离散。湍流模型采用标准的k-ε模型。形成的代数方程组用带有预条件器的共轭梯度平方法（CGS）求解。最后利用所给出的算法对一实验用的旋转离心叶轮进行计算。计算结果和实验数据基本相符,说明该算法合理,编程正确,可作为流场预测的理论工具。     

带气膜冷却的火焰筒壁温的数值分析

 采用数值方法，预测了由两通道两股气流构成的二维平缝气膜冷却的火焰筒壁温，获得了整个流场、温度场以及相应的流动和换热特性。在计算中，对紊流采用“ｋ ε”两方程模型处理；在固壁附近运用壁面函数法；运用热流法计算了高温气体的热辐射；对两种不同介质的交接面的导热系数运用调和平均法。最后与实验结果相比较，得到良好的一致性。     

k-ω SST两方程湍流模型中参数影响的初步分析

Menter的k-ωSST两方程湍流模型在流体力学计算中有良好表现。利用均匀试验设计方法对k-ωSST模型中八个参数在小攻角亚声速、跨声速和大攻角三种典型的翼型绕流流场计算中的影响规律进行了分析。分析结果表明在附着流中参数a1是对流场的影响最大的参数;而在分离流中,参数a1对流场的影响与湍流产生与耗散等项的作用相当,模型中的四个参数:σω2、β2、β*和a1的取值都会对计算结果有较大的影响。     

非线性 k-ε 紊流模型在平面叶栅紊流计算中的应用

 采用Speziale的非线性k-ε模型计算了平面叶栅内的紊流流动,并且与标准的k-ε模型的计算结果作了对比,表明非线性模型同样可以应用于具有周期性边界条件的复杂流动,在紊流量的计算结果上表现出了明显的各向异性。