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

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

基于CFD耦合DEM并行算法的旋翼飞行器桨叶沙盲现象模拟分析

旋翼飞行器处于地效（IGE）悬停状态易吸入沙尘从而引起沙盲现象。为了研究出现沙盲现象时的悬停流场与沙尘的运动规律,采用基于雷诺平均 N-S 方程及 k-ω（SST）湍流模型,结合计算流体力学（CFD）与离散元（DEM）耦合并行算法,使用 Fluent 软件进行流场连续相计算;通过 API 传递流场信息至离散相,计算直升机旋翼的悬停流场,并与可得到的 PIV 测量结果进行对比;再将流场动量数据传输给离散相,计算沙尘的完整运动轨迹和颗粒分布。结果表明：在对流场精确模拟的同时,采用离散元方法对沙尘颗粒受力运动以及空间分布较好地模拟,能够多参数化地了解地效流场和细微沙尘颗粒受力运动规律,直观呈现沙盲的发展过程。     

直接模拟蒙特卡罗羽流模拟的两相作用模型

为了数值模拟包含气体分子和凝相颗粒的两相羽流场,使用自主开发的羽流计算软件(PWS)对两相作用模型进行研究。计算均匀气流对颗粒的加速过程,结果与理论解一致,证明气体分子对颗粒的单向作用模型合理。计算气体分子在颗粒表面的反射时,将全局坐标系转换到颗粒坐标系确定其法线,再转换到法线坐标系进行镜面反射或漫反射。将这种间接法与使用反射向量的分布函数的直接法进行了对比。对同一个气体分子在颗粒表面的反射进行1×106次数值试验,两种方法所得镜面反射的速度向量均符合空间均匀分布,漫反射的速度向量分布相互一致。使用两种方法模拟150N液体火箭发动机的两相羽流场,所得气相羽流场几乎相同,间接法计算效率较高。液滴颗粒对气流场的影响主要集中在轴线附近,与无液滴的气流场相比,有液滴时会使气流温度升高、压强增大、马赫数减小。     

气固紊流剪切流中颗粒弥散的拉格朗日模拟

本文提出了一种对于气固两相紊流剪切流中圆形固体颗粒弥散的拉格朗日拟计算方法，其中考虑了颗粒间的磁撞对流体相和颗料相的影响，应用该方法对一气固紊流剪切流场进行了模拟计算，并对有、无颗粒间磁撞情况下的模拟计算结果与Lavieville用大涡模拟方法的研究结果进行了比较，并进行了讨论。     

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

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

固体火箭发动机塞式喷管两相流场与性能分析

为了解固体发动机塞式喷管中两相流场的特点和性能的情况,采用欧拉-拉格朗日两相方法对二维全长塞式喷管的两相流场进行了计算.将纯气相流场和两相流场进行了比较,结果表明由于颗粒相的存在使气流的比热比减小,使得低空环境中,塞锥表面的压强峰后移;高空环境中,塞锥表面的压强会升高.大直径颗粒的运动轨迹受外界反压和塞锥气流变化的影响小,小直径颗粒的运动轨迹受外界反压变化和塞锥气流变化的影响大;在两相流性能方面塞式喷管与钟型喷管类似,微粒直径较大的情况下,塞式喷管的效率较低;塞式喷管的效率随着颗粒含量的升高而降低.     

双工况氢氧发动机燃烧与传热数值分析

应用三维湍流N-S方程以及颗粒轨道模型描述双工况氢氧发动机内部喷雾两相燃烧流动过程。两相之间的质量、能量交换由液滴蒸发模型计算，气相化学反应速率由Arrhenius公式计算。通过耦合求解气液两相的模型方程，对发动机转工况前后的三维流场进行了数值计算，并耦合计算了燃气与壁面之间的传热以确定壁面的温度和热流分布。另外还对分别采用同轴离心式喷嘴和直流式喷嘴得到的燃烧流场与燃烧效率进行了比较。计算结果表明转工况前的壁面温度与热流都比转工况后大。离心式喷嘴的雾化混合效果与燃烧效率都比直流式喷嘴好。     

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

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

两相流喷管扩散段型面优化计算

本文论述了用两相喷管流动计算程序按直接优化法寻找两相流喷管最佳扩散段型面的方法,着重讨论了两相超音速流场计算中的网格结点、颗粒阻力系数和传热系数。算例中给出了抛物线、圆弧线和三次曲线喷管的计算结果。算例指出;喷管扩散段采用不同的曲线族,其最佳型面的初始膨胀角是不同的。      

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

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

面对面呼吸飞沫传播和防护的流体力学初步分析

针对SARS-CoV-2疫情,围绕病人和戴防飞沫面罩的健康人面对面呼吸状态下飞沫的流动传播问题,利用计算流体力学方法进行了数值模拟。首先采用二维非稳态不可压缩雷诺平均N-S方程描述平均流场,计算区域采用H型结构网格拓扑结构,将多相流的欧拉模型用于描述空气-飞沫流动,湍流脉动采用RNG k-ε湍流模型。计算获得了呼吸飞沫喷射速度、粒径和入口含量对飞沫流动传播的影响规律,分析了水平传播线和垂直传播线的飞沫浓度分布。研究表明:呼吸飞沫喷射速度、入口含量和流动传播速度、浓度分布呈正相关性,粒径主要影响飞沫流动传播的浓度峰值大小;呼吸飞沫跟随空气流动先以水平传播为主,后转变为以垂直扩散为主,扩散过程会形成“猫眼”状飞沫云团;呼吸飞沫在气流作用下侵入到防飞沫面罩内形成旋涡,导致防飞沫面罩的失效,因此,防飞沫面罩必须结合口罩共同使用。     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

环管燃烧室热态流场的数值计算及性能分析

借助流体计算软件对某环管燃烧室内的热态流场进行了数值分析。气流的湍流流动采用Realizab lek-ε湍流模型进行模拟,用拉格朗日法跟踪油滴的轨迹,非预混PDF燃烧模型模拟气液两相的燃烧过程,辐射传热通过离散坐标模型来考虑。分析了气流流动规律、燃油液滴的运动轨迹、温度场的分布,获得了总压恢复系数、燃烧效率、燃烧室出口的径向温度分布系数等性能参数。计算结果可为环管燃烧室的优化设计提供参考。     

雾化燃料在超声速气流中横向喷射混合的数值模拟

本文用双流体模型对液态雾化燃料在直简形燃烧室中横向喷射问题进行了数值研究,气相用二阶迎负ＴＶＤ格式求解三维全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,液相用二阶预估、校正ＮＮＤ格式求解三维Ｅｕｌｅｒ方程。描述相间相互作用的常微分方程用预估、校正的二阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求解。分别考察了喷口处气相压力、相间速度滑移、液滴直径和燃烧室入口处预置的超声速流向涡对雾化燃料混合效果的影响。结果表明：气相较液相的扩     

液氢/液氧火箭发动机燃烧过程二维数值模拟分析

本文应用数值模拟方法对液氢/液氧火箭发动机燃烧过程进行了计算,给出了发动机燃烧室内燃气参数分布。计算中通过随机喷雾模型确定液氧滴进入燃烧室的初始条件,应用高压蒸发模型计算液氧蒸发速率及Arrhenius公式计算氢氧气相化学反应速率。贴体网格由TTM方法生成,数值方法采用一步隐式预测、两步显式校正的PISO算法,控制方程在非交错网格下离散并且对流项采用了二阶迎风格式。      

隐式耦合求解燃烧室内喷雾两相流场

在任意非结构网格下,采用开源程序对某型涡喷发动机环形燃烧室内的喷雾燃烧过程进行数值模拟.首先选用k-ε双方程模型、EBU(eddy break up)-Arrhenius紊流燃烧模型、颗粒随机轨道模型等模拟流场中的物理过程.然后用隐式耦合算法求解气相流动、拉格朗日方法跟踪液滴相在流动区域的运动,并考虑气相与液滴相的完全双向耦合作用.仿真实验得到的 燃烧室内温度分布、马赫数分布、速度分布、压力分布均较合理,温度分布与文献中的基本一致,这表明隐式耦合方法能够较好地模拟燃烧室内流场,反映其内部流动特点.      

过冷大水滴动力学特性对结冰影响数值研究

对比常规水滴,过冷大水滴结冰对飞机气动性能和飞行安全具有更大的危害性。在常规水滴结冰过程数值模拟基础上,针对过冷大水滴条件下各种动力学特性及其对成冰过程的影响进行了数值研究。采用结构化网格和中心有限体积法求解 N-S 方程获得空气流场,用拉格朗日法求解水滴流场,基于改进的 Messinger 热力学模型完成结冰过程模拟;对于过冷大水滴,采用泰勒类比理论和水滴碰撞模型,对比分析了水滴动力学特性及其对结冰过程的影响;通过数值模拟获得了水滴 Langmuir D 分布的撞击特性及结冰计算结果,研究了水滴多尺度分布对成冰过程的影响规律。通过一些典型结冰算例的数值模拟,并与参考文献结果及实验数据的对比分析,说明了本文计算方法及计算模型模拟过冷大水滴结冰是可行和正确的。     

胶质微推进液滴加速过程的数值模拟

采用拉格朗日粒子跟踪法计算液滴的运动轨迹,对胶质推进器中带电液滴在加速区的运动过程进行了数值模拟,获得液滴速度和密度的空间分布。使用两种溶剂(TBP3,甲酰胺)作为推进剂,分析了电压、流量变化对喷雾扩散角、液滴轴向速度及对推进器性能的影响。模拟结果表明,随着电压的增加,液滴轴向速度增大,喷雾扩散角减小;随着流量的增大,液滴轴向速度和喷雾扩散角均减小;使用甲酰胺推进剂比使用TP3推进剂得到更高的推力和比冲;模拟得到的推进器性能能良好符合实验结果。     

环形燃烧室两相燃烧数值研究

根据环形燃烧室设计需要,运用计算机模拟技术研究模型环形燃烧室三维有旋液雾燃烧流场,研究两相流动与混合,液雾蒸发对燃烧过程的影响。两相流动模型采用颗粒群轨道模型,气相紊流输运方程求解采用双方程k-ε模型和EBU-Arrhenius紊流燃烧模型,热辐射采用热通量法模型,本文提供计算方法可以用来研究不同性能参数对环形燃烧室两相燃烧特性影响。     

基于双流体模型的超声速两相横向喷射流场研究

对两相超声速横向喷射试验进行了数值模拟.采用双流体两相流模型,在不同直径液滴和不同动压比条件下分别对两相超声速横向喷射流场进行了计算.结果表明:雾化液滴直径和喷射动压比对液体的横向穿透深度影响很大;雾化液滴直径的大小或者动压比的大小对喷口下游流场的影响区域是不同的.数值模拟结果与试验结果符合良好,为进一步的研究打下良好基础.