sonicFoam #include[s] rhoEqn.H...
...but the rhoEqn.H cannot be found in the solver directory (version 1.6). I was curious when I was using the solver before, it gave results for everything, but the rho was always calculated as 0. Does anybody know the purpose of this? Is this intentional? Can I incorporate some other rhoEqn.H easily?

It's in the folder OpenFOAM1.6.x/src/finiteVolume/cfdTools/compressible.

Revisiting this thread after some time:
I am running a case at the moment, using sonicFoam. The case is pretty much the same as /home/username/OpenFOAM/OpenFOAM1.6/tutorials/compressible/sonicFoam/ras/prism/ but instead of "prismWall", I have "wall" which is the wall of my domain. While running, in each time step, I get: Code:
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0 My questions are:

Hi,
Any solution for this problem. 
Probably you need to supply an equation of state in your case dir? that is dont you need a thermoPhysicalProperties dict? I cant say for shure as I am looking at engineFoam and rhoPimpleFOAM at the moment, but that seems to be what's happening there. rhoEqn.h is just the mass conservation equation.

Quote:
But still don't know how to get the rho values of the fields, two options: 1. state equation 2. phi=rho*velocity*A 
Solved?
Do you guys solved this issue? Why is rho equation not being calculated?

Hi Maders
Try to add these extra sentence at the end of controlDict, and the rho will output in the calculation folder. Quote:

Worked!
Quote:
that worked! I have already done it by inserting write() functions at runTime, but I found your way more straightforward.. Thanks!:) 
Actually the answer is..
Actually the main answer iin this topic was about the rho calculation.. That seemed not to be doing by the "0" values of iterations and residuals. Here is the aswer, which is very logic by the way: http://www.cfdonline.com/Forums/ope...opisofoam.html
Regards, Lisandro 
Additional note
I would just like to add to this that by adding the rho function I was able to use the R utility to compute Reynolds stress when using the sonicFoam solver. This was very helpful if you are using some other turbulence model at the start of the solution with the intent to switch to LRR model after the flow field has been initialized.

All times are GMT 4. The time now is 19:40. 