CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

OF7: strain rate invariant computation in Epsilon Wall Function

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   October 25, 2021, 16:00
Default OF7: strain rate invariant computation in Epsilon Wall Function
New Member
Join Date: May 2020
Posts: 22
Rep Power: 4
Bibill is on a distinguished road
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!

Last edited by Bibill; October 26, 2021 at 03:19.
Bibill is offline   Reply With Quote


openfoam-7, programming

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 13 March 7, 2022 17:22
Function object for strain rate tensor? blebon OpenFOAM Post-Processing 18 June 24, 2021 03:50
decomposePar problem: Cell 0contains face labels out of range vaina74 OpenFOAM Pre-Processing 37 July 20, 2020 05:38
What is wall shear stress and shear strain rate? warlocklw Main CFD Forum 2 May 21, 2012 06:16
Wall function in adverse pressure gradients stephane baralon Main CFD Forum 11 September 2, 1999 04:05

All times are GMT -4. The time now is 00:45.