[Phicau]

Hi all,

I am currently working in porous media flow with free surface using a customized interFoam (OF2.1.0). So far the model has proven to be stable, but it runs for laminar turbulence model. Now I have to modify kEpsilon to add the closure terms needed. In order to do so I have to make kEpsilon to accept another field (porosity) at the time of creation, but that is kind of a secondary question for now.

As a first step I am trying to eliminate the unwanted functions e.g. R(), as I will only be using the model for myInterFoam. I changed the name of the model to kEpsilonPor and it compiles flawlessly. Then I removed all the references to R() in kEpsilonPor.H and kEpsilonPor.C and it did not compile.

To proceed I changed "files" and "options" from Make folder:
files

kEpsilonPor.C

LIB = $(FOAM_USER_LIBBIN)/libPorousRASModels

options

EXE_INC = \
    -I$(LIB_SRC)/finiteVolume/lnInclude \
    -I$(LIB_SRC)/meshTools/lnInclude

LIB_LIBS = \
    -lincompressibleTurbulenceModel \
    -lfiniteVolume \
    -lmeshTools

and copied several files to the folder, in order to modify them: RASModel.C, RASModel.H, transportModel.C, transportModel.H, turbulenceModel.C, turbulenceModel.H, backwardsCompatibilityWallFunctions.C, backwardsCompatibilityWallFunctions.H and backwardsCompatibilityWallFunctionsTemplates.C.

Everything compiles well, getting rid of the mentioned references to R() in RASModel.C, RASModel.H and turbulenceModel.H.

The problem is that when I run any case using the new turbulence model I keep getting:

kEpsilonPorCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    sigmaEps        1.3;
}

time step after time step, not only the first one, as expected.

Any hints on what could be happening and how to fix it?

Thanks

Pablo