Proposal of A Linked Mirror Configuration for Magnetic Confinement Experiment Zhichen Feng, Guodong Yu, Peiyou Jiang, and GuoYong Fu∗ Institute for Fusion Theory and Simulation and Department of physics, Zhejiang University, Hangzhou 310027, China Abstract A new linked mirror device for magnetic confinement experiment is proposed. The new linked mirror device consists of two straight magnetic mirrors connected by two half-torus. The structure of the configuration as a whole is three dimensional because the two linear mirror sections are not parallel. The angle between the two mirror sections generates rotational transform which results in good magnetic confinement of toroidally processing passing particles. In this way the usual loss cone of the traditional linear mirror machines is eliminated. The single particle confinement is similar to that of tokamaks with most of particles well confined. The calculated neoclassical confinement is very good and is even better than that of an equivalent tokamak. Potentially the proposed linked mirror configuration can allow MHD stable high beta equilibria with good plasma confinement suitable for neutron sources and magnetic fusion reactors. arXiv:2103.09457v1 [physics.plasm-ph] 17 Mar 2021 ∗corresponding author’s Email: [email protected] 1 I. INTRODUCTION Magnetically controlled fusion energy research has been pursued since early 1950’s. Many types of magnetic confinement devices have been invented, such as magnetic mirrors[1, 2], stellarators[3, 4] and tokamaks. One of earliest and simplest magnetic confinement devices is magnetic mirror. In early days the magnetic mirror devices suffered from two major defects: the fast end loss and strong MHD interchange instability[2]. In particular the end loss of passing particles is very fast leading to poor plasma confinement. Since then various schemes have been developed to mitigate the fast end loss, such as tandem mirror approach[5] and the more modern approach of GDT[6, 7]. These approaches have indeed improved the plasma confinement and MHD stability considerably. However some residual end loss remains and as a result plasma confinement is still not good enough for fusion reactors. In this paper a new linear-torus hybrid device of Linked Mirror Configuration (LMC) has been proposed. The main idea of the new configuration was inspired by the concepts of linear magnetic mirror and the figure-8 stellarator[3]. The hybrid configuration consists of two straight mirror sections connected by two half-torus sections. The two linear mirror sections are not parallel but with a finite angle between them. This angle results in three dimensional (3D) toroidal structure of magnetic field leading to finite rotational transform, just as in the early figure-8 stellarator. Correspondingly the toroidally transiting passing particles are well confined and the end loss of a linear mirror is eliminated. Indeed our simulation results show that the single particle confinement is very good and the neoclassical transport is very low, even lower than that of an equivalent tokamak. The linked mirror approach has been considered previously. The earlier ELMO Bumpy Torus (EBT) consists of many mirrors linked toroidally[8]. Later in 1981 the DRAKON closed fusion trap was proposed by Glagolev et al.[9, 10], consists of two long linear mirror sections joined by two curvilinear sections (CRELs). More recently, Mazzucato proposed a fusion reactor scheme of two linear mirrors joined by two toroidal sections with helical coils[11]. Conceptually our linked mirror approach is similar to those of the DRAKON[9, 10] and Mazzucato’s[11]. The main difference is that in our approach the rotational transform is generated by the angle between the two non-parallel mirror sections rather than the two 3D half-torus sections. As will be shown later, our LMC posses following good properties: (1) there is no end loss; (2) neoclassical confinement is very good and better than that of equivalent tokamaks; (3)All the current coils have the simple geometry of circular planar shape. Finally the device is potentially MHD stable order-of-unity plasma beta. Because of these good properties, the proposed LMC is potentially suitable for neutron sources and fusion reactors. 2 The paper is organized as follows. Section II introduces the basic structure of the LMC. Section III describes the rotational transform properties. Section IV analyses the quasi- omnigeneity and single particle confinement of LMC. Section V presents the numerical results of neoclassical transport. Section VI briefly discusses MHD stability. Section VII gives conclusions. II. BASIC STRUCTURES OF THE LINKED MIRROR CONFIGURATION (a)Top view of the linked mirror configuration (b)Side view of the linked mirror configuration FIG. 1: Basic structure of the linked mirror configuration The basic structure of LMC is shown in Fig. 1 with cartesian coordinates (x,y,z). From the top view (Fig. 1(a)), the configuration shape looks like a race track. The configuration consists of two straight mirror sections and two half-torus sections. However the two straight sections are not parallel as seen from the side view in Fig. 1(b). The angle between the two straight sections is defined as 2Θ with Θ being the angle between one of the straight section and the horizontal plane (z=0 plane). We will show in next section that this angle is responsible for generating rotational transform of the whole configuration. The magnetic 3 0.6 0.4 0.2 0 0.3 z (m) 0.2 -0.2 0.1 -0.4 0 Z (m) -0.1 -0.6 -0.2 -0.3 2.5 3 3.5 4 8.6 8.8 9 9.2 9.4 y (m) x (m) (a)Cross section of magnetic (b)Cross section of magnetic flux surfaces at flux surfaces at φ = 0 plane φ = π/2 plane FIG. 2: Cross section of magnetic flux surfaces of the linked mirror configuration, where φ is the toroidal angle field is provided by simple circular planar coils of equal radius and different currents. For the two mirror sections, the coils forms a circular cylinder. On the other hand, for the half-torus sections, the circular coils are placed in such a way that the centers of the coils form a half circle connecting the two ends of the mirror sections. Obviously, the connections between the straight and half-torus sections can not be tangent, but the magnetic field lines and the magnetic flux surfaces are always smooth, as will be shown later. Besides the angle Θ, other key configuration parameters include L the length of the straight mirror section, d the half distance between the two mirror sections, r the radius of the coils, and Ii the current in each coil. Fig. 1 plots the coils in blue, one of the magnetic surfaces in orange, and a magnetic field line on the flux surface. The parameters used are: Θ = π/20, L = 12m, d = 3m, r = 1m, Itorus/Itube =5I0 :1I0, which means the currents in the coils of the half-torus sections are 5 times of that in the straight mirror tubes. It should be pointed out that the currents in the four end coils in each mirror section are 5I0 too. Fig. 2 plots the cross-section of the magnetic flux surfaces of LMC obtained from field line tracing for the φ = 0 plane (a) and the φ = π/2 plane (b) respectively. Here φ is the toroidal angle about the z-axis through the center of the whole configuration. It shows that by this type of simple connections, nice magnetic flux surfaces are formed with finite rotational transform. The magnetic flux surfaces plotted in red indicate that the flux surfaces go out of the cylindrical tube within coils, as shown in Fig. 3. They are still closed magnetic flux 4 surfaces and can potentially be used as divertors. FIG. 3: The last closed magnetic flux surface with the coils III. ROTATIONAL TRANSFORM PROPERTIES The LMC as just defined above is a special toroidal magnetic confinement configuration with finite external rotational transform generated by the coils. The mechanism for the rotational transform is the 3D magnetic field lines induced by finite Θ just as in the early figure-8 stellarator. In this section, we will show the dependence of the rotational transform on key configuration parameters. 0.17836 0.17834 0.17832 0.1783 0.17828 0.17826 0.17824 0.17822 9.05 9.1 9.15 9.2 x (m) FIG. 4: Rotational transform versus x, where x is the starting point of each field line with z = 0 at φ = 0 plane 5 Fig. 4 plots the rotational transform profile against radius variable x for the LMC plotted in Fig. 1, where x is the major radius at φ = 0 plane. We observe that the rotational transform profile is almost constant with little magnetic shear. 0.6 0.3 0.5 0.25 0.4 0.2 0.3 0.15 0.2 0.1 0.1 0.05 0 0 0 0.1 0.2 0.3 0.4 0.5 6 8 10 12 14 16 2 / L (m) (a)Rotational transform versus Θ (b)Rotational transform versus L 0.3 0.25 0.2 0.15 0.1 0.05 0 2 2.5 3 3.5 4 4.5 5 5.5 6 d (m) (c)Rotational transform versus d FIG. 5: The dependence of rotational transform on the main parameters Fig. 5(a), 5(b) and 5(c) plot the dependence of the rotational transform on Θ, L and d respectively. For each plot, the following parameters are held fixed unless it is varied: Θ = π/20, L = 12m, d = 3m, r = 1m, Itours/Itube = 5I0 : 1I0. Fig.