[bgoeppner]

Hi there,

that's maybe a very simple question, but where can I define which fields OpenFOAM will write? I solved a case with sonicFoam and wanted to have a look at the density, but there is no rho in the time-directories...

Can I make OpenFOAM to calculate and write rho without re-running the hole job?

BTW: Am I right that the pressure is the dynamicPressure? Because then I could calculate my desity from velocity...

Thanks a lot,
Ben

[CedricVH]

The field p is the static pressure (in OpenFOAM 1.6). However, you can calculate the total pressue (static + dynamic pressure) by running the command ptot on the solved case. From there on, you can use foamCalc to calculate other parameters.

[bgoeppner]

Thanks a lot. I tried that but as my p-Field is in Pa (kg/ms²) not in m²/s² ptot wants to have a rho-file... And that is exactly what I'm not having ;-)

[aruv]

Hi Benedikt,
I am facing with the same problem with sonicFoam.. how to calculate ptot ... have u found any solution to the problem ???

[nakul]

Hi,

You may change definition of rho in createFields.H file of sonicFoam locaed in "$applications/solvers".

You may define rho the way U is defined and the recompile your solver. In this way sonicFoam will output rho during runtime.

[aruv]

Hi,
I made the changes in the createFields.H and recompiled the solver, finally it works. thanks for the suggestion

Using Sampledict one can calculate p, T, U etc . do u have any idea how to calculate variables other than the standard ones for Example (isentropic Mach Number) around a compressor profile

thank u

[nakul]

Hi,

To get local Mach No. in the domain you may run the command "Mach".

Also refer to the utilities section of the User Guide to find out what all you can get as data from OF.

[bgoeppner]

Hi aruv,

sorry for my late answer...

You can also add the following to the controlDict file:

functions
{
   rhofunc
   {
      type                 writeRegisteredObject;
      functionObjectLibs   ("libIOFunctionObjects.so");
      outputControl        outputTime;
      outputInterval       1;
      objectNames
      (
         "rho"
         "psi"
      );
   }
}

This will work with every solver and gives you all files you want to have... (In my case this is rho and psi)

Have fun!

[farhhad]

I was just reading your posts and leant lots. There is one thing that seems to be forgotten and It is the dimensions of p and p_rgh in OF 1.7 and OF 1.6. The dimension of p is m^2/s^2 which is the square of velocity dimension. Basically, for p OF skips the density in its calculations. On the other hand, the dimension of p_rgh is kg/ms^2 meaning that density is included. Therefore, if p is the addition of dynamic and static pressure, the actual pressure has yet to be obtained multiplying the density to p.