关于多物理场耦合仿真的相关思考

关于多物理场耦合仿真，写点自己的思考。

1 核心本质

多物理场耦合仿真，听起来是个挺高大上的名词。不少人被各种名词创新搞得云里雾里，不知所谓。

实际上，多物理场耦合仿真理解起来并不算复杂。搞清楚了本质，做多物理场耦合仿真就从根源上理顺了思路。

所谓仿真思路，就是处理三个问题：

1. 分析什么对象的什么问题？
2. 分析的问题可用什么方程求解？
3. 有哪些已知信息？

仿真本质都是解方程，仿真不同问题则为求解不同方程。各种算法，如有限元法、有限体积法、格子玻尔兹曼法等，都是求解方程的方法。

解方程，需要三类输入：

1. 方程的表达式
2. 求解区域
3. 边界条件

仿真的已知信息，包括边界条件和方程中的项与系数。

一图总结上述内容：

所谓多物理场仿真，简单讲即：由于需要分析的问题缺失部分重要输入信息且缺失的输入信息受到另外一个物理场的影响，因此需要另外一个物理场的求解结果来作为输入信息。

由此，可将多物理场仿真涉及到的物理场分为主物理场和次物理场。主物理场是需要分析的问题所直接涉及的物理场，次物理场为补充缺失的输入信息而引入。

举一个很常见的例子：风力作用下的结构振动，即风致振动。分析风致振动问题是高层建筑、桥梁、飞机机翼等领域必需事项。风致振动分析的问题为结构的振动规律，风力是影响结构振动的一个重要因素。

风致振动引发 1940 年美国塔科马海峡大桥坍塌（图源：geopop.it）

在结构线弹性变形假设下，结构振动的方程为：

linear format: m∙(∂^2 x ⃗)/(∂t^2 )+c∙(∂x ⃗)/∂t+k∙x ⃗=F ⃗(t)

方程右侧的外力项中，风力的变化规律很难预先得知，若人为指定为某个函数，仿真可信度必然成疑。对此，需要先利用 CFD 仿真获取风力的变化规律，即上述振动方程右侧的外力作用项信息。

2 耦合方式

不同物理场的求解对象和计算域通常是不同的，但是一定会有公共部分，如公共面、公共实体等。数据传递仅能传递两个物理场公共部分，且为接收方可作为输入信息的数据。

如上述的风致振动问题，CFD 仿真仅传递两者的公共部分（结构和空气的接触面）的压力数据，作为结构的作用力。结构其他部分不和空气接触，不受风力的作用。CFD 仿真的流动速度等物理量也并非结构振动仿真所需要的输入信息。

根据主物理场和次物理场之间是否存在相互影响，可分为单向耦合和双向耦合两类：

• 单向耦合：主物理场的变化不影响次物理场，仅需要次物理场的数据传递到主物理场
• 双向耦合：主物理场的变化会影响次物理场，需要考虑主物理场和次物理场的相互影响

如上述的风致振动问题，若不考虑结构振动对空气流动的影响，即为单项耦合；若考虑结构振动对空气流动的影响，即为双向耦合。