读论文《Modeling of fast ionization waves in pure nitrogen at moderate pressure》PSST 2021

摘要

采用耦合二维等离子体流体模型，研究了在27 mbar毛细管中纯氮中由正高压纳秒脉冲引发的快速电离波放电。对基于局部平均能量近似的二维流体代码进行了验证和应用，提出了纯氮的扩展三指数亥姆霍兹光电离模型。分析了快速电离波FIW的发展和结构，计算了电场和电流，并与实验测量结果进行了比较；分析了FIW发展阶段和余辉阶段电子N2(C3)的径向分布演化，电子密度在FIW发展阶段的径向分布呈“空心”结构，N2(C3)的时空演变主要受N2(A3)池化反应和电子猝灭；讨论了光电离对FIW径向分布的影响，提出了氮气放电时光电离的等效背景电子密度为10¹³m⁻³。得到了比能沉积(SED)、快速燃气加热(FGH)能量和温升的空间分布，加热效率随电场E/N和SED的变化而变化，高SED时加热效率趋于10%。在近轴区域，N2(A3)和N2(B3)的池化反应贡献了高达80%的FGH能量，氮分子的电子碰撞解离贡献了约10%，e-N₂⁺复合和N2分子对N(2D)原子的猝灭贡献了其余部分。

引言

介绍了如下内容：
1、纳秒脉冲放电的应用（典型的非平衡等离子体）：流动控制、点火助燃、能源转换、燃料重整、生物医疗等。
2、FIW优点：放电可重复性高、等离子体分布均匀、高约化电场、高沉积能量、高传播速度、高电离度、强化学活性。
3、其他文献已做的快速电离波的实验诊断和数值仿真研究
4、本文目的：通过耦合代码模拟纯氮的FIW，以获得详细动力学信息。模拟条件与实验相同，通过与实验结果对比验证模型。

模型与实验

模型：PASSKEy-流体模型，漂移-扩散-反应方程，平均电子能方程，电场的泊松方程，光电离的亥姆霍兹方程和流体动力学的欧拉方程耦合。

模拟几何图如下

结果与讨论

运用2D流体代码PASSKEy研究了中等气压下的快速电离波(FIW)，计算了p=27 mbar、初始气体温度Tgas=300 K、直径1.5 mm的纯氮毛细管放电，验证了基于LMEA的流体模型和扩展三指数Helmholtz光电离模型。

首先讨论了氮气FIW的传播，在近10 ns传播阶段，电离波头部的电子密度可达10²⁰ m^-3，平均传播速度为0.7 cm/ns，头部电场超过10000 Td，产生的电离速率在10³⁰ m³/s以上，电子密度的分布呈现出空心结构，径向电子密度峰值偏离对称轴。

将计算所得电场和电流与电参数测量对比，计算和实验在数值和趋势上定量一致，随后分析了起辉阶段和通道封闭前造成计算和实验差异的原因。

分析了电子和N2C3的时空演化，在FIW传播阶段，N2C3密度的径向分布重复电子密度的径向分布趋势，峰值密度偏离对称轴。在后辉光阶段，电子径向分布趋势和数值上改变缓慢，且由于N2A3和N2a1结合电离而重复上述两种组分的分布趋势；N2C3径向分布趋势在后辉光阶段变动较大，其时间演化可归结为N2A3汇集反应和电子熄灭N2C3反应的竞争机制。N2A3分裂对正确计算N2C3的时空分布非常关键。

研究了光电离和分步解离对氮气FIW的影响，通过比较电子碰撞电离和光电离速率，给出了氮气放电的等效背景电子密度范围4∼6×10¹³ m^-3，较大气压空气流注放电低两个数量级。通过开关Helmholtz模型研究了光电离对电子和N2C3径向分布的影响，背景电子密度越小，电离通道越细。分步解离——分子激发态N2 (A3, B3, C3)的电子碰撞解离在高SED下对N2C3、N4D、N2D影响较大，在降低N2C3密度至原来一半的同时可将N原子密度产率提升至原来的3倍。

分析了能量特性和快速气体加热(FGH)，放电总沉积能量为12.7 mJ，其中约9.4%用于FGH。加热效率是电场E/N和SED的函数，随SED增加其值趋向于10%。归纳了对FGH贡献较大的反应，在管中心区域r<0.45 mm，N2A3和N2B3的汇集反应贡献率可达80%，电子碰撞解离贡献10%，e-N2^+复合和N2熄灭N2D贡献了剩余的能量。温度分布趋势与电子密度相同，250 ns内FGH导致的局部最大温升可达1200 K。

论文获取

Chen Xiancong, Zhu Yifei, Wu Yun, et al. Modeling of fast ionization waves in pure nitrogen at moderate pressure. Plasma Sources Science and Technology, 2021, 30(6): 065002