[faiazk]

Hello All,

I have been using kepsilon turbulence model so far and was able to get it running both with simpleFoam and pisoFoam before. Now, for my particular case I need to use RSM (LRR) turbulence model and I am planning to implement it in simpleFoam. Right now I am trying to implement it in an empty domain. Initially, I was having few errors while implementing LRR which I resolved following instructions from CFD online. But, right now I am having a 'time continuity error' right after going through 6 time steps. I am attaching the error here. Please help me with this.

Thank you everyone.

Time = 10.05

DILUPBiCGStab: Solving for Ux, Initial residual = 1.02274e-13, Final residual = 1.02274e-13, No Iterations 0
DILUPBiCGStab: Solving for Uy, Initial residual = 4.08144e-14, Final residual = 4.08144e-14, No Iterations 0
DILUPBiCGStab: Solving for Uz, Initial residual = 1.24302e-13, Final residual = 1.24302e-13, No Iterations 0
DICPCG: Solving for p, Initial residual = 0.874444, Final residual = 0.00868986, No Iterations 65
DICPCG: Solving for p, Initial residual = 4.34332e-22, Final residual = 4.34332e-22, No Iterations 0
DICPCG: Solving for p, Initial residual = 4.34332e-22, Final residual = 4.34332e-22, No Iterations 0
time step continuity errors : sum local = 0.000316474, global = -4.65325e-05, cumulative = -4.65324e-05
DILUPBiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.0422663, No Iterations 1
bounding epsilon, min: -1.28744e+14 max: 6.48512e+19 average: 4.96778e+15
DILUPBiCGStab: Solving for Rxx, Initial residual = 5.02367e-11, Final residual = 1.37406e-16, No Iterations 1
DILUPBiCGStab: Solving for Rxy, Initial residual = 6.12384e-09, Final residual = 5.08093e-17, No Iterations 1
DILUPBiCGStab: Solving for Rxz, Initial residual = 3.86031e-09, Final residual = 1.21732e-15, No Iterations 1
DILUPBiCGStab: Solving for Ryy, Initial residual = 1.6043e-10, Final residual = 1.38089e-16, No Iterations 1
DILUPBiCGStab: Solving for Ryz, Initial residual = 6.69048e-09, Final residual = 2.86888e-16, No Iterations 1
DILUPBiCGStab: Solving for Rzz, Initial residual = 1.35119e-10, Final residual = 1.46835e-16, No Iterations 1
ExecutionTime = 32.55 s ClockTime = 33 s

fieldAverage fieldAverage1 write:
Calculating averages

Time = 10.06

DILUPBiCGStab: Solving for Ux, Initial residual = 7.05192e-08, Final residual = 6.13968e-09, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 2.11358e-07, Final residual = 4.31941e-09, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 3.85413e-07, Final residual = 9.98109e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00982227, No Iterations 225
DICPCG: Solving for p, Initial residual = 0.000615733, Final residual = 5.6734e-06, No Iterations 230
DICPCG: Solving for p, Initial residual = 1.43324e-05, Final residual = 1.37545e-07, No Iterations 258
time step continuity errors : sum local = 3.2069e+32, global = 8.42778e+30, cumulative = 8.42778e+30
DILUPBiCGStab: Solving for epsilon, Initial residual = 1, Final residual = 0.0319122, No Iterations 1
bounding epsilon, min: -3.67017e+46 max: 2.42397e+47 average: 4.13154e+42
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::lduMatrix const&) at ??:?
#4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam:ILUPreconditioner>::New(Foam ::lduMatrix::solver const&, Foam::dictionary const&) at ??:?
#6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foam::dictionary const&) at ??:?
#7 Foam::PBiCGStab::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<Foam::SymmTensor<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#9 Foam::fvMatrix<Foam::SymmTensor<double> >::solve(Foam::dictionary const&) at ??:?
#10 Foam::fvMatrix<Foam::SymmTensor<double> >::solve() at ??:?
#11 Foam::SolverPerformance<Foam::SymmTensor<double> > Foam::solve<Foam::SymmTensor<double> >(Foam::tmp<Foam::fvMatrix<Foam::SymmTensor<double > > > const&) at ??:?
#12 Foam::RASModels::LRR<Foam::IncompressibleTurbulenc eModel<Foam::transportModel> >::correct() at ??:?
#13 ? at ??:?
#14 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#15 ? at ??:?
Floating point exception

Thanks

Regards,
Faiaz