在玻璃上沉积一条蛇形电阻膜，给两端施加电压，就能做成一个隔着玻璃加热流体的小装置。真正建模时，难点不只是把电压设成12 V：薄膜只有10 μm厚，玻璃却有2 mm厚。电流在薄膜内产生的热，应该作为三维域热源输入，还是作为玻璃表面的热通量？如果结果变量已经是 W/m²，还要不要再乘膜厚？
这篇围绕 COMSOL 6.3 的 Heating Circuit（模型 ID 465），把薄膜表示方法、焦耳热耦合和功率积分连起来。完整案例还包含 Solid Mechanics 与热膨胀，用于分析应力和变形；这里聚焦电热部分，不展开整套热应力设置。
一、先看几何：薄膜不是一块“没有厚度”的材料
装置由玻璃板、蛇形镍铬电阻膜与两端银焊盘组成。玻璃板尺寸为130 mm×80 mm×2 mm；镍铬膜宽5 mm、厚10 μm；银焊盘为10 mm×10 mm，厚度同样是10 μm。沉积回路的一侧接触空气，玻璃背面接触被加热的流体，玻璃侧边则按热绝缘处理。
图1：加热装置与空气、玻璃、流体的位置关系。来源：COMSOL 6.3 Heating Circuit 教程图1，第2页。
玻璃在这里提供电路与流体之间的隔离，并承担导热作用。薄膜产生的热既可以传给空气，也可以进入玻璃，再通过背面与流体换热。最终温度由整个热通路决定，不能只看电阻膜发了多少热。
图2：蛇形镍铬回路与两端银焊盘的几何。来源：同版本教程图2，第3页。
若把10 μm薄膜完整拉伸为三维实体，再与毫米级玻璃一起划分网格，会出现明显的尺度差异。这里采用层状壳表示电流，在固体传热中使用 Thin Layer 表示薄膜。它把薄膜放在几何边界上，通过材料层数据保留厚度和物性，从而避免必须沿真实膜厚建立三维网格。
这不等于把厚度设为零。电流导通能力、面发热和薄层导热仍依赖真实的 d_layer。建模中省去的是显式三维薄膜实体，保留的是薄层的物理作用。
二、Thermally thin 近似究竟简化了什么
在 Heat Transfer in Solids 下添加 Thin Layer，边界选择 Circuit，Layer type 采用 Thermally thin approximation。这里的近似针对薄膜厚度方向的温差，并不要求整条蛇形回路温度相同。沿回路及平面方向仍然可以存在温度梯度，转弯和中央位置也可以呈现不同温度。
判断这种近似是否适用，要看厚度方向温差是否足够小，以及研究目标是否需要解析膜内的厚度温度分布。仅仅看到“10 μm”就认定任何工况都可用薄层近似，并不充分。对于厚度方向温差显著、接触热阻重要或具有复杂多层结构的装置，需要重新检查近似及层间热传递表示方法。
电学使用 Electric Currents in Layered Shells，标识为 ecis，也选择 Circuit。它描述薄膜中的电流传导；玻璃作为三维固体域参与传热。两者作用在不同几何层级，不能把薄膜壳接口误当成已经在整个玻璃体积中求解导电。
在 Definitions 中，Circuit 是边界层级的显式选择，包含边界6—8。把它同时用于壳电流与传热薄层，有助于保持两者作用范围一致。边界编号只对应这份几何，自己重建图形后应核对位置，不应只复制数字。
三、材料与电边界：银焊盘和镍铬不能混成一层
参数表中设置 V_in = 12 V、d_layer = 10 μm，银的电导率 sigma_silver = 6.3×10⁷ S/m，镍铬的电导率 sigma_nichrome = 9.3×10⁵ S/m。这些电导率采用常数，不包含电阻随温度改变的反馈。
玻璃体积使用 Silica glass。薄膜材料采用 Single Layer Material：先建立 Silver Layer，选择整个 Circuit，Position 设为 Bottom side on boundary；再建立 Nichrome Layer，只选择边界7，并同样定义层位置。镍铬材料在蛇形回路区域覆盖银材料分配，两端保留银焊盘。不要把这个操作误解成在整条回路上叠了一层银和一层镍铬。
两种薄层材料中的 Thickness 都使用 d_layer。厚度若在材料处和其他手工表达式中被重复处理，后面的功率就可能出现数量级错误。热导率也不同：银为420 W/(m·K)，镍铬为15 W/(m·K)，玻璃为1.38 W/(m·K)，因此材料既影响发热分配，也影响热量在结构中的扩散。
电势施加在壳的边缘：Edge 43 使用 Ground，Edge 10 使用 Electric Potential，V0 输入 V_in。这里是边缘电条件，不是把玻璃某个三维面设成电极。随后添加 Continuity，Source 选 Silver Layer，Destination 选 Nichrome Layer，使两种材料层在连接位置满足相应连续条件。
材料选区、厚度、电极位置与层间连续性要一起检查。仅仅看到两端有12 V和0 V，还不足以证明电流已经沿预期的银焊盘—镍铬蛇形膜通路流动。
四、最关键的单位转换：W/m³ 怎样变成 W/m²
Q_DC = J·E = σ|∇tV|²，单位 W/m³。
J 是电流密度，单位 A/m²；E 是电场强度，单位 V/m；∇tV 表示沿薄层切向的电势梯度。J·E 描述每单位体积产生多少热功率，还不是整个回路的功率。
q_prod = d × Q_DC，单位 W/m²。
其中 d 是膜厚，单位 m；这个简式对应厚度方向采用相应均匀近似的情况。可以直接检查量纲：m × W/m³ = W/m²。薄膜被等效到边界后，这份热量以面发热的形式传给热模型，物理上并没有少算薄膜体积。
更一般地，应把体积损耗沿厚度积分；只有本例这种均匀厚度、采用相应薄层近似的情形，才可直接写成厚度与体积密度的乘积。对于厚度不均或损耗沿厚度显著变化的薄膜，不能用一个固定厚度代替实际积分。
但不能把这条公式机械地套到每一个变量上。ecis 提供的 ecis.Qsh 是 Surface loss density, electromagnetic，单位已经是 W/m²。它已经属于面功率密度，计算回路总发热时应对 Circuit 做面积积分：
P_gen = ∫Circuit ecis.Qsh dA，单位 W。
如果先将 ecis.Qsh 再乘 d_layer，再做面积积分，量纲会变成 W·m，而不是功率。这是重复计算厚度，不是让结果“更完整”。另一方面，若使用的是 Q_DC 这种 W/m³ 体积密度，直接做面积积分又缺少一个长度因子。判断该不该乘厚度，先看变量定义和单位，再看当前几何维度。
模型通过 Electromagnetic Heating, Layered Shell 耦合，把壳接口计算的电阻损耗接入传热。已有这项耦合时，不应再把同一份 q_prod 或 ecis.Qsh 手动添加一次，否则可能把焦耳热重复计入。热源总量不合理时，应同时检查厚度与热源是否重复，而不是先修改换热系数去压低温度。
热量首先在薄层中产生，随后向空气与玻璃两侧分流。等效为边界上的发热，不代表所有电功率都进入玻璃；玻璃接收到的热量，也不代表最终全部进入流体。热源、跨界面传热与外边界净散热，是能量路径上的不同量。
五、空气与流体是两组换热边界
空气侧 Heat Flux 选择边界4、6—8，Flux type 采用 Convective heat flux，h 输入 h_air = 5 W/(m²·K)，外部温度输入 T_air。参数表将 T_air 定义为20 °C，即293.15 K；教程模型说明中的293 K是取整表达。
流体侧 Heat Flux 只选择背面边界3，h 输入 h_fluid = 20 W/(m²·K)，外部温度输入 T_fluid = 353 K。其余侧边按热绝缘处理。两个换热系数的单位都是 W/(m²·K)，不能当成空气或流体的导热系数 W/(m·K)。
这里的353 K是边界输入，不是计算得到的出口温度。模型用它和换热系数概括流体对玻璃背面的热作用，没有求解流体速度场、沿程温升或出口混合温度。因此，得到向流体传热的功率，也不能直接说已经获得某一流量下的加热效果。
约定向外散热为正时，边界热通量可理解为 q_out = h(T−T_ext)。局部固体温度高于外部温度时，热量向外传递；若低于外部温度，该区域也可能吸热。积分应保留这种方向性，不能把所有局部热通量取绝对值后相加，再称为净传热量。
六、从面损耗图走到总功率和温度
图3绘制的是 ecis.Qsh，单位 W/m²。在12 V工况下，回路转弯内侧具有较高的电流密度，因此局部面发热更强。它说明哪里单位面积发热集中，不能把色标最大值直接读成整个装置的瓦数。
图3：电阻薄膜的电磁面损耗密度 ecis.Qsh，单位 W/m²。来源：COMSOL 6.3 Heating Circuit 教程图3，第5页。
后处理可复制解数据集，增加 Boundary Selection 并选择 Circuit；在三维绘图组中添加 Surface，表达式选择 ecis.Qsh。计算总量则在 Derived Values 中添加 Surface Integration，选 Circuit，积分同一个变量，得到总发热约13.8 W。面积积分的选区必须覆盖整条回路，而不是只选最热转角。
温度分布显示，最高温度约428 K，位于回路中央区域。注意，最强局部发热在转角内侧，而最高温度出现在中央，两者并不矛盾。温度除了受局部热源影响，还受邻近区域导热、边界散热和到散热面的热通路影响，因此热点位置不能只靠损耗图判断。
图4：回路与玻璃板的稳态温度分布，单位 K。来源：同版本教程图4，第6页。
流体侧边界3上的净传热量大小约8.5 W。这里比较的是净传给流体的功率与总电发热，并不是把8.5 W当作某一点的热通量。若后处理选择 ht.q0，变量名称是 Inward heat flux，应先核对它的向内正号约定；采用向外为正汇总能量时，需要统一转换符号，不能直接把带方向的积分值当作无方向数值。
两项功率的差额不能称为“被玻璃消耗掉的热量”。玻璃的热阻改变温度分布与热量分配，热量本身不会因为经过热阻而消失。在侧边绝热、稳态且没有其他热源的设置下，还需要把空气侧净散热纳入平衡，才能判断总能量是否闭合。
七、网格与核查应该围绕这个建模选择
网格在边界4、6—8上建立 Free Triangular，对 Circuit 设置最大单元尺寸2 mm；玻璃体积使用 Swept，并在厚度方向分成3层。薄膜厚度由层数据承载，不应误以为玻璃的3层扫掠网格就是把10 μm电阻膜划成了3层实体单元。
转角电流集中与局部温度梯度值得进一步检查。做网格对照时，应比较同一位置的温度、同一选区上的总发热及面损耗变化，不能只凭图像平滑认定网格足够。教程中的1e-3设置针对位移自由度的 Scaling，用来改善完整模型的求解表现，并不是电热求解的相对容差，更不能当成结果误差。
总功率是空间积分量，局部峰值则对转角附近的分辨率更敏感。一个指标变化较小，不能自动证明另一个指标也稳定。需要先确定设计关注的是总供热还是局部过热，再选择对应的比较指标，避免用平均表现掩盖局部风险。
可以把核查分成三个问题：面发热积分是否具有 W 的正确单位；电热耦合是否只计入一份损耗；稳态总发热是否与空气侧和流体侧的净散热相匹配。核对能量时要统一选区与符号，侧边绝热也应与几何选择一致。这里是稳态平衡，不含随时间变化的储能项；若改成瞬态，就不能省略材料吸热储能。
对于自己的器件，还要重新判断电导率温度依赖、膜与玻璃之间的接触热阻、真实流体换热条件等假设。完整案例中的结构场可以分析热应力，但温度和功率正确并不自动证明薄膜不会脱层，也不能据此保证长期可靠性。
这类模型最值得记住的，是“变量单位—几何层级—积分方式”必须一致。薄层可以保留真实厚度而不建立三维薄膜域；体积焦耳热经过厚度积分成为面发热；已经以 W/m² 给出的壳损耗，再做面积积分才得到 W。把这条链理顺，后续的温度图和能量比较才有清楚的物理基础。
