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

数据准备：
（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+计算结果是一致的：