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

[yfclark]

特征时间尺度，一般用在时间尺度分析上，实际的时间尺度和特征时间尺度之间还差一个系数，这么分析能行吗？

[Albert Lee]

@yfclark 这只是时间尺度的估计，能得到一个基本的数量级，具体问题的话还要看条件。在这种毛细力主导的流动中，该特征时间与毛细波的振荡频率相对应，至少保持在同一数量级。

[yfclark]

fluent的东西不太懂，level-set和vof耦合有专门的算法，要保证level set估算的网格体积分数和vof基本一致，不知道fluent是如何实现的，另外你给出的图的虚假流动太大了，把计算域扩大，边界条件用一对平行的wall，一对平行的slipwall试一下

[李东岳]

我在出差，回去用openfoam测试下给你反馈

[Albert Lee]

@yfclark 谢谢！我也试过了将计算域从原来的1cm x 1cm 扩大到1m x 1m，虚拟流动比原来的略有减小，但是没有量级上的变化。fluent中耦合level-set的方法只是用来计算相界面法向量和曲率，体积分数仍然由vof方法确定。

[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里面