线圈电流密度计算

Comsol油浸式变压器电磁-温度-流体多物理场耦合仿真；可以得到变压器稳定运行时内部热点温度及油流速度分布

油浸式变压器那滋滋的电流声背后，藏着一场肉眼看不见的物理狂欢。今天咱们直接上硬菜，用COMSOL还原变压器稳定运行时内部的电磁-热-流三重奏。先甩个结果镇楼：仿真完成后能看到115℃的局部热点正被2.3m/s的油流疯狂降温，这可比看温度计刺激多了。

先来搭电磁场的基本框架。在AC/DC模块里选磁场接口，线圈设置有个小技巧——直接给绕组截面加载电流密度更省事：

J_peak = (I_rms * N_turns) / (A_coil * fill_factor) # 峰值电流密度

这里fill_factor取0.7是实测经验值，太满会影响散热。边界条件记得勾选磁绝缘，别让磁场乱跑。求解频率设50Hz，瞬态求解器跑三个周期基本就能稳定。

电磁场算出的焦耳热是后续温度场的火种。耦合方式很巧妙——直接把体积损耗密度q_dot传到固体传热接口：

// 电磁损耗到热源的传递 ht.q0 = emw.Qh/emw.rho // 单位体积发热量

注意这里的材料属性必须用温度相关函数，特别是绕组电阻率随温度变化的系数取0.004/K，这个数值偏差1%能让最终温度差5℃以上。

Comsol油浸式变压器电磁-温度-流体多物理场耦合仿真；可以得到变压器稳定运行时内部热点温度及油流速度分布

温度场和流场的双向耦合才是重头戏。打开非等温流动接口，设置自然对流时别忘了重力方向（Y轴负方向）。粘度模型用Carreau方程处理变压器油的非牛顿特性：

% 变压器油粘度模型 mu_inf = 0.01; % 无限剪切粘度 mu_0 = 0.12; % 零剪切粘度 lambda = 0.5; % 松弛时间 n = 0.3; % 幂律指数 mu_eff = mu_inf + (mu_0 - mu_inf)*(1 + (lambda*gamma_dot)^2)^((n-1)/2)

求解器配置有讲究：先用稳态求解器跑两轮，等流场基本稳定了再切到瞬态。时间步长从0.1秒开始，自适应调整更保险。遇到不收敛别慌，把流场初始速度设为层流分布能破局。

后处理阶段得玩点花样：在绕组表面创建探针组监测热点温度，用粒子追踪看油流路径。有个骚操作是在结果表中设置温度超过105℃自动标红，配合流线箭头图，整个散热过程跟看动作大片似的。

最后验证别忘了对比油顶温升：仿真结果比实测高3-5℃算正常，毕竟实际运行中绝缘纸的老化会略微影响散热。要是误差超过10℃，赶紧检查流场边界条件——八成是忘了设置散热片处的对流传热系数，或者泵循环的入口流速设成了固定值。