A Continuum-Based Fully-Lagrangian Model for Ice-Water Dynamics

Recognition of the behavior of the river ice dynamics is of great importance, especially in cold regions, as they have the potential to trigger the so-called ice jamming phenomena. Since these phenomena are highly dynamic, its dynamics can only be captured by proper numerical methods. Numerical methods for ice dynamic problems are either based on the discrete or continuum descriptions of the ice floes. While the discrete-based approaches, such as the Discrete Element Method (DEM), are very accurate, they are computationally expensive for large-scale problems, as they deal with individual ice floes. On the other hand, the continuum description offers scalability and computational affordability, by considering the assembly of ice parcels as a body of continuum and solving the conservation equations. In combination with the mesh-free Lagrangian (particle) numerical methods, such as Smoothed Particle Hydrodynamic (SPH) and Moving Particle Semi implicit (MPS), the continuum description can also offer accuracy and flexibility comparable to the discrete description. In this study, we propose an MPS Lagrangian numerical method for simulation of multiphase, multiphasic ice-water dynamics, based on continuum description of both ice and water, considered as non-Newtonian and Newtonian fluids, respectively.