interFoam求解器的液滴的奇怪现象（虚假速度？）

[Albert Lee]

@麦迪文 虚拟流动方面现在一些开源代码如Gerris等都已经做得很好了，貌似level set效果也好一些。时间充裕的话可以看看，但是对于自己的问题又可能不太适用，修改起来挺麻烦。

[麦迪文]

（热毛细对流反向加在自由面后，自由面上速度被扭转回来）

[李东岳]

你这个设置壁面接触角了么

[麦迪文]

我这个没设接触角，不过在openfoam的tutorials
OpenFOAM-5.x/tutorials/multiphase/interFoam/laminar/capillaryRise/
这个case里，已经设置了接触角，然后将液相入口设为壁面条件，入口压力也设为fixedFluxPressure，结果如下

自由面附近出现了较大的速度

[李东岳]

你这个图看起来需要设置接触角？尤其是C图

[麦迪文]

@东岳 下面这个实验的图效果好点

下面是fluent设置接触角 下壁面设置的110度，上壁面设置70度时数据

工作站有点问题，暂时用不了openfoam，过两天贴下of加接触角的结果

[李东岳]

不是水？液体不润湿？

[Albert Lee]

@麦迪文 看了一下您这个模拟也是将热毛细力加到气液界面上了吧？我想请教一下udf里面是先用adjust求出梯度，用udm纪录后再在将热毛细力添加到动量源项里？我是算一个带有热毛细的液体晃动（不做液桥方向），感觉跟这里类似，所以请教一下您

[麦迪文]

这个在三维就还复杂些，在二维，实际简单得多，热毛细力就是在帮助里的温度梯度*（方向扭转90度的相分数梯度）再乘一个系数，这些在帮助里都能查到

[aaaachenjian]

@麦迪文 我做过些液桥震动的，不过用的fluent。你这个考虑了温度可能会不太一样。
但是我感觉，你模拟的这个液桥，体积和实验应该不一样吧，看着模拟里面好像体积大些。还有其他面你设置的什么边界条件呢？

[麦迪文]

@aaaachenjian 模拟是在微重力环境下，上下壁面有个温差，在自由面上加个热毛细力，已经放弃用VOF方法了，像这种变形较小的直接挪网格，这是算法天生问题不太好解决