pisoFoam不收敛



  • 各位老师,我修改了下pisoFoam程序,跑一个带湍流的模型,CFL数一直没有超过0.5,算到0.00014s时就报错了,而且U , p ,epsilon都计算了,不知道是不是算到k的时候不收敛了。请大家帮忙看看,计算过程中的信息如下:

    Create time
    
    Create mesh for time = 0
    
    
    PISO: Operating solver in PISO mode
    
    Reading field p
    
    Reading field U
    
    Reading field IBForce
    
    Reading/calculating face flux field phi
    
    Selecting incompressible transport model Newtonian
    Selecting turbulence model type RAS
    Selecting RAS turbulence model kEpsilon
    RAS
    {
        RASModel        kEpsilon;
        turbulence      on;
        printCoeffs     on;
        Cmu             0.09;
        C1              1.44;
        C2              1.92;
        C3              0;
        sigmak          1;
        sigmaEps        1.3;
    }
    
    No MRF models present
    
    No finite volume options present
    
    Starting time loop
    
    Time = 2e-05
    
    Courant Number mean: 1.69015e-05 max: 0.00679925
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 2.63502e-08, No Iterations 8
    smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 3.55693e-08, No Iterations 9
    DICPCG:  Solving for p, Initial residual = 1, Final residual = 0.000986401, No Iterations 843
    time step continuity errors : sum local = 3.33434e-08, global = -7.61897e-12, cumulative = -7.61897e-12
    DICPCG:  Solving for p, Initial residual = 0.000781366, Final residual = 9.93475e-07, No Iterations 795
    time step continuity errors : sum local = 2.74641e-08, global = 1.4937e-10, cumulative = 1.41751e-10
    DICPCG:  Solving for p, Initial residual = 8.01697e-05, Final residual = 9.97195e-08, No Iterations 784
    time step continuity errors : sum local = 2.3913e-09, global = -1.70151e-11, cumulative = 1.24736e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.0208534, Final residual = 8.99048e-08, No Iterations 7
    smoothSolver:  Solving for k, Initial residual = 1, Final residual = 7.82183e-08, No Iterations 12
    ExecutionTime = 16.87 s  ClockTime = 17 s
    
    Time = 4e-05
    
    Courant Number mean: 0.0328074 max: 0.504716
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.940406, Final residual = 9.91358e-08, No Iterations 81
    smoothSolver:  Solving for Uy, Initial residual = 0.347834, Final residual = 9.54873e-08, No Iterations 110
    DICPCG:  Solving for p, Initial residual = 0.00165164, Final residual = 1.6306e-06, No Iterations 848
    time step continuity errors : sum local = 3.83469e-08, global = -9.46615e-11, cumulative = 3.00746e-11
    DICPCG:  Solving for p, Initial residual = 0.0055641, Final residual = 5.49123e-06, No Iterations 802
    time step continuity errors : sum local = 2.57625e-08, global = -1.17822e-10, cumulative = -8.77479e-11
    DICPCG:  Solving for p, Initial residual = 0.0079625, Final residual = 9.61041e-08, No Iterations 918
    time step continuity errors : sum local = 3.00268e-11, global = 1.8938e-13, cumulative = -8.75585e-11
    smoothSolver:  Solving for epsilon, Initial residual = 0.210254, Final residual = 9.1651e-08, No Iterations 77
    smoothSolver:  Solving for k, Initial residual = 0.474029, Final residual = 9.53898e-08, No Iterations 182
    ExecutionTime = 34.73 s  ClockTime = 35 s
    
    Time = 6e-05
    
    Courant Number mean: 0.0327892 max: 0.506609
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0633166, Final residual = 9.62574e-08, No Iterations 101
    smoothSolver:  Solving for Uy, Initial residual = 0.0486054, Final residual = 9.55871e-08, No Iterations 94
    DICPCG:  Solving for p, Initial residual = 0.112246, Final residual = 0.000110043, No Iterations 795
    time step continuity errors : sum local = 1.10285e-08, global = -7.3844e-11, cumulative = -1.61402e-10
    DICPCG:  Solving for p, Initial residual = 0.125497, Final residual = 0.000125075, No Iterations 769
    time step continuity errors : sum local = 6.75731e-09, global = -7.69153e-11, cumulative = -2.38318e-10
    DICPCG:  Solving for p, Initial residual = 0.0423941, Final residual = 9.79415e-08, No Iterations 931
    time step continuity errors : sum local = 2.29889e-12, global = 2.2377e-14, cumulative = -2.38295e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.0961413, Final residual = 9.11179e-08, No Iterations 91
    smoothSolver:  Solving for k, Initial residual = 0.0817304, Final residual = 9.26786e-08, No Iterations 126
    ExecutionTime = 51.84 s  ClockTime = 53 s
    
    Time = 8e-05
    
    Courant Number mean: 0.0328029 max: 0.355755
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0213574, Final residual = 9.90208e-08, No Iterations 55
    smoothSolver:  Solving for Uy, Initial residual = 0.0222823, Final residual = 9.94334e-08, No Iterations 61
    DICPCG:  Solving for p, Initial residual = 0.514316, Final residual = 0.000503296, No Iterations 711
    time step continuity errors : sum local = 2.01553e-08, global = -7.32031e-11, cumulative = -3.11499e-10
    DICPCG:  Solving for p, Initial residual = 0.415154, Final residual = 0.0004039, No Iterations 723
    time step continuity errors : sum local = 7.62055e-09, global = -6.56086e-11, cumulative = -3.77107e-10
    DICPCG:  Solving for p, Initial residual = 0.0798351, Final residual = 9.7236e-08, No Iterations 930
    time step continuity errors : sum local = 1.49244e-12, global = -4.25276e-15, cumulative = -3.77111e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.00348876, Final residual = 8.96214e-08, No Iterations 23
    smoothSolver:  Solving for k, Initial residual = 0.00931361, Final residual = 9.23941e-08, No Iterations 64
    ExecutionTime = 67.81 s  ClockTime = 69 s
    
    Time = 0.0001
    
    Courant Number mean: 0.0328098 max: 0.346648
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0089072, Final residual = 9.20369e-08, No Iterations 73
    smoothSolver:  Solving for Uy, Initial residual = 0.0109404, Final residual = 9.92794e-08, No Iterations 95
    DICPCG:  Solving for p, Initial residual = 0.531085, Final residual = 0.000523354, No Iterations 780
    time step continuity errors : sum local = 8.3745e-09, global = 1.16496e-10, cumulative = -2.60615e-10
    DICPCG:  Solving for p, Initial residual = 0.315397, Final residual = 0.00031492, No Iterations 758
    time step continuity errors : sum local = 6.16328e-09, global = -8.10869e-11, cumulative = -3.41702e-10
    DICPCG:  Solving for p, Initial residual = 0.0714293, Final residual = 9.83643e-08, No Iterations 949
    time step continuity errors : sum local = 1.62371e-12, global = 3.36497e-15, cumulative = -3.41699e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.0028404, Final residual = 9.96149e-08, No Iterations 52
    smoothSolver:  Solving for k, Initial residual = 0.0108154, Final residual = 9.70358e-08, No Iterations 98
    ExecutionTime = 84.14 s  ClockTime = 85 s
    
    Time = 0.00012
    
    Courant Number mean: 0.0328065 max: 0.33996
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.00834367, Final residual = 9.57569e-08, No Iterations 65
    smoothSolver:  Solving for Uy, Initial residual = 0.0164384, Final residual = 9.00942e-08, No Iterations 84
    DICPCG:  Solving for p, Initial residual = 0.741021, Final residual = 0.00073043, No Iterations 776
    time step continuity errors : sum local = 2.01386e-08, global = -3.08812e-10, cumulative = -6.50511e-10
    DICPCG:  Solving for p, Initial residual = 0.340083, Final residual = 0.000331943, No Iterations 752
    time step continuity errors : sum local = 1.55979e-08, global = 1.91446e-10, cumulative = -4.59065e-10
    DICPCG:  Solving for p, Initial residual = 0.0827811, Final residual = 9.85134e-08, No Iterations 962
    time step continuity errors : sum local = 3.73664e-12, global = -5.7697e-15, cumulative = -4.5907e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.00349239, Final residual = 8.94009e-08, No Iterations 46
    smoothSolver:  Solving for k, Initial residual = 0.0421945, Final residual = 8.85794e-08, No Iterations 92
    ExecutionTime = 100.37 s  ClockTime = 101 s
    
    Time = 0.00014
    
    Courant Number mean: 0.0328501 max: 0.333415
    
    Create IB points
    
    smoothSolver:  Solving for Ux, Initial residual = 0.0207804, Final residual = 9.78132e-08, No Iterations 168
    smoothSolver:  Solving for Uy, Initial residual = 0.0496483, Final residual = 9.85341e-08, No Iterations 225
    DICPCG:  Solving for p, Initial residual = 0.791384, Final residual = 0.000787745, No Iterations 800
    time step continuity errors : sum local = 4.27946e-08, global = 2.6387e-10, cumulative = -1.95201e-10
    DICPCG:  Solving for p, Initial residual = 0.325651, Final residual = 0.000320259, No Iterations 785
    time step continuity errors : sum local = 4.60358e-08, global = -4.45217e-10, cumulative = -6.40418e-10
    DICPCG:  Solving for p, Initial residual = 0.0700157, Final residual = 9.82096e-08, No Iterations 981
    time step continuity errors : sum local = 1.0812e-11, global = -8.00759e-14, cumulative = -6.40498e-10
    smoothSolver:  Solving for epsilon, Initial residual = 0.0447545, Final residual = 9.80666e-08, No Iterations 170
    #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::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, long) at ??:?
    #4  Foam::GaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, long) const at ??:?
    #5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
    #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
    #7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
    #8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
    #9  Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<double> > const&) at ??:?
    #10  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() at ??:?
    #11  ? in ~/OpenFOAM/saumuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoam
    #12  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
    #13  ? in ~/OpenFOAM/saumuel-v1806/platforms/linux64GccDPInt64Opt/bin/pisoIBMFoam
    
    


  • pisoIBMFoam没有pisoFoam稳健

    你试试用PBiCGStab求解器算k和epsilon能不能收敛



  • @东岳 好的,我试试


Log in to reply