Robust and stable discrete adjoint solver development for shape optimisation of incompressible flows with industrial applications