全文HTML
--> --> -->当将电源频率提高到甚高频乃至特高频频段时, 引入脉冲调制机制, 这时放电会出现不同于通常射频(13.56 MHz)放电的新现象, 比如在电压开启的第一个周期内会出现一次大电流放电现象[17,22], 即电压开启期间的第一个周期的放电电流峰值远远大于其稳定后的电流脉冲幅值; 并且基于粒子模拟的结果可知, 该过程同时会出现具有50 eV左右电子能量的高能电子[17], 这就意味着该类型放电具有潜在的巨大应用价值, 比如这些高能电子有可能打开某些化学键, 这样为大气压等离子体在医学领域的应用带来极大的便利[23-26].
但是, 当在电极上覆盖介质后, 脉冲调制是否对介质阻挡放电有重要的影响还没有深入研究[27]. 一般来说, 当电源频率增加至甚高频乃至特高频频段的时候, 基于传统的放电理论, 很多放电现象与极板表面是否覆盖介质已经关系不大, 因为快速变化的电场导致放电空间的带电粒子很难到达极板表面, 绝大多数带电粒子被束缚在放电空间内[3,4,28,29]. 但是已有的实验与计算研究表明, 介质的引入还是能够有效调控射频放电过程, 比如调控放电模式的转换等[30-32]. 因此当引入脉冲调制之后, 在甚高频乃至特高频频段下, 介质阻挡放电是否呈现出新的特点还需要进一步深入的研究.
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M1.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M2.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M3.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M4.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M5.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M6.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M7.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M8.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M9.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M10.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M11.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M12.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M13.png)
放电空间的总电流密度
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M14.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M15.png)
为了得到电子温度, 电子能量方程也需要自洽地求解,
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M16.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M17.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M18.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M19.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M20.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M21.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M22.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M23.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M24.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M25.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M26.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M27.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M30.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M29.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M28.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M30.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M29.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M28.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M30.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M29.png)
![](https://wulixb.iphy.ac.cn/cn/article/doi/10.7498/11-20191853_M28.png)
本文所使用的介质阻挡放电模型中, 放电间隙为0.05 cm, 介质板的厚度为0.03 cm. 射频频率固定为500 MHz, 该频率已经处于特高频频段. 实际上研究表明[22], 当放电频率提高至几百兆赫兹后, 脉冲调制放电中出现的大电流脉冲现象会变得极为显著, 因此在本文的研究中放电频率取为500 MHz. 模拟中所使用的气体为纯氦气, 其中所考虑的化学反应及反应系数均来自文献[34], 二次电子发射作为重要的边界条件在模拟中应加以考虑, 为简化计算, 这里二次电子发射系数取为0.03, 同时电子、离子及其他亚稳态粒子的初始密度取为1.0 × 108 /cm3, 在一定范围内, 初始粒子密度的取值对最终计算结果并没有影响. 本文所采用的计算程序的有效性已经通过与实验数据的对比得到了验证[3,4].
由于介质的引入, 在电压开启期间的第一个周期内, 如图1所示, 可以看到在正负半周期内均出现了一次大电流脉冲现象, 也就是出现“双电流脉冲”现象, 这与裸电极条件下只在电压正半周期出现一个大电流脉冲的放电现象是不同的[22-24]. 同时, 由图1可知, 正半周期的大电流脉冲幅值达到了23.4 A/cm2, 而负半周期的大电流脉冲幅值为23.1 A/cm2, 即正半周期大电流幅值比负半周期稍微大一点. 而随着外加电压的持续施加, 放电稳定后的电流脉冲正负半周期幅值相等, 均为19.2 A/cm2. 而且, 由于在裸电极条件下大电流脉冲中会同时产生大量的高能电子[22,23], 考虑类似的机理, 这些高能电子应该也会在介质阻挡放电中出现, 而且会在正负半周期各产生一次, 这就使得该类型的放电具有较大的应用优势与潜力. 然而, 考虑到本文所使用的流体模拟是基于局域场近似的, 此处并不深入分析双电流脉冲的产生机理及高能电子的产生过程(高能电子的问题将另文基于粒子模拟进行讨论), 而是更多地关注双电流脉冲随放电参数的演化过程, 这也更能充分体现流体模型计算速度快、遍历参数范围大且计算曲线光滑等优势.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-1_mini.jpg)
Figure1. Temple evolution of current densities in DBDs with a modulation frequency of 6.25 MHz, voltage amplitude of 800 V and duty cycle of 60%.
2
3.1.占空比的影响
在给定调制频率条件下, 随着占空比的逐渐增加, 电压开启时间越来越长, 相应的电压关断时间越来越短, 乃至当占空比为100%的时候放电将逐步转化为通常的连续性放电. 为保证数值模拟结果稳定性的同时, 能够更明显地观察双电流脉冲现象, 本节将射频电压幅值设定为800 V, 调制频率设定为6.25 MHz, 电压调制比为0. 图2分别给出了占空比为10%, 30%, 50%, 70%, 90%及100%的情况下的脉冲调制电流密度波形. 此时一个脉冲调制周期持续160 ns, 含射频周期数约80个. 在调制频率不变的情况下增加占空比, 从图2中可以直观地看出电压开启期间的射频周期数明显增加. 占空比为30%时, 电压开启期间的射频周期数约为24个, 电压关断时间为112 ns, 双电流脉冲的正峰值达到了21.6 A/cm2, 而稳定后的电流脉冲峰值为18.8 A/cm2, 两者差值为2.8 A/cm2. 占空比为70%时, 电压开启期间的射频周期数为56个, 电压关断时间为48 ns, 双电流脉冲的正峰值达到了23.6 A/cm2, 而稳定后的电流脉冲峰值为19.3 A/cm2, 两者差值为4.3 A/cm2. 从图2整体上来看, 双电流脉冲峰值与稳定后的电流幅值相比差值是比较大的, 而且随着占空比的逐渐增大, 双电流脉冲峰值与稳定后电流脉冲峰值的差值的变化规律也是比较复杂的, 基于计算数据, 将其变化规律总结在图3中.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-2_mini.jpg)
Figure2. Temporal evolution of current densities at a given modulation frequency of 6.25 MHz and voltage amplitude of 800 V for various duty cycles from 10% to 100%.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-3_mini.jpg)
Figure3. Peak values of current densities as a function of duty cycle at a given modulation frequency and voltage amplitude.
图3详细给出了双电流脉冲正负峰值及稳定后的电流脉冲峰值随占空比的变化曲线. 由图3可知, 随着占空比的增加, 双电流脉冲的正峰值总体呈现出先增大后减小的趋势, 并且大约在占空比为70%时达到最大值, 该值接近23.6 A/cm2. 而双电流脉冲的负峰值与正峰值的变化趋势相同, 只是负半周期的电流密度幅值略微小一些. 稳定之后的电流脉冲峰值随着占空比的增加而逐渐增大, 最后在占空比达到100%(即放电连续)时电流脉冲峰值接近20.2 A/cm2, 与双电流脉冲的正负峰值重合.
为了进一步分析说明上述变化趋势出现的内在机理, 图4给出了第一个电压周期的正半周期电流峰值时刻的电场强度空间分布情况, 图5给出了相同时刻的电子与离子密度空间分布, 其中两图均为左侧为瞬时阳极, 右侧为瞬时阴极. 其中的放电条件中固定电压为800 V, 调制频率为6.25 MHz, 调节占空比从10%—100%. 由图4可以看出, 当占空比为100%, 即为连续放电的时候, 阳极与阴极电场方向为异号, 其中阳极和阴极附近的电场强度分别突变为–38 kV/cm和75 kV/cm, 电场的空间分布为正常的射频放电情况下的电场分布, 这已经有了较为充分的讨论[2-4]. 当占空比小于100%时, 电场分布规律基本一致, 均是在阴极和阳极附近出现相同极性的电场, 这与占空比为100%的情况是完全不同的. 阴极鞘层电场的最大值在40 kV/cm左右, 而阳极的电场数值较小, 最大在17.5 kV/cm左右, 但是已经能比较明显地在阳极附近再次加速电子进而形成较大的放电电流, 并能提高电子能量形成一定数量的高能电子.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-4_mini.jpg)
Figure4. Spatial distribution of electric fields at the moment when the positive current density reaches the top value for various duty cycles.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-5_mini.jpg)
Figure5. Spatial profiles of electron density (solid line) and ion density (dash line) at the instant when the positive current density reaches the peak value for various duty cycles.
一般说来, 当电压关闭后, 等离子体在放电空间主要以扩散运动为主, 特别是在极板表面附近, 由于电子运动的较快, 离子运动的较慢, 最终形成相对稳定的双极扩散运动, 极板前面存在正离子区域以形成双极电场. 当电压开启后, 就阳极附近而言, 一般情况下, 随着电压的升高, 大量电子在电场的作用下到达极板并离开放电空间, 阳极附近形成较高密度的正离子区域, 从而在阳极附近形成较强的与外加电场反向的电场, 在电压连续施加的情况下也是这样, 如图5中的占空比为100%的情况. 然而, 如果外加电压变化较快, 比如电源频率为500 MHz, 由于外加电压在极短时间内上升到较大的数值, 导致在发生击穿的时候, 阳极极板附近的电子还没有完全响应电场的变化, 依旧有大量电子停留在极板附近, 如图5中占空比为0.1和0.4的情况, 从而在阳极与电子区域之间形成较强的正电场, 如图4所示.
图5描述了正脉冲峰值时刻电子密度与离子密度的空间分布. 由图5可知, 电子与离子的密度空间变化规律基本一致, 仅在极板鞘层附近的分布有明显的差别. 在占空比小于100% (即放电仍为不连续)时, 随着占空比的逐渐增大, 介质板附近的等离子体鞘层区宽度基本不变, 阳极鞘层宽度约为110 μm, 阴极鞘层宽度约为125 μm (结合图4). 而电子密度与离子密度均逐步增大, 在介质板之间的放电空间内, 电子密度和离子密度均从3.8 × 1013 /cm3 (占空比为10%)升高至7.1 × 1013 /cm3(占空比为40%). 而在占空比达到100%时, 电子密度与离子密度在介质板之间的空间内继续增大, 均达到了1.3 × 1014 /cm3. 但是在介质板附近的鞘层区却呈现出了减小的趋势, 这与图4中鞘层的变化趋势大体一致. 此时介质板附近离子密度要大于电子密度, 两者的密度梯度差随占空比的变化决定了图4中的电场分布变化.
当然, 结合电场(图4)与密度(图5)的分布, 也可以更好地理解图3中电流密度的变化, 简单地说电流密度正比于电场与等离子体密度的乘积, 图4可以看出阳极鞘层的电场是逐渐减弱的, 而与此同时图5表明等离子体密度是逐渐增加的. 随着占空比的增加, 电子密度的增加对电流的贡献较大, 大电流脉冲逐渐上升, 而另一方面当阳极电场减弱到一定程度的时候, 虽然粒子密度依旧增大, 但是电流值已不再增加, 逐渐下降.
到了外加电压的负半周期, 两极板电压极性反转, 特别是沉积在介质表面的壁电荷产生的电场比较显著, 可以影响放电空间中, 特别是介质板附近的电场, 这也是介质阻挡放电的重要特性[5,35]. 如图6所示, 当占空比为100%的时候, 电场的分布类似于正常的射频放电的电场分布, 并且其分布情况与正半周期除了阴极与阳极反转之外也完全相同. 但是当占空比小于100%的时候, 电场的分布呈现出新的特点, 主要体现在阳极附近(图6右侧), 从左至右先出现与阴极同向的电场, 该电场主要是由于离子区域的存在与表面电荷共同作用的结果, 这个电场将可以有效地再次加速电子; 然后是与阴极反向的电场, 其中阳极附近的电子区域起了重要作用.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-6_mini.jpg)
Figure6. Spatial distribution of the electric fields at the moment when the negative current density reaches the top value for various duty cycles.
同时可以看到, 此时在阴极附近的电场值得到了极大的增强, 要明显大于上半周期的数值, 这显然是表面电荷极大地增强了极板表面的电场, 这也体现了介质阻挡放电的特点.
可以概括地说, 电压开启前的剩余等离子体密度空间分布以及由此引起的电场的空间分布, 是影响电压开启后放电行为的重要因素. 如果电源频率非常高, 即电压变化非常快, 放电空间的剩余等离子体来不及响应外加电场的变化(如图5所示), 放电时阴极与阳极均会有较高的正电场(如图4), 从而在电压开启后第一个周期形成较大的电流密度脉冲.
2
3.2.调制频率的影响
调制频率对放电也有重要的影响, 在给定占空比的情况下, 提高调制频率意味着电压关断时间与开启时间均缩短. 图7给出了占空比为60%且外施电压峰值为600 V时, 改变脉冲调制频率(从6.25 MHz开始逐步增加至50 MHz为止)所得到的电流密度脉冲波形图. 当调制频率较低的时候, 比如6.25 MHz, 如图7(a)所示, 电压开启期间与关断期间的时间都相对较长, 在电压开启期间施加的电压周期数较多, 达到了48个电压周期, 产生的等离子体密度也较高, 同时较长的电压关断时间也能使等离子体的空间输运过程比较充分. 当调制频率较高的时候, 电压开启时间与关断时间均较短, 如图7(f)所示, 当调制频率为50 MHz的时候, 在电压开启期间仅有6个电压周期, 电压关断时间也非常短, 仅为8 ns. 由图7可知, 随着脉冲调制频率的增大, 双电流脉冲的正负峰值均渐渐减小, 而稳定后的电流脉冲峰值几乎不变. 当脉冲频率达到较高的50 MHz的时候, 双电流脉冲现象几乎消失.![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-7_mini.jpg)
Figure7. Temporal evolution of current densities at a given voltage amplitude of 600 V and duty cycle of 60% for various modulation frequencies from 6.25 MHz to 50 MHz
图8是通过采集图7各个图像中双电流脉冲的正负峰值及稳定后的电流脉冲峰值得出的电流脉冲峰值变化曲线. 随着调制频率的增加, 双电流脉冲的正峰值从16.1 A/cm2减小至13.6 A/cm2, 负峰值从15.4 A/cm2减小至13.3 A/cm2, 两者均呈现单调减小的趋势, 而稳定后的电流脉冲峰值始终保持在13.5 A/cm2左右. 这意味着当调制频率大于25 MHz后, 负半周期的电流脉冲峰值已经小于稳定后的电流峰值了.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-8_mini.jpg)
Figure8. Peak values of current densities as a function of modulation frequency at a given applied voltage of 800 V and duty cycle of 60%.
从计算结果可以得出, 调制频率的增加不会影响稳定后的电流脉冲峰值, 但是对于双电流脉冲的正负峰值有较大影响. 主要原因在于随着脉冲调制频率的增加, 一个脉冲调制周期的时间越来越短, 由于放电稳定后的电流基本不变, 而电压关断时间的变短, 却会影响电压关断期间带电粒子的运动(主要以扩散运动为主), 这就不可避免地影响了下一次电压开启前空间电场的分布, 进而影响电压开启后第一个电压周期的放电电流.
2
3.3.电压调制比的影响
通过逐渐增加电压调制比, 即逐渐增加电压关断期间内的电压幅值, 从另一角度实现了非连续性放电向连续性放电的转化过程, 有利于更好地从理论上分析电压关断期间的等离子体输运过程, 特别是其最终的空间分布对后继放电行为的影响.图9表示电压固定为800 V, 调制频率固定为6.25 MHz, 电压调制比从0 (对应电压为0 V)上升到1.0 (对应电压为800 V)情况下脉冲调制电流密度波形. 从图9中可以直观地看到, 随着电压调制比的逐渐增大, 电压开启与关断期间的电压差值在逐渐缩小. 在电压开启期间稳定后的电流脉冲峰值大小基本维持在19.5 A/cm2左右, 这也是电压调制比为1.0, 即完全连续放电的情况下(图9(f))的放电电流值. 而双电流脉冲的大小却受电压调制比的影响较大, 这点从图10中也可以看出.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-9_mini.jpg)
Figure9. Temporal evolutions of current densities at a given modulation frequency of 6.25 MHz and voltage amplitude of 800 V for various voltage modulated rates from 0 to 1.0.
![](https://wulixb.iphy.ac.cn/fileWLXB/journal/article/wlxb/2020/11/PIC/11-20191853-10_mini.jpg)
Figure10. Peak values of current densities as a function of voltage modulated rates at a given modulate frequency of 6.25 MHz and voltage amplitude of 800 V.
图10通过采集图9各个图像中双电流脉冲的正负峰值及稳定后的电流脉冲峰值得到了电流脉冲峰值的变化曲线. 从图10中可以看出, 双电流脉冲正负峰值的最大值出现在电压调制比为0时, 且正负电流脉冲峰值均达到23.1 A/cm2以上. 随着电压调制比的增加, 双电流脉冲幅值先逐渐减小, 甚至在电压调制比为0.35左右后小于稳定后的电流脉冲峰值, 在电压调制比为0.6左右达到最小值18.5 A/cm2, 然后逐渐增加, 直到电压调制比为1.0时, 各电流值趋于一致, 达到19.8 A/cm2.
从计算结果可以看出, 电压调制比的增加对于稳定后的电流脉冲峰值几乎没有影响, 对电压开启后第一个周期产生的双电流脉冲有较大影响. 其主要原因在于, 随着电压调制比的增加, 电压关断期间的外施电压幅值逐渐升高, 此时施加的外加电场可以有效调控电压关断期间等离子体的空间分布, 导致在电压开启前无法在阳极附近形成较强的正电场(如图4所示), 这样在电压开启后也就无法形成显著的大电流脉冲现象(大电流脉冲逐渐下降). 甚至随着电压调制比的逐渐增大, 在电压开启前阳极附近逐渐形成类似于连续放电的的负电场, 大电流脉冲现象不再出现, 且由于此时外加电压较小, 导致电压开启后放电形成的电流也比较小. 而随着电压调制比的进一步增大, 电压开启后形成的电流也逐渐增大, 直至当电压调制比为1的时候, 电压关断与电压开启期间形成统一的放电电流密度值.