做高压放电实验的人应该都对着高速相机拍下的负流注画面纠结过针尖明明只是一根几十微米曲率半径的钨针板上电压也就是负几千伏可放电通道居然能在几十纳秒里跨过几毫米间隙头部还会拐弯、分叉形状跟正流注完全是两回事。为了把“负极性流注为什么能持续前进”这个问题看清楚我在COMSOL里用等离子体模块搭了一个二维轴对称的针-板空气负流注放电模型把电场演化、空间电荷分布和电子密度增长整个跑了一遍。这个模型能帮的忙不少验证针尖附近的Laplace电场是否达到击穿阈值、观察负流注头部的空间电荷层如何自持推进、估算流注跨越间隙的速度以及为实验里的ICCD影像和放电电流提供数值参照。想上手气体放电模拟、或者已经被收敛性问题折磨到怀疑人生的朋友这篇内容可以直接当避坑手册用。1. 先弄懂负流注为何难模拟针-板结构里的两套物理1.1 负流注与正流注的电场增强机制差别很多人刚接触流注放电时以为正流注和负流注只是电压极性调个头模型里改个符号就行。这可能是模拟里第一个大坑。正流注的头部是正离子空间电荷区这些正离子会在头部前方叠加一个与原电场同向的附加电场于是在传播方向上形成“越走越强”的电场增强而负流注走的是另一条路它的头部前方是电子堆积形成的负空间电荷层真正推动流注前进的电场增强区很薄只存在于这个负电荷层的前沿。这个差别直接决定了初始化和边界条件的设置方式。做负流注模拟时种子电子从哪来是个绕不开的问题。正流注可以依赖气体内部的背景电离加上针尖前方局部增强的电场迅速启动负流注从阴极出发没有二次电子发射或者预置电离层的话初始阶段会非常难产甚至干脆放不出来。所以理解机理不只是理论问题它直接决定你往模型里填什么参数。从数值角度看负空间电荷层在传播过程中会形成极高的局部电场梯度。这个梯度区和正流注头部的正离子区相比厚度更薄、电子密度峰值更陡对网格分辨率的要求也就更高。很多人在正流注模型里跑得好好的网格策略一换到负流注就发散发疯根子就在这里。1.2 针-板几何为什么是负流注研究的标准配置针-板结构本质上是一个“电场不均匀度可控”的放电平台。针尖曲率半径决定了局域电场增强倍数板间距决定了流注加速成长的空间长度这两者刚好对应流注模拟里最关心的两个问题起始条件和发展距离。模拟里我常看到有人把针尖画得特别尖锐觉得越尖越好其实这不是纯粹几何问题。针尖曲率半径如果小到和电子平均自由程同一量级COMSOL里默认的流体近似就开始不靠谱了因为连续性方程描述的是大量粒子的统计行为这种尺度下局部场近似LFA的适用性会动摇。经典的经验是让最小特征尺寸保持在1微米以上并且远大于德拜长度这样漂移-扩散方程才站得住。另外针体本身不是无限细的。真实的针电极是一根有直径的圆柱针体侧面的电场虽然远低于针尖但会影响流注在近阴极区域的分布形态。建模时不能省掉这段针体否则初始Laplace场分布会偏掉一块后续流注轨迹也跟着偏。2. 从零搭几何轴对称针-板模型的关键尺寸与网格纪律2.1 几何参数怎么设在COMSOL里第一步就是选二维轴对称空间维度。三维针-板也能算但计算代价非常大而且负流注头部的薄层结构在三维网格下很难喂饱。只要保证几何、边界条件和激励都是轴对称的二维轴对称得到的结果在工程分析层面就足够了。几何参数我建议直接用全局参数定义后面扫参不用到处改。一套常见的起始参数大概是这样的参数含义典型取值R_tip针尖曲率半径20~100 μmR_rod针杆半径0.5~1 mmd针-板间隙3~8 mmR_dom计算域半径2~3倍板间距p气压1 atmT气体温度300 K针尖曲率半径直接用圆弧代替圆弧与针杆相切连接。这个圆弧决定了最大电场值的量级不能为了好画网格给它偷偷倒角加粗那是自欺欺人。板电极放在间隙底部作为接地边界。针电极顶面接激励电压。计算域外侧不是电极实体而是开放空间边界半径至少取到板间距的两倍以上让外围电场能够自然衰落到接近零避免人为边界反射干扰流注通道。2.2 网格只在针尖和流注路径上较劲网格策略在这类模型里几乎决定着生死。负流注头部是一个强梯度薄层电子密度在微米尺度内可以从背景值跳到10^19 m^-3量级网格不够细的话数值扩散会把薄层“抹”成一片模糊流注速度算出来虚高一倍都不奇怪。我的做法是在针尖表面和板之间的轴线上做一条矩形加密带。这条带子的宽度大约覆盖流注预期径向发展范围轴向网格尺寸控制在1微米左右。如果条件允许针尖附近的网格可以压到0.2~0.5微米。这个加密带不需要铺满整个计算域外围区域用自由三角形网格逐渐过渡就行。但要提醒一点不要刚开跑就直接全局极细化先跑一个中等网格看流注头部大概走哪条路径再针对性加密。负流注在针-板间隙里基本沿轴线走但这不代表它不会在中途出现径向膨胀头部通过区域的径向网格也要留好余量。边界的网格同样不能含糊阴极针尖表面需要加边界层网格来解析阴极鞘层。即使采用漂移-扩散近似而不是完全求解玻尔兹曼方程边界处的电子通量条件仍然会受到网格质量影响粗糙的边界网格会让二次电子发射的数值实现变得很脏。2.3 初始条件背景电离密度和种子电子的给法负流注模拟里初始条件不是随便填个0就行。真实空气里始终存在宇宙射线和背景辐射维持的微弱电离电子数密度通常在10^11~10^13 m^-3量级。这个背景电离在模型里扮演的是“点火底火”的角色当针尖附近电场达到击穿条件后背景电子被加速、碰撞电离电子密度像滚雪球一样涨起来。实际操作中我在全计算域设置均匀背景电子密度和正离子密度相等保持初始电中性典型值取1×10^13 m^-3。光有背景还不够为了让起始过程更快更稳可以在针尖前方一个很小的半球区域里叠加一个高斯分布的种子电子团峰值密度给到1×10^15 m^-3左右。这里有个经验值得写下来负流注模型对背景电离密度相当敏感背景给得太低流注启动延迟长甚至可能在整个模拟时间内都处于“勉强生长”状态给得太高流注起始位置会提前速度也会偏快。所以最好把背景电离密度作为一个扫描参数跑10^11、10^12、10^13三档看启动时间差多少再和实验对照定标。2.4 边界条件电极、间隙边缘的兜底选择边界条件配置直接对应实验里的电极系统。针电极本体是电位边界施加随时间变化的负高压波形。板电极是接地电位并且要选为电子和离子的收集边界这样后续可以直接积分板电流和实验电流波形做对照。计算域外围边界要合理选择。我一般用开放边界配合无限元域或者简单的高电位接地边界关键是让外围区域不要积累空间电荷干扰内部电场。如果计算域半径只比板间距大一丁点边界附近电场会明显被压缩变形流注走到接近边界时还会出现虚假的加速这个很多人容易忽略。阴极表面的二次电子发射参数也在这个阶段就要填好后面物理场章节再说。这其实是一个“几何建模时就要思考物理边界”的例子等物理场搭完再回头补边界条件往往已经忘了自己当初画的边界入口到底对应哪个实体现容易把发射项设错位置。3. 等离子体模块设置把空气放电化学反应精简到可算3.1 放电物理场接口选型与漂移-扩散方程COMSOL等离子体模块里的核心接口是“等离子体放电”Plasma, Discharge它默认采用漂移-扩散近似来描述电子和离子输运再用静电接口求解泊松方程。这个组合不是唯一选择但针对空气流注这种中等气压放电它兼顾了精度和计算量是当前最主流的做法。漂移-扩散近似的核心思想是带电粒子的运动由两部分构成一是电场作用下的整体漂移二是浓度梯度造成的扩散。电子密度和平均电子能量是两个独立求解的量电子输运系数迁移率、扩散系数不是常数而是通过局部约化场强E/N插值得到的。这一近似隐含了“电子动能分布能快速跟随电场变化”的假设在纳秒量级流注时间尺度下通常是成立的。我选择额外求解电子平均能量方程而不是直接用局域平衡近似LFA。两者差别在于LFA假设电子能量完全由当地电场决定算得快但会把流注头部能量弛豫过程抹掉求解能量方程虽然增加一个输运变量但对负流注这种强非平衡的区域是值得的。物理场设置上记得把电子输运系数切换成“从碰撞截面计算”的模式。如果只是用固定迁移率E/N变化范围很大时结果会偏差明显。气体温度对约化场强有影响最好在材料属性里定义成随空间可变的理想气体密度用300 K常温和1 atm气压计算背景数密度。3.2 空气简化反应集电离、附着、复合三条主线空气放电的完整化学反应集可以复杂到上百个反应但流注模拟里最核心的只有三条主线电子碰撞电离产生新电子、氧分子的电子附着消耗电子、正负离子的复合消除电荷。在这三条主线上我建议先把氮氧相关的激发态反应尽量压缩等基础模型跑通后再逐步加复杂。一个能用且好收敛的简化反应集大概是这样的反应类型作用e N2 → N2 2e电离主要电子增殖通道e O2 → O2 2e电离次要电离通道e O2 → O2-两体附着电子损耗通道e O2 N2 → O2- N2三体附着高压下附着主导通道e N2 → N2复合正离子消除e O2 → O2复合正离子消除反应速率常数从哪来是新手最容易卡壳的地方。最好的来源是权威数据库里发布的碰撞截面数据然后通过COMSOL的“反应”设置导入。如果你是第一次搭直接用模块内置的空气放电反应库起步也行。但有个细节不同的内置模型里“Air”的物种和反应集合差异很大有的只含氮氧电离有的含了激发态你要看清自己选的哪一套再对照文献确认电子能量损失方程里是否包含了对应的非弹性碰撞损耗项。物种定义上电子是必须的正离子N2、O2以及负离子O2-也要创建。负氧离子的存在对负流注很关键因为氧的附着效应会不断抽走自由电子抑制流注过度发展。如果不把附着反应放进去流注头部电子密度会虚高通道增长偏快。3.3 阴极二次电子发射负流注得以维持的钥匙我必须专门给这个参数写一节。负流注和正流注的模拟在边界上的最大区别就是阴极的二次电子发射。负流注的电子从阴极出发针尖表面的正离子轰击会释放出新的电子这些电子在强场区被立即加速并发生碰撞电离。关闭二次发射整个负流注起步会变得极其艰难很多情况下甚至根本形不成稳定的流注头部。在COMSOL放电接口里阴极边界的二次电子发射通常用发射系数γ来定义表示每个入射正离子在表面释放的平均电子数。这个值与电极材料、表面状态、气体种类都有关系工程上常见取值在0.01到0.1之间。我在模型里先取0.05做基准再扫描0.01和0.1看结果差异。仿真时二次发射带来的电子会作为边界源项注入到等离子体域内同时也带来能量边界条件。被发射出来的电子初始能量很小大约是几电子伏特COMSOL里可以在边界设置平均初始能量这个细节别漏掉。如果漏了发射电子在阴极附近就会携带异常高的能量导致起始阶段电场畸变。3.4 瞬态研究配置和激励电压波形研究类型选瞬态。电压波形和实验保持一致的思路是针电极从0开始在几个纳秒内斜坡下降到目标负高压然后保持。直接给一个阶跃负高压的话初始电场突变会在模型里甩出强烈的数值冲击电子密度容易出现负值或振荡。我用平滑斜坡0到-10 kV用5 ns之后维持恒定。具体电压幅值怎么定一个重要检验环节是先用稳态静电求解器算一次纯Laplace场看针尖表面最大电场是否在击穿阈值附近。1 atm空气中的击穿场强约在3 MV/m量级针尖表面的几何增强效应通常能把场放大到几十MV/m。如果静电场算出来针尖场强根本到不了这个水平就别急着跑瞬态先把几何和电压调对再说。施加电压的边界是针电极顶部的电位边界电压函数里用带平滑斜坡的表达式通过全局参数控制幅值和上升时间。我这里给的参考波形是一个起点实际实验如果用的是纳秒脉冲源还要考虑脉冲的上升沿过冲和脉宽限制这些都要在波形函数里如实体现。4. 求解与稳定性负流注模拟里的收敛驯服记录4.1 耦合求解器与时步策略等离子体放电接口和静电接口是强耦合的带电粒子分布改变了空间电场电场反过来又驱动粒子运动。COMSOL处理这种多物理场耦合时我一般用全耦合求解器并启用“按时间分步”的依赖。相比分离式迭代全耦合在负流注这种强反馈场景下更容易稳定代价是每步计算量稍大。时步控制是这个模型最容易出问题的地方。流注头部电子漂移速度可以到10^5~10^6 m/s量级配合细网格CFL条件给出的时间步长小得惊人。比如网格最小尺寸是1微米、电子迁移速度考虑高场下1×10^5 m/s理论时步上限大约1×10^-11 s量级实际还要更保守。COMSOL自动时步未必能识别到这么细的限制我通常在求解器设置里手动加上最大时步取值约在几×10^-12 s到1×10^-11 s之间。直接跑完整模拟时间可能很漫长。解决办法是分阶段先跑一个较短的时间窗口比如10 ns验证起步正常检查一次流注头部位置和形状确认无异常后再把总时间延伸到100 ns。每一步都保留结果帧后面提取流注轨迹要用。4.2 三个典型翻车现场和对应处置把我在这个模型上踩过的坑直接列出来比给一百个参数都管用。翻车现场一负流注死活启动不了电子密度在针尖附近纹丝不动。查下来通常是两个原因叠加背景电离密度给得太低阴极二次发射系数设成了0。处置办法是把背景电离提到1×10^13 m^-3以上同时给阴极边界指定γ0.05。注意看起始阶段的电子密度云图如果出现针尖表面一个很薄的高密度层但向外没有发展说明二次电子是有在发射的只是发射能量或边界网格出了问题。翻车现场二电子密度出现负值或沿轴锯齿形振荡。这几乎可以肯定是网格分辨率不够、时步太大或数值扩散综合作用的结果。处置顺序是先加密流注头部要经过的区域网格把局部尺寸从1微米降到0.5微米再把手动最大时步缩小到原来的四分之一重新跑。如果锯齿只出现在特定位置通常是那个位置的网格过渡太突然把加密带边界修平滑。翻车现场三流注传播速度虚高远远偏离文献或实验量级。常见原因是背景电离过高、施加电压过高或者忽略了附着反应。前者容易排查后者则容易被忽略因为流注头部电子密度很高时电离占主导附着好像不起作用。但流注尾部区域其实一直在靠附着削减自由电子附着缺席的话整个通道的导电性会被高估尾部电荷泄漏不了头部推进自然偏快。第五个容易踩到的问题是计算域边界反射。如果外围边界设置成导体接地且距离间隙很近流注发展到接近边界时电场会被边界“吸”一下速度突然加速。我把计算域半径从5 mm加大到15 mm后才稳定下来。这类边界效应问题在云图里往往很难一眼看出来最好沿轴线提取电场分布观察边界附近是否存在异常的场强抬升。4.3 性能与网格依赖评估不管结果显示得多漂亮都得做一遍网格无关性验证。我一般跑三套网格基础网格头部尺寸1微米、加密网格0.5微米、超细网格0.25微米。对比的是同一时刻下流注头部位置和最大电子密度。如果头部位置在不同网格间漂移超过10%说明网格还没收敛需要继续加密如果加密后结果变化在5%以内基础网格就可以接受。同时留意计算资源的分配。二维轴对称模型看起来小但因为时间步长极短总步数可能上万每一步的求解规模又不小跑一次完整参数扫描可能要数小时到一天。我通常会先跑短时间窗口、粗网格做参数筛选确定合理参数范围后再做精细计算省下的时间非常可观。5. 后处理判读什么才算“流注建立”5.1 电子密度与净空间电荷的演化特征流注有没有建立不是看“有没有放电”而是看空间电荷是否形成了自持的头部结构。判读时先在云图里看电子密度分布一簇从针尖向板电极方向伸展的高密度通道头部保持一个清晰的球形或蘑菇形边界这就是流注形态。负流注另一个判据是净空间电荷分布。在针尖附近和流注头部前方净空间电荷应为负电子堆积超前于正离子。而在流注通道内部往后地带净空间电荷应以正电荷为主这正是电子离去后遗留的正离子成分。这种正负交替的空间电荷格局如果清晰出现说明模型跑出了名副其实的负流注而不只是阴极附近一团随便扩散的电离云。如果云图里看到的是一整片电子密度均匀膨胀、没有明确头部边界那大概率是背景电离过高预电离直接把整个间隙“点燃”了流注形态被淹没这时候降低背景电离密度或者减小种子电子区域重新观察。5.2 提取流注头部速度和电场增强倍数流注速度是模拟结果和实验对照最常用的物理量。提取办法不复杂沿轴线做一条截线在不同时刻导出电子密度剖面找到电子密度峰值或固定阈值比如1×10^18 m^-3对应的轴向位置位置差除以时间差就是该时间段的平均速度。多个时刻连成一条位置-时间曲线后可以估算流注是匀速、加速还是减速。另一个值得提取的是头部最大电场与背景电场之比。在负流注发展过程中头部前方的电场增强区应该显著高于Laplace场背景值典型增强倍数可以从几倍到十几倍。具体数值依赖电压、间隙和气压。如果电场增强区完全不见说明空间电荷没有有效累加起来模型还是停留在“电子漂移电离”而不是“空间电荷主导的流注机制”。判断流注是否进入稳态传播阶段可以看头部速度是否趋于稳定。如果是这时的速度才有和实验对照的价值如果从头到尾一直在加速说明电压可能过高流注接近击穿贯通不再是典型的由空间电荷自持传播过程。5.3 从模拟回到实验的对照做法最后聊点接地气的。模拟结果不能只在软件里自嗨要和实验互相印证。ICCD相机拍到的发光层位置对应模拟中电子密度高且激发反应强烈的区域你可以在云图里把高电子密度等值面的位置和实验影像里的发光前锋对齐粗略验证传播速度。板电流的波形对照也很有价值。在板电极边界上对电子和离子通量做积分得到的就是放电电流它随时间的变化包含了起始阶段、传播阶段、到达板极阶段的丰富信息。很多文献里的电流波形都带有典型的阶梯状上升特征模拟电流如果形状对得上整个模型的置信度会大幅提升。我自己迭代这套模型时的体会是负流注模拟不是一次跑成功就结束的事它更像是在“模型自洽”和“实验合理性”两头反复校准的过程。参数灵敏度分析、网格无关性验证、和实验波形对照这三件事每一项都比单纯把云图调好看重要。先跑通一套基础反应集和初始条件再逐步把背景电离、二次发射系数、电压波形这些因素逐个扫过去积累起来的数据会远比一张单张云图更能说明问题。后面如果你在这个模型上遇到具体问题或者想知道某个反应项具体怎么导入、某个边界条件怎么设欢迎按你手头的参数和情况来聊。模拟里没有放之四海而皆准的参数但有通用的排查思路把思路理顺大部分坑都是可以绕过去的。