在COMSOL里跑通一个沿面放电模型我大概花了将近两周坦白讲难点根本不在物理概念而在于怎么把“沿面放电”这个物理图像准确翻译成软件的边界条件和求解设置。网上能直接抄的案例太少了大部分教程到针板放电路径或者介质阻挡放电就停了真要把放电过程放到绝缘子表面去算同时拿到可靠的电子数密度和电势分布很多帖子要么语焉不详、要么默认你懂一大堆等离子体物理背景。这篇文章把我实际调试过程中的完整思路、参数配置和踩坑记录整理出来重点讲清楚电子数密度和电势分布这两个量是怎么一步步获取的模型又该怎么验证适合正在用COMSOL做沿面放电、绝缘表面闪络或高压绝缘结构电场分析的工程师和研究生参考。1. 沿面放电模型在建什么两个输出量的物理意义1.1 与纯气隙放电的本质区别沿面放电这个词字面上看是气体放电但真正让它区别于纯气隙放电的是绝缘介质表面的参与。纯气隙放电的电场由电极几何和空间电荷决定而沿面放电多了一个与放电通道直接接触的固-气界面这个界面的介电常数、表面电荷积累能力和二次电子发射特性都会改变放电路径的发展。典型场景是高压开关设备里的环氧树脂绝缘子、GIS母线支撑件或者PCB板上爬电距离不够时发生的沿面闪络。绝缘子表面在切向电场作用下放电通道倾向于沿介质表面发展而不是直接击穿空气间隙。原因也很直白介质表面的电场畸变比纯空气间隙更强而且在放电过程中表面会积累带电粒子积累的电荷反过来会扭曲空间电场形成正反馈最终导致闪络。这个反馈回路恰恰是沿面放电模型与普通放电模型最大的区别点。你在COMSOL里做针板电极气体放电边界条件基本上就是两个电极加开放边界绝缘子一进去就必须额外考虑绝缘表面的电荷连续性、二次电子发射和介电边界的电场折射这也是很多从气隙放电转过来的人一开始最容易忽略的差异。1.2 电子数密度与电势分布为什么是“必取”输出做过放电实验的人都有体会实验里能测的往往是外加电压、放电电流、或者高速ICCD相机拍到的发光通道形态但电子数密度的空间分布基本没法直接测电势分布在放电过程中更是几乎不可能用探针去测因为探针一伸进去就会扰动放电过程本身。所以模型里取这两个量本质上是在做实验给不了的事情。电子数密度的重要性在于它就是放电发展程度最直接的量化指标。等离子体中电子密度高的区域意味着电离反应剧烈、放电通道的主干在哪儿流注头部发展到什么位置全部可以从电子数密度云图上读出来。电势分布则是用来判断空间电荷如何影响电场的。当放电进入流注阶段空间电荷造成的电场畸变往往比外加电场还要强只看初始静电场分布根本不能预测放电后续路径这个时候电势分布云图能直观告诉你电场增强区域往哪儿移动了、绝缘表面的电位反转出现在什么时段。这两个量一个反映放电强度一个反映电场结构结合在一起基本上就能把沿面放电的演化过程讲清楚了。后面的模型校验、合理性判断包括和实验对比也都是围绕这两个量来展开的。1.3 模型的适用边界需要提前说明的是我们这里讨论的是流体近似下的漂移-扩散模型它适用于环境压力下比如常压空气、氮气、氩气的弱电离等离子体。弱电离意味着电子密度远低于中性粒子密度电子与中性分子的碰撞非常频繁平均自由程远远小于放电通道的特征尺度这种情况下电子能量分布可以用局域近似的迁移率和扩散系数来描述。如果你的场景是低压条件或者特征尺寸小到微米级、平均自由程接近特征尺寸那流体模型就会失真需要考虑PIC/MCC这类动理学方法。不过对工程上绝大多数沿面放电问题漂移-扩散框架是够用的。我在COMSOL里建模型的时候默认了常压空气、300K、放电间隙从毫米到厘米这个尺度范围所有参数和结论都建立在这个前提上。2. 方程与参数漂移扩散框架下要填的几组关键数字2.1 核心方程组不是只有电子连续性方程和泊松方程COMSOL的等离子体模块本质上帮你把漂移扩散近似的偏微分方程组集成到了图形界面里。电子输运由连续性方程控制[ \frac{\partial n_e}{\partial t}\nabla\cdot(-D_e\nabla n_e - n_e\mu_e \mathbf{E}) S_e ]其中 (n_e) 是电子数密度(D_e) 是电子扩散系数(\mu_e) 是电子迁移率(\mathbf{E}) 是电场强度(S_e) 是由电离、附着、复合等化学反应产生的净源项。电势 (V) 则由泊松方程决定[ -\nabla\cdot(\varepsilon_0\varepsilon_r\nabla V)e(n_ - n_- - n_e) ]这里需要提醒一下对于沿面放电绝缘介质区内部的电势同样要参与求解所以几何里要包含绝缘子本体并且给不同的域分配不同的相对介电常数 (\varepsilon_r)空气区 (\varepsilon_r1)环氧树脂大概在3.5~4.5陶瓷在6~9之间。因为介电常数不连续绝缘表面的电场法向分量会发生跳变这正是沿面放电中表面电场增强的根源之一。如果你需要更精细地考虑电子温度的非局域效应还应该加上电子能量守恒方程。这个时候需要额外的变量 (n_\varepsilon n_e\bar{\varepsilon})也就是电子能量密度方程形式为[ \frac{\partial n_\varepsilon}{\partial t}\nabla\cdot(-D_\varepsilon\nabla n_\varepsilon - n_\varepsilon\mu_\varepsilon \mathbf{E}) -e\mathbf{\Gamma}_e\cdot\mathbf{E} - n_e\sum_k \bar{\varepsilon}_k k_k n_k ]实际工程中如果放电主要靠外加电压和空间电荷场来驱动且电子温度的空间分布相对均匀可以暂时不求解能量方程直接把输运系数设成约化场强 (E/N) 的函数这样能省去不少非线性迭代的麻烦。但代价是丢失了电子温度在放电通道内、特别是流注头部快速变化的信息这个取舍要做清楚。2.2 输运参数和反应系数怎么填电子迁移率和扩散系数从哪里来我通常用Bolsig这类玻尔兹曼方程求解软件先算好再以插值表的形式导入COMSOL也可以直接用等离子体模块自带的气体数据文件。常压空气下的典型量级是电子迁移率 (\mu_e)约 (4\times10^{-2}) m²/(V·s) 量级随约化场强变化电子扩散系数 (D_e)约 (0.05) m²/s 量级爱因斯坦关系 (D_e/\mu_e k_B T_e/e) 可以粗略自洽温度流注头部电子温度通常3~5 eV强场区更高如果用空气做放电介质最简化的反应集至少包含三类反应直接电离(e N_2 \rightarrow 2e N_2^)附着(e O_2 \rightarrow O_2^-)复合(e O_2^ \rightarrow O_2)这三个反应分别决定了电子产生、损耗和整体的倍增平衡。反应速率系数的数值强烈依赖平均电子能量在COMSOL里通常以表格形式输入“平均电子能量—反应速率系数”的曲线。如果你用的是氩气反应集更简单基本只需要考虑电离和两体/三体复合初次调试模型的时候用氩气比空气省很多求解时间。2.3 能量方程到底解不解我的建议是初始模型阶段先不解能量方程用常数或简单插值的输运参数把框架跑通确认边界条件和几何没有问题后再打开能量方程。理由是能量方程引入了额外的非线性项和对流项求解器收敛难度会明显上升一旦同时出了多个错误排查起来很痛苦。但有一个例外如果你关注的是流注在绝缘表面的加速或停滞能量方程就必须解因为电子温度决定了电离反应速率而电离速率决定了流注头部能不能持续维持高场强。我曾经试过只在模型里放一个简单的电离反应集、不加能量方程算出来的电子密度分布虽然也呈现通道状但流注头部不往前走问题就出在电离反应速率系数被设置成常数无法反映头部电子的非局域加热效应。加了能量方程并改用基于平均电子能量的反应速率插值后结果才合理。3. 几何和边界条件电子数密度分布可信度的源头3.1 电极-绝缘子的几何搭建与尺寸选择沿面放电模型的几何并不复杂但尺寸选择直接影响放电起始位置和通道形态。我常用的结构有两种针-板电极加平板绝缘子适用于研究放电起始和沿面发展机制棒-板或者板-板电极加绝缘子更接近实际开关设备的绝缘结构。一个典型参考尺寸针电极尖端曲率半径0.2~0.5mm针板间距5mm绝缘子厚度2~4mm放置在接地电极上。针尖到绝缘表面的距离可以调这个距离决定了电场的切向分量和法向分量的比例。如果针尖离绝缘表面太近放电几乎立刻贴表面发展不利于观察起始阶段的连续演化过程离得太远放电在空气间隙中发展沿面特征又不明显。这里有个比较实用的建议一开始建立几何时不要追求和工程图纸完全一致用简化几何把模型跑通再逐步逼近真实结构。COMSOL中从CAD导入复杂绝缘子模型时偶尔会报“转换为CAD内核时不支持的拓扑”这通常意味着模型里有细小的倒角、曲面缝合问题简化几何可以绕开这类麻烦。仿真研究本质上是先抓住主导物理再考虑工程细节。3.2 表面电荷累积处理沿面放电和纯气隙放电的边界条件最大的不同就在绝缘子表面。放电产生的带电粒子到达表面后不会消失而是形成表面电荷表面电荷会反过来改变电场分布。这个负反馈回路是沿面放电最核心的物理过程边界条件必须把表面电荷的积累过程算进去。在COMSOL中最简单可靠的做法是自己定义表面电荷密度变量 (\sigma_s)它的时间演化可以写成[ \frac{d\sigma_s}{dt} e\left(\sum_i Z_i \mathbf{\Gamma}_i\cdot\mathbf{n} - \mathbf{\Gamma}_e\cdot\mathbf{n}\right) ]其中 (\mathbf{n}) 是表面法向。这个ODE和等离子体方程组是双向耦合的因为 (\sigma_s) 会改变表面的电场边界条件从而影响后续到达表面的粒子通量。实际处理中我把 (\sigma_s) 作为额外变量加在边界上同时修改导体表面附近电场计算时的边界条件。另一个办法是利用等离子体模块自带的“表面电荷积累”边界条件。不过我对默认的符号约定有过踩坑经历建议还是在第一次跑模型时单独设置一个探针监视表面电荷密度确认电荷符号和积累速率在物理上是合理的。比如正离子主导的表面电荷应该让表面浮电位上升如果算出来反而下降那就是符号方向搞反了。3.3 二次电子发射和壁面损失绝缘子表面的二次电子发射对沿面放电形态影响很大尤其对于短间隙、高场强的工况。二次电子发射系数γ一般取值在0.01~0.1之间具体取决于材料、表面状态、离子能量和入射角度。COMSOL里可以把γ设成常数也可以用表达式让它随入射离子能量变化。做绝缘子表面闪络研究的时候报道里常用0.05这个典型数值如果你手头没有实验数据从0.05起步是比较稳妥的。壁面处的电子损失也很关键。COMSOL中等离子体模块默认在壁面施加热流边界条件对于电子连续性方程壁面处通常采用零通量或对称条件。对于绝缘表面我的处理方式是让电子通量等于“到达壁面的热运动通量减去二次电子发射产生的电子通量”公式虽然看起来啰嗦但物理上是正确的。我实际调试中发现不少收敛问题恰恰出在壁面条件上。比如把金属电极的边界条件复制到绝缘表面忘了设置表面电荷积累点导致计算到十几微秒时电场分布乱掉又比如二次电子发射项在强场区域产生非常陡峭的源项导致时间步长退化得非常小。遇到这类情况先把二次电子发射关掉看看是否只是数值问题再加入回来。4. 调参全记录网格、求解器和收敛问题的实战处理4.1 网格策略放电路径附近绝不能偷懒沿面放电模型跨尺度很明显放电通道的特征宽度可能是几十微米而整个计算域可能是几厘米。如果网格在整个域上均匀划分单元数量会爆炸。我的做法是先用粗略的均匀网格计算初始电场分布找清楚电场增强最强烈的区域通常是针尖电极附近的绝缘表面然后在这一带做局部加密。具体设置针尖附近最大单元尺寸控制在0.05mm以内放电通道可能经过的沿面区域设置一层边界层网格层数6~8层最薄层厚度0.02mm左右。绝缘子内部网格可以很粗因为内部的电场是准静态的没有空间电荷效应不需要那么高的分辨率。网格这块有一个很容易被忽视的点在求解过程中流注头部会移动如果初始网格分辨率只在针尖附近流注发展到中段以后那里的网格太粗会导致电子密度数值扩散看起来通道变宽、传播速度变快。解决方式要么是把沿面路径全部加密要么启用自适应网格重构。COMSOL的物理场控制网格在等离子体问题里往往不够聪明我倾向于手动在沿面区域设置尺寸域让网格沿着放电可能路径均匀分布。4.2 求解器选择与时间步长沿面放电方程组是高度刚性的原因在于电子输运和反应特征时间非常短纳秒甚至亚纳秒而放电的整体演化持续微秒到毫秒量级。这时候显式求解器基本没法用必须用隐式BDF求解器配合自适应时间步进。COMSOL中默认的时步设置基本可靠但要留意最大时间步长我通常限制在纳秒量级防止求解器为了缩短计算时间跨过放电起始阶段。相对容差设置上默认0.01往往不够电子数密度这种跨越多个数量级的变量建议把容差压到 (10^{-4}) 到 (10^{-5})不然算出来的密度云图会有一圈圈的数值振荡。当然代价是求解时间变长。我通常在调试阶段用0.01容差快速跑通流程确认物理合理后再把容差收紧做最终计算。4.3 两个最典型的收敛失败案例第一个案例电子数密度出现了负值。原因通常是初始电子密度设得太低加上非线性求解过程中数值振荡把密度压到了0以下。解决办法是给初始电子密度一个合理的底值比如 (10^{10}) m⁻³量级而不是从零开始同时在密度方程中启用对数格式。COMSOL里可以设置因变量的对数缩放这会让 (n_e) 始终为正代价是方程形式复杂一点但值得。第二个案例“未能求出一致初始值”这个报错在等离子体瞬态求解中非常常见。根本原因是初始值不能同时满足泊松方程和连续性方程特别是初始电场和初始电子密度分布不自洽。我的处理办法分三步先求解纯静电模块得到稳定的电势分布再把这个结果作为等离子体模块的初始值或者先固定一个极低的电子密度只求泊松方程收敛最后再激活全部反应源项做瞬态计算。用辅助扫描逐步升高外加电压也是一个屡试不爽的办法相当于在物理上给求解器一个连续的“冷启动”过程。5. 后处理与数据导出让结果变成能直接用的图5.1 云图可视化电子数密度必须用对数刻度电子数密度在放电通道内部可能达到 (10^{20}) m⁻³而背景区域可能低到 (10^{10}) 以下差四五个数量级。云图如果用线性刻度最终看到的只是一条窄窄的亮线周围全是深色根本看不出放电结构的细节。我出图的时候一定选对数刻度并把色标范围截断到有物理意义的区间比如 (10^{12}) 到 (10^{21}) m⁻³背景的极低密度区域用一种颜色统一显示就好。电势分布则用线性刻度云图加等势线叠加可以比较清楚看到电势降集中在放电通道头部。在COMSOL的二维绘图组里我的组合套路是表面图显示电势等势线叠加再把电子数密度用半透明的第二层表面图画上去这样就可以在一张图上同时看两个量的空间关联。5.2 一维截线数据提取与二次绘图很多论文需要沿绝缘表面的物理量分布曲线比如“沿面距离-电子数密度”和“沿面距离-电势”。在COMSOL里最简单的方法是画一条沿绝缘表面的截线然后用一维绘图组输出这个路径上的数据。这里有一个小技巧截线不要太贴紧表面。因为壁面附近数值上存在边界层效应电子数密度可能因为壁面损失条件出现非物理的低值取距表面0.1mm的平行线既能代表沿面放电通道内的真实值又不会被壁面振荡干扰。导出到CSV之后我习惯再用Origin或者Matplotlib重新画图把坐标轴、字体、线宽统一调整成论文标准直接在COMSOL里面出图虽然方便但排版自由度还是差一些。5.3 用探针追踪关键位置的动态过程云图只能看某一时刻的快照放电路径的动态演化过程还需要探针来定量追踪。我通常在针尖附近、绝缘表面每隔1mm等距放置几个探针记录电子数密度和电位的时间曲线。这些曲线能给出流注沿面传播的速度相邻探针密度峰值出现的时间差以及表面电位翻转的时序。比如我跑的一个典型的针板沿面放电模型里探针数据显示流注头部在经过绝缘表面距针尖2mm的位置时电子数密度在几个纳秒内从 (10^{16}) 冲到了 (10^{20}) m⁻³同时表面电位从正几十伏翻转到了负几百伏这个翻转点是表面电荷积累的量化体现。这些动态数据是云图给不了的却是论文和工程报告中非常有说服力的内容。6. 结果有没有算错几个必须做的自检6.1 量级自检拿到计算结果的第一件事不是看云图长得好不好看而是检查量级靠不靠谱。常压空气流注通道的电子数密度通常在 (10^{19}\sim10^{21}) m⁻³ 量级如果算出来只有 (10^{14})多半是反应速率系数单位错了或者初始密度太低如果算出来超过 (10^{22})电离反应集可能设置得过度或者忽略了附着损耗。电势分布自检相对简单导体表面的电位应该等于你设置的边界电压绝缘表面电位不能超过外加电压的极值如果出现超幅振荡网格分辨率或时间步长大概率有问题。还有一个小技巧检查放电前后初始静电场分布放电起始前模型应该退化成纯静电问题绝缘子表面的电场法向分量跳变应该符合 (\varepsilon_1 E_{n1} \varepsilon_2 E_{n2}) 的关系这一步就验证了介质边界的处理是否正确。6.2 典型报错排查清单报错/异常可能原因排查与解决未能求出一致初始值初始电场和电子密度不自洽先求解纯静电逐步升压辅助扫描达到最大牛顿迭代次数网格太粗、反应源项太刚硬加密放电路径网格关掉二次电子发射再跑电子密度负值振荡初始值过低、线性求解器过冲启用对数缩放提高初始密度到 (10^{10}) m⁻³时间步长退化到极小某一点上的源项或通量剧烈变化检查壁面条件清理几何尖角松弛边界条件几何转换支持问题导入CAD复杂拓扑失败简化几何删除细碎倒角用内置几何工具重建这张表里的内容都是我自己跑沿面放电模型时真实遇过的其中“时间步长退化”那一条排查起来最花时间因为问题常常不在很明显的位置而是在你以为无关紧要的几何尖角处。等离子体求解对网格质量极其敏感几何建模阶段花一分钟处理掉的尖角能在求解阶段节省几小时。6.3 与实验对照的基本思路最后仿真结果必须经得起实验的对照。最直观的是放电通道形态高速相机拍到的发光通道位置在仿真里对应电子数密度大于某个阈值比如 (10^{18}) m⁻³的区域两者应该对得上。其次是放电的发展速度实验上用ICCD不同延时的照片可以估计流注传播速度仿真输出探针曲线也可以算出来两个数量级一致就是很好的验证。我个人的体会是仿真和实验很少有完美重合的时候关键是偏差方向要可控。比如计算得到的流注传播速度比实验快20%~30%这个偏差可能是忽略表面吸附、水汽或者绝缘子表面粗糙度导致的但如果速度差了好几倍那更可能是模型本身有结构性问题而不是参数微调能解决的。这时候回头检查边界条件和反应集比继续调网格有意义得多。
在COMSOL里跑通一个沿面放电模型我大概花了将近两周坦白讲难点根本不在物理概念而在于怎么把“沿面放电”这个物理图像准确翻译成软件的边界条件和求解设置。网上能直接抄的案例太少了大部分教程到针板放电路径或者介质阻挡放电就停了真要把放电过程放到绝缘子表面去算同时拿到可靠的电子数密度和电势分布很多帖子要么语焉不详、要么默认你懂一大堆等离子体物理背景。这篇文章把我实际调试过程中的完整思路、参数配置和踩坑记录整理出来重点讲清楚电子数密度和电势分布这两个量是怎么一步步获取的模型又该怎么验证适合正在用COMSOL做沿面放电、绝缘表面闪络或高压绝缘结构电场分析的工程师和研究生参考。1. 沿面放电模型在建什么两个输出量的物理意义1.1 与纯气隙放电的本质区别沿面放电这个词字面上看是气体放电但真正让它区别于纯气隙放电的是绝缘介质表面的参与。纯气隙放电的电场由电极几何和空间电荷决定而沿面放电多了一个与放电通道直接接触的固-气界面这个界面的介电常数、表面电荷积累能力和二次电子发射特性都会改变放电路径的发展。典型场景是高压开关设备里的环氧树脂绝缘子、GIS母线支撑件或者PCB板上爬电距离不够时发生的沿面闪络。绝缘子表面在切向电场作用下放电通道倾向于沿介质表面发展而不是直接击穿空气间隙。原因也很直白介质表面的电场畸变比纯空气间隙更强而且在放电过程中表面会积累带电粒子积累的电荷反过来会扭曲空间电场形成正反馈最终导致闪络。这个反馈回路恰恰是沿面放电模型与普通放电模型最大的区别点。你在COMSOL里做针板电极气体放电边界条件基本上就是两个电极加开放边界绝缘子一进去就必须额外考虑绝缘表面的电荷连续性、二次电子发射和介电边界的电场折射这也是很多从气隙放电转过来的人一开始最容易忽略的差异。1.2 电子数密度与电势分布为什么是“必取”输出做过放电实验的人都有体会实验里能测的往往是外加电压、放电电流、或者高速ICCD相机拍到的发光通道形态但电子数密度的空间分布基本没法直接测电势分布在放电过程中更是几乎不可能用探针去测因为探针一伸进去就会扰动放电过程本身。所以模型里取这两个量本质上是在做实验给不了的事情。电子数密度的重要性在于它就是放电发展程度最直接的量化指标。等离子体中电子密度高的区域意味着电离反应剧烈、放电通道的主干在哪儿流注头部发展到什么位置全部可以从电子数密度云图上读出来。电势分布则是用来判断空间电荷如何影响电场的。当放电进入流注阶段空间电荷造成的电场畸变往往比外加电场还要强只看初始静电场分布根本不能预测放电后续路径这个时候电势分布云图能直观告诉你电场增强区域往哪儿移动了、绝缘表面的电位反转出现在什么时段。这两个量一个反映放电强度一个反映电场结构结合在一起基本上就能把沿面放电的演化过程讲清楚了。后面的模型校验、合理性判断包括和实验对比也都是围绕这两个量来展开的。1.3 模型的适用边界需要提前说明的是我们这里讨论的是流体近似下的漂移-扩散模型它适用于环境压力下比如常压空气、氮气、氩气的弱电离等离子体。弱电离意味着电子密度远低于中性粒子密度电子与中性分子的碰撞非常频繁平均自由程远远小于放电通道的特征尺度这种情况下电子能量分布可以用局域近似的迁移率和扩散系数来描述。如果你的场景是低压条件或者特征尺寸小到微米级、平均自由程接近特征尺寸那流体模型就会失真需要考虑PIC/MCC这类动理学方法。不过对工程上绝大多数沿面放电问题漂移-扩散框架是够用的。我在COMSOL里建模型的时候默认了常压空气、300K、放电间隙从毫米到厘米这个尺度范围所有参数和结论都建立在这个前提上。2. 方程与参数漂移扩散框架下要填的几组关键数字2.1 核心方程组不是只有电子连续性方程和泊松方程COMSOL的等离子体模块本质上帮你把漂移扩散近似的偏微分方程组集成到了图形界面里。电子输运由连续性方程控制[ \frac{\partial n_e}{\partial t}\nabla\cdot(-D_e\nabla n_e - n_e\mu_e \mathbf{E}) S_e ]其中 (n_e) 是电子数密度(D_e) 是电子扩散系数(\mu_e) 是电子迁移率(\mathbf{E}) 是电场强度(S_e) 是由电离、附着、复合等化学反应产生的净源项。电势 (V) 则由泊松方程决定[ -\nabla\cdot(\varepsilon_0\varepsilon_r\nabla V)e(n_ - n_- - n_e) ]这里需要提醒一下对于沿面放电绝缘介质区内部的电势同样要参与求解所以几何里要包含绝缘子本体并且给不同的域分配不同的相对介电常数 (\varepsilon_r)空气区 (\varepsilon_r1)环氧树脂大概在3.5~4.5陶瓷在6~9之间。因为介电常数不连续绝缘表面的电场法向分量会发生跳变这正是沿面放电中表面电场增强的根源之一。如果你需要更精细地考虑电子温度的非局域效应还应该加上电子能量守恒方程。这个时候需要额外的变量 (n_\varepsilon n_e\bar{\varepsilon})也就是电子能量密度方程形式为[ \frac{\partial n_\varepsilon}{\partial t}\nabla\cdot(-D_\varepsilon\nabla n_\varepsilon - n_\varepsilon\mu_\varepsilon \mathbf{E}) -e\mathbf{\Gamma}_e\cdot\mathbf{E} - n_e\sum_k \bar{\varepsilon}_k k_k n_k ]实际工程中如果放电主要靠外加电压和空间电荷场来驱动且电子温度的空间分布相对均匀可以暂时不求解能量方程直接把输运系数设成约化场强 (E/N) 的函数这样能省去不少非线性迭代的麻烦。但代价是丢失了电子温度在放电通道内、特别是流注头部快速变化的信息这个取舍要做清楚。2.2 输运参数和反应系数怎么填电子迁移率和扩散系数从哪里来我通常用Bolsig这类玻尔兹曼方程求解软件先算好再以插值表的形式导入COMSOL也可以直接用等离子体模块自带的气体数据文件。常压空气下的典型量级是电子迁移率 (\mu_e)约 (4\times10^{-2}) m²/(V·s) 量级随约化场强变化电子扩散系数 (D_e)约 (0.05) m²/s 量级爱因斯坦关系 (D_e/\mu_e k_B T_e/e) 可以粗略自洽温度流注头部电子温度通常3~5 eV强场区更高如果用空气做放电介质最简化的反应集至少包含三类反应直接电离(e N_2 \rightarrow 2e N_2^)附着(e O_2 \rightarrow O_2^-)复合(e O_2^ \rightarrow O_2)这三个反应分别决定了电子产生、损耗和整体的倍增平衡。反应速率系数的数值强烈依赖平均电子能量在COMSOL里通常以表格形式输入“平均电子能量—反应速率系数”的曲线。如果你用的是氩气反应集更简单基本只需要考虑电离和两体/三体复合初次调试模型的时候用氩气比空气省很多求解时间。2.3 能量方程到底解不解我的建议是初始模型阶段先不解能量方程用常数或简单插值的输运参数把框架跑通确认边界条件和几何没有问题后再打开能量方程。理由是能量方程引入了额外的非线性项和对流项求解器收敛难度会明显上升一旦同时出了多个错误排查起来很痛苦。但有一个例外如果你关注的是流注在绝缘表面的加速或停滞能量方程就必须解因为电子温度决定了电离反应速率而电离速率决定了流注头部能不能持续维持高场强。我曾经试过只在模型里放一个简单的电离反应集、不加能量方程算出来的电子密度分布虽然也呈现通道状但流注头部不往前走问题就出在电离反应速率系数被设置成常数无法反映头部电子的非局域加热效应。加了能量方程并改用基于平均电子能量的反应速率插值后结果才合理。3. 几何和边界条件电子数密度分布可信度的源头3.1 电极-绝缘子的几何搭建与尺寸选择沿面放电模型的几何并不复杂但尺寸选择直接影响放电起始位置和通道形态。我常用的结构有两种针-板电极加平板绝缘子适用于研究放电起始和沿面发展机制棒-板或者板-板电极加绝缘子更接近实际开关设备的绝缘结构。一个典型参考尺寸针电极尖端曲率半径0.2~0.5mm针板间距5mm绝缘子厚度2~4mm放置在接地电极上。针尖到绝缘表面的距离可以调这个距离决定了电场的切向分量和法向分量的比例。如果针尖离绝缘表面太近放电几乎立刻贴表面发展不利于观察起始阶段的连续演化过程离得太远放电在空气间隙中发展沿面特征又不明显。这里有个比较实用的建议一开始建立几何时不要追求和工程图纸完全一致用简化几何把模型跑通再逐步逼近真实结构。COMSOL中从CAD导入复杂绝缘子模型时偶尔会报“转换为CAD内核时不支持的拓扑”这通常意味着模型里有细小的倒角、曲面缝合问题简化几何可以绕开这类麻烦。仿真研究本质上是先抓住主导物理再考虑工程细节。3.2 表面电荷累积处理沿面放电和纯气隙放电的边界条件最大的不同就在绝缘子表面。放电产生的带电粒子到达表面后不会消失而是形成表面电荷表面电荷会反过来改变电场分布。这个负反馈回路是沿面放电最核心的物理过程边界条件必须把表面电荷的积累过程算进去。在COMSOL中最简单可靠的做法是自己定义表面电荷密度变量 (\sigma_s)它的时间演化可以写成[ \frac{d\sigma_s}{dt} e\left(\sum_i Z_i \mathbf{\Gamma}_i\cdot\mathbf{n} - \mathbf{\Gamma}_e\cdot\mathbf{n}\right) ]其中 (\mathbf{n}) 是表面法向。这个ODE和等离子体方程组是双向耦合的因为 (\sigma_s) 会改变表面的电场边界条件从而影响后续到达表面的粒子通量。实际处理中我把 (\sigma_s) 作为额外变量加在边界上同时修改导体表面附近电场计算时的边界条件。另一个办法是利用等离子体模块自带的“表面电荷积累”边界条件。不过我对默认的符号约定有过踩坑经历建议还是在第一次跑模型时单独设置一个探针监视表面电荷密度确认电荷符号和积累速率在物理上是合理的。比如正离子主导的表面电荷应该让表面浮电位上升如果算出来反而下降那就是符号方向搞反了。3.3 二次电子发射和壁面损失绝缘子表面的二次电子发射对沿面放电形态影响很大尤其对于短间隙、高场强的工况。二次电子发射系数γ一般取值在0.01~0.1之间具体取决于材料、表面状态、离子能量和入射角度。COMSOL里可以把γ设成常数也可以用表达式让它随入射离子能量变化。做绝缘子表面闪络研究的时候报道里常用0.05这个典型数值如果你手头没有实验数据从0.05起步是比较稳妥的。壁面处的电子损失也很关键。COMSOL中等离子体模块默认在壁面施加热流边界条件对于电子连续性方程壁面处通常采用零通量或对称条件。对于绝缘表面我的处理方式是让电子通量等于“到达壁面的热运动通量减去二次电子发射产生的电子通量”公式虽然看起来啰嗦但物理上是正确的。我实际调试中发现不少收敛问题恰恰出在壁面条件上。比如把金属电极的边界条件复制到绝缘表面忘了设置表面电荷积累点导致计算到十几微秒时电场分布乱掉又比如二次电子发射项在强场区域产生非常陡峭的源项导致时间步长退化得非常小。遇到这类情况先把二次电子发射关掉看看是否只是数值问题再加入回来。4. 调参全记录网格、求解器和收敛问题的实战处理4.1 网格策略放电路径附近绝不能偷懒沿面放电模型跨尺度很明显放电通道的特征宽度可能是几十微米而整个计算域可能是几厘米。如果网格在整个域上均匀划分单元数量会爆炸。我的做法是先用粗略的均匀网格计算初始电场分布找清楚电场增强最强烈的区域通常是针尖电极附近的绝缘表面然后在这一带做局部加密。具体设置针尖附近最大单元尺寸控制在0.05mm以内放电通道可能经过的沿面区域设置一层边界层网格层数6~8层最薄层厚度0.02mm左右。绝缘子内部网格可以很粗因为内部的电场是准静态的没有空间电荷效应不需要那么高的分辨率。网格这块有一个很容易被忽视的点在求解过程中流注头部会移动如果初始网格分辨率只在针尖附近流注发展到中段以后那里的网格太粗会导致电子密度数值扩散看起来通道变宽、传播速度变快。解决方式要么是把沿面路径全部加密要么启用自适应网格重构。COMSOL的物理场控制网格在等离子体问题里往往不够聪明我倾向于手动在沿面区域设置尺寸域让网格沿着放电可能路径均匀分布。4.2 求解器选择与时间步长沿面放电方程组是高度刚性的原因在于电子输运和反应特征时间非常短纳秒甚至亚纳秒而放电的整体演化持续微秒到毫秒量级。这时候显式求解器基本没法用必须用隐式BDF求解器配合自适应时间步进。COMSOL中默认的时步设置基本可靠但要留意最大时间步长我通常限制在纳秒量级防止求解器为了缩短计算时间跨过放电起始阶段。相对容差设置上默认0.01往往不够电子数密度这种跨越多个数量级的变量建议把容差压到 (10^{-4}) 到 (10^{-5})不然算出来的密度云图会有一圈圈的数值振荡。当然代价是求解时间变长。我通常在调试阶段用0.01容差快速跑通流程确认物理合理后再把容差收紧做最终计算。4.3 两个最典型的收敛失败案例第一个案例电子数密度出现了负值。原因通常是初始电子密度设得太低加上非线性求解过程中数值振荡把密度压到了0以下。解决办法是给初始电子密度一个合理的底值比如 (10^{10}) m⁻³量级而不是从零开始同时在密度方程中启用对数格式。COMSOL里可以设置因变量的对数缩放这会让 (n_e) 始终为正代价是方程形式复杂一点但值得。第二个案例“未能求出一致初始值”这个报错在等离子体瞬态求解中非常常见。根本原因是初始值不能同时满足泊松方程和连续性方程特别是初始电场和初始电子密度分布不自洽。我的处理办法分三步先求解纯静电模块得到稳定的电势分布再把这个结果作为等离子体模块的初始值或者先固定一个极低的电子密度只求泊松方程收敛最后再激活全部反应源项做瞬态计算。用辅助扫描逐步升高外加电压也是一个屡试不爽的办法相当于在物理上给求解器一个连续的“冷启动”过程。5. 后处理与数据导出让结果变成能直接用的图5.1 云图可视化电子数密度必须用对数刻度电子数密度在放电通道内部可能达到 (10^{20}) m⁻³而背景区域可能低到 (10^{10}) 以下差四五个数量级。云图如果用线性刻度最终看到的只是一条窄窄的亮线周围全是深色根本看不出放电结构的细节。我出图的时候一定选对数刻度并把色标范围截断到有物理意义的区间比如 (10^{12}) 到 (10^{21}) m⁻³背景的极低密度区域用一种颜色统一显示就好。电势分布则用线性刻度云图加等势线叠加可以比较清楚看到电势降集中在放电通道头部。在COMSOL的二维绘图组里我的组合套路是表面图显示电势等势线叠加再把电子数密度用半透明的第二层表面图画上去这样就可以在一张图上同时看两个量的空间关联。5.2 一维截线数据提取与二次绘图很多论文需要沿绝缘表面的物理量分布曲线比如“沿面距离-电子数密度”和“沿面距离-电势”。在COMSOL里最简单的方法是画一条沿绝缘表面的截线然后用一维绘图组输出这个路径上的数据。这里有一个小技巧截线不要太贴紧表面。因为壁面附近数值上存在边界层效应电子数密度可能因为壁面损失条件出现非物理的低值取距表面0.1mm的平行线既能代表沿面放电通道内的真实值又不会被壁面振荡干扰。导出到CSV之后我习惯再用Origin或者Matplotlib重新画图把坐标轴、字体、线宽统一调整成论文标准直接在COMSOL里面出图虽然方便但排版自由度还是差一些。5.3 用探针追踪关键位置的动态过程云图只能看某一时刻的快照放电路径的动态演化过程还需要探针来定量追踪。我通常在针尖附近、绝缘表面每隔1mm等距放置几个探针记录电子数密度和电位的时间曲线。这些曲线能给出流注沿面传播的速度相邻探针密度峰值出现的时间差以及表面电位翻转的时序。比如我跑的一个典型的针板沿面放电模型里探针数据显示流注头部在经过绝缘表面距针尖2mm的位置时电子数密度在几个纳秒内从 (10^{16}) 冲到了 (10^{20}) m⁻³同时表面电位从正几十伏翻转到了负几百伏这个翻转点是表面电荷积累的量化体现。这些动态数据是云图给不了的却是论文和工程报告中非常有说服力的内容。6. 结果有没有算错几个必须做的自检6.1 量级自检拿到计算结果的第一件事不是看云图长得好不好看而是检查量级靠不靠谱。常压空气流注通道的电子数密度通常在 (10^{19}\sim10^{21}) m⁻³ 量级如果算出来只有 (10^{14})多半是反应速率系数单位错了或者初始密度太低如果算出来超过 (10^{22})电离反应集可能设置得过度或者忽略了附着损耗。电势分布自检相对简单导体表面的电位应该等于你设置的边界电压绝缘表面电位不能超过外加电压的极值如果出现超幅振荡网格分辨率或时间步长大概率有问题。还有一个小技巧检查放电前后初始静电场分布放电起始前模型应该退化成纯静电问题绝缘子表面的电场法向分量跳变应该符合 (\varepsilon_1 E_{n1} \varepsilon_2 E_{n2}) 的关系这一步就验证了介质边界的处理是否正确。6.2 典型报错排查清单报错/异常可能原因排查与解决未能求出一致初始值初始电场和电子密度不自洽先求解纯静电逐步升压辅助扫描达到最大牛顿迭代次数网格太粗、反应源项太刚硬加密放电路径网格关掉二次电子发射再跑电子密度负值振荡初始值过低、线性求解器过冲启用对数缩放提高初始密度到 (10^{10}) m⁻³时间步长退化到极小某一点上的源项或通量剧烈变化检查壁面条件清理几何尖角松弛边界条件几何转换支持问题导入CAD复杂拓扑失败简化几何删除细碎倒角用内置几何工具重建这张表里的内容都是我自己跑沿面放电模型时真实遇过的其中“时间步长退化”那一条排查起来最花时间因为问题常常不在很明显的位置而是在你以为无关紧要的几何尖角处。等离子体求解对网格质量极其敏感几何建模阶段花一分钟处理掉的尖角能在求解阶段节省几小时。6.3 与实验对照的基本思路最后仿真结果必须经得起实验的对照。最直观的是放电通道形态高速相机拍到的发光通道位置在仿真里对应电子数密度大于某个阈值比如 (10^{18}) m⁻³的区域两者应该对得上。其次是放电的发展速度实验上用ICCD不同延时的照片可以估计流注传播速度仿真输出探针曲线也可以算出来两个数量级一致就是很好的验证。我个人的体会是仿真和实验很少有完美重合的时候关键是偏差方向要可控。比如计算得到的流注传播速度比实验快20%~30%这个偏差可能是忽略表面吸附、水汽或者绝缘子表面粗糙度导致的但如果速度差了好几倍那更可能是模型本身有结构性问题而不是参数微调能解决的。这时候回头检查边界条件和反应集比继续调网格有意义得多。
