1. 项目概述：空气流注放电模型的物理意义与应用场景
空气流注放电是气体放电物理中一种重要的非平衡等离子体现象，广泛存在于高压绝缘、闪电形成、工业电晕处理等场景。当外加电场超过空气的击穿阈值时，电子雪崩会发展成丝状导电通道，这个动态过程涉及复杂的电子碰撞电离、空间电荷效应和光电离机制。
在Comsol Multiphysics中模拟这一现象，本质上是在求解耦合的麦克斯韦方程组、流体力学方程和化学反应动力学方程组。等离子体模块提供了现成的接口，让我们能够避免从零开始推导这些偏微分方程。我选择Comsol进行这项研究，主要基于三个实际考量：一是其多物理场耦合能力可以自然处理放电过程中的电磁-流体-化学相互作用；二是内置的Bolsig+接口能自动计算电子输运系数；三是后处理可视化功能强大，便于分析复杂的空间分布数据。
2. 模型搭建的关键技术环节
2.1 几何建模与网格划分策略
创建一个典型的针-板电极几何结构，针尖曲率半径设为50μm，极间距5mm。这种结构会产生高度不均匀的电场分布，有利于流注起始。网格划分需要特别注意：
- 针尖附近采用边界层网格，最小单元尺寸0.1μm
- 放电通道区域使用自适应网格细化
- 外围区域可以用较粗网格节省计算资源
重要提示：必须开启"变形几何"功能，因为流注发展会导致等离子体通道形变
2.2 物理场接口配置
核心物理场包括：
- 等离子体接口：选择"流注放电"预设，包含电子密度、正负离子密度三个因变量
- 静电接口：计算空间电荷修正后的电场分布
- 流体流动接口：模拟放电产生的气体加热和流动效应
- 化学反应工程：定义关键的23个电子碰撞反应（包括电离、激发、附着等）
关键参数设置示例：
comsol复制// 电子传输系数
electron_mobility = bolsig(Te)
electron_diffusion = bolsig(Te)
// 初始种子电子
ne_init = 1e12*exp(-(x^2+y^2)/1e-6^2) [1/m^3]
2.3 材料属性与边界条件
空气混合物定义：
- N2: 78%
- O2: 21%
- H2O: 1% (湿度影响很重要)
边界条件特殊处理：
- 电极表面使用二次电子发射系数γ=0.01
- 开放边界设置空间电荷限制条件
- 温度边界设为300K恒温（除非专门研究热效应）
3. 求解器配置与计算技巧
3.1 多物理场耦合策略
采用顺序耦合求解策略：
- 先求解稳态静电场（作为初始条件）
- 开启瞬态研究，时间步长从1ns逐步增加到100ns
- 每步先更新电场，再求解等离子体输运方程
3.2 非线性求解器调参
遇到收敛困难时尝试：
comsol复制// 调整阻尼因子
solvers.stepdampingfactor = 0.7
// 启用雅可比矩阵更新
solvers.jacobianupdate = on
// 限制变量变化幅度
solvers.maxdeltator = 0.1
3.3 高性能计算设置
对于大规模模型：
- 使用分布式内存求解器
- 开启矩阵预处理(ILU)
- 保存频率设为每10步存一次
4. 典型结果分析与验证
4.1 流注发展动态可视化
通过切片图展示：
- 电子密度分布（log10标度）
- 电场增强因子(E/E0)
- 空间电荷密度分布
动画输出建议设置为100fps，可以清晰观察到流注先导的跳跃式发展特征。
4.2 定量验证方法
-
与Peek公式对比起始电压：
math复制V_start = 31δd(1 + 0.308/√(δr))δ为相对空气密度，d为间隙距离，r为针尖半径
-
流注速度测量：
comsol复制// 通过最大电子密度位置随时间变化计算 v_streamer = d(position(ne_max))/dt
5. 常见问题排查指南
| 问题现象 | 可能原因 | 解决方案 |
|---|---|---|
| 初始不放电 | 种子电子不足 | 添加人工背景电离(1e6 m^-3) |
| 计算发散 | 时间步长过大 | 采用自适应步长，初始1ps |
| 内存不足 | 网格太密 | 先试算2D轴对称模型 |
| 物理场不耦合 | 求解顺序错误 | 检查变量依赖关系图 |
6. 进阶应用方向
- 混合模型开发：将流体模型与蒙特卡洛碰撞(MCC)方法结合，在关键区域使用粒子模拟
- 工业应用扩展：
- 电晕处理表面改性（需添加表面反应）
- 臭氧发生器优化（增加O3生成反应）
- 多尺度耦合：将微观放电参数传递给宏观设备级仿真
我在实际建模中发现，湿度对流注分支形态的影响比理论预期更显著。当相对湿度超过60%时，需要在反应体系中额外添加10个与水分子相关的碰撞反应截面，这可以通过Bolsig+的交叉截面数据库导入实现。