VOF添加斥力模型

学流体的小明

李老师，想请教一下源项的添加的问题。
比如我现在在有一个repulsion源项，要想植入，应该是有三个地方要添加。

一，速度预测方程，即UEqn.H中的fvMatrix UEqn。但我不清楚这里是否要加上源项。
加的理由是：fvOptions这种人工添加的源项都在。
不加的理由是：写的那个icoFoam解析，以及这里实际上只保留了对流项和粘度项，压力、表面张力和重力都没有加，这个斥力在形式上是和表面张力有关的，就像表面张力那样处理。

fvVectorMatrix UEqn
    (
        fvm::ddt(rho, U) + fvm::div(rhoPhi, U)
      + MRF.DDt(rho, U)
      + turbulence->divDevRhoReff(rho, U)
     ==
        fvOptions(rho, U)
      + repulsion_Vector // add
    );

二，真实速度方程，一般的想法都是加上这个源项，就是最后的速度修正项是

但是在reactingTwoPhaseEulerFoam中，我看到这些源项都是直接加到最后的速度上的，不知道这里是您省略它们的具体形式了，还是这些$M$确实就不除以$ A_{d,p} $

就是这两种方式应该是哪一种？还是两种都可以？

三，上面两种不同方式的速度方程，对应不同的压力泊松方程，也就差个$ A_{d,p} $

李东岳

速度预测方程，即UEqn.H中的fvMatrix UEqn。但我不清楚这里是否要加上源项。

不需要加

源项你直接1）加到phi上面组建压力方程，2）最后的速度更新要加上源项

我那个M都忘记除掉A了，是要除掉，我有遗漏，我周末更新过去
http://dyfluid.com/reactingTwoPhaseEulerFoam.html 这个里面41-46之前写的一直没更新 29之前是我最近更新的没问题

学流体的小明

不需要加

源项你直接1）加到phi上面组建压力方程，2）最后的速度更新要加上源项

这种方法算得的结果就是那些有压力振荡不稳定的。

另一种方法：仅添加到UEqn当中，压力方程不做处理，斥力的贡献依靠UEqn.H参与到后续的计算当中。结果就没有压力振荡了。
目前初步的结果是这样的，后续可能还有其它问题。

非常感谢李老师这么久这么多的指点，帮助太大了。

李东岳

另一种方法：仅添加到UEqn当中，压力方程不做处理，斥力的贡献依靠UEqn.H参与到后续的计算当中。结果就没有压力振荡了。

Interesting！

李东岳

@学流体的小明大佬后来你们多相流直接模拟的工作发出来没

学流体的小明

@李东岳谢谢李老师关心呀。算例模拟完了，文章😂我还在写。
算出来的结果不是很好，气泡槽道流收敛性还比较差，几个参数之间的规律不是很明显，故事不太好讲。

李东岳

不好讲故事是比较困境了，哎..

CFD中文网

VOF添加斥力模型