去年做换流变阀侧绝缘可靠性评估的时候被一个问题反复纠缠如果绕组里掉出几颗铜屑或者绝缘纸板边缘剥落了几根细小纤维这些颗粒在绝缘油里到底怎么跑油泵启停、负载周期性波动带来的油流振荡正好对应了往复运动这个工况。为了把这个问题彻底摸清我用COMSOL搭了一个简化模型把铜颗粒和纤维颗粒放进绝缘油流场里对比它们在往复流动下的轨迹、相位跟随性和沉积规律。这篇文章把从建模取舍、受力排序、边界条件设置到后处理读图的完整过程记录下来给同样在研究绝缘油颗粒污染物运移的朋友当一份可复现的参考。1. 换流变和油浸设备里为什么要死盯铜颗粒和纤维颗粒1.1 颗粒从哪来又为什么危险油浸式电力设备内部的绝缘系统从来不是绝对干净的。铜颗粒最常见的来源是绕组焊接毛刺脱落、导体在电动力或机械振动下的微动磨损以及电化学腐蚀产物。纤维颗粒则来自绝缘纸、纸板或压板的边缘磨损还有制造和干燥过程中残留的碎屑。这两种颗粒一旦进入油隙问题就不是脏不脏那么简单了。在强电场环境下导电的铜颗粒会畸变局部电场造成电场集中容易诱发局部放电甚至直接搭桥导致击穿。纤维颗粒虽然本身半绝缘但吸湿后电导率上升同样会成为放电起始点而且纤维容易在油流和电场作用下聚集形成小桥这就是所谓的小桥效应。国内外不少绝缘油色谱分析和设备故障归因里颗粒污染都被列为重要诱因。我在这类可靠性分析里最关心的不是有没有颗粒而是颗粒在特定运行工况下会去哪里、做什么运动。因为不同运动形式直接决定了颗粒是停留在油流中被循环过滤掉还是反复撞击绝缘件表面、卡在狭窄油道里、或者在电极附近来回穿梭。厂商给的油颗粒度测试只能测出浓度和粒径分布测不出运动形态所以必须靠仿真补位。1.2 往复运动不是拍脑袋想出来的工况很多人一看到往复会以为是人为设置的理想化条件其实这套工况在真实设备里到处都是。第一类来源是油泵与冷却系统的周期性调节。强油循环冷却的变压器油泵并不是一直恒定转速运行温度控制会触发启停或变频调节每次切换都会在油道里产生一次压力波油流速度和方向都会出现周期性变化。第二类来源是负荷波动与绕组振动。负载电流日波动、换流阀周期性触发都会让绕组和铁心产生振动这个振动通过结构传递到油流中形成受迫振荡。第三类来源是交变电场本身。极性反转或交流电压下半绝缘颗粒和带电金属颗粒会受到周期性的电动力也会形成往复迁移。那为什么不能只研究单向流动因为单向流是扫过去就完了颗粒要么被带走要么沉积在某处而往复流是颗粒在一个区域里来回穿梭有反复积累、反复碰撞、反复穿越危险区的可能。比如一个铜颗粒在单向流里可能一次就被冲到油箱底部但在往复流里可能每半个周期都被带回电极间隙附近这种涡流式反复扫掠造成的放电风险比单向流严重得多。研究往复运动形式本质上是研究颗粒在振荡流场中的响应特性。1.3 这个模拟到底想回答哪几个问题在做模型之前我先列了项目真正想回答的问题清单全仿真过程都围绕这些问题展开在给定的油流振荡幅度和频率下颗粒是跟随油流同步往返还是出现明显的相位滞后铜颗粒和纤维颗粒谁更容易沉积谁更容易被带离危险区域大颗粒和小颗粒的运动形式差异有多大能不能按粒径分选出不同行为颗粒会不会在某个周期里反复撞击电极或绝缘障碍物影响局放风险带着问题去看模型结果就不容易陷入做了一堆漂亮轨迹图但不知道说明什么的误区。2. 动手建模前先把受力主次、量级和耦合方式定下来2.1 把真实工况抽象成可计算的物理场景真实变压器里的油道几何非常复杂绕组、垫块、撑条到处都是没有必要一上来就建全尺寸模型。我的做法是先做一个能反映核心物理过程的简化几何平行板电极油隙电极间距取典型绝缘油道尺度比如5毫米或10毫米。流场方向沿板间长度方向入口在一侧出口在另一侧油流以正弦速度振荡。这个二维轴对称或二维平面模型足以对比颗粒在往复流中的运动规律。如果你要研究的是局部结构比如某条窄油道入口处的颗粒行为也可以换成三维短通道模型但计算量会增加一个量级。我的建议是先用二维把物理规律跑通再决定要不要升级到三维。COMSOL里几何改动比较方便网格和边界条件重设很快把二维阶段的参数和经验直接带过去就行。2.2 颗粒受力的主次排序不能一股脑全开COMSOL的粒子追踪模块里可以激活一大堆力曳力、重力/浮力、电场力、压力梯度力、虚拟质量力、Saffman升力、Brownian力、Basset历史力等等。工程上最忌讳的就是把能开的力全打开计算慢不说还容易数值不稳定。实际项目里要根据颗粒尺寸和工况判断哪些是主控力哪些可以忽略。我通常遵循这样的优先顺序曳力Stokes曳力几乎所有颗粒运动场景的主控项100%要开。重力与浮力对铜颗粒这种高密度差颗粒必须开对纤维颗粒密度只比油高一点点作用弱一些但也不能完全省因为往复流中的净沉降会影响沉积位置。电场力如果研究的是带电荷的金属颗粒或极性反转电场中的带电颗粒电场力是核心。但要注意中性颗粒在梯度电场里受到的介电泳力不等于简单的库仑力需要单独去算。虚拟质量力与压力梯度力在往复振荡流场里流体加速度不可忽略这两个力对颗粒跟随性的影响很明显建议保留。Basset历史力高频往复和高密度比时会有影响但数值上很难收敛如果频率只在几赫兹以内可以先不开用结果对比决定要不要补。Saffman升力与速度梯度有关在壁面附近剪切强的地方影响大如果关注近壁运动可以开。Brownian力对几十微米以上的颗粒几乎没有影响只有亚微米颗粒才需要考虑一般绝缘油颗粒研究可以不开。2.3 先用斯托克斯数判断颗粒响应再谈建模难度这里有一个快速估算颗粒跟随性的办法弛豫时间。对于球形颗粒Stokes弛豫时间可以写成τ_p ρ_p * d_p² / (18 * μ)其中ρ_p是颗粒密度d_p是粒径μ是油液动力粘度。这个时间表示颗粒在流体速度阶跃后追上流体所需要的时间尺度。典型变压器油40℃时动力粘度大约0.008 Pa·s油密度约890 kg/m³。我按这个估算了几组数颗粒类型粒径密度弛豫时间铜颗粒10 μm89606.2e-7 s铜颗粒100 μm89606.2e-5 s铜颗粒1 mm89606.2e-3 s纤维颗粒10 μm13509.4e-8 s纤维颗粒100 μm13509.4e-6 s再取流动特征时间往复运动下取1/f。比如频率1Hz流动特征时间约0.16秒。把弛豫时间除以流动特征时间得到斯托克斯数Stk。上面所有数值算下来Stk都在0.04以下说明粒径在100微米以内的颗粒对油流振荡的跟随性都非常好基本是油流到哪颗粒到哪。但不代表它们的行为完全一致因为铜颗粒密度摆在那里重力引起的沉降速度差异巨大这在后续结果里体现得很明显。如果你处理的工况频率很高比如几百赫兹的振动激励那Stk会迅速增大颗粒开始跟不上流场运动形式会完全不同。这时候Basset力、虚拟质量力的权重都会上来建模策略也需要调整。2.4 单向耦合还是双向耦合先看颗粒浓度还要决定颗粒与流场的相互作用。我们关心的是稀颗粒问题颗粒体积分数远小于0.1%颗粒对油流反作用可以忽略直接采用单向耦合先解流场再让颗粒在流场中运动。COMSOL里完全可以在一个瞬态研究中同时求解层流和粒子追踪粒子追踪读取流场结果但不会反过来改变流场。如果颗粒浓度很高比如油液已经严重劣化、颗粒聚集或者你在做滤芯过滤性能分析那就要考虑双向耦合让颗粒拖拽对流体动量方程产生额外源项。双向耦合的计算代价明显更高而且需要特别注意颗粒数量统计与源项光滑化不然压力场会出现剧烈的局部波动。我这个项目里浓度很低单向耦合完全够了。3. 在COMSOL里把往复运动模型搭起来3.1 物理场接口与多物理场耦合关系我使用的物理场接口组合如下层流接口spf描述绝缘油的流动选用瞬态研究。如果最大流速低、几何尺度小到雷诺数远小于1可以直接改蠕动流接口计算更稳定但往复流瞬时速度会变化不少时段雷诺数还是能到几十几百所以直接保留层流更省心。流体流动颗粒追踪接口fpt释放铜颗粒和纤维颗粒设置各受力项。如果考虑电场对带电颗粒的影响再加一个静电场接口es给电极施加载荷然后在粒子追踪中通过表达式施加电泳力。多物理场耦合的本质是把流场速度场和电场结果传给粒子追踪模块。在COMSOL 6.x里层流接口的默认变量是spf.U粒子追踪接口可以直接引用spf.U作为流体速度电力计算则通常要靠自定义表达式把电场强度空间分量传给颗粒。这类模型的物理关系可以用一句话概括流体运动通过曳力带动颗粒电场通过电荷作用力改变颗粒轨迹颗粒作为稀相不反过来改变流场和电场这就是整个耦合的闭环。3.2 几何、油品和颗粒参数的基准值我建议所有参数都用全局定义管理这样后面的参数化扫描会很方便。基准参数如下参数名值说明L50 mm油道长度H5 mm电极间隙rho_oil890 kg/m³绝缘油密度mu_oil0.008 Pa·s动力粘度U00.05 m/s往复速度幅值f1 Hz往复频率rho_cu8960 kg/m³铜颗粒密度rho_fiber1350 kg/m³纤维颗粒密度d_p10-100 μm颗粒直径N_p100-500每批释放数量初始释放位置我一般选在入口附近均匀分布或者直接在油域里随机均匀分布前者能看入口效应后者能更直观反映全局行为。颗粒数量不要一上来就放几千个一方面是后处理轨迹图会糊成一团另一方面是粒子追踪的计算时间与颗粒数量强相关先用100-200个颗粒摸清趋势再增加。3.3 三种往复驱动方式的实现与选择往复运动在COMSOL里有三种典型的驱动实现我三种都试过各自适用面不同。第一种是入口速度正弦振荡也是我推荐大多数人优先尝试的方案。在层流入口边界上直接写U0 * sin(2*pi*f*t)这种设置简单直接求解稳定后油域内的速度场就会呈往复振荡。要注意的是第一个计算步如果直接从0跳到U0压力场容易产生数值冲击我习惯乘一个短时间的斜坡过渡后面会在踩坑部分细说。第二种是移动网格模拟机械振动。比如你想模拟电极或绝缘件在振动信号作用下的微位移给某个边界加周期性位移条件。这需要增加变形几何接口让网格随着边界运动而变形。往复频率、振幅都直接作用在几何上颗粒在移动边界和变形油域中运动。这种方案最贴近机械振动场景但动网格对网格质量和时间步要求很高边界振幅稍大网格就可能在拐角处翻转。第三种是交变电场驱动。静电场接口中给电极施加随时间正弦变化的电位颗粒如果有净电荷就会受到周期性的电场力从而往复运动。这种方式适合研究带电颗粒在电场作用下的往复迁移但要注意真实油浸设备里颗粒的带电状态并不稳定电荷量是仿真里最难标定的参数通常要做参数扫描。我从实际操作角度说如果你关心的是油流振荡对颗粒的影响用第一种就够了加了电场耦合后再叠加交变电场组合成既有油流振荡又有电应力的双重往复能看的东西更多但也要做好参数敏感性分析。3.4 颗粒边界条件和时间步设置的原则粒子追踪的壁边界条件决定了颗粒碰到电极表面之后怎么处理。COMSOL里提供了冻结、反弹、消失、吸附等选项。我一般把电极表面设成冻结或吸附因为一旦颗粒到达电极表面就代表它脱离了油流主体后续运动不再是这个模型关心的重点让它停下来反而方便统计沉积量。如果是研究颗粒撞击频率那就设成反弹同时设定恢复系数钢-油-铜这类刚性碰撞的恢复系数我一般取0.8左右。另一个必须设置的是粒子最大时间步。COMSOL的粒子追踪求解器允许多时间步推进即全局流场求解器走一个大步粒子可以在每个大步内再细分更多子步。这个机制处理快速受力变化时非常关键能有效避免颗粒在单个时间步里穿越未解析的薄网格层。我通常把最大粒子子步设成流场网格的特征长度除以颗粒最大可能速度的几分之一宁小勿大。4. 轨迹图之外铜颗粒和纤维颗粒的运动形式到底差在哪4.1 相位滞后、振幅衰减与净漂移等仿真跑完首先看颗粒水平位移随时间变化的曲线。把油流入口速度曲线叠加上去对比就能直观看出颗粒对往复流场是否忠诚执行指令。我实测下来10微米的纤维颗粒基本上贴在流线上的震荡路径上位移曲线与油流速度波形几乎同步相位滞后肉眼几乎不可见。这是因为它弛豫时间只有约1e-7秒远小于流动周期。而100微米铜颗粒就明显不一样了位移曲线相对油流速度出现了明显的相位滞后而且往复振幅也偏小。原因是颗粒惯性在每一个换向时刻都试图维持原来的运动状态振荡流中表现为振幅衰减相位滞后频率越高这两个效应越明显。更有意思的是颗粒的净漂移。往复流场理论上是零净通量的但颗粒如果在垂直方向存在重力沉降就会和水平的振荡流动产生耦合形成一种类似走楼梯的轨迹每往复数周期颗粒就朝某个方向偏移一点。铜颗粒因为沉降速度快这种楼梯式净漂移非常明显纤维颗粒沉降弱净漂移就小得多基本在流线附近来回振荡。这个现象直接指导了危险区域风险评估铜颗粒更容易从油流中被甩到壁面附近长期滞留纤维颗粒则容易被油流循环带走。4.2 沉积位置与壁面碰撞频率把粒子轨迹按颜色区分后能看到非常清晰的模式铜颗粒几乎都倾向聚集在油道底部电极表面特别是在速度幅值低、重力主导的时段沿壁面形成一条沉积带。而在往复流换向的瞬时速度极低点颗粒几乎只受重力作用有足够时间沉降到壁面。纤维颗粒则更加分散它们倾向于追随流线在无束缚的区域来回穿梭只有遇到网格涡或壁面旋涡时才会被卷到近壁区停留。如果对壁面碰撞次数做统计纤维颗粒在单位时间内的碰撞次数往往比同粒径铜颗粒更高因为它们在流场里待的时间更长、运动范围更大自然更容易碰到电极和障碍物。这个结论对局部放电风险评估很重要不一定所有沉积颗粒都危险恰恰是那些反复碰撞、在电极附近来回穿梭的颗粒更容易激发微放电。4.3 把结果量化而不是只看动画轨迹图很直观但项目汇报需要定量结论。我建议至少输出以下几类数据不同粒径颗粒在一个周期结束时的净位移向量画成粒径-净位移曲线能快速看出哪个粒径范围颗粒最难被油流带走。电极表面颗粒沉积数量随时间的变化曲线斜率直接对应沉积速率。颗粒在指定危险区域比如电极间隙中心的平均逗留时间这个可以用粒子追踪的时间积分表达式统计。如果开了电场耦合再统计高场强区的颗粒出现频率和局放风险直接挂钩。以10微米纤维和100微米铜颗粒为例我的典型结果是100微米铜颗粒在电极附近的沉积比例远高于10微米纤维但单个纤维颗粒在油域中的存活时间和穿越频率更高。这说明两种颗粒的危险模式完全不同一个偏沉积型一个偏穿梭型评估时不能混为一谈。5. 运行中反复踩的坑和我的修正经验5.1 颗粒莫名其妙穿出油域或穿进壁面第一次跑通模型时我发现部分颗粒会在接近壁面时直接穿过好像壁面不存在。排查后确定是粒子时间步太大流场网格在近壁边界层非常薄颗粒在一个子步内的位移比边界层厚度还大穿透就在所难免。修正办法有两条一是在粒子追踪设置里限制最大粒子子步我一般按网格最小尺寸的1/3除以颗粒最大可能速度来估二是提高近壁网格分辨率确保边界层至少有两到三层网格托住颗粒运动。还要在结果里统计粒子数守恒COMSOL有粒子数量统计功能如果穿出壁面的颗粒占比超过几个百分点基本可以判定步长或网格有问题。5.2 入口速度突变导致压力震荡直接把入口速度设成U0*sin(2*pi*f*t)t0时速度是0按理说没有阶跃但压力和速度场在第一个时间步很容易出现非物理的高频振荡。原因在于初始场假设静止而速度边界开始变化后压力波需要一段时间传遍整个计算域如果初始步长过大压力修正跟不上。更稳妥的办法是在初始半个到几个周期内加斜坡系数比如写成U0*sin(2*pi*f*t)*min(t/0.05, 1)或者用U0*sin(2*pi*f*t)*(1-exp(-t/0.02))。让入口速度从0软启动到目标振荡状态能明显改善压力场初期的数值噪声。类似问题在交变电场中也一样电压也建议用软启动函数不让电场在第一步就满负荷加载。5.3 纤维颗粒用球形近似失真COMSOL默认颗粒是球形Stokes曳力公式也基于球形假设。纤维颗粒实际上是细长体长径比可能到几十用同体积球形去模拟曳力和旋转响应都会偏差。我在这块走过弯路直接用了球形模型跑纤维结果颗粒跟随性偏强、沉积率偏低和实验现象对不上。后来改用等效直径加形状因子的修正方案在曳力公式里引入动态形状因子对于纤维状颗粒形状因子可取1.5到3视长径比而定。COMSOL里可以在粒子追踪的曳力设置中选用户定义手动添加形状因子修正。如果你有细长体多段链模型也可以用多个球体组合成纤维但那对网格和计算量要求高得多一般项目没必要。5.4 流场和粒子追踪的求解时间步不一致导致的假振荡瞬态求解器采用自适应步长流场大时间步时粒子追踪虽然可以内部细分步长但如果全局时间步进过大粒子在每个时间步内看到的流场变化被线性化处理高速往复下会丢失真实的速度变化细节轨迹看起来就像假振荡或者死板直线。解决方法要么是把整个求解器时间步限制得足够小比如固定为流动周期的1/100到1/200要么是在粒子追踪中明确设置允许粒子子步和最大子步数。必须做网格无关性和时间步无关性验证把时间步减半如果颗粒位移曲线变化很小说明时间分辨率够了。5.5 没有实验数据时怎么给结果背书现场没有专门的透明油道颗粒可视化实验任何人都会问一句你的仿真凭什么可信。我的经验是至少做三层自检第一层是网格无关性验证粗中细三套网格比对颗粒轨迹和沉积率确认结果不依赖网格密度。第二层是极限工况对比把往复流频率降到接近0模型应该退化为稳态流动此时颗粒轨迹应该近似Stokes沉降解可以用沉降末端速度的解析公式直接对比把流场关掉只留重力颗粒自由沉降速度也应该与理论值一致。第三层是物理守恒检查统计粒子总数、质量、能量变化确保没有数值质量凭空消失。这三层自检过关后即使没有实验数据仿真结果在工程决策里也够用了。去年我就是依靠这套模型把何种粒径颗粒会在换流变阀侧油道中滞留并靠近高场强区的结论输送给绝缘设计团队后续油样实测中发现的颗粒粒径分布区间与模拟预判基本吻合算是一次比较成功的仿真加实测闭环。最后分享一个我固定使用的习惯把往复频率、速度幅值、粒径这三个参数做成扫描一次性跑完一组参数化研究再把不同条件下的净位移和沉积率汇总成表。这个动作看起来不起眼但做参数敏感性分析时省掉的重复操作时间非常可观。真正做这类颗粒运动仿真一两个孤立的算例说服力有限能说明在什么参数区间内颗粒行为会发生转折的结论才最有工程价值。