2. 多场耦合模型的可信度挑战

单物理场模型的可信度验证，核心问题通常是“网格够不够细、材料参数准不准、边界条件对不对”。一旦进入多物理场耦合，验证的难度会发生质变：误差不再只沿单一求解链传播，而是会在物理场之间来回反射、相互放大，甚至出现“每个单场都验证通过，耦合结果却不可信”的情况。理解这种质变的来源，是建立多场模型可信度评估体系的前提。

一、多物理场耦合的五种基本类型

从耦合的数学结构和实现方式看，COMSOL 中的多物理场耦合大致可归为五类，它们对可信度的威胁方式各不相同。

• 单向耦合：场 A 的解作为场 B 的输入，B 不影响 A。例如焦耳热中电流产生焦耳热，热场反过来对电导率的影响被忽略。可信度风险主要来自“忽略反向影响是否合理”，而非数值耦合本身。
• 双向耦合：两个场互相影响，必须联立求解或迭代求解。例如热膨胀中温度场改变结构变形，变形又改变传热路径。误差会在两个场之间往返传播。
• 强耦合（全耦合）：耦合项直接进入同一方程组，如压电效应中应力与电场通过本构矩阵同时求解。耦合项误差直接进入雅可比矩阵，影响收敛与精度。
• 弱耦合（分离式）：各场分别求解，通过外迭代交换数据。例如流固耦合中流体求解器与结构求解器交替迭代。可信度风险集中在迭代是否真正收敛、交换量是否守恒。
• 间接耦合：通过中间量或第三方场传递，如电磁—热—结构三级链。误差在链上逐级累积，且中间场的验证往往被跳过。

判断一个模型属于哪一类，不能只看物理现象，而要看实际求解设置：是同一个雅可比矩阵，还是分离步；是单向传递，还是双向迭代。同一物理问题在不同求解策略下，可信度风险完全不同。

二、耦合项误差传播：误差如何被放大

耦合项是连接两个物理场的“桥梁”，也是误差最容易放大的位置。以热—结构耦合为例，热应变项 ε_th = α(T)(T − T_ref) 中，热膨胀系数 α 本身随温度变化，若在模型中将其简化为常数，则温度每升高 100 K，α 的取值偏差可能带来 5%～15% 的应变偏差；而应变又通过结构刚度影响位移，位移再通过几何非线性影响传热边界，形成闭环放大。

更隐蔽的是耦合项的符号与量纲一致性。在多场方程中，源项、汇项、通量项的方向约定必须统一。例如在电化学—热耦合中，若电流密度方向与热源符号约定不一致，焦耳热可能被错误地写成负值，导致温度场反向变化。这类错误在单场验证中完全看不出来，只有做耦合守恒检查才会暴露。

三、参数不确定性在耦合中的非线性放大

多场模型中，参数不确定性不是简单相加，而是通过耦合关系非线性放大。以温度相关材料为例：热导率 k(T) 的不确定性影响温度分布，温度分布又影响热膨胀，热膨胀影响结构应力，应力通过压阻效应影响电导率，电导率再影响焦耳热。这条链上，每个参数的敏感度都会被前一级的梯度放大。

实际验证中，可以用局部敏感度 S_i = ∂Y/∂p_i · p_i/Y 来量化。若某耦合路径上敏感度乘积大于 1，说明该路径是误差放大器，应优先做参数标定或不确定性量化。对于强非线性耦合，建议至少对关键参数做 ±10% 的扰动扫描，观察输出响应是否线性；若明显非线性，则不能用简单误差传递公式估计可信区间。

四、边界条件交互与尺度差异

多场模型的边界条件往往不是独立的。例如流固耦合中，流体侧的压力边界与结构侧的位移边界在交界面上必须满足力和位移的连续性；若两侧分别设置，容易出现“过约束”或“欠约束”。过约束会导致求解器报错或收敛到非物理解，欠约束则会让刚体位移混入结果。

尺度差异是另一个常被低估的挑战。电磁场中的趋肤深度可能在微米量级，而结构热变形在毫米量级，两者相差三个数量级。若用同一套网格同时解析，网格量会爆炸；若用不同网格并通过映射耦合，则映射误差会引入额外的守恒误差。验证时应检查：映射前后总通量是否守恒、界面上的物理量是否连续、粗网格侧是否出现非物理振荡。

五、面向可信度的验证策略

面对上述挑战，验证不能停留在“单场残差小”的层面，而应建立三层检查。第一层是耦合守恒检查：在交界面上分别计算两侧的通量（热流、电流、力），确认差值在网格收敛后趋于零。第二层是耦合项敏感度扫描：对关键耦合参数做 ±10% 扰动，记录输出变化，识别误差放大路径。第三层是降阶对比：用简化解析解或单场近似解作为参照，确认耦合结果在物理上合理。

一个可操作的检查清单如下，可直接用于多场模型的可信度自检：

需要强调的是，多场模型的可信度不是“一次验证通过”就一劳永逸的。当几何、材料、工况发生变化时，耦合路径的敏感度会改变，原先合格的判据可能失效。因此，可信度验证应作为模型迭代的一部分，随模型演进而重新执行关键检查项。