流化进气口数量对粉末燃料供给特性的影响_邓哲.pdf

第 43 卷第 3 期2023 年 6 月弹箭与制导学报Journal of Projectiles，ockets，Missiles and GuidanceVol.43 No.3Jun.2023DOI:10 15892/j cnki djzdxb 2023 03 018收稿日期:2023 02 22基金项目:航空科学基金(20200001056001);江西省微小航空发动机重点实验室开放基金(Ef202180433)资助作者简介:邓哲(1987)，男，副研究员，博士，研究方向:航空宇航科学与技术。通信作者:孙海俊(1987)，男，副教授，博士，研究方向:粉末燃料发动机。引用本文:邓哲，任冠龙，沈文，等 流化进气口数量对粉末燃料供给特性的影响J 弹箭与制导学报，2023，43(3):115-124DENG Zhe，EN Guanlong，SHEN Wen，et al Effects of the number of fluidization inlets on powder fuel supply characteristics J Journal of Projectiles，ockets，Missiles and Guidance，2023，43(3):115-124流化进气口数量对粉末燃料供给特性的影响邓哲1，任冠龙2，沈文2，蔡平2，薛瑞3，孙海俊2(1西北工业大学宁波研究院，浙江 宁波 315103;2南昌航空大学飞行器工程学院江西省微小航空发动机重点实验室，江西 南昌330063;3西安交通大学航天航空学院机械结构强度与振动国家重点实验室，陕西 西安710049)摘要:粉末燃料供给特性是影响粉末发动机工作性能的关键因素之一，为探究粉末燃料的流化和输送过程，通过用户自定义函数实现运动壁面运动，构建了气体 粉末 运动壁面多作用耦合模型，基于欧拉双流体模型，开展了不同进气口数量对粉末燃料供给特性的数值模拟研究，并总结了动壁作用下的稠密气固两相流动规律。研究结果表明:在供粉初期，不同进气口数量下的固相流型演化过程相似，但粉末储箱上部的粉末体积分数随进气口数量增加而降低;在供粉中后期，随进气口数量的增加，粉床高度降低，粉末储箱下部的粉末体积分数增大。当进气口数量为 6 时，相比数量 4 和 2 时，粉末流量波动和压力波动分别降低了约 7 6%和 30 1%，即适当增加进气口数量，有利于提高粉末供给稳定性。关键词:粉末发动机;粉末供给特性;进气口;气固两相流;准二维结构;数值模拟中图分类号:TJ763;V238文献标志码:A文章编码:1673-9728(2023)03-0115-10Effects of the Number of Fluidization Inlets onPowder Fuel Supply CharacteristicsDENG Zhe1，EN Guanlong2，SHEN Wen2，CAI Ping2，XUE ui3，SUN Haijun2(1Ningbo Institute of Northwestern Polytechnical University，Ningbo 315103，Zhejiang，China;2Jiangxi Key Laboratory of Micro Aeroengine，School of Aircraft Engineering，Nanchang Hangkong University，Nanchang 330063，Jiangxi，China;3State Key Laboratory for Strength and Vibration of Mechanical Structures，School of Aerospace Engineering，Xi an Jiaotong University，Xi an 710049，Shaanxi，China)Abstract:Powder fuel supply characteristics were one of the key factors that affected the performance of powder engines Inorder to explore the powder fuel fluidization and conveying process，the moving wall movement was realized through user-de-fined functions，and the gas-powder-moving wall multi-action coupling model was constructed The numerical simulation studyon the powder supply characteristics with different numbers of gas inlets was carried out using the Eulerian two-fluid model，andthe law of dense gas-solid two-phase flow under the action of a moving wall was summarized The results show that the evolutionprocess of the solid-phase flow pattern under different numbers of gas inlets is similar at the beginning of the powder supply，butthe powder volume fraction in the upper part of the powder storage tank decreases as the number of gas inlets increases;in themiddle and late stages of the powder supply，the powder bed height decreases as the number of gas inlets increases，and thepowder volume fraction at the lower part of the powder storage tank increases When the number of gas inlets is six，the powderflow rate fluctuation and pressure fluctuation decrease by about 7.6%and 30.1%，respectively，compared with four and twogas inlets Properly increasing the number of gas inlets is conducive to improving the stability of the powder supplyKeywords:powder engine;powder supply characteristics;gas inlet;gas-solid two-phase flow;quasi-two-dimensional struc-ture;numerical simulation弹 箭 与 制 导 学 报第 43 卷0引言粉末发动机是一种以粉末为燃料，固体、气体或液体作为氧化剂的新型发动机。由于燃料自身的特殊性，使得粉末发动机具有高能量、高密度、高比冲等优点，可以更容易实现燃料的稳定输送和自身流量的灵活调节。此类发动机具有多脉冲起动和可调推力的功能，正在逐步得到人们的重视1。目前，已经开发出了多种类型的粉末发动机，如粉末火箭发动机2-5、粉末冲压发动机6-7 以及 Mg/CO2粉末发动机8-9 等，粉末发动机的发展正趋于多元化。在粉末发动机发展过程中，身为其核心部件的推进剂供给系统也在同步发展。推进剂供给装置的作用是将粉末燃料输出，实现粉末的不断供给。然而，粉末燃料自身不具备流体的性质，必须借助相应的载体(气体或液体)使其具备流动的能力，在此载体作用下，粉末经历起动、悬浮以及输运等过程，最终以气固或气液两相流的形式输出。有关粉末发动机中粉末供给过程的研究，最早是由 Fricke 等10 所开展，设计了一种流化床式粉末供给装置，但由于结构较为复杂，也未对该供粉装置进行更多的研究。在此基础上，Meyer 等11 对 Al/Mg 发动机中供粉装置的运动壁面结构进行了简化，并将气孔与运动壁面相连，既能够实现流化气的均匀分配，还有效避免了粉末回流。随后，Foote 等12 设计了一种可使运动壁面运动更平稳的容积式流化床(positive displacement fluidizedbed，PDFB)的粉末推进剂供给装置，但由于气缸结构的存在，增大了发动机的总质量。Miller 等13 对供粉装置做了进一步的改进，通过供气管将流化气体输送至粉末床，既能很好的调节气体流量，也能平衡粉末储箱内压强，不足之处在于软管结构的存在，降低了粉末储箱内粉末的初始装填率。从国内相关进展来看，很多科研团队对供粉装置也进行了深入的研究，使得粉末燃料供给装置得到了进一步的发展。韩超14、申慧君等15 对粉末冲压发动机的燃料供给装置进行了实验和数值仿真研究，认为固气比、粒径、装置结构等参数对电机驱动式供粉特性和发动机燃烧性能会产生较大的影响。张虎等多个团队对结构更为简便的气压驱动式粉末燃料供给装置开展了大量研究16-21，也验证了气压驱动式粉末供给装置使用在部分种类粉末发动机上的可行性。由此表明，有关粉末燃料供给装置的探究热度逐渐递增，供粉装置依旧存在很大的研究前景和应用空间。在粉末发动机技术中，最难也是最亟需解决的问题是如何实现供粉装置中粉末燃料的稳定输送。粉末燃料从静止到完成整个输运过程，其间涉及十分复杂的稠密气固两相流动，气固两相流动状态以及固相流型变化能够反映出粉末流态化性能，并直接影响粉末供给稳定性，可看出粉末流化特性是影响供粉装置性能的关键。任冠龙等22 研究了准二维粉末储箱构型下的粉末流态化和输送特性，总结了进气流量变化对粉末供给特性的影响规律。霍东兴等23 通过改变气体温度、压强、粉末粒径等参数，对粉末沉降速度和最小流化速度进行了探究，并总结了相应的变化规律。苏晓杰等24 对爆轰发动机中粉末燃料供给特性开展了研究，结果表明在预压压力条件下，粉末流化性能较好，且粉末流化能力与活塞速度成正相关。孙海俊等对处于高压状态下粉末起动特性和粉末供给特性开展了研究25-27，认为高压环境下，粉末能够形成稳定的气固分界型面，极大提高了粉末输送稳定性。en 等28 探究了不同粉末储箱内初始工作压力对粉末流化和输送特性的影响，认为高压环境下的粉末输送状态更为稳定。Liang 等29 对高压下的密相气力输送过程进行了数值研究和分析，并着重探究了流态化数和粉末粒径对送料稳定性的影响。从目前研究来看，几乎都是探究流化参数(流化气量、压强环境、粉末燃料粒径等)对粉末流化和输送特性的影响，而有关不同进气结构或者进气的布置方式等对粉末供给过程产生的影响的研究较少，再加上粉末在重力作用下极易沉降，气体入口位置对粉末运动行为的影响并不明确，尤其是在气体 粉末 运动壁面三者耦合作用机制下的粉末流态化过程尤为模糊，且有关此方面的研究也少之甚少。另外，有关粉末供给特性的研究主要集中在实验方面，所表现出的固相宏观现象不足以充分理解详细的气固运动细节，借助数值仿真手段能够有效完善对粉末流化和输送过程的分析内容，这对于深入了解气固两相流动行为和固相流型的演化过程十分重要，因此有必要开展相关方面的研究。基于以上研究，文中应用数值计算方法，在构建的气体 粉末 运动壁面相互耦合作用模型下，采用欧拉双流体来描述稠密气固两相流动行为，着重讨论不同进气口数量对粉末流化和输送特性的影响，以揭示详细的气固两相流动细节，为后续供粉装置的改进和工程应用提供一定的理论基础和设计支撑。1几何模型与数值计算方法1 1粉末储箱构型在前期工作中，已经验证了准二维粉末储箱构型611第 3 期邓哲等:流化进气口数量对粉末燃料供给特性的影响的可行性22，相关文献也表明准二维结构在减少计算量的前提下，所表现出的流动特性与三维结构下的流动状态相似30-31。因此继续选用准二维粉末储箱构型，图 1 为准二维粉末储箱构型。不同粉末储箱构型的差异主要是进气口位置和进气口数量，设计的核心思路是在保证入口表面积不变的条件下，改变进气口数量，即在收敛段上分别均匀布置若干个气体入口。表 1 给出了详细的入口数量与工况的对应关系。表 1计算工况Table 1Simulation casesCases123456Total number ofair intakes24681012另外，图 1 所示的粉末储箱水平放置，其结构主要由运动壁面、流化腔以及输送管道 3 部分组成。其中，运动壁面始终沿着 X 轴正方向匀速运动。在粉末储箱出口附近设计喉道结构，以便控制出口粉末流量。具体的结构参数如图 1 所示。图 1准二维粉末储箱构型Fig 1Powder storage configuration1 2数值计算方法和边界条件采用双流体模型(two-fluid model，TFM)来描述气固两相流动，采用 NG k-湍流模型和颗粒动力学模型(kinetic theory of granular flow，KTGF)32-33 分别描述气相运动和颗粒行为，气固相间作用由 Gidaspow曳力模型34 来描述，PC-SIMPLE(phase coupled SIM-PLE)算法用于求解压力 速度耦合，动量、湍流和体积分数等方程中的对流项采用二阶迎风差分格式。流化气体和粉末燃料分别采用空气和铝粉末，气体入口和气固出口分别选用质量流量入口和压力出口边界条件。通过用户自定义函数(user definedfunction，UDF)实现运动壁面的匀速运动。另外，入口流化气体流量由理论粉末流量的百分比给出，而理论粉末流量为17:mpt=packingvmwAmw=pvmwAmw(1)式中:mpt为理论粉末流量;packing为初始粉末装填密度;vmw为运动壁面速度;Amw为运动壁面面积;p为粉末密度;为初始粉末装填率。表 2粉末属性和操作条件Table 2Powder properties and operating conditionsParameterValuePowder paiticle size/mm0 02Inlet gas mass flow rate/(g/s)0 054Initial powder filling rate0 55Max powder filling rate0 63Growitational acceleration/(g/(m/s2)9 81Powder viscosity/(g/(Pas)1 72 105Powder density/(kg/m3)2 719Operating pressure/MPa0 101Moving wall surface velocity/(m/s)0 05Time step/ms0 5Iteration step202网格无关性验证采用商业软件 CFD-ICEM 对粉末储箱构型进行结构化六面体网格划分，以 Case 1 为例，六面体网格划分结果如图 2 所示。采用 1 2 节数值计算方法和边界条件对粉末供给过程进行数值模拟，并选取网格数为 3 万、10 万、20 万以及 40 万的网格计算域以验证网格无关性。图 3 为不同网格数下 X=0 090 m 处粉末轴向速度对比结果。由图可知，网格数为 3 万的计算结果与其余数量网格相比时，差异较大。而网格数为 10 万的计算结果与 20 万和 40 万相比时，粉末速度曲线具有很好的吻合形，但在 Y=0 055 m 附近，计算结果差异较大，综合考虑计算精度和计算成本，最终选择 20 万左右网格用于全文的数值计算。711弹 箭 与 制 导 学 报第 43 卷图 2粉末储箱的六面体网格Fig 2Structured hexahedral meshes of the powder storage tank图 3网格无关性验证结果对比Fig 3Comparison of grid independence verification results3计算结果与分析3 1计算模型验证为验证所选计算模型的准确性，将构型 1 的计算结果与实验结果26 进行对比，其结果如图 4 所示。图 4计算模型验证结果对比Fig 4Comparison of calculation model validation results在粉末供给过程中，粉末储箱内压降率是研究粉末流动行为的重要参数，其含义是单位时间内粉末储箱内压力变化量与时间的比值，其计算公式为:p=p1 p0t(2)式中:p1为初始总压;p0为终止压强;t 为时间间隔。由图 4 可知，数值计算结果与实验结果具有较好的一致性，最大相对误差在 13%以内，由此认为所选择的计算方法和边界条件用于模拟粉末的供给过程是合理且准确的。3 2粉末供给过程粉末流态化过程涉及复杂的气固流动，其流态化性能对粉末输出有直接影响。图 5 为 Z 向中心截面(Z=0 003 m)处不同时刻下粉末流化和供给过程。由图可知，当流化时刻 t1=0 01 s 时，不同工况下的气相和固相作用范围均呈现上下对称分布。在流化气体卷吸和运动壁面推动的共同作用下，粉末开始运动并不断输出粉末储箱。当流化时刻 t2=0 10 s 时，上一时刻的气固对称分布现象消失，小气泡逐渐演变成了大气泡，呈现出收敛段上部处气相范围大于下部的现象，这主要是因为粉末在自身重力和运动壁面的推动作用下，不断挤压下部进气口，导致下部气体流动受限，而上部气相范围受到粉末重力的影响较弱，因此上部处气相范围较大。另外，从图中还能看出，入口数量的改变对锥形段上部的气相分布范围影响不大，而下部处气相范围随着入口数量的增加而增大，形成的气固分界面结构稳定性随入口数量的增加而降低，尤其在 Case 2 Case 4 中，出现了大范围粉末低浓度区域。分析原因如下:当只存在两个进气口时(即上下各一个)，上部气体流入粉末储箱的运动方向与下部处气体相反，上部处气体不断向粉末储箱内部流动，下部处气体在粉末破挤压作用下大多会流出粉末储箱，此时气体与气体之间的湍动作用较弱。然而，随着入口数量的增多，从不同位置入口处流入的气体会相互影响和冲击，导致气相湍流脉动增大，再加上粉床不断挤压下部处粉末，在强烈的气体 气体、气体 固体以及固体 固体相互作用下，气固分界面结构以及粉末分布受到较大的影响。另外，由于此流化时刻还处于粉末储箱的启动初期，气体卷吸输出的粉末量与运动壁面推动的粉末量还未达到很好的匹配关系，导致粉末缺额现象的发生，形成了粉末低浓度区域。当流化时刻 t3=0 40 s 时，首先粉末储箱上部的气相范围不断增大，演变成了大气腔结构，而且不同进气口数量下的气固分界型面结构不稳定，整体呈现波浪形分布，且轴向范围 0 10 m，0 15 m 内的气固分界面厚度高于轴向范围 0 05 m，0 10 m，主要原因是流化气体卷吸输出的粉末量与运动壁面推动的粉末量未达到很好的一致性，即流化气体卷吸输出粉末所留下的空隙应及时被运动壁面推动的粉末填充，粉末输出和粉末装填过程不断博弈，完成整个粉末供给过程。然而，图中大量气腔结构的存在表明流化气811第 3 期邓哲等:流化进气口数量对粉末燃料供给特性的影响体卷吸输出的粉末量大于运动壁面推动的粉末量，粉末不能及时填补空隙，使得气相不断扩展。此时，在气体湍动作用下，向运动壁面方向运动的气体卷吸形成的回流旋涡与从上部入口流入粉末储箱的气体相互作用和冲击，导致粉末被卷吸扬起和悬浮，形成了较厚的气固分界面。另外，还能发现随着进气口数量的增加，锥形段下部壁面处会形成类似条状的粉末分布，这主要是因为入口与入口之间存在间隙，间隙中的粉末在重力和气体力的影响下受力平衡，形成上述粉末分布现象。图 5粉末流化和输送过程Fig 5The powder fluidization and conveying process当流化时刻 t4=1 00 s 时，不同工况下粉末储箱内的气固分界面近乎水平分布，相比上一时刻，气固分界面厚度明显增加，这主要归因于气体的提升作用，粉末储箱上部的粉末被气体完全卷吸输出，出现明显的上下气固分层现象。此时，粉床上部的粉末受到气体力和自身重力的影响，悬浮在气流中，粉末与911弹 箭 与 制 导 学 报第 43 卷粉末间距增大，粉末浓度降低，最终表现在气固分界面的厚度上。当流化时刻 t5=1 60 s 和 t6=2 90 s 时，不同工况下粉末储箱内的气固分界面水平分布，且气固分界面厚度减小，表明此时粉末填补量与粉末输出量处于动态平衡，流动状态趋于稳定，由此形成了结构相对稳定的气固分界型面。此外，相比前几个流化时刻，当前时刻下存在一个明显的固相现象，即随着进气口数量的增加，在中心轴线上无法形成稳定的气固两相通道，尤其是在 Case 6 中，气体流动成为了气固两相流动的主流，原因在于随着供粉过程的不断进行，粉末储箱内残余的粉末量降低，导致粉床对气相的挤压作用减弱。当入口数量为 2 时，下部进气口对粉末储箱底部的粉末影响较小，此时粉末流化性能较差，在运动壁面的推动作用下不断堆积在收敛段下部。随着进气口数量的增加，堆积在收敛段下部的粉末会被气体夹带而流出粉末储箱，但由于粉末无法及时填补气体输出粉末留下的间隙，导致气相流动占据主导地位，再加上粉末的持续输出，多个进气口粉末储箱内的粉床高度也较低。3 3粉末体积分数变化特性由上述粉末流化和输送过程可知，进气口数量的变化显著改变了固相流型的演变行为。为此，针对粉末空间分布采取量化处理，这对于深入理解粉末运动行为十分重要。图 6 给出了中心截面(Z=0 003 m)与轴向位置 X=0 150 m 处截面交线处的粉末体积分数分布。由图可知，t1 t44 个流化时刻下的粉末体积分数曲线分布趋势不尽相同，需要逐一分析。在流化时刻 t1=0 03 s，不同进气口数量对粉末体积分数的影响主要体现在 Y 轴范围 0 05 m，0 06 m 内，其余 Y 轴范围内的粉末体积分数在 0 56 左右变化，这是因为粉末储箱刚启动，该范围内的粉末在自身重力和运动壁面的推动作用以及粉末间的相互挤压作用下处于自由堆积状态，流化气体入口变化对该位置处粉末运动产生的影响微乎其微。在粉末储箱顶部 Y轴范围 0 05 m，0 06 m 内，不同工况的粉末体积分数下降，且下降速率随入口数量的增多而加快，其原因在于随着进气口数量的增多，从收敛段尾部入口流入的气体对粉末有冲击和卷吸作用，该作用随进气口数量增多而增强，导致粉末间距增大，粉末浓度降低。而在 Case 1 中，气体入口位于收敛段的中心位置，气体对收敛段尾部处粉末影响较小，故粉末体积分数无明显变化。图 6不同时刻下粉末体积分数分布Fig 6The distribution of powder volume fraction atdifferent moments在流化时刻 t2=0 17 s，不同工况的粉末体积分数在不同 Y 轴位置发生了较大的变化。在 Y 轴范围 0 m，0 04 m 内，不同进气口数量下粉末体积分数无明显变化，均在057 左右稳定变化。在Y 轴范围 004 m，0 06 m 内，不同工况下的粉末体积分数都存在两个下降阶段，第 1 个下降阶段处于 Y 轴范围 0 04 m，021第 3 期邓哲等:流化进气口数量对粉末燃料供给特性的影响0.05 m 内，粉末体积分数从 0 56 下降至 0 28 左右。其中，Case 1 的粉末体积分数减小速率高于其余工况，且同一 Y 轴位置，Case 1 的粉末体积分数低于其余工况。其原因为:相比多个气体入口，从单个进气口流入粉末储箱时的气体动能更为集中，气体卷吸粉末的能力较强，此时，粉床中粉末受到的气体提升力与粉末重力博弈，弥散在粉床中，增大了粉末与粉末间距，降低了粉末体积分数。随着进气口数量的增加，上述博弈程度减弱，粉末趋于沉降，造成粉末体积分数增大。第 2 个下降阶段处于 Y 轴范围 0 05 m，0 06 m 内，粉末体积分数从 0 28 下降至 0 14 左右，对应的曲线下降速率低于前一阶段。原因是这一阶段的粉末体积分数是由气体卷吸扬起的粉末作用的，粉末浓度势必低于粉床中粉末浓度。另外，由于 Case1 中只有单个气体入口，相比多个进气口，气体与粉末之间的摩擦耗散减少，其充足的动能可扬起更多的粉末，使得粉末体积分数高于其余工况。在流化时刻 t3=1 10 s 和 t4=2 90 s，不同进气口数量下的粉末体积分数无明显差异，但相比前两个时刻，该两时刻下的粉末体积分数曲线下降拐点提前，拐点分别处于 Y=0 035 m 和 Y=0 030 m 附近，这主要是因为该两时刻下粉末储箱上部几乎全部被气腔结构充满，只有零星粉末弥散和悬浮在气流中，导致粉末浓度大幅度降低。另外，从图中还能发现，在拐点之前，Case 1 的粉末体积分数大于其余工况，表明在粉末流化过程的中后期，多个进气口结构的存在对压实状态下的粉末扰动作用较大。3 4粉末输出特性粉末输出特性可由出口粉末流量来表征。图 7为不同进气口数量下粉末流量分布。由图可见，不同工况下粉末流量均会出现大范围波动和相对稳定变化阶段。在流化时段 0 1 5 s 内，不同工况下的粉末流量出现大范围波动，且持续时间较长，表明有大量粉末喷出粉末储箱，十分不利于粉末的稳定输送。另外，随着进气口数量的增加，粉末流量波动幅度降低，对应的流量曲线更加平滑，表明进气口数量的增加能够削弱气固相间作用，有利于改善初始阶段粉末的不稳定喷动现象。除此之外，还能发现不同工况在流化时刻 t=1 0 s 左右都存在一个较小的粉末流量峰值，这主要归因于气固分界面结构的改变，在此时刻之前，气固分界面呈现波浪形分布，在此时刻之后，气固分界面呈现水平分布，此时，流化气体卷吸输出的粉末量较多，再加上运动壁面的推动作用，最终导致粉末输出时再次形成一流量峰值。随后，在流化时段1 5 2 9 s 内，粉末流量曲线虽有波动，但整体相对平稳，其中 Case 1 的流量曲线呈现高频振荡。另外，在稳定输送阶段，各工况的瞬时粉末流量曲线与理论粉末流量曲线间存在较大差距，且均低于理论粉末流量，这主要是因为在 0 1 3 s 时间段内，大量粉末喷出粉末储箱，导致粉末储箱内粉末量减少，再加上运动壁面推送的粉末量远小于流化气体卷吸输出的粉末量，所以导致稳定阶段输出的粉末流量低于理论粉末流量。121弹 箭 与 制 导 学 报第 43 卷图 7不同工况下粉末流量分布Fig 7The distribution of powder flow rate in different cases为了定量分析不同进气口数量下粉末流量的波动程度，选用均方差(standard deviation，STD)来表征瞬时粉末流量偏离平均值的程度26。图 8 展示了不同进气口数量下瞬时粉末流量的均方差分布。可看出，Case 3 的均方差最小，相比最大均方差值即 Case2，粉末流量波动降低了约 7 6%，表明当在锥形段上下各布置 3 个进气口时，粉末流量波动较小。图 8不同工况的粉末流量均方差对比Fig 8Comparison of powder flow rate standard deviation indifferent cases总的来说，在整个粉末流化和供给过程中，出口粉末流量存在剧烈波动，这种流量现象符合之前的实验结果17，25。通过对粉末流化过程和粉末流量进行分析，可看出进气口数量对出口粉末流量的影响较大，尤其是当进气口数量为 4 时，粉末流量波动最大，这就要求在实际工程领域中设计合适的进气口数量来提高粉末输送稳定性。虽然在采用准二维粉末储箱构型模拟粉末供给过程时与三维流动区域相比存在误差，但在计算成本降低的基础上，所表现出的气固流动规律亦能很好的为后续二维粉末储箱的应用以及三维粉末储箱的设计和发展提供思路。3 5粉末储箱内压力变化特性图 9 展示了不同进气口数量下粉末储箱内压力分布对比结果。整体来看，不同工况下压力变化趋势相似，在经历两个压力峰值之后开始小幅度稳定变化，但在压力波动幅度上，不同进气口数量下略微不同，尤其是 Case 1 的压力波动幅度较大，体现在整个压力变化过程中。这主要是因为粉末以气固两相流的形式流出粉末储箱，Case 1 中的粉末流量波动幅度较大，直接导致气体输出过程不稳定，此时，气体流量波动也随之增大，恒定的气体输入量与变化的气体输出量不断博弈，导致粉末储箱内压力波动幅度较大。图 9不同工况下粉末储箱内压力分布Fig 9The distribution of pressure in the powder storagetank in different cases图 10 为不同进气口数量下粉末储箱内压力的均方差分布。可看出，均方差的最大值和最小值分别对应于 Case 1 和 Case 3，相比之下，Case 3 的压力均方差降低了约 30 1%，表明当进气口数量为 6 时，粉末储箱内压力波动幅度最小。图 10不同工况下粉末储箱内压力均方差分布Fig 10The standard deviation distribution of pressure in thepowder storage tank in different cases221第 3 期邓哲等:流化进气口数量对粉末燃料供给特性的影响4结论针对准二维粉末储箱构型，在构建的气体 粉末 运动壁面多物理耦合作用模型上，采用 TFM 探究了不同进气口数量对粉末燃料供给特性的影响，主要结论如下:1)在不同进气口数量下，气固分界面从波浪形分布逐渐演变成水平分布，并最终形成明显的上下气固分层现象。当粉末供给过程接近尾声时，粉床高度随进气口数量增加而降低。2)在粉末供给初期，粉末储箱上部的粉末体积分数随进气口数量增加而降低。在粉末供给后期，粉末储箱下部的粉末体积分数随进气口数量增加而增大。3)当进气口数量为 4 时，粉末流量均方差最大。当进气口数量为 6 时，粉末流量均方差最小。相比前者，后者的粉末流量波动降低了约 7 6%。4)当进气口数量从 2 增加至 6 时，粉末储箱内压力波动降低了约 30 1%。通过适当增加进气口数量的方式，能够改善粉末流化性能，提高粉末输运稳定性。参考文献:1 董新刚，霍东兴，张强，等 粉末发动机技术研究现状及展望 J 固体火箭技术，2021，44(2):166 178DONG X G，HUO D X，ZHANG Q，et al esearch prog-resses and prospect of powdered fuel engine technology J Journal of Solid ocket Technology，2021，44(2):166 178 2 LOFTUS H J，MASHALL D，MONTANINO L N Powderrocket feasibility evaluation M New York:Bell Aerosys-tems Co，1973 3 李悦，胡春波，胡加明，等 粉末火箭发动机研究进展 J 推进技术，2018，39(8):1681 1695LI Y，HU C B，HU J M，et al Progress of powder rocketengine technologyJ Journal of Propulsion Technology，2018，39(8):1681 1695 4 武冠杰，任全彬，李超，等 Al/AP 粉末火箭发动机推力调节性能试验J 固体火箭技术，2018，41(5):543 548WU G J，EN Q B，LI C，et al Experimental study onthrust regulation characteristics of aluminum/ammoniumperchlorate powder rocket moto J Journal of Solid ocketTechnology，2018，41(5):543 548 5 李超，胡春波，武冠杰，等 喷注压降对 Al/AP 粉末火箭发动机工作特性影响研究 J 推进技术，2018，39(12):2873 2880LI C，HU C B，WU G J，et al Investigation on effects ofinjection pressure drop on Al/AP powder rocket motor oper-ation characteristics J Journal of Propulsion Technology，2018，39(12):2873 2880 6 李慧强，徐旭，朱清波，等 以粉末燃料冲压发动机为动力的火星巡航飞行器方案初步研究 J 载人航天，2021，27(3):334 338Li H Q，XU X，ZHU Q B，et al Study on preliminaryscheme of mars cruise vehicle powdered by powder fuelramjetJ Manned Spaceflight，2021，27(3):334 338 7 任蒙飞，席文雄，罗世彬，等 粉末燃料冲压发动机头部组织掺混流动数值模拟J 火箭推进，2020，46(5):35 41EN M F，XI W X，LUO S B，et al Numerical simulationof mixing flow in the head of powder fuel ramjet J Jour-nal of ocket Propulsion，2020，46(5):35 41 8 朱小飞，胡春波，胡加明，等 氧化剂气量配比对 Mg/CO2 发动机性能影响实验J 航空动力学报，2019，34(2):479 485ZHU X F，HU C B，HU J M，et al Experiment on influ-ence of oxygen gas intake plan on performance of Mg/CO2engineJ Journal of Aerospace Power，2019，34(2):479 485 9 WEI G，HU C B，WU F Z，et al Heat-transfer charac-teristics of CO2 boiling flow in the regenerative coolingchannel of an Mg/CO2 powder rocket engine for Mars mis-sions J Acta Astronautica，2021，189:43 54 10 FICKE H D，BU J W，SOBIENIAK M G Fluidizedpowders-A new approach to storable missile fuelsC Denver-12th JANNAF Liquid Propulsion Meeting，1970 11 MEYE M L Powdered Aluminum and oxygen rocket pro-pellants:Subscale combustion experiments:NASA-TM-106439 Washington:NASA，1993 12 FOOTE J P，LITCHFOD J Powdered magnesium:Carbon dioxide combustion f