[Andrea1984]

Hello,

I am trying to define a source term for the epsilon equation using a scalarCodedSource in fvOptions.

I am using OF v6 and my solver is simpleFoam.

In my code I can easily access attributes of my mesh as, for instance:

const vectorField& C = mesh_.C();

to access the coordinates of the cell centres.

The source term I am trying to implement depends also on k and epsilon. Any idea on how I can access the turbulence variable fields from fvOptions?

I have tried both

const scalarField& eps = epsilon_;

and

const scalarField& eps = turbulence->epsilon();

but without any success.

Any help would be greatly appreciated.

Andrea

[jherb]

Try something like (see https://github.com/OpenFOAM/OpenFOAM...nsfer.C#L105):

const compressible::turbulenceModel &turb = mesh_.lookupObject<compressible::turbulenceModel>
        (
            turbulenceModel::propertiesName
        );
const volScalarField &eps = turb.epsion();

Perhaps you have to change it to a non-compressible turbulence class depending on your solver.

[Andrea1984]

Hi Joachim

thanks for your reply. I have tested your approach but it does not work because compressible::turbulenceModel (or incompressible::turbulenceModel for incompressible cases) is not a valid type to declare within fvOptions. Maybe this can be fixed with a few include statements but I could not find a quick fix.

However, I have discovered that the following syntax does the trick:

const scalarField& eps = mesh_.lookupObject<volScalarField>("epsilon");

this can be used without including any additional header files since mesh_ is already available in fvOptions.

Andrea