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

[李东岳]

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]

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