Numerical methods for matrix flows on Grassmann manifold