Abstract: As an important component transporting resources such as oil and mineral ores mixture from the seabed to the surface in ocean engineering, vortex-induced vibration (VIV) of flexible risers can be encountered when the risers are subjected to the external environmental conditions. As VIV can lead to structural fatigue for the riser system, which threatens to the facility safety during deepsea resource exploitation, it is of great significance to investigate VIV mechanism and dynamics. Therefore, VIV dynamics of a flexible fluid-conveying riser undergoing external shear current is studied based on the combination of the Euler-Bernoulli beam theory and the semi-empirical hydrodynamic model. The finite element method and Newmark-β method are adopted to discretize and solve the governing equation. The model is firstly validated by comparing with the experimental data in order to examine the accuracy of the present model. Subsequently, cross-flow (CF) VIV response of the fluid-conveying riser is mainly examined and analyzed while various internal flow velocity and fluid density are considered and changed. The results show that when the flexible riser is subjected to both internal flow and shear current, there appears multi-frequency response for CF VIV. And the CF vibrating frequency and the CF root mean square (RMS) displacement are evidently influenced by the internal flow velocity and fluid density. With the increase of the internal flow velocity and fluid density, the CF vibrating frequency decreases while the RMS displacement shows an increasing trend in CF direction. Furthermore, in addition to the variation of the CF vibrating frequency and RMS displacement, the change of internal flow densities can cause notable mode and frequency transitions.