An accurate and robust axisymmetric SPH method based on Riemann solver with applications in ocean engineering