COMSOL等离子体模块电子碰撞反应速率系数与Bolsig+计算结果差异及修正

最新推荐文章于 2025-04-25 22:17:35 发布

落叶_小唱

最新推荐文章于 2025-04-25 22:17:35 发布

阅读量8.1k

点赞数 13

分类专栏： COMSOL 文章标签： comsol 等离子体 bolsig+ 碰撞截面数据 afivo-streamer

本文链接：https://blog.csdn.net/ouening/article/details/125881321

版权

COMSOL 专栏收录该内容

1 篇文章

订阅专栏

本文探讨了如何通过comsol的等离子体仿真模型，结合air_chemistry_v3.txt中的运输数据和BOLSIG+碰撞截面数据库，进行电子碰撞反应及化学反应速率系数的计算。重点介绍了数据转换和工具应用，以及在不同软件间的单位一致性处理。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

数据准备：
（1）comsol等离子体仿真模型（ 包含N2和O2电子碰撞反应以及化学反应）
（2）air_chemistry_v3.txt : 该数据是afivo-streamer等离子体仿真项目中提供的输运数据（transport data），如下图所示：

在这里插入图片描述
对查表法而言，例如C25 N2 Ionization 15.6 eV，在文件air_chemistry_v3.txt中展示如下：

理论上也可以去网站https://nl.lxcat.net/data/download.php下载Phelps碰撞截面数据库，然后用bolsig+计算求解，保存计算结果（注意要用Combined table格式）

然后再用afivo-streamer/tools提供的python转换脚本转换：
python .\bolsig_convert.py .\bolsig_combined_output.dat out_bolsig.txt
其结果是一致的。

COMSOL中的速率系数展示的是电子平均能量和速率系数的关系，如下图所示：
在这里插入图片描述
而在air_chemistry_v3.txt文件中是约化场强(E/N）与速率数据之间的关系，需要进行简单转换，注意到air_chemistry_v3.txt中有约化场强和电子平均能量的关系，根据对应关系转换即可。另外还需注意的是由bolsog+计算的速率系数表达式为：
在这里插入图片描述
单位量纲和comsol的不一致，如果在comsol中要用到bolsig+计算的速率系数还需乘上阿伏伽德罗常数NA=6.022e+23/mol. matlab代码可视化如下所示：

最后需要注意的是comsol等离子体模型的电子密度和能量中电子传递属性设置查找表时，其迁移率和扩散系数的单位和bolsig+计算结果是一致的：
在这里插入图片描述