added edge cases in k for div(u)
Different output and another sign change
corrected non-dimensionaliezed version

and corrected pressure control