[Bibill]

Dear all,

I would like to compute both the strain rate and the vorticity invariants in the calculate() function of the Epsilon Wall function in OF7.

Here is the structure of the code:

void Foam::epsilonWallFunctionFvPatchScalarField::calculate
(
    const turbulenceModel& turbModel,
    const List<scalar>& cornerWeights,
    const fvPatch& patch,
    scalarField& G0,
    scalarField& epsilon0
)
{

    const tmp<volScalarField> tk = turbModel.k();
    const volScalarField& k = tk();

    const tmp<volScalarField> te = turbModel.epsilon();
    const volScalarField& epsilon = te();

    const volVectorField& Ut = turbModel.U();

    /* Strain rate invariant */
    volSymmTensorField S = symm(fvc::grad(Ut));

    /* Vorticity invariant */
    volTensorField W = skew(fvc::grad(Ut));

...

}

The code is compiled but I have an MPI_Waitall error when I try to run my test case.

Does anyone have an idea to solve my problem ?

Thanks in advance!