1. 为什么要用 COMSOL 做棒板电极流注放电仿真棒板电极流注放电是气体放电研究里非常经典也很有代表性的一类非均匀场问题。高压电气工程里人们关心它的击穿特性等离子体应用领域关心它能不能稳定地形成放电通道诊断研究里又困于实验很难拍到几个纳秒内流注头部的发展细节。我最早接触这类仿真是为了弄明白一个很小的问题同样的间隙距离为什么换成正棒负板之后放电起始电压明显不同。单纯靠实验去测得到的是一个总电流和电压的结果中间的暂态过程完全是黑箱。后来转向氩气等离子体仿真用 COMSOL 软件做解算和后处理才逐渐把流注头部的高场增强、空间电荷壁障、通道建立这些物理过程一个一个从仿真结果里拎出来看明白。这篇内容没有高深的数学推导重点是把一个最小可用的流注放电仿真模型讲清楚物理模型怎么选棒板电极几何怎么建氩气的反应机理怎么搭网格怎么剖求解器怎么调以及结果怎么看。内容主要面向两类读者一类是刚开始用多物理场仿真做气休放电的人另一类是已经有了静电或稳态等离子体仿真基础想进一步触碰瞬态流注发展的人。我会尽量按“实际跑通一个模型”的顺序来写并且把踩过的坑放在对应环节里提而不是憋到最后才讲。1.1 流注放电到底难在哪里流注放电本质上是电子雪崩在极不均匀场中失稳发展的结果。棒板结构里棒电极附近电场极高电子在这里获得足够能量电离产生的大量电子又会形成新的雪崩当空间电荷的畸变电场足够强时前方出现一个自持的发射区域称为流注头。流注头以每秒十的五次方到六次方米的速度向对面电极发展后方则会留下一条近似等离子的导电通道。难就难在这个过程是多个物理场耦合在一起的电场由带电粒子空间电荷决定而粒子的生成率和迁移率又反过来依赖于电场。用解析公式去估算起始电压还好但想还原流注从起始到贯穿的完整过程必须要做时空上的瞬态求解。普通稳态电晕仿真做不了这件事因为它默认空间电荷效应已达到平衡实际流注发展过程中电子密度和电场在几十纳秒内剧烈变化求解器稍有不当就会出现负密度、振荡甚至发散。1.2 为什么选择氩气作为放电气体很多流注实验在空气中进行但仿真初建模型时我更推荐氩气。原因很实际空气里有氧气和水汽电子附着效应很强负离子参与化学反应反应体系动辄几十个组分数值刚性问题非常头疼。氩气是惰性气体负离子基本可以忽略最小模型只需要考虑电子、氩离子和亚稳态氩原子就能呈现出流注放电的核心特征。氩气还有一个好处就是碰撞截面数据比较成熟。国际上有 LXCat 等开源数据库可以直接查到电子碰撞激发、电离和弹性散射截面用 BOLSIG 这类玻尔兹曼求解器离线计算出约化场 E/N 对应的电子迁移率和反应速率系数再以插值表的形式导进 COMSOL整个参数链很清晰。相较之下如果你一开始就上六氟化硫或者潮湿空气数据不全是一回事反应硬编码起来也很难调试。至于选 COMSOL 而不是自己写有限差分程序我的理由比较朴素COMSOL 的 PDE 接口可以自己定义控制方程同时又有非常成熟的求解器、网格生成和后处理模块能让我把精力放在物理机制而不是程序开发上。尤其是 2D 轴对称模型内存消耗相对可控参数扫描很方便结果可视化也直观。如果你要做的只是单次对称针板放电用 COMSOL 的等离子体模块或者自定义 PDE 方案完全可以跑出现象来。2. 棒板电极模型搭建与物理接口选择2.1 几何、极性与关键参数先看几何。棒板电极在实验中一般是一根有一定曲率半径的棒垂直于一块大平板放置。仿真中没必要建三维用二维轴对称就可以把对称轴设为 z 轴棒电极看作对称轴上的细长圆柱端部做成圆弧或小平头板电极则是一段带半径的水平线段。我常用的一组示意参数如下参数建议值备注棒半径0.5 mm曲率越小电场越集中流注起始越容易棒长度10 mm主要影响远端边界尽量留长一点板电极半径20 mm要明显大于放电区域避免外部边界干扰间隙距离5 mm以实际实验为准这里只作示例气压1 bar氩气温度 300 K施加电压5–10 kV建议从小到大做斜坡加压极性选择上正棒负板是比较典型的正流注场景。流注首先在阳极棒尖附近形成然后向阴极板方向贯穿负极性情形下放电往往从阴极发射起始机制略有不同仿真边界条件也要相应调整。如果你刚开始做建议先从正极性入手因为正流注现象更规整头部发展路径也容易和实验照片对照。2.2 等离子体模块还是自定义 PDE在 COMSOL 里建流注仿真有两条路一是用等离子体模块里的直流放电接口二是用系数形式 PDE 自己搭方程。我两个都试过总体感受是内置等离子体接口在低气压、大面积放电上比较成熟比如辉光放电或介质阻挡放电但对棒板电极流注这种强不均匀、强瞬态的问题内置接口的默认反应模型和边界设置反而限制了发挥我最后更倾向于自定义 PDE 加静电接口的组合。从物理上说流注放电的流体模型由三块组成电子连续性方程、离子连续性方程含可能的亚稳态粒子、静电泊松方程。以电子为例漂移扩散近似下方程是$$\frac{\partial n_e}{\partial t} abla \cdot \left(-n_e \mu_e \mathbf{E} - D_e abla n_e\right) S_e$$右侧 $S_e$ 是电离、激发、复合等反应源项。泊松方程右侧是空间电荷$$\rho e \left(n_i - n_e\right)$$每个粒子的迁移率 $\mu$ 和扩散系数 $D$ 可以定义为约化电场的函数也可以借助爱因斯坦关系由迁移率反推扩散系数。COMSOL 的“系数形式 PDE”接口允许你直接输入这些方程每个粒子的密度作为一个因变量再加上静电接口里的电势变量 V就能完成耦合。这样做的好处是你能清楚地看到每个源项怎么加、边界怎么设出了数值问题也更好定位。2.3 边界条件电极和开放边界的设置棒板电极仿真里边界条件比几何更容易踩坑。先看电极边界棒阳极设置固定电势板阴极设置为 0 V。除此之外电子和离子在电极表面的处理会直接影响二次电子发射和鞘层结构。常见的做法是给电子设置一个“电子通量”边界包含漂移扩散通量加二次电子发射项。二次电子发射系数 $\gamma$ 不是随便给的它应该取决于入射离子能量和电极材料的实际二次发射特性。对于氩气放电$\gamma$ 取 0.05–0.1 是常见的经验范围但如果你只在流注传播阶段做仿真、不关心阴极鞘层的微小差异也可以先把 $\gamma$ 设为零减少收敛难度。开放边界设在板电极外侧和棒后端要距离放电核心足够远。常见的错误是模型外边界离流注通道太近导致等电位线被外部人工边界压缩结果场强偏低。经验做法是外边界离轴线的距离至少取间隙距离的两倍以上必要时再用无限元域把人造反射效应压下去。3. 氩气等离子体化学反应体系与输运参数3.1 最小反应机理几个方程就够做流注仿真不等于机理越全越好。氩气放电流程中有很多微观过程但对棒板电极流注发展影响最大的也就是电子碰撞电离、亚稳态激发和亚稳态参与的附加电离。我建议第一版模型只保留四个粒子电子 e氩正离子 Ar氩亚稳态 Ar*背景中性氩原子 Ar相应的核心反应如下反应类型说明e Ar → e Ar弹性散射交换能量并影响电子迁移e Ar → e Ar*电子激发形成亚稳态阈值约 11.5 eVe Ar → 2e Ar直接电离阈值约 15.8 eV流注头部的主要电离来源e Ar* → 2e Ar分级电离亚稳态再电离有效降低电场门槛Ar* Ar* → e Ar Ar潘宁电离两个亚稳态碰撞产生额外电子注意实际速率常数不能直接在常数表里抄一个数因为放电过程中约化场 E/N 一直在变。常见做法是用 BOLSIG 求解玻尔兹曼方程得到不同 E/N 下的 EEDF 电子能量分布函数再与碰撞截面做积分得到各反应速率系数和电子迁移率。把这些结果做成插值表在 COMSOL 里作为变量函数引用。3.2 如何把速率系数导入 COMSOL具体操作上我会在外部脚本里把 BOLSIG 算出来的若干组 E/N 值比如从 10 Td 到 1000 Td共 40 到 80 个点导出成两列或三列的文本文件。然后在 COMSOL 的“全局定义”里用插值函数加载导入。比如定义函数kiz(EbyN)表示直接电离速率系数单位用m^3/s。此后在 PDE 源项里直接写kiz(En)*ne*nAr其中En是约化场变量normE*1e21/ nArnormE是电场模nAr是中性粒子数密度。用插值函数的好处是求解器在迭代时能自动求导收敛特性比用 if 判断写好得多。不要因为嫌麻烦就把速率系数设成常数否则你大概率会发现流注根本不发展或者起始电压偏离实验太多。3.3 迁移率、扩散系数与电子能量近似氩气里电子迁移率同样强烈依赖于 E/N。不能把它当成常数否则流注头部传播速度会失真。建议也和反应速率系数一样做成插值函数。离子迁移率变化相对不大但也别直接抄常温下的单值先确认你的气压范围是否合适。扩散系数方面电子扩散系数可以用 Einstein 关系估算但注意电子温度不是常数。更严谨的做法是再加一个平均电子能量方程用它来给电子迁移率和反应速率提供本地电子温度。也就是说模型变成五变量电子数密度、电子能量密度、Ar 数密度、Ar* 数密度、电势。我建议第一版先只做电子数密度加局域场近似也就是认为电子能量主要由当地 E/N 决定不额外算能量方程模型跑稳以后再升级成非局域模型。这个顺序可以避免一上手就被电子能量方程的小时间常数折磨得寸步难行。4. 网格、求解器和时间步进怎么把模型跑稳4.1 网格策略棒尖附近必须细密棒板电极流注仿真对网格的要求和前几章说的普通电晕仿真不太一样。流注头部尺度很小往往在几十微米量级如果网格太粗空间电荷产生的电场峰值被平均掉流注起不来更谈不上传播。实际操作中我会在棒尖端和轴线附近专门做一个细网格区棒电极端部圆弧处网格尺寸取 0.01 mm 甚至更小流注路径所在的轴线区域网格尺寸控制在 0.02–0.05 mm往外再逐渐放宽到 0.5 mm 甚至几个毫米。同时在电极表面加边界层网格因为电极附近鞘层区的密度梯度很大。网格质量直接影响时间步长。有时候放电物理上是稳定的但网格太粗时数值黏性会把流注头部磨平。定期检查一下剖分质量不要只看“最小质量”报告重点关心棒尖附近的单元长宽比。4.2 初始条件、电压斜坡和稳定性流注仿真需要给一个背景电子密度否则雪崩无法起始。这个值可以从宇宙射线本底电离或实验预电离来估计通常取 10^10 到 10^12 m^-3。如果你设成零绝大多数求解器会直接给零除以零的崩溃。也不要试图用一次特别大的随机扰动来代替背景密度网格不均会造成伪放电。电压不要一步直接加满。物理上实验系统有电容和上升沿突然跳变的电压会在模型里激发高频振荡。建议设为斜坡函数V Vin * min(1, t/t_rise)上升时间取 10–50 ns与实验中高压脉冲沿大致匹配。即便如此初始阶段还是要自适应时间步长把电场变化放缓否则泊松方程解的变化太快瞬态求解器很容易发散。4.3 求解器调参与参数扫描对这样的多变量瞬态模型我很少用“全耦合求解器”一揽子代。分离式求解器更可控先把泊松方程算一轮再依次算电子密度、离子密度和电子能量再回到泊松方程循环。每一轮内部可以限制迭代次数和阻尼因子。COMSOL 6.4 里切换分离式和全耦合求解器的地方很直观但我建议直接选择分离式。时间步进则交给自适应算法设置相对容差在 1e-3 左右。太严苛的容差会让求解器在小时间步上反复挣扎太宽松又会把空间电荷细节抹掉。参数扫描可以交给 COMSOL 的辅助扫描用气压、电压或间隙距离作为扫描参数。如果扫描量大就配合 LiveLink for Python 或 MATLAB 批量修改模型参数把每组仿真的关键输出导出成 CSV后处理再汇总。像“棒半径 0.2 mm vs 0.8 mm 对起始电压的影响”这种问题用参数扫描比手动改一次跑一次要高效得多。5. 典型仿真结果与流注发展分析5.1 如何判断模型是否产生流注跑完以后不要急着看放电通道的漂亮彩图先检查几个无量纲指标。最直接的是电场增强幅度棒电极表面上未加空间电荷时的“几何电场”和加入空间电荷后的“实际电场”往往有明显区别。流注头部前方会形成额外场增强电子密度峰值迅速抬升到背景值的几个数量级以上。如果等离子体只在棒尖附近一直呆着形成一圈稳态电晕那说明电压低于流注起始阈值或者电离源项不足。真正的流注起始电离会在空间电荷反馈下自行加速头部电子密度增加空间电荷场再增强再刺激更多电离。仿真中你会看到某个时间点后电子密度曲线突然不再向四周均匀扩散而是顺着轴线方向拉出一条狭窄的高密度路径这就是流注通道的雏形。5.2 电场分布与空间电荷屏蔽画一条沿对称轴的电场模曲线你会发现流注没起来时电场最高点在棒尖表面沿轴线单调下降。流注发展起来后头部前缘会有一个明显的电场尖峰而通道内部电场反而被压低这是因为通道里电子和正离子大致中和形成了准中性等离子体。这里有一个很经典的物理读图技巧看等势线是不是被“压成扁椭圆”。等离子体通道接近导体会把外部电场挤压到头部前方造成更高的有效场强。仿真后处理里把电势分布和电子密度叠加在一起能非常直观地看到通道屏蔽外部场的过程。这也是流注能够低速起振、高速贯穿的根本原因。5.3 电流波形与放电阶段划分电极电流也是重要结果。在 COMSOL 里可以定义一个探测边界把电极表面的总电流密度对面积积分再加进位移电流项。流注发展时电流波形一般会有几个阶段先是指数上升的雪崩阶段电流很小但增速极快随后是通道向阴极推进的阶段电流出现较大的上升脉冲一旦通道贯穿间隙电流会直接进入导通阶段同时通道电阻迅速下降。有了电流波形就可以定量讨论放电模式。同样的施加电压低气压下可能只出现慢速正流注高气压下则会转成快速流注甚至击穿。做参数扫描时我把这些电流峰值和起始电压汇总成一张表再和实验对照就能发现模型哪里需要调整。6. 常见问题排查与避坑经验6.1 电子密度出现负值怎么办这是新手最常遇到的问题。根因通常是漂移扩散项的对流占优越过了稳定域或者反应源项的时间尺度比网格快。流注头部电子梯度很高如果不给密度的数值解加下限保护迭代时很容易在局部出现负密度。我建议两条腿走路一是在源项里对电子密度引用做下限截断比如写成max(ne, 1e8)避免负密度参与反应计算二是适当缩小时间步长让电子密度变化不要跨越几个数量级。也可以考虑把电子密度的因变量设成对数形式但这样会带来源项里很多非线性项调试成本更高非必要不推荐。6.2 不收敛或长期振荡如果求解器报错先不要怪时间步。把电压先降一半或者把背景电子密度调高一点看看是不是模型本身物理上太激进。常见原因是反应速率系数表在 E/N 上限处外推太猛建议把 E/N 范围设置得比实际需要宽一点并让插值曲线平滑外推为常数。再一个常见原因是边界条件刚度的零冲突棒尖附近网格特别细但开放边界离得不够远导致电势外边界对解约束过强。此时电场在外部边界出现“反弹”流注头部会像撞到玻璃板一样停止生长。解决方法是增大计算域或者改用无限元域吸收外部场。6.3 网格敏感性与结果不可重复如果结果每次调整网格都变说明网格没有收敛。最简单的判断方法把棒尖网格尺寸减半看最大电场和流注到达时间是否变化超过百分之十如果变化明显就要继续加密。注意细网格会带来很大的计算量建议先在一个短的间隙距离上验证网格收敛再放大正式工况。6.4 COMSOL 仿真耗时太长放电时间是几十纳秒但因为过程高度非线性仿真耗时长也很正常。减少耗时的几个实用招数把轴对称模型选成 2D 而不是 2D 轴对称去偷工减料这是错误示范所有计算都应该开轴对称关闭不必要的后处理存储只在探测点上保存数据分离式求解器里限制每一牛顿步的迭代次数。还有一个更有效的办法是先用很少的氩气反应机理把 mesh 和时间方案调对再逐步加反应。6.5 实验对照时模型还缺什么如果你的仿真和你自己的实验总是差一口气别急着怀疑求解器。先检查实验里气体纯度氩气里掺一点空气会让流注形态发生明显改变而你的模型还是纯氩气。必要时加入氮气或氧气作为杂质反应。另一个原因可能是电极表面粗糙度和二次电子发射系数不一致这在小间隙放电里尤其明显。我自己的经验是流注仿真永远是一个“物理优先、数值其次”的建模过程先把实验中的电压波形、间隙距离和气压三个输入看准再去追求复杂反应机理。模型里最重要的验证指标是最小起始电压和传播速度只要这两个量能和实验对上说明整个模型框架没有原则性错误。7. 一个可以继续深入的方向从单次流注到重复脉冲与参数控制每次评完这些坑都有人问我后面还能做点什么。我觉得比较值得投入的方向是把单次流注仿真扩展到重复脉冲场景。因为很多实际的等离子体应用比如材料处理、点火辅助用的是几十到几百千赫兹的重复脉冲脉冲间隔里前一脉冲留下的带电粒子和亚稳态还会影响下一次放电。COMSOL 里做这个扩展不算难。先跑完第一次流注把最终解作为下一次脉冲的初始条件再在电极边界上重新施加电压斜坡观察残余等离子体如何降低后续脉冲的起始电压。实践下来你会发现重复脉冲的起始电压往往比第一个脉冲低得多这就是记忆效应。配合 LiveLink for MATLAB 或 Python把多周期脉冲扫描自动化还能画出“起始电压随脉冲重复率变化”的曲线。另一个值得探索的方向是针板结构里的共轭电场与流体耦合。如果气压很低或者你想研究放电产生的高温气体膨胀就需要把层流单元耦合进来。不过先把纯静电和粒子输运问题理解透彻再考虑流动才能保证每一步结果的可靠性。棒板电极加氩气的最小模型是我认为平衡了物理丰富度和仿真可行性的最佳起点。