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并无太大区别。
    所以我很好奇:

    1. 以上求解器是否求解了d(rho)/dt+div (rhoU)=0?
    2. 如果求解了,其代码段在何处?
    3. 如果有Foamer愿意指点下关于不可压两相流的连续性方程的一些基础问题,我也很感兴趣。

    希望大家一起讨论,希望大神指导赐教



  • 连续性方程 -> 压力方程,因此压力方程就是你说的这个方程


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2