做仿真这些年我陆陆续续被问过好多次关于负电晕放电模型的问题尤其是Trichel脉冲电流这个东西。说起来这个现象在气体放电里算是老前辈了1938年就有研究者观察到并做了描述但直到今天能用仿真把它稳定复现出来依然是不少人头疼的事。COMSOL里做负电晕放电仿真模型本身看着不难针-板电极、空气域、加上等离子体模块设置好边界条件就能跑但真正做起来你会发现Trichel脉冲这个电流波形特别“娇气”网格、时间步、边界条件稍微有点不合适脉冲就没了或者干脆发散不收敛。这篇文章我会把自己跑负电晕放电模型的经验完整梳理一遍重点放在Trichel脉冲电流的产生机制、COMSOL里的模型搭建细节以及我最常踩的几个坑上希望能帮正准备做这块仿真的朋友省点时间。1. 先搞清楚负电晕放电和Trichel脉冲到底是什么1.1 负电晕的物理图景与放电特征先说清楚咱们在仿真的对象。负电晕放电指的是高压电极接阴极、低压电极或者接地平板接阳极的放电构型典型结构就是一根曲率半径很小的针尖对着一个平面。针尖附近电场强度极高空气在这个局域强场区里被电离形成放电等离子体。这里有个关键点得先明白负电晕和正电晕的放电形态差别非常大。正电晕放电通常比较“温和”电流是连续的辉光均匀包围在针尖周围而负电晕的电流却有明显的周期性振荡这就是我们说的Trichel脉冲。这个脉冲的频率范围很宽从几十kHz到几MHz不等取决于电压、气压、针尖半径、电极间距这些参数。Trichel脉冲的波形特征是窄的尖峰电流脉冲上升沿极短只有纳秒量级下降沿稍长一点脉冲与脉冲之间有一个相对平静的间歇期。初始电压刚超过起晕电压时脉冲频率低、幅值大随着电压升高脉冲频率越来越高幅值反而会有所下降。如果电压继续升高到一定程度脉冲逐渐变得不规则最终过渡到连续放电模式。1.2 Trichel脉冲产生机制与关键物理过程Trichel脉冲的产生机制简单概括就是负空间电荷的“产生—屏蔽—消散”循环。针尖附近电离产生的电子很快被阳极吸走因为电子迁移率极高而正离子迁移率低得多滞留在针尖附近形成正空间电荷层这个层会增强针尖前面的电场促进电离。但真正给出Trichel机制核心解释的是负离子的贡献附着过程产生的负离子O₂⁻等会聚集在针尖下游它们减弱了针尖前方的电场抑制了电离过程导致电流降低。这里有个容易混淆的点需要单独强调。在负电晕放电过程中正离子虽然产生得多、迁移率低、在针尖附近形成正离子云但是正离子云反而会让针尖表面电场更高让放电增强而真正在电流脉冲下降和熄灭中起作用的是负离子。电子从针尖表面出来、被电场加速、与氧气分子发生三体附着形成负离子这些负离子向平板迁移得非常慢在针尖前方形成一个负空间电荷团直接削弱了针尖前的电离区电场从而把放电“掐灭”。之后负离子继续下行空间电荷效应逐渐减弱针尖电场恢复新一轮电离开始第二个脉冲产生。所以Trichel脉冲的重复频率本质上是由负空间电荷的产生、输运和消散时间决定的。这一点想在仿真里复现出来数值上必须能捕捉到负离子的产生、堆积和撤退这个过程尤其是在时间尺度上频率几十kHz对应的时间是几十微秒但单个脉冲的上升沿可能只有几纳秒这对时间步长控制提出了很苛刻的要求。1.3 为什么要用等离子体模块来做这件事COMSOL里做气体放电很多人第一反应是只用AC/DC模块做静电场或者电流场计算。这样当然能算出电场分布也能大致判断针尖电场是否达到起晕条件但你算不出电流的振荡。因为Tricheel脉冲本质上是粒子密度时空演化的结果不是静电问题。这时候就需要等离子体模块Plasma Module出场了。等离子体模块基于漂移-扩散近似把各种带电粒子的连续性方程、动量方程简化成漂移扩散通量、以及泊松方程耦合在一起求解。对于空气放电来说涉及的主要粒子包括电子e、正离子N₂⁺、O₂⁺负离子O₂⁻还可能包括激发态粒子但如果我们重点看Trichel脉冲其实电子、一种正离子、一种负离子这个简化模型就能跑出基本脉冲形态。完整的三维针-板放电模型不是不能做但计算量极大一般个人电脑很难承受。实际工作中绝大多数文献和工程实践采用的都是二维轴对称模型。针尖是回转体平板是无限大平面整个放电结构绕中轴线旋转对称用二维轴对称坐标系建模可以把三维问题降到二维网格数量和计算时间直接下降一到两个数量级。这是做负电晕仿真最优的建模方式没有之一。2. COMSOL模型构建从几何到物理条件的完整设置2.1 几何建模与电极尺寸选择几何建模听起来简单实际上尺寸选择直接决定了电场分布的准确性。针尖曲率半径是一个极端重要的参数Trichel脉冲的频率、幅值、起晕电压都对针尖半径极敏感。文献报道中常见的针尖半径在10微米到100微米这个范围仿真中通常用抛物线旋转体近似针尖形状。我自己的经验是针尖半径取50微米左右比较适合初学者原因有两个一是这个半径的电场强度不会高到让数值求解特别困难二是网格划分时的曲率约束相对容易满足。如果半径取得太小比如10微米针尖附近就需要非常细的网格才能捕捉到电场集中效应网格数量会急剧膨胀。电极间距方面针尖到平板的距离通常在3到10毫米之间。间距越小起晕电压越低但间距过小也会让放电空间很小负离子输运距离不够脉冲可能不明显。我通常取5毫米作为初始值这个距离下负离子的输运时间在几十微秒量级对应的Trichel脉冲频率落在几十到几百kHz的典型范围内。平板电极就简单了一个矩形域宽度取间距的2倍以上就够了比如10毫米宽这样平板边缘效应对中心区域的影响可以忽略。2.2 多物理场耦合与等离子体接口选择在COMSOL里负电晕放电需要“等离子体”接口和“静电”接口或者直接使用等离子体内置的电场耦合。具体步骤是这样的添加“等离子体Plasma”物理场选择“漂移扩散”近似时间相关研究然后需要确保“静电Electrostatics”接口被自动添加或者手动添加因为等离子体接口依赖静电接口求解泊松方程。空气放电的化学反应集合是模型的灵魂。最粗糙的模型也要包含电子碰撞电离、附着和复合这三类过程。以空气的五种粒子简化模型为例电离反应e N₂ → 2e N₂⁺对应电子碰撞电离系数α附着反应e O₂ O₂ → O₂⁻ O₂这是三体附着反应负离子生成的绝对主力复合反应e N₂⁺ → N₂e O₂⁺ → O₂以及正负离子复合O₂⁻ N₂⁺ → 产物要注意的是反应速率常数不是随便填的。在漂移扩散框架里电离和附着反应速率通常写成电子能量或者约化电场E/N的函数COMSOL内置了部分气体的反应数据但空气放电建议使用Boltzmann方程求解器比如Bolsig或者COMSOL内置的Boltzmann模块预先算出速率常数随E/N的变化然后导入为插值函数。这是仿真结果能否和实验对上的关键。2.3 边界条件与初始值设定边界条件这块我单独拎出来说因为大部分模型跑飞都是边界条件给错了。针尖电极是阴极设置为金属壁边界。电子在金属表面的行为比较复杂电场方向会决定二次电子发射的方向。在负电晕中电场从针尖指向等离子体也就是电子受到的电场力朝向针尖电子带负电受力方向与电场方向相反正离子被加速轰击针尖表面产生二次电子发射。常用的简化处理是固定二次电子发射系数γ例如取0.01~0.1之间。γ的取值对Trichel脉冲的形状和频率都有影响但实际仿真中很多人忽视了它的作用默认取零或者一个常数。我建议至少取0.05左右试试如果脉冲完全出不来调整γ是一个重要手段。平板电极是阳极设置为接地金属壁电子会在这里被吸收正离子在表面发生中和反应不再返回等离子体。这个边界处理相对简单。初始值设置上电子密度和离子密度的初始值通常在一个很低的背景水平比如1e10 m⁻³模拟环境里总是有微量的背景电离存在。如果初始值设成绝对零放电起不来数值上也会有除法除零的风险。电压加载方式也需要注意不要一开始就加载完整的直流电压最好从零开始线性爬升到目标值爬升时间取1微秒左右这样一个温和的激励能显著降低初始阶段的数值冲击。3. 网格、求解器与时间步决定Trichel脉冲能否出来的三驾马车3.1 针尖附近的网格细化策略我在这个项目上花时间最多的地方就是网格。可以说网格不合适模型其他部分再怎么设置都没用。负电晕放电的放电区集中在针尖前方极薄的区域内电离区厚度通常只有几微米到几十微米。这个区域的网格尺寸必须小到能解析电子密度的空间梯度。我的做法是针尖表面第一层网格厚度取0.5~1微米针尖曲率半径范围内至少分布10到20个网格点向外逐渐加粗用“边界层”网格控制针尖附近的法向分辨率而远离针尖的区域用自由三角形网格。放电区域的体网格大小也需要跟得上前面的边界层网格由于漂移扩散方程中离子的特征输运尺度远大于电子德拜长度和电离层的厚度通常的做法是在紧贴针尖的区域用非常细的网格然后设置一个最大单元增长率让网格尺寸缓慢膨胀。我常用的参数是紧邻针尖的圆环域内最大单元尺寸2~5微米向外过渡到30~50微米底部平板附近放宽到100微米以上。有一个我踩过的坑需要特别提醒网格不仅要在针尖表面细化还要在针尖到平板之间整个中轴区域保留足够的纵向分辨率。虽然电离和附着只发生在高压区附近但是负离子云会沿着轴线向平板运动这条输运路径上的网格如果太粗负离子云的“疏散”过程在数值上会被涂抹掉导致脉冲间歇期重现不了最终看到的是单个脉冲而不是一串均匀的脉冲序列。3.2 求解器配置与非线性控制COMSOL的等离子体漂移扩散模型求解的是一个强非线性、强耦合的偏微分方程组默认求解器设置通常不够用需要手动调整。时间步长控制建议改成“手动”或者“中间Intermediate”并设定一个最大时间步长经验值在0.1纳秒到1纳秒之间。这个看起来很苛刻因为总仿真时长远超几百微秒按这个时间步长推下去步数会上十万但这是没办法的事。Trichel脉冲的上升沿只有几纳秒你如果不强制时间步足够小求解器会“跳过”脉冲。不过不需要全程保持这么小的时间步放电初期和脉冲之间的间歇期控制器会自动把步长放大真正细密的步长集中在脉冲发生时。求解器方面直接求解器是必须的。MUMPS或PARDISO都可以我习惯用PARDISO多核并行效率更高内存占用也更稳定。非线性方法选“牛顿Newton”阻尼因子设成0.9左右可以改善强非线性方程的收敛性。等离子体中电子密度动态范围极大可能跨越10个数量级默认的代数求解器有时候会处理不过来如果出现“找不到更小的阻尼因子”这类报错试着把相对容差放宽到0.01或者反过来把绝对容差收紧到1e-8量级。3.3 初始瞬态与激励方式模型是从零场开始的。直接加载高压直流电压数值上会产生一个初始时刻的巨大位移电流很容易把求解器打崩。我建议在入口电压上添加一个平滑的斜坡函数例如U(t) U0 · (1 - exp(-t/τ))τ取0.1~0.5微秒这样电压在几微秒内平滑过渡到目标值。过程中放电会自然点燃而不是被数值扰动强行驱动起来。另一个更进阶的做法是先用稳态求解器算一个不含放电的纯静电场分布以此为初始条件再跑瞬态这样等离子体的点燃过程会更接近物理实际。不过实际操作上我一般不用这个直接从低背景密度加斜坡电压也能跑出正确的结果反而更省事。4. 结果提取与后处理把Trichel脉冲波形稳稳拿出来4.1 电流提取方法与积分技巧仿真跑完之后怎么把Trichel脉冲电流提取出来COMSOL里计算边界电流有多种方式这里有一个需要注意的坑如果直接输出针尖表面的总电流密度积分得到的电流包含位移电流分量。在脉冲上升沿极短的情况下位移电流对总电流的贡献不可忽略甚至可能掩盖传导电流的真实形态。比较稳妥的做法是取针尖电极表面“总电流密度”变量COMSOL中通常会提供类似ec.nJ这样的边界法向总电流密度包含位移电流也可以用ec.nJc提取传导电流密度。具体提取传导电流时还要注意所有带电粒子通量相加。为了和实验电流波形做对比总电流更合适因为实验测到的电流本身就是位移电流和传导电流的总和。电流波形提取后用“一维绘图组”绘制电流随时间的变化曲线然后在结果节点里使用“全局计算”求平均周期、脉冲最大值、基频等。频率分析可以用“快速傅里叶变换”节点直接对电流时间序列做FFT得到频谱峰值对应的频率就是Trichel脉冲重复频率。4.2 典型脉冲波形的解读一个成功的仿真电流波形应该是这样的开始阶段有一个非常高的初始尖峰这是点火阶段的瞬态响应随后经过短暂的调整进入周期性稳定的脉冲串。每个脉冲的宽度通常在几十纳秒量级脉冲峰值电流可达毫安级甚至更高取决于电压和几何尺寸而相邻脉冲的间隔在几微秒到十几微秒之间。查看电子密度和负离子密度的空间分布可以看到脉冲相位的关键特征脉冲上升期针尖前电子密度快速增长形成高密度的电子团脉冲下降期这个电子团下游形成一大团负离子云负离子云的密度峰值逐渐远离针尖间歇期负离子云继续向平板慢慢移动针尖前的电场逐步恢复直到下一次电离雪崩重新发生。这个循环过程如果能在动画中展示出来理解和汇报效果都很好。建议做两个图一个是电流波形图另一个是针尖到平板轴线上的电子密度、负离子密度和电场强度的分布随时间变化图几个关键时间点截出来放在一起对比。4.3 参数扫描与验证模型跑通后我建议做一个简单的参数扫描来验证模型物理上是否可靠。固定其他参数只改变施加电压观察脉冲频率和幅值的变化规律电压越高脉冲频率越高脉冲幅值可能会有略微下降的趋势。这个趋势和文献中实验观测是吻合的如果仿真结果出现频率随电压升高反而降低的情况那就要回头检查边界条件和反应系数了。如果手头有实验数据可以把仿真电流波形的周期和脉冲幅值和实验值做对比。不过在对不上数值的时候也不用太纠结因为反应速率、二次电子发射系数、几何尺寸的实际误差都会影响结果能对上一个数量级已经算是不错的模型了。重点是趋势一致、物理过程正确这比数值完全吻合更重要。5. 实操中常见的坑与排查思路5.1 等离子体点不着火电流始终为零这是最常见的故障。先看计算得到的电子密度如果电子密度一直停留在初始背景值说明电离反应没有发生。排查顺序是网格是否捕捉到了针尖附近的高场区电场强度是否达到击穿阈值空气中约3e6 V/m电离反应速率常数是否合理施加电压是否超过了起晕电压。针尖半径太大会导致表面电场不够强这是几何因素上最容易忽视的原因。举个例子间距5毫米电压5kV如果针尖半径是200微米表面电场可能只有1e6 V/m量级根本达不到起晕条件。这种情况不是模型算错而是电晕放电在物理上就不会发生。5.2 只出现一个尖峰无法形成重复脉冲这个问题通常出在负离子产生不足或负离子输运过快上。Trichel脉冲的本质是负空间电荷的周期性消长如果附着反应速率设置得太低负离子无法在针尖前方形成有效屏蔽放电就会一直持续下去看到的电流是一条直线而不是脉冲序列。另外一个隐蔽的原因是计算域的高度不够。如果平板距离针尖太近负离子云被阳极快速吸收来不及在间隙中形成“撤退”的过程脉冲也会消失。试试把电极间距加大到10毫米或者把附着反应速率提高1~2个数量级观察波形变化。5.3 发散与不收敛求解器发散在最开始阶段最容易发生。处理方法很直接时间步长最大值收紧、绝对容差降低、网格在针尖区域加密、电压斜坡时间拉长逐个尝试。还有一个我常用的技巧是先把二次电子发射系数设成0跑一个粗结果再逐步把这个系数加回去相当于用分段方式逼近完整物理模型。另外要注意COMSOL的等离子体模块在默认情况下会对归一化使用对数尺度有时候负的密度值会干扰求解。这时候把“密度比饱和度”相关的选项打开限制重粒子的密度下限能显著提高稳定性。具体选项在不同版本里名字略有差异但关键词都是“saturation”或者“logarithmic”。5.4 参数合理范围的参考下面给出一个我常用的初始参数表以针尖半径50微米、电极间距5毫米、空气间隙、电压5kV直流负极性为基准参数推荐初始值备注针尖半径50 μm越小电场越强起晕电压越低电极间距5 mm间距影响脉冲频率和输运时间施加电压5 kV负极性扫描范围可到3~8 kV电子背景密度1e10 m⁻³初始种子电子正负离子背景密度1e10 m⁻³等离子体电中性初始条件二次电子发射系数γ0.05调整范围0.01~0.1气体温度293.15 K影响迁移率和扩散系数气体压力1 atm常压条件电离反应速率由玻尔兹曼解得出的插值函数不要手工填常数附着反应速率同上O₂三体附着速率与E/N强相关这套参数跑下来一般能得到频率在几十到几百kHz的Trichel脉冲串电流幅值在微安到毫安量级。当然不同的反应数据来源、不同的网格策略会让数值有一些差异实属正常。经验之谈与后续方向做负电晕仿真这件事我最大的感受是物理理解比软件操作重要得多。COMSOL的界面和建模流程并不复杂真正难的是你心里有没有那幅物理图像——你知道电流的下降是负离子干的知道脉冲的频率是由输运时间决定的知道网格解析不了梯度的时候一切结果都是假的。有了这层认识模型里各种参数设置就不需要死记硬背而是能顺着物理图像推理出来。如果你只是想把Trichel脉冲跑出来作为演示用我上面说的五粒子简化模型就够了。如果你要做更细致的研究可以往模型里加入电子能量平衡方程非局域热平衡近似、更多的氧氮激发态反应、光电离效应甚至是三维非对称几何下的流注演化。进阶方向里还有一个很有意思的领域就是用负电晕放电做静电除尘器的数值优化电场、空间电荷和颗粒物充电过程耦合在一起物理更丰富发文章和做工程都很有价值。仿真这条路就是一步一个脚印走出来的先把一个现象复现出来再慢慢把它做得更像物理。
