In most previous ice–ship interaction studies involving fluid effects, ice was taken as unbreakable. Building breakable level ice on water domain is still a big challenge in numerical simulation. This paper overcomes this difficulty and presents a numerical modeling of a ship moving in level ice on the water by using a one-way CFD-DEM (computational fluid dynamics-discrete element method) coupling method. The detailed numerical processes and techniques are introduced. The ice crack propagation process including radial and circular cracks have been observed. Numerical results are compared with previous experimental data and good agreement has been achieved. The results show that water resistance is an order of magnitude smaller than ice resistance during the ice-breaking process. Ice resistance shows strong oscillation along with ice failure process, which are affected by ship speed and ice thickness significantly.

In this work, the finite element method (PD-FEM) coupling strategy is used to simulate ship-ice interaction. Two numerical benchmark tests are selected to validate the coupling approach and its program. During the ice-breaking process simulation, the generation and propagation of radial and circular cracks in level ice are modeled and phenomena such as the shedding of wedge ice, flipping of brash ice, and cleaning of the channel are observed to be broadly consistent with experimental observation. The influence of ship speed and ice thickness on the ice load are investigated and analyzed. The ice load obtained from the numerical simulations is in general agreement with that given by Lindqvist’s empirical formula. The boundary effect on the crack path can also be avoid with the current coupling method.

When an airgun releases high-pressure gas underwater below an ice plate, it is observed that a bubble is formed rapidly while the ice plate is broken fiercely. In order to study the ice-water-gas interaction during this transient and violent phenomenon, a set of laboratory-scale devices was designed and a series of icebreaking experiments were carried out. High-speed photography was used to capture the evolution of the bubble and the ice plate. It was found that the airgun bubble had a unique ‘pear’ shape compared with the spherical bubble generated by electric sparking. The pressure induced by the pulsation of the airgun bubble near a rigid wall was measured by the pressure sensor. The initial shockwave, oscillatory pressure peaks caused by the directional fast air injection, secondary shockwave, and pressure peak caused by the bubble jet impact were clearly recorded. Three damage patterns of ice plates were observed and corresponding reasons were analyzed. The influence of dimensionless parameters, such as airgun-ice distance $H$ and ice thickness $T$, was also investigated. The physical mechanism of ice-water-gas interaction was summarized.

Icebreaking by using underwater explosion bubbles and compressed high-pressure gas bubbles has gradually become an effective icebreaking method. In order to compare the damaging effect of these two methods on the ice body, a fluid–structure coupling model was established based on the arbitrary Lagrangian–Eulerian (ALE) method and a series of calculations were carried out. The morphological changes of underwater explosion bubbles and compressed gas bubbles at the same energy under the free surface; the changes of flow load near the rigid wall; and the damage caused to the ice plate were studied and compared. The damage effect of the ice plate was analyzed by detecting the number of failure elements of the ice plate, and the optimum standoff distance was found. For an ice plate with a radius of 0.19 m and a thickness of 0.15 m, the optimum standoff distance of the compressed gas bubbles with 120 J is 0.03 m, and the optimum standoff distance of the TNT with 120 J is 0.02875 m. The similarities and differences of the two sources of bubbles on ice plate damage were summarized.

A free-rising buoyant sphere can break an ice plate floating above it. The problem is when the light sphere breaks the ice plate most, or the optimal relative density of the sphere which can break the ice plate the most severely. This experimental study was done to answer this problem. A set of experimental devices were designed, and a high-speed camera system was adopted to record the whole dynamic process, including the free-rising of the sphere, the collision between the sphere and the ice plate, the crack initiation and propagation, as well as the breakup of the ice plate. The failure mode of the ice plate under impact load was analyzed. It was found that conical cracks were formed under the reflected tensile wave at the top surface of the ice plate. On this basis, the influences of ice thickness, the initial submergence depth, and the relative density of the sphere on icebreaking were further investigated. An optimal relative density of the sphere was found when the sphere was released at a certain initial submergence depth, at which point the ice was damaged the most severely. For example, when the dimensionless initial submergence depth of the sphere was 2.31, the optimal relative density of the sphere was close to 0.4, with the probability of the ice plate breakup as high as 91.7%. It was also found from the experiments that the degree of damage to the ice plate correlated well with the kinetic energy of the sphere just before collision. Results showed that the optimal relative density can be estimated by theoretical analysis of the kinetic energy of the sphere, which will provide a reference for potential icebreaking applications in the future.