棒板电极是气体放电里最经典的非均匀间隙结构。一根针尖对着一块平板施加几十kV的脉冲电压几纳秒内就能看到流注从尖端长出来然后沿着轴线一路发展、分叉、甚至击穿。这个现象背后既有“电子雪崩如何启动”这种基础物理又有“设备绝缘缺陷怎么评估”这种工程问题。这次想说的就是在 COMSOL Multiphysics 里用等离子体模块把棒板空气流注放电完整搭起来化学反应部分采用一套 22 反应的综合模型既算电子密度分布也算电势、电场畸变和活性产物生成。这篇文章适合三类人一是刚接触气体放电仿真、想在 COMSOL 里跑通第一个流注模型的同学二是做高压绝缘或放电实验、想用数值结果对照发光照片和电流波形的研究者三是想了解 22 反应这个规模到底够不够用、要不要往上加反应集的从业者。我会把模型怎么设计、反应怎么选、边界怎么给、求解器怎么调、常见坑怎么避都讲一遍内容按“能从零复现”的标准来写参数可以直接抄。1. 为什么做棒板流注放电模型从现象到数值1.1 棒板结构到底解决了什么棒板结构在气体放电研究里地位很高是因为它极好地制造了“非均匀电场”。尖端附近的电场强度远高于间隙平均值电场集中系数可以达到几十倍。这种局部增强让放电不需要很高的外加电压就能先发生电晕再发展为流注。工程上输电线下面的导体毛刺、GIS 内部尖端缺陷、绝缘子表面附着颗粒本质上都可以抽象成某种“棒-板”或“针-板”结构。所以拿棒板间隙做流注研究不只是为了看热闹而是能直接映射到绝缘评估、电晕放电等方面。从仿真角度还有一个非常实际的好处棒板沿公共轴线对称用二维轴对称模型就能描述整个间隙计算成本比三维模型低一到两个量级。流注在发展初期也基本沿轴向传播二维轴对称下的结果已经能解释很多实验现象。当然流注到后期会出现螺旋和分支二维轴对称天然抹掉了角度方向的分支自由度这一点后面会专门说。但作为机理研究和参数扫掠二维轴对称棒板模型是性价比最高的起点。1.2 22反应体系要回答的问题流注放电的本质是电子在强电场中加速、碰撞电离、产生更多电子同时空间电荷畸变电场、反过来强化头部电离。这个过程里如果只关心“电子从哪里来”一个电离反应加一个附着反应就够了。但空气放电的实际情况要复杂得多电子能量不光用来电离还用来激发氮气和氧气的分子态这些激发态又会通过退激发、潘宁电离、能量转移等方式影响电子密度和活性产物。选择 22 个反应是一个典型的“精度与成本平衡”决定。反应太少比如只用五六步反应会丢激发态能量通道和重粒子化学流注传播速度、发光特性、臭氧浓度都算不准反应太多比如拉一份上百步的完整空气化学集非线性刚性和计算量会急剧上升尤其是在纳秒级瞬态里网格加密后每一步都要面对几十个输运方程的强耦合调试周期会拖得特别长。22 步这个规模刚好能把流注头部最关键的电子碰撞电离、附着、激发、退激发都覆盖住同时把尾巴上氮氧化物和臭氧的生成路径也搭起来是一个“既能跑出好看的流注又能拿到活性粒子数据”的折中点。2. 22反应体系怎么搭等离子体化学的平衡2.1 反应分类与代表性反应我在搭建这个模型时把 22 个反应按物理作用切成了八组。分组不是为了好看而是为了后面排查问题电子密度崩了先查电离和附着发光对不上先查激发和退激发产物浓度不对先查重粒子反应每一类问题都能快速定位到对应反应组。反应类别数量作用代表性反应弹性散射2电子动量交换影响迁移率与扩散系数e N2 → e N2e O2 → e O2电子碰撞激发5分子电子态激发能量损失和发光诊断e N2 → e N2(A3Σu)e O2 → e O2(a1Δg)电子碰撞电离与离解4电子雪崩核心二次电子来源e N2 → 2e N2e O2 → e O O电子附着2消耗自由电子形成负离子e O2 N2 → O2- N2e O2 O2 → O2- O2超弹性退激发2激发态将能量还给电子影响通道尾部电子温度e N2(A3Σu) → e N2离子-离子复合2放电通道内电荷中和尾部衰减N2 O2- → N2 O2O2 O2- → 2O2电荷转移与潘宁电离2离子物种转化、额外电子源N2 O2 → N2 O2N2(A3Σu) O2 → N2 O2 e重粒子生成与臭氧路径3臭氧、氮氧化物等活性产物的生成与损耗O O2 M → O3 MO O3 → 2O2这个表里电离、激发、附着、超弹性退激发是最核心的“电子动力学模块”它们决定了流注头部的电子增长和传播速度。重粒子生成与臭氧路径是“产物输出模块”它们对流注本身的传播影响很小但决定了放电产生了多少臭氧和活性氮是后续环境污染控制和等离子体应用研究最关心的量。2.2 为什么电离、附着和激发三者必须同时存在流注能不能传播本质上由三个竞争的电子碰撞过程决定电离产生自由电子、附着消耗自由电子、激发抢走电子能量。空气里氧气是附着性很强的气体氧气分子在低能区对电子有较大的附着截面很容易形成 O2-跟纯氮气放电差别巨大。如果你在模型里把氧气附着反应删了流注传播速度会明显偏快通道电流也偏高。所以做空气放电电子碰撞离解和附着这两个反应组必须保留不能偷懒。激发态的作用经常被简化模型忽略但实验上最常用的可是激发态发光。比如氮气 C 态跃迁回 B 态时发出 337 nm 附近的紫外光是 ICCD 相机拍流注照片时最常采集的信号。在仿真里电子碰撞激发反应 N2(C3Πu) 的速率空间分布可以近似代表发光强度的分布。这样就能把仿真云图和实验照片直接对照验证模型算出来的流注形态对不对。这也是为什么我坚持在 22 反应里保留 5 个激发反应而不是只留一两个。超弹性退激发则经常被不懂等离子体化学的人漏掉。N2(A3Σu) 这类长寿命激发态能量可以通过碰撞“给回”低能电子相当于一种加热机制。在流注通道尾部电场已经很低了但电子温度并没有立刻掉下来很大程度就是超弹性碰撞的功劳。如果你希望模型能算流注拖尾期间的电子输运系数这两步反应不能省。2.3 重粒子反应和活性产物怎么接进来流注通道里产生的大量 O 原子主要来自氧气分子的电子碰撞解离。有了 O 原子才能通过三体反应 O O2 M → O3 M 生成臭氧。这个三体反应在 1 个大气压下非常重要因为第三体 M 的浓度很高反应进行得比想象中快得多。做空气放电、废水处理这类应用研究时臭氧浓度是必须输出的指标所以 22 反应集里专门保留了臭氧的生成和损耗路径。氮氧化物的路径则跟潘宁电离和电荷转移密切相关。N2 在大气压空气中寿命极短因为 N2 O2 → N2 O2 这一步电荷转移反应速率非常快几乎刚生成就转化成了 O2。潘宁电离 N2(A3Σu) O2 又提供了额外的 O2 电子来源。这些反应决定了正离子到底是 N2 为主还是 O2 为主也影响了后续 NO、NO2 等含氮产物的生成路径。模型里只设置了正负离子复合反应没有把精细的 NOx 基元反应铺开是因为 22 反应这个规模不适合追求完整的含氮化学。如果研究目标是氮氧化物产物的精确浓度后面需要单独扩展反应集但用来判断“放电产生了多少氮氧化物趋势”已经够用。3. COMSOL等离子体模块建模从几何到边界3.1 几何尺寸和物理场接口怎么选几何上我用的是二维轴对称模型棒电极从计算域顶部中心向下伸出棒半径 0.5 mm棒尖与板电极之间的距离设为 10 mm。板电极就是计算域的最底边设置为接地电位。计算域不需要把左右两侧拉得很远流注主要活动区域在轴线周围径向设置 4 mm 左右就能保证边界效应不明显。实际图里计算域是一个宽度 4 mm、高度略大于 10 mm 的矩形区域棒电极位于上部间隙部分就是空气域。物理场选择上直接使用 COMSOL 的“等离子体”模块里的漂移-扩散接口它会自动耦合静电方程。漂移-扩散近似对大气压流注来说是在精度和计算量之间最合理的选项。完整的 Boltzmann 方程输运求解非常贵而且在大气压条件下流注头部尺度远大于电子平均自由程漂移-扩散近似的适用性是可以接受的。特别要注意的是一定要开启电子能量平衡方程让电子温度作为一个独立变量参与计算。这样才不是纯局域场近似才能算出发光强度分布。3.2 反应速率数据和输运参数怎么处理COMSOL 的等离子体接口支持两种方式输入电子碰撞反应数据一种是直接给反应速率常数随电子温度变化的曲线或表格另一种是给碰撞截面由软件利用两项近似求解电子能量分布函数。对于 22 反应模型我更推荐后一种因为空气碰撞截面数据库已经很成熟电子迁移率、扩散系数、电离速率这些量都能从截面数据自洽推出来不容易出现“电离速率和迁移率来自不同文献、对不上”的问题。输运参数里电子迁移率和扩散系数可以通过截面数据计算得到离子则按大气压条件下的典型迁移率来设定。N2、O2、O2- 的迁移率直接填常数即可扩散系数用爱因斯坦关系估算。电场强度对迁移率的非线性影响在流注头部强场区比较明显所以重粒子的输运参数不需要做得很精细电子部分才是这个模型的灵魂。针对反应速率要提醒一句所有反应速率常数都必须保持单位一致。COMSOL 里浓度单位如果选成 mol/m³速率常数也要按摩尔浓度换算如果选成 m^-3那速率常数单位就是 m³/s、m⁶/s 这种形式。我建议直接用粒数密度 m^-3这样跟截面数据和文献值对起来最省事反应速率常数的数值量级也不会因为摩尔常数产生 10^-3 级别的偏移。3.3 边界条件和种子电子怎么给棒电极表面设置为高电压电极板电极接地。等离子体边界条件的核心处理是电极对电子是吸收边界如果没有二次电子发射就把二次发射系数设为零。这个假设对裸金属电极在常温下是合理的。如果后续研究电极材料对放电的影响再去调整二次发射系数比如考虑冷阴极发射效果时这个系数就会变成一个关键参数。初始条件是新手最容易忽略的地方。纯空气背景里自由电子浓度为 0数值上放电永远启动不了必须人为给一个“种子电子”分布。我在棒尖附近给了一个峰值 1e13 m^-3、高斯半径 200 μm 的电子密度分布同时在全域设置背景电离密度 1e11 m^-3。背景电离本质上是模拟光致电离给放电前方提供的“备用电荷”。种子的具体峰值可以调但数量级不能太小太小了电离发展太慢数值上也容易被数值噪声淹没。电极电压不建议直接给理想阶跃因为瞬间跳变的电场会带来初始时刻的数值冲击。我用的是上升沿 5 ns 的斜坡电压让电压从零线性爬升到目标值这样流注启动过程更真实求解器也好收得多。3.4 网格尺度和加密策略流注头部是电场最强、密度梯度最陡的地方网格不加密算出来就是一团模糊的波前甚至会出现密度振荡。对这个 10 mm 间隙的模型我在棒尖正下方的轴上把网格加密到 10 μm 左右流注传播路径上保持 20 μm 左右到计算域边缘再逐渐放松到 0.5 mm。为什么要这么密因为流注头部厚度在几十微米量级至少要有三五个网格点跨过这个尺度才能分辨头部结构。网格一粗头部等效厚度变大流注速度就会失真电场增强也体现不出来。时间和网格是绑在一起的。流注传播速度约在 0.1~1 mm/ns按 10 μm 网格算库朗条件要求时间步长至少在 0.01~0.1 ns 量级。再加上电离反应时间尺度极短初始步长我会直接设置成 1e-12 s然后交给求解器自动增长。总模拟时间设置在 100~200 ns这个时间基本够正流注穿过 10 mm 间隙又不至于让离子输运和中性粒子变化把计算拖到完全跑不动。4. 求解设置与参数清单照着填就行4.1 推荐的可复现参数表下面这套参数是我在多次调试后觉得“开箱即用”的一组典型值。具体实验电压和间隙要根据你的研究目标调整但数量级和设置思路可以直接用。参数数值备注间隙距离10 mm棒尖到板电极棒电极半径0.5 mm尖端控制电场增强倍数峰值电压20~30 kV正极性典型值负极性另调电压上升沿5 ns避免阶跃冲击气压1 atm标准大气压气体温度300 K背景气体密度约 2.5e25 m^-3N2:O2 比例79:21干空气近似种子电子峰值1e13 m^-3高斯分布半径 200 μm背景电子密度1e11 m^-3等效光致电离计算域半径4 mm轴对称总模拟时间100~200 ns流注穿隙时间相对容差1e-3瞬态求解绝对容差1e-6浓度单位 m^-3初始时间步长1e-12 s自动可调电压值怎么定流注起始本质上要求尖端附近电子雪崩发展到临界尺寸即电离增长积分达到十几到二十几。对于棒半径 0.5 mm、间隙 10 mm 的结构20 kV 以下往往只能看到局部电晕很难形成贯穿流注超过 30 kV 后流注发展太快容易快速击穿对观测和网格分辨都不利。所以推荐先从 25 kV 开始试观察流注能否在几十纳秒内形成稳定传播。4.2 求解器怎么选分离式还是全耦合等离子体模块的瞬态求解是一个强对流、强反应源项的刚性方程系统。我用的是全耦合牛顿求解打开“自动缩放”功能变量缩放优先考虑电子密度和电势因为这两个量的量级跨了好几个数量级缩放做不好非常影响收敛。矩阵求解器我习惯用 PARDISO对二维轴对称这个问题规模来说速度很快。如果全耦合一直难以收敛可以改成分离式求解先解静电再解输运交替进行。分离式每步算得快但时间步长往往会因为拆分误差被限制得很小。我自己在 22 反应模型上还是推荐全耦合因为它对物理反馈的处理更自洽长时间推进反而总耗时更低。研究步骤上我强烈建议先建一个“稳态”研究只求解静电方程拿到无空间电荷时的初始电场分布。然后建“瞬态”研究把稳态解作为初始值并开启等离子体反应源项。这样流注演化是从一个合理的电场基线起步的而不是从全零势场硬算。4.3 结果怎么读自己算完怎么验证对不对把时间步推进到几十纳秒第一件看的事情是电子密度云图分布。典型正流注应该是棒尖附近出现一个高密度核向外延伸出一条细通道通道最前沿是一个密度陡增的“弹头状”头部头部后面通道的电子密度维持在高位但电场很低。如果云图变成一大团鼓包或者沿边界乱串网格或边界条件一定有问题。第二件看电场和电势分布。流注通道内部应该接近等势体因为通道电导率高到可以把电位“短路”掉外部电压几乎全落在头部前方一个极窄的区域内形成电场增强。如果你看到电场还在电极尖端而通道内部电场迟迟不下降说明通道电导没有建立起来电子密度不够放电还没进入流注阶段。第三件是对照实验信号。把 N2(C3Πu) 激发反应速率画出来这个量的空间分布近似发光强度跟 ICCD 发光的锥形分布去比。再把电极电流积分出来看电流波形有没有几十到几百毫安量级、随时间快速上升的形态。如果计算电流和实验电流差出一个数量级先回头检查初始种子密度和电压上升沿这两个参数对流注发展影响最大。5. 常见问题与排查技巧流注模型最容易踩的坑5.1 电子密度变成负值怎么办负电子密度是流注仿真初期最容易遇到的问题原因基本逃不开三个网格太粗导致对流项数值过冲、时间步长过大导致反应源项刚性误差、容差设置不合理导致迭代未收敛。排查顺序也是先网格、再时间、后容差。网格方面把流注头部区域的网格从 20 μm 加密到 10 μm甚至 5 μm如果负密度明显减少就是分辨率问题。时间方面把最大时间步长限制在 0.01 ns强制别让求解器“跳步”。容差方面相对容差从 1e-3 收紧到 1e-4绝对容差从 1e-6 收紧到 1e-8。但要注意绝对容差太小会大大增加迭代次数建议先试 1e-7。如果以上都调过还是有负值检查反应速率常数里有没有把正向和反向方向写反尤其是三体附着反应 e O2 M → O2- M它的速率常数通常随 E/N 下降而增大在低场区很强。一个不合适的附着速率会在低密度区把电子“抽空”导致浓度越过零线变成负值。处理这个问题的长远办法是给输运方程开启单调上风格式牺牲一点数值扩散换来非负性。5.2 流注起不来或者速度慢得离谱流注起不来的第一嫌疑是种子电子密度不足。我见过有同学把初始电子密度设成 1e7 m^-3算了一整天流注始终不出现就是因为种子太弱电离雪崩需要指数级增长增长初期信号完全被数值噪声淹没。直接改成 1e13 m^-3 后几步就启动了。第二嫌疑是电压不够。判断方法是后处理里查看最大约化场强 E/N。在棒尖附近E/N 至少要达到 100 Td 以上N2 和 O2 的电离反应速率才会有明显贡献。如果最大 E/N 只有几十 Td说明电压偏低或者棒尖曲率半径偏大要么升压要么把尖端半径调小。第三嫌疑是缺少光致电离等效源项。正流注传播依赖前方区域的背景电荷没有背景电子流注头部推进会非常吃力甚至停在半途。均匀初始背景密度 1e11 m^-3 是为了兜住这个物理过程如果你改了背景密度之后速度还是慢可以适当提高到 1e12 m^-3。但注意这是工程近似最好用实验流注速度标定一下别为了追求速度乱拉高背景密度。5.3 求解器不收敛的治理顺序模型一开算就报“达到最大迭代次数”先别急着改求解器按下面这个顺序排查第一步检查几何和边界。有没有悬空域、有没有把电极设成求解域、有没有漏掉接地边界。这些低级问题会导致初始电位分布异常进而让整个瞬态从起点就发散。第二步检查网格质量。尖端处的最大单元质量如果低于 0.1几何奇点处电场计算会非常离谱。在尖端点附近画一个微小的倒角或者圆角过渡比强行细分网格更有效。第三步降低初始阻尼。在最开始的一两个时间步里把牛顿阻尼因子设到 0.01给求解器一个“软启动”。等离子体问题初始猜测往往离解很远阻尼太小会振荡阻尼太大又可能直接拒绝迭代。第四步还是不行就拆开来试。先把所有重粒子反应禁用只剩下电离、附着、激发这几步电子碰撞反应跑通流程。跑通了再逐组把反应加回来。每加一组反应对比一下电子密度空间分布有没有跳变。这样定位问题反应非常高效我从这个操作里受益很多。5.4 流注传播速度和导电通道电流对不上实验如果流注形态挺好看但电流波形和实验差距大多半是通道电导率算高了。大气压电弧温度很高真实流注通道里气体被加热、密度下降电导率会显著升高但漂移-扩散模型默认背景气体密度恒定没算气体升温所以通道电导和电流往往会偏低。这个限制在 22 反应模型里没办法直接解决属于简化假设的代价。另一个原因是通道内长寿命激发态的能量回馈。如果你把超弹性退激发反应删了通道尾部电子温度会掉得过快通道残余电导下降电流拖尾也会跟实验对不上。所以做电流波形对比时这两步退激发反应必须保留。想进一步修正可以后续耦合一个简单的气体加热方程把 N2 转动温度或平动温度随时间变化算进去再反过来修正背景密度这属于模型扩展的范畴但至少你知道偏差来自哪里。6. 模型能拿来做什么工程价值与扩展方向从工程角度看这个棒板流注模型能回答的问题很直接给定一个尖端缺陷电压多高会起流注、流注多久能穿过间隙、通道电流有多大。这对绝缘子局部放电风险评估、电缆终端内部气隙放电研究都有参考价值。你也可以把模型改一改成棒-棒或板-板结构几何调整很便宜关键是放电物理没变还是那套漂移扩散加等离子体化学。从应用学科看空气流注放电的活性产物是臭氧和含氮氧化物的来源这个模型算出的 O3、NO、NO2 空间分布可以用于放电等离子体水处理、空气净化器效率模拟、人体等离子体医学的剂量估算。虽然 22 反应对精密化学计量还不够但能提供一个大体可信的空间分布趋势。真要精确算产物浓度可以在这个模型基础上再扩展反应集和输运方程。几个我认为值得做的扩展方向一是显式加入光致电离源项把“背景密度”这个近似升级为从氧分子激发态辐射计算的非局部源更接近真实流注分叉物理二是加入介质阻挡层把棒电极或板电极上覆盖一层介质材料研究表面电荷积累对流注发展的抑制三是把几何改成三维来研究流注分支代价很高但对机理研究很有帮助四是耦合简单的气体加热方程把放电后期的热膨胀和通道电阻下降补上。最后分享几点实战经验我自己第一次搭这个模型时一口气把 22 个反应全塞进界面结果跑通之前先花了两周调收敛。后来学乖了先用一个 6 反应的最小集把流注跑出来确认网格、边界、种子电子都没问题然后逐步把激发、附着、重粒子反应一组一组加回去。每加一组就对比一次电子密度基线和流注速度这样一旦结果异常马上知道是哪组反应在捣乱。这个“渐进搭建”的习惯救了我很多次。另外一个容易踩的小坑是反应速率常数的数据来源。电子碰撞反应速率如果是从截面数据算的温度相关曲线通常没问题重粒子反应速率常数大多来自文献不同文献给出的数值可能差两三倍。臭氧生成和三体反应这种对第三体浓度敏感的反应一定要确认速率常数的浓度参照系是空气、N2 还是纯 O2弄错了臭氧浓度会差一个量级。最后建议算完记得把 N2 激发发光分布导出来跟实验照片叠图。我试过叠加之后你才会真正意识到云图和照片的亮度分布哪怕有一点偏差都能暴露模型里某个反应给多了或给少了。这个模型做到那一步才算真正“活”了。