总是算到一定时间步,计算停止报错,有哪位大神知道这样的报错信息是问题出在哪吗?是发散的问题吗



  • kEqnCoeffs
    {
        Ce              1.048;
        Ck              0.094;
    }
    
    
    Reading g
    
    Reading hRef
    Calculating field g.h
    
    No MRF models present
    
    No finite volume options present
    
    GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 8.8753743e-06, No Iterations 8
    GAMGPCG:  Solving for pcorr, Initial residual = 0.024255991, Final residual = 2.1939681e-06, No Iterations 3
    time step continuity errors : sum local = 4.7008361e-15, global = -6.386216e-17, cumulative = -6.386216e-17
    Courant Number mean: 0.001130128 max: 63.68082
    
    Starting time loop
    
    Courant Number mean: 0.001130128 max: 63.68082
    Interface Courant Number mean: 7.3912381e-05 max: 9.8596523
    Time = 0.0021
    
    PIMPLE: iteration 1
    MULES: Solving for alpha.water
    Phase-1 volume fraction = 0.50002593  Min(alpha.water) = -2.5864875e-05  Max(alpha.water) = 1.0000001
    MULES: Solving for alpha.water
    Phase-1 volume fraction = 0.50002492  Min(alpha.water) = -0.010606339  Max(alpha.water) = 1.0000001
    MULES: Solving for alpha.water
    Phase-1 volume fraction = 0.50002392  Min(alpha.water) = -0.0074242149  Max(alpha.water) = 1.0000001
    MULES: Solving for alpha.water
    Phase-1 volume fraction = 0.50002292  Min(alpha.water) = -0.0078825052  Max(alpha.water) = 1.0000001
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::sigFpe::sigHandler(int) at ??:?
    #2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
    #3  Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:?
    #4  Foam::DICSmoother::DICSmoother(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:?
    #5  Foam::lduMatrix::smoother::addsymMatrixConstructorToTable<Foam::DICSmoother>::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&) at ??:?
    #6  Foam::lduMatrix::smoother::New(Foam::word const&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::dictionary const&) at ??:?
    #7  Foam::GAMGSolver::initVcycle(Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::lduMatrix::smoother>&, Foam::Field<double>&, Foam::Field<double>&) const at ??:?
    #8  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
    #9  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
    #10  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
    #11  ? at ??:?
    #12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #13  ? at ??:?
    浮点数例外 (核心已转储)
    

    求助,跪谢!


  • 自由表面模型副教授 OpenFOAM教授

    @regina #3 Foam::DICPreconditioner::calcReciprocalD(Foam::Field<double>&, Foam::lduMatrix const&) at ??:?应该是GAMG代数方程求解器的DIC预处理器遇到值为0的主对角元素了,除0错误。需要检查一下算例。


Log in to reply