大密度比两相流发散



  • 最近在用interThermalPhaseChangeFoam算气泡冷凝的验证算例,利用1.3bar水和水蒸气的参数进行模拟,不过总是发散,个人认为是密度比太大了?有没有人做过这方面的研究啊?


  • 网格教授 OpenFOAM教授 管理员

    interThermalPhaseChangeFoam这个求解器我还没看到过…



  • @东岳 https://www.sciencedirect.com/science/article/pii/S2352711016300309 这个不是官方的,个人开发并开源的


  • 网格教授 OpenFOAM教授 管理员

    还没看提前插一句,softwareX这个期刊citeScore很高啊,10点多,但是没有影响因子。关注这个期刊很久了。不知道业内怎么看。好像知道这个期刊的很少。

    一会我看看全文。



  • @东岳 我对这个期刊也不熟。只是在CFDOnline上看到了这个solver。拿来用一下,结果老是发散。

    
    Time = 0.0001
    
    MULES: Solving for alpha1
    Liquid phase volume fraction = 0.999782  Min(alpha1) = 0  Max(alpha1) = 1
    GAMG:  Solving for p_rgh, Initial residual = 1, Final residual = 9.17326e-07, No Iterations 32
    GAMG:  Solving for p_rgh, Initial residual = 4.56128e-07, Final residual = 4.56128e-07, No Iterations 0
    time step continuity errors : sum local = 1.15564e-10, global = 4.67937e-12, cumulative = 4.67937e-12
    GAMG:  Solving for p_rgh, Initial residual = 0.000258322, Final residual = 7.79975e-07, No Iterations 4
    GAMG:  Solving for p_rgh, Initial residual = 7.79882e-07, Final residual = 7.79882e-07, No Iterations 0
    time step continuity errors : sum local = 1.96234e-10, global = 4.41825e-12, cumulative = 9.09762e-12
    GAMG:  Solving for p_rgh, Initial residual = 3.261e-06, Final residual = 6.33275e-07, No Iterations 1
    GAMGPCG:  Solving for p_rgh, Initial residual = 6.33275e-07, Final residual = 9.32048e-09, No Iterations 5
    time step continuity errors : sum local = 3.247e-11, global = -1.80463e-12, cumulative = 7.29299e-12
    ****
    ****Pressure range: 313.455 Pa
    ****Max velocity: 0.168515 m/s
    ****Phase change energy: 0 W
    ****Volume change: 0 m^3/s
    DILUPBiCG:  Solving for H, Initial residual = 0.00100367, Final residual = 7.26908e-09, No Iterations 9
    DILUPBiCG:  Solving for H, Initial residual = 0.000932397, Final residual = 7.34408e-09, No Iterations 9
    DILUPBiCG:  Solving for H, Initial residual = 0.000911756, Final residual = 7.5305e-09, No Iterations 9
    ExecutionTime = 15.77 s  ClockTime = 17 s
    
    Courant Number mean: 2.65292e-05 max: 0.283422
    Interface Courant Number mean: 0 max: 0
    Fourier number mean: 0.00271598 max: 1.81326
    Time = 0.0002
    
    MULES: Solving for alpha1
    Liquid phase volume fraction = 0.999782  Min(alpha1) = -1.24343e-09  Max(alpha1) = 1
    GAMG:  Solving for p_rgh, Initial residual = 0.158707, Final residual = 9.59991e-07, No Iterations 33
    GAMG:  Solving for p_rgh, Initial residual = 9.90688e-07, Final residual = 9.90688e-07, No Iterations 0
    time step continuity errors : sum local = 2.63908e-08, global = -1.91596e-08, cumulative = -1.91524e-08
    GAMG:  Solving for p_rgh, Initial residual = 0.00214601, Final residual = 9.69321e-07, No Iterations 12
    GAMG:  Solving for p_rgh, Initial residual = 9.70692e-07, Final residual = 9.70692e-07, No Iterations 0
    time step continuity errors : sum local = 2.63875e-08, global = -1.9196e-08, cumulative = -3.83484e-08
    GAMG:  Solving for p_rgh, Initial residual = 0.000147321, Final residual = 6.99269e-07, No Iterations 4
    GAMGPCG:  Solving for p_rgh, Initial residual = 6.99262e-07, Final residual = 4.28891e-09, No Iterations 6
    time step continuity errors : sum local = 2.62568e-08, global = -1.91768e-08, cumulative = -5.75252e-08
    ****
    ****Pressure range: 314.046 Pa
    ****Max velocity: 0.158209 m/s
    ****Phase change energy: -0.0533324 W
    ****Volume change: -4.60434e-10 m^3/s
    DILUPBiCG:  Solving for H, Initial residual = 0.0112104, Final residual = 7.17193e-09, No Iterations 11
    DILUPBiCG:  Solving for H, Initial residual = 0.0204525, Final residual = 6.41804e-09, No Iterations 11
    DILUPBiCG:  Solving for H, Initial residual = 0.0172126, Final residual = 4.83197e-09, No Iterations 11
    ExecutionTime = 30.71 s  ClockTime = 32 s
    
    Courant Number mean: 3.50588e-05 max: 0.283714
    Interface Courant Number mean: 1.07527e-05 max: 0.225319
    Fourier number mean: 0.00259948 max: 1.81296
    Time = 0.0003
    
    MULES: Solving for alpha1
    Liquid phase volume fraction = 0.999788  Min(alpha1) = -4.86607e-10  Max(alpha1) = 1
    GAMG:  Solving for p_rgh, Initial residual = 0.248811, Final residual = 8.889e-07, No Iterations 56
    GAMG:  Solving for p_rgh, Initial residual = 4.8953e-08, Final residual = 4.8953e-08, No Iterations 0
    time step continuity errors : sum local = 7.55439e-06, global = -5.52413e-06, cumulative = -5.58166e-06
    GAMG:  Solving for p_rgh, Initial residual = 0.000592748, Final residual = 8.74594e-07, No Iterations 21
    GAMG:  Solving for p_rgh, Initial residual = 8.86839e-07, Final residual = 8.86839e-07, No Iterations 0
    time step continuity errors : sum local = 7.55687e-06, global = -5.52428e-06, cumulative = -1.11059e-05
    GAMG:  Solving for p_rgh, Initial residual = 2.6548e-05, Final residual = 6.06255e-07, No Iterations 2
    GAMGPCG:  Solving for p_rgh, Initial residual = 6.06256e-07, Final residual = 9.35864e-09, No Iterations 5
    time step continuity errors : sum local = 7.55429e-06, global = -5.52414e-06, cumulative = -1.66301e-05
    ****
    ****Pressure range: 488.866 Pa
    ****Max velocity: 0.275259 m/s
    ****Phase change energy: -15.3631 W
    ****Volume change: -1.32635e-07 m^3/s
    DILUPBiCG:  Solving for H, Initial residual = 0.048899, Final residual = 9.83152e-09, No Iterations 11
    DILUPBiCG:  Solving for H, Initial residual = 0.0988781, Final residual = 9.59266e-09, No Iterations 11
    DILUPBiCG:  Solving for H, Initial residual = 0.0302533, Final residual = 4.60914e-09, No Iterations 11
    ExecutionTime = 49.86 s  ClockTime = 51 s
    
    Courant Number mean: 0.00111238 max: 0.451657
    Interface Courant Number mean: 2.23494e-05 max: 0.451657
    Fourier number mean: 0.00249786 max: 1.81736
    Time = 0.0004
    
    MULES: Solving for alpha1
    Liquid phase volume fraction = 0.999796  Min(alpha1) = -0.0745499  Max(alpha1) = 1.0008
    GAMG:  Solving for p_rgh, Initial residual = 0.0207549, Final residual = 9.63291e-07, No Iterations 41
    GAMG:  Solving for p_rgh, Initial residual = 3.93396e-06, Final residual = 8.35811e-07, No Iterations 8
    time step continuity errors : sum local = 9.83598e-06, global = -6.79663e-06, cumulative = -2.34267e-05
    GAMG:  Solving for p_rgh, Initial residual = 0.00157222, Final residual = 8.93124e-07, No Iterations 24
    GAMG:  Solving for p_rgh, Initial residual = 9.05354e-07, Final residual = 9.05354e-07, No Iterations 0
    time step continuity errors : sum local = 9.83602e-06, global = -6.79665e-06, cumulative = -3.02234e-05
    GAMG:  Solving for p_rgh, Initial residual = 0.000103908, Final residual = 5.83835e-07, No Iterations 3
    GAMGPCG:  Solving for p_rgh, Initial residual = 5.83854e-07, Final residual = 7.0116e-09, No Iterations 5
    time step continuity errors : sum local = 9.83541e-06, global = -6.79658e-06, cumulative = -3.70199e-05
    ****
    ****Pressure range: 350.415 Pa
    ****Max velocity: 0.218803 m/s
    ****Phase change energy: -18.9019 W
    ****Volume change: -1.63186e-07 m^3/s
    DILUPBiCG:  Solving for H, Initial residual = 0.1639, Final residual = 9.47499e-09, No Iterations 12
    DILUPBiCG:  Solving for H, Initial residual = 0.272875, Final residual = 5.27952e-09, No Iterations 12
    DILUPBiCG:  Solving for H, Initial residual = 0.214969, Final residual = 6.45167e-09, No Iterations 12
    ExecutionTime = 68.7 s  ClockTime = 70 s
    
    Courant Number mean: 0.00134918 max: 0.480577
    Interface Courant Number mean: 3.19172e-05 max: 0.404985
    Fourier number mean: 0.00247607 max: 1.8178
    Time = 0.0005
    
    MULES: Solving for alpha1
    Liquid phase volume fraction = 0.999804  Min(alpha1) = -0.251971  Max(alpha1) = 1.04115
    [4] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    [4] #1  Foam::sigFpe::sigHandler(int) at ??:?
    [4] #2   in "/lib64/libc.so.6"
    [4] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
    [4] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
    [4] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
    [4] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
    [4] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
    [4] #8  
    [4]  at ??:?
    [4] #9  __libc_start_main in "/lib64/libc.so.6"
    [4] #10  
    [4]  at /usr/src/packages/BUILD/glibc-2.11.3/csu/../sysdeps/x86_64/elf/start.S:116
    APPLICATION TERMINATED WITH THE EXIT STRING: Floating point exception (signal 8)