[jock]

Hi foamers
I have some questions about divSchemes , when I use MRFSimpleFoam to simulate the centrifugal pump .
If I use below code ,I can simulate successfully.

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

But if I use two order accuracy limitedLinear ,it will be divergent.

divSchemes
{
default none;
div(phi,U) Gauss limitedLinearV 1;
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

My version of OpenFOAM is foam-extend-4.0. And its run log as follows,it only can run one step

Create time

Create mesh for time = 0

Initializing the GGI interpolator between master/shadow patches: GUIDE-OUTLET/IMPELLER-INLET
Initializing the GGI interpolator between master/shadow patches: IMPELLER-OUTLET/VOLUTE-INLET

SIMPLE: no convergence criteria found. Calculations will run for 1400 steps.

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting RAS turbulence model RNGkEpsilon
RNGkEpsilonCoeffs
{
Cmu 0.0845;
C1 1.42;
C2 1.68;
alphak 1.39;
alphaEps 1.39;
eta0 4.38;
beta 0.012;
sigmak 0.71942;
sigmaEps 0.71942;
}

Creating MRF for cell zone IMPELLER. rpm = 2900
Excluded patches: 2(3 4)

Starting time loop

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00436496, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.00451732, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.00450457, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0359147, No Iterations 6
GAMG: Solving for p, Initial residual = 0.0362698, Final residual = 0.000442889, No Iterations 2
time step continuity errors : sum local = 5.86704, global = 0.471927, cumulative = 0.471927
smoothSolver: Solving for epsilon, Initial residual = 0.07982, Final residual = 0.000307204, No Iterations 2
bounding epsilon, min: -1.76856 max: 52.4243 average: 2.25425
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00876813, No Iterations 2
ExecutionTime = 13.49 s ClockTime = 14 s

Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.168279, Final residual = 0.0013948, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.215064, Final residual = 0.00194542, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.185276, Final residual = 0.00162896, No Iterations 2

Please help me,thank you
jack