openfoam的twomixingLiquidfoam求解器改成稳态求解器
-
各位老师好,我最近想把openfoam的twomixingLiquidfoam求解器改成稳态求解器,不知道怎样实现最简单?我的初步想法是把求解速度场和压力场的pimple求解器先改成simple求解器,然后把求解相alpha
1的方程中非稳态项(t的导数项)去掉,这样做可以吗?MULES的相关求解部分需要进行修改吗?或者各位老师有什么更简便的建议和方法? -
简单来说不太好弄。主要是MULES那面是按照时间步推进的。目前OpenFOAM的做法是使用局部时间步,就是local time step。还涉及到一个耦合问题。比如单相流是pU耦合。但是可压缩稳态就设计到rho p U耦合。那你多相流就是alpha p U耦合。这个迭代流程可能需要自己设计一下。以上是我的想法。
不过我知道商软好像可以做稳态VOF,稳态EE,这个需要研究下。你可以找找相关的sci,看看怎么处理的。
你有进展也可以保持更新,我跟进一下,挺有意思。
-
@李东岳 好的,谢谢老师,我目前主要是想对多组分气体(氧气和氮气)的稳态流动扩散进行求解,想着基于twomixingLiquidfoam简单修改修改,因为我对MULES研究的不多,如果是稳态计算的话,我可否直接在相方程里面添加对流项fvm::div(phi,alpha1),将求解alpha1的方程写成:
fvm::div(phi,alpha1) - fvm::laplacian
(
volScalarField("Dab", Dab + alphatab*turbulence->nut()),
alpha1
)
从而删去MULES部分呢?