最近在帮着调一套氩气介质阻挡放电射流模型电极是典型的双环结构两个金属环贴在石英管外壁上模型建在COMSOL里用的二维轴对称坐标。折腾下来最大的感触是这类等离子体射流仿真真正难的不在于把模型搭起来而在于让放电物理、流场和边界条件在同一个模型里自洽地活起来。这篇把完整建模过程、关键参数和踩过的坑都写出来给打算做双环电极射流仿真的朋友一个参考。COMSOL做等离子体射流最常见的是Ar气在石英管中流动双环电极提供外部电场激励工作气压为1个大气压。这种结构在低温等离子体灭菌、材料表面改性、生物医学处理里很常见好处是电极不接触被处理对象等离子体以射流形式从管口喷出温度低、浓度高、活性粒子丰富。另一个常被问到的点是为什么不直接用更简单的针-环电极因为针电极会扰动气体流场管中心静电势集中也容易过渡到弧光放电双环电极的电场更均匀径向分布也更接近实验里观察到的环状放电结构。仿真层面二维轴对称模型就是把三维圆柱结构在r-z平面内求解网格数量从几百万降到几万到十几万计算时间从几天缩到几小时是这类问题最合理的折中方案。标题里后半段没写完按国内做这类仿真常用的路线一般就是“采用等效电路法处理介质层与放电间隙的电容耦合”或者“采用等离子体流体模型漂移扩散近似”。我这套模型用的是后者也就是Plasma模块默认的建模框架把电子密度、平均电子能、离子密度作为场变量和静电场方程联立求解。这个选择的依据很直接流体模型能给出空间和时间分辨的电子密度、电场、电离速率分布这对分析射流长度、传播速度、活性粒子浓度来说足够用计算量也控制在个人工作站能承受的范围内。如果改用粒子模拟PIC-MCC二维轴对称下千万级宏粒子的推进代价实在太高射流整个传播过程要算几十微秒甚至上百微秒粒子的运动步长却在皮秒量级算完基本等于一次小规模超算中心的机时消耗。1. 双环电极射流模型的整体设计思路1.1 为什么选双环电极它的电场结构有什么特点双环电极也叫ring-ring电极就是在介质管外壁套两个金属环一个接高压和接地之间加交流激励另一个悬空或直接接地。两个环之间隔着一段距离这段管壁内壁表面就会形成切向电场分量这恰恰是射流能沿着气流方向喷出管口的关键。实验室里用铜箔或者铝箔缠绕在石英管外就能做成这种电极制作简单、重复性好不需要把金属电极伸进放电管内避免了电极材料和气体直接接触带来的污染问题对生物医学应用这个优点尤其关键。从电场分布上看双环电极的高场强区域主要集中在两个电极之间的管段和管口附近。高压电极边缘电场最强因为金属边缘曲率大电荷密度集中会产生局部的电场增强效应。当外加电压是正负交替的正弦波时每个半周期内的极性效应并不完全对称——正半周期电子向电极迁移负半周期相反最终宏观表现就是放电看起来在每个周期里都有两次脉冲一次是主放电一次是残余电荷造成的反向放电只不过第二个通常弱一些。仿真结果里需要熬过最初的几个周期才能看到这样的周期性规律如果只算一两个周期就停很容易下错误的结论。双环电极另一个优势是放电体积可控。两个环的宽度、间距和管径决定放电区的尺寸管口喷出的射流长度基本由气流速度、电压幅值、频率和气体种类共同决定。仿真模型里调哪个参数能影响射流长度这是很快会上手的但想算得准就必须把几何、放电物理和流场三者绑在一起看。1.2 二维轴对称模型的建立逻辑与坐标选取二维轴对称是一个很巧妙但也很容易用错的模型假设。它的前提是几何结构绕中心轴完全旋转对称且边界条件、初始条件在方位角方向上没有任何变化。对双环电极射流这种结构来说石英管、金属环、气流入口全都在x-y平面里绕同一根轴转对称所以取r-z坐标系只需要计算一个纵截面包含中心轴的场分布再按对称性旋转展开就能反映整个三维空间的物理状态。但要注意一个隐含假设等离子体在对称轴附近的行为也必须旋转对称。在实际射流中放电可以出现径向模不稳定性、螺旋波结构或由外场不均匀导致的方位向不对称。实验里经常能看到射流并不是完美的直杆状而是有一定倾斜或分叉。如果模型目标是复现这种三维不稳定性二维轴对称模型是做不到的它天生认为一切物理量关于哪一个θ都相同任何方位角方向上的对称破缺都表达不了。所以在本文场景里二维轴对称适用于结构稳定、气流稳定、电压幅值适中、射流基本呈线状的放电状态。坐标选取上我把z轴定义为石英管的轴向r轴为径向对称轴放在r0的位置。COMSOL里就是在“几何”中选择二维轴对称组件2D Axisymmetric然后直接画矩形和圆环截面就行。几何建模时管壁内气体区域、管壁介质区域、管外开放空间区域要分开建因为它们的材料属性、参与计算的物理场接口不同。1.3 “采用等……”方案流体模型还是等效电路模型标题那半句没写完但国内做等离子体射流仿真的一般路线可以从项目背景里猜个大概。频率在几十k赫兹到几百kHz范围内的介质阻挡放电激励波长远大于模型特征尺寸所以可以用准静电近似不需要完整求解Maxwell方程组只要解静电场泊松方程就可以。对放电间隙和介质阻挡层常见的两种做法是把介质层等效为电容叠加在外部回路上或者把介质层的电荷积累作为边界条件显式地算进去。我在模型里采用的是后者——用COMSOL Plasma模块自带的介质表面电荷积累边界条件。这块算得准不准直接决定了放电电压和电流波形对不对。介质阻挡放电本质上是电荷在介质表面积累形成内建电场反过来削弱外加电场最终熄灭放电。等效电路模型适合快速估算放电功率、电流幅值这类宏观量但拿不到空间分布也分析不了射流形态所以对标实验ICCD图像的空间结构时必须走流体模型全耦合路线。整体建模框架分三层先是流体计算给定入口流速算出稳态速度场然后把速度场作为已知条件输送给等离子体模块放电产生的焦耳热和气流加热影响暂且忽略最后等离子体模块内耦合电子输运、重离子输运、静电场、介质表面电荷四个物理过程。一句话概括先算流场再算放电双向耦合放后面优化。这套方案对大多数射流仿真需求都足够稳定跑起来也省心。2. 几何建模与参数设置里的门道2.1 等离子体模块方程与关键物理设置COMSOL Plasma模块采用的是漂移-扩散近似电子用双矩方程描述电子数密度方程和电子平均能方程。电子通量由迁移项和扩散项组成系数来自外加电场和背景气体参数电场则通过泊松方程求解。重离子正离子、负离子、激发态中性粒子的输运方程同样的是迁移-扩散型式。这个近似在低气压到常压的介质阻挡放电中适用得很广前提是粒子碰撞频繁、平均自由程远小于特征尺寸。放电区气体的离子化学反应在氩气里不算复杂但要算得准反应列表也不能太省。我的Ar反应集包含了以下几种电子碰撞激发与去激发Ar e - Ar* e阈值能量11.5 eV左右电子碰撞电离Ar e - Ar 2e阈值能量15.76 eV二体与三体复合反应Ar e - Ar hν以及Ar 2Ar - Ar2 Ar潘宁电离Ar* Ar* - Ar Ar e这个反应对常压氩气放电的点火和传播贡献很大电子与背景气体的弹性碰撞和激发碰撞截面数据可以从LXCat库导出再通过Boltzmann求解器计算电子迁移率、扩散系数和反应速率常数作为电场的函数然后在等离子体模块中用一个插值表直接调用。这么做比用恒定迁移率靠谱得多因为常压放电中平均电子能在110 eV之间变化迁移率随电场强度变化可能相差好几倍。2.2 双环电极几何建模与参数表以一个典型的氩气射流放电管为例我从实际项目里拿来的尺寸参数是这样的参数名称取值说明石英管内径4 mm放电气体通道石英管外径6 mm管壁厚度1 mm高压电极宽度5 mm铜环贴于管外壁接地电极宽度5 mm铜环贴于管外壁电极间距8 mm两环之间距离高压电极前缘到管口10 mm射流在管口外形成计算域长度60 mmz方向包含管外空间计算域径向半径20 mm管外开放空间电极间距值得多聊两句。间距太小放电集中在两个电极之间的狭窄区域射流很难有效往管口推送间距太大放电区离管口太远活性粒子在输运过程中损失严重到了管口基本衰减殆尽。实验上最优间距通常在515 mm模型里可以先按这个范围扫几个值。要注意的是电极间的电场是相对均匀的而管口外空间由于介质管端的几何突变电场分布变得非常不均匀所以射流在管口外是依靠电场梯度自持传播的。输出远场等离子体参数的准确性对网格密度极其敏感后面会细说。高压电极的驱动电压我给出的典型值是8 kV峰值、13.56 kHz正弦波这是RF频段的低频端也是常压等离子体射流实验常用参数。如果换用10 kHz甚至50 kHz射流形态和放电模式会发生明显变化所以参数表里的频率不是随便填的要和实验条件对应。2.3 边界条件设置的几个关键细节边界条件是这类仿真里最容易出错的地方也是最难排查的地方。逐条说高压电极表面设置为电压边界条件V0*sin(2*pi*f*t)V08000 Vf13.56 kHz。同时这个边界需要开启“金属接触”属性让电荷能在电极表面积聚并对放电产生影响。接地电极设V0。如果实际实验用的是浮动电位比如接一个大电阻到地那么模型里需要用“浮动电位”边界条件去描述悬浮电极的电位自洽变化两种设置得到的放电电流和电压电流相位差差别不小。介质表面设置“表面电荷积累”边界条件。这是模拟介质阻挡放电的核心。石英管壁的表面电荷密度σ作为额外变量由流过介质边界的光通量和总电流决定。我试过用等效电容法替代结果是电流波形里的窄脉冲消失只有平滑的位移电流和实验电流波形完全对不上。所以这一步千万别省。管口和外部开放边界管外空间边界设为接地或开放边界让电场自然收敛到远场零电位。注意不要让任何导电边界直接和高压电极靠得太近否则模型会放电到边界上产生一堆非物理现象。入口气体边界气体入口的等离子体密度设为本底值电子含量极低比如1e12 m^-3这是为了给放电提供初始种子电子。这个值的设定参考了宇宙射线和背景电离对大气放电的贡献实际实验管内气体源源不断初始密度对稳态周期放电的影响经过几个周期后基本会被抹平但太小会让第一个放电脉冲出现时间明显延后。出口边界等离子体组分离子或高能电子碰到开放边界时一般设置为流出漫反射或零梯度如果设置的边界条件和流场耦合不对会导致气体在管道内累积压力模型报错。3. 实战操作从网格划分到结果后处理3.1 网格划分的策略与必要性射流放电的物理特征尺度跨度非常大这是网格划分最难的地方。放电通道流注头部直径可能在百微米量级而整个计算域长度却是60 mm差了三个数量级。如果网格均匀加密到能分辨流注头部节点数直接爆炸如果不加密流注根本走不动或者扩散得太厉害。我的做法是分区域差异化划分放电区两个电极之间以及管口向外约10 mm、半径约5 mm的区域用最大单元尺寸0.1 mm的三角形网格靠近管壁内表面和电极边缘处再加两层边界层网格第一层厚度约0.02 mm增长率1.2。这是整个模型最要的地方。常压放电中的鞘层厚度在几十微米量级网格太粗相当于人为增大了扩散系数电子密度会被严重低估。管壁介质区域网格可以放松一点最大0.2 mm因为介质内部只算泊松方程没有等离子体反应。管外远场区域离管口超过15 mm、半径超过8 mm网格可以用0.51 mm这里放电已经变得很弱主要是捕获活性粒子输运趋势。统计下来总网格数大约8万到12万单元对二维轴对称模型来说非常轻松。COMSOL自动网格生成器给出的默认网格在这个问题里偏粗一定要手动调整。如果你跑完发现电子密度在管口沿轴向出现锯齿形振荡先怀疑网格而不是怀疑物理。3.2 求解器设置与时间步长控制等离子体放电模型的瞬态求解是硬仗因为物理过程横跨多个时间尺度电子在电场中响应时间是皮秒量级流注传播速度是10^410^5 m/s在60 mm空间走完也就微秒量级而流体对流的时间尺度是毫秒量级。所以时间步长控制是仿真的核心难点。我的做法是瞬态求解器BDF向后差分公式方法最大阶数2因为等离子体化学体系的刚性很强高阶BDF容易振荡发散。初始步长给8e-13 s是的皮秒级最大步长限制在1e-8 s。在每一个正弦周期13.56 kHz下约73.7 μs内至少采集200个点否则电流波形的窄脉冲是完全看不到的。最开始我试过默认的自适应步长结果COMSOL把步长拉到微秒级放电脉冲的时间分辨率完全丢失电流曲线几乎是一条直线。求解器选择PARDISO直接求解器。模型未知数不多二维轴对称下约20万自由度的非线性方程组直接求解器比迭代求解器稳定得多。等离子体模块中强非线性反应项频繁造成Jacobian矩阵异常迭代求解器很容易发散这也是为什么很多人在COMSOL中跑这种模型会莫名报错。模拟总时间我一般设置包含5到10个电压周期前两个周期用来让系统从初始状态过渡到准周期状态后面几个周期做统计和处理。如果只算2个周期放电电流幅值和相位可能与稳定状态差别很大做参数扫描时错误特别明显。3.3 结果处理与常用输出量输出的核心量包括电子密度分布这是判断放电是否成功起辉、射流是否喷出的直接依据。正常成功放电时电子密度峰值应该在10^17~10^19 m^-3量级放电区域集中在两个电极之间以及管口外围的等离子体羽流中。电场分布重点分析管口外电场是否足够强。射流的本质之一是流注头部自持传播——头部的高电场电离前方气体产生新的带电粒子维持传播。观察电场分布能看到管口外的电场增强区域是否沿着轴向延伸。电离速率这个量可以直观地看到放电在哪里发生、强度如何。对比不同电压下的电离速率能确定击穿阈值。电流波形积分电极上电流密度得到气体放电电流。波形里如果有代表微放电的窄脉冲说明模型把介质阻挡放电的脉冲特性抓住了。这个判断模型正确性的关键指标之一。后处理时把电子密度在轴向的分布画出来也就是取r0附近的值做一维曲线这能直接展示射流长度。射流长度可以用电子密度衰减到某个阈值比如1e16 m^-3处的轴向位置来量化。用“全局计算”功能提取每周期放电功率P 1/T ∫ V(t) * I(t) dt这是和实验功耗对比的常用参数。4. 常见问题与排查技巧实录4.1 模型一直不起辉或者起辉太弱这大概是被问到最多的问题之一。检查顺序是初始电子密度给太低放电需要很长的孕育时间才能从本底电离发展到击穿。把初始电子密度从1e10改成1e131e14 m^-3通常在第一个电压半周期内就能看到电子密度显著增长。这不是作弊放电气体的本底电离受宇宙射线和环境中放射性物质影响确实存在一个残余电子密度而实验里上一轮放电产生的残余电荷也会促进下一次击穿。电压幅值不够或者频率不对劲。8 kV峰值在13.56 kHz下通常能击穿1 atm氩气但如果你把频率降到50 Hz同样的电压幅值就不一定够。频率升高每个半周期内参与放电的有效离子次数减少击穿阈值会变化。次级电子发射系数设太低。介质表面吸附的电子被正离子轰击打出这些二次电子对维持常压DBD至关重要。很多默认模型用γ0.01我实际调下来氩气石英管体系γ取0.1更接近实验。这个参数不好直接测量但可以通过对比放电电流脉冲间隔来调。4.2 电子密度出现负值或严重振荡如果方程解出的电子密度在局部区域变负数说明数值解已经不稳定。负密度在物理上毫无意义但在数值上很容易出现因为漂移-扩散方程对流项在网格太粗或时间步长过大时会产生数值振荡。解决办法把最大时间步长缩小到5e-10 s级别在放电前沿处增加网格密度把等离子体模块的“对流稳定”选项打开。另外检查BDF方法阶数把阶数固定到1或2别让求解器自己切换因为高阶格式在刚性体系里很容易跳过物理极值点。4.3 电流波形完全没有窄脉冲特征有人跑完模型电流曲线就是一条正弦形状的位移电流曲线没有微放电叠加特征这说明放电过程没有真正发生或者介质表面电荷积累模拟得不充分。首先检查放电区最大电子密度是否达到10^17量级如果只有10^14说明只是弱电离放电没有完全发展。其次检查介质表面电荷积累边界条件——有时候忘选“边界条件类型电荷积累”导致介质表面永远不带电放电持续进行电流变成纯电容性电流。最后检查电压频率和数值时间步长是否匹配如果你每个周期只采了二十来个点那种纳秒级微放电脉冲是不可能被分辨出来的。4.4 计算速度慢到无法接受这个模型虽然只有二维轴对称但如果网格在放电区搞得太密比如统统0.02 mm那节点数也会到上百万非线性瞬态求解可能一整天跑不出结果。我自己试过一套加快策略组合先用轻度粗网格跑通模型流程确认没有物理设置错误再逐步加密并对比结果。 放电区用边界层各向异性网格而不是各向同性加密。 把等离子体模块里不太重要的中性激发态粒子暂时关掉只算电子和其他离子结构确认后再加回来。 如果流场和放电时间尺度差太多就用单向耦合不用双向耦合。 重离子输运方程的时间常数比电子长得多可以测试一下把重离子的扩散项关闭只保留迁移项有些情况下对结果影响很小但计算速度能翻倍。4.5 管口外射流形态和实验对比总对不上这个问题的原因往往很综合。射流在管口外的传播方向会受气流速度分布、电极极性、电场分布三个因素共同影响。如果实验里射流长度为15 mm模型算出只有5 mm优先排查电压幅值是否不够激发的活性粒子消耗太快气体中是否有杂质混入氮气或氧气会对氩气放电的电子能量分配产生显著影响模型里用纯氩气实验里如果用的是钢瓶氩气比如99.99%纯度就有微量的杂质管口气体速度是否足够氩气流速从1 m/s调到10 m/s射流长度变化明显仿真里也要扫。4.6 模型在某个特定时间点突然发散这个问题我在调双环电极间距时遇到过。几何参数变化导致放电模式突变从均匀DBD过渡到丝状DBD数值上表现为局部电子密度尖峰和极端的电场畸变方程组Jacobian矩阵变得奇异。这时候首先要确认“发散”是不是由网格奇点导致的。我在电极矩形拐角处加了一个0.1 mm的圆角把金属边角上的电场奇异去掉问题往往就解决了。模拟实验里真实的金属环边缘不可能是理想直角几何带圆角比直角更接近物理实际数值稳定性也会提高。如果圆角后仍在发散减小时间步长到1e-11 s试一次。等离子体模块在放电击穿瞬态电子密度在纳秒内可能变化6个数量级步长太大Jacobian更新根本跟不上。5. 仿真结果如何与实验验证对接很多人跑通模型就满足了但做工程应用的人肯定是想看模型对不对。我强烈建议在电脑仿真之前或者同步做小实验验证至少一个宏观指标比如放电电流波形或者射流长度。实验电流测量用电流探头罗氏线圈或皮尔逊探头夹在接地回路上得到的电流波形对比仿真电流波形。匹配度高就说明放电功率、脉冲幅值、相位关系都对。常压DBD电流脉冲一般在几十毫安到几百毫安量级如果仿真平均电流比实验高一个量级多半是边界条件设置太理想比如把电极到介质管之间接触电阻忽略不计了。射流长度对比用高速相机或ICCD拍射流形态量出射流从管口延伸的长度和仿真中电子密度降到背景值的距离对比。这个对网格密度、反应集完整度、流场准确度都极其敏感是一个很好的“真伪测试”。如果没有实验条件做验证那就至少做网格无关性分析和参数灵敏度分析在关键参数附近做一组粗、中、细网格结果对比确认预测值随网格加密趋于稳定然后再谈结论。这看起来是老生常谈但在等离子体仿真里真有不少人把网格相关误差当成物理结果在写文章的。6. 双环电极射流仿真还能怎么扩展双环电极模型的基础框架搭好后往各个方向扩展其实很顺这也是我当初看好这个模型的原因。常见的扩展方向包括加气体成分在氩气中混入少量空气或含氧成分研究放电产生的活性氧物种。需要在反应集里加入O2/N2电子的激发和离解反应并扩展电子碰撞截面表计算量会明显上升可控性还行。加气溶胶或液体液滴研究等离子体射流与液滴的相互作用需要开启相场或颗粒追踪接口并把等离子体化学和流体表面反应耦合起来复杂度和难度是另外一档。但双环电极能提供稳定的射流搭这种多相流模型是可行的。换电极结构比如从双环改成针-环、环-环间距扫描、改变介质管材质都能快速用现有模型做参数化扫描找到最优放电结构。比每次做实验试错快十倍不止。流场与放电双向耦合如果需求是分析射流加热或热对流影响那要将放电产生的焦耳热通过能量方程耦合给流体包含温升对气体密度、输运系数带来的反作用。这个会让模型计算成本增加一个量级我目前只在需要深入了解气压降和温度分布时才做。多个射流管并列仿真研究射流阵列之间的相互影响比如多个双环电极管并排放电需要用三维全模型这时候计算资源需求会非常大一般先用二维轴对称单管模型估算单管放电参数再建立简化的三维阵列模型。这套双环电极二维轴对称模型我现在已经跑通了每次调完几何或电压参数基本能在一轮仿真内看到合理的趋势变化。最后再分享一个小技巧建模时候千万别把电极和介质管壁的接触简化成完全导电现实里电极是缠绕在管子外壁的微观表面有粗糙度和微小气隙这些在快速高压放电过程会产生额外压降。如果实验电流波形和仿真对不上先往电极接触处加一点接触电阻往往能改善不少。做仿真这行要紧的从来不是把事情算得越复杂越好而是用最小的模型代价把事物的核心物理抓住。