compressibleInterFoam中植入drift flux的过程是否正确

[Lubing]

最近在基于compressibleInterFoam添加漂移通量，想请大佬们帮忙看看过程上是否存在明显的问题：组分输运方程如下：

参照compressibleInterFoam中的处理方法，可以整理为：

两式相加得到div v_m:

在第三个方程中添加±α1 div v_m，相对速度的影响拆分成两部分处理，一部分作为源项参照dgdt的处理方法，一部分作为通量处理：

按照这样的处理方法是否存在明显的问题呢？有没有更好的处理方案呢~ 请各位大佬提提意见~

[Lubing]

这里α1表示分散相，α2表示连续项

[李东岳]

vm v1m为什么不写在一起

[Lubing]

可压的场景下，vm和v1m前面的系数不同，我计划参照compressibleInterFoam中的方程的形式来处理组分输运方程，组分输运方程最终的输运形式如下：

完成这个方程的处理之后，按照compressibleInterFoam中的方法，能够得到alpha1Phi和alpha2Phi，刚好能够用于pEqn中的参数构造。想请问下这样是否存在问题~

[Lubing]

采用这样处理的原因而不是直接使用$\mathrm{\partial }(\alpha \ast \rho )/\mathrm{\partial }t$是想直接套用compressibleInterFoam中的MULES方法，保证alpha的有界性

[李东岳]

我看你第一个方程前面都是alpha1rho1

[Lubing]

是的，所以我参照着compressibleInterFoam中的方法对方程做了调整以使最终形式编程\partial alpha/ \partial t+ alpha\div vm的形式

[李东岳]

参照着compressibleInterFoam中的方法

没看明白，compressibleInterFoam没有飘逸对流项，你要用MULES处理有界，你要把你的对流项都用MULES处理。我看你第一个图把对流项放在源项那面去了。所以为什么不把vm v1m放在一起形成一个对流项用MULES处理。

[Lubing]

我大概理解李老师的意思了，我对OpenFOAM中的MULES的原理理解的不透彻，李老师的意思应该是直接基于第一个方程： 直接求解此方程，漂移对流项放到前面去。
我这里处理的时候是想同时使用这两个方程：
最终构造成
，把漂移的影响分开到了源项和对流项中，需要同时修改alphaEqn和alphaSuSp。您觉得这样处理有必要么~
73348540-9b27-449d-9432-26ec0c6fe88b-alphaEqn.H d6ed1b9c-4078-46e3-b212-f7ed1662c6c0-alphaSuSp.H

[李东岳]

你最后一个图，里面还是出现了2个对流项，右边还有个密度梯度（错的），还不如划到一起更简单。

你这个属于求解不同的控制方程，都是描述一套物理，后面的方法更复杂，我不觉得会更好。不过你的PDE越复杂，潜在的所需要的数值方法越多。

[Lubing]

了解了，谢谢李老师~为什么右边的密度梯度是错的呢~

[李东岳]

密度梯度是3个分量，你的alpha方程是1个

[Lubing]

了解了，密度梯度那里是个点乘，v1m·grad(rho1)，在这种情况下是不是也有可能出现有激波时的震荡啊，因为直接使用了密度的梯度~

[李东岳]

那就不好说了 :-)

[Lubing]

我先用一维气液激波管试试看看会不会有显著的震荡，谢谢李老师详细的解答~