前阵子做氩气常压流注放电的仿真被一个“粗通道”问题折腾了两周。COMSOL里Ar流注放电的模型大家不陌生但网上能查到的案例大多是小针板间隙里细丝状的流注真正面向毫米级粗通道、棒电极对平板电极的案例很少。正好这几天把Ar棒板粗通道流注放电从反应机理、几何建模到数值求解完整跑通了写一篇实操笔记给同样在COMSOL里跟等离子体仿真较劲的朋友一个参考。这篇文章会按我建模型的顺序讲先解释流注的物理背景再讲几何和边界条件然后重点拆解Ar反应体系怎么设置参数最后谈求解稳定和结果判读。如果你是刚接触气体放电仿真也能从里面找到可直接照抄的套路。1. 流注放电和“粗通道”到底是个什么状态1.1 先搞清楚棒板结构里的电场长什么样棒板电极是气体放电里最经典的强非均匀场结构。高压棒电极曲率半径小低压板电极是大平面两者之间的电场分布极不均匀电场线在棒头高度集中离棒头越远场强衰减越快。实际计算里有个很好用的近似公式对球头棒电极半径R对无限大平板间隙d棒头顶端电场强度大概是E_max ≈ 2V / [R * ln(1 4d/R)]这里V是加载电压。举个例子R0.5mmd10mmV20kV那么E_max算出来大约在440kV/cm量级而间隙平均场强只有2kV/cm。差了200倍。这个电场畸变是流注启动的根本条件高电场区域集中在棒头附近电子在这里获得足够能量碰撞电离产生初始等离子体团之后空间电荷场进一步畸变电场等离子体通道向下发展。需要提醒的是这里不能直接套用平行板的均匀场概念。棒板结构里外加电场和空间电荷产生的电场方向相反——通道内部的电场被屏蔽高场区被压缩到流注头部薄层里。这种正负电荷分离产生的自增强场才是流注能自持传播的核心机制。1.2 粗通道流注与细流注、先导的区别低压或短间隙下流注通道半径往往只有几十微米可以用细流注模型描述。但当我做高气压、长间隙、较大电流的情况时放电通道会明显变粗这就是标题里的“粗通道”。粗通道不是说形状被故意画大了而是物理行为发生了质变。细流注头部薄、空间电荷层强传播主要靠头部极高场强粗通道内部等离子体密度高、电导率高通道整体更接近准中性空间电荷层变厚头部半径不断扩张。在Ar里粗通道现象特别明显原因是Ar没有电子附着过程。空气或者SF6里电子容易被O2或含电负性气体捕获形成负离子把自由电子消耗掉Ar是单原子惰性气体几乎没有附着截面电子只能通过复合和扩散损失所以放电通道能维持很高的电导率通道不容易熄灭反而继续扩张变粗。这个性质对雷电先导研究、脉冲放电等离子体辅助燃烧、材料表面处理都有用。如果你之后想从细流注过渡到粗通道第一步就要意识到不能用太短的间隙几百微米那种出来全是细丝放电也不能在Ar里硬塞附着反应——物理上就不对。我做仿真时从紫外灯实验数据里看到Ar中常压流注通道半径可以到2~5mm量级传播速度在10^5~10^6 m/s范围。这个量级可以作为仿真结果合理的参考锚点。2. 在COMSOL里搭Ar棒板放电模型从几何到物理场2.1 几何简化为什么用二维轴对称而不是三维棒板结构天然绕中心轴对称直接建二维轴对称模型计算量大减。三维模型对Ar流注几乎没有必要除非你要模拟多个流注分枝、气压极低导致三维波动明显或者电极本身不对称。我的几何做法是这样的计算域用矩形域代表气体区域宽度取间隙长度的2~3倍。比如间隙10mm区域径向宽度至少20mm避免边界影响通道演化。棒电极简化为直径1mm的半圆头圆柱轴向长度取5mm左右伸入气体区域。板电极直接用模型底部边界代表不额外建实体。接地边界板电极处设V0。开放边界计算域外边界设定为对称/开放条件避免电荷堆积反射。这里有个很关键的细节气体区域的径向范围一定要够大。一开始我图省事区域半径只设为间隙长度的一半结果通道发展到一半时径向电场受到人工边界限制流注头部出现奇怪的扁平化。2.2 物理场选择等离子体模块还是自己手写方程COMSOL 6.x的Plasma模块Plasma, Time Explicit或Plasma, Electrostatic内置了多组分漂移扩散方程、泊松方程和反应源项。流注放电的典型尺寸是毫米级流注头部尺度远大于电子平均自由程常压下电子平均自由程约1~10微米所以漂移扩散近似完全够用不需要粒子模拟PIC。我在这个模型里用的是“The Plasma Module”的静电接口把电子、离子和中性粒子都作为稀薄组分处理。如果你手头没有Plasma模块理论上可以手动添加多个“系数型PDE”接口写漂移扩散方程再耦合静电接口。但我建议不要这么搞耦合项的符号和单位很容易出错光调试就够你喝一壶。等离子体物理场下默认包含三类方程泊松方程电场由电位梯度得到空间电荷密度由粒子密度差分决定。组分输运方程电子和离子的连续性方程包含迁移、扩散、反应产生和消耗。电子能量方程如果启用非局域能量平衡COMSOL会额外求解电子温度T_e。对于流注问题建议保留电子能量方程。Ar的激发阈值和电离阈值差近4.3eV局域场近似Local Field Approximation会把电子温度直接当成约化场强E/N的函数在流注头部这种强梯度区域会引入不小误差。2.3 边界条件和初始种子电子怎么处理边界条件的工作量主要在棒电极和壁面上高压棒电极设置电位V(t)我一般用1ns的线性上升沿从0升到目标电压避免瞬态冲击导致初始电场过强而炸数。板电极电位恒为0。壁面反应Ar放电中壁面主要发生表面复合可以设二次电子发射系数γ0.01~0.1。这个参数对流注传播影响不大因为流注主要靠气体电离维持但对放电熄灭阶段有影响。初始种子电子密度我取1e12~1e13 m^-3放在棒头附近的高斯分布标准偏差取0.2~0.5mm。初始密度太低模型会花很长时间产生第一个电子雪崩太高则会跳过流注成核阶段看不到完整的物理过程。一个很实用的经验先跑一个没有种子电子、只有背景电离的模型测试电压高于击穿阈值时的电场分布。再把种子电子放进去跑完整流注。分两步走可以快速排查几何和边界设置的问题。3. Ar的反应机理与“填表”工程3.1 一套能用的Ar简化反应集Ar反应集是整个模型里最容易踩坑的部分。大而全的反应集动辄几十个反应参数查起来费劲数值刚性也大。流注放电本身是微秒以下尺度的非平衡过程没必要追求热平衡态的完整化学反应。我用的这套简化反应集适合常压粗通道流注物种只有五种e、Ar、Ar*这里指4s亚稳态把两个相近亚稳态合并、Ar2、Ar2*。反应包含编号反应类型备注R1e Ar - e Ar弹性碰撞主要损失电子能量速率常数来自截面数据R2e Ar - e Ar*激发阈值11.5~11.7eVR3e Ar - e Ar e直接电离阈值15.76eVR4e Ar* - e Ar e逐级电离阈值约4.2eV对放电维持很关键R5Ar* Ar* - Ar Ar e潘宁电离亚稳态浓度高时是重要电离通道R6Ar 2Ar - Ar2 Ar三体离子转换决定Ar2占比R7e Ar2 - Ar* Ar复合通道熄灭阶段重要这个反应集剔除了大部分激发态的辐射跃迁和精细能级分裂在流注传播阶段误差可接受。3.2 速率常数、迁移率和扩散系数从哪里来Ar的电子碰撞截面数据可以从LXCat官网的Phelps数据库或者Itikawa数据库拿。下载到截面之后用BOLSIG求解玻尔兹曼方程得到约化场强E/N从10到1000 Td范围内的速率常数和迁移率。COMSOL里可以用插值函数把速率常数写成约化场强E/N的函数也可以拟合乘Arrhenius形式。我建议你把BOLSIG算出来的表格导出在COMSOL里用“插值”函数读进去这样比手搓Arrhenius公式稳得多也方便调参时对比。离子迁移率方面Ar在Ar中的约化迁移率大约在1.5~2.0 cm²/(V·s)量级Ar2略高一些。注意单位换算COMSOL使用SI单位需要把cm²/(V·s)乘以1e-4转成m²/(V·s)。扩散系数用爱因斯坦关系D μ * k_B * T / e估算离子温度取300K电子温度用求解出的T_e。这里我多说一句不要偷懒把所有离子的迁移率填成同一个值。Ar和Ar2的迁移率差异虽然不大但三体转换反应R6在高压下非常快如果迁移率设置不当离子成分比例会出错后期复合阶段的结果完全不可信。3.3 那些容易填错的单位和系数这类模型90%的数值发散和“假物理”结果根源都在单位。双体反应速率常数单位是cm³/s还是m³/sCOMSOL里如果以mol为单位速率常数是m³/(mol·s)如果用粒子数密度就是m³/s。我的经验是直接以粒子数密度为基础全过程使用“m³/s”别混用。三体反应速率常数单位是cm⁶/s换算成m⁶/s要乘以1e-12。一旦填错放电通道内部Ar2浓度会差几个数量级。能量阈值单位激发阈值11.5eV电离阈值15.76eV在电子能量方程里要把eV换算成J乘1.602e-19。还有一个常被忽略的点R2输入的反应速率如果是BOLSIG算出来的随E/N变化的数据导入COMSOL时横坐标E/N的单位通常是Td1 Td 1e-21 V·m²。如果你做插值函数时横坐标直接填了V/m数据点会整体错位放电电压响应会变得完全失真。4. 瞬态求解的数值稳定性控制4.1 时间尺度和求解区间怎么确定流注传播速度一般取5e5 m/s量级10mm间隙的穿越时间大概在20ns左右。但流注启动时的弧光发展、以及到达对极板后的电荷消散光靠20ns不够观察。我的做法是总时长设100ns输出步长0.5ns电压上升沿1ns。这样既能完整看到流注形成、传播、到达板电极又能覆盖通道扩张的早期阶段。如果目标是粗通道的后期热化或者收缩时间要拉长到微秒级模型就要加上气体加热和密度扰动方程已经不是纯电模型范畴了。4.2 网格细化与CFL限制Ar流注头部厚度在常压下约0.01~0.1mm这是电场梯度和电荷分离最剧烈的区域。网格必须在这个尺度内充足分辨。我在轴向从棒头到板电极方向用边界层网格棒头附近和流注路径上最小网格0.02mm其余区域最大网格0.5mm。径向同样在轴线附近加密。整体网格量大约几十万自由度在COMSOL里属于中等规模。网格太粗时流注头部会“糊”。如果发现电子密度云图沿轴向一片模糊没有清晰的头部前沿基本就是网格分辨率不够。这里有个特别管用的自检方法把网格尺寸减半看流注传播速度和通道半径是否变化超过5%。如果变化明显说明网格没收敛必须继续加密。时间步长由COMSOL自适应控制BDF方法能处理等离子体的刚性源项。不过我建议手工设一个最大步长限制在0.2ns避免头部高速移动时数据点跨过了整个头部薄层。4.3 求解器配置与收敛判断流注问题本质是强非线性耦合问题电荷分离产生极强电场梯度同时电离源项对电场呈指数敏感。这种问题我强烈建议用直接求解器比如PARDISO或MUMPS别用迭代求解器。迭代求解器在这种病态矩阵下很容易收敛慢甚至发散。相对容差设置1e-4即可太紧会让时间步长小到无法接受。绝对容差按各组分分设电子密度绝对容差设1e10 m^-3左右比初始电子密度低两个量级即可电位容差0.1V。看收敛性不能只看最后结果图。要养成习惯每跑几步看求解器日志里的时间步长序列。如果步长一路走低到1e-16量级说明出现了剧烈不稳定多半是反应源项或者网格出了问题。正常的时间步长应该在1e-13到1e-10之间波动。5. 结果里看什么放电通道、电场与速度5.1 电子密度演化与通道扩张跑通之后最有信息量的是电子密度云图。你会看到棒头的种子电子被强电场加速形成电子雪崩电子密度从初始的1e12 m^-3迅速升高到1e20 m^-3以上高密度区域呈一个个“气泡”状向下推进。粗通道流注的特征在云图上很明显等离子体通道带而不是细线。如果你看到的是一条宽度几百微米不到的细丝那说明你的条件还是细流注模式或者Ar反应集里电离速率设置偏低。5.2 流注头部速度和电场增强因子的估算取沿轴线的电子密度剖面找到密度最大值的位置或者密度上升到某个阈值的前沿位置记录它随时间的位移差分就是流注头部瞬时速度。正常情况下Ar常压粗通道流注速度在10^5~10^6 m/s。如果你的结果高到10^7 m/s量级十有八九是初始条件给的种子电子密度太高或者电压上升斜率太陡流注还没充分成核就被强行推出去了。电场增强因子也是重要判据取流注头部的最大电场值除以间隙平均场强V/d。细流注的增强因子可以到20~50粗通道通常在10~25之间。这个数字可以用来对比文献判断你的模型处于合理的放电模式区间。5.3 二维云图和沿轴分布的常用出图方式COMSOL里报告结果遵循三个原则电位和电场用2D云图色标范围要手动设避免被棒头的局部高场“拉爆”色标。电子密度沿轴线用一个一维截点图方便看头部前沿和通道内部准中性程度。放电电流随时间变化用全局积分I ∫J_z dS对高压电极边界做积分可以得到放电电流波形。粗通道流注的电流波形通常有快速上升沿和缓慢衰减尾巴对得上实验示波器基本可以确认物理行为正确。6. 仿真翻车现场常见问题与排查思路6.1 电子密度出现负值或震荡电子密度出现负值几乎是每个等离子体仿真新手的必经之路。原因通常是头部电子密度梯度太大而网格太粗或时间步长过大导致对流项过冲。先检查网格再检查最大时间步长这两步能解决九成问题。如果还不行给组分输运方程打开人工扩散或流线扩散稳定数值稳定参数取0.1左右即可太高会引入不可忽略的虚假扩散把流注头部磨平。6.2 流注速度高得离谱或根本不动流注速度异常偏向两个极端。速度极高检查是不是电压上升时间太短或者初始电离度太高。速度几乎为零大概率是Ar反应集里电离反应没生效——常见错误是R3的速率常数表只覆盖了低约化场强范围而流注头部实际E/N远超插值函数上界插值默认返回边界值或0。我吃过一次亏BOLSIG算的速率表上界只到500Td流注头部局部约化场强跑到800Td电离速率直接被插值函数钉在上边界值速度偏差巨大。解决方法是把截面数据在BOLSIG里算到1000~2000Td再导入。6.3 计算时间过长与内存爆炸粗通道流注因为通道半径大网格需求比细流注高但没必要全计算域线性加密。把网格细化的范围限制在轴心附近直径2~3倍通道半径的区域里外围改用粗网格自由度能降一半以上。另一个值得注意的优化是精简反应集。有时候觉得反应加得越多越准确但其实不少反应在流注传播阶段贡献极小。把R4和R5拿掉分别跑一次对比如果电子密度分布和速度都差不多就可以放心把次要反应删掉减少刚性源项数量求解器会舒服很多。跑流注仿真我现在的固定流程是先0.5mm间隙小模型验证反应集合理性再切到10mm粗通道模型做网格和时间收敛性测试最后才放开电压和间隙做大参数扫描。这套流程看着慢实际上比一股脑全尺寸跑完后发现机理错误要省出几个通宵。Ar粗通道流注的魅力就在于只要反应集对了、数值稳住了COMSOL云图上的那一团团蔓延的等离子体通道比任何公式都更能让人理解气体放电究竟是怎么“爬”过去的。
COMSOL氩气棒板粗通道流注放电仿真建模与数值求解实战
前阵子做氩气常压流注放电的仿真被一个“粗通道”问题折腾了两周。COMSOL里Ar流注放电的模型大家不陌生但网上能查到的案例大多是小针板间隙里细丝状的流注真正面向毫米级粗通道、棒电极对平板电极的案例很少。正好这几天把Ar棒板粗通道流注放电从反应机理、几何建模到数值求解完整跑通了写一篇实操笔记给同样在COMSOL里跟等离子体仿真较劲的朋友一个参考。这篇文章会按我建模型的顺序讲先解释流注的物理背景再讲几何和边界条件然后重点拆解Ar反应体系怎么设置参数最后谈求解稳定和结果判读。如果你是刚接触气体放电仿真也能从里面找到可直接照抄的套路。1. 流注放电和“粗通道”到底是个什么状态1.1 先搞清楚棒板结构里的电场长什么样棒板电极是气体放电里最经典的强非均匀场结构。高压棒电极曲率半径小低压板电极是大平面两者之间的电场分布极不均匀电场线在棒头高度集中离棒头越远场强衰减越快。实际计算里有个很好用的近似公式对球头棒电极半径R对无限大平板间隙d棒头顶端电场强度大概是E_max ≈ 2V / [R * ln(1 4d/R)]这里V是加载电压。举个例子R0.5mmd10mmV20kV那么E_max算出来大约在440kV/cm量级而间隙平均场强只有2kV/cm。差了200倍。这个电场畸变是流注启动的根本条件高电场区域集中在棒头附近电子在这里获得足够能量碰撞电离产生初始等离子体团之后空间电荷场进一步畸变电场等离子体通道向下发展。需要提醒的是这里不能直接套用平行板的均匀场概念。棒板结构里外加电场和空间电荷产生的电场方向相反——通道内部的电场被屏蔽高场区被压缩到流注头部薄层里。这种正负电荷分离产生的自增强场才是流注能自持传播的核心机制。1.2 粗通道流注与细流注、先导的区别低压或短间隙下流注通道半径往往只有几十微米可以用细流注模型描述。但当我做高气压、长间隙、较大电流的情况时放电通道会明显变粗这就是标题里的“粗通道”。粗通道不是说形状被故意画大了而是物理行为发生了质变。细流注头部薄、空间电荷层强传播主要靠头部极高场强粗通道内部等离子体密度高、电导率高通道整体更接近准中性空间电荷层变厚头部半径不断扩张。在Ar里粗通道现象特别明显原因是Ar没有电子附着过程。空气或者SF6里电子容易被O2或含电负性气体捕获形成负离子把自由电子消耗掉Ar是单原子惰性气体几乎没有附着截面电子只能通过复合和扩散损失所以放电通道能维持很高的电导率通道不容易熄灭反而继续扩张变粗。这个性质对雷电先导研究、脉冲放电等离子体辅助燃烧、材料表面处理都有用。如果你之后想从细流注过渡到粗通道第一步就要意识到不能用太短的间隙几百微米那种出来全是细丝放电也不能在Ar里硬塞附着反应——物理上就不对。我做仿真时从紫外灯实验数据里看到Ar中常压流注通道半径可以到2~5mm量级传播速度在10^5~10^6 m/s范围。这个量级可以作为仿真结果合理的参考锚点。2. 在COMSOL里搭Ar棒板放电模型从几何到物理场2.1 几何简化为什么用二维轴对称而不是三维棒板结构天然绕中心轴对称直接建二维轴对称模型计算量大减。三维模型对Ar流注几乎没有必要除非你要模拟多个流注分枝、气压极低导致三维波动明显或者电极本身不对称。我的几何做法是这样的计算域用矩形域代表气体区域宽度取间隙长度的2~3倍。比如间隙10mm区域径向宽度至少20mm避免边界影响通道演化。棒电极简化为直径1mm的半圆头圆柱轴向长度取5mm左右伸入气体区域。板电极直接用模型底部边界代表不额外建实体。接地边界板电极处设V0。开放边界计算域外边界设定为对称/开放条件避免电荷堆积反射。这里有个很关键的细节气体区域的径向范围一定要够大。一开始我图省事区域半径只设为间隙长度的一半结果通道发展到一半时径向电场受到人工边界限制流注头部出现奇怪的扁平化。2.2 物理场选择等离子体模块还是自己手写方程COMSOL 6.x的Plasma模块Plasma, Time Explicit或Plasma, Electrostatic内置了多组分漂移扩散方程、泊松方程和反应源项。流注放电的典型尺寸是毫米级流注头部尺度远大于电子平均自由程常压下电子平均自由程约1~10微米所以漂移扩散近似完全够用不需要粒子模拟PIC。我在这个模型里用的是“The Plasma Module”的静电接口把电子、离子和中性粒子都作为稀薄组分处理。如果你手头没有Plasma模块理论上可以手动添加多个“系数型PDE”接口写漂移扩散方程再耦合静电接口。但我建议不要这么搞耦合项的符号和单位很容易出错光调试就够你喝一壶。等离子体物理场下默认包含三类方程泊松方程电场由电位梯度得到空间电荷密度由粒子密度差分决定。组分输运方程电子和离子的连续性方程包含迁移、扩散、反应产生和消耗。电子能量方程如果启用非局域能量平衡COMSOL会额外求解电子温度T_e。对于流注问题建议保留电子能量方程。Ar的激发阈值和电离阈值差近4.3eV局域场近似Local Field Approximation会把电子温度直接当成约化场强E/N的函数在流注头部这种强梯度区域会引入不小误差。2.3 边界条件和初始种子电子怎么处理边界条件的工作量主要在棒电极和壁面上高压棒电极设置电位V(t)我一般用1ns的线性上升沿从0升到目标电压避免瞬态冲击导致初始电场过强而炸数。板电极电位恒为0。壁面反应Ar放电中壁面主要发生表面复合可以设二次电子发射系数γ0.01~0.1。这个参数对流注传播影响不大因为流注主要靠气体电离维持但对放电熄灭阶段有影响。初始种子电子密度我取1e12~1e13 m^-3放在棒头附近的高斯分布标准偏差取0.2~0.5mm。初始密度太低模型会花很长时间产生第一个电子雪崩太高则会跳过流注成核阶段看不到完整的物理过程。一个很实用的经验先跑一个没有种子电子、只有背景电离的模型测试电压高于击穿阈值时的电场分布。再把种子电子放进去跑完整流注。分两步走可以快速排查几何和边界设置的问题。3. Ar的反应机理与“填表”工程3.1 一套能用的Ar简化反应集Ar反应集是整个模型里最容易踩坑的部分。大而全的反应集动辄几十个反应参数查起来费劲数值刚性也大。流注放电本身是微秒以下尺度的非平衡过程没必要追求热平衡态的完整化学反应。我用的这套简化反应集适合常压粗通道流注物种只有五种e、Ar、Ar*这里指4s亚稳态把两个相近亚稳态合并、Ar2、Ar2*。反应包含编号反应类型备注R1e Ar - e Ar弹性碰撞主要损失电子能量速率常数来自截面数据R2e Ar - e Ar*激发阈值11.5~11.7eVR3e Ar - e Ar e直接电离阈值15.76eVR4e Ar* - e Ar e逐级电离阈值约4.2eV对放电维持很关键R5Ar* Ar* - Ar Ar e潘宁电离亚稳态浓度高时是重要电离通道R6Ar 2Ar - Ar2 Ar三体离子转换决定Ar2占比R7e Ar2 - Ar* Ar复合通道熄灭阶段重要这个反应集剔除了大部分激发态的辐射跃迁和精细能级分裂在流注传播阶段误差可接受。3.2 速率常数、迁移率和扩散系数从哪里来Ar的电子碰撞截面数据可以从LXCat官网的Phelps数据库或者Itikawa数据库拿。下载到截面之后用BOLSIG求解玻尔兹曼方程得到约化场强E/N从10到1000 Td范围内的速率常数和迁移率。COMSOL里可以用插值函数把速率常数写成约化场强E/N的函数也可以拟合乘Arrhenius形式。我建议你把BOLSIG算出来的表格导出在COMSOL里用“插值”函数读进去这样比手搓Arrhenius公式稳得多也方便调参时对比。离子迁移率方面Ar在Ar中的约化迁移率大约在1.5~2.0 cm²/(V·s)量级Ar2略高一些。注意单位换算COMSOL使用SI单位需要把cm²/(V·s)乘以1e-4转成m²/(V·s)。扩散系数用爱因斯坦关系D μ * k_B * T / e估算离子温度取300K电子温度用求解出的T_e。这里我多说一句不要偷懒把所有离子的迁移率填成同一个值。Ar和Ar2的迁移率差异虽然不大但三体转换反应R6在高压下非常快如果迁移率设置不当离子成分比例会出错后期复合阶段的结果完全不可信。3.3 那些容易填错的单位和系数这类模型90%的数值发散和“假物理”结果根源都在单位。双体反应速率常数单位是cm³/s还是m³/sCOMSOL里如果以mol为单位速率常数是m³/(mol·s)如果用粒子数密度就是m³/s。我的经验是直接以粒子数密度为基础全过程使用“m³/s”别混用。三体反应速率常数单位是cm⁶/s换算成m⁶/s要乘以1e-12。一旦填错放电通道内部Ar2浓度会差几个数量级。能量阈值单位激发阈值11.5eV电离阈值15.76eV在电子能量方程里要把eV换算成J乘1.602e-19。还有一个常被忽略的点R2输入的反应速率如果是BOLSIG算出来的随E/N变化的数据导入COMSOL时横坐标E/N的单位通常是Td1 Td 1e-21 V·m²。如果你做插值函数时横坐标直接填了V/m数据点会整体错位放电电压响应会变得完全失真。4. 瞬态求解的数值稳定性控制4.1 时间尺度和求解区间怎么确定流注传播速度一般取5e5 m/s量级10mm间隙的穿越时间大概在20ns左右。但流注启动时的弧光发展、以及到达对极板后的电荷消散光靠20ns不够观察。我的做法是总时长设100ns输出步长0.5ns电压上升沿1ns。这样既能完整看到流注形成、传播、到达板电极又能覆盖通道扩张的早期阶段。如果目标是粗通道的后期热化或者收缩时间要拉长到微秒级模型就要加上气体加热和密度扰动方程已经不是纯电模型范畴了。4.2 网格细化与CFL限制Ar流注头部厚度在常压下约0.01~0.1mm这是电场梯度和电荷分离最剧烈的区域。网格必须在这个尺度内充足分辨。我在轴向从棒头到板电极方向用边界层网格棒头附近和流注路径上最小网格0.02mm其余区域最大网格0.5mm。径向同样在轴线附近加密。整体网格量大约几十万自由度在COMSOL里属于中等规模。网格太粗时流注头部会“糊”。如果发现电子密度云图沿轴向一片模糊没有清晰的头部前沿基本就是网格分辨率不够。这里有个特别管用的自检方法把网格尺寸减半看流注传播速度和通道半径是否变化超过5%。如果变化明显说明网格没收敛必须继续加密。时间步长由COMSOL自适应控制BDF方法能处理等离子体的刚性源项。不过我建议手工设一个最大步长限制在0.2ns避免头部高速移动时数据点跨过了整个头部薄层。4.3 求解器配置与收敛判断流注问题本质是强非线性耦合问题电荷分离产生极强电场梯度同时电离源项对电场呈指数敏感。这种问题我强烈建议用直接求解器比如PARDISO或MUMPS别用迭代求解器。迭代求解器在这种病态矩阵下很容易收敛慢甚至发散。相对容差设置1e-4即可太紧会让时间步长小到无法接受。绝对容差按各组分分设电子密度绝对容差设1e10 m^-3左右比初始电子密度低两个量级即可电位容差0.1V。看收敛性不能只看最后结果图。要养成习惯每跑几步看求解器日志里的时间步长序列。如果步长一路走低到1e-16量级说明出现了剧烈不稳定多半是反应源项或者网格出了问题。正常的时间步长应该在1e-13到1e-10之间波动。5. 结果里看什么放电通道、电场与速度5.1 电子密度演化与通道扩张跑通之后最有信息量的是电子密度云图。你会看到棒头的种子电子被强电场加速形成电子雪崩电子密度从初始的1e12 m^-3迅速升高到1e20 m^-3以上高密度区域呈一个个“气泡”状向下推进。粗通道流注的特征在云图上很明显等离子体通道带而不是细线。如果你看到的是一条宽度几百微米不到的细丝那说明你的条件还是细流注模式或者Ar反应集里电离速率设置偏低。5.2 流注头部速度和电场增强因子的估算取沿轴线的电子密度剖面找到密度最大值的位置或者密度上升到某个阈值的前沿位置记录它随时间的位移差分就是流注头部瞬时速度。正常情况下Ar常压粗通道流注速度在10^5~10^6 m/s。如果你的结果高到10^7 m/s量级十有八九是初始条件给的种子电子密度太高或者电压上升斜率太陡流注还没充分成核就被强行推出去了。电场增强因子也是重要判据取流注头部的最大电场值除以间隙平均场强V/d。细流注的增强因子可以到20~50粗通道通常在10~25之间。这个数字可以用来对比文献判断你的模型处于合理的放电模式区间。5.3 二维云图和沿轴分布的常用出图方式COMSOL里报告结果遵循三个原则电位和电场用2D云图色标范围要手动设避免被棒头的局部高场“拉爆”色标。电子密度沿轴线用一个一维截点图方便看头部前沿和通道内部准中性程度。放电电流随时间变化用全局积分I ∫J_z dS对高压电极边界做积分可以得到放电电流波形。粗通道流注的电流波形通常有快速上升沿和缓慢衰减尾巴对得上实验示波器基本可以确认物理行为正确。6. 仿真翻车现场常见问题与排查思路6.1 电子密度出现负值或震荡电子密度出现负值几乎是每个等离子体仿真新手的必经之路。原因通常是头部电子密度梯度太大而网格太粗或时间步长过大导致对流项过冲。先检查网格再检查最大时间步长这两步能解决九成问题。如果还不行给组分输运方程打开人工扩散或流线扩散稳定数值稳定参数取0.1左右即可太高会引入不可忽略的虚假扩散把流注头部磨平。6.2 流注速度高得离谱或根本不动流注速度异常偏向两个极端。速度极高检查是不是电压上升时间太短或者初始电离度太高。速度几乎为零大概率是Ar反应集里电离反应没生效——常见错误是R3的速率常数表只覆盖了低约化场强范围而流注头部实际E/N远超插值函数上界插值默认返回边界值或0。我吃过一次亏BOLSIG算的速率表上界只到500Td流注头部局部约化场强跑到800Td电离速率直接被插值函数钉在上边界值速度偏差巨大。解决方法是把截面数据在BOLSIG里算到1000~2000Td再导入。6.3 计算时间过长与内存爆炸粗通道流注因为通道半径大网格需求比细流注高但没必要全计算域线性加密。把网格细化的范围限制在轴心附近直径2~3倍通道半径的区域里外围改用粗网格自由度能降一半以上。另一个值得注意的优化是精简反应集。有时候觉得反应加得越多越准确但其实不少反应在流注传播阶段贡献极小。把R4和R5拿掉分别跑一次对比如果电子密度分布和速度都差不多就可以放心把次要反应删掉减少刚性源项数量求解器会舒服很多。跑流注仿真我现在的固定流程是先0.5mm间隙小模型验证反应集合理性再切到10mm粗通道模型做网格和时间收敛性测试最后才放开电压和间隙做大参数扫描。这套流程看着慢实际上比一股脑全尺寸跑完后发现机理错误要省出几个通宵。Ar粗通道流注的魅力就在于只要反应集对了、数值稳住了COMSOL云图上的那一团团蔓延的等离子体通道比任何公式都更能让人理解气体放电究竟是怎么“爬”过去的。