Coupled level set-VOF方法

Shihang Chen

@Zhujh 您好，也采用了与您基本相同的CLSVOF方法进行计算，然而对于毛细张力主导的问题，这个方法表现出了更强的虚假流动（寄生流动）的问题，如无重力状态下的水滴。请问您遇到过类似的问题么？同时我也参考fluent里面的两种抑制虚假流动方法（密度和H函数），但是收效甚微，想问下您有什么建议吗

EricLiu

@Zhujh @Shihang-Chen 两位大佬好，请问你们用到的CLSVOF 方法有相关的代码资源吗，我最近也打算结合CLS和VOF方法，但是没有找到合适的切入点，openfoam这块属于刚开始学

Shihang Chen

@EricLiu 不是大佬。这个讨论顶部文章附带一些代码，或者直接github搜索clsvof也能有很多收获。

EricLiu

@Shihang-Chen 十分感谢，我去查一下

capillaryFix

@Shihang-Chen 关于spurious currents/velocities，我们此前在OpenFOAM里植入了一个非常简单但是可行的方法，代码量与植入难度不大，具体请参考：https://www.sciencedirect.com/science/article/pii/S1877750323002557

Shihang Chen

@capillaryFix 太好了，感谢感谢，我会好好看看。我目前一直使用的消除寄生流的方法是基于sharpen surface model来的，这个效果也是很好。

李东岳

@capillaryFix 这个流弊，老铁

落花风

@capillaryFix ，您好。我仔细研究了您的文章，有两点疑惑，向您请教：

最后一步是关于capillary flux的过滤是基于原始的$F_{st,f}$，过滤结束的$\phi^{filt}_{cf}$是放在这里吗？

solve
        (
            UEqn
         ==
            fvc::reconstruct
            (
                (
                   - ghf*fvc::snGrad(rho)
                    - fvc::snGrad(p_rgh)
                ) * mesh.magSf()+phi_filt_cf
            ),
            mesh.solver(U.select(finalIter))
        );

和

surfaceScalarField phig
    (
        (
          - ghf*fvc::snGrad(rho)
        )*rAUf*mesh.magSf()+ phi_filt_cf*rAUf
    );

这步$F^{f}{st,f} = F{st,f} − F^{filt}_{st,f}$不知道怎么植入进方程，是不是关于capillary flux的过滤应该基于这步过滤的结果

谢谢、谢谢

capillaryFix

@落花风 1中你的理解是对的，关于2，需要首先根据Eqn. 30,计算出p_{st}，有了p_{st}，就能求解Eqn. 34了呀。祝好运！

落花风

@capillaryFix 是的，我求解了Eqn.30 -> Eqn.34 ->Eqn.33，但是Eqn.33的结果用在哪里了？

CFD中文网

Coupled level set-VOF方法