[gu1]

Hi,

I would like to know if it is possible to add a term as highlighted below (original PISO code):

fvVectorMatrix UEqn
(
fvm::ddt(U)
+ fvm::div(phi, U)
+ turbulence->divDevReff(U)
+ myTerm->divMyTerm(U)
==
fvOptions(U)
);

How could I do that?
I'm trying to add a term to the kEpsilon model, but unfortunately I'm having trouble understanding how the pisoFoam solver looks for this information contained in turbulence->divDevReff(U) and that are making it impossible for me to add another term.

And I have a second question, also related to the kEpsilon model. Based on original file: kEpsilon

On lines 132 and 133 I could substitute for:

bound(k_, kMin_);
bound(epsilon_, epsilonMin_);

for

bound(k_, SMALL);
bound(epsilon_, SMALL);

and lines 259 and 275:

bound(epsilon_, epsilonMin_);
bound(k_, kMin_);

for

bound(epsilon_, SMALL);
bound(k_, SMALL);

I could not understand if kMin and epsilonMin are updated every iteration. Could someone explain to me what their function is? Would I disrupt the efficiency of the kEpsilon model if I removed them?

OBS: I understand that SMALL is the default value of this function.