Kinetic model solution for axisymmetric flow by the method of discrete ordinates