The multistage centrifugal pump is the critical component of mineral resources lifting in deep-sea mining. The reflux of nodules in the lifting pipe caused by the emergency pump stop can easily cause the pump to clog. In this paper, coupled Computational Fluid Dynamics and Discrete Element Method (CFD-DEM co-simulations) are used to clarify the solid-liquid two-phase flow in two-stage centrifugal pumps under different particle sizes (10–20, 20–30, 30–40, 40–50 mm) with constant particle concentration. The movement and accumulation behaviour of particles in different flow fields (pipeline to pump, the first to the second pump stage) is investigated. Meanwhile, the effect of particle size and particle reflux velocity on the blockage of the flow channel in the pump was investigated. Particle accumulation in the pump was observed to determine the key factors affecting the pump’s reflux capacity. The residual mass of particles in the pump at different particle sizes was counted. Simultaneously, the percentage of residual mass of 10–20 mm particles in the pump was compared between the experiment and the simulation with an acceptable tolerance of within 10%. In addition, pressure changes in the blockage-prone section were also investigated. A comparison between experiments and simulations verifies the consistency of the trend on the pump inlet pressure when clogged with 50 mm particles. It was found that larger particles in the range of 10–30 mm can better ensure the pump’s reflux performance.

Six-stage centrifugal pumps are used in deep-sea mining lifting systems and are required to convey slurry containing coarse particles. A six-stage centrifugal pump suitable for operation in a natural mining system was manufactured. High-flow and full-scaled slurry conveying experiments at a 5% and 9% volume concentration of particles was carried out at a large modified test site with artificial nodules. CFD-DEM simulations were carried out to obtain slurry transport characteristic curves, particle transport and distribution characteristics, where the simulation method was validated by the experiment data. A clarified two-stage pump can be used instead of a multi-stage pump for simplified simulation calculations with acceptable accuracy. Local agglomeration of particles caused by backflow was found at the outlet of the diffuser, and such agglomeration decreased with increasing flow rates. It was found that particles are transported non-uniformly, particles transport in diffusers in strands. Particles are transported in a pulse-like mode within the pump, with the latter stage showing similar particle characteristics to those transported in the previous pump stage.

The seafloor soil is characterized by high water content, strong compressibility, and low shear strength. Deep-sea mining vehicles (DSMV) are prone to sinking when walking on the surface of the soil, which will cause significant reduction in traction performance. Therefore, it is necessary to study the sinkage performance. The track is usually considered the travelling mechanism of the DSMV, and the track plate is an important part of the movement system. The study of the interaction between the track plate and the soil is of great significance to the study of the DSMV’s sinkage performance. In this study, firstly, based on the in situ seafloor soil samples of 1000 m in a region of the South China Sea collected by a box sampler, the physical and mechanical parameters of soil were measured by indoor geotechnical instruments. Secondly, an elastoplastic soil numerical model similar to that of in situ soil was established. Based on coupled Eulerian-Lagrangian (CEL) method, a numerical model of the interaction between the track plate and soil was established. Considering the dynamic process, the structure of the track plate and the physical and mechanical properties of the soil, the numerical simulation were carried out under different conditions, such as different dynamic loading, the plate structural parameters and the soil physical and mechanical properties. It is found that the plate-sinkage curve were significantly influenced by these factors. The findings are as follows, firstly, with the increase in the pressure loading rate, the soil sinkage decreasing at the same pressure. On the other hand, with the increase in velocity, soil flow was accelerated, and the nonlinear relationship between resistance and velocity became more obvious; the L/B ratio of different track plates affects the variation law of the curve, and the maximum sinkage gradually decreases as the ratio of L/B increases; with the increase in the grouser height, the maximum sinkage gradually decreases, and the pressure-sinkage curve changes obviously with the grouser type; and different soil physical and mechanical properties affect the variation of pressure-sinkage curve. Innovatively, the heterogeneous soil stress distribution mode was obtained through the fitting function and Python secondary development. This study can provide a reference for studying the sinkage performance of the DSMV.

The vortex-induced vibration (VIV) dynamics of commercial-scale deep-sea mining risers with complex component arrangements (pumps, buffer stations, buoyancy modules) remain insufficiently explored, especially for 6000 m systems with nonlinear tension. This study investigates VIV control strategy by adjusting tension for a nonlinear riser system using the Iwan-Blevins wake oscillator model integrated with Morison equation-based analysis. An analytical model incorporating four typical current profiles was established to quantify the dynamic response under different flow velocities, internal flow density, and structural parameters. Increased buffer station mass effectively suppressed drift distance (over 35% reduction under specific conditions) by regulating axial tension. Dynamic comparisons demonstrated distinct VIV energy distribution patterns under different current conditions. Spectral analysis revealed that the vibration follows Strouhal vortex shedding lock-in principles. Spatial modal differentiation was observed due to nonlinear variations in velocity profiles, pipe diameters, and axial tension, accompanied by multi-frequency resonance, coexistence of standing and traveling waves, and broadband resonance with amplitude surges under critical velocities (1.75 m/s in Current-B). This study proposes to control the VIV amplitude by adjusting internal flow density and buffer mass, which is proved effective for reducing vibrations in upper (0–2000 m) risers. It validates vibration amplitude and frequency control through current velocity, buffer mass and slurry density regulation in a nonlinear riser system.