Numerical Treatment of Cylindrical Coordinate Centerline Singularities