这几天被沿面放电仿真搞得有点上头。起因是一个表面处理项目需要预测干燥空气下放电沿介质表面的发展形态我就在COMSOL里搭了一个介质阻挡放电DBD模型。实话讲这活儿比想象中要酸爽得多——核心就是电场、化学反应、边界条件三者的强耦合电场决定电子往哪跑、跑多快化学反应决定电子能量怎么被消耗、产生哪些离子和激发态粒子介质表面上的电荷积累又会反过来把电场扭曲到面目全非。每个环节单独拆开看都不算复杂但一旦合在一起就是一台高度非线性、强刚性、极易发散的精密机器。这篇文章把整个建模过程完整梳理一遍包括物理模型怎么简化、几何和边界条件怎么设、反应列表怎么选、求解器怎么调。主要写给正在入坑等离子体流体仿真、尤其是想把沿面DBD算出个门道的人。如果你只是想把电子密度云图跑得漂亮一点也能从后边的排查部分捞到几个关键设置。1. 物理模型拆解电场、化学反应、边界条件到底在吵什么1.1 沿面放电不是普通的体相DBD很多人第一次接触DBD脑子里浮现的是两块平行板电极中间夹一层气隙、介质板隔开电极的经典结构。那种叫体相DBD放电通道基本沿电场方向贯穿整个气隙表面电荷虽然存在但更多是作为“阻挡层”的一部分对放电形态的影响是全局性的。我这里要搭的是沿面DBD结构完全不同高压电极紧贴介质板上表面接地电极在介质板另一侧放电气隙实际是介质板上方靠近表面的那层气体。放电从高压电极边缘起始沿着介质表面向外围发展而不是垂直穿透气隙。这种形态下表面电荷的作用被放大了好几个量级——它直接影响放电如何一步步“爬”过介质表面也决定了每一轮放电何时熄灭、何时再触发。这个区别不是概念层面的咬文嚼字。在COMSOL里体相DBD模型往往可以用较粗的网格和较少的表面处理凑合过去沿面DBD不行。你必须在气体与介质交界面上把电荷积累、二次电子发射、电场不连续三个事情同时处理对漏掉任何一个放电形态都会跑偏。1.2 三者是同一件事的三个侧面为什么说“相爱相杀”我们可以把整个过程看成一个循环加压后高压电极边缘局部电场强度极高远超过空气击穿阈值背景种子电子在这个强场区被加速获得足够能量后碰撞电离中性分子产生电子雪崩电子、离子在电场中分离到达介质表面后沉积下来形成表面电荷表面电荷产生附加电场削弱电极方向的场强同时沿切向改变介质表面的电场分布电场畸变导致局部区域场强反而增强放电可能熄灭也可能沿表面继续推进。这个循环里电场是驱动者化学反应是燃料和刹车表面电荷是记分牌三者互相制约、互相依赖。你在仿真里不能把它们拆开处理必须同时在控制方程里求解。这也是为什么沿面DBD仿真比纯粹的电场求解或者单纯的化学动力学计算都要难一个数量级。1.3 建模前一定要做的简化假设开始建模之前先想清楚哪些地方必须简化和舍弃。我的做法是四个维度切一刀粒子模拟直接放弃。大气压DBD中电子自由程极短、碰撞频率极高如果上PIC粒子模拟一个毫米尺度的区域跑一个脉冲就需要天文数字的粒子数和时间步。流体模型漂移-扩散近似把电子按密度和平均能量处理虽然牺牲了部分微观细节但对宏观放电过程已经足够是公认性价比最高的路线。忽略流体动力学。放电过程发生在纳秒到微秒尺度气体分子根本来不及产生宏观流动。温度变化也不显著能量过程主要由电场输入和碰撞损失覆盖。除非你要研究气流辅助放电或者热耗散否则不要一开始就把流场耦合进去。气体成分简化。真实空气是复杂混合物但一个初始模型中只需要氮氧四比一水分、二氧化碳、稀有气体全部砍掉。这不是偷懒是为了把化学刚度控制住。后面想研究臭氧产量或者NOx生成再往反应列表里补。介质表面简化。电子显微镜下介质表面是坑坑洼洼的还有吸附水、灰尘。初始模型把这些因素全部忽略假设表面干净平整处理起来干净利落。这些简化做完了模型的主要物理特征仍然保留强电场驱动、电离附着竞争、表面电荷反馈。足够回答工程上的核心问题。2. 几何搭建与边界条件配置把“相爱相杀”落实成设置项2.1 几何结构、网格与尺寸表沿面DBD的几何在三维里是旋转对称的所以直接用2D轴对称模型能把自由度数量压到很低普通笔记本就能跑起来。具体尺寸我给一个参考结构件参数数值计算域半径×高度15 mm × 10 mm高压薄膜电极半径×厚度5 mm × 0.2 mm紧贴介质上表面电极边缘倒角半径0.05 mm介质层厚度1 mmAl2O3相对介电常数9.8接地电极覆盖区域介质层下表面全部气隙区域高度介质层以上10 mm电极实体画出来但在等离子体物理场中把电极域排除只作为静电端终端边界处理。电极边缘的倒角很关键不能省否则尖角处场强趋近无穷场分布完全失真放电起始位置也会被错误地钉在尖角上。网格划分策略上核心思想是“沿面方向加密远处粗放”。电极边缘附近场梯度极大网格最大尺寸控在0.05 mm以下放电沿介质表面向外传播的一段区域r5~15 mm网格控制在0.1 mm左右远处气体区域可以放宽到2 mm。介质层内部建议用扫掠网格或者边界层至少5层保证电位移场过渡平滑。提示沿面放电的空间尺度关键就在介质表面附近几十微米的薄层之内。网格太粗会把表面电荷的局部效应平均掉放电形态会变成一团模糊的高密度区而不是一支清晰的放电通道。2.2 电压激励与初始电场验证电压激励直接用50 Hz工频正弦幅值5 kVV(t)5000·sin(2π·50·t)。真正开始跑等离子体化学之前先单独用静电模块跑一遍稳态电场这一步几乎零成本但能发现大量几何和边界问题。检查两个事一是高压电极边缘的峰值场强。空气在大气压下的击穿场强约在30 kV/cm量级电极边缘有场增强效应峰值应该能到每米数百万伏特也就是几十kV/cm以上。如果峰值只有几千V/cm说明几何有问题或者电压太低放电物理根本起不来。二是沿介质表面方向的电场分布。沿面放电的特点是切向电场分量沿表面向外延伸这个分量的衰减快慢直接决定了放电能跑多远。用一条剖面线从电极边缘拖到计算域边界看电场怎么衰减心里就有数了。2.3 表面电荷积累和二次电子发射重头戏在介质表面这个界面上。绝缘介质与气体交界处电位移矢量的法向分量有一个突跳ε_gas E_n,gas − ε_diel E_n,diel ρ_s这里ρ_s就是表面电荷密度。静电场方程里这个不连续自动成立但如果等离子体物理场里没有把电荷积累方程挂上去ρ_s永远是零沿面放电的物理就缺了半边。表面电荷随时间的变化由电子和离子到达表面的通量决定dρ_s/dt e(Γ_i·n − Γ_e·n)在COMSOL的等离子体接口中这个操作是在介质表面的边界条件节点里勾选表面电荷积累选项把电子通量和离子通量的法向分量作为源项。注意法向方向定义为从介质指向气体符号搞反会导致电场自我增强模型立即发散。二次电子发射也要在这个边界上设置。γ系数取0.05发射电子初始能量给1~2 eV。二次电子发射的意义在于放电脉冲期间正离子轰击介质表面会产生新的电子进入气隙这些电子在强场区迅速加速成为下一轮雪崩的种子。如果没有二次电子发射DBD单个放电脉冲之后就会因为电子附着和复合而彻底熄灭无法形成周期性的微放电序列。注意我踩过最深的坑是表面电荷方程的符号问题。由于COMSOL的边界法向在不同几何中会随构建方向变化建议建完几何后先检查一下法向方向再确认电荷积累项的符号否则一跑就发散的模型八成是这个问题。3. 化学反应与输运参数反应列表才是模型的灵魂3.1 选择“漂移-扩散”等离子体接口COMSOL等离子体模块里空间放电仿真应该用漂移-扩散接口而不是零维的全局反应堆模型。这个接口求解的方程组包括电子密度连续性方程、电子平均能量方程、重物质输运方程、静电场的泊松方程。电子通量的形式是典型的漂移-扩散分解Γ_e −μ_e n_e E − D_e∇n_e电子迁移率μ_e和扩散系数D_e可以通过能量分布函数计算。COMSOL内部自带一组常用气体的截面数据和输运数据空气的电子输运属性可以直接从材料库调用不需要自己查表。如果你追求更准确的E/N依赖关系也可以用BOLSIG预跑一张迁移率随归约场强变化的表再导入COMSOL做插值但这不是新手阶段该做的事默认数据足够跑通。电子能量方程里的碰撞能量损失项会自动根据每个电子碰撞反应的能量阈值和反应速率系数计算不需要人手动去估计每个反应消耗多少能量。这也是用COMSOL做这件事比手搓方程省心太多的地方。3.2 最小可用空气反应组反应列表是整个模型的核心资产。这里必须克制——反应越多方程越刚收敛越难。一个能跑出沿面放电基本形态的模型我把反应压缩到12个编号类型反应式数据来源/说明1弹性碰撞e N2 → e N2动量交换截面2激发e N2 → e N2(A3Σu)氮气电子态激发3激发e N2 → e N2(B3Πg)放电发光的主要通道4电离e N2 → 2e N2氮气主电离通道5弹性碰撞e O2 → e O2动量交换截面6激发e O2 → e O2(a1Δg)氧气单线态激发7附着e O2 → O O−解离附着8电离e O2 → 2e O2氧气主电离通道9三体附着e O2 M → O2− M电子附着M为N2/O210离子转化N2 O2 → O2 N2重物质反应常量速率11电子离子复合O2 e → O2体复合12激发态淬灭N2(A3Σu) O2 → N2 O2激发态能量释放电子碰撞反应1~9号在COMSOL里通过碰撞截面数据自动计算速率系数也就是从截面做能量积分重物质反应10~12号直接输入速率常数。这个最小反应组已经能够覆盖放电启动、电子雪崩、附着熄灭、脉冲后粒子回收的全过程。为什么砍掉了几十个更精细的反应因为沿面放电的宏观形态主要靠电离与附着之间的竞争来平衡这个平衡抓住了模型的骨架就立住了。氮氧化物、臭氧的高精度量化是后续工程问题不是模型跑通的前提。每增加一个电子碰撞反应通道化学刚性的复杂度就涨一截收敛难度成倍增加。我一直坚持一个原则先最小集跑通再逐个塞进新反应每次只改一个变量。3.3 初始值与种子电子的作用放电不可能凭空启动。真实环境里宇宙射线和辐射本底不断产生稀疏的电子离子对仿真里用一个均匀的背景电子密度和离子密度来模拟这个效果。初始电子密度和正离子密度都设为1×10^12 m⁻³负离子可以忽略。这个数量级能保证在强场区快速出现雪崩但又不至于人为地提前触发放电。把初始密度调到1×10^10可能一个电压周期过去了放电还没起调到1×10^14放电会提前发生脉冲频率也偏快。所以初始值不是一个可以随意拍脑袋的细节——它直接影响放电的相位和脉冲形态。提示如果第一步只是想验证边界条件和化学模型没有明显错误可以临时把初始电子密度调到1×10^14让放电快速触发跑通流程后再降回真实量级。不要把这个临时值当成最终设置否则放电时序会严重失真。4. 求解策略与收敛调教实录4.1 求解器设置BDF、容差、时间步沿面DBD模型的刚性问题主要来自工频慢变化和放电纳秒快变化的对比。一个5 kV、50 Hz的正弦电压半周期10 ms但放电脉冲宽度只有几十纳秒。如果要求解器用均匀步长把整个过程都解析出来那得跑多少个时间步COMSOL里选用BDF向后差分求解器允许自适应步长。关键设置有三个初始时间步建议1×10⁻¹² s。不要指望求解器自动给一个靠谱的初值放电触发瞬间的突变如果发生在第一步内直接就看不到。给一个足够小的初始步长让求解器先“站稳”再自行放大步长。最大时间步放电事件发生的窗口附近建议限制在1×10⁻⁸~1×10⁻⁹ s远离放电窗口可以放宽。如果整个时间域都限制在1 ns步长计算量会让电脑直接罢工。一个务实做法是分两段跑先跑0~10 μs的短窗口看单个脉冲再跑长窗口看周期序列。容差相对容差用1×10⁻³就够绝对容差按变量量级来。电子密度从1×10¹²跳到1×10¹⁹绝对容差给到1×10¹⁰~1×10¹¹比较稳妥电势的绝对容差给1~10 V。绝对容差设得太小求解器会在每个小变化上死磕步长缩到飞起几小时算不完一个周期。我的操作顺序是先用固定电压比如5 kV直流跑几十纳秒确认放电形态、电子密度峰值、表面电荷增长都是合理的再切换到工频正弦激励观察不同相位下的放电行为。这把“放电物理”和“电压波形调度”两类问题隔离开排查起来轻松太多。4.2 典型发散问题与排查速查问题1初始几步残差就爆了。大概率是初始时间步太大或者初始电子密度渲染出局部过大的电场。先减小初始时间步到1×10⁻¹² s把绝对容差调大一个数量级。如果还爆把二次电子发射临时关掉看是不是表面电荷反馈太强。问题2表面电荷振荡、放电形态呈锯齿状。网格太粗或者表面电荷方程符号反了。沿介质表面的网格不要超过0.2 mm重点区域0.1 mm。符号问题按前面说的检查边界法向方向。问题3放电迟迟不触发。先查高压电极边缘峰值场强是否达到几十kV/cm再查初始电子密度是不是设低了最后确认介质表面有没有误设成接地或者高电位边界。问题4计算量爆炸一天跑不完一个周期。反应列表太多、时间窗太长、网格太密。把时间窗先从μs级开始反应组砍到最小集远离放电区域的网格放宽到毫米量级。跑出一个能完算的模型比跑出一个永远算不完的精致模型有价值得多。问题5电子密度云图呈弥散状没有清晰的放电通道。漂移-扩散近似下这是正常现象。流体模型天然会抹掉部分丝状放电的精细结构大气压空气中的流注也因此只能看到大致的亮区轮廓。别急着调整模型先看电流脉冲波形和时序是否与实验常识对得上。4.3 结果提取与合理性校验沿面DBD模型跑完重点看四个输出量电子密度云图用对数色标看峰值应该在1×10¹⁸~1×10²¹ m⁻³量级。如果峰值只有1×10¹⁵说明电离反应没有真正起主导如果超过1×10²²说明网格或初始条件有问题。介质表面电荷密度沿半径方向做一条线图观察电荷积累随位置和时间的变化。量级在μC/m²量级是合理的。不同位置电荷极性的空间分布能直接反映放电沿表面推进的形态。电场分布从电极边缘沿介质表面拉剖面线看电场切向分量如何衰减可以判断表面电荷是否有效屏蔽了前端电场。放电电流对电极边界做积分得到电流波形。工频激励下通常能看到每半周期一组脉冲簇脉冲宽度几十纳秒到几百纳秒电流幅值在毫安到几十毫安量级。这四个量对不上任何一个回去查模型问题比继续调参更重要。最后说点不中听的实话。沿面DBD仿真几乎每个人都栽在同一个坑一上来就想要完整空气化学模型立马跑出多周期微放电结果就是模型越来越复杂、残差越来越难看、机器越来越烫。我给的建议是先把反应砍到最小集先用固定电压激发单个脉冲再慢慢加物理对象。一个能收敛的粗模型远胜一个堆满网格和反应却永远跑不动的精致模型。我自己每次做新结构都坚持先几何后静电、先短窗后长窗、先少反应后多反应的顺序每一步都盯住输出量做对比确认不跑偏再往下走。这套流程看起来琐碎实际上才是最短路径。