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

[李东岳]

你为什么不生成这样一个网格（红色的线）？旁边黑色的计算域有什么用？

同一个网格点不同的传输速度那个算法太复杂了 请参考：
https://www.jianguoyun.com/p/DUDlftMQ9s3ZBhjeq8EC
https://www.jianguoyun.com/p/DQeqIUMQ9s3ZBhjqq8EC

[麦迪文]

是按照实验来的，周围是空气，上下是俩圆柱

[李东岳]

Fluent那面结果如何？

[麦迪文]

之前的数据没找到，又重新画网格算了一下，也有这个问题。。。。但是，fluent把热毛细力的UDF插入进去后，虚假速度影响会减弱，自由面上速度方向会变化，但openfoam用fvoption方法加入热毛细力，没法抑制虚假速度（以上两图是没加热毛细力，纯自然状态的形状）

晚上会把fluent加热毛细力使自由面的速度反向图附上

[Albert Lee]

@麦迪文 您好！我想在计算域里面的气液界面加热毛细条件的udf，不知道如何添加。请问您是怎么做的呢？方便给我一份您的热毛细udf参考一下吗？非常感谢您！我的邮箱是lijichengbitao@163.com

[麦迪文]

首先，拿到温度梯度和界面梯度，两个变量做一点数学计算就得到沿界面的热毛细力了，一旦我发文章了会把代码贴文章后面的，数学公式也很简单，UDF不是很麻烦哈

[Albert Lee]

@麦迪文 您的热毛细力是加在边界吧？我是要加在计算域里头，所以有点麻烦，热毛细力的公式是$\mu \nabla u = \sigma_T \nabla T$吧？

[麦迪文]

我之前用的VOF法，只能加在计算域里，如果用普通方法，就加在边界上

[Albert Lee]

@麦迪文 我也是用vof方法，因为还有气相，所以热毛细力只能加在气液界面上。但是微重力条件下虚拟流动很厉害，跟热毛细速度在一个量级上，所以对结果准确性影响比较大

[麦迪文]

是是是，我也因为这个没啥解决办法，好在我导师不怎么管我，自己定方向，研究了一阵子没突破出来就整稀相流去了

[麦迪文]

@Albert-Lee 我想到的一种解决办法是界面变形不大的情况下或许可以用动网格的方法解决

[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度的相分数梯度）再乘一个系数，这些在帮助里都能查到