2. 多物理场控制方程与耦合形式

多物理场问题的本质，是多个物理场在同一空间区域内同时演化，并通过共享的未知量、材料参数或界面条件相互影响。与单物理场仿真相比，多物理场建模的难点并不在于单个方程本身，而在于耦合项如何进入方程、耦合强度如何影响求解策略、以及界面条件如何保证物理量守恒。本节围绕偏微分方程耦合系统、强耦合与弱耦合、单向与双向耦合、界面耦合条件与整体求解框架展开。

一、偏微分方程耦合系统的基本形式

设求解区域为 Ω，未知场向量记为 u = [u₁, u₂, …, u_n]^T，每个分量代表一个物理场，例如温度 T、位移 u、压力 p、电势 φ、浓度 c 等。多物理场控制方程可以统一写成如下耦合系统：

L₁(u₁; u₂, u₃, …) = f₁   在 Ω₁ 内
L₂(u₂; u₁, u₃, …) = f₂   在 Ω₂ 内
…
B(u₁, u₂, …) = g         在界面 Γ 上

其中 L_i 是第 i 个物理场的微分算子，f_i 是源项或外载荷，B 表示界面耦合算子，g 为界面给定条件。关键点在于：L_i 中出现了其他场的未知量，这正是耦合的来源。例如热弹性问题中，温度场通过热膨胀项进入力学平衡方程，而力学变形又通过热弹耦合项影响能量方程；流固耦合中，流体压力作为面载荷作用于固体，固体位移又改变流体计算域。

从数学结构看，耦合系统可分为三类：体积耦合（耦合项出现在区域内部方程中，如焦耳热、热膨胀）、界面耦合（耦合通过交界面上的通量或位移连续条件实现，如流固界面、接触界面）、全局耦合（通过电路、外回路等集中参数元件实现，如电磁-电路耦合）。

二、强耦合与弱耦合

强耦合与弱耦合描述的是耦合项在离散方程组中的处理方式，而不是物理耦合的强弱。

强耦合（ monolithic / fully coupled ）：将所有物理场的未知量组装进同一个大型方程组，一次性求解。以热-结构耦合为例，离散后得到：

[ K_TT   K_Tu ] [ ΔT ]   [ R_T ]
[ K_uT   K_uu ] [ Δu ] = [ R_u ]

其中 K_Tu 和 K_uT 为耦合雅可比块。强耦合的优点是数值稳定性好，适合强非线性、强相互作用问题；缺点是矩阵规模大、内存与计算成本高，且需要为耦合块提供一致的雅可比矩阵。

弱耦合（ partitioned / staggered ）：将各物理场分开求解，通过交换界面数据实现耦合。典型流程是：先求解场 A，把结果作为边界条件或源项传给场 B，再求解场 B，必要时迭代若干次。弱耦合的优点是模块化、可复用已有单场求解器、内存占用低；缺点是当耦合很强时可能收敛慢甚至发散，需要松弛因子或 Aitken 加速。

三、单向耦合与双向耦合

单向耦合指一个场影响另一个场，但反向影响可忽略。例如：已知温度分布，计算热应力，而不考虑变形对温度的反作用。单向耦合只需一次传递，计算量小，适合热应力、压电传感等弱反作用场景。

双向耦合指两场相互影响，必须迭代或联立求解。例如：流固耦合中，流体压力使固体变形，固体变形又改变流道形状和流场；电磁-热耦合中，电流产生焦耳热，温度又改变电导率。双向耦合的收敛判据通常取界面物理量的相对变化：

|| u_Γ^(k+1) - u_Γ^(k) || / || u_Γ^(k+1) || < ε

其中 ε 常取 1e-4 ~ 1e-6，具体取决于工程精度要求。

四、界面耦合条件

界面耦合条件是多物理场在交界面上的“物理契约”，必须同时满足运动学连续和动力学平衡。以流固耦合界面 Γ_FSI 为例：

• 运动学条件：流体速度等于固体速度，v_f = ∂d_s/∂t；
• 动力学条件：流体应力与固体应力满足 σ_f·n = σ_s·n，即牵引力连续；
• 几何条件：界面位置由固体位移更新，网格需随动或重划。

热-结构界面则要求温度连续 T₁ = T₂ 和热流连续 q₁·n = q₂·n。电磁-热界面要求电势连续和电流密度法向连续。界面条件写错，往往表现为能量不守恒、界面温度跳变或计算发散。

五、整体求解框架

一个完整的多物理场求解框架通常包含以下层次：几何与网格层负责区域划分和界面标记；物理场层定义各场控制方程与本构关系；耦合层负责界面数据映射、耦合项组装与收敛控制；求解层选择强耦合或弱耦合策略；后处理层输出场分布与耦合指标。

下面用一张结构图展示多物理场耦合系统的整体求解框架，重点表达各物理场、耦合层与求解器之间的数据流与界面条件传递关系。

该框架的核心思想是：物理场层定义“各自怎么算”，耦合层定义“彼此怎么传”，求解层定义“整体怎么解”。在实际建模中，先明确哪些界面存在耦合、耦合是单向还是双向、耦合强度是否要求强耦合，再选择求解策略，可以显著减少试错成本。

总结来看，多物理场控制方程与耦合形式的掌握，需要同时回答四个问题：方程里耦合项写在哪、界面条件是否守恒、耦合方向是单向还是双向、求解时是联立还是分步。这四个问题决定了后续模型降阶与快速求解的可行路径——例如弱耦合结构天然适合分块降阶，而强耦合系统则需要保留耦合雅可比块的低秩近似。