pisoFoam不收敛

[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

[李东岳]

pisoIBMFoam没有pisoFoam稳健

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

[Samuel-Tu]

@东岳 好的，我试试