In the cold fusion experiment, there observed vortex patterns on a thick Pd electrode surface during long-term electrolysis in 0.1M LiOD. To understand the peculiar phenomenon, we had proposed N-cycle model, which is composed of four sequential processes including the CF reaction. In this study, numerical simulation for the analysis of the vortex patterns was performed to elucidate the relation between the vortex formation and the CF reaction. Supposing that the hypothetical particles mass evolved due to the CF reaction energy, we performed the numerical simulation using discretization method for the analysis of 3D motion of the hypothetical particles mass. The experimental results showed that the ejected particles mass from one side to another in a given rectangular box (simulate the sub-surface layer) exhibited a helical trajectory as a function of magnetic field. The traces of the ejected particles mass on another side composed vortex pattern. There has been remained more accurate simulation to examine the structure of the layer and the influence of the magnetic field.