来都来了总得带点COMSOL干货再走吧~

来源公众号：COMSOL仿真交流

多物理场强耦合

边界非线性问题

我们从单场仿真说起。单场仿真中最典型的边界非线性问题，是力学分析中的接触问题。接触上就有约束，没接触上就没有约束。是否接触上显然与形变的位移有关。

事实上，边界非线性的概念可以被理解为描述某种反馈控制，所以在其他物理场分析中也非常常见。

比如加热管道的瞬态传热分析，边界上的热源不是一个恒定值，而是要根据管道中间某点的实测温度来实时调整，边界条件是温度的函数。另外，传热问题中我们常用的自然换热边界：热通量Flux = h*(T -Text)，这其实也会为计算带来非线性。

在多物理分析中这种现象更为普遍。

比如，入口流速可控的反应器中溶质的浓度分布问题。这是一个对流扩散的过程，入口处存在初始浓度，在扩散作用和对流作用下，随着时间的增长，浓度会在反应器中形成一定的分布。我们在反应器中某位置放置一个探测针，用于读取该点浓度值，并要求入口流速可以不断调整，以使探针处浓度值恒定在指定值。显然，入口的流速边界条件此时与浓度分布有关，而浓度分布受对流作用影响，流体的计算也依赖于入口边界的设定。

下图展示的就是一个浓度监控的流动入口速度PID控制模拟。我们可以想象入口速度将是一个先增大后减小，再增大再减小，振荡一会 后趋向稳定的过程，所以显然必须用瞬态分析才能获得正确的结果。图中绘制的是3s瞬时状态下的浓度分布和流速流线。

再如，电热耦合问题中，如果我们需要根据温度分布调整工作电流，也是非常类似的边界非线性问题。