We propose a multi-phase, multi-resolution SPH method for fluid/solid interaction with multi-GPU implementation and dynamic load balancing following the movement of the refinement regions. The primary design goal of this framework is to maintain the efficiency of the single-resolution SPH model running on a single GPU. To this end, a multi-background mesh is introduced, and the domain is regarded as a nested multi-domain with different resolutions. Validation using both a δ-SPH and Riemann SPH model is shown, and applications to the simulation of the water entry of a projectile with a high Froude number are considered, with comparisons to experimental data from three challenging test cases, showing the proposed model's ability to correctly reproduce the free surface evolution on water entry, the motion of the projectile, and the formation and evolution of multiple cavities depending on entry angle and velocity. An analysis of the computational performance and resolutions achieved (up to 120 million particles) is also provided across several test cases. •A multi-GPU, multi-node implementation is developed using CUDA for the GPU execution.•The efficiencies of the multi-GPU and multi-resolution model are presented and validated.•The two-phase flow of three water entry cases of slender body is simulated using the model.•The SPH model can provide a sharp interface between the air and water for water entry problems.