Numerical verification of Orbit and Nemato codes for magnetic topology diagnosis

We present the results of a benchmark study involving the Orbit and Nemato codes. The two codes have been used to compare magnetic structure in a reversed-field pinch (RFP), where conserved magnetic structures/islands appear both in the core (dominated by m=1 modes) and in the edge (dominated by m=0 modes). As input, a snapshot of a 3D nonlinear MHD visco-resistive simulation (produced by the SpeCyl code) has been used. The first test is given by the reconstruction via Poincare surface of section plot of an island generated by a single mode. In this case, the magnetic field topology corresponds to a time-independent Hamiltonian and shows conserved flux-surfaces used as a reference. Both codes successfully yield field lines which follow the same flux surfaces, in both the m=1 and m=0 cases. The benchmark between the codes has then been extended to a more complex configuration with chaotic magnetic field, using as input a fully 3D multiple mode RFP condition, characterized by the typical chain of edge magnetic islands providing a transport barrier. Finally, a quantitative benchmark has been performed, using the same 3D input, by estimating the correlation length of the magnetic field line in a bounded stochastic domain. The conclusion is that both codes yield consistent Poincare plot structure on one hand, and very good quantitative agreement in correlation length estimate. This gives confidence to the application of both codes to magnetic topology in the RFX-mod device, for which they are used routinely, as well as for a generic fusion device.