CFD Online Discussion Forums

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)

turbulentwakes December 4, 2011 15:05

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.

sahas February 24, 2012 02:22

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;
    }


sahas March 21, 2012 04:02

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


All times are GMT -4. The time now is 00:26.