2. 多物理场耦合界面与边界类型解析
在多物理场仿真中，真正决定结果可信度的往往不是各物理场内部的求解精度，而是耦合界面上的物理量是否被正确传递、边界条件是否被正确施加。边界一致性检查的核心，就是确认"界面两侧说的是同一件事"。本节先厘清常见耦合界面的物理量传递关系，再区分各类边界条件的数学本质，最后给出每类边界在一致性检查中的关注重点。
一、耦合界面的物理量传递关系
耦合界面本质上是两个物理场共享的几何面，其一致性取决于"传递什么量、以什么方向传递、是否守恒"。常见四类耦合界面的传递关系如下：
- 流固耦合（FSI）：流体侧向固体传递压力
p（面力），固体侧向流体传递位移u与速度∂u/∂t。这是典型的双向传递，检查重点是界面网格是否匹配、压力方向是否沿界面法向、固体位移是否引起流体域网格变形。 - 热固耦合：温度场向结构传递温度
T（产生热应变与热应力），结构一般不再向温度场回传力学量（单向耦合）；若考虑摩擦生热或变形做功，则需回传热流密度。检查重点是温度在界面上的连续性、热膨胀参考温度的取值是否两侧一致。 - 电磁热耦合：电磁场向温度场传递焦耳热/介质损耗
Q（体热源，单位 W/m³），温度场向电磁场回传电导率σ(T)等温变材料参数。检查重点是热源是否落在正确的体区域而非界面、材料参数的温度依赖是否双向启用。 - 声固耦合：声压
p作用于结构表面，结构法向加速度aₙ回传给声场作为声源项。检查重点是法向方向约定是否一致、声学域是否为封闭腔体、界面是否设为"声-结构耦合"而非普通壁面。
需要强调一个贯穿全文的约定：本文所有界面法向 n 均定义为由主侧指向从侧，压力以压入界面为正。切换软件或切换主从侧时必须重新核对符号，否则会出现"压力方向反了但残差仍收敛"的隐蔽错误。
二、边界条件的数学类型与物理对应
耦合界面之外，各物理场自身的边界条件同样需要一致性检查。按数学形式可分为四类，它们在有限元中对应不同的施加方式，检查重点也完全不同。
| 类型 | 数学形式 | 物理实例 | 一致性检查重点 |
|---|---|---|---|
| Dirichlet（第一类） | 指定场变量值 u = u₀ | 固定约束、给定温度、壁面速度、声压节点 | 是否与相邻面重复约束同一自由度；给定值单位与量纲；约束方向是否与坐标系一致 |
| Neumann（第二类） | 指定法向通量 q·n = q₀ | 热流密度、面力、声质点速度、电流密度 | 通量方向与法向符号；零通量（绝热/自由面）是否被误设为默认；单位面积还是总面积 |
| Robin（第三类） | 通量与场值线性组合 q·n = h(u∞ − u) | 对流换热、辐射线性化、声阻抗边界 | 换热系数与参考温度是否两侧一致；辐射是否已线性化；阻抗实部是否为正 |
| 接触/界面 | 两侧场量连续或跳跃关系 | 热接触热阻、结构接触、薄层界面 | 接触对主从面选择；初始间隙；接触热阻是否重复计入 |
这里最容易出错的是 Robin 边界与对流/辐射的混用。以热分析为例，对流边界 q = h(T∞ − T) 与辐射边界 q = εσ(T∞⁴ − T⁴) 在数学上都属于 Robin 型，但辐射是非线性的，若软件中误将辐射按常系数对流处理，低温段误差可达数倍。检查时应确认：辐射是否启用了非线性求解或等效辐射换热系数，且该系数是否随温度更新。
三、一致性检查的关注重点与操作路径
把上述类型落到检查动作上，可按"几何—物理量—数值"三层推进。以下给出一个可执行的检查顺序，示例参数均为教学假设，用于说明判断依据。
- 几何层：确认耦合面在两侧模型中为同一组面（面积偏差建议小于 1%）。若流固界面面积在流体侧为 0.02 m²、固体侧为 0.0198 m²，说明存在 1% 的网格偏差，需检查是否漏选了边界面。
- 物理量层：逐项核对传递量的单位与符号。例如压力传递若流体侧单位为 Pa、固体侧误设为 MPa，则数值放大 10⁶ 倍，表现为"收敛但变形异常"。符号上按本文约定，压入界面为正，若固体侧反号，则变形方向整体反转。
- 数值层：检查松弛因子与收敛判据。双向耦合中，若界面残差在 1e-3 附近振荡而不下降，常见原因是松弛因子过大（如取 1.0），可降至 0.3～0.5 后重试；若残差单调下降但物理量不守恒，则应回到物理量层复查。
常见陷阱：同一几何面上同时施加了 Dirichlet 与 Neumann 条件（如既固定温度又加热流），求解器通常只取其一而不报错，结果看似合理却违背物理。检查时应导出边界条件清单，按面逐一比对，确认每个面只有一类主约束。
快速自查口诀：面要同、量要同、号要同、类要清。"面同"指耦合面几何一致，"量同"指单位量纲一致，"号同"指法向与正负约定一致，"类清"指每个面的边界类型唯一且与物理意图相符。四者齐备，边界一致性检查即可闭环。
需要说明的是，不同求解器对"接触"与"界面"的处理粒度不同：有的把接触热阻作为界面属性，有的作为薄层单元。检查前应先确认当前软件把该物理效应归入哪一类，再决定是查边界条件表还是查材料/单元属性，避免在错误的位置反复排查。