[gu1]

Hello,

I would like to know if it is possible, within the kEpsilon model to implement a boundary condition for Epsilon, similar to the example below:

epsilon_
(
IOobject
(
IOobject::groupName("epsilon", U.group()),
U.time().timeName(),
U.mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
U.mesh(),
epsilon.bondaryCondition()

I need to implement the epsilonWallFunction for my incompressible model. However, because of my code, when I compile the BC (epsilonWallFunction) and call it in my solver, the answer we get is that it is necessary to use a turbulence model (kEpsilon ... LamBr ...) and my algorithm despite being based on the OpenFOAM turbulence models, is an algorithm of its own, there is no need for the turbulenceProperties and transportProperties files.
So, I'm looking for ways to circumvent this, since completely redoing the epsilonWallFunction condition for my solver does not seem to me very simple. I wanted to know if only the most important parts is possible to put in the file '' kEpsilon.C ''.