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,
https://www.openfoam.com/documentation/guides/latest/api/bound_8C_source.html
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.
Regards!