Physical Design of MW-class Steady-state Spherical Tokamak, QUEST

Similar documents
Steady-state operation scenario and the first experimental result on QUEST

Noninductive Formation of Spherical Tokamak at 7 Times the Plasma Cutoff Density by Electron Bernstein Wave Heating and Current Drive on LATE

Recent results on non-inductive startup of highly overdense ST plasma by electron Bernstein wave on LATE

Electron Bernstein Wave (EBW) Physics In NSTX and PEGASUS

Recent Development of LHD Experiment. O.Motojima for the LHD team National Institute for Fusion Science

Advanced Tokamak Research in JT-60U and JT-60SA

Overview of Spherical Tokamak Research in Japan

Study of High-energy Ion Tail Formation with Second Harmonic ICRF Heating and NBI in LHD

Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks

Introduction to Nuclear Fusion. Prof. Dr. Yong-Su Na

Plasma and Fusion Research: Regular Articles Volume 10, (2015)

Production of Over-dense Plasmas by Launching. 2.45GHz Electron Cyclotron Waves in a Helical Device

Toward the Realization of Fusion Energy

Der Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk

Microwave Spherical Torus Experiment and Prospect for Compact Fusion Reactor

Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations

ICRH Experiments on the Spherical Tokamak Globus-M

Innovative fabrication method of superconducting magnets using high T c superconductors with joints

Physics of fusion power. Lecture 14: Anomalous transport / ITER

D- 3 He HA tokamak device for experiments and power generations

Heating and Current Drive by Electron Cyclotron Waves in JT-60U

Evolution of Bootstrap-Sustained Discharge in JT-60U

Recent results from lower hybrid current drive experiments on Alcator C-Mod

DT Fusion Ignition of LHD-Type Helical Reactor by Joule Heating Associated with Magnetic Axis Shift )

Reduced-Size LHD-Type Fusion Reactor with D-Shaped Magnetic Surface )

Design window analysis of LHD-type Heliotron DEMO reactors

RELATIVISTIC EFFECTS IN ELECTRON CYCLOTRON RESONANCE HEATING AND CURRENT DRIVE

Comparison of Ion Internal Transport Barrier Formation between Hydrogen and Helium Dominated Plasmas )

D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706

Studies of Lower Hybrid Range of Frequencies Actuators in the ARC Device

Conceptual design of an energy recovering divertor

Progress of Confinement Physics Study in Compact Helical System

Extension of High-Beta Plasma Operation to Low Collisional Regime

Design concept of near term DEMO reactor with high temperature blanket

2D full wave analysis of wave structure by TASK/WF2

PHYSICS OF CFETR. Baonian Wan for CFETR physics group Institute of Plasma Physcis, Chinese Academy of Sciences, Hefei, China.

Overview of Pilot Plant Studies

Modelling of JT-60U Detached Divertor Plasma using SONIC code

Non-Solenoidal Plasma Startup in

Study of Current drive efficiency and its correlation with photon temperature in the HT-7 tokomak

Computational Study of Non-Inductive Current Buildup in Compact DEMO Plant with Slim Center Solenoid

Observations of Counter-Current Toroidal Rotation in Alcator C-Mod LHCD Plasmas

Physics and Engineering Studies of the Advanced Divertor for a Fusion Reactor

Advancing Local Helicity Injection for Non-Solenoidal Tokamak Startup

HIGH PERFORMANCE EXPERIMENTS IN JT-60U REVERSED SHEAR DISCHARGES

On Electron-Cyclotron Waves in Relativistic Non-Thermal Tokamak Plasmas

Dynamics of ion internal transport barrier in LHD heliotron and JT-60U tokamak plasmas

Driving Mechanism of SOL Plasma Flow and Effects on the Divertor Performance in JT-60U

and expectations for the future

Highlights from (3D) Modeling of Tokamak Disruptions

Influence of ECR Heating on NBI-driven Alfvén Eigenmodes in the TJ-II Stellarator

Scaling of divertor heat flux profile widths in DIII-D

MHD. Jeff Freidberg MIT

RECENT PROGRESSES ON HIGH PERFORMANCE STEADY-STATE PLASMAS IN THE SUPERCONDUCTING TOKAMAK TRIAM-1M

ICRF Mode Conversion Flow Drive on Alcator C-Mod and Projections to Other Tokamaks

Formation of An Advanced Tokamak Plasma without the Use of Ohmic Heating Solenoid in JT-60U

1. Motivation power exhaust in JT-60SA tokamak. 2. Tool COREDIV code. 3. Operational scenarios of JT-60SA. 4. Results. 5.

Possibilities for Long Pulse Ignited Tokamak Experiments Using Resistive Magnets

Additional Heating Experiments of FRC Plasma

Full-wave Simulations of Lower Hybrid Wave Propagation in the EAST Tokamak

Lower Hybrid Current Drive Experiments on Alcator C-Mod: Comparison with Theory and Simulation

Evolution of Bootstrap-Sustained Discharge in JT-60U

Magnetic Confinement Fusion and Tokamaks Chijin Xiao Department of Physics and Engineering Physics University of Saskatchewan

0 Magnetically Confined Plasma

Ion Heating Experiments Using Perpendicular Neutral Beam Injection in the Large Helical Device

Heating and current drive: Radio Frequency

D.J. Schlossberg, D.J. Battaglia, M.W. Bongard, R.J. Fonck, A.J. Redd. University of Wisconsin - Madison 1500 Engineering Drive Madison, WI 53706

STEADY-STATE EXHAUST OF HELIUM ASH IN THE W-SHAPED DIVERTOR OF JT-60U

Plan of Off-axis Neutral Beam Injector in KSTAR

ª 10 KeV. In 2XIIB and the tandem mirrors built to date, in which the plug radius R p. ª r Li

Energetic Ion Confinement and Lost Ion Distribution in Heliotrons

Fluctuation Suppression during the ECH Induced Potential Formation in the Tandem Mirror GAMMA 10

Magnetic Field Configuration Dependence of Plasma Production and Parallel Transport in a Linear Plasma Device NUMBER )

HT-7U* Superconducting Tokamak: Physics design, engineering progress and. schedule

High Beta Discharges with Hydrogen Storage Electrode Biasing in the Tohoku University Heliac

Ripple Loss of Alpha Particles in a Low-Aspect-Ratio Tokamak Reactor

Configuration Optimization of a Planar-Axis Stellarator with a Reduced Shafranov Shift )

Integrated Heat Transport Simulation of High Ion Temperature Plasma of LHD

STELLARATOR REACTOR OPTIMIZATION AND ASSESSMENT

Neutral beam plasma heating

Current Drive Experiments in the HIT-II Spherical Tokamak

Chapter IX: Nuclear fusion

Spherical Torus Fusion Contributions and Game-Changing Issues

The Path to Fusion Energy creating a star on earth. S. Prager Princeton Plasma Physics Laboratory

Fast ion physics in the C-2U advanced, beam-driven FRC

Atomic physics in fusion development

Progressing Performance Tokamak Core Physics. Marco Wischmeier Max-Planck-Institut für Plasmaphysik Garching marco.wischmeier at ipp.mpg.

1 EX/P7-12. Transient and Intermittent Magnetic Reconnection in TS-3 / UTST Merging Startup Experiments

Fusion Development Facility (FDF) Mission and Concept

Heating and Confinement Study of Globus-M Low Aspect Ratio Plasma

(Inductive tokamak plasma initial start-up)

1 EX/P6-5 Analysis of Pedestal Characteristics in JT-60U H-mode Plasmas Based on Monte-Carlo Neutral Transport Simulation

Toroidal confinement devices

(a) (b) (c) (d) (e) (f) r (minor radius) time. time. Soft X-ray. T_e contours (ECE) r (minor radius) time time

Non-inductive plasma startup and current profile modification in Pegasus spherical torus discharges

Plasma Wall Interactions in Tokamak

ASSESSMENT AND MODELING OF INDUCTIVE AND NON-INDUCTIVE SCENARIOS FOR ITER

DEMO Concept Development and Assessment of Relevant Technologies. Physics and Engineering Studies of the Advanced Divertor for a Fusion Reactor

GA A23713 RECENT ECCD EXPERIMENTAL STUDIES ON DIII D

Fast ion generation with novel three-ion radiofrequency heating scenarios:

Transcription:

Physical Design of MW-class Steady-state Spherical Tokamak, QUEST FT/P3-25 K.Hanada 1), K.N.Sato 1), H.Zushi 1), K.Nakamura 1), M.Sakamoto 1), H.Idei 1), M.Hasegawa 1), Y.Takase 2), O.Mitarai 3), T.Maekawa 4), Y.Kishimoto 4), M.Ishiguro 5), T.Yoshinaga 6), H.Igami 6), S.Kawasaki 1), H.Nakashima 1), A.Higashijima 1), Y.Higashizono 1), A.Ando 7), N.Asakura 8), A.Ejiri 2), Y.Hirooka 6), A.Ishida 9), A.Komori 6), M.Matsukawa 8), O.Motojima 6), Y.Ogawa 2), N.Ohno 10), Y.Ono 2), M.Peng 11), S.Sudo 6), H.Yamada 6), N.Yoshida 1), Z.Yoshida 2). 1) Research Institute for Applied Mechanics, Kyushu University 2) University of Tokyo 3) Tokai University 4) Kyoto University 5) IGSES, Kyushu University 6) National Institute for Fusion Science 7) Tohoku University 8) Japan Atomic Energy Agency 9) Niigata University 10) Nagoya University 11) Fusion Energy Division, ORNL E-mail:hanada@triam.kyushu-u.ac.jp QUEST (R=0.68m, a=0.4m) focuses on the steady state operation of the spherical tokamak (ST) by controlled PWI and electron Bernstain wave (EBW) current drive (CD). The QUEST project will be developed along two phases, phase I: steady state operation with plasma current, I p =20-30kA on open divertor configuration and phase II: steady state operation with I p =100kA and β of 10% in short pulse on closed divertor configuration. Feasibility of the missions on QUEST was investigated and the suitable machine size of QUEST was decided based on the physical view of plasma parameters. Electron Bernstein wave (EBW) current drive are planned to establish the maintenance of plasma current in steady state. Mode conversion efficiency to EBW was calculated and the conversion of 95% will be expected. A new type antenna for QUEST has been fabricated to excite EBW effectively. The situation of heat and particle handling is challenging, and W and high temperature wall is adopted. The start-up scenario of plasma current was investigated based on the driven current by energetic electron and the most favorable magnetic configuration for start-up is proposed. 1. Introduction It is important to conduct the academic basics research supporting high beta and steady state operation approaches. The QUEST project focuses on the steady state operation of the spherical tokamak (ST) which has the capability to attain high β rather than conventional tokamaks. A final target of the project is the steady state operation of ST with relatively high β (10%) under controlled plasma wall interaction (PWI). The program will be developed in increment step such as, I. low β steady state operation in limiter configuration, II. low β steady state operation in divertor configuration, III relatively high β steady state operation in closed divertor configuration. The specific purposes in phase I and II are: (1) To examine the steady state current drive and plasma generation by EBW. (2) To comprehensively establish recycling control based on wall temperature control, advanced wall control under high plasma performance. (3) To improve diverter concepts and to establish the way of controlling particles and heat loads during long duration operation.

FT/P3-25 Fig. 1 Left: QUEST (R=0.68m, a=0.4m) and Right: first plasma. fiscal year 05 items H17 06 07 08 09 10 11 12 13 14 H18 H19 H20 H21 H22 H23 H24 H25 H26 further construction High ξ Plasma start-up >10% (1sec) RF+OH 20% RF+OH+NBI Current drive RF RF(8.2GHz)+NBI 16GHz PWI W W, high Temp. wall Divertor fueling open Gas puff Control of Recycling closed advanced CT injection pellet Table 1 Time table of the important items of physics issues of QUEST. Vertical blue line show the end of construction phase and two red lines show the end of phase I and II. Present status of QUEST is shown in Fig. 1 and the time table of the important items is shown in table. 1. One of the most important items is the investigation of PWI. QUEST will be used no carbon in the vacuum vessel and only metal wall, that is mainly W, will be available. NBI of 2MW is planned, however it can be operated in short pulse (<1s). NBI can used for the plasma start up and high β experiment in the short pulse. In steady state operation of QUEST. RF of 1MW will work mainly to maintain and heat the plasma. Especially plasma

start-up is one of the key issues and only RF and RF+ohmic heating (OH) will be tested. The available magnetic flux for OH is less than 0.3 Vs, which is not enough to obtain the plasma current up to 0.3MA. Therefore the help of RF is quite significant in QUEST. Open divertor is installed in phase I to investigate heat and particle flux to the divertor plate. Based on the experimental results in Phase I, closed divertor will be designed and tested in Phase II. Efficient fueling is important in steady state operation. Advanced fueling is planned on QUEST, which is compact toroid (CT) injection. The fueling using CT injection is challenging, however it is suitable way in low magnetic field device such as ST. 2. Purpose of the QUEST project Pinj τ dis /R [MW s/m] 10 5 10 4 10 3 10 2 TRIAM-1M QUEST LHD JT-60U Power Plant ITER ne Vpτ D /( τ p R) [1/m] 10 27 10 26 10 25 10 24 10 23 TRIAM-1M LHD QUEST JT-60U Power Plant ITER 10 1 0.01 0.1 1 10 P inj /R [MW/m] 10 22 10 20 2 4 6 10 21 2 4 6 10 22 2 4 6 10 23 n e V p /(τ p R) [1/m s] Fig. 1 Operation region of QUEST on heat (left) and particle (right) handling. Left: The vertical axis is injected power (MW) x discharge duration (s) / major radius (R), which is approximately proportional to heat load to divertor and the horizontal axis is injected power (MW) / major radius (R), which is approximately proportional to heat flux to the divertor. Right: The vertical axis is density (m -3 ) x plasma volume (m 3 ) x discharge duration (s) / major radius (R) /particle confinement time (s), which is approximately proportional to fluence to the divertor and the horizontal axis is density (m -3 ) x plasma volume (m 3 ) / major radius (R) /particle confinement time (s), which is approximately proportional to particle flux, where we assume particle confinement time equals energy confinement time. One of the main purposes of QUEST is an achievement of long duration discharge with MW-class injected power. Operation region of QUEST in view of heat handling is shown in Fig. 2. The operation of QUEST in heat handling is so challenging and PWI control on QUEST plays an essential role in realization of a power plant. Moreover, particle handling is crucial to obtain the steady state operation of magnetic confinement devices. As shown in Fig. 3, QUEST is at a dispensable position towards DEMO in operation region on particle handling. Especially, wall phenomena on PWI, which are co-deposition, T retention, dust and so on, strongly depends on fluence rather than particle flux. In this meaning, QUEST has a competitive edge over the other devices. QUEST will equip full-scale closed divertor to control particle flux around divertor region. First wall of QUEST will be completely surrounded by W and keeps high temperature (>600K) during discharge at least to obtain recycling ratio, R=1. A steady state operation under R=1 marks a milestone for realization of steady state magnetic fusion power plant.

3. Feasibility of the missions on QUEST QUEST has two different types of missions, which are high β Τ (>10%) in short pulse (< 1s) and low β T in long pulse (>1h). A plasma size should be decided by the physical requirement from the missions. The value of toroidal field of QUEST, which is 0.25 T at plasma center was decided in view of the position of fundamental electron cyclotron resonance (ECR) and 2 nd one for the frequency of the present RF power sources (f=8.2ghz). Figure 4 shows a feasibility window on minor and major radii considering about the abovedescribed missions. The solid and dotted lines show the line for q 95 =4 and 3, respectively. Generally speaking, it is difficult to operate on q 95 <3 in tokamaks and the solid line provide a low q limit. The solid and dotted line shows the heat load limitation (10MW/m 2 ) to a divertor plate under the condition of the fraction of radiation, P rad =40% and 50% in single null configuration. In steady state, the heat load to any plasma facing components should be kept down on less than 10MW/m 2. This is also the requirement for ITER. When the major radius is less than the solid line, the heat flux to divertor plate goes beyond the boundary of the ITER requirement. The solid and dotted lines show the aspect ratio, A=1.4 and 1.3, respectively. The reduction of aspect ratio leads to the difficulty of the construction of the machine. The solid and dotted lines show the constant β of 10% and 5% based on the ITER 89P L-mode scaling under the condition of I P =300kA, B T =0.25T, P inj =3MW, respectively. The selected plasma size β=5% κ=2.55 δ=0.68 f rad =40% f rad =50% β=10% Fig. 4 Plasma size determination. (a point on the figure) has the candidate to attain the mission. 4. Preparation for EBWCD in steady state operation In Phase I, the most important issue of the project is to achieve steady state operation. Non-inductive current drive is indispensable for steady state operation. The different way to drive the plasma current from that in conventional tokamaks should be developed for ST.

Electron Bernstein wave current drive (EBWCD) and higher harmonics fast wave current drive (HHFWCD) are candidates to establish the maintenance of plasma current in steady state by use of RF. As we do not have any heating source for HHFWCD in the view of frequency at present, EBWCD is suitable way to apply to the machine at first. The current drive efficiency of 0.5 x 10 18 A/W/m 2 at considerable low electron temperature is an overarching goal for current drive in this phase. At first, two excitation scenario of EBW (O-X-B and X-B) should be discussed, because injection way of RF depends on the scenario of excitation. The O-X-B scenario is utilized a mode conversion from O-mode to X-mode at the cut-off region for O mode. The conversion rate from O-mode to X-mode is well-understood and depends on the injection angle to the cut-off layer. When we are willing to adopt this scenario, the adjustment of the injection angle is crucial. The converted X-mode will transfer to the upper hybrid resonance (UHR) from high field side and it will convert to EBW. EBW propagates inwards again and is absorbed bulk electrons effectively. By the effect of the magnetic shear, EBW has the single directed momentum in specified situation and it delivered to plasma. Accordingly plasma is provided single directed momentum from EBW and plasma current can be driven. The X-B scenario is the same in delivery and receipt of the momentum, however the mode conversion process is different. The RF will injected perpendicular to the magnetic field and 3 wave coupling play an essential role in mode conversion process and the mode conversion efficiency significantly depends on the scale length of electron density at UHR. To convert to EBW, we should control the density gradient at UHR. Figure 5 shows the conversion map of EBW [1]. More than 95% mode conversion rate to EBW can be obtained in the three types of parameter regions, the Perp-X, the g1 and the O-X-B, where O-X-B means O-mode wave converts to X-mode at plasma cutoff and then the converts to EBW and Perp-X means X-mode penetrates right hand cutoff and converts to EBW at upper hybrid resonance, and g1 means the optimal mixture of the O and the X-modes. Increasing the peak density, the optimal injection mode for EBW varies from the O-mode to the g1 mode. Therefore the injected RF should be modified during a shot. To realize the scenario, we designed and fabricated a new phased-array antenna [2] and it will be installed on QUEST. Experimental observations of EBWCD are obtained in various devices. Driven current of 100 ka at 60GHz 600kW was achieved in COMPASS-D on the O-X-B scenario. Full non-inductive plasma current up to 20 ka in a sequence of plasma start-up on CS-less configuration was achieved on LATE [3]. The X-B scenario was executed on TST-2 by appropriating RF power source of 8.2 GHz, 200kW [4]. Plasma heating could be observed in the case of application to ohmic heated plasma. Full non-inductive plasma current of 4kA for 0.3 sec could be sustained by only-rf on TST-2. A simulation was done for high β discharge on NSTX and accordingly non-inductive plasma current of 30 ka will be obtained by 1MW injection. Fokker-Planck simulations were executed and some results were presented by B.Jones and C.B.Forst for NSTX and Pegasus, respectively. According to their paper, expected current drive efficiencies were η CD =2.7X10 18 for NSTX and 1.5x10 18 for Pegasus. These values are not far from the target value of QUEST. Fokker-Planck simulation was carried out by TASK/FP [5].

Fig.5 The conversion map of EBW on a slab model. Ln/λ is the ratio of the density scale-length at the plasma cutoff to the vacuum wavelength and ω/ω is the ratio of the wave frequency to the electron cyclotron frequency. The numbers in this figure correspond to the peak values of the given density profiles of QUEST. The red shaded region indicates the area where the mode-conversion rate is more than 95%, while the green meshed region indicates the area where the mode-conversion rate is less than 10%. 5. Start-up scenario by use of RF ST has small aspect ratio (A = R 0 /a, where R 0 is the major radii and a is the minor radii of tokamak plasma) and its central space is limited [1]. Therefore, inductive electric field in ST is weaker than that in conventional tokamak, and it is required that plasma current start-up is conducted by non inductive method such as radio frequency (RF) wave injection or ohmic method supported by RF injection (RF + OH) in very low inductive electric field. In the RF injection method, it is pointed that energetic electrons, which is tail electrons, confined in open field contribute to the plasma current initiation [6]. Thus, plasma current initiation with RF injection has quite different physics from the normal OH initiation with bulk electrons, which is conducted in conventional tokamak. This peculiar condition applies to superconductive tokamak because its inductive electric field is also weak and it needs to use RF injection to initiate plasma. Therefore we can use the knowledge of ST plasma ramp-up to conduct current initiation in operation of superconductive tokamak, like ITER. In terms of it, the study of ST current ramp-up is important to realize future fusion reactor. Figure 7 shows that the orbits of electrons in the opened field. Electron orbits are classified three patterns. They are named lost, banana, and co-moving respectively. The electrons which are not confined in the opened field are classified into lost. The electrons classified into banana are trapped in opened field by magnetic mirror because of toloidal field dependence on major radii R, B=B 0 R 0 /R and poroidal field which transports the electrons to higher magnetic field strength side. The pattern named co-moving is for the electrons confined in the open field without magnetic mirror effect. Electrons in tokamak configuration have two poroidal velocities. One is the poroidal component of velocity along

the magnetic field line and the other is drift velocity generated by toloidal field curvature and gradient. When these two velocities have opposite direction each other and almost same value, electrons hardly move to poroidal direction and they are confined in small area in opened field. Which orbits a electron traces depends on its initial velocity, pitch angle, poroidal field shape and strength. Thus poroidal field configuration is important to confine electrons in open field. Especially, N-index on mid plane should be positive. Expected plasma current driven by energetic electrons, I es, is calculated in the following way. At first we assume that the velocity distribution function of electron is Maxwell distribution and loss electrons do not carry any plasma current. The start point of orbit calculation sets on the ECR layer on the midplane. Density and temperature of energetic electrons are 1x10 17 m -3 and 1keV, respectively. Figure 8 shows I es as the function of poloidal magnetic field, Bp in the magnetic configuration on QUEST. From this result, we expect the most suitable magnetic configuration for start-up of plasma current. Fig.7 Orbits of electrons in open field. Solid lines show orbits projection and wall of vacuum chamber, and dashed line show magnetic line projection on poroidal plane. The magnetic field strength in left side in the figure is stronger than that in right side. Ies[kA] 1.5 1.0 0.5 2 4 6 Bp[T] 8 10 x10-3 Fig. 8 Expected plasma current driven by energetic electron in open magnetic configuration as the function of poloidal magnetic field, B p.

Summary One of the main purposes of QUEST is an achievement of long duration discharge with MW-class injected power in ST. The QUEST project will be developed along two phases, phase I:steady state operation of 20-30kA on open divertor configuration and phase II: steady state operation of 100kA and β of 10% in short pulse on closed divertor configuration. Time table for important physical items of QUEST project is proposed. Feasibility of the missions on QUEST was investigated and the suitable machine size of QUEST was decided based on the physical view of plasma parameters. Non-inductive current drive is indispensable for steady state operation. Electron Bernstein wave current drive (EBWCD) are promising way to establish the maintenance of plasma current in steady state. Mode conversion efficiency to EBW was calculated and 95% conversion will be expected in narrow window. New antenna system was designed to sophisticate the injected wave. The start-up scenario of plasma current was investigated based on the driven current by energetic electron and the most favorable magnetic configuration is proposed. Acknowledgements This work is performed with the support and under the auspices of the NIFS Collaboration Research Program (NIFS05KUTR014, NIFS08KUTR022). This work is partially supported by a Grant-in-aid for Scientific Research from Ministry of Education, Science and Culture of Japan. References [1] H.Igami et al. PPCF 48 (2006) 573-598. [2] H. Idei et al. Proc. 32nd International Conf. on Infrared and Millimetre Waves, 2007. [3] T.Yoshinaga et al. PRL 96, 125005 (2006). [4] S.Shiraiwa et al. PRL 96, 185003 (2006) [5] H.Idei et al. ICCP2008 BEH P2-157. [6] A.Ejiri and Y.Takase Nuclear Fusion, Vol. 47 (2007) 403-416.