[wiking69]

Hi all,

I have searched the code but I can't seem to find the line where the boundary conditions for the k- or R-field are enforced.

There are plenty such lines for epsilon / omega / nut such as

this->nut_.correctBoundaryConditions();
or: omega_.boundaryFieldRef().updateCoeffs();
or: omegaEqn.ref().boundaryManipulate(omega_.boundaryF ieldRef());

but none for k / R.

I had a case where R on a wall-boundary started to go wild and didn't rest at fixedValue uniform (0 0 0 0 0 0)

To prevent questions, I'm using a LRR/SSG-omega Reynods-Stress model that is SST-like and is supposed to work without wall functions for R.

Gladly waiting for your hints.

Hendrik