diff --git a/openfoam2012/freeCylinderFlow/system/fvSchemes b/openfoam2012/freeCylinderFlow/system/fvSchemes index 888aa0933f49188bcab48746b97f4a625a2b417d..5a805ab74e371e4e30df9c6113a4611701498c00 100644 --- a/openfoam2012/freeCylinderFlow/system/fvSchemes +++ b/openfoam2012/freeCylinderFlow/system/fvSchemes @@ -17,7 +17,8 @@ FoamFile ddtSchemes { - default Euler; + /* default Euler; */ + default backward; } gradSchemes @@ -29,7 +30,8 @@ divSchemes { default none; - div(phi,U) Gauss linearUpwind grad(U); + /* div(phi,U) Gauss linearUpwind grad(U); */ + div(phi,U) Gauss limitedLinearV 1; div(phi,T) Gauss limitedLinear 1; turbulence Gauss limitedLinear 1; @@ -37,7 +39,7 @@ divSchemes div(phi,epsilon) $turbulence; div(phi,omega) $turbulence; // error was that this was missing. Just do the same as with the other div(phi,...)... - div((nuEff*dev2(T(grad(U))))) Gauss linear; + div((nuEff*dev2(T(grad(U))))) $turbulence; // to have the same everywhere Gauss linear; } laplacianSchemes