COMSOL仿真绝缘油中铜/纤维颗粒往复运动轨迹与沉积规律

去年做换流变阀侧绝缘可靠性评估的时候，被一个问题反复纠缠：如果绕组里掉出几颗铜屑，或者绝缘纸板边缘剥落了几根细小纤维，这些颗粒在绝缘油里到底怎么跑？油泵启停、负载周期性波动带来的油流振荡，正好对应了往复运动这个工况。为了把这个问题彻底摸清，我用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³。我按这个估算了几组数：

再取流动特征时间，往复运动下取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 几何、油品和颗粒参数的基准值

我建议所有参数都用全局定义管理，这样后面的参数化扫描会很方便。基准参数如下：

初始释放位置我一般选在入口附近均匀分布，或者直接在油域里随机均匀分布，前者能看入口效应，后者能更直观反映全局行为。颗粒数量不要一上来就放几千个，一方面是后处理轨迹图会糊成一团，另一方面是粒子追踪的计算时间与颗粒数量强相关，先用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)，t=0时速度是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沉降解，可以用沉降末端速度的解析公式直接对比；把流场关掉只留重力，颗粒自由沉降速度也应该与理论值一致。第三层是物理守恒检查，统计粒子总数、质量、能量变化，确保没有数值质量凭空消失。

这三层自检过关后，即使没有实验数据，仿真结果在工程决策里也够用了。去年我就是依靠这套模型，把"何种粒径颗粒会在换流变阀侧油道中滞留并靠近高场强区"的结论输送给绝缘设计团队，后续油样实测中发现的颗粒粒径分布区间与模拟预判基本吻合，算是一次比较成功的仿真加实测闭环。

最后分享一个我固定使用的习惯：把往复频率、速度幅值、粒径这三个参数做成扫描，一次性跑完一组参数化研究，再把不同条件下的净位移和沉积率汇总成表。这个动作看起来不起眼，但做参数敏感性分析时省掉的重复操作时间非常可观。真正做这类颗粒运动仿真，一两个孤立的算例说服力有限，能说明"在什么参数区间内颗粒行为会发生转折"的结论才最有工程价值。