CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    pisoFoam不收敛

    OpenFOAM
    2
    3
    554
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • S
      Samuel-Tu 最后由 编辑

      各位老师,我修改了下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
      
      
      1 条回复 最后回复 回复 引用
      • 李东岳
        李东岳 管理员 最后由 编辑

        pisoIBMFoam没有pisoFoam稳健

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

        CFD高性能服务器 http://dyfluid.com/servers.html
        2023年,线下CFD课,预热一下 http://dyfluid.com/class.html

        S 1 条回复 最后回复 回复 引用
        • S
          Samuel-Tu @李东岳 最后由 编辑

          @东岳 好的,我试试

          1 条回复 最后回复 回复 引用
          • First post
            Last post