Particle-resolved simulations of ice slurry flows in a square duct by the thermal immersed boundary–lattice Boltzmann method