CFD Online Discussion Forums - laminar compressible solver for rhoCentralFoam (OF-2.0.x)

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)

- - laminar compressible solver for rhoCentralFoam (OF-2.0.x) (http://www.cfd-online.com/Forums/openfoam-solving/95051-laminar-compressible-solver-rhocentralfoam-2-0-x.html)

laminar compressible solver for rhoCentralFoam (OF-2.0.x)

Hi,

I've installed OpenFOAM-2.0.1 few days back and started using it. In all the tutorials for rhoCentralFoam, I see mu =0 in thermophysical properties. If I want to change it's value for laminar case, it returns error like this

----------------------------------------------------------------------------------------
--> FOAM FATAL IO ERROR:
keyword e is undefined in dictionary "/home/rajesh/OpenFOAM/rajesh-2.0.1/run/tutorials/compressible/rhoCentralFoam/forwardStep/system/fvSolution::solvers"

file: /home/rajesh/OpenFOAM/rajesh-2.0.1/run/tutorials/compressible/rhoCentralFoam/forwardStep/system/fvSolution::solvers from line 22 to line 38.

From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 461.

FOAM exiting
---------------------------------------------------------------------------------------

I don't understand how to proceed, can somebody help? Also in programmers guide I see mention of a tutorial with non-zero mu for forward step case (sonicFoam case), but I don't see the same in the installed folder. Also, could somebody give some reference where laminar/turbulent cases at high reynolds no. have been solved using rhoCentralFoam.

You should specify solver for e in fvSolution.
For example, instead of

Code:

   U

    {

        solver          smoothSolver;

        smoother        GaussSeidel;

        nSweeps         2;

        tolerance       1e-09;

        relTol          0.01;

    }

in forwardStep tutorial use

Code:

   "(U|e)"

    {

        solver          smoothSolver;

        smoother        GaussSeidel;

        nSweeps         2;

        tolerance       1e-09;

        relTol          0.01;

    }

Error in Prandtl number

I've just found that Prandtl number in rhoCentralFoam is always 1 (independent of what you have set in constant/thermophysicalProperties).
See bug http://www.openfoam.org/mantisbt/view.php?id=475