参考文献仍然是An accurate adaptive solver for surface-tension-driven interfacial flows
标准的2D规则网格曲率计算:
- 以欲计算曲率的格子$(i,j,k)$为中心,搞一个$3\times 7$或者$7\times 3$的stencil。具体是哪种取决于之前重建界面的时候MYC算出来的法向朝哪个方向。
- 在这个stencil里面把每一列的VOF值相加,得到高度函数$y=h(x)$或者$x=h(y)$.
- 使用有限差分计算曲率$\kappa=\frac{h”}{(1+h’^2)^{3/2}}$.
这样的问题是$3\tim[……]