你如果要搞电弧模拟第一反应多半是“这不就是CFD里加个热源嘛”。真上手了才知道电弧这玩意儿是流体、传热、电磁场、相变甚至化学反应的混合体典型的多物理场强耦合问题。Fluent在这类场景下的尴尬之处在于它自带的功能跑个普通流动收敛没问题但你要把焦耳热、洛伦兹力、辐射散热这些源项塞进控制方程就必须跟UDF“深度绑定了”。我自己第一次用Fluent做电弧仿真是在一个钨极氩弧焊的简化模型上。当时天真地以为只要在能量方程里加个热源就行结果算出来的温度场完全失真电弧形态也根本不是那个“钟罩形”。后来才明白电弧内部电流产生的焦耳热是分布式的电流密度本身又由电场控制而电场又受电导率温度的函数调制——这三者必须同步求解缺一不可。这就是UDF的用武之地通过用户自定义标量UDS把电势方程“塞”进Fluent再用源项把电磁力、焦耳热耦合回动量方程和能量方程才算把电焊工在虚拟世界里真正“造”出来。这篇文章我打算把整个思路从物理模型建立、UDF代码骨架到调参、收敛、排查完整过一遍。适合哪些人看刚接触电弧仿真、被多物理场耦合搞得头皮发麻的研究生或工程师还有那些已经把Fluent跑熟但还没碰过UDF、想朝多物理场更进一步的人。我不会给你堆公式推导重点讲“Fluent里怎么实现”和“我踩过的那些坑”。1. 方案选型为什么不能只用Fluent自带模块硬算1.1 自带MHD模块的局限性Fluent的MHD模块在很多型号里是默认关闭的你需要通过Scheme命令手动激活。就算激活了它适用的场景也极其有限三维电流、感应磁场计算还凑合但电弧模型里的边界条件、源项往往高度非线性MHD模块没法让你自由定义“电导率随温度突变”这种极端行为。更麻烦的是MHD模块的电流方程求解基于磁场方程矢量势法而大部分电弧模拟更适合直接用标量电势法。标量电势的好处是方程少、边界条件直观——阳极、阴极可以直接给定电流或电压边界。矢量势在三维里有分量、有规范选择问题对一般工程人员来说调试成本太高。所以我的结论很直接用UDS解电势方程再用UDF写源项和物性是目前精确可控、自由度最高的方案。1.2 UDS方程组的设计思路电弧的电场控制方程本质是电流连续性方程∇·(σ∇φ) 0σ是电导率φ是电势。这个方程在Fluent里可以用一个UDS来解标量命名为“potential”扩散系数设为σ对流项为零无瞬态项。只要把UDS的扩散系数设置成和电导率挂钩求解器就能在当前温度场下解出电势分布。这里有个关键设计电导率是温度的函数而温度场会改变电导率电导率又反过来决定焦耳热的大小焦耳热再改变温度场。所以整个求解过程必须外迭代循环先固定温度求电势算焦耳热和洛伦兹力再更新温度和流场反正中间的物性也要同步刷新。这个耦合闭环必须处理好否则会出现典型的数值失稳。1.3 为什么推荐“温度上限截断”的物性建模电弧温度动辄上万开尔文空气或保护气体的物性数据在这个区间内是不能查常温表的。常见做法是直接从文献里引数据写插值函数比如空气在300K到30000K的电导率、粘度、比热、密度。这些数据点很多论文里都有拿过来做线性插值即可。这里很容易犯的错是电导率在低温段比如5000K以下近似为零若直接用插值源项里可能会除零或产生极端非线性。实践上我习惯给电导率加一个下限比如1e-2或者更小防止数值爆炸同时如果温度超过数据库上限直接按上限值截断。别觉得这“不够物理”工程仿真里这种处理很常见否则收敛性没法保证。2. 核心细节电势UDS、源项UDF和物性控制三板斧2.1 电势方程UDS的Fluent配置在Fluent里添加一个UDS需要进入User-Defined Scalars面板。你需要把UDS-1的名字改成“potential”对它的扩散项通过UDF指定。具体操作是在UDS的Diffusivity那里选user-defined然后关联到你的DEFINE_DIFFUSIVITY函数。这个扩散系数的函数很简单返回一个实数值即当前单元温度对应的电导率。注意Fluent的UDS默认是无量纲的而你在C代码里写的所有系数都必须是国际单位制。电导率的单位是西门子每米S/m你需要在物性插值函数里确保返回的值是这个单位。边界条件设置上如果是电流控制型电弧比如钨极氩弧焊通常恒流控制阴极表面设UDS边界为固定电流密度实际上Fluent里UDS的边界条件里没有“电流”这个概念你得通过“通量”边界或源项处理。常用的技巧是给阴极边界设一个弱约束电势值或者干脆在近阴极单元加一个体积源项等效注入电流。这个方法操作起来有点绕但保证总电流大致恒定。阳极边界一般设零电势作为参考地。如果你用电压驱动型比如电弧炉的相电压那就更简单直接阴极、阳极各给一个固定电压。不过那种情况下电流分布要靠收敛自然确定初始迭代阶段容易震荡建议用较小的欠松弛。2.2 能量源项焦耳热的正确写法焦耳热功率密度是q σ |∇φ|²在UDF里你需要取当前网格单元梯度。Fluent提供了C_UDSI_G(cell, thread, index)宏可以直接拿到UDS梯度。这个梯度向量是三维的计算模的平方即可。但有个坑必须提醒C_UDSI_G取出的是迭代过程中存储的梯度如果你用的是基于格林-高斯的梯度重构某些单元交界处会出现梯度振荡。这会导致焦耳热源项在温度梯度和电势梯度都很陡的电弧弧柱区出现非物理尖峰。解决的办法是对源项做限制clipping比如超过某一最大功率密度的值强行改成该最大值。具体阈值需要根据你的网格密度和电流大小估计我一般取理论最大值的1.2倍左右。这个限幅操作看虽然不“严格”但实际计算稳定性和收敛速度都会好很多。你可以在代码里用Message宏把超限的次数打出来如果次数很少几万分之一记说明限幅影响可忽略如果频繁被触发说明你的网格和电流设置可能不合理需要回去改。2.3 动量源项洛伦兹力的三维分量处理电流密度是矢量J -σ∇φ负号源自电场与电势梯度的关系磁场B由电流分布感应产生。完整的三维磁场需要通过另一组UDS磁矢量势或毕奥-萨伐尔积分实现。但对轴对称的电弧模型洛伦兹力的径向分量有一个经典简化公式磁场只有环向分量Bθ大小可以通过安培环路定理从电流分布算出。简化成公式就是Jz 轴向电流密度r 半径该半径内通过的总电流I_enc则 Bθ μ0 I_enc / (2πr)洛伦兹力fr -Jz × Bθ方向挤压电弧fz Jr × Bθ。在轴对称模型里Fluent的坐标轴中心线是x轴或z轴你需要根据你的模型方向写分量。如果你做的是三维模型建议还是稳妥地用磁场UDS方法否则二维简化在三维里的误差不好估计。我在C代码里一般用C_UD_MI(cell, thread)存储每个单元的中心坐标然后基于全局电流积分算I_enc。这个“积分”操作需要遍历网格单元所以UDF里要写一个循环所有单元的宏用NV_V等宏逐点累加。这是很费时的部分我建议每时间步或每迭代步都更新一次即可不要多线程同时写共享数组否则会遇到并行计算的race condition问题。2.4 物性插值电导率、比热、黏度的温度依赖物性这一块Fluent自带的材料库基本帮不上忙。你需要自己把温度-物性数据表做成二维数组在UDF里用线性插值函数返回对应值。比如电导率的典型输入是温度开尔文输出S/m。数据来源可以用公开发表的空气热物理性质表或者从电弧仿真经典论文里提取曲线。写插值函数有几个细节一是温度区间要覆盖你仿真可能出现的所有范围比如300K到25000K超出上限就返回上限点的值低于下限就返回下限点的值别让插值函数自己外推二是要注意数据点的单调性电导率在某个温度区间可能非单调先升后降又升如何处理不影响你焦耳热源项的正确性但会影响收敛。比热和粘度、热导率同理直接做插值。密度的话如果不考虑电离引起的组分变化理想气体状态方程即可。对于保护气氛氩气、氦气物性注意换成对应气体的数据空气的数据不适用。3. 实操过程与核心环节实现3.1 几何建模与网格划分的注意事项电弧模拟的几何往往是轴对称的简化成二维可以大幅降低计算量。但二维轴对称和三维的差别一定要心里有数二维里没法考虑电弧摆动、偏吹等实际现象如果你关注的是电弧炉或开关电弧等强三维效应场景建议直接上三维网格。网格划分的重点区域是弧柱区。弧柱的温度梯度极大从中心上万K到边界几千K电导率变化剧烈这里的网格尺寸不能太大。我一般用渐近网格弧柱中心区域加密到0.1mm以下对于氩气电弧往外逐渐变疏。网格质量上是老生常谈的偏斜率尽量低于0.8负体积绝对不能有。电极区域的处理也可以在这步决定。钨极和铜阳极是固体域但Fluent默认是流固共轭传热CHT吗如果选共轭传热固体域同样要网格划分好且必须设置好内部界面interface。这里很容易出bug界面两边的网格节点不对齐会导致插值误差和收敛变慢。建议用节点对齐conformal网格来处理固液交界面。有个小技巧把阴极表面附近几个单元的网格设成楔形或边界层网格这样能更精确捕捉近极区的压降和高温梯度。电弧的阳极压降和阴极压降虽然实际物理机制复杂但在简化的CFD模型里常常被建模成很薄的近壁高温区如果网格太粗这个薄层的温度梯度和电势梯度都捕不到可能导致整体电流偏小。3.2 UDF代码框架从初始化到迭代循环我习惯把UDF按功能拆成几个独立的源文件再在Fluent里用编译型UDFcompiled而不是解释型UDF。解释型UDF慢且对C语言的限制多做电弧这种复杂问题不合适。代码的骨架大致是DEFINE_PROPERTY(electrical_conductivity, cell, thread) { real temp C_T(cell, thread); return interpolate_sigma(temp); } DEFINE_DIFFUSIVITY(potential_diffusivity, cell, thread, i) { return electrical_conductivity(cell, thread); } DEFINE_SOURCE(energy_source, cell, thread, dS, eqn) { real sigma electrical_conductivity(cell, thread); real grad[3]; C_UDSI_G(cell, thread, 0, grad); real q sigma * (grad[0]*grad[0] grad[1]*grad[1] grad[2]*grad[2]); dS[eqn] 2.0 * sigma * (grad[0]*...); // 对温度的导数近似留空或简化 return q; }注意源项的dS[eqn]是雅可比项写得好能改善收敛性。焦耳热对温度的导数严格说需要电导率对温度的函数求导这一步比较复杂。我自己的经验是初期先把这个导数设为零让源项显式处理对稳态问题残差收敛慢但能接受。如果实在不收敛可以试试用数值差分近似dS[eqn]。迭代运行过程中最好每隔一定步数打印一下总电流、最高温度、弧柱半径等关键量。在UDF里用Message宏输出即可。这样你能在收敛之前就判断结果是不是靠谱的避免算了几天发现温度根本不在合理范围内。3.3 磁场处理的两种实现方案对比磁场这一块我试过两条路线一是额外的UDS求磁矢量势三维,二是轴对称简化公式。对比下来轴对称简化方案实现简单、代码量少、计算速度快但前提是你必须确认你的问题可以降维。如果你在后处理想要展示磁场分布云图简化方案也能出图因为从Bθ公式可以直接计算每个网格的磁场。三维的UDS方案需要多解三个标量方程Ax, Ay, Az。磁矢量势方程是∇²A μ0 J矢量势UDS的扩散系数设为恒定值1/μ0源项为电流密度分量。这个方法求解稳定性还可以但金属域里的电流被固体区域截止需要特别处理固液界面处的电流连续性代码复杂度高不少。我建议你做第一版先跑二维简化把耦合逻辑跑通了再扩展三维。别一上来就啃最难的那种。做仿真最忌贪大求全先有个能出结果的东西再说。3.4 求解器设置与松弛因子的实战选择Fluent求解器里压力-速度耦合算法用Coupled还是SIMPLE我的经验SIMPLE比较稳Coupled收敛快但在强源项情况下容易发散。电弧这种问题源项高度非线性我习惯用SIMPLE加一阶迎风起步。温度方程、电势方程、动量方程都要设欠松弛因子。能量方程一般可以0.9以上电势UDS的欠松弛建议从0.5起步。如果你发现残差锯齿形波动剧烈第一反应不是动时间和步长而是把电势UDS的松弛因子降到0.2左右看看。时间步长如果是瞬态或伪瞬态步长也很有讲究。电弧稳态模拟可以用伪瞬态pseudo transient方法Fluent有开关键。伪瞬态步长设置个1e-5秒左右慢慢推进比定常求解稳健得多。这个技巧很多人不知道强烈推荐试试。4. 常见问题与排查技巧实录4.1 电势场不更新或始终为零这个现象通常分两类一是UDS方程本身没被激活求解检查一下Solution Methods里有没有勾选UDS方程求解二是扩散系数设置错误。一个很隐蔽的坑是如果你在DEFINE_DIFFUSIVITY里返回的是一个常数而不是随温度变化的插值结果那么温度场的反馈循环就被切断了。排查方法很简单在扩散系数UDF里加一条Message打印当前温度范围跑几步看输出。如果打印出的温度始终是初始值说明UDS压根没有参与迭代。4.2 电弧温度虚高或飞温飞温是电弧模拟中最常见的问题。虚高的原因一般是焦耳热源项没有限幅或者电导率对温度的插值在高温端太“平”。什么意思呢当温度升高到一定程度后电导率增长缓慢但电源功率电流固定还在持续输入能量辐射散热跟不上温度就一路飙上去。解决方法有三步排查检查能量方程里是否包含了辐射损失项——电弧模拟不能忽略辐射尤其大电流和高压情况下辐射散热占比可能达到30%以上。如果你设的是恒流边界条件还要检查出口边界温度是否允许过冲。最后就是限幅源项了直接把异常尖峰砍掉。4.3 残差震荡不收敛残差收敛曲线出现持续的正弦波震荡且不随时间衰减几乎可以断定是电磁场和流场耦合过强的问题。应对策略按优先级排序降低电势UDS欠松弛因子改用一阶迎风格式调小伪瞬态时间步长。另一个有意思的排查方向是动网格问题。如果你用的是动网格模型模拟钨极进给网格更新时电弧形态会跟着变化旧时间步的场数据在网格变形后会剧烈振荡。动网格设置里建议把remeshing的尺寸函数范围控制好避免网格畸变导致UDS梯度错误。Fluent里CPU绑定核数对UDF的并行计算稳定也有影响我试过单机用8核比用16核更稳具体原因没深究大概率是共享内存的race condition问题。如果你发现并行计算结果和串行不一致第一件事就是把UDF里的全局数组和公共变量检查一遍。4.4 一个容易被忽略的“无量纲”陷阱Fluent的UDS方程在求解时后台会对变量做无量纲化。如果你在UDF里做的计算涉及物理单位务必搞清楚Fluent的参考值Reference Values设置。默认的参考长度、参考密度等可能会让UDS梯度输出一个异常大的数值再乘上电导率焦耳热源项可能直接溢出。这个坑我当年查了很久。Gradient取出来在数值上是无量纲的需要转换回有物理量纲的值才能用于源项。麻烦的是这个转换因子Fluent没有单独提供给UDF你需要从网格尺寸和标量参考值自己推导。例如电势标量的单位是伏特如果参考长度设为1米梯度单位就是V/m但如果参考长度是毫米梯度值可能放大了1000倍。最省事的规避办法把所有UDS相关的物理量都手动控制在SI单位然后Fluent里把参考值全部改成和网格几何一致比如尺度单位用米。运行几个迭代步后输出电流密度检查一下是否在合理范围。5. 一点经验总结非套话电弧模拟这活儿本质上是一半物理、一半数值技巧。物理模型选对了剩下的全在和收敛性较劲。我做了几个月的电弧仿真最大的感悟是别迷信高端数值格式一阶迎风格式配合合理的欠松弛因子往往比什么QUICK、高阶格式都够用。还包括那些看起来很聪明的限幅、截断处理——这些对不上文献的时候你可能觉得心虚但算得上算得稳比“严格”但发散的结果有意义得多。如果你刚开始做我建议先跑一个非常简单的空气电弧模型把UDS电势方程和焦耳热耦合拎顺了再加母线力、辐射散热和动网格。每加一个物理场至少留一周的调参时间。这个过程很磨人但等你把第一个电弧形态算出来了那种成就感真的很难替代。最后再分享一个小技巧Fluent的收敛判定不要全看残差电弧模拟中全局电流、最高温度、弧柱半径这些物理量是否稳定才是更重要的指标。有的算例残差到1e-4甚至到1e-3就基本稳定了过分追求残差反而浪费时间。把后处理里的监测曲线加好比盯着迭代日志有用得多。
