system/fvSchemes: ok but epsilon blows slowly.

FoamFile{version 2.0; format ascii; class dictionary; location "system"; object fvSchemes;}

snGradSchemes {
    default         limited 0.7;  // limited 0.7; for bad meshes;
                                // corrected; for good meshes;
                                // must fit to laplacianSchemes
    }

laplacianSchemes {default none;
   //laplacian(nuEff,U) Gauss linear corrected;          // for good meshes
   laplacian(nuEff,U) Gauss linear limited 0.7;      // for bad meshes

    //laplacian((1|A(U)),p) Gauss linear corrected;      // for good meshes
    laplacian((1|A(U)),p) Gauss linear limited 0.7;  // for bad meshes
//  for turbulent flow only:
//    laplacian(DkEff,k) Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
    laplacian(DREff,R) Gauss linear corrected;
//    laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; // SA only
}

fluxRequired {
    default         no;
    p              ;
    phi             ;
    }

gradSchemes { // Gradient calculation.  Gauss linear first, cell Limited Gauss linear 1; for better stability with poor meshes.
    default  Gauss linear;  // Gauss linear or leastSquares
    grad(p)         faceLimited leastSquares 0 1;  
    grad(U)         Gauss linear;              
    }

interpolationSchemes {default  linear;     // linear is fine
   U               linear;
//   interpolate(U) reconCentral phi cellLimited leastSquare 1.0; // new tet mesh interpolator. Use with reconCentral in divSchemes
}

divSchemes {default none; div(phi,U)      Gauss upwind;               // stable: Gauss upwind
//div(phi,U)      Gauss linearUpwind Gauss; // faster linearUpwind Gauss
//    div(phi,U)      Gauss reconCentral cellLimited leastSquares 1.0;    // new, tets mesh setting
//    div(phi,k)      Gauss upwind;             // for turbulent flow only
    div(phi,epsilon) Gauss upwind;            // for turbulent flow only
    div(phi,R)      Gauss upwind;             // for turbulent flow only
    div(R)          Gauss linear;             // for turbulent flow only
//    div(phi,nuTilda) Gauss upwind;            // for turbulent SA flow only
div((nuEff*dev(grad(U).T()))) Gauss upwind; // this scheme is not used for
//div((nuEff*dev(grad(U).T()))) Gauss linear; // this scheme is not used for
                                                // steady and laminar flow, but
                                                // it is read by the solver, so
                                                // it must be defined.
//div((nuEff*dev(T(grad(U))))) Gauss linear; // not advised, taken from motorbike case
div((nuEff*dev(T(grad(U))))) Gauss upwind; // not advised, taken from motorbike case
}

ddtSchemes {default steadyState;}    // no time dependence her