氧等离子体中臭氧形成过程数值模拟.doc

氧等离子体中臭氧形成过程数值模拟欧阳建明，邵福球，林明 东国防科技大学理学院，长沙，410073通讯地址：国防科技大学物理系，410073E-mail 地址：oyjmoyjm@摘要：本文利用数值模拟程序对不同条件下氧等离子体和大气等离子体中臭氧的形成过程进行数值模拟结果发现：臭氧的形成与游离态氧原子密切相关；臭氧的最大数密度与初始电子数密度有关；臭氧最大数密度持续的时间与脉冲放电大气等离子体中放电时间选取相关，最佳的放电间隔是游离态氧原子达到最大的时间一. 引言臭氧具有极强氧化能力，是强氧化剂、 杀菌消毒剂、催化剂、脱色除臭剂由于臭氧本身的不稳定性，经过 一段时间后能自行分解为氧气，不会造成二次污染随着人们对环境质量要求的提高，臭氧将是人们理想的绿色环保药剂早在 1857年，西门 子公式就开始利用臭氧进行消毒 [1] 目前人工产生臭氧主要是采用沿面放电、介质阻挡放电等方法，利用氧等离子体或大气等离离子体产生臭氧，相关理论和实验研究已经取得了很大的进展 [1-4]等离子体中臭氧的形成 过程数值模拟是目前理论研究的一个重要方向，也是较为直接的一种方法本文以 较为完善的大气等离子体物理化学过程模拟程序 [5]为基础，充分考虑与臭氧形成 过程相关的组分和反应式，对程序进行改进。

利用改进程序数值模拟不同条件下氧等离子体和大气等离子体中臭氧的形成过程二. 物理模型和数值模拟方法模型考虑了 41 种组分，其中包含有中性粒子（包括激发态氧分子和原子），带正电粒子（正离子、水合正离子），带负电粒子（电 子、 负离子、水合负离子）以及与这 41 种组分相关的 249 个反应式，包括二体反应，如 ，三体eOe2反应，如 各反 应式由资料查得 [2,6-8]232OO模型为零维模型，即只考虑各组分随时间的关系，而不考虑扩散、漂移等空 间过程，这 是由于形成臭氧的等离子体一般都被局限在一个小的空间当中，空间过程对臭氧的产生过程的影响不大模拟的方程组如下 [5]：),.(21lii ytfdy 41,2imjimjijijii yknknf4141其中 为第 i 组分的数密度， 是与第 i 种组 分反应相关的函数 和iyif ijk分别是 i、j 组分之间的二体反应的反应速率和 i、j、m 组分之间的三体反应的ijmk反应速率； 和 为系数，当反应生成 i 组分是取 1，消耗 i 组分时取-1ijnijm数值模拟程序采用拟稳态近似法（QSSA） [9]，中心思想是将不同的刚性度的方程分类，采用不同的公式 进行近似计算， 该方法在大气物理的模拟中广泛采用[10]。

模 拟程序采用了质量守恒技术 为了加快计算 时间，程序的时间步长可以自动调整同 时为了保证计算精度，程序还采用精度控制方法三. 模拟条件模型假设等离子体由气体放电产生，带电粒子主要为电子 和正粒子 ，和e2O，为 保持电中性条件，正负电 荷相同不同条件下各带电粒子的初始数密度不2N同各方应式的反应系数也由资料查得 [2,6-8]，其中有些反应系数与大气温度、电子温度有关，大气温度设为 300K，电子温度假设为 2ev，约合 23200K初始的时间步长为 t0 为 s，程序采用可 变步长 方法，在计算过程中，随150着离子数密度的减小步长自动增大，结合精度控制方法，在保证精度的前提下提高计算效率四. 结果分析利用数值模拟程序，模拟了氧等离子体和大气等离子体（大气中同样包含大量氧分子）中的物理化学过程1 氧等离子体中臭氧形成氧等离子体中初始带电粒子假设为电子和 图 1 给出了在初始电子数密2O度 条件下，与臭氧形成过程相关的几种组分： ， ， ， 随时间3120cme43O演化图，其中横坐标为数密度，纵坐标为反应时间，均取对数坐标从图中可以看出，电 子在反应开始的短 时间内生成 离子。

游离态氧原子 的生成与电子4O有关游离态氧原子产生的主要反应： ，该反应与电子相关反 应ee2开始后，由于电子具有较高的数密度，游离态氧原子 迅速产生随着 电子的衰减， 生成减缓，再 0.05μs 时刻达到最大值 ，在 1μs 后氧原子迅速O310cm衰减形成臭氧 游离态氧原子 生成的同时，氧原子与氧分子及其它组分反应生O成臭氧 在 产生氧等离子体的同时会形成大量的激发态氧分子 ，激3 )(12gaO发态氧分子不稳定，最终生成氧分子图 2 给出了不同电子数密度条件下臭氧和氧原子数密度随时间演化图从图中可以看出，随着初始电子数密度的增加，氧原子的最大数密度和臭氧的最大数密度都增加，并且与电子初始数密度近似成正比2 大气等离子体中臭氧形成大气等离子体是假设在大气中放电产生等离子体，大气中同样有大量的氧分子初始带电粒子为电子、 和 图 3 给出了在初始 电子数密度2ON条件下，与臭氧形成过程相关的几种组分随时间演化图图形与图 13120cm具有很大差别，电子衰减过程中首先生成大量的负离子 ， 在 0.1μs 时刻左2O右占优，而后才是 占优 [5]游离态氧原子迅速产生，在 0.1μs 达到最大值4O，达到最大的时间比同条件下的氧等离子体长，最大值比氧等离子体3102cm大约一个数量级，这也导致臭氧最大值较氧等离子体中大约 1 个数量级。

氧原子形成的同时产生臭氧，氧原子在 10μs 后迅速衰减，形成臭氧臭氧分子数密度在 0.1ms 达到最大值 ，此后保持一段时间，在 0.1 秒后迅速衰减激3102cm发态氧分子的过程与氧等离子体的基本相同图 4 给出了不同初始电子数密度条件下臭氧和氧原子数密度随时间演化过程从图 中可以看出，随着初始电子数密度的增加，氧原子和臭氧数密度最大 值增加但与初始电子数密度不成正比，初始电子低时较 氧等离子体大，初始 电子数密度高时与氧等离子体相同氧等离子体中臭氧形成后较为稳定，但大气等离子体中形成的臭氧不稳定，臭氧最大数密度保持一段时间后，臭氧衰减初始 电子数密度越大，保持的时间 越长初始 电子数密度 条件下可以保持到3140cm10s图 2 氧等离子体 不同电 子初始数密度 和 随时间演化 图O3图 1 氧等离子体中主要组分随时间演化图3 脉冲放电大气等离子体中臭氧形成前面模拟的是单脉冲放电形成等离子体中臭氧的形成过程，臭氧的最大数密度有限为了加大生成臭氧的最大数密度，可以采用脉冲放电多次产生等离子体，形成高数密度臭氧图 5 和图 6 是初始电子数密度为 条件下，脉冲放电 10 次产生臭3120cm氧的过程。

我们关心的是臭氧的形成过程，因此假 设 每次放电就形成等离子体，并且电子数密度都是 为了进行对比分析，我们选取两个放电间隔时3120c间，图 5 间隔时间是臭氧最大数密度时间 0.1ms；图 6 间隔时间是游离态氧原子最大数密度时间 1μs对比两幅图不难发现，臭氧的最大数密度都是 ，但是达到最大的时间3120cm以及持续的时间不同放电间隔 0.1ms 时，在 1ms 达到最大值，但持续的时间较短，在 10ms 后急剧衰减放电间隔为 1μs 时，在 10μs 秒达到最大值，最大值也为 最大值持续较长时间，一直持续到约 0.1s 才明显衰减因此实际3120cm的臭氧发生器放电间隔应当选取在氧原子数密度达到极大值的时间图 3 大气等离子体中主要组分随时间演化图图 4 大气等离子体 不同电子初始数密度 和 随时间演化图O3五. 结论本文以大气等离子体物理化学数值模拟程序为基础，充分考虑氧等离子体和大气等离子体中的组分， 对模型和模拟程序进行改进模型中包含 41 种组分及相关的 249 个反应式利用改进的数值模拟程序模拟不同条件下氧等离子体和大气等离子体中臭氧的形成过程，给出了臭氧及与臭氧形成相关的组分数密度随时间的变化过程，得到了以下结论：1，臭氧的形成与游离态氧原子的数密度密切相关，主要由游离态氧原子反应形成。

2，氧等离子体中臭氧最大密度持续时间较长；臭氧的最大数密度与初始电子数密度有关，近似成正比3，大气等离子体中同样可以大量形成臭氧；在低初始电子数密度条件下，臭氧最大数密度较氧等离子体中的高，高电子数密度下基本一致；臭氧保持一段时间后迅速衰减，初始电子数密度越高，持续时间越长 4，臭氧最大数密度持续的时间与脉冲放电大气等离子体中放电时间选取相关，最佳的放电间隔是游离 态氧原子达到最大时间参考文献图 5 脉冲放电 大气等离子体中 O 和O3 随时间演化图（脉冲间隔 10-4s）图 6 脉冲放电 大气等离子体中 O 和O3 随时间演化图（脉冲间隔 10-6s）[1] Takashi Kimura, Yosuke Hattori ，Akinori Oda，Japanese Journal of Applied Physics, Vol. 43（2004）, No. 11A, P: 7689[2] Atsushi Yamatake, Koichi Yasuoka and Shozo Ishii, Japanese Journal of Applied Physics, Vol. 43（2004）, No. 9A, P: 6381[3] Ashraf Yehia, Akira Mizuno and Mohamed A. M. EL-Osealy, Japanese Journal of Applied Physics, Vol. 43(2004), No. 8A, P: 5558[4] P.A. Balakrishnana, A. Arunagiria, P.G. Rao, Journal of Electrostatics, Vol. 56 (2002), P:77[5] Ou-yang J M, Guo W, Wang L and Shao F Q, 2004, Chinese Physics, Vol.13, P: 2174[6] Vidmar R J, 1990, IEEE Transactions on Plasma Science, Vol.18, 4733[7] Brasseur G and Solomon S, 1984, Aeronomy of the Middle Atmosphere (Boston: D. Reidel Publishing Company)[8] Heicklen J, 1976, Atmospheric Chemistry (New York: Academic Press)[9] Hesstwedt E．et al, Internal.J.Chem.Kinet ,Vol.10(1978), P:971[10] Wang Tiiian, li Zongkai, Quarterly Journal of Applied Meteorology, Vol.7(1996)， No.4, P: 466[王体健，李宗恺，应用气象学报，Vol.7(1996)，No.4, P: 466]。