The Helical Fusion group (Helical Fusion Co., Ltd. in Japan and Helical Fusion USA, Inc.) is developing a 50 MWe-class steady-state helical fusion reactor, which is a kind of stellarator called heliotron composed of two continuous helical coils similar to the large helical device and is operated without plasma current. HESTIA plays the role of the fusion pilot plant. The first-of-a-kind fusion power plant that would follow HESTIA will be a 100 MWe-class steady-state helical fusion reactor. After more than one year of continuous operation, maintenance will be completed within 3 months, aiming for an availability larger than 80%. High-temperature superconducting magnet coils are adopted to generate 8 T of the magnetic field at the helical coil center. HESTIA is a deuterium-tritium fusion reactor, where tritium is self-produced using liquid metal blanket systems. The first wall of the liquid metal blanket including the divertor strike zone is covered and protected by liquid metal free-surface flow, and therefore, individual divertor systems are not required in HESTIA. Electron cyclotron heating is adopted for plasma heating. Since plasma current drive is not required, HESTIA can be operated at a low fusion gain of ∼13, and steady-state operation is principally possible on the order of a year. After a few years of individual development phases, a prototype device is planned to be constructed and operated for the integrated demonstration before starting the construction of HESTIA.

Helical Fusion Co., Ltd. in Japan was established at the end of 2021 aiming for the implementation of a commercial helical fusion reactor for base load electric power to society, and Helical Fusion USA, Inc. is a wholly owned subsidiary of Helical Fusion Co., Ltd. (jointly referred to as the HF group). The HF group is developing a commercial helical fusion reactor based on the design studies of the force free helical reactor (FFHR) series,1 which have similar magnetic configurations to that of the large helical device (LHD),2 incorporating the latest technologies of high-temperature superconducting (HTS) magnets and liquid blankets. Here, we choose the LHD-type heliotron configuration generated by two continuous helical coils because of the capability for steady-state long pulse operation, high density limit, and large apertures for maintenance ports between magnet coils. LHD is the only facility in the world to date that has achieved ion temperature higher than 10 keV,3 electron density higher than 1 × 1021 m−3,4 and plasma duration time longer than 3000 s,5 although these were achieved individually. These achievements are comparable to or better than those of large tokamaks. In particular, the operable density limit and plasma duration time are far superior to those of tokamaks with similar device size and magnetic field strength.6 The necessity of plasma current is the fundamental difference between tokamak and helical. It is not necessary to drive a plasma current inside the helical plasmas since the nested magnetic surfaces confining the high-temperature plasma are completely generated by external magnet coils. Helical plasma is quite stable and does not require sophisticated plasma position control as in tokamaks, even in plasma startup and/or shut-down phases. Helical plasma is free from harmful plasma current disruptions.

The first-of-a-kind (FOAK) fusion power plant envisioned by the HF group is a 100 MWe-class steady-state helical fusion reactor. The electricity output is relatively small compared with conventional fusion reactor designs. This is one of the distinguished points of the HF group's vision. The small-output fusion reactor can be commercialized not only as a small-scale distributed power source but also as a large-scale centralized power source by combining multiple reactors together. The small-output fusion reactor is applicable as, for example, a local power source for factories such as aluminum or titanium smelters, remote islands, or large ships. The construction cost of the FOAK power plant is estimated to be around USD 5 billion. After more than one year of continuous operation, maintenance is to be completed within 3 months to achieve a high availability of more than 80%.

Before proceeding to the construction of the FOAK power plant, a fusion pilot plant (FPP) is planned to be constructed and operated to demonstrate steady-state net-electricity generation. At present, two options named FFHR-HESTIA (or simply, HESTIA, HElical reactor with the surface flow of TIn alloy) and FFHR-SOARHER (or SOARHER, Superconducting Optimally ARranged HElical Reactor) are being considered as candidates for the FPP. HESTIA adopts the modular-type Liquid Metal (LM) blanket7 and flexible HTS WISE (Wound and Impregnated Stacked Elastic tapes) conductors for magnet coils.8 The first wall of the LM blanket including the divertor strike zone is covered and protected by LM free-surface flow, and therefore, individual divertor systems are not required in HESTIA. Electron cyclotron heating (ECH) is adopted for plasma heating. Neutral beam injection (NBI) and ion cyclotron heating (ICH) are not required. By skipping the development of individual divertor systems, NBI, and ICH, which are widely recognized as the most difficult issues in fusion reactor development, a fast-track plan is possible in the case of HESTIA. On the other hand, SOARHER explores a more conservative development strategy by employing either the molten-salt1 or ceramic pebble blanket,9 the tungsten divertor,10,11 and rigid HTS STARS (Stacked Tapes Assembled in Rigid Structure) conductors for magnet coils.11,12 Both HESTIA and SOARHER have a device size twice that of LHD, a magnetic field strength of ∼3 times that of LHD, and a maximum electrical output of 50 MWe.

In this manuscript, we would focus on HESTIA. Section II provides an overview of HESTIA. Section III describes the development plan for HESTIA. Discussions on strategies to secure fusion fuels are given in Sec. IV. Finally, these are summarized in Sec. V. While the development plan discussed here provides a rough sketch of a commercial helical fusion reactor, they are only the envisioned goals that must ultimately be achieved to demonstrate the economic feasibility of fusion. Our current focus is on technical feasibility, i.e., showing that a fusion power plant can indeed produce net electricity.

A schematic view of HESTIA is shown in Fig. 1(a). The diameter of the main body of HESTIA is about 30 m [Fig. 1(b)] and a building of 60 × 160 m2 of floor space, and 100 m height is large enough to install all apparatus needed for reactor operation [Fig. 1(c)]. The helical coil major radius, Rc, is 7.8 m (twice that of LHD), and the magnetic field strength, Bc, at Rc is 8.0 T (∼3 times that of LHD). The helical coil minor radius, ac, is defined by the minor radius of the helical coil current center. In LHD with layer-wound helical coils, ac is variable according to the ratio of the currents through the three layers of helical coils. On the other hand, HESTIA plans to adopt double pancake coil winding, and, therefore, ac is basically determined by the helical coil geometry. The trajectory of the helical coil guiding center in the LHD-type heliotron is described as follows:
(1)
(2)
(3)
where θ, ϕ, m, , and α are the poloidal angle, toroidal angle, toroidal pitch number, poloidal pitch number (the number of helical coils), and the helical coil pitch modulation, respectively. The helical coil pitch parameter, γc, is defined by
(4)
FIG. 1.

(a) Schematic of plasma with two helical coils and four vertical field coils, (b) the main body of 50 MWe-class steady-state helical fusion reactor, HESTIA, and (c) the reactor building layout.

FIG. 1.

(a) Schematic of plasma with two helical coils and four vertical field coils, (b) the main body of 50 MWe-class steady-state helical fusion reactor, HESTIA, and (c) the reactor building layout.

Close modal

In the case of the LHD, m = 10,  = 2, γc is variable in the range of 1.13–1.33 by changing the current ratio in the three helical coil layers, and α = 0.1. In HESTIA-Primary and HESTIA, m and are the same as LHD, but γc = 1.20 and α = 0.0 are adopted as the first candidate of magnetic configuration, because of the larger blanket space and the better confinement property.13,14

HESTIA is a deuterium–tritium (DT) fusion reactor, where tritium is self-produced using liquid metal (LM) blanket systems. The maximum net electrical output is 70 MWe, with the following core plasma parameters: central ion temperature, Ti0, of ∼8 keV, central ion density, ni0, of ∼4 × 1020 m−3, central plasma beta, β0 (ratio of central plasma pressure to magnetic pressure), of ∼3%, and fusion power, Pfusion, of ∼260 MW. ECH is supposed to be applied with the heating power of PECH = 20 MW, and therefore Q ∼13. The HTS magnet coils require a current density, jc, of ∼85 A/mm2. The stored magnetic energy of the HTS magnet system, Wmag, is 66.2 GJ. These and other design parameters of HESTIA are summarized in Table I, together with the reference parameters obtained in LHD and those expected for a prototype device named HESTIA-Primary and a FOAK power plant after HESTIA, which will be discussed in Sec. III.

TABLE I.

Parameters in LHD, HESTIA-Primary, HESTIA, and the FOAK power plant.

LHD (#115772, t = 3.533 s) HESTIA-Primary (prototype) HESTIA (FPP) FOAK power plant
High-T High-nTτ
Rc (m)  3.9  2.73  ←  7.8  ← 
ac (m)  0.936  0.655  ←  1.87  ← 
Bc (T)  0.91  4.0  ←  8.0  6.6 
R0 (m)  3.60  2.49  ←  7.1  ← 
a0 (m)  0.56  0.39  ←  1.13  ← 
B0 (T)  1.0  4.39  ←  8.79  7.25 
α  0.1  0.0  ←  ←  ← 
γ  1.20  ←  ←  ←  ← 
ı2/3  0.79  0.77  0.79  0.77  0.76 
jc (A/mm2 12.4  122  ←  85.1  70.2 
Ic (MA)  1.80  5.46  ←  31.2  25.7 
Wmag (GJ)  0.09  0.67  ←  66.2  45.6 
ne0 (1020 m−3 0.154  0.1  1.41  3.90  3.35 
Te0 (keV)  1.88  30.7  9.81  8.41  10.6 
Ti0 (keV)  N.A.  2.19  5.61  8.36  10.3 
β0 (%)  2.35  0.68  4.5  3.4  5.4 
(Pdep / Pdep1)avg  0.536  0.852  0.853  0.710  0.711 
γDPE  1.32  ←  1.18  ← 
H  1.3  ←  ←  ← 
τE (s)  0.025  0.011  0.073  5.26  2.40 
τEISS95 (s)  0.017  0.0058  0.025  1.57  0.71 
τEISS04v3 (s) (fren = 0.93)  0.024  0.0080  0.038  2.59  1.15 
nedge (1020 m−3 0.50  0.06  0.86  2.38  2.04 
nSudo (1020 m−3 0.59  4.40  4.26  1.21  1.71 
ΓFW,avg (MW/m2 ←  ←  0.65  0.85 
Γdiv,avg (MW/m2 3.4  31  ←  4.3  7.8 
PECH (MW)  0 (PNBI = 6.8)  30  ←  20  31 
Pfusion (MW)  N.A.  ←  260  340 
Q  N.A.  ←  ←  13.1  11.1 
fα  N.A.  ←  ←  ← 
ηα  N.A.  85  ←  ←  ← 
ηEG  N.A.  53  ←  ←  ← 
ηECH  30  30  ←  ←  ← 
ηcryo  1.5  ←  ← 
Pnet (MW)  ←  ←  70.4  103 
Pgross (MW)  ←  ←  139  181 
Qeng  N.A.  ←  ←  2.0  2.3 
Cdirect ($M)  1,000  480  ←  5,000  ← 
τlife (year)  N.A.  ←  ←  6.4  > 5.0 
Pnet τlife (GWh)  N.A.  ←  ←  4,096  4,527 
Cdirect/(Pnet τlife) ($/kWh)  N.A.  ←  ←  1.22  1.10 
LHD (#115772, t = 3.533 s) HESTIA-Primary (prototype) HESTIA (FPP) FOAK power plant
High-T High-nTτ
Rc (m)  3.9  2.73  ←  7.8  ← 
ac (m)  0.936  0.655  ←  1.87  ← 
Bc (T)  0.91  4.0  ←  8.0  6.6 
R0 (m)  3.60  2.49  ←  7.1  ← 
a0 (m)  0.56  0.39  ←  1.13  ← 
B0 (T)  1.0  4.39  ←  8.79  7.25 
α  0.1  0.0  ←  ←  ← 
γ  1.20  ←  ←  ←  ← 
ı2/3  0.79  0.77  0.79  0.77  0.76 
jc (A/mm2 12.4  122  ←  85.1  70.2 
Ic (MA)  1.80  5.46  ←  31.2  25.7 
Wmag (GJ)  0.09  0.67  ←  66.2  45.6 
ne0 (1020 m−3 0.154  0.1  1.41  3.90  3.35 
Te0 (keV)  1.88  30.7  9.81  8.41  10.6 
Ti0 (keV)  N.A.  2.19  5.61  8.36  10.3 
β0 (%)  2.35  0.68  4.5  3.4  5.4 
(Pdep / Pdep1)avg  0.536  0.852  0.853  0.710  0.711 
γDPE  1.32  ←  1.18  ← 
H  1.3  ←  ←  ← 
τE (s)  0.025  0.011  0.073  5.26  2.40 
τEISS95 (s)  0.017  0.0058  0.025  1.57  0.71 
τEISS04v3 (s) (fren = 0.93)  0.024  0.0080  0.038  2.59  1.15 
nedge (1020 m−3 0.50  0.06  0.86  2.38  2.04 
nSudo (1020 m−3 0.59  4.40  4.26  1.21  1.71 
ΓFW,avg (MW/m2 ←  ←  0.65  0.85 
Γdiv,avg (MW/m2 3.4  31  ←  4.3  7.8 
PECH (MW)  0 (PNBI = 6.8)  30  ←  20  31 
Pfusion (MW)  N.A.  ←  260  340 
Q  N.A.  ←  ←  13.1  11.1 
fα  N.A.  ←  ←  ← 
ηα  N.A.  85  ←  ←  ← 
ηEG  N.A.  53  ←  ←  ← 
ηECH  30  30  ←  ←  ← 
ηcryo  1.5  ←  ← 
Pnet (MW)  ←  ←  70.4  103 
Pgross (MW)  ←  ←  139  181 
Qeng  N.A.  ←  ←  2.0  2.3 
Cdirect ($M)  1,000  480  ←  5,000  ← 
τlife (year)  N.A.  ←  ←  6.4  > 5.0 
Pnet τlife (GWh)  N.A.  ←  ←  4,096  4,527 
Cdirect/(Pnet τlife) ($/kWh)  N.A.  ←  ←  1.22  1.10 

In Table I, R0 and a0 are the major radius and the minor radius of the magnetic axis in vacuum, respectively, B0 is the magnetic field strength at R0, ı2/3 is the rotational transform ı = ι/(2π) = 1/q (q is the safety factor in tokamaks) at the normalized minor radius ρ = r/a = 2/3, where a is the toroidally averaged plasma minor radius in the magnetic equilibrium determined by the plasma beta profile, Ic is the helical coil current, ne0 is the electron density at the plasma center, Te0 is the electron temperature at the plasma center, (Pdep/Pdep1)avg is the line-averaged value of the volume-integrated heat deposition profile, Pdep(ρ), normalized by the total heating power, Pdep1 = Pdep(1), γDPE is the confinement improvement factor given by ((Pdep/Pdep1)avg,reactor/(Pdep/Pdep1)avg,exp)0.6 (Ref. 15), H is the additional confinement improvement factor expected in the optimized magnetic field configuration, τE is the energy confinement time, τEISS95 and τEISS04v3 are the energy confinement times predicted by empirical scaling laws of ISS9516 and ISS04 with the renormalization factor fren = 0.93,17 respectively, nedge is the edge electron density at ρ = 1, nSudo is the Sudo density limit,18,19 ΓFW,avg (MW/m2) is the averaged neutron heat flux on the blanket first wall, Γdiv,avg (MW/m2) is the averaged divertor heat flux, PNBI is the NBI heating power used in the reference discharge, fα is the contamination ratio of helium ash, ηα is the ratio of the confined alpha particle, ηEG is the conversion efficiency of electricity generation from the thermal energy absorbed by the blanket, ηECH is the conversion efficiency from the wall-plug electricity to PECH, ηcryo is the conversion efficiency from the wall-plug electricity to the cryogenic power for HTS coils, Pnet is the net electric generation, Pgross is the gross electric generation, Qeng is the engineering fusion gain defined by the ratio of Pnet to total electric power for reactor operation, Cdirect is the direct cost for construction other than land, cite preparation, and licensing, and τlife is the reactor life at full power operation. Note that the Cdirect in Table I is estimated by our original cost evaluation model. This cost model is based as much as possible on the actual performance of LHD and ITER at the end of the 1990s, but it does not take into account inflation and other factors over the past 20 years. In fact, prices in Japan have hardly changed at all over the past 20 years. On the other hand, prices have more than doubled over the past 20 years in other countries. Therefore, if this inflation is taken into account, Cdirect must be modified by a factor of 2 or more. The parameters of HESTIA-Primary, HESTIA, and the FOAK power plant in Table I are estimated using the systems code HELIOSCOPE20 and the integrated transport analysis code TASK3D.21 

To realize HESTIA and the FOAK power plant, however, we need to overcome the technology gaps that remained in the categories of (1) reactor core plasma design, including enhancements in energy confinement in heliotron-type stellarators, (2) HTS magnet, including more flexible HTS conductor suitable for construction of complex coil shapes, (3) LM blanket, including LM that can provide a single solution for both plasma exhaust and tritium breeding, non-ferromagnetic structural metals, and efficient LM circulation pumps, (4) ECH, including high-power, high-frequency gyrotrons, (5) fueling and evacuation, including high-frequency, high-speed pellet injectors, and direct recycling system, and (6) electricity and H2 generation, needed for a novel startup strategy including self-generation of startup tritium using the electricity generated from green H2.

Detailed discussions on reactor core plasma design, HTS magnet, LM blanket, ECH, fueling and evacuation, electricity generation, and H2 generation are given in the Subsections II A–II F.

In stellarator/heliotron devices, including LHD and W7-X,22 the parameter dependence of energy confinement time is known to follow the predictions of a physical model called the gyro-Bohm model, which is also true for tokamaks.23 In the gyro-Bohm model, the energy confinement time is proportional to the 0.6th power of the density, the 0.8th power of the magnetic field strength, and so on. Miyazawa et al. have developed the DPE (direct profile extrapolation) method, which extrapolates the experimentally observed plasma density and temperature profiles in the minor radius direction to the reactor condition using the gyro-Bohm property.24,25 This enables highly accurate prediction of fusion power in fusion reactor core plasmas.13,26

In Fig. 2, radial profiles in HETSIA estimated using the DPE method are shown together with the reference profiles of LHD from which the extrapolation was made. The reference data were taken from an NBI-heated hydrogen discharge of #115772 (t = 3.533 s). Radial profiles of electron density, electron temperature, and plasma beta are shown in Figs. 2(a)–2(c), respectively, where the factors fn, fT, and fβ shown on the right-top of each panel denote the ratio of the corresponding parameters in HESTIA to those in the reference data of LHD. The NBI power deposition profile in LHD, pNBILHD(ρ), is shown in Fig. 2(d). In Fig. 2(e), shown are the profiles of ECH power deposition, pECH(ρ), the alpha heating power deposition, pα(ρ), the bremsstrahlung loss, pB(ρ), and the total heat deposition, ptot(ρ) = pECH(ρ) + pα(ρ) − pB(ρ), in HESTIA. The neoclassical heat flux, QtotneoS = QeneoS + QineoS, where subscripts e and i denote the contributions of electrons and ions, respectively, passing through the magnetic surface of an area S(ρ) = dV(ρ)/dρ at each ρ, where V(ρ) is the volume inside ρ, and the volume-integrated total heating power, Ptot(ρ), given by
(5)
is shown in Fig. 2(f).
FIG. 2.

Radial profiles of (a) electron density, (b) electron temperature, (c) plasma beta in HESTIA (solid lines) and the reference data of LHD (open circles with smoothed broken lines). The NBI power deposition profile in the reference data of LHD is shown in (d). Profiles of ECH power deposition (open squares), alpha heating power deposition (open circles), bremsstrahlung loss (open diamonds), and the total heating power deposition (closed circles) in HESTIA are shown in (e). In (f), the volume integrated total heating power (bold solid line) is compared with the neoclassical heat flux (bold broken line) that is the sum of the neoclassical heat fluxes in the electron (thin broken line) and ion (thin solid line) channels, respectively.

FIG. 2.

Radial profiles of (a) electron density, (b) electron temperature, (c) plasma beta in HESTIA (solid lines) and the reference data of LHD (open circles with smoothed broken lines). The NBI power deposition profile in the reference data of LHD is shown in (d). Profiles of ECH power deposition (open squares), alpha heating power deposition (open circles), bremsstrahlung loss (open diamonds), and the total heating power deposition (closed circles) in HESTIA are shown in (e). In (f), the volume integrated total heating power (bold solid line) is compared with the neoclassical heat flux (bold broken line) that is the sum of the neoclassical heat fluxes in the electron (thin broken line) and ion (thin solid line) channels, respectively.

Close modal

The Sudo density limit18 for the case of HESTIA is also shown in Fig. 2(a). Note that nSudo corresponds to the edge density limit.19 Although nedge exceeds nSudo in HESTIA, as in Table I, the difference seems to be not critical in Fig. 2(a). If this difference between nedge and nSudo appears to be critical, it will become necessary to add an off-axis ECH to increase nSudo without largely modifying the beta profile. It should be noted that the beta profile changes if the heating power is increased, and the magnetic axis moves outward due to the Shafranov shift. In the LHD-type heliotron configuration, the outward shift of the magnetic axis results in the increase in the neoclassical transport. Therefore, the off-axis ECH that has a weaker impact on the beta profile compared with the on-axis ECH will be favorable.

In designing the core plasma parameters of HESTIA shown in Table I, we assumed a modest confinement improvement factor of H = 1.3, together with the heating profile effect of ECH and the alpha heating profile localized at the plasma center. It has been observed in LHD that the energy confinement apparently improves in the cases with central heating compared with the cases with broad or hollow heating profiles. This heating profile effect has been already formulated and incorporated into the DPE method using (Pdep/Pdep1)avg.15 As can be seen in Fig. 2(e), the heat deposition profile in HESTIA is extremely peaked at the plasma center. This results in a large (Pdep/Pdep1)avg of 0.710, which is larger than that of the reference data of 0.536. The confinement improvement due to this center-peaked heating profile is taken into account using a factor of γDPE = (0.710/0.536)0.6 ∼1.18. Although the extremely center-peaked heating profile has a possibility to trigger the formation of an internal transport barrier (ITB) that is not taken into account at this moment, because the physics scenario of ITB formation is still unknown.

Instead of assuming confinement improvement due to the ITB formation, we are now trying to find out the optimized magnetic configuration to achieve an additional confinement improvement of H = 1.3, other than γDPE, by slightly deforming the helical coils from those of LHD.27 It is known that the inward-shifted magnetic configuration with R0 ∼ 3.6 m in LHD conforms to a σ-optimized field,28 where the magnetic field strength at helical ripple bottoms is similar along a field line.29,30 However, there is still ample potential for optimizing the magnetic configuration of LHD further,14 because we have gained experimental and theoretical knowledge of plasma physics and have established numerical analysis tools during the 25 years of LHD operation. We are now searching for an optimized magnetic configuration that can achieve high MHD stability and reduced neoclassical transport at the same time. Although it is expected that the neoclassical heat flux of QtotneoS can be compensated with Ptot as shown in Fig. 2(f), it is marginal at ρ ∼ 0.5. Therefore, it is favorable to reduce neoclassical transport further to secure a margin for turbulent transport. If the neoclassical transport cannot be reduced and the turbulent transport is not negligible at ρ ∼ 0.5, additional heating power should be added inside ρ < 0.5 to compensate for the turbulent heat flux. The off-axis ECH without largely modifying the beta profile will be favorable for this, as was explained already.

As for the high-energy particle confinement, we assume ηα = 85%, according to the simulation results in the former study,13 in which it is shown that the ratio of confined alpha particles in the case of the magnetic configuration similar to HESTIA is larger than 85% at ρ < 0.7, as long as alpha particle loss boundary is set on the blanket first wall, and the reentering effect is taken into account. Note that pα at ρ > 0.7 is negligibly small as shown in Fig. 2(e). It is expected that the high-energy particle confinement property can be also improved in the optimized magnetic configuration with reduced neoclassical transport.

Our reactor core plasma design does not rely on the non-thermal ion population. Instead, we plan to use only electron heating by ECH and alpha heating. The so-called energy equipartition by Coulomb collision between electrons and ions will take place, and the ion temperature is expected to become similar to the electron temperature. HESTIA and the FOAK power plant will be operated at high density as shown in Table I, which promotes electron-ion collisions. The equipartition effect is calculated using a theoretical model and taken into account in our reactor core plasma design.31 

Although LHD is the world's first large superconducting helical device and our helical reactor design largely relies on the knowledge accumulated in LHD, we need further R&D. Especially, there is a large difference that LTS (low-temperature superconducting) magnet coils cooled to 4 K by liquid helium that are adopted in LHD, while HTS magnet coils cooled to 20 K by helium gas that are adopted in our helical reactors. There are several reasons why HTS magnets must be used. First, the critical magnetic field is higher with HTS magnets than with LTS magnets. Second, the coil current density can be higher than that of LTS magnets. A high coil current density is essential to ensure a large blanket space. Third, the amount of helium required as a coolant can be significantly reduced compared to the LTS magnets. In recent years, the global supply of helium has been depleted. In our helical reactors, gas helium, of which the density at 20 K is 1/50 of liquid helium at 4 K, will be used to cool superconducting magnets to address this problem. The 20 K helium gas will be produced using a helium refrigerator, and the inlet temperature could be precisely kept at 20 K using the heat exchange with liquid hydrogen that is generated by solar power as will be discussed in Subsection II F. A cooling efficiency of 2% is assumed for the cryogenic system for HTS magnets in HESTIA. The cooling efficiency is supposed to be increased from 1.5% to 2% in the step-by-step development from the prototype to HESTIA.

The maximum magnetic field strength on the magnet coils is expected to increase from 6 T in LHD to ∼20 T in HESTIA. The high current density of 85 A/mm2 is also required for the HTS conductor of HESTIA, which is higher than that of LHD of 34 A/mm2 at the maximum. The HF group is developing a new HTS WISE conductor8 using stacked REBCO tapes. The WISE conductor is quite flexible to be wound to fabricate the HTS coils with complicated shapes, such as the helical coils. After the winding procedure is completed, the volume and the WISE conductor inside the coil body are impregnated by a low-melting point alloy to fix and reinforce the WISE conductor. It is supposed to fabricate the HTS helical coils for HESTIA without insulation. Such a no-insulation HTS coil is expected to have a high resilience to the quench accident, although the time constant of the coil increases.8,32 If a normal transition takes place locally in the no-insulation HTS coil, the coil current flowing the normal part is automatically bypassed to the neighboring HTS tapes which are still superconducting, and further progress of normal transition can be avoided.

A schematic view of the LM blanket in HESTIA is shown in Fig. 3(a). The LM blanket is composed of nine main modules as shown in Fig. 3(b) and one test blanket module at each of ten similar sections.7 Our LM blanket is designed aiming at easy construction and maintenance. Although 90 blanket modules in total are used in HESTIA or the FOAK power plant, it is possible to replace them within 1 month, if three or four modules are replaced in one day. Because fusion reactors have high radiation levels immediately after the operation, access to the inside of the reactor is difficult. All LM blanket cartridges in HESTIA are designed to be maintained from the upper port using a three-point suspension crane. It should be noted that no pipe cutting or welding is necessary for our LM blanket.7 These eliminate the need for a robot to enter the reactor during maintenance. After 1 year or more of continuous operation, a maintenance period of about 3 months is supposed. Then, the availability is larger than 80%. The maintenance period consists of 1 month for decreasing the magnetic field and for cooling the radioactivity and tritium degassing of blanket materials, 1 month for inspection and replacement of the blanket as a necessity, and 1 month for increasing the magnetic field and for preheating of blanket modules for the next operation.

FIG. 3.

Schematics of (a) the LM blanket with the first wall covered by LM free-surface flow, (b) nine types of the main LM blanket modules in one section, and (c) the pressurized-gas-driven LM circulation system.

FIG. 3.

Schematics of (a) the LM blanket with the first wall covered by LM free-surface flow, (b) nine types of the main LM blanket modules in one section, and (c) the pressurized-gas-driven LM circulation system.

Close modal

The blanket first wall is surely one of the key components in magnetic fusion reactors. It is reported from the long pulse discharge experiment in LHD that the particle recycling changes with a slow time constant of a few tens of minutes.33 The slow variation of particle recycling can cause difficulty in density control in a fusion reactor. In HESTIA, this problem is solved by adopting the LM first wall, i.e., the first wall of our blanket is covered and protected by the LM free-surface flow. Since the first wall is always refreshed by flowing LM, the hydrogen recycling condition is expected to be kept stable. To implement the LM first wall, the plasma-facing sides of each blanket module are composed of porous material. The LM leaks out from the porous material and forms an LM free-surface flow on the first wall.

Since the free surface of the LM is exposed to the plasma, the vapor pressure of the LM must be low enough not to affect the plasma performance. Tin is one of the low melting point materials and has a very low vapor pressure.34 Therefore, we plan to develop a tin–lead–lithium alloy LM breeder material for fusion reactor blankets, which is composed mainly of tin mixed with lead for neutron multiplication and lithium for tritium production. Here, lead can be substituted by other materials if only it results in a better tritium breeding ratio (TBR) without largely affecting the LM characteristics such as the melting point. Since the hydrogen retention in the pure tin is as low as ∼0.1 at.%,35–39 hydrogen isotopes including the DT molecules hitting the LM and T atoms generated inside the LM will be immediately emitted to the vacuum vessel (VV) and exhausted by the cryopumps. Hydrogen retention in the tin-lead-lithium alloy, which may differ from pure tin, will be investigated in future studies. Tin is highly corrosive to steel and therefore special consideration for corrosion protection is required for structural materials and the porous material for the blanket of HESTIA. Research on corrosion protection is also considered in our research plan.

Contamination of the main plasma by heavy impurities such as tin and lead should be avoided. At this point, we expect the impurity shielding effect of the ergodic layer surrounding the last-closed-flux-surface (LCFS) in the LHD-type heliotron configuration. In LHD, it has been observed that the friction forces inside the ergodic layer that push the impurity ions in the direction from LCFS to divertor become dominant over the thermal forces that work vice versa.40 It was also reported that this impurity shielding effect of the ergodic layer is enhanced at higher density and for heavier impurities.

On the other hand, we are also planning to develop new structural material. Although RAFM (reduced activation ferritic martensitic) steels are widely adopted in the fusion community as the structural material for fusion reactors, they are ferromagnetic. In magnetic confinement fusion devices, the use of ferromagnetic materials that produce error fields should be avoided. Therefore, we are going to develop a new non-magnetic material. One of the possible candidates is a non-magnetic reduced activation high-manganese steel. In the past, high-manganese steel was considered for use in fusion reactors, but the first choice was given to RAFM because the high-manganese steel has problems of large decay heat and swelling caused by neutron irradiation.41 These weaknesses of high-manganese steel, however, can be greatly mitigated in low-power fusion reactors with fusion power of several tens of MW, such as the HESTIA and the FOAK power plant we are envisioning.

Additive manufacturing or other 3D printing method using pure titanium, titanium alloy, or the new high-manganese steel will be adopted to fabricate the porous material equipped on the blanket's first wall. In the case of titanium porous material, its surface will be already covered by a strong titanium oxide layer that will be effective for corrosion protection. In the case of high-manganese porous material, it will be necessary to form the anti-corrosion oxide layer on its surface. In both cases, the first wall side of the porous material should have a high wettability to form the LM free-surface flow. Research on wettability is also planned in our research plan.

Independent divertor devices are not needed in HESTIA or the FOAK power plant of the HESTIA-type, because the divertor strike zones of the blanket first wall are already covered and protected by the LM surface flow. No special heat exchange equipment for the divertor is needed either since all plasma energy reaching the divertor is absorbed by the LM of the blanket. The averaged divertor heat load in HESTIA is expected to be ∼4 MW/m2. This can be easily managed by the LM free-surface flow while keeping the metal vapor pressure low enough. These factors have contributed greatly to the simplification of the reactor system and the reduction of the development program, construction time, and total cost.

We also need a reliable LM circulation system. Because the mechanical pump with high-rotation speed has low reliability, and the electromagnetic pump has only poor driving efficiency, we will newly develop an original reciprocating pump driven by pressurized gas as shown in Fig. 3(c). This system is composed of multiple reservoir tanks containing LM. High-pressure gas is alternately introduced to each tank to push the LM out to HESTIA. The LM is introduced to blanket modules through the top of the VV. The LM is heated inside the blanket modules and exhausted to the lower pool. The high-temperature LM then flows into the heat exchanger and finally returns to the reservoir tanks again. It should be noted that the electric power needed for the new LM circulation system is quite unknown at this moment. In Table I, the electricity for LM circulation is estimated assuming mechanical pumps. We are expecting that our new LM circulation system can be operated with smaller electricity compared with mechanical or electromagnetic pumps. The new LM circulation system will be applied to the prototype and then to HESTIA.

ECH is adopted in HESTIA as the only auxiliary heating method to enable a year-long continuous operation without maintenance. As shown in Table I, PECH = 20 MW is required in HESTIA. The resonance frequency of ECH in HESTIA and the FOAK power plant is as high as 250 GHz, due to the high magnetic field of around 9 T at the plasma center. It becomes difficult to design a gyrotron as the oscillation frequency increases. In LHD, 154 GHz gyrotrons were developed and applied in the experiment. The 154 GHz gyrotron could put 1 MW of ECH power for only a short period of shorter than 2 s and the output power becomes limited to lower than 0.5 MW at continuous operation. No gyrotron working at 250 GHz with 1 MW of output power at CW (Continuous Wave) is available at this moment. Therefore, we will develop a 250 GHz-1MW-CW gyrotron by ourselves. Although it is designed to operate continuously at 1 MW, we will prepare 60 gyrotrons of 250 GHz-1MW-CW for HESTIA to ease the requirements for the gyrotron. These gyrotrons will be operated alternately to inject 20 MW of ECH power continuously. The 60 gyrotrons are evenly distributed at ten outer ports of VV, six at a port, as shown in Fig. 4. ECH is vertically injected from the outside of the horizontal cross section of the helical plasma.

FIG. 4.

Schematics of (a) all ten ECH systems, and (b) the close-up view of two ECH systems, where each system includes six 1 MW-CW gyrotrons. Waveguides are connected to gyrotrons one by one.

FIG. 4.

Schematics of (a) all ten ECH systems, and (b) the close-up view of two ECH systems, where each system includes six 1 MW-CW gyrotrons. Waveguides are connected to gyrotrons one by one.

Close modal

The 250 GHz-1MW-CW gyrotrons we develop would have a high wall-plug efficiency of 50%, i.e., 40 MW of the wall-plug electricity is required to inject 20 MW of ECH power into HESTIA. Although the wall-plug efficiency of ECH in LHD is around 30%, it is possible to increase the efficiency to larger than 50% depending on the design of the gyrotron.

As for the fueling apparatus, we will develop 30-barrel pipe-gun DT-ice Pellet Injectors (PIs) together with a direct fuel gas recycling system.42,43 It will be required in HESTIA to inject high-speed DT pellets at a few Hz of injection frequency. To obtain high-speed pellets, the pipe-gun pellet injector is the best candidate among the various pellet injection systems. On the other hand, the difficulty of pipe-gun pellet injectors is that it requires at least one minute to form a pellet in the pipe. This is the reason why we need to develop 30-barrel PIs. We will set ten PIs at all ten outer ports of HESTIA as shown in Fig. 5, i.e., 300 barrels in total. If a pellet is formed within 100 s, then we can inject pellets at an injection frequency higher than 3 Hz. In LHD, a 20-barrel PIs has been developed and used in the experiment. We are aiming to increase the barrel number from 20 to 30.

FIG. 5.

Schematics of (a) the pellet injection system incorporated with the direct fuel gas recycling system consisting of cryopumps and DT gas conditioning systems and (b) the close-up view of a 30-burrel pipe-gun DT-ice pellet injector.

FIG. 5.

Schematics of (a) the pellet injection system incorporated with the direct fuel gas recycling system consisting of cryopumps and DT gas conditioning systems and (b) the close-up view of a 30-burrel pipe-gun DT-ice pellet injector.

Close modal

The PIs are to be coupled with the direct fuel gas recycling system, where the DT gas exhausted from the VV is directly supplied to the pellet injectors after being purified and pressurized in the DT gas conditioning system. By this, it becomes possible to suppress the tritium site inventory. Multiple cryopumps as shown in Figs. 1(b) and 5 will be mainly used for vacuum pumping in HESTIA. Cryopumps are accumulating vacuum pumps and must be periodically regenerated to expel the accumulated hydrogen isotopes. Steady-state operation is achieved by alternately repeating pumping and regeneration operations using multiple cryopumps.

In HESTIA, a supercritical CO2 (SCO2) gas turbine generator44 will be adopted to achieve an electricity generation efficiency higher than 50%. The working temperature will be 800–1200 K. High energy conversion efficiency and relatively small size compared with ordinary steam turbines are the merits of SCO2 gas turbines. We will develop an SCO2 gas turbine power generation facility for HESTIA.

Large electricity is required at the startup phase of a fusion reactor, to cool the superconducting magnets and to preheat the blanket modules. We will use H2 produced by solar power and an H2 gas turbine generator to provide this startup power. We will also develop a solar-powered H2 production facility and an H2 gas turbine power generation facility. These systems can be also applied to HESTIA-Primary which is planned to operate without DT fusion reaction or power generation. The volume of H2 can be reduced by ∼1/730 times by liquifying the H2 gas. Furthermore, the cold heat of the 20 K liquid H2 is available to stabilize the inlet temperature of the HTS magnets.

HESTIA is supposed to operate in a steady state for 1 year. As shown in Table I, the fusion output of HESTIA is 260 MW. The gross electric power generated from this fusion output will be ∼140 MW. Subtracting a circulating electric power required for reactor operation, the net electricity of 70 MW can be supplied to the power grid or electricity consumers. The majority of the circulating electric power is used in the cryogenic system for HTS magnet cooling and the ECH system for plasma heating.

As discussed in the Subsection II A, the HF group will adopt an optimized magnetic configuration. To verify the improved plasma performance, it is planned to construct and operate a prototype reactor named HESTIA-Primary. It is possible to examine if the energy confinement is improved or not, without the DT fusion reaction. Therefore, HESTIA-Primary will be operated in the non-nuclear condition, where only H2 gas is used without D or T. If the parameters listed in Table I for HESTIA-Primary, which are estimated using the DPE method with the confinement improvement factor of 1.3, are achieved, that can be strong evidence showing that the improved energy confinement is achieved in the optimized magnetic configuration. Then, we will apply these results to HESTIA, and the parameters listed in Table I for HESTIA that satisfy the FPP requirements will surely be achieved with high prediction accuracy.

The HF group aims to demonstrate the functionality of the key elements, such as the HTS magnet and LM blanket, by the early 2020s, and to integrate them into HESTIA-Primary, which is expected to be constructed within roughly 3 years and to begin operation by the end of the 2020s. This is a fairly difficult schedule, but not necessarily an impossible one, given the budget and human resources. Of course, if the necessary budget cannot be secured, the schedule will be pushed back. HESTIA-Primary will be about 1/3 the size of HESTIA to minimize construction costs as low as USD 480 million. As was discussed in Sec. II, this cost evaluation is based on the prices at the end of the 1990s, and if inflation during these 20 years is taken into account, the construction cost must be modified by a factor of 2 or more. The size of HESTIA-Primary is large enough to simultaneously demonstrate the functionality of the key elements and the high plasma performance. The construction and operation of HESTIA-Primary as a prototype before HESTIA can reduce the development risk.

After the successful operation of HESTIA-Primary, HESTIA will be constructed. The construction period is expected to be 5–6 years. Therefore, in the fastest case, we can start the operation of HESTIA by the early 2030s. In HESTIA, a net electricity generation of 50 MWe at the maximum and a steady-state operation for more than 1 year will be demonstrated. Thereafter, maintenance will be completed within 3 months, and the next operation will be resumed to demonstrate the availability of more than 80%. Our FOAK fusion plant assumes net electrical output of 100 MWe in a device of the same size as HESTIA. The thermal/neutron load handling of the in-vessel components and the energy extraction from the reactor core can be realized by similar systems as HESTIA. On the other hand, the neutron flux to the HTS coils will increase. Therefore, it is necessary to improve the neutron shielding performance by, for example, expanding the blanket space by further optimization of the magnetic field configuration and/or developing efficient neutron shielding materials, to achieve a sufficient plant lifetime from the viewpoint of plant economy. As listed in Table I, the central beta of the plasma needs to be increased from 3.4% in HESTIA to 5.4% in the FOAK power plant to increase the fusion output from 260 to 340 MW and the net electricity from 70 to 103 MW. A high central beta larger than 5% must be demonstrated in HESTIA even for a short time of, for example, 3 h, which is short enough not to influence the HESTIA's life significantly.

Based on the knowledge obtained in HESTIA and further progresses in the magnetic configuration optimization and material studies, the FOAK power plants of 100 MWe net electricity will become available by the late 2030s to the early 2040s.

HESTIA is a DT fusion reactor and strategies to secure the initial tritium inventory at startup and enriched 6Li supplies for reactor operation are quite important. Our strategies for these issues include a so-called Deuterium–Deuterium (DD) startup operation and a precise 3D neutron transport simulation.

The initial tritium inventory is generated by the HESTIA itself by the DD startup operation, where the DD fusion reaction is used to produce fast neutrons to produce tritium in the blanket via a nuclear spallation reaction of 6Li atoms by 2.45 MeV DD neutron. The ratio of D at the very initial stage of the DD startup for the first operation of HESTIA is 100%. The ratio of T gradually increases with reactor operation. It is possible to shorten the period of DD startup operation by mixing the tritium atoms produced in the blanket into the fuel gas. After the DT ratio reaches 50:50, the excess T is separated, recovered, and stored for the second and subsequent operations of HESTIA or for the FOAK power plant. A cryogenic distillation system similar to that of ITER45 will be adopted for tritium separation. Although detailed and accurate estimations of the DD startup scenario will be performed in the very early phase of our development plan, it is roughly estimated that a few months of DD operation is enough to produce the tritium inventory needed for a few hours of DT operation.

The use of isotope-enriched lithium is known to be effective to increase the TBR of the blanket in a fusion reactor. Also in our case, the isotope-enriched lithium that includes 80 at. % 6Li and 20 at. % 7Li is supposed. The accurate enrichment factor of 6Li from the natural isotope ratio of ∼8 at. % 6Li and ∼92 at. % 7Li will be determined at the initial stage of our project by the 3D neutron transport simulation using an accurate 3D-CAD model of HESTIA.

The FPP concept we would propose is a 50 MWe-class steady-state helical fusion reactor named HESTIA. The major radius of helical coils is ∼8 m, and the magnetic field strength at the plasma center is ∼9 T. HESTIA is a DT fusion reactor, where tritium is self-produced using LM blanket systems. Since plasma current drive power is unnecessary, HESTIA can be operated at a low fusion gain of ∼13, and steady-state operation is principally possible on the order of a year. The construction cost of HESTIA is estimated to be USD 5 billion. The availability will be larger than 80%.

HF group has started research and development on magnetic configuration optimization, HTS conductor, LM blanket equipped with the LM first wall, 250 GHz-1MW-CW gyrotrons, 30-burrel pipe-gun DT-ice pellet injector, SCO2 gas turbine generator, solar-powered H2 generation system, H2 gas turbine generator, and HTS magnet cooling system using liquid H2. These will be integrated into a prototype device named HESTIA-Primary. After the integrated technology demonstration by HESTIA-Primary, HESTIA will be constructed to demonstrate 70 MWe of net electricity generation by the fusion reaction. A long pulse operation for more than 1 year and a demonstration of high availability of > 80% will be also conducted in HESTIA. Finally, in the fastest case, the FOAK power plants of 100 MWe net electricity will become available by the late 2030s.

The parameters, technological issues, and schedules discussed here are our targets that must ultimately be achieved to demonstrate the economic feasibility of fusion. Our current focus is to demonstrate the technical feasibility that a fusion power plant can indeed produce net electricity.

The author is grateful to members of the HF group (T. Taguchi, N. Yanagi, T. Shimozuma, A. Nishimura, and A. Sagara) and collaborators in the United States (D. Andruczyk (Univeristy of Illinois), [A. Khodak and R. Maingi (PPPL), C. Kessel, S. Smolentsev, B. Pint, L. Baylor (ORNL), and R. Yamaguchi (PARC)] and in Japan (Y. Sugawara (Aoyama Gakuin University), M. Sakama (Tokushima University), R. Kasada (Tohoku University), Y. Narushima, Y. Hamaji, and G. Motojima (NIFS).

The authors have no conflicts to disclose.

Junichi Miyazawa: Conceptualization (lead); Data curation (equal); Formal analysis (equal); Funding acquisition (lead); Investigation (lead); Methodology (equal); Project administration (lead); Resources (equal); Software (supporting); Supervision (lead); Validation (equal); Visualization (equal); Writing – original draft (lead); Writing – review & editing (lead). Takuya Goto: Data curation (equal); Formal analysis (equal).

The data that support the findings of this study are available within the article.

1.
A.
Sagara
,
J.
Miyazawa
,
H.
Tamura
,
T.
Tanaka
,
T.
Goto
,
N.
Yanagi
,
R.
Sakamoto
,
S.
Masuzaki
, and
H.
Ohtani
,
FFHR Design Group
.
Nucl. Fusion
57
,
086046
(
2017
).
2.
A.
Komori
,
H.
Yamada
,
S.
Imagawa
,
O.
Kaneko
,
K.
Kawahata
,
K.
Mutoh
,
N.
Ohyabu
,
Y.
Takeiri
,
K.
Ida
,
T.
Mito
et al,
Fusion Sci. Technol.
58
,
1
11
(
2010
).
3.
M.
Osakabe
,
H.
Takahashi
,
H.
Yamada
,
K.
Tanaka
,
T.
Kobayashi
,
K.
Ida
,
S.
Ohdachi
,
J.
Varela
,
K.
Ogawa
,
M.
Kobayashi
et al,
Nucl. Fusion
62
,
042019
(
2022
).
4.
J.
Miyazawa
,
S.
Masuzaki
,
R.
Sakamoto
,
B. J.
Peterson
,
N.
Tamura
,
M.
Goto
,
M.
Kobayashi
,
M.
Shoji
,
T.
Akiyama
,
H.
Yamada
, and
LHD Experiment Group
.
Fusion Sci. Technol.
58
,
200
(
2010
).
5.
T.
Mutoh
,
R.
Kumazawa
,
T.
Seki
,
K.
Saito
,
H.
Kasahara
,
Y.
Nakamura
,
S.
Masuzaki
,
S.
Kubo
,
Y.
Takeiri
,
T.
Shimozuma
et al,
Nucl. Fusion
47
,
1250
(
2007
).
6.
J.
Miyazawa
,
S.
Masuzaki
,
H.
Yamada
,
R.
Sakamoto
,
B. J.
Peterson
,
M.
Shoji
,
N.
Ohyabu
,
A.
Komori
,
O.
Motojima
, et al,
Fusion Sci. Technol.
50
,
192
(
2006
).
7.
J.
Miyazawa
,
T.
Goto
,
Y.
Hamaji
, and
M. I.
Kobayashi
,
Nucl. Fusion
61
,
126062
(
2021
).
8.
S.
Matsunaga
,
Y.
Narushima
,
Y.
Onodera
,
Y.
Terazaki
,
J.
Miyazawa
, and
N.
Yanagi
,
IEEE Trans. Appl. Supercond.
30
,
4601405
(
2020
).
9.
K.
Tobita
,
R.
Hiwatari
,
Y.
Sakamoto
,
Y.
Someya
,
N.
Asakura
,
H.
Utoh
,
Y.
Miyoshi
,
S.
Tokunaga
,
Y.
Homma
,
S.
Kakudate
et al,
Fusion Sci. Technol.
75
,
372
(
2019
).
10.
T.
Hirai
,
S.
Panayotis
,
V.
Barabash
,
C.
Amzallag
,
F.
Escourbiac
,
A.
Durocher
,
M.
Merola
,
J.
Linke
,
T.
Loewenhoff
,
G.
Pintsuk
et al,
Nucl. Mater. Energy
9
,
616
(
2016
).
11.
N.
Yanagi
,
T.
Goto
,
J.
Miyazawa
,
H.
Tamura
,
Y.
Terazaki
,
S.
Ito
,
H.
Hashizume
, and
A.
Sagara
,
J. Fusion Energy
38
,
147
(
2019
).
12.
N.
Yanagi
,
Y.
Terazaki
,
S.
Ito
,
H.
Tamura
,
S.
Hamaguchi
,
T.
Mito
,
H.
Hashizume
, and
A.
Sagara
,
Cryogenics
80
,
243
(
2016
).
13.
J.
Miyazawa
,
Y.
Suzuki
,
S.
Satake
,
R.
Seki
,
Y.
Masaoka
,
S.
Murakami
,
M.
Yokoyama
,
Y.
Narushima
,
M.
Nunami
,
T.
Goto
et al,
Nucl. Fusion
54
,
043010
(
2014
).
14.
T.
Goto
,
K.
Ichiguchi
,
H.
Tamura
,
J.
Miyazawa
,
S.
Satake
,
H.
Yamaguchi
, and
N.
Yanagi
,
FFHR Design Group
.
Plasma Fusion Res.
16
,
1405085
(
2021
).
15.
J.
Miyazawa
,
T.
Goto
,
R.
Sakamoto
,
G.
Motojima
,
C.
Suzuki
,
H.
Funaba
,
T.
Morisaki
,
S.
Masuzaki
,
I.
Yamada
,
S.
Murakami
et al,
Nucl. Fusion
52
,
123007
(
2012
).
16.
U.
Stroth
,
M.
Murakami
,
R. A.
Dory
,
H.
Yamada
,
S.
Okamura
,
F.
Sano
, and
T.
Obiki
,
Nucl. Fusion
36
,
1063
(
1996
).
17.
H.
Yamada
,
J. H.
Harris
,
A.
Dinklage
,
E.
Ascasibar
,
F.
Sano
,
S.
Okamura
,
J.
Talmadge
,
U.
Stroth
,
A.
Kus
,
S.
Murakami
,
M.
Yokoyama
,
C. D.
Beidler
,
V.
Tribaldos
,
K. Y.
Watanabe
, and
Y.
Suzuki
,
Nucl. Fusion
45
,
1684
(
2005
).
18.
S.
Sudo
,
Y.
Takeiri
,
H.
Zushi
,
F.
Sano
,
K.
Itoh
,
K.
Kondo
, and
A.
Iiyoshi
,
Nucl. Fusion
30
,
11
(
1990
).
19.
J.
Miyazawa
,
R.
Sakamoto
,
S.
Masuzaki
,
B. J.
Peterson
,
N.
Tamura
,
M.
Goto
,
I.
Yamada
,
K.
Narihara
,
K.
Tanaka
,
T.
Tokuzawa
et al,
Nucl. Fusion
48
,
015003
(
2008
).
20.
T.
Goto
,
Y.
Suzuki
,
N.
Yanagi
,
K. Y.
Watanabe
,
S.
Imagawa
, and
A.
Sagara
,
Nucl. Fusion
51
,
083045
(
2011
).
21.
T.
Goto
,
J.
Miyazawa
,
R.
Sakamoto
,
R.
Seki
,
C.
Suzuki
,
M.
Yokoyama
,
S.
Satake
, and
A.
Sagara
,
FFHR Design Group
.
Nucl. Fusion
55
,
063040
(
2015
).
22.
R. C.
Wolf
,
A.
Ali
,
A.
Alonso
,
J.
Baldzuhn
,
C.
Beidler
,
M.
Beurskens
,
C.
Biedermann
,
H.-S.
Bosch
,
S.
Bozhenkov
,
R.
Brakel
et al,
Nucl. Fusion
57
,
102020
(
2017
).
23.
U.
Stroth
,
Plasma Phys. Controlled Fusion
40
,
9
(
1998
).
24.
J.
Miyazawa
,
T.
Goto
,
T.
Morisaki
,
M.
Goto
,
R.
Sakamoto
,
G.
Motojima
,
B. J.
Peterson
,
C.
Suzuki
,
K.
Ida
,
H.
Yamada
et al,
Fusion Eng. Des.
86
,
2879
(
2011
).
25.
J.
Miyazawa
,
T.
Goto
,
R.
Sakamoto
,
A.
Sagara
, and
FFHR Design Group
.
Nucl. Fusion
54
,
013014
(
2014
).
26.
A.
Sagara
,
H.
Tamura
,
T.
Tanaka
,
N.
Yanagi
,
J.
Miyazawa
,
T.
Goto
,
R.
Sakamoto
,
J.
Yagi
,
T.
Watanabe
,
S.
Takayama
, and
FFHR Design Group
.
Fusion Eng. Des.
89
,
2114
(
2014
).
27.
H. E.
Mynick
,
T. K.
Chu
, and
A. H.
Boozer
,
Phys. Rev. Lett.
48
,
322
(
1982
).
28.
S.
Murakami
,
A.
Wakasa
,
H.
Maaßberg
,
C. D.
Beidler
,
H.
Yamada
,
K. Y.
Watanabe
, and
LHD Experimental Group
.
Nucl. Fusion
42
,
L19
(
2002
).
29.
S.
Murakami
,
H.
Yamada
,
A.
Wakasa
,
H.
Inagaki
,
K.
Tanaka
,
K.
Narihara
,
S.
Kubo
,
T.
Shimozuma
,
H.
Funaba
,
J.
Miyazawa
et al,
Fusion Sci. Technol.
51
,
112
(
2007
).
30.
H.
Yamaguchi
,
S.
Satake
,
M.
Nakata
,
A.
Shimizu
,
Y.
Suzuki
, and
W7-X Team
,
Nucl. Fusion
61
,
106004
(
2021
).
31.
T.
Goto
,
J.
Miyazawa
,
N.
Yanagi
,
H.
Tamura
,
T.
Tanaka
,
R.
Sakamoto
,
C.
Suzuki
,
R.
Seki
,
S.
Satake
,
M.
Nunami
et al,
Plasma Phys. Controlled Fusion
60
,
074001
(
2018
).
32.
S.
Hahn
,
D. K.
Park
,
J.
Bascunan
, and
Y.
Iwasa
,
IEEE Trans. Appl. Supercond.
21
,
1592
(
2011
).
33.
G.
Motojima
,
S.
Masuzaki
,
M.
Tokitani
,
H.
Kasahara
,
Y.
Yoshimura
,
M.
Kobayashi
,
R.
Sakamoto
,
T.
Morisaki
,
J.
Miyazawa
,
T.
Akiyama
et al,
J. Nucl. Mater.
463
,
1080
(
2015
).
34.
M.
Kondo
and
Y.
Nakajima
,
Fusion Eng. Des.
88
,
2556
(
2013
).
35.
J. P. S.
Loureiro
,
H.
Fernandes
,
F. L.
Tabarés
,
G.
Mazzitelli
,
C.
Silva
,
R.
Gomes
,
E.
Alves
,
R.
Mateus
,
T.
Pereira
,
H.
Figueiredo
, and
H.
Alves
,
Nucl. Mater. Energy
12
,
709
(
2017
).
36.
A.
Cremona
,
E.
Vassallo
,
E.
Alves
,
F.
Causa
,
S.
De luliis
,
R.
Dondè
,
G.
Giacomi
,
G.
Gervasini
,
G.
Granucci
,
M.
Iafrati
et al,
Nucl. Mater. Energy
17
,
253
(
2018
).
37.
W.
Ou
,
R. S.
Al
,
J. W. M.
Vernimmen
,
S.
Brons
,
P.
Rindt
, and
T. W.
Morgan
,
Nucl. Fusion
60
,
026008
(
2020
).
38.
M.
Shimada
and
K.
Tobita
,
Plasma Fusion Res.
15
,
1401011
(
2020
).
39.
K.
Tamura
,
H.
Suzuki
,
J.
Miyazawa
,
S.
Masuzaki
, and
H.
Toyoda
,
Fusion Eng. Des.
170
,
112532
(
2021
).
40.
S.
Morita
,
C. F.
Dong
,
M.
Kobayashi
,
M.
Goto
,
X. L.
Huang
,
I.
Murakami
,
T.
Oishi
,
E. H.
Wang
,
N.
Ashikawa
,
K.
Fujii
et al,
Nucl. Fusion
53
,
093017
(
2013
).
41.
H.
Attaya
,
Fusion Technol.
19
,
1331
(
1991
).
42.
A.
Sagara
,
T.
Goto
,
J.
Miyazawa
,
N.
Yanagi
,
T.
Tanaka
,
H.
Tamura
,
R.
Sakamoto
,
M.
Tanaka
,
K.
Tsumori
,
O.
Mitarai
et al,
Fusion Eng. Des.
87
,
594
(
2012
).
43.
C.
Day
and
T.
Giegerich
,
Fusion Eng. Des.
88
,
616
(
2013
).
44.
S.
Ishiyama
,
A.
Sagara
,
H.
Chikaraishi
, and
N.
Yanagi
,
Plasma Fusion Res.
17
,
1405103
(
2022
).
45.
I.
Cristescu
,
I. R.
Cristescu
,
L.
Dörr
,
M.
Glugla
, and
D.
Murdoch
,
Fusion Sci. Technol.
52
,
667
(
2007
).