多孔介质相变 发散



  • 各位老师好,我目前的算例是模拟一个池空间的蒸发现象,下方是多孔介质(eps=0.5),多孔介质内充满液体,上方自由空间是气体。整个空间的温度均为饱和温度下的373.15K,下方壁面温度为378.15K,相变模型采用Lee模型。结果在跑了几步之后就出现发散现象。我把算例文件上传在下方,麻烦各位老师看看哪里出了问题~

    Phase-1 volume fraction = -2.39612e+89  Min(alpha.liquid) = -3.11025e+94  Max(alpha.liquid) = 5.09398e+91
    smoothSolver:  Solving for alpha.liquid, Initial residual = 1, Final residual = 2.10491e-91, No Iterations 1
    Phase-1 volume fraction = 1.2399e+51  Min(alpha.liquid) = -1.12053e+52  Max(alpha.liquid) = 2.23433e+56
    MULES: Correcting alpha.liquid
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in /lib/x86_64-linux-gnu/libpthread.so.0
    #3  ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/twoPhaseChangePorousFoam
    #4  ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/twoPhaseChangePorousFoam
    #5  ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/twoPhaseChangePorousFoam
    #6  ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/twoPhaseChangePorousFoam
    #7  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
    #8  ? in ~/OpenFOAM/dyfluid-v2006/platforms/linux64GccDPInt32Opt/bin/twoPhaseChangePorousFoam
    Floating point exception (core dumped)
    

    test2.zip



  • 补充一下~算例用的求解器是自己在原求解器基础上修改的,这个求解器已经可以确定没有语句或者逻辑上的问题了。之前用这个求解器已经正确跑过算例了,结果也是和物理相符合的。只是前两天给虚拟机扩容,操作失败,然后所有的算例文件还没来得及保存。这两天新设的这个算例总是发散,我觉得是边界条件的问题,但是又不太懂,所以还是得麻烦各位老师~


Log in to reply
 


CFD中文网 | 东岳流体学术 | 东岳流体商业 | 吉ICP备20003622号-1