Fluent VOF 耦合Level-Set方法时的Density correction 抑制虚拟流动

[Albert Lee]

@东岳 真的十分感谢！算例上传上来了，见附件。spuriousCurrentArc4-1.zip

[刘晓刚]

@Albert-Lee 请问您知道fluent中Vof耦合level set如何进行初始化吗？简单液滴撞平壁算例，初始化时要对液滴区域给定level set function函数值吗？ 我用的是heaviside方法。

麻烦看一下这个提问，http://www.cfd-china.com/topic/3500/fluent-vof-和level-set-耦合

[Albert Lee]

@刘晓刚 直接全场给0或者1皆可，一个迭代之后自然会算出相应的levelset值。

[刘晓刚]

@Albert-Lee 那您看我这个level set function算的是不是不对啊？
http://www.cfd-china.com/topic/3500/fluent-vof-和level-set-耦合

[刘晓刚]

正常level set function=0 时不应该为两相界面吗？

[Albert Lee]

@刘晓刚 这个耦合了vof，即使你初始化levelset为0，后面计算会结合vof自动赋值的

[刘晓刚]

再麻烦问一下，您知道我这个计算大概哪出问题了吗？vof 单算和vof+level set的结果不一样啊，明显vof的交界面更加光滑。现在通过level set function=0 也判断不了相界面啊？

[Albert Lee]

@刘晓刚 我回复至你那个问题里面了。单个因素试着来，逐一排查，应该会找到原因的

[anubis]

你现在解决虚拟流动的问题了吗？

[martinliu]

@Albert-Lee 我跟你遇到了相同的问题，你的测试和我的结果也相差不大。目前我也在寻找一个解决方法，我看到的在UDF中获取Level set function 是使用 C_LSF(c,t)宏，其中c = cell index, t = mixture phase thread，原链接：https://www.cfd-online.com/Forums/fluent/144991-access-level-set-function-udf.html

注意：需要额外include "sg_ls.h" 因为他没有包含在udf.h里面