[Frithjof]

Hello fellow openfoam-users,

I'm trying to modify wallGradU to write the sn-gradient of scalar T at the wall calculated by scalarTransportFoam.

It seems to be very straightforward, it compiles and runs but I want to make sure the results are correct.
Here is what I did:

Im creating an IOObject,
and an object wallGradT:

            volScalarField wallGradT
            (
                IOobject
                (
                    "wallGradT",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                mesh,
                dimensionedScalar
                (
                    "wallGradT",
                    T.dimensions()/dimLength,
                    0
                )
            );

and calculating the surface-normal gradient just as in wallGradU:

wallGradT.boundaryField()[patchi] = -T.boundaryField()[patchi].snGrad();

Is this the right way?
I want to make sure I'm getting correct results..

Thanks in advance!

[boger]

That looks believable. You can find almost the same logic in applications/utilities/postProcessing/wall/wallGradU/wallGradU.C