Ships sailing through cold regions frequently encounter floe ice fields. An air-bubble system that reduces friction between the hull and ice floes is thus considered useful for the reduction of ice-induced resistance. In this study, a numerical analysis procedure based on coupled finite volume method (FVM) and discrete element method (DEM) is proposed to simulate complicated hull-water-gas-ice interactions for ice-going ships installed with air-bubble systems. The simulations reveal that after turning on the air-bubble system ice floes in contact with the hull side wall are pushed away from the hull by the gas-water mixture, resulting in an ice-free zone close to the side hull. It is found that the drag reduction rate increases with the increase of ventilation, while the bow ventilation plays a deciding role in the overall ice-resistance reduction. The proposed procedure is expected to facilitate design of new generations of ice-going ships.