二维气固混合层中拟序结构对 固粒运动影响的研究＊

本文对固粒在混合层涡卷起和配对过程中的运动进行了研究。对流场用拟谱方法直接数值模拟,对固粒用颗粒－轨道模型。计算结果给出了涡卷起和涡配对过程中固粒的运动情况,说明本流场Ｂａｓｓｅｔ力对固粒运动的影响甚微;拟序结构演变过程的不同阶段造成不同的固粒分布,该结果与流场显示的实验结果相符;在相同流场条件下,Ｓｔｏｋｅｓ数是影响固粒在拟序结构流场中运动的控制参数。文中所得结论对工程应用具有指导意义。     

稀薄流过渡区气固两相喷流的建模与数值模拟

从描述粒子运动的细观层次出发,运用DSMC方法,构造了稀薄流过渡区气固相互作用模型和颗粒碰撞模型,数值模拟了有气固两相横向喷流干扰的二维有限平板高超声速流场,分析了流场结构和物面气动系数分布等特性.结果表明,对于不同的来流Knudsen数和颗粒质量载荷等参数,喷流干扰流场存在差异,同时反映了基于DSMC方法的气固相互作用模型是研究稀薄及过渡条件下的气固两相流动的有效方法.     

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

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

固体火箭喷管颗粒尺寸分级两相跨音速流场计算

本文中对固体火箭喷管颗粒尺寸分级的两相跨音速流场作了计算.气相控制方程采用隐式近似因子分解法求解,尺寸分级的颗粒控制方程采用特征线求解,然后,二者进行充分的耦合,可以获得固体火箭发动机含有任意颗粒质量分数和不同颗粒尺寸时轴对称喷管跨音速流场的参数分布.文中讨论了不同颗粒半径和质量分数对流场的影响,对单一颗粒尺寸和颗粒尺寸分级的参数进行了比较.两相耦合计算的迭代收敛速度取决于气相,本文中气相方程求解的格式除部分边界外是隐式的,CFL数可取至6左右,收敛速度快.特别是对颗粒尺寸分级的计算,得益更大,其得益的倍数为颗粒的分级数.     

入口温度剖面对喷管流场结构的影响

应用质量平均的Navier-Stokes方程和B－L代数湍流模型，对超燃冲压发动机尾喷管的流场进行了数值模拟研究，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散，通过数值模拟，对比研究了温度非均匀性，三维效应对尾喷管的流场结构的影响。     

星孔型装药发动机三维两相流场的数值模拟

为了研究颗粒在星孔型装药固体火箭发动机燃烧室和喷管中的运动轨迹以及颗粒与发动机壁面的碰撞情况,针对可压两相流动,采用了高雷诺数下的k ε湍流模型和欧拉 拉格朗日两相流模型,用全速度SIMPLE方法对方程组进行求解,并用PSIC方法进行气固耦合计算。计算得出了流场内两相的速度、温度等参数的分布及多种情况下固体颗粒的运动轨迹。在燃气生成量确定的情况下,从距离喷管较近的某些位置进入流场的颗粒比较容易撞击壁面;颗粒的尺寸和局部产生的旋涡对颗粒的轨迹和碰撞也会产生较大的影响。     

导弹垂直发射燃气排导系统两相流场数值分析

根据燃气排导系统内两相流场的分布特点,利用可压缩雷诺 N-S平均方程对流场进行计算。湍流模型采用标准 k-ε二方程模型,采用颗粒轨道模型对颗粒相进行模拟,利用 FLUENT就发射高度及燃烧室压强对三维流场参数分布的影响进行了重点分析,为优化发射结构设计提供了重要的理论依据。     

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

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

冲压发动机伴随蒸发过程的两相流流场计算

为了计算固液混合式火箭冲压发动机补燃室内带有燃料液滴的流场,用块隐式法计算气相流场,用分离流连续介质模型计算颗粒相的流动与蒸发过程,在控制方程组的联合求荽中,采用特定的算法来体现颗粒相的质量守恒律与尺寸的变化。算例表明,迭代过程能够快速收敛,计算结果定性正确。     

前缘对进气道性能影响的数值模拟

应用质量平均的N-S方程和B-L代数湍流模型，在不同压缩面前缘半径条件下，以高超声速进气道的流场进行了数值模拟，对比研究了进气道压前缘半径对流场结构及性能的影响，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散。     

不同密度材料颗粒对两相湍流调制的实验

 在颗粒雷诺数较小范围内,针对不同密度材料颗粒对气相湍流的调制规律进行了实验研究。采用相多普勒颗粒分析仪(PDPA)测量了气粒两相圆湍射流中气相速度和湍流度的分布,比较了相同形状、相同粒径和体积分数近似相同时玻璃微珠和聚苯乙烯颗粒分别存在时的两相射流流动的气相速度和湍流度。实验发现:不同密度材料的颗粒对气相平均速度的调制作用差别不明显,但对气相湍流度的调制作用显著不同。在近场区域,密度大的玻璃微珠对气相湍流度的增强作用高于密度小的聚苯乙烯颗粒;在远场区域,气相湍流度被削弱,并且密度大的玻璃微珠的削弱作用强于密度小的聚苯乙烯颗粒。该研究为进一步提出合理的两相湍射流的控制方法提供了理论依据和指导。     

气固两相湍流场纳米颗粒演变特性综述

含纳米颗粒的气固两相湍流场在包括航空等众多领域中很常见,以单体、聚集体和团聚体不同形式存在的纳米颗粒在流场中经过生成、对流、扩散、凝并、破碎等过程,其数密度、尺度、尺度分散度等将发生变化。本文就以上相关研究状况进行了回顾,说明颗粒生成是气相化学反应产生的可冷凝蒸汽物质因表面冷却、绝热膨胀或混合、湍流混合或化学过程产生的过饱和所导致;导致颗粒凝并的原因包括布朗运动、湍流剪切、速度梯度、差异沉降;颗粒的凝并取决于颗粒的尺度和流场的特性,并受初始颗粒分布及湍流扩散控制;湍流场对颗粒凝并的影响除了湍流强度的因素外,还体现在由湍流脉动所引发的颗粒数密度的脉动;颗粒凝并后形成尺度较大的团聚体容易在流场剪切和其他因素作用下发生破碎;剪切破碎是导致颗粒破碎的主要因素,有效破碎系数取决于剪切率和颗粒的体积分数;颗粒的沉降取决于颗粒尺度、形状和流体性质等因素;导致颗粒沉降的因素有重力、扩散、惯性撞击、电场和热迁移等;当存在温度梯度时,热泳力对颗粒沉降也起到重要作用。本文最后提出了有待进一步研究的若干问题。     

航空发动机尾喷管固体颗粒流动特性研究

为了通过发动机尾喷管内固体颗粒的荷电水平来判断发动机的故障程度,对发动机尾喷管内固体颗粒的流动特性进行研究。选取某型发动机尾喷管为研究对象,利用ICEM软件进行建模与网格划分,通过Fluent模拟仿真分析颗粒直径和初始速度对颗粒分布的影响,探讨固相颗粒产生绳流、均匀弥散流、环流及层流等流型的条件。研究表明:局部注入的颗粒分布为绳流;对于面均匀注入方式,随着固体颗粒直径增大,其分布由均匀弥散流逐渐变为环流;随着颗粒直径和速度的增加,颗粒与尾喷管壁面碰撞加剧,部分区域出现层流;当颗粒直径和速度达到一定程度时,颗粒的分布规律保持稳定。     

激波与固体颗粒群相互作用实验研究

为了研究FAE武器、温压武器中燃料抛撒问题,在垂直激波管里进行了气-固两相流的实验研究。主要研究激波（Ma=1．96）与3种不同堆密度的石英砂颗粒群相互作用,用高速摄影仪对颗粒在抛撒过程进行图像捕捉,得到颗粒在抛撒阶段的速度。研究颗粒群在超声速气流下的运动规律以及激波对颗粒群加速的影响规律,分析激波通过颗粒群时气-固两相流中的动力学过程。定量分析颗粒堆密度对激波的影响,得出了固体颗粒群受到激波作用后抛撒过程中纵向和横向颗粒云团体积增加,颗粒云团平均浓度减少;颗粒群的堆密度与颗粒粒径、抛撒初速度、反射激波的强度都成反相关关系。     

动能BGK算法在近连续流模拟中的应用

首先比较详细地介绍了动能BGK算法的基本思想,然后应用该算法,数值模拟了克努森数(Knudsen number)为0.0001～1.56,马赫数为1.96绕圆柱的二维粘性流动,研究了滑移边界条件和无滑移边界条件对圆柱阻力的影响,通过与实验结果[5]进行比较,验证了文献[1]关于由连续流到自由分子流的这种流动模型的划分.计算结果表明:采用滑移边界条件,直到克努森数为0.3,计算得到的阻力与实验都符合得很好.因此,采用BGK算法加滑移边界条件,能够模拟由连续流到近连续区的流动.     

Tracking characteristics of tracer particles for PIV measurements in supersonic flows

The tracking characteristics of tracer particles for particle image velocimetry (PIV) mea-surements in supersonic flows were investigated. The experimental tests were conducted at Mach number 4 in Multi-Mach Wind Tunnel (MMWT) of Shanghai Jiao Tong University. The motion of tracer particles carried by the supersonic flow across shockwaves was theoretically modelled, and then their aerodynamic characteristics with compressibility and rarefaction effects were evalu-ated. According to the proposed selection criterion of tracer particles, the PIV measured results clearly identified that the shockwave amplitude is in good agreement with theory and Schlieren visu-alizations. For the tracer particles in nanoscales, their effective aerodynamic sizes in the diagnostic zone can be faithfully estimated to characterize the tracking capability and dispersity performance based on their relaxation motion across oblique shockwaves. On the other hand, the seeding system enabled the tracer particles well-controlled and repeatable dispersity against the storage and humidity.     

固体发动机喷管内两相流的有限元数值模拟

本文给出了固体发动机喷管内两相流的数值模似。从非守恒型两相流甚本方程出发,采用加权积分技术,导出了离散点的流动物理量的有限元数值计算公式。数值模拟分别与JPL喷管的两相流壁面静压分布测量结果以及与差分法计算结果作了比较,比较结果表明:二者是吻合的。本文还以中型喷管为例,计算了粒子的轨迹,粒子的速度滞后和燃气的温度滞后,计算结果表明;入口粒子冲击角对粒子轨迹有明显的影响。     

长尾喷管故障诊断中的两相流动计算

为了给长尾喷管故障诊断提供理论依据并提供改进措施,利用颗粒轨道模型和有限体积的Jameson格式计算了长尾喷管两相流流场,得到了流场中温度和马赫数的分布和不同燃烧室总压下不同尺寸的粒子对喷管壁面的撞击情况,结果表明大尺寸粒子对直管段壁面的撞击位置存在一个集中点,这是造成喷管烧穿的重要原因。得出了总压越高粒子撞击壁面的次数越多、尺寸越大粒子撞击壁面的次数越多的结论,给出了提高长尾喷管热防护结构可靠性的一些措施。      

高速流动中PIV示踪粒子松弛特性研究

示踪粒子的随流能力是PIV技术在高速流动中应用的关键点之一。在上海交通大学变马赫数高速风洞中开展了PIV实验研究,重点提出一种评价示踪粒子随流能力的松弛特性分析模型。在马赫数4条件下尖锥、尖劈等模型PIV实验研究中,可以准确分析粒子的松弛特性,粒径分析结果与实验吻合较好,并验证了高速流动PIv的测试精度和示踪粒子布撒能力。     

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

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