CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Can't output Reynolds stress in pisoFoam (http://www.cfd-online.com/Forums/openfoam/113676-cant-output-reynolds-stress-pisofoam.html)

bhattach February 25, 2013 03:47

Can't output Reynolds stress in pisoFoam
 
I am trying to run a channel flow using the k-epsilon model, with the pisoFoam solver. While the result seems to be converging, and I can see outputs for U, k, epsilon, nut, p, phi, I don't see any output files for R (Reynolds stress). Is there any way I can make the solver output R as well ? I can of course calculate R from k, nut and grad U, but it would be more convenient to just get R outputted directly.

I am attaching the contents of the 0/R file below. Any help will be appreciated.

regards

Amitabh

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volSymmTensorField;
object R;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform (0 0 0 0 0 0);

boundaryField
{
inlet
{
type zeroGradient;
}

outlet
{
type zeroGradient;
}

fixedWalls
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}

frontAndBack
{
type empty;
}
}

// ************************************************** *********************** //

chegdan February 25, 2013 13:06

You could use the R utility that is already there in OpenFOAM. Just go to your case directory and type

Code:

R -help
If you want to add this to pimpleFoam, go check out the R utility in

Code:

$FOAM_APP/utilities/postProcessing/turbulence/R
for inspiration.

bhattach February 25, 2013 13:17

Thanks a lot, Dan, that worked great for me.


All times are GMT -4. The time now is 02:58.