我正在寻找易于理解的用于跟踪阶段接口的水平集方法的示例。我想使用 FEM 解决它,因为我的解决方案是基于第二 Fick 定律的 FEM 解决方案。初始化步骤很简单——它只是每个节点与接口的距离。但我对第二步有问题 - 界面速度方程的解:
如果我理解正确,我需要一个微分方程的弱公式,其中 F 是界面的线速度。解决方案会给我一个新的界面位置。这是正确的吗?
我正在寻找易于理解的用于跟踪阶段接口的水平集方法的示例。我想使用 FEM 解决它,因为我的解决方案是基于第二 Fick 定律的 FEM 解决方案。初始化步骤很简单——它只是每个节点与接口的距离。但我对第二步有问题 - 界面速度方程的解:
如果我理解正确,我需要一个微分方程的弱公式,其中 F 是界面的线速度。解决方案会给我一个新的界面位置。这是正确的吗?
TL;DR:是的,总是在以后找到接口。请参阅下面的更多细节:
一个水平集函数,,在二维中,描述了一个函数,其零级集描述了您的界面。准确地说,一个水平集函数划分您的域分为三个不相交的子域:
假如说 总是表示接口,显然接口必须满足:
此外,由于在,这个方程可以写成只涉及法线方向 的速度
这是您的方程式版本. 它有时被称为法线方向的运动。
总而言之,零水平集当您求解水平集方程(以原始形式或正常运动形式)时,始终定位界面。我建议您参考以下两个标准文本以获取更多详细信息和实际实现:
第一个包括关于解决方案的弱公式和 FEM 方法的一章,而第二个主要涉及 FDM 方法。