多孔介质相变发散

hongjiewang

各位老师好，我目前的算例是模拟一个池空间的蒸发现象，下方是多孔介质（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

hongjiewang

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

多孔介质相变 发散

多孔介质相变发散