An enhanced implicit viscosity ISPH method for simulating free-surface flow coupled with solid-liquid phase change.

Free-surface flow coupled with solid-liquid phase change can be observed in many industrial applications, and it has a great influence on many industrial processes. In this paper, the Incompressible Smoothed Particle Hydrodynamics (ISPH) method is extended with solid-liquid phase change model. A novel implicit viscosity solver is adopted in which the implicit viscosity solver is moved behind the calculation of pressure force to eliminate the time step limitation and the numerical creeping of high viscosity. To reduce the compression caused by the delayed incompressibility correction, accurate wall boundary conditions are applied and an enhanced implicit viscosity Laplacian operator is proposed. Besides, a dynamic solid boundary criterion is proposed to deal with the changing solid-liquid interface. In addition, the enthalpy-based heat transfer model and the mushy zone viscosity model are introduced. The inter-particle conductivity between liquid and solid is adopted to consider the large contact thermal resistance between fluid and solid. As will be shown, our ISPH method is verified in its performance in simulating free-surface flow coupled with solid-liquid phase change by four numerical examples and has good advantages in simulation quality over the existing implicit viscosity algorithm, which, hopefully, is to provide applications in engineering. • ISPH method is developed for free-surface flow coupled with solid-liquid phase change. • Enhanced implicit viscosity algorithm is proposed for high viscosity flow. • Dynamic solid boundary criterion is proposed for the changing solid-liquid interface. • The contact thermal resistance model and the mushy zone viscosity model is introduced. • Numerical examples verify the method and demonstrate the improvements. [ABSTRACT FROM AUTHOR]