[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

[David 1985]

I am having the same question.
Why should it not be necessary to reinforce the boundary condition for k?


I found out that adding the terms

k_.boundaryFieldRef().updateCoeffs();

and later

k.ref().boundaryManipulate(k_.boundaryFieldRef());

improves convergence in the v2f model in the Obi diffuser case.


So I ask myself, if there is a reason behind leaving the corresponding equations for k out.

[dlahaye]

turbulence.correct() solves turbulence eqns, including boundary conditions handling.