割线法求解过程_浸入边界法求解流固耦合问题

作者：王敏

浸入边界法是近些年发展起来求解流固耦合问题的新算法[1-5]，其核心思想是将固体区域简化为固体边界，将固体与流体之间的相互作用以分布力和分布热的方式加以考虑，从而完成流体与固体的流固耦合。

具体来说，采用欧拉节点描述流体的运动，采用拉格朗日节点描述固体边界的运动，两种节点之间通过狄拉克函数传递数据，如下图所示。

浸入边界法网格系统示意图

浸入边界法的主要优点包括：

1. 整个计算过程都在结构化网格上进行，不需要依照固体形状生成贴体网格或非结构化网格。

2. 当固体处于运动状态或固体在流体作用下存在变形时，采用浸入边界法计算不需要在各个时层上重新生成计算网格，这将显著降低计算量，明显提高计算效率。

到目前为止，浸入边界法已经广泛应用于血液流动、动物游动、圆柱和翼型绕流等领域。

包含浸入边界法的自然对流控制方程为：

其中，f和q分别为欧拉节点