In this paper, a coupled computational fluid dynamics-discrete element (CFD-DEM) approach is used to investigate the performance of abrasive slurry jets. Firstly, the CFD method was used to simulate the VES slurry jet flow characteristics. Secondly, an orthogonal design method was used to analyse the sensitivity of the parameters of zero shear viscosity , relaxation time and flow index of the Cross fluid, in which zero shear viscosity had the most significant effect on the flow field. As the zero-shear viscosity increases, the velocity of the flow field gradually increases, but the core section of the jet becomes less concentrated and the core section of the jet decreases. Based on the VES fluid formulation used in this project, the optimum zero-shear viscosity, relaxation time and flow index were optimised to 2.5 Pa-s, -0.1s and 0.2 respectively. The results showed that the addition of VES surfactant enhanced the erosion rate of the specimens and that the erosion rate increased with increasing OTAC concentration.