In order to minimize numerical diffusion PIC/FLIP hybrid method which mixes the perspectives of solving the system from a particle point of view (Lagrangian) and solving the system from a grid point of view (Eulerian) is used. Particle displacement is calculated after each individual dike injection, while heat conduction is solved on a fixed grid.