In the cold fusion experiment, there observed vortex patterns on a thick Pd electrode surface during longterm electrolysis in 0.1M LiOD. To understand the peculiar phenomenon, we have proposed N-cycle model, which is
composed of four sequential processes including the cold fusion (CF) reaction. In this study, a numerical simulation for the analysis of the vortex patterns was performed to elucidate the relation between the vortex formation and CF. Supposing that the hypothetical particles mass evolved due to CF reaction energy, two numerical simulations of the motion of the hypothetical particles mass were performed:
Cellular automata (CA) and discretization method. Results of both simulations have shown generation of a vortex or a cascade of vortices, which were identical to those obtained experimentally. However, there still exist differences between the experimental and theoretical patterns in scale and form. Then, it motivates us to build a more sophisticated model. On the other hand, Miesch et al. performed the numerical simulation of the convection structure in the solar convection zone, where the appearance of vortex was deduced from the helioseismic data. The numerical simulation showed that the strong
downflow was characterized by the helicity reversal near the base of the convection zone. Reviewing the numerical
calculation method and the restricted conditions inspired us the important term of helicity change under the conservation, relaxation and transportation across a boundary.