水下爆轰燃气泡形态与激波传播过程研究*(爆震专刊)

针对脉冲爆轰发动机在水下工作过程中形成的燃气射流问题,搭建了水下爆轰燃气实验系统,研究了第一个爆轰循环在水下的燃气泡发展变化过程。建立了基于气液两相双流体模型的脉冲爆轰发动机水下喷射模型,采用时-空守恒元和求解元方法,模拟了爆轰波退化为激波在水中的传播及衰减过程。研究结果表明：燃气泡前期受外界水环境阻滞作用呈现“豌豆状”形态,充分发展阶段气液交界面逐渐失稳,在达到最大尺寸后开始收缩并在中心轴线位置出现凸出的射流；水下爆轰燃气射流发展过程中同时存在脱离燃气泡的水中前导激波和管口燃气泡内的高压区两部分,水中前导激波在传播过程中压力迅速衰减至常压量级,而管口燃气泡内则一直保持较高压力；中心轴线区域气液交界面处反射激波的回传使管口附近出现回击现象,并导致前导激波波阵面上压力峰值逐渐出现在30°方向上。     

狭缝内乙烯/氧气预混气体爆轰几何极限的实验

为研究狭缝内爆轰波传播极限,实验得到了不同初始压力(0.004~0.04MPa)下化学恰当比的乙烯/氧气预混气体在狭缝高度为1.0~4.0mm狭缝内的爆轰性能.采用烟膜板记录爆轰波运行轨迹,高速摄影捕捉火焰面.结果表明:狭缝高度越小,爆轰极限对应的临界初始压力越高.对于近极限条件内不稳定爆轰传播模式,包括"结巴"爆轰和驰振爆轰,其对应的初始压力范围随着狭缝高度的降低而变宽.考虑初始和边界条件,将水力直径与胞格宽度之比作为合理的爆轰敏感性参数描述预混气体爆轰特性,得出在不同狭缝通道内爆轰极限范围为,即该比值在0.326~0.403之间.      

壁面的表面粗糙度引起的Stokes层亚临界不稳定性

用数值模拟的方法研究了二维壁面的表面粗糙度下Stokes层的非线性亚临界不稳定性问题.发现当粗糙度高度极小时,响应系数曲线与线性情况就会产生较大偏离.随着粗糙度高度的增加,扰动1阶谱会出现亚谐波的成分,粗糙度高度的进一步增加使扰动1阶谱进入混乱阶段,显示出亚临界失稳的过程.根据粗糙度高度与扰动1阶谱演化的特征关系,定义了临界粗糙度高度,并给出临界粗糙度高度与雷诺数的关系曲线.结果表明:临界粗糙度高度随着雷诺数增大而减小.雷诺数为300左右时,微米量级的粗糙度高度就可能引起Stokes层的亚临界失稳,发生转捩,由此也可以给出实验中观测到的转捩通常都发生在雷诺数为300附近的原因.      

直流喷射首次破碎的形变过程研究

为深入了解直流喷射过程中射流柱的细部结构以及其脱落过程,利用LES结合VOF的方法,对静止大气中垂直射流现象进行模拟。模拟得到的液相喷雾结构和试验结果能够很好地吻合。通过计算观察得到射流柱进入到大气中后由于Rayleigh-Taylor(RT)不稳定性迅速形成伞状的头部,头部的边缘在不断变薄失稳的过程中脱落形成液带,液带脱落产生的液滴具有等距性,间距为0.106mm,并在伞状边缘的下方形成一个气涡,气涡与上游的射流柱相互作用,促进射流柱表面的脱落。射流柱表面由于Kelvin-Helmholtz(K-H)表面波的作用呈现鱼鳞状的结构,并导致射流柱整体断裂、破碎,其表面波波长由初始的0.26mm迅速增长到0.78mm。     

狭缝内爆轰波传播模式的实验研究

为研究狭缝内爆轰波起爆及各种传播模式的特点,实验得到了不同初始压力下(p0=10~60k Pa)化学当量比的乙烯/氧气预混气体在高度1.0mm,宽度20mm狭缝内的爆轰性能。采用烟膜片记录爆轰波运行轨迹,高速摄影捕捉火焰面。结果表明:随着初始压力的降低,实验依次得到稳定爆轰模式、"结巴"爆轰模式、驰振爆轰模式和低速稳定传播模式,中间两者为不稳定模式。不稳定爆轰传播模式具有强烈的速度震荡,在一个周期内,烟膜板中周期性的出现过驱爆轰的胞格结构,随后"结巴"爆轰以单头形式与侧壁面发生碰撞并向前传播,驰振爆轰短暂熄灭。低速稳定传播速度为45%的C-J速度,烟膜板中无胞格结构出现。狭缝内爆轰波速度亏损值为10%~15%,高于大尺寸通道中的速度亏损。     

文氏管出口张角对旋流杯流动不稳定性的影响

为探究旋流杯内流动不稳定现象,针对带不同出口张角文氏管的旋流杯冷态流场进行了大涡模拟,并结合PIV试验结果进行了验证。研究结果表明,文氏管出口张角对时均流场的影响主要在旋流杯内部区域,对下游恢复区域基本无影响;旋流杯瞬态流场中存在旋进涡核（PVC）和小涡旋结构,出口张角对PVC的运动频率和形态基本无影响,对小涡旋结构影响较大,相同Q值下,56°张角旋流杯方案的小涡旋结构更丰富;经本征正交分解方法（POD）分析,旋流杯的强脉动结构为PVC和小涡旋结构。56°张角旋流杯方案的小涡旋结构能量更强,脉动幅度较大,能更好地对燃油进行剪切破碎,有利于强化雾化效果。     

富氢燃气旋转爆轰波传播特性实验研究

为研究富氢燃气旋转爆轰波传播特性,利用氢气与氧气预燃烧产生的富氢燃气作为燃料,空气为氧化剂,开展了旋转爆轰实验研究。对富氢燃气旋转爆轰压力变化、时频特性及传播速度等参数进行了分析,研究了不同传播模态下富氢燃气旋转爆轰波传播特性。研究表明：本文实验条件下,富氢燃气与空气旋转爆轰的传播模态主要受当量比影响,当量比高于1.06时呈现单波模态,随着当量比减小,旋转爆轰波呈现单波-双波过渡模态,即同一工况下,单波模态和双波模态交替出现,当量比减小到0.68左右时,基本呈现复杂的双波模态;在270 g/s的空气流量下,当量比增大,旋转爆轰波在环形燃烧室内的传播速度随之提高,但当量比到达临界点以后,传播速度提高不明显;在相同当量比下,当空气流量增大到370 g/s 时,旋转爆轰波的传播速度会进一步提高;空气流量越大,临界点对应的当量比越低,其中270 g/s空气流量对应临界当量比为1.32,370 g/s空气流量对应临界当量比1.16;达到临界当量比以后,传播速度受当量比和空气流量影响不大。     

亚声飞行Ma数变化对脉冲爆震发动机性能影响

采用经验数据和数值计算相结合的方法研究亚声飞行Ma数变化对带气动阀的脉冲爆震发动机(Pulsedetonation engine,简称PDE)的性能影响。分别计算了飞行Ma数在0.52,0.71,0.82和0.93时对PDE工作频率、推力和油耗的影响。计算表明当飞行Ma数从0.82增加到0.93时,频率上升11.1%,推力上升25.0%,油耗下降6.2%;而Ma数从0.82减小到0.52时,频率下降31.0%,推力下降43.6%,油耗增加14.8%。因此飞行Ma数变化对预估飞行器在亚声阶段的加速性能有较大的参考价值。     

增设收缩截面后的爆轰驱动激波管性能研究

为提高前向爆轰驱动激波风洞的性能，采用在驱动段靠近主膜处增设一个收缩喉道的方法，以此来产生较为均匀的驱动气源。实验结果表明此方法有效地削弱了爆轰波后Taylor稀疏波的影响。     

电离层中SEE现象的静电粒子模拟

利用网格粒子云（PIC）模型。在一定条件下，对电离层加热实验中静电波的激发和参量不稳定性的产生，作一维静电粒子模拟，模拟结果表明，SEE的特征与加热泵波频率ω0和电子回旋频率nωce之间的比值有关。入射泵波的电场与地磁场之间的夹角也是影响SEE的一个因素。