Numerical simulation of the linear and second-order surface flows around circular cylinders in random waves