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

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

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

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

for inspiration.

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