随机颗粒轨道模型在长尾喷管发动机流场计算中的应用

应用随机轨道模型对长尾喷管发动机长尾段内的轴对称两相流动进行了数值模拟,并对确定轨道和随机轨道获得的计算结果进行了对比分析,研究了10、30、50、80、100μm不同直径颗粒及3种不同颗粒直径分布下颗粒的运动轨迹和颗粒在壁面聚集浓度的区别。研究表明,随着颗粒直径的增加,颗粒的随流性变差,惯性增强,湍流脉动对颗粒的作用减小,不同颗粒轨道模型获得的计算结果区别明显,尤其是喷管扩张段的差异;随机轨道相对于确定轨道获得的结果,流场中几乎不存在无颗粒区域;不同颗粒直径的分布形式对计算结果影响也较明显;颗粒和壁面碰撞的形式主要为小直径颗粒主要受湍流脉动的作用,大直径颗粒主要是由于惯性的作用。     

射流推力矢量喷管中粒子运动轨迹的数值模拟

用颗粒轨道模型对基于激波控制的二维收缩-扩张(2-DCD)喷管中的气固两相流动进行了数值模拟,研究了1～60 μm不同直径颗粒的运动轨迹.研究表明,颗粒尺寸越小,粒子的随流性越好;颗粒尺寸越大,主流中的粒子与收缩段壁面碰撞越剧烈,无粒子区越大,二次流中的粒子对下壁面的碰撞越剧烈.     

潜入式摆动喷管两相内流场数值模拟

气相采用三维薄层近似N-S方程，固相采用拉格朗日坐标系下颗粒轨道模型，计算方法采用显式ENO差分格式和分区多块网格技术，建立了固体火箭发动机内流场模拟软件。对轴对称JPL喷管进行模拟，计算结果与文献进行比较，符合较好。推广到实际三维装药型面时，发现颗粒数目增加导致计算效率降低。采用软件中的气相计算模块对潜入式摆动喷管内流场进行数值模拟，比较了喷管摆动角在发动机喉道附近引起流动变化，研究了非对称流动对喷管内表面压强的影响。     

"长二丙改"上面级CPKM发动机喷管两相粘性流动数值模拟研究

对“长二丙改”火箭上面级CPKM发动机喷管的高空性能进行了预示。气相采用显式MacCormack差分格式、颗粒相采用特征线法,数值求解采用轴对称二维两相粘性流动模型。对CPKM发动机喷管进行数值模拟,所得的发动机平均推力与实测性能数据相比较,相对误差为2．3％。研究表明,所采用的流动模型和数值模拟方法,对高空大面积比喷管性能进行的预示是有效的。     

Al/AP粉末发动机燃料喷射流场的数值研究

为研究Al/AP粉末发动机燃料输送方式对粉末分散及混合效果的影响,利用气相控制方程、湍流模型、欧拉模型,对金属粉末发动机粉末燃料喷射掺混过程进行了详细的数值模拟。数值计算结果表明:固相体积分数和喷射入口有无锥体结构对流场参数具有较大影响。随着Al/AP粉末入口固相质量分数增加,粉末混合区域由燃烧室尾部向燃烧室头部移动。入口处安装锥体有利于粉末在燃烧室内分散,燃烧室头部颗粒相质量分数增加,轴线上平均固相质量分数降低,喷管前粉末局部积聚现象消失。     

固体火箭发动机喷管扩张段粒子冲刷流场分析

某翼柱形药柱固体火箭发动机喷管扩张段出口部位在试验后出现了与药柱翼槽位置相对应的冲刷痕迹,为了研究Al2O3粒子对喷管扩张段的冲刷规律,对喷管型面改进提供依据,对比了不同湍流模型、颗粒轨道模型对形成冲刷痕迹的影响,分析了发动机喷管扩张段两相流场特征,确定了形成冲刷痕迹的粒径范围,判断了冲刷痕迹的形成时间,提出了喷管型面改进方案。结果表明,喷管扩张段的冲刷痕迹形成于发动机工作的15 s时刻之前,主要由药柱后翼燃烧产物中颗粒粒径分布为10~16μm区间的粒子造成,改进后的喷管型面可有效降低粒子对喷管扩张段的冲刷。     

铝镁推进剂固冲发动机两相燃烧数值模拟

针对铝镁推进剂中心进气固冲发动机,湍流模型采用Reynolds应力方程模型,气相燃烧采用涡耗散模型,两相流采用颗粒随机轨道模型,铝颗粒燃烧采用Brooks燃烧模型,对二次燃烧和流动进行了三维、两相和化学反应流场数值模拟,对比和分析了冲压空气与一次燃气无旋、同向及反向旋转3种工况下补燃室燃烧流动特性及燃烧效率,并对不同工...     

基于混合正交法的环喉式膨胀偏流喷管性能影响因素研究

利用混合正交法设计数值模拟实验对环喉式膨胀偏流喷管(ATEDN)的工作过程进行研究，以获得ATEDN在不同工作模态的流场特征并分析入口压力、扩张比和塞锥基底半径三个影响因素共同作用下喷管性能的变化规律及各因素对喷管性能影响的显著性。采用FLUENT商用软件对ATEDN在0～31 km高度范围内的工作状态进行模拟。模拟结果表明，ATEDN在不同工作高度下存在开放-固壁反射、开放-自由边界反射和闭合三种基本工作模态。根据喷管流动经历的不同工作模态特征，将喷管的工作过程进行分类，具有不同工作过程的喷管推力性能存在显著差异。结果对比发现，在所选的参数水平下，塞锥基底半径是影响高度积分平均比冲的最显著因素，其次是扩张比，入口压力影响最小。随着喷管扩张比的增加，高度积分平均比冲降低；随着喷管入口压力和塞锥基底半径增加，高度积分平均比冲增加。     

高过载条件下固体发动机内流场及绝热层冲蚀研究

针对Ф315mm实验发动机和某发动机工作状态及结构特点，进行了发动机燃烧室内三维两相流动数值模拟和内绝热层的炭化冲蚀规律研究。计算中应用了颗粒轨道模型和二阶迎风有限体积方法，对纵、横加速度载荷下的两相流动进行了模拟，分析了纵、横向载荷对两种发动机燃烧室内粒子场和聚集带的影响；应用绝热层炭化冲蚀和两相流粒子热增量模型，分析了实验发动机在多种纵、横向过载作用下的绝热层冲蚀规律。在与Ф315mm实验发动机结果对照后，修正了炭化冲蚀计算所需参数，进一步预示了发动机的三维两相流场和炭化冲蚀率，并在35gn过载下分析了推进剂含铝量对粒子聚集密度和炭化冲蚀的影响。计算结果为发动机故障分析和实验发动机试车结果评定提供了理论基础。     

过氧化氢／煤油发动机推力室气液燃烧数值模拟

采用离散相模拟颗粒随机轨道方法处理两相流动,采用SST后刊湍流模型,利用燃料的单步总反应模型模拟煤油蒸汽燃烧,针对过氧化氢／煤油发动机推力室,进行三维两相燃烧数值模拟。对推力室的总体流场特征及性能进行了初步研究,为开展进一步的数值研究奠定基础。     

微型固体姿控发动机微喷管内气粒两相流动规律的CFD-DSMC研究

微型固体姿控发动机在航天领域具有广泛的应用前景。以基于MEMS技术的微喷管为研究对象,首先通过计算微喷管中的克努森数,得到了微喷管中的气相流动状态;然后,采用CFD-DSMC方法,模拟了微喷管中的气粒两相流动,并研究了颗粒相质量分数和粒径对气相流动的影响。结果表明,在所研究的来流条件下,微喷管中的连续介质假设是成立的;气相与颗粒相间的动量和能量交换,导致气相马赫数降低、温度升高,同时也导致颗粒相速度增加、温度降低;颗粒相质量分数和粒径均能显著影响气相的马赫数和温度。     

粒子流量可调的喷管烧蚀试验方法

发展了一种保持燃气参数不变的情况下能实现粒子流量可调的喷管烧蚀试验方法,并研制了试验装置。该试验方法是将两相流燃气中的一部分粒子收集起来,以减少流经喷管的粒子流量,通过改变收集孔和收敛角的大小来调节粒子流量。采用该方法开展了变粒子流量的喷管烧蚀试验,试验结果验证了该方法是有效的,试验条件下喷管喉部平均线烧蚀率随粒子流量减小而降低。     

粒子尺寸分级的喷管两相流场计算

应用计算流体力学软件FLUENT，对固体火箭发动机喷管内的二维两相粘性湍流进行数值模拟。用拉格朗日方法对粒子的速度、温度、轨迹进行处理。5组级别粒子尺寸的气流被用来实验和计算，计算结果表明，粒子尺寸分布对两相流场有显著影响。     

粉末燃料冲压发动机燃烧室两相流数值模拟

采用颗粒轨道模型对镁粉燃料冲压发动机的两相流场进行了三维数值模拟,目的是为进一步的实验研究提供指导和参考。结合理论性能分析,提出了一种发动机构型,通过数值模拟分析了颗粒粒径、产物相态等因素对该发动机燃烧效率的影响。结果表明,粉末燃料的粒径较小时点火延迟时间短,颗粒在发动机中的滞留时间长,燃料燃烧效率较高。     

两相流环缝塞式喷管理想型面的设计方法

目前固体火箭发动机塞式喷管没有成熟的理论设计方法,设计方法需考虑两相流因素和极限粒子流线的几何约束。在常滞后两相流假设下提出改进的Angelino理想型面法设计两相流环缝塞式喷管,证明了常滞后两相流中的普朗特-迈耶函数关系式,并给出了最终设计公式。用F luent软件计算了改进法设计的喷管型面性能。算例结果表明,与未考虑两相流效应的纯气相理想型面相比,该法设计的型面长度缩短近33%,推力增大约1%,证明了提出常滞后两相流假设的合理性及改进法设计两相流环缝塞式喷管的有效性。     

双燃速星孔药柱长尾喷管发动机三维两相流场数值模拟

采用欧拉-拉格朗日两相方法模拟了带长尾喷管的固体火箭发动机三维两相内流场。在同位网格基础上用有限体积法离散N-S方程,采用不完全LU分解预处理B iCGStab算法求解线性代数方程组,通过可压缩SIMPLE算法求解气相流场。用PSIC方法进行气粒耦合计算,得到了气粒两相的速度、温度等参数的分布以及不同尺寸的A l2O3粒子运动轨迹。计算结果表明,在发动机工作初期气粒两相流场呈现强三维特征,长尾管中后段、喷管收敛段以及燃烧室后封头等部位是烧蚀最严重的部位。     

高温固体颗粒气动阻力系数数值计算

在发动机喷管的气粒两相流场中,金属颗粒阻力的准确计算对整个流场的模拟有着重要作用.而高温颗粒与周围气相的传热会改变气相的性质,进而影响颗粒的气动阻力.采用k-ε湍流模型和强化壁面处理方法,计算了雷诺数Re在0～10 000范围内,颗粒与气相不同温差情况下球形单颗粒的气动阻力系数.颗粒温度升高,颗粒所受的气动阻力随之增加,但是增加的幅度随Re的增加而减小;根据计算结果,拟合了以颗粒雷诺数为自变量的高温颗粒气动阻力系数的计算公式.通过与低Re下的实验数据的对比,验证了所拟合公式的适用性和可信性.