Bound does not work in some boundaryFields
Hello folks we've found an issue with the bound method, particularly in the line,
vsf.boundaryFieldRef() = max(vsf.boundaryField(), lowerBound.value());
Line 58 here,
In the attached case (run it with pimpleFoam, v2006) an inletOutlet BC is set in the outlet for omega with inletValue 0. Bound should override this value with 1E-5 as set in constant/turbulenceProperties and in fact that happens in the RHS of line 58, but after "=" operator evaluation the original BC inletValue = 0 is set (or left). Due to this 0 value in omega solver crashes immediately after bound.