October 25, 2021, 16:00
|
OF7: strain rate invariant computation in Epsilon Wall Function
|
#1
|
New Member
Join Date: May 2020
Posts: 23
Rep Power: 6
|
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:
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!
Last edited by Bibill; October 26, 2021 at 03:19.
|
|
|