3. 验证前的模型审查

在COMSOL多物理场建模中，模型一旦求解完成，结果往往看起来“有云图、有曲线”，但这并不等于模型可信。验证工作的第一步不是对比实验数据，而是对模型本身做一次系统审查。审查的目标是回答一个根本问题：这个模型所求解的方程、所用的参数、所设的边界与初始条件，是否真实反映了我们要研究的物理对象？如果这一层存在结构性错误，后续任何网格收敛性分析或实验对比都失去意义。本节围绕五个审查维度展开，并给出可操作的检查方法。

一、控制方程与物理场选择合理性

COMSOL的每个物理场接口背后都对应一组特定的控制方程。选择接口时，首先要确认方程形式与物理过程匹配。例如，传热问题中，若固体内部以导热为主、流体区域存在明显流动，就应使用“共轭传热”而非单纯“固体传热”；若流速较高、惯性效应不可忽略，则非等温流动接口比层流接口更合适。审查时建议逐条核对：

• 方程是否包含所研究的主导效应（对流、扩散、辐射、相变、多孔介质渗流等）；
• 是否引入了不必要的强耦合导致求解困难或虚假稳定；
• 因变量与材料属性是否在同一坐标系和参考系下定义；
• 多物理场耦合节点是否真正传递了正确的物理量（如温度→黏度、压力→密度）。

一个常见错误是：在自然对流问题中使用了“层流”接口但未开启重力体积力，导致浮升力缺失，流场完全错误。审查时应检查体积力、源项是否被显式加入。

二、材料参数来源审查

材料参数是模型可信度的“输入质量”。审查时需为每个关键参数建立来源档案，至少区分以下层级：

特别要注意温度依赖参数。例如硅钢片磁导率随温度变化、聚合物黏度随剪切速率变化，若在模型中设为常数，需在审查记录中明确说明该简化的适用温度区间。对于示例性教学模型，若采用“假设导热系数为0.2 W/(m·K)”这类设定，应在文档中标注“示例假设值，非实测”。

三、几何简化假设审查

几何简化是建模中最容易被忽视的环节。常见的简化包括：忽略倒角与圆角、将三维轴对称问题降为二维轴对称、用多孔介质等效复杂内部结构、截断远场区域。审查时需逐项确认：

• 简化是否改变了主导传热/传质路径？例如去掉散热翅片间的细小间隙可能显著改变对流换热面积；
• 对称面是否真实存在？若边界上存在非对称热源或流动，使用对称简化会强制错误结果；
• 截断边界距离是否足够远？一般要求远场边界距扰动源至少5～10倍特征长度；
• 薄层结构是否用“薄壳”或“过渡边界条件”替代实体网格，并检查其数学一致性。

四、边界条件物理意义审查

边界条件不是数学装饰，而是物理环境的抽象。审查时应对每个边界逐一追问：它代表真实世界中的什么？例如：

• 热绝缘边界：代表无热量交换，还是仅仅因为“不知道设什么”而默认？若实际存在自然对流散热，绝缘假设会高估温度；
• 固定温度边界：代表恒温浴、大热容基底还是理想散热器？需确认其热阻是否可忽略；
• 流出/流入边界：在流体中，压力出口与速度出口的物理含义不同，前者允许回流，后者强制方向；
• 对称边界：同时约束了法向通量为零和切向梯度为零，是否与真实对称条件一致。

一个实用的审查方法是：把每个边界条件翻译成一句自然语言，写在模型文档中。如果翻译后与实验描述矛盾，则边界条件需要修改。

五、初始条件一致性检查

初始条件对瞬态问题尤为关键。审查要点包括：

• 初始值是否与边界条件相容？例如初始温度设为300 K，而入口边界立即设为500 K，会在入口附近产生非物理的陡峭梯度，需通过渐变加载或延长入口段缓解；
• 初始场是否满足约束方程？例如压力初始场若与密度、温度不满足状态方程，求解器可能在第一步就报错；
• 多物理场初始值是否相互一致？流场初始速度为零而温度场已有强浮升力时，可能引发瞬态振荡。

对于稳态问题，初始条件只影响收敛路径，但仍建议用物理上合理的初值以改善收敛。审查时可先求解一个简化物理场（如只算导热）作为初值，再耦合完整物理场。

审查记录与执行建议

上述五个维度建议以检查表形式固化到项目文档中。每个维度至少记录：检查项、结论（通过/不通过/待确认）、依据或假设、责任人。对于教学示例模型，尤其要区分“实测参数”与“假设参数”，避免读者误将示例值当作真实材料数据。

实际操作中，可按以下顺序推进：先确认控制方程与物理场接口，再核对材料参数来源，然后审查几何简化是否影响主导物理过程，接着逐边界翻译物理意义，最后检查初始条件相容性。任何一项不通过，都应返回修改并重新审查，而不是带着已知缺陷进入验证对比。只有审查通过后，网格收敛性分析和实验/解析对比才具有可信度基础。