惰性颗粒抑爆过程的数值模拟

对高温火球诱导的爆炸及惰性颗粒的抑爆过程进行了数值研究。由于扰动特征时间 ,两相弛豫特征时间和化学平衡特征时间的差异 ,这种带激波的两相化学反应流守恒方程具有刚性。利用处理包括组元守恒方程的全耦合TVD格式和可处理方程源项的Lax Wendroff Rubin格式分别求解气相和颗粒相 ,并根据各特征时间的比值用分步法处理方程的刚性。计算结果反映了均温燃烧火球在可燃环境中形成加速火焰和激波 ,进而成长为爆轰波的过程 ,以及惰性颗粒对爆轰波的抑制过程 ,反映了激波、化学反应和惰性颗粒之间的相互作用和耦合。计算结果表明 ,仅当颗粒浓度大于某值时才能有效抑制爆轰 ,否则无论多大的颗粒散布区域 ,皆不能使爆轰波完全抑制 ,爆轰波或是在其中发展为两相爆轰波 ,或是穿越该区域后重新成长为气相爆轰。