On a rational magnetic surface, an isochronous bifurcation transforms one island chain into another chain with the same winding number. This transformation has been the subject of recent studies in tokamak plasmas. Namely, visco-resistive magnetohydrodynamic simulations of NSTX-U and DIII-D plasmas showed the onset of bifurcations with new magnetic isochronous islands for two competing helical perturbations on the same rational magnetic surface. To investigate these bifurcations, we use a cylindrical plasma model, with first-order correction for toroidicity, subject to externally applied magnetic perturbations, generated by a pair of resonant helical windings (RHWs) on the external wall and superposed to a helical current sheet (HCS) located on a rational plasma surface. We numerically integrate the magnetic field line equation and show that isochronous islands emerge when the perturbation created by the HCS increases. We present examples of such bifurcations on primary and secondary magnetic surfaces for different RHW configurations.

In the phase space of twist Hamiltonian systems, the emergence of different isochronous island chains on the same rational surface is observed when multiple resonant perturbations with different wave numbers and the same winding numbers are applied.1–3 All the concerned islands in the chains, called isochronous islands,2 have the same winding number, which means they have the same average speed around an invariant circle. The transformation from one of these chains to another occurs through bifurcations called isochronous. Such isochronous islands and bifurcations have been found in several systems—see Refs. 1–3 and references cited therein.

In tokamaks, the magnetohydrodynamics equations for the magnetic field lines can be, in general, described as a twist Hamiltonian system after defining a suitable field line Hamiltonian.4 Thus, the emergence of isochronous bifurcations can be expected for multiple resonant perturbations with the same winding number, as considered in Refs. 5 and 6. In a Poincaré section, the magnetic island chains can be obtained by integrating the field line equation.

Recently, studies with visco-resistive magnetohydrodynamic (MHD) simulations have investigated the role of resonant magnetic perturbations (RMPs) on the onset of magnetic islands in the isochronous bifurcations at the NSTX-U tokamak.5,6 In Ref. 5, a series of Poincaré sections were calculated by following magnetic field lines for two different approaches: (i) superposing the equilibrium axisymmetric field with the externally applied magnetic perturbation, but without considering the plasma response to this perturbation (the so-called vacuum approach); and (ii) superposing the equilibrium axisymmetric field with both the externally applied magnetic perturbation and the field due to the response of a visco-resistive single-fluid plasma. Perturbations create the reported islands due to the application of resonant external coils on the plasma equilibrium perturbed by an internal mode described by a linear resistive MHD approximation.5,6 These calculations show the existence of island bifurcations. Evans et al.,6 however, have shown that island elongation is followed by a sequence of isochronous bifurcations when the perturbation current is increased.

Experimental observation of isochronous island chains during DIII-D plasma discharges has also been reported, reinforcing that further studies on the impact of RMPs are required to control these heteroclinic structures.7 These studies described a bifurcation from heteroclinic to homoclinic topology in DIII-D plasmas.

To study these effects on magnetic topology, we use a model that can generate isochronous bifurcations similar to the ones above when two RMPs are superposed for a cylindrical plasma with first-order correction on toroidicity. The first perturbation is created by a pair of resonant helical windings (RHWs) located on the external wall. It is usually applied to reduce wall erosion by generating a chaotic layer in the plasma edge.8–13 The second perturbation is created by a helical current sheet (HCS) located on a rational flux surface that acts as a plasma current perturbation.14,15 Since a first-order correction on toroidicity is applied to this cylindrical plasma model, an RMP field will result in the formation of magnetic islands at the resonant rational surface (the primary mode associated with the RMP) and also on the neighboring rational surfaces (the secondary modes associated with the RMP due to poloidal mode coupling induced by toroidicity). Similarly to the study reported in Ref. 5, the resulting magnetic topology around a rational surface in this work is a combination of primary and secondary modes caused by RHWs and HCS. Here, we estimate the current amplitude required for the onset of isochronous islands around rational surfaces.

Additionally, the analytical expressions for the magnetic field are obtained for a cylindrical plasma in equilibrium with first-order corrections on toroidicity included and with the superposition of perturbations due to a pair of RHWs16 and an HCS.17–20 Solving the field line equations allows us to calculate the Poincaré maps for a plasma cross section, considering different combinations of parameters for specified perturbations. We have identified isochronous bifurcations located around rational flux surfaces, all of them modifying the distribution of field lines similarly to the bifurcations observed in visco-resistive MHD simulations5,6 and in experiments.7 

This paper is structured as follows: Sec. II presents the MHD equilibrium for a cylindrical plasma with first-order correction due to toroidicity. The magnetic field generated by a pair of RHWs and an HCS is calculated in Secs. III and IV, respectively. Section V shows the Poincaré sections of field lines, while conclusions are presented in Sec. VI. The  Appendix details the numerical procedure used to calculate the Poincaré sections and to identify the isochronous bifurcations of Sec. V.

The static MHD equilibrium for a low beta cylindrical plasma is described by the Grad–Shafranov equation in cylindrical coordinates (r, θ, z),17,18
(1)
where Ψ is the magnetic poloidal flux, μ 0 is the vacuum magnetic permeability, and Jz is the current density. Within this approximation, Eq. (1) can be obtained from the Grad–Shafranov equation as the lowest order approximation in local coordinates, expanding in powers of inverse aspect ratio.21 
In this work, we suppose that Jz peaks at the plasma column center (r = 0) and vanishes at the plasma edge,
(2)
where Ip is the plasma current, a is the plasma radius, and γ is a constant. Since the plasma is confined inside a cylindrical chamber (radius b), we also consider that a < b.
The equilibrium magnetic field ( B 0) for a cylindrical plasma can be written as
(3)
where B 0 θ and B 0 z are the angular and axial components, expressed in T/m and T, respectively, and the vectors e r, e θ, and e z form a non-normalized basis that satisfies the relation e z × e r = e θ / r. The equilibrium magnetic field is calculated by solving Eq. (1) with the current profile shown in Eq. (2), resulting in the following expression:
(4)
where Ie is an equivalent current generating the toroidal magnetic field, and 2 π R 0 is the periodic length of the cylindrical plasma along the z axis.
In this work, we apply a toroidal correction to the axial component of the equilibrium magnetic field,19 
(5)
The safety factor is defined as16 
(6)
where qc(r), the safety factor for a cylindrical plasma, is calculated by the following expression:
(7)

The parameters used in Secs. III–VI are related to the operation of the TCABR tokamak:22  Ip = 100 kA, Ie = 4 MA, a = 0.2 m, R0 = 0.6 m, and γ = 3. Figure 1 shows the safety factor radial profile considering Eq. (6) and the TCABR parameters. This figure also shows that radial positions (r/a) 0.610, 0.800, and 0.925 are approximately associated with safety factor values of 2, 3, and 4, respectively. In the following sections, trajectories of magnetic field lines around these radial positions are analyzed.

FIG. 1.

Safety factor radial profile across the plasma column considering a cylindrical plasma with toroidal correction.

FIG. 1.

Safety factor radial profile across the plasma column considering a cylindrical plasma with toroidal correction.

Close modal

A pair of RHWs generates an external perturbation that is superposed to the equilibrium magnetic field. It consists of two helical wires located on the surface r = b separated by 180° and with current Ih in opposing directions, as indicated in Fig. 2.16,23,24

FIG. 2.

Schematic diagram of RHWs on the surface r = b. The red and blue wires conduct current Ih in opposing directions.

FIG. 2.

Schematic diagram of RHWs on the surface r = b. The red and blue wires conduct current Ih in opposing directions.

Close modal
We define the winding law as
(8)
where mh and nh are the applied field's poloidal and toroidal mode numbers, respectively. The RHWs with mode (mh,nh) are resonant on the surface, where q r h = m h / n h.
The current density J h is defined as16–18 
(9)
where I l h = I h / ( 2 π b ) is the perturbation current density, considering b = 0.22 m, and e hel is a vector defining the helix direction.16–18 

Although the current distributions in this work are located on a surface—in Eq. (9), for instance, the surface is r = b—considering a current layer of finite thickness would result in more realistic values for the resonance amplitudes and the corresponding Poincaré maps. However, these alterations should not change the resonance sequence and the observed bifurcations.

In order to calculate the magnetic field due to the RHWs ( B h), we use the fact that · B h ( r , θ , z ) = 0 and write B h ( r , θ , z ) = ϕ h ( r , θ , z ) in the region r b, resulting in the following equation:16–18 
(10)

In the low beta plasma considered in this article, the solution of Eq. (10) determines approximately the perturbing magnetic field in the vacuum approximation, obtained by neglecting the effect of viscous and resistive plasma reaction.

The solution of this Laplace equation is25,26
(11)
where C k z , k θ is a constant that can be calculated from boundary conditions and I k θ k z β r is the modified Bessel function of the first kind.18 
Thus, the perturbing magnetic field B h r , θ , z , created by the helical windings, can be calculated from the potential ϕ h for the region r < b,
(12)
(13)
(14)

This perturbation excites the main mode (mh,nh) on the magnetic rational surface associated with the safety factor q = mh/nh and other sideband modes due to the toroidal correction.

In this work, a perturbation to the plasma current is represented as an RMP generated by a HCS located at the equilibrium rational surface r = rs, with rs < b and q(rs) = ms/ns (ms and ns here are the mode numbers associated with this perturbation).15  Figure 3 shows a schematic diagram of the surfaces r = b and r = rs with a representation of the helical direction in which the current density is applied.

FIG. 3.

Schematic diagram of the HCS highlighting surfaces r = b (external wall, shown in black) and r = rs (HCS, shown in red). The single helix wound around the surface r = rs represents the direction in which the surface current density is applied.

FIG. 3.

Schematic diagram of the HCS highlighting surfaces r = b (external wall, shown in black) and r = rs (HCS, shown in red). The single helix wound around the surface r = rs represents the direction in which the surface current density is applied.

Close modal
The current density related to the HCS is described by14,15,19
(15)
where j (in A/m) is the amplitude of the current density, ϕ p is a phase angle, and e s is the helical vector. To compare the amplitude j with parameters such as the plasma current Ip, we also define the current I s = j 2 π r s. The helical vector e s is defined as
(16)

In Eq. (16), the parameter α s = B 0 θ / B 0 z defines the helical direction of the equilibrium magnetic field on the radial position r = rs.

The magnetic field due to the current sheet ( B s) is calculated by solving · B s ( r , θ , z ) = 0, considering that B s r , θ , z = ϕ s ( r , θ , z ) for r r s, which results in the following Laplace equation:18,
(17)
whose solution is18,25,26
(18)
where C k z , k θ i and C k z , k θ e are constants calculated from the boundary conditions and I k θ k z β r and K k θ k z β r are the modified Bessel functions of the first and second kind, respectively.18 
Thus, the perturbing magnetic field B s r , θ , z , created by the internal current density, can be calculated from the potential ϕ s in the region r < rs as follows:
(19)
(20)
(21)
In the region r > rs, the perturbing field components are
(22)
(23)
(24)

This perturbation excites the main mode (ms,ns) on the magnetic rational surface associated with the safety factor q = ms/ns and other sideband modes due to the toroidal correction.

In this section, we analyze the effect of the superposition of two perturbations—one generated by a pair of RHWs around the plasma and the other one by an HCS, as modeled in Eq. (12)—on the magnetic field topology by solving the magnetic field line equation,
(25)
where d l is an infinitesimal displacement along the field line.

The pair of RHWs is resonant to mode (mh,nh) = (2,1) and has a perturbation parameter Ih/Ip = 0.4%, while the plasma current perturbation is resonant to mode (ms,ns) = (8,2) and has amplitude j. The chosen wave numbers are among those observed in usual tokamak discharges. The toroidal correction of our cylindrical approach introduces secondary sideband resonant modes interacting with the dominant modes.

In this article, the “vacuum approximation” has been used to calculate the magnetic field structure, although plasma response effects may result in a screening (or even an amplification) of the external perturbations. However, for low beta plasma and non-marginally stable modes, as considered here, this effect is not so noticeable.27,28

A set of initial values of radial and angular coordinates is chosen at the plane z = 0, and Eq. (25) is integrated numerically to obtain the Poincaré sections showing the intersection of field lines with this plane. The chosen wave numbers are among those observed in usual tokamak discharges. The  Appendix provides a detailed explanation of the numerical procedure used to solve Eq. (25), generate the Poincaré maps, and identify the bifurcations.

A complete toroidal description would make the identification of bifurcations more difficult than our cylindrical approach with toroidal corrections, which introduces the sideband modes but underestimates their amplitudes. In fact, for large-amplitude perturbations, the chaotic area increases while the island size is reduced, making the identification of bifurcations more difficult.

In the Poincaré sections, we consider the isochronous island chains (islands with the same winding number) created by the primary and secondary modes. Thus, we analyze the island chains on the q = mh/nh magnetic surface due to the primary mode (mh,nh) generated by the external RHWs and the secondary mode (ms,ns) generated by the internal plasma current perturbation. The competition between these primary and secondary modes gives rise to a sequence of isochronous bifurcations, on the considered rational surfaces, as the mode amplitude changes. These bifurcations modify the number of island chains formed by isochronous islands. As identified in Secs. V B–V E, the changes in the island chains occur due to local bifurcations29 that change fixed point stabilities and create new fixed points,3 changing the number of island chains.

Initially, we analyze the distribution of field lines across the z = 0 plane when the perturbation is generated only by RHWs (in this first case, the perturbation due to the HCS is not considered) and then by an HCS (in this second case, the perturbation due to the RHWs is not considered).

Figure 4(a) shows the Poincaré section associated with the first case, considering Ih/Ip = 0.4% and Is = 0. Two primary magnetic islands associated with the RHWs primary mode are formed at the q = 2 surface (r/a = 0.610). The toroidal correction applied to the axial component of the equilibrium magnetic field causes smaller secondary islands to appear on the neighboring rational surfaces: four magnetic islands are located at the q = 4 surface (r/a = 0.925), while three islands are observed at the q = 3 surface (r/a = 0.800).

FIG. 4.

Poincaré section of field lines for perturbation generated by (a) RHWs with (mh,nh) = (2,1) for Ih/Ip = 0.4% and (b) HCS with (ms,ns) = (8,2) for Is/Ip = 0.116%.

FIG. 4.

Poincaré section of field lines for perturbation generated by (a) RHWs with (mh,nh) = (2,1) for Ih/Ip = 0.4% and (b) HCS with (ms,ns) = (8,2) for Is/Ip = 0.116%.

Close modal

When only the HCS is present, as shown in Fig. 4(b), eight primary magnetic islands are formed at the q = 4 surface (r/a = 0.925). Secondary islands also form on different rational surfaces due to toroidicity induced poloidal mode coupling. Note that two sets of isochronous islands appear on both the q = 4 surface (r/a = 0.925—green and magenta islands), associated with the HCS primary mode, and on the q = 3 surface (r/a = 0.800—red and blue islands), associated with an HCS secondary mode. In fact, as indicated in Fig. 4(b), field lines located on the green islands do not intercept the magenta ones at the q = 4 surface, while blue islands do not intercept the red ones at the q = 3 surface.

When the superposition of perturbations due to RHWs and HCS is considered, as shown in the Poincaré section in Fig. 5, two main magnetic islands are formed at the q = 2 surface (r/a = 0.610). Smaller islands are formed on neighboring rational surfaces, indicating that the trajectory of field lines is similar to the case without the HCS perturbation, as shown in Fig. 4(a). The magnetic islands located around r/a = 0.925 and r/a = 0.800, associated with safety factor values 4 and 3, respectively, are only slightly distorted compared to the Poincaré section in Fig. 4(a). However, as will be shown in the next figures, these islands bifurcate, generating new islands if the amplitude Is is large enough. In these two rational surfaces, we observe isochronous islands when the perturbation is generated only by an HCS located at the q = 4 surface (rs/a = 0.925), as shown in Fig. 4(b).

FIG. 5.

Poincaré section of field lines for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4% and plasma current perturbation modeled as an HCS with (ms,ns) = (8,2) and Is/Ip = 0.116%.

FIG. 5.

Poincaré section of field lines for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4% and plasma current perturbation modeled as an HCS with (ms,ns) = (8,2) and Is/Ip = 0.116%.

Close modal

In this section, we consider the primary modes (mh,nh) = (2,1) and (ms,ns) = (8,2) and present examples of topological bifurcations that are generated when the current Is is increased and Ih is kept constant. In all cases considered here, the perturbation is generated by RHWs and HCS, with the phase angle ( ϕ p) equal to zero. The examples correspond to bifurcations observed at two different rational surfaces with q = 3 and q = 4.

When the Is/Ip is increased from 0.116% to 0.232%, while Ih/Ip is kept constant at 0.4%, a bifurcation is observed at the q = 3 surface (r/a = 0.800), as shown in the Poincaré section in Fig. 6(a). In this case, each of the three main magnetic islands (represented in red), which also appear in Figs. 4(a) and 5, now contains two additional smaller isochronous islands (green and blue), with the magnetic topology resulting in a combination of secondary modes associated with the RHWs (three red islands) and HCS [six isochronous islands in blue and green, also seen in Fig. 4(b)]. Field lines on the green region do not access or intercept the blue one. The observed elongation and bifurcation of magnetic islands is consistent with isochronous bifurcations observed experimentally7 and numerically investigated with resistive MHD simulations.5,6 Figure 6(b) shows an enlarged view of the Poincaré section around r/a = 0.800. When Is/Ip is increased to 4.65%, as shown in Fig. 6(c), the green and blue islands become larger and, instead of the main magnetic islands observed in Figs. 6(a) and 6(b), there is a thick chaotic layer around r/a = 0.800, highlighting the predominant role of the HCS secondary mode (six isochronous islands) in shaping the magnetic topology around this rational surface.

FIG. 6.

Poincaré sections for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and HCS with (ms,ns) = (8,2) and (a) Is/Ip = 0.232%, showing a region around rh/a = 0.610, associated with safety factor q(rh) = 2; (b) Is/Ip = 0.232%, showing a region around r/a = 0.800, associated with safety factor q(r) = 3; and (c) Is/Ip = 4.65%, showing a region around r/a = 0.800.

FIG. 6.

Poincaré sections for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and HCS with (ms,ns) = (8,2) and (a) Is/Ip = 0.232%, showing a region around rh/a = 0.610, associated with safety factor q(rh) = 2; (b) Is/Ip = 0.232%, showing a region around r/a = 0.800, associated with safety factor q(r) = 3; and (c) Is/Ip = 4.65%, showing a region around r/a = 0.800.

Close modal

Figure 7 shows a series of Poincaré sections of field lines around the q = 4 surface (r/a = 0.925), considering that Is/Ip is modified in the range from 0.291% to 0.407%, while the Ih/Ip is kept at 0.4% for all cases. In Fig. 7(a), four magnetic islands (represented in magenta) are located around r/a = 0.925, indicating that the RHW secondary mode (four islands) is predominant on topology. The distortion of field lines observed around r/a = 0.925 is similar to the modification of magnetic islands around rational surfaces observed in resistive MHD simulations (see Fig. 3 of Ref. 5).

FIG. 7.

Poincaré sections of field lines in the region around r/a = 0.925, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.291%, (b) 0.314%, (c) 0.320%, and (d) 0.407%. The magenta islands correspond to those originally generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

FIG. 7.

Poincaré sections of field lines in the region around r/a = 0.925, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.291%, (b) 0.314%, (c) 0.320%, and (d) 0.407%. The magenta islands correspond to those originally generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

Close modal

When Is/Ip is increased to 0.314%, as shown in Fig. 7(b), a bifurcation is observed in the q = 4 surface: two small islands (green and blue) are formed inside the main ones (magenta) around r/a = 0.925, which is also the region where two sets of isochronous islands are located when only the perturbation due to the HCS is considered, as shown in Fig. 4(b). In Fig. 7(c), it is shown that, when Is/Ip is modified to 0.320%, another bifurcation is observed: the secondary (green and blue) islands are completely separated from the main islands (magenta), and they are delimited by a magnetic surface (red) that does not connect to the main islands. Thus, we have one chain with 4 islands and the other with 8, a total of 12 islands. However, when the Is/Ip is increased to 0.407%, only two sets of isochronous islands are observed around r/a = 0.925, represented in magenta (4 islands) and red (4 islands) in Fig. 7(d). In this case, the influence of the HCS secondary mode, characterized by two sets of 4 isochronous islands, is predominant on the magnetic topology [as observed in Fig. 4(b) for Ih = 0], without the secondary islands created by the RHWs in Fig. 4(a) for Is = 0.

Considering that a bifurcation is observed around the rational surface with q = 3 when Is/Ip = 0.232% and that another set of isochronous islands is observed around the surface with q = 4 when Is/Ip  0.314%, our results are consistent with those reported in Ref. 5, in which isochronous bifurcations are observed on different rational surfaces as the perturbation parameter was increased.

Figure 8 shows an enlarged view of the Poincaré sections in Fig. 7 in the region 0.923 ≤ rs/a ≤ 0.929 and 0.25 ≤ θ/(2π) ≤ 0.35, considering the same perturbation parameters of Fig. 7. The described bifurcation process can be confirmed; in other words, Figs. 8(a)–8(d), clarifying the observed local bifurcations, are magnifications of Figs. 7(a)–7(d) around a fixed point.

FIG. 8.

Poincaré sections of field lines in the region 0.923 ≤ r/a ≤ 0.929 and 0.20 ≤ θ/(2π) ≤ 0.40, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.291%, (b) 0.314%, (c) 0.320%, and (d) 0.407%. The magenta islands correspond to those initially generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

FIG. 8.

Poincaré sections of field lines in the region 0.923 ≤ r/a ≤ 0.929 and 0.20 ≤ θ/(2π) ≤ 0.40, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.291%, (b) 0.314%, (c) 0.320%, and (d) 0.407%. The magenta islands correspond to those initially generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

Close modal

In Sec. V B, we analyzed the effect of increasing Is on the distribution of field lines, supposing that the phase angle in Eq. (12) was ϕ p = 0. We now study the effect on the magnetic topology of a perturbation created by an HCS, located at the surface r/a = 0.925 with ϕ p = π and (ms,ns) = (8,2), superposed to a perturbation created by a pair of RHWs located at r = b with (mh,nh) = (2,1). While Ih/Ip is kept constant at 0.4%, we change the parameter Is/Ip from 0.116% to 0.232% and analyze the Poincaré sections of field lines, as shown in Fig. 9.

FIG. 9.

Poincaré sections of field lines in the region around r/a = 0.925, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2), finite phase angle ( ϕ p = π) and Is/Ip equal to (a) 0.116% and (b) 0.232%. The magenta islands correspond to those initially generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

FIG. 9.

Poincaré sections of field lines in the region around r/a = 0.925, for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2), finite phase angle ( ϕ p = π) and Is/Ip equal to (a) 0.116% and (b) 0.232%. The magenta islands correspond to those initially generated by the RHWs around r/a = 0.925, while the blue and green islands are formed when the current associated with the HCS is increased.

Close modal

Figure 9(a) shows that four secondary islands (associated with the RHWs) are formed, when Is/Ip = 0.116%, at r/a = 0.925. When Is/Ip is modified to 0.232%, a bifurcation is observed, as shown in Fig. 9(b): small islands (represented in green and blue) are formed inside the main magnetic islands (magenta). There are, however, two main differences between the bifurcations observed in Figs. 7(b) (null phase angle) and 9(b) (finite phase angle): first, the two magnetic islands represented as green and blue surfaces in Fig. 9(b), for finite phase angle, appeared when the current Is was set at a value approximately equal to 74% of the one associated with Fig. 7(b); second, the small islands (green and blue) in Fig. 7(b) are located outside the main magnetic islands [represented in magenta in Fig. 7(b)], while in Fig. 9(b), for a finite phase angle, they are formed inside the main islands.

This example shows that the phase angle between the two perturbing currents, which can be adjusted by modifying the RHW position, alters the bifurcations in the plasma.

The previous results were calculated considering that the RHWs were placed outside the plasma at r = b with the current described by Eq. (9). To investigate the effect of an HCS on field line topology when the placement of electric wires around the external surface is modified, the following current density for the RHWs is now considered:
(26)

The current expressed in Eq. (26) has a reversed direction in comparison to the one shown in Eq. (9). Considering that Ih/Ip = 0.4%, and for RHWs with (mh,nh) = (2,1), Fig. 10 shows the Poincaré section of field lines for the z = 0 plane, initially without the effect of the HCS. In Fig. 10(a), two islands are formed around r/a = 0.610, and smaller islands are generated around the radial positions with rational safety factor values due to the inclusion of toroidicity correction on the axial component of the magnetic equilibrium field. These islands are similar to the results shown in Fig. 4(a), with the difference that the positioning of magnetic islands is shifted along the angular direction due to the modification of electric wires on the external wall. In Fig. 10(b), we show that four magnetic islands, represented in magenta, are formed around rs/a = 0.925. In the following section, we analyze an isochronous bifurcation observed around this region when the HCS is considered.

FIG. 10.

Poincaré sections of field lines for reversed RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, showing the region (a) around rh/a = 0.610, associated with safety factor q(rh) = 2, and (b) an enlarged view around r/a = 0.925, associated with safety factor q(r) = 4. The magenta islands correspond to the originally generated islands by the RHWs around rs/a = 0.925.

FIG. 10.

Poincaré sections of field lines for reversed RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, showing the region (a) around rh/a = 0.610, associated with safety factor q(rh) = 2, and (b) an enlarged view around r/a = 0.925, associated with safety factor q(r) = 4. The magenta islands correspond to the originally generated islands by the RHWs around rs/a = 0.925.

Close modal

The influence of two perturbations on the distribution of magnetic field lines is analyzed in this section: the first one is caused by a pair of RHWs with (mh,nh) = (2,1) whose current is described by Eq. (26), while the second one is generated by a plasma current perturbation modeled as an HCS with zero phase angle ( ϕ p = 0) and (ms,ns) = (8,2).

Figure 11 shows the Poincaré plot of field lines for two cases: in the first one, shown in Fig. 11(a), Is/Ip is equal to 0.116%, while in the second one, shown in Fig. 11(b), Is/Ip is equal to 0.174%. In Fig. 11(a), the magnetic islands are distorted around r/a = 0.925 in comparison to the Poincaré section shown in Fig. 10(b), with the RHW secondary mode shaping the distribution of magnetic islands around this rational surface. When Is/Ip is increased to 0.174% as shown in Fig. 11(b), the effect of the HCS secondary mode on magnetic topology is amplified, and a new set of four red magnetic islands around r/a = 0.925 is observed. This type of bifurcation, caused by the HCS, is similar to the ones reported in numerical and experimental studies.5–7 

FIG. 11.

Poincaré sections of field lines for reversed RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.116% and (b) 0.174%.

FIG. 11.

Poincaré sections of field lines for reversed RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and for HCS with (ms,ns) = (8,2) and Is/Ip equal to (a) 0.116% and (b) 0.174%.

Close modal

In plasmas confined in large aspect ratio tokamaks, represented as a cylindrical geometry with toroidicity included as a first-order correction, we considered magnetic islands created by a pair of resonant helical windings (RHWs), located on the external wall, and a helical current sheet (HCS), representing a plasma current perturbation near the plasma edge. Separately, these perturbations created resonant islands on the plasma's rational magnetic surfaces. When the two perturbations acted together, the resultant island's configuration changed through bifurcations as the current amplitudes were varied. In this work, we investigated isochronous bifurcations on a given rational magnetic surface for which the previous and final islands are isochronous, i.e., have the same winding number.

Poincaré sections of magnetic field lines were calculated by integrating the field line equation, and topological bifurcations were observed on sets of magnetic islands when one perturbing current was modified. Isochronous bifurcations were triggered when the current amplitude associated with the HCS was large enough while the current on the RHWs was kept constant, reorganizing field lines on magnetic islands near the radial position where the HCS is located. In these bifurcations, the number of island chains changes from the one determined by the RHW resonance to the other determined by the HCS resonance.

Isochronous bifurcations were also observed when the current in the HCS was increased and displaced 180° along the angular direction, resulting in pairs of small magnetic islands appearing inside the main magnetic islands where the HCS is located. These bifurcations also appeared when the current on the RHWs was reversed; in this case, sets of isochronous islands were observed around the surface where the HCS is located.

The isochronous bifurcations presented in this work are similar to the ones observed in visco-resistive MHD simulations of NSTX-U plasmas and observed in DIII-D plasmas. In further research, our results may yield additional insight into how plasma transport is affected by isochronous bifurcations in these and other machines.

The authors thank the São Paulo Research Foundation (FAPESP, Brazil) under Grant Nos. 2022/08699-2 and 2018/03211-6 for financial support. The authors acknowledge Dr. Marisa Roberto for her invaluable contributions to this work.

The authors have no conflicts to disclose.

Andre C. Fraile Jr.: Conceptualization (equal); Data curation (lead); Formal analysis (equal); Methodology (equal); Software (lead); Writing – original draft (lead); Writing – review & editing (equal). Marisa Roberto: Conceptualization (equal); Formal analysis (equal); Methodology (equal); Writing – original draft (supporting); Writing – review & editing (equal). Gustavo P. Canal: Formal analysis (equal); Writing – review & editing (equal). Ibere L. Caldas: Conceptualization (equal); Formal analysis (equal); Methodology (equal); Writing – original draft (supporting); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

In this section, we present the numerical procedure used to calculate the Poincaré maps shown in Sec. V. The displacement along a magnetic field line was numerically calculated by solving Eq. (25) with the Dormand and Prince method30 with order 5(4) and an adaptive step. In order to calculate the Poincaré map for the plasma cross section z = 0, we choose a set of NIC initial conditions (r, θ, z) = (ri,θi,0), with i = 1, 2, …, NIC. Then, Eq. (25) is solved NIC times from z = 0 to z = 2πR0Nt, where Nt is the number of turns along the periodical length 2πR0. Each time the magnetic field line crosses the z = 0 section, the corresponding position (r, θ) is computed, and the Poincaré map is plotted with the set of NtNIC points.

As an example of this calculation, Fig. 12 shows the Poincaré map calculated with NIC = 22 and Nt = 2000 considering the superposition of two perturbations: a pair of RHWs with (mh,nh) = (2,1) for Ih/Ip = 0.2% and an HCS with (ms,ns) = (8,2) for Is/Ip = 0.058%.

FIG. 12.

Poincaré section of field lines for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.2%, and HCS with (ms,ns) = (8,2) and Is/Ip = 0.058%.

FIG. 12.

Poincaré section of field lines for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.2%, and HCS with (ms,ns) = (8,2) and Is/Ip = 0.058%.

Close modal

As shown in Fig. 12, magnetic islands are formed around radial positions with rational safety factor values. After a Poincaré map such as that in Fig. 12 was obtained, we identified isochronous bifurcations (as those presented in Sec. V) by increasing the parameter Is/Ip and observing the distribution of magnetic islands around rational surfaces while the parameter Ih/Ip is kept constant.

While the results presented in Sec. V were obtained for Ih/Ip = 0.4%, similar isochronous bifurcations were observed for Ih/Ip ≤ 1.0%. However, Ih/Ip = 0.4% was chosen in this work because it is close to values commonly observed in discharges, and higher values of this parameter would result in chaotic regions, which are not the focus of this study. In fact, for Ih/Ip = 0.4% and ϕ p = 0 (with Is/Ip ≤ 0.6%, as shown in Sec. V B), we can estimate numerically, from Poincaré maps, a stochastic parameter (η) similar to the Chirikov parameter:31 considering the magnetic island chains associated with q(r4) = 4 and q(r9/2) = 9/2, we estimate their half-widths δ4 and δ9/2 and their radial positions r4 and r9/2, respectively. We define
(A1)

Figure 13 shows this stochastic parameter estimated for Ih/Ip = 0.4% and ϕ p = 0.

FIG. 13.

Numerical estimative of stochastic parameter for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and HCS with (ms,ns) = (8,2) and variable Is/Ip.

FIG. 13.

Numerical estimative of stochastic parameter for RHWs with (mh,nh) = (2,1) and Ih/Ip = 0.4%, and HCS with (ms,ns) = (8,2) and variable Is/Ip.

Close modal

Since η < 1 in the range Is/Ip ≤ 0.6%, neighboring separatrices do not intercept for the magnetic islands chains associated with q(r4) = 4 and q(r9/2) = 9/2.

1.
M.
Mugnaine
,
B. B.
Leal
,
A. M. O.
de Almeida
,
R. L.
Viana
, and
I. L.
Caldas
, arXiv:2312.12552 (
2023
).
2.
M. C.
de Sousa
,
I. L.
Caldas
,
A. M.
Ozorio de Almeida
,
F. B.
Rizzato
, and
R.
Pakter
,
Phys. Rev. E
88
,
064901
(
2013
).
3.
B. B.
Leal
,
I. L.
Caldas
,
M. C.
de Souza
,
R. L.
Viana
, and
A. M.
Ozorio de Almeida
,
Phys. Rev. E
109
,
014230
(
2024
).
4.
R. L.
Viana
,
M.
Mugnaine
, and
I. L.
Caldas
,
Phys. Plasmas
30
,
090901
(
2023
).
5.
W.
Wu
,
T. E.
Evans
,
G. P.
Canal
,
N. M.
Ferraro
,
B. C.
Lyons
, and
D. M.
Orlov
,
Nucl. Fusion
59
,
066010
(
2019
).
6.
T. E.
Evans
,
W.
Wu
,
G. P.
Canal
, and
N. M.
Ferraro
,
Phys. Rev. E
103
,
013209
(
2021
).
7.
L.
Bardóczi
and
T. E.
Evans
,
Phys. Rev. Lett.
126
,
085003
(
2021
).
8.
S. S.
Abdullaev
,
Magnetic Stochasticity in Magnetically Confined Fusion Plasmas
, 1st ed. (
Springer
,
Berlin
,
2014
), pp.
227
262
.
9.
M.
Lehnen
,
S. S.
Abdullaev
,
W.
Biel
,
S.
Brezinsek
,
K. H.
Finken
,
D.
Harting
,
M.
von Hellermann
,
M.
Jakubowski
,
R.
Jaspers
,
M.
Kobayashi
,
H. R.
Koslowski
,
A.
Krämer-Flecken
,
G.
Matsunaga
,
A.
Pospieszczyk
,
D.
Reiter
,
T.
Van Rompuy
,
U.
Samm
,
O.
Schmitz
,
G.
Sergienko
,
B.
Unterberg
,
R.
Wolf
,
O.
Zimmermann
, and
TEXTOR Team
,
J. Nucl. Mater.
337–339
,
171
(
2005
).
10.
K. H.
Finken
,
B.
Unterberg
,
Y.
Xu
,
S. S.
Abdullaev
,
M.
Jakubowski
,
M.
Lehnen
,
M. F. M.
de Bock
,
S.
Bozhenkov
,
S.
Brezinsek
,
C.
Busch
,
I. G. J.
Classen
,
J. W.
Coenen
,
D.
Harting
,
M.
von Hellermann
,
S.
Jachmich
,
R. J. E.
Jaspers
,
Y.
Kikuchi
,
A.
Krämer-Flecken
,
Y.
Liang
,
M.
Mitri
,
P.
Peleman
,
A.
Pospieszczyk
,
D.
Reiser
,
D.
Reiter
,
U.
Samm
,
D.
Schega
,
O.
Schmitz
,
S.
Soldatov
,
M.
Van Schoor
,
M.
Vergote
,
R. R.
Weynants
,
R.
Wolf
,
O.
Zimmermann
, and
TEXTOR Team
,
Nucl. Fusion
47
(
7
),
522
(
2007
).
11.
O.
Schmitz
,
T. E.
Evans
,
M. E.
Fenstermacher
,
E. A.
Unterberg
,
M. E.
Austin
,
B. D.
Bray
,
N. H.
Brooks
,
H.
Frerichs
,
M.
Groth
,
M. W.
Jakubowski
,
C. J.
Lasnier
,
M.
Lehnen
,
A. W.
Leonard
,
S.
Mordijck
,
R. A.
Moyer
,
T. H.
Osborne
,
D.
Reiter
,
U.
Samm
,
M. J.
Schaffer
,
B.
Unterberg
, and
W. P.
West
(DIII-D and TEXTOR Research Teams)
,
Phys. Rev. Lett.
103
,
165005
(
2009
).
12.
T. E.
Evans
,
Plasma Phys. Controlled Fusion
57
,
123001
(
2015
).
13.
E. C.
da Silva
,
I. L.
Caldas
,
R. L.
Viana
, and
M. A. F.
Sanjuán
,
Phys. Plasmas
9
(
12
),
4917
(
2002
).
14.
M. S. T.
Araújo
,
A.
Vannucci
, and
K. A.
Oliveira
,
Braz. J. Phys.
32
(
1
),
131
(
2002
).
15.
P.
Chrisman
,
J.
Clark
, and
J.
Rome
,
Oak Ridge Report No. ORNL-TM-4501
,
1974
.
16.
E. C.
de Silva
,
I. L.
Caldas
, and
R. L.
Viana
,
IEEE Trans. Plasma Sci.
29
(
4
),
617
(
2001
).
17.
A. C.
Fraile
, Jr.
,
M.
Roberto
, and
I. L.
Caldas
,
Braz. J. Phys.
48
,
426
(
2018
).
18.
A. C.
Fraile
, Jr.
,
M.
Roberto
,
I. L.
Caldas
, and
C. G. L.
Martins
,
IEEE Trans. Plasma Sci.
45
(
11
),
2906
(
2017
).
20.
A.
Vannucci
,
I. C.
Nascimento
, and
I. L.
Caldas
,
Plasma Phys. Controlled Fusion
31
(
2
),
147
(
1989
).
21.
J.
Wesson
,
Tokamaks
, 4th ed. (Oxford University Press,
Oxford
,
2011
).
22.
E. C. D.
Silva
,
I. L.
Caldas
, and
R. L.
Viana
,
Braz. J. Phys.
32
(
1
),
39
(
2002
).
23.
F.
Karger
,
H.
Wobig
,
S.
Corti
,
J.
Gernhard
,
O.
Klüber
,
G.
Lisitano
,
K.
McCormick
,
D.
Meisel
, and
S.
Sesnick
, in
6th International Conference on Plasma Physics and Controlled Nuclear Fusion Research 1974 (Tokyo)
(
IAEA
,
Vienna
,
1975
), Vol.
1
, p.
207
.
24.
A.
Vannucci
,
O. W.
Bender
,
I. L.
Caldas
,
I. H.
Tan
,
I. C.
Nascimento
, and
E. C.
Sanada
,
Nuovo Cim. D
10
(
10
),
1193
(
1988
).
25.
A. S.
Fernandes
and
I. L.
Caldas
, Campos Ressonantes Criados por Correntes Elétricas Helicoidais em Tokamaks (IFUSP, São Paulo,
1985
).
26.
A. S.
Fernandes
,
M. V. A. P.
Heller
, and
I. L.
Caldas
,
Plasma Phys Control Fusion
30
,
1203
(
1988
).
27.
T. E.
Evans
,
R. A.
Moyer
, and
P.
Monat
,
Phys Plasmas
9
,
4957
(
2002
).
28.
H.
Frerichs
,
D.
Reiter
,
O.
Schmitz
,
P.
Cahyna
,
T. E.
Evans
,
Y.
Feng
, and
E.
Nardon
,
Phys Plasmas
19
,
052507
(
2012
).
29.
A. J.
Lichtenberg
and
M. A.
Lieberman
,
Regular and Chaotic Dynamics
(
Springer Science & Business Media
,
New York
,
1992
), Vol.
38
.
30.
J. R.
Dormand
and
P. J.
Prince
,
J. Comput. Appl. Math.
6
,
19
(
1980
).
31.
G.
Corso
,
G.
Oda
, and
I. L.
Caldas
,
Chaos, Soliton Fractals
8
(
12
),
1891
(
1997
).