[srivatta]

Hello,
I am trying to write out stream function values on a 2-d surface. Once the case is done running, I invoke the streamFunction utility. I make sure that the streamFunction file exists in each time directory.
Then I invoke the sample utility. The sampleDict is given here. (I got rid of the headers here)

surfaceFormat vtk;
interpolationScheme cellPoint;
surfaces
(
testplane
{
type plane;
basePoint (0 0 0.001);
normalVector (0 0 1);
}
);
fields (

streamFunction

);

After the sample utility is run, the post-processing/surface/time folders show up, but there are no files inside the time directories. Is the "streamFunction" field recognized by the sample utility? If I replace streamFunction with U, p etc, they work.

Thanks,
Srivathsan

[alexeym]

Hi,

streamFunction utility creates pointScalarField (while p is volScalarField). While sample utility can't sample these kind of fields

Once I've made streamFunctionV (https://bitbucket.org/mrklein/streamfunctionv) which first creates pointScalarField and then interpolates point values into cell values. Basically it's streamFunction with the following addition:

            streamFunctionV = dimensionedScalar("zero", phi.dimensions(), 0.0);
            weight = 0.0;
            const labelListList& cellPoints = streamFunctionV.mesh().cellPoints();
            forAll(streamFunctionV, cellI)
            {
                const labelList& curCellPoints = cellPoints[cellI];
                forAll(curCellPoints, cellPointI)
                {
                    streamFunctionV[cellI] += 
                        streamFunction[curCellPoints[cellPointI]];
                    weight[cellI] += 1.0;
                }
            }

            streamFunctionV /= (weight + 1.0);
            streamFunctionV.boundaryField() = 0.0;
            streamFunctionV.write();

[srivatta]

Thank you for your quick reply. I will try this and let you know.
Sri

[srivatta]

This idea works the way I want it to.
Thank you very much.