ADN基推力器热回浸对毛细管微尺度流动特性影响的数值模拟研究

小推力ADN基推力器在工作中,其毛细管内的推进剂容易在壁面传热作用下发生相变,进而影响推力器的正常工作。为了深入理解推进剂在毛细管内的相变和流动特性,采用三维数值模拟方法对毛细管微尺度流动和相变特性进行计算。计算考虑了毛细管与喷注器内的流固耦合传热、推进剂相变过程。推进剂相变采用Lee模型,气液体积分数的求解和气液界面捕捉重构采用VOF-CSF方法。本文首次采用VOF模型耦合Lee相变模型计算了ADN基推进剂在毛细管内的相变过程,并结合气泡空间分布仿真结果得到了质量流量和热回浸温度对流动特性的影响规律。计算结果显示,受到ADN基推进剂的冷却作用,毛细管内壁面温度要略低于外表面。毛细管内的气泡形成于弯管处,其体积沿着毛细管轴向下游逐渐增加,采取散热措施可减少并推迟气泡的形成。随着ADN基推进剂质量流量的降低或下游热回浸温度的增加,毛细管内的气泡均显著增加,且形成区域更接近上游区域。当热回浸温度从800K增加至1100K时,毛细管内的气泡体积从25.6mm~3增加至58.7mm~3。     

S型进气道沿程结冰参数变化的数值模拟

采用欧拉-欧拉法,对低来流马赫数和大气结冰条件下S型进气道沿程的空气-过冷水滴两相流场进行了数值计算,获得了进气道沿程和出口面上结冰参数的分布和变化特征.计算结果表明:当来流马赫数从0.1增加到0.5时,进气道入口处附近相对于来流的气动降温从6K增大到17K左右,进气道出口面上的温度畸变指数逐渐增大,液态水含量畸变指数先急剧增大后缓慢减小,在研究范围内来流马赫数为0.2时液态水含量畸变值达到最大.      

高速来流发动机射流预冷及压气机湿压缩特性一体化分析方法

为了研究高空高速来流条件下发动机射流预冷及压气机湿压缩特性,本文针对发动机进气道部分和压气机部分的气液两相流动过程分别建立相应的计算模型,发展了一体化分析方法。对于进气道喷水后的气液两相流动,采用基于微元段思想的气动效应-蒸发冷却修正两步法建立一维计算模型;对于压气机湿压缩过程,采用基于CFD技术的欧拉-拉格朗日方法对气液两相流动进行建模。利用所建立的一体化分析方法,考虑真实高空高速来流条件,同时保证发动机在喷水前后具有等换算转速调节规律,对发动机射流预冷及压气机湿压缩特性进行分析,结果表明：射流预冷能够显著降低高马赫数来流下的压气机进气温度,使得相同最高物理转速下压气机换算转速提高,使压气机具有更好的气动性能;若液滴在进入压气机前蒸发完全,混入工质中的过量水蒸气会降低压气机的无量纲转速,使得同换算转速下压气机压比较不喷水时降低;小粒径液滴更快的蒸发作用有助于对高温气流的快速降温,当液滴粒径3μm和喷湿量为10%时,压气机进气温度较喷水前降低247.13K,同换算转速下物理转速仅为不喷水时的76%。     

气液两相流中气泡速度的电磁互相关测量

气液两相流作为一种工业领域中普遍存在的流态,在测量时需考虑气相和液相间的速度差。相较于流动规律平稳的气相而言,研究不稳定的气泡流动状态具有更高的技术难度。本文针对空气-水为工质的两相流中一项重要参数——气泡速度进行测量,引入基于互相关运算的电磁检测方法,通过提高系统激励频率得到较大的检测信号。实验设计了高频(1 MHz)电磁检测系统,选用垂直上升管道,并在管道的2个平行截面分别安装电磁传感器,每组传感器均包含一个激励线圈及一个接收线圈,对2组接收线圈上的相位信号进行互相关运算,从而求得气泡的速度。实验中,对3种不同速度气泡采集到的信号进行了对比分析,相对误差控制在10%以内,提供了一种完全非接触非侵入的测量气泡速度的方法。此方法可以进行后续补充完善,用于其他工业场合,如金属液中气泡参数的测量等。     

气液两相流在180°弯头上T形管处相分离的实验研究

利用180°弯头上的T形管考察气液两相流的相分离行为,以加深对两相流在T形管处相分离机理的理解。以空气和水为两相流工作介质,在鼓泡流和环状流条件下,入口主管分别以垂直向上和垂直向下2种方式放置,通过改变气液入口流型和流速,考察了180°弯头上T形管处合力的大小和流型的变化,测定了该处的相分离数据。实验结果表明:入口主管垂直向上时,在鼓泡流条件下,气体主要受液体的浮力作用,液相主要受重力作用,侧支管以气相采出占优,增加气液流速对气相采出不利;而在环状流条件下,液体离心力占主导地位,侧支管以液相采出占优为主,气液两相流速增加对相分离有利;主管垂直向下时,在环状流条件下,以液体向下的离心力和重力占主导地位,侧支管中液相采出占优,增大气体流速或者液体流速,不利于液相在侧支管的采出。利用T形管处的合力大小、入口流型和两相入口动量能有效解释相分离结果的变化规律。     

动载对管内汽水两相流流阻、空隙率和传热的影响

利用旋转平台对动载作用下单相水和沸腾两相流流动及传热特性进行了实验研究.通过改变旋转速度、入口温度、液体流量等参数,获得了静止和旋转2种状态下的单相水及沸腾两相流实验数据.对旋转状态下的管内两相流型进行可视化拍摄,得到了动载作用下的两相流型图.结果表明,过载方向与流动方向相反时,过载越大,管内压力和流阻越大,液相流量、空隙率和流体得热量越小.可见,动载阻碍流体向前流动,强化了管外散热,削弱了管内流体的相变换热能力.     

动载对管内沸腾两相流传热特性的影响

对旋转状态下的单相水及加热单相水至沸腾的两相流进行了对比实验.通过对热耗散、流型、流阻及传热的分析与计算,研究了动载对流体流动与传热的影响.结果表明:过载强烈影响两相流动(尤其是流阻)和热传递方式,增大了流体的流阻和热损失,并削弱了流体的得热能力.动载引起的冲击交混流改变了管壁温度的均匀性和变化规律,即随着过载的增大,在一定范围内,管壁底部温度降低,顶部温度升高.      

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.      

回流燃烧室燃烧过程的三维数值模拟

在三维任意曲线坐标系下数值模拟回流燃烧室火焰筒内两相燃烧过程,采用RNG k-ε模型模拟紊流粘性,EBU-Arrhenius模型模拟燃烧反应速率、离散坐标模型以及六通量模型考虑辐射传热,液相采用颗粒轨道模型,气相采用SIMPLE算法求解,并用PSIC算法考虑气液两相之间的相互作用的影响,计算得到燃烧室内速度、温度等各气流参数分布.通过将计算与实验结果对比表明,计算方法可靠,离散坐标模型优于六通量辐射模型,更适用于模拟火焰筒内两相燃烧流场.      

N2O/HTPB固液火箭发动机喷管两相流计算

利用二维轴对称N-S方程对选用氧化亚氮/丁羟基燃料推进剂的固液混合火箭发动机的喷管两相流进行了计算.计算采用MacCormack时间推进预报校正二步格式,采用了Baldwin-Lomax代数湍流模型和两相平衡流模型.计算了三种氧燃比下4个不同喷管的喷管流场参数,并计算了喷管性能,通过比较两相流和气相流的计算结果,分析了不同氧燃比和喷管形状对喷管性能的影响,认为固液火箭发动机的性能主要受氧燃比的影响,为固液混合火箭发动机的设计提供了依据.