interFoam求解器的液滴的奇怪现象（虚假速度？）
-
这个稳态simpleFoam不可能出来，必须是瞬态的，瞬态流场耦合dpm，单向耦合，不考虑颗粒对流场的作用，就能出来三角形的PAS
水平太低，不能理解 “同一个网格点不同的传输速度” 指的是啥，我现在的简化思路有两个
1.就是在连续相的动量方程中去掉曳力项
2.就是用pimlefoam ，外挂两个方程
固相动量方程
固相分数也就是想用欧拉数值模化拉格朗日模型
感觉上可行，还没具体开始操作。。。。这两个想法在fluent里实现特别难（得用巨量的udf），也然而，在openfoam里实现却没有那么的麻烦，这也是openfoam比fluent的魅力所在啊 -
还有就是我觉得颗粒作为杂质也可能会有一些聚并与破碎现象，cfd-PBM改造后可能可以应用在这一领域上，但是需要一些实验数据支撑，实验台我们组已经搭好，就是实验着实难做，而且对openfoam-7中的pbm模型也是一头雾水
-
需加速度主要产生于表面张力呈主要贡献的工况，尤其是使用CSF模型的时候，稍微一点点曲率计算不正确都会导致压力梯度和表面张力失衡，进而需要虚假速度去补偿。楼主后来这个解决了么
看Zaleski对虚假速度的讨论 https://www.jianguoyun.com/p/DRs785wQ9s3ZBhj6yr8C
-
不好意思，回复晚了，没有什么进展，已经用了fluent解决此问题, 如李老师所说，应该是openfoam曲率计算不准导致，当在openfoam把重力设置为0时，表面不再变形，虚假速度也相继大幅度降低基本消失，虽然猜测是这个问题导致，但是由于数学基础不好，所以没有什么改进的措施
上传的case是在openfoam5_x版本，直接运行interFOAM即可（由于研究的是热毛细对流，其特性是使自由面上出现沿自由面的切向速度，刚好与虚假速度方向相同或相反，从而使计算结果误差巨大）
-
附：
在OF5中 interfoam计算曲率在：
OpenFOAM-5.x/applications/solvers/multiphase/multiphaseInterFoam/multiphaseMixture/multiphaseMixture.C中的414行-----530行 -
回复上上个话题，同一个网格点的不同传输速度不会处理所以没有处理。。。。。。在网格内的颗粒速度认为一致。。。。。请问泰勒涡的这个不同算法对比的case的文章发表了吗？能不能发个这个文章的连接学习一下
-
你为什么不生成这样一个网格（红色的线）？旁边黑色的计算域有什么用？
同一个网格点不同的传输速度那个算法太复杂了 请参考：
https://www.jianguoyun.com/p/DUDlftMQ9s3ZBhjeq8EC
https://www.jianguoyun.com/p/DQeqIUMQ9s3ZBhjqq8EC
-
是按照实验来的，周围是空气，上下是俩圆柱
-
之前的数据没找到，又重新画网格算了一下，也有这个问题。。。。但是，fluent把热毛细力的UDF插入进去后，虚假速度影响会减弱，自由面上速度方向会变化，但openfoam用fvoption方法加入热毛细力，没法抑制虚假速度（以上两图是没加热毛细力，纯自然状态的形状）
晚上会把fluent加热毛细力使自由面的速度反向图附上
-
@麦迪文 您好！我想在计算域里面的气液界面加热毛细条件的udf，不知道如何添加。请问您是怎么做的呢？方便给我一份您的热毛细udf参考一下吗？非常感谢您！我的邮箱是lijichengbitao@163.com
-
首先，拿到温度梯度和界面梯度，两个变量做一点数学计算就得到沿界面的热毛细力了，一旦我发文章了会把代码贴文章后面的，数学公式也很简单，UDF不是很麻烦哈
-
@麦迪文 您的热毛细力是加在边界吧？我是要加在计算域里头，所以有点麻烦，热毛细力的公式是
-
我之前用的VOF法，只能加在计算域里，如果用普通方法，就加在边界上
-
@麦迪文 我也是用vof方法，因为还有气相，所以热毛细力只能加在气液界面上。但是微重力条件下虚拟流动很厉害，跟热毛细速度在一个量级上，所以对结果准确性影响比较大
-
是是是，我也因为这个没啥解决办法，好在我导师不怎么管我，自己定方向，研究了一阵子没突破出来就整稀相流去了
-
@Albert-Lee 我想到的一种解决办法是界面变形不大的情况下或许可以用动网格的方法解决
-
@麦迪文 虚拟流动方面现在一些开源代码如Gerris等都已经做得很好了，貌似level set效果也好一些。时间充裕的话可以看看，但是对于自己的问题又可能不太适用，修改起来挺麻烦。
18/38