openfoam的twomixingLiquidfoam求解器改成稳态求解器

[1064168551]

各位老师好，我最近想把openfoam的twomixingLiquidfoam求解器改成稳态求解器，不知道怎样实现最简单？我的初步想法是把求解速度场和压力场的pimple求解器先改成simple求解器，然后把求解相alpha
1的方程中非稳态项（t的导数项）去掉，这样做可以吗？MULES的相关求解部分需要进行修改吗？或者各位老师有什么更简便的建议和方法？

[李东岳]

简单来说不太好弄。主要是MULES那面是按照时间步推进的。目前OpenFOAM的做法是使用局部时间步，就是local time step。还涉及到一个耦合问题。比如单相流是pU耦合。但是可压缩稳态就设计到rho p U耦合。那你多相流就是alpha p U耦合。这个迭代流程可能需要自己设计一下。以上是我的想法。

不过我知道商软好像可以做稳态VOF，稳态EE，这个需要研究下。你可以找找相关的sci，看看怎么处理的。

你有进展也可以保持更新，我跟进一下，挺有意思。

[1064168551]

@李东岳 好的，谢谢老师，我目前主要是想对多组分气体（氧气和氮气）的稳态流动扩散进行求解，想着基于twomixingLiquidfoam简单修改修改，因为我对MULES研究的不多，如果是稳态计算的话，我可否直接在相方程里面添加对流项fvm::div(phi,alpha1)，将求解alpha1的方程写成：
fvm::div(phi,alpha1) - fvm::laplacian
(
volScalarField("Dab", Dab + alphatab*turbulence->nut()),
alpha1
)
从而删去MULES部分呢？