mixutre模型连续性方程
-
最近在做关于mixutre模型的二次开发。遇到了一个疑惑。是关于连续性方程
d(rho)/dt+div (rhoU)=0 (http://dyfluid.com/twoLiquidMixingFoam.html中的方程6)
的求解代码。
例如，在interFoam中，不需要求解此方程。而是求解alpha方程，和不可压缩形式的连续性方程
div(U)=0。事实上，我觉得这个方程也不是严格成立的。因为他对纯水和纯气的区域当然成立。但是在液面的两相混合处，我觉得理论上来说，描述方程并不可如此简化。但是毕竟液面范围很小，所以这种近似可以接受。
但是对于两相可以混合的twoLiquidMixingFoam和driftFluxFoam而言，我觉得上述方程d(rho)/dt+div (rhoU)=0则必须求解。但是我对比了下这三个求解器的代码，感觉其中的UEqn和pEqn并无太大区别。
所以我很好奇：- 以上求解器是否求解了d(rho)/dt+div (rhoU)=0？
- 如果求解了，其代码段在何处？
- 如果有Foamer愿意指点下关于不可压两相流的连续性方程的一些基础问题，我也很感兴趣。
希望大家一起讨论，希望大神指导赐教
