Assessment of Global Stellarator Confinement: Status of the International Stellarator Confinement Data Base

Similar documents
Confinement Study of Net-Current Free Toroidal Plasmas Based on Extended International Stellarator Database

Effect of rotational transform and magnetic shear on confinement of stellarators

Confinement Study of Net-Current Free Toroidal Plasmas Based on Extended International Stellarator Database

1 EX/P5-9 International Stellarator/Heliotron Database Activities on High-Beta Confinement and Operational Boundaries

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

1) H-mode in Helical Devices. 2) Construction status and scientific objectives of the Wendelstein 7-X stellarator

Exploration of Configurational Space for Quasi-isodynamic Stellarators with Poloidally Closed Contours of the Magnetic Field Strength

Impurity issues frequently mentioned

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

A THEORETICAL AND EXPERIMENTAL INVESTIGATION INTO ENERGY TRANSPORT IN HIGH TEMPERATURE TOKAMAK PLASMAS

Significance of MHD Effects in Stellarator Confinement

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

BRIEF COMMUNICATION. Near-magnetic-axis Geometry of a Closely Quasi-Isodynamic Stellarator. Greifswald, Wendelsteinstr. 1, Greifswald, Germany

Der Stellarator Ein alternatives Einschlusskonzept für ein Fusionskraftwerk

GA A23736 EFFECTS OF CROSS-SECTION SHAPE ON L MODE AND H MODE ENERGY TRANSPORT

Shear Flow Generation in Stellarators - Configurational Variations

H-mode in Helical Devices

Active MHD Control Needs in Helical Configurations

Confinement and edge studies towards low ρ* and ν* at JET

Multifarious Physics Analyses of the Core Plasma Properties in a Helical DEMO Reactor FFHR-d1

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

On the physics of shear flows in 3D geometry

Integrated Bayesian Experimental Design

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

TH/P8-4 Second Ballooning Stability Effect on H-mode Pedestal Scalings

Stellarators. Dr Ben Dudson. 6 th February Department of Physics, University of York Heslington, York YO10 5DD, UK

- Effect of Stochastic Field and Resonant Magnetic Perturbation on Global MHD Fluctuation -

Ion orbits and ion confinement studies on ECRH plasmas in TJ-II stellarator

Impact of Toroidal Flow on ITB H-Mode Plasma Performance in Fusion Tokamak

Physics Considerations in the Design of NCSX *

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

Electron Transport and Improved Confinement on Tore Supra

Finite-Orbit-Width Effect and the Radial Electric Field in Neoclassical Transport Phenomena

Spatial, temporal and spectral structure of the turbulence-flow interaction at the L-H transition

Configuration Control Experiments in Heliotron J

Numerical calculation of the Hamada basis vectors for three-dimensional toroidal magnetic configurations

UNDERSTANDING TRANSPORT THROUGH DIMENSIONLESS PARAMETER SCALING EXPERIMENTS

Characterization of the Perpendicular Rotation Velocity of the Turbulence by Reflectometry in the Stellarator TJ-II

Scaling of divertor heat flux profile widths in DIII-D

Long-Range Correlations and Edge Transport Bifurcation in Fusion Plasmas

Integrated Heat Transport Simulation of High Ion Temperature Plasma of LHD

Quasi-Symmetric Stellarators as a Strategic Element in the US Fusion Energy Research Plan

ISOTOPE EFFECTS ON CONFINEMENT AND TURBULENCE IN ECRH PLASMA OF LHD

Energetic Ion Confinement and Lost Ion Distribution in Heliotrons

Confinement in Wendelstein 7-X Limiter Plasmas

Progress of Confinement Physics Study in Compact Helical System

STELLARATOR REACTOR OPTIMIZATION AND ASSESSMENT

Particle transport results from collisionality scans and perturbative experiments on DIII-D

19 th IAEA- Fusion Energy Conference FT/1-6

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

Role of Magnetic Configuration and Heating Power in ITB Formation in JET.

Presented on the 15th International Stellarator Workshop, Madrid, 3-7 Oct Submitted to Fusion Science and Technology Revised 20 March 2006

Reduction of Neoclassical Transport and Observation of a Fast Electron Driven Instability with Quasisymmetry in HSX

Triggering Mechanisms for Transport Barriers

Innovative Concepts Workshop Austin, Texas February 13-15, 2006

Particle Transport Measurements in the LHD Stochastic Magnetic Boundary Plasma using Mach Probes and Ion Sensitive Probe

Columbia Non-neutral Torus completes construction. In this issue... and starts operation

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

Extension of High-Beta Plasma Operation to Low Collisional Regime

Initial Experimental Program Plan for HSX

Direct drive by cyclotron heating can explain spontaneous rotation in tokamaks

Integrated Transport Modeling of High-Field Tokamak Burning Plasma Devices

Design window analysis of LHD-type Heliotron DEMO reactors

Dimensionless Identity Experiments in JT-60U and JET

Flow, current, & electric field in omnigenous stellarators

Overview of Tokamak Rotation and Momentum Transport Phenomenology and Motivations

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

Self-consistent modeling of ITER with BALDUR integrated predictive modeling code

MHD Simulation of High Wavenumber Ballooning-like Modes in LHD

Pellet injection experiments in LHD

Development of TASK3D-a, and its linkage to Stellarator-Heliotron Database (ISH-DB)

Advances in stellarator gyrokinetics

The effect of the radial electric field on neoclassical flows in a tokamak pedestal

Electrode Biasing Experiment in the Large Helical Device

and expectations for the future

The Initial Program of Wendelstein7-X on the Way to a HELIAS Fusion Power Plant

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

Study of neoclassical transport and bootstrap current for W7-X in the 1/ν regime, using results from the PIES code

Sizing Up Plasmas Using Dimensionless Parameters

Neoclassical transport

First Observation of ELM Suppression by Magnetic Perturbations in ASDEX Upgrade and Comparison to DIII-D Matched-Shape Plasmas

The EPED Pedestal Model: Extensions, Application to ELM-Suppressed Regimes, and ITER Predictions

Power Deposition Measurements in Deuterium and Helium Discharges in JET MKIIGB Divertor by IR-Thermography

Erosion and Confinement of Tungsten in ASDEX Upgrade

Results From the International Collaboration on Neoclassical Transport in Stellarators (ICNTS)

Correlation Between Core and Pedestal Temperatures in JT-60U: Experiment and Modeling

Progress in characterization of the H-mode pedestal

2. STELLARATOR PHYSICS. James F. Lyon James A. Rome John L. Johnson

Inter-linkage of transports and its bridging mechanism

Comparison of ITER Performance Predicted by Semi-Empirical and Theory-Based Transport Models

3D analysis of impurity transport and radiation for ITER limiter start-up configurations

The performance of improved H-modes at ASDEX Upgrade and projection to ITER

Observation of Neo-Classical Ion Pinch in the Electric Tokamak*

MHD instability driven by supra-thermal electrons in TJ-II stellarator

TURBULENT TRANSPORT THEORY

C-Mod Core Transport Program. Presented by Martin Greenwald C-Mod PAC Feb. 6-8, 2008 MIT Plasma Science & Fusion Center

Introduction to Integrated Data Analysis

GA A23114 DEPENDENCE OF HEAT AND PARTICLE TRANSPORT ON THE RATIO OF THE ION AND ELECTRON TEMPERATURES

MHD. Jeff Freidberg MIT

ELM Suppression in DIII-D Hybrid Plasmas Using n=3 Resonant Magnetic Perturbations

Transcription:

Assessment of Global Stellarator Confinement: Status of the International Stellarator Confinement Data Base A.Dinklage (1,*), E.Ascasibar (2), C.D.Beidler (1), R. Brakel (1), J. Geiger (1), J.H. Harris (3), A. Kus (1), S. Murakami (4), S. Okamura (5), R. Preuss (1), F. Sano (4), U. Stroth (6), Y. Suzuki (5), J. Talmadge (7), V.Tribaldos (2), K.Y.Watanabe (5), A. Weller (1), H.Yamada (5), M.Yokoyama (5) (1) Max-Planck-Institut für Plasmaphysik, EURATOM-Association, Greifswald, Germany (2) Laboratorio Nacional de Fusión, EURATOM-CIEMAT, 28040 Madrid, Spain (3) Australian National University, Canberra, Australia (4) Kyoto University, Kyoto, Japan (5) National Institute for Fusion Science,Toki,Gifu 509-5292, Japan (6) Institut für Plasmaforschung, Universität Stuttgart Germany (7) University of Wisconsin, Madison, USA Corresponding Author: A. Dinklage, Max-Planck-Institut für Plasmaphysik, EURATOM Association, Wendelsteinstr.1, 17491 Greifswald, Germany Tel.: +49 38 34 88 23 28, Fax: +49 38 34 88 20 09, e-mail: dinklage@ipp.mpg.de 28 pages, 6 figures, 0 tables

Assessment of Global Stellarator Confinement: Status of the International Stellarator Confinement Data Base (a) A.Dinklage (1,*), E.Ascasibar (2), C.D.Beidler (1), R. Brakel (1), J. Geiger (1), J.H. Harris (3), A. Kus (1), S. Murakami (4), S. Okamura (5), R. Preuss (1), F. Sano (4), U. Stroth (6), Y. Suzuki (5), J. Talmadge (7), V.Tribaldos (2), K.Y.Watanabe (5), A. Weller (1), H.Yamada (5), M.Yokoyama (5) (1) Max-Planck-Institut für Plasmaphysik, EURATOM-Association, Greifswald, Germany (2) Laboratorio Nacional de Fusión, EURATOM-CIEMAT, 28040 Madrid, Spain (3) Australian National University, Canberra, Australia (4) Kyoto University, Kyoto, Japan (5) National Institute for Fusion Science,Toki,Gifu 509-5292, Japan (6) Institut für Plasmaforschung, Universität Stuttgart Germany (7) University of Wisconsin, Madison, USA (*) dinklage@ipp.mpg.de Abstract Further exploration of confinement trends in the evolving International Stellarator Confinement Database is reported. The impact of configurations on confinement is confirmed and the performance close to operational limits discussed. Model comparison techniques allow for tests of physical models against data. Correlations of configuration-describing parameters against an empirical confinement enhancement factor are investigated.

1. Introduction Because of their three-dimensional nature, stellarator/heliotron devices cover a large magnetic configuration space. Since the ultimate goal of stellarator research is to develop an alternative fusion reactor concept, the exploration of the most promising configurations requires comparative assessments of the plasma performance and how different aspects of a 3D configuration influence it. To this end, the International Stellarator Confinement Database (ISCDB) project was been re-initiated in 2004 to extend the ISS95 [1] database to roughly 3300 discharges from nine different devices. The revision of a data set restricted to comparable scenarios led to the ISS04 scaling law [2] which confirmed the general tendencies of ISS95, but also revealed clearly the necessity of incorporating configuration-descriptive parameters. In other words, an extension beyond the set of regression parameters used for ISS95/ISS04 appears to be necessary, and candidates, such as the flux-surface elongation are investigated. Since grouping of data is a keyissue for deriving ISS04, basic assumptions are revisited in this paper, e.g. the dependence on the heating scheme. Furthermore, the statistical approaches used are assessed with respect to their impact on the scaling. A crucial issue is the weighting of data groups which is discussed in terms of error-in-variable techniques. Bayesian model comparison is employed for testing scaling hypotheses based on scaling invariance principles thus allowing the assessment of applicability of theory-based scaling laws on stellarator confinement. Again, subsets of ISCDB are used in order to ensure comparability of the data under investigation.

2. The ISCDB dataset and the ISS04 stellarator scaling Plotting the data from the new ISCDB data set against the original ISS95 scaling relation (Fig. 1) shows subgroup clusters around an overall trend. Systematic offsets are visible between the various groups of data: for example, the LHD data indicate enhanced performance compared to the ISS95 prediction. This motivates to confirm specifically the gross scaling of regression parameters (minor and major radius (a (m), R(m)), mean density n (10 19 m -3 ), absorbed power P (MW), magnetic field B (T) and rotational transform ι). The tokamak data [3] in Fig. 1 are shown for comparison. This comparison is qualified because of differences in the set of regression parameters and assumptions on the rotational transform profiles [1] as follows. From the tokamak database the tok tok elongation corrected minor radius ( a κ ) and 2 ι / 3 = 1/ q2 / 3 was calculated using an 4 2 ( ) ι = a ρ where ρ is the effective radius. ι a ( ρ = 1) was 2 ι -profile as ( ρ) ι 1 ( 1 ρ ) κ 2 derived from I P a 2 = 2.5 B t R 1 ( 1 + ) a ι which relates the plasma current I P to the edge rotational transform and the toroidal magnetic field B t. The ITER prediction is consistent with ELMy-H mode data. Although a number of non-validated assumptions enter the comparison, the stellarator data join the high confinement data from tokamaks in the range of confinement times covered. Hence, there is no indication of a degraded confinement in particular from the data of the largest device, LHD. A second conclusion can also be drawn from Fig. 1: Data from different devices form clusters whose coarse tendency may differ from the general trend of ISS95. A closer inspection even shows subgroups within discharges of a single device. This observation led to the ad hoc assumption of ISS04: stellarator configurations are

different, but a common scaling in the leading regression parameters is to be investigated. Technically, ISS04 assigns to each subgroup a configuration parameter which is determined as the ratio of experimentally determined confinement times and a reference scaling (here ISS95 was chosen, f is normalized to the W7-AS (limiter) subgroup) [2]. The regression of renormalized data results in: τ ISS 04 E = 0.134a τ Bohm ρ 2.28 R 0.79 * ν 0.64 P * 0.00 0.61 β n 0.54 e 0.19 a 0 * B ι 0.84 0.41 2 / 3 (1) where a * reflects the minor radius dependencies not contained in the dimensionless scaling parameters (normalized collisionality ν*, normalized gyro-radius ρ * and plasma β). A vanishing a * dependence indicates the scaling to be dimensionally correct. The RMSE value of the regression is 0.0267 if the weighted renormalized confinement time is employed as regression data. Using the original confinement time data gives an RMSE of 0.2366. The scaling law ISS04 confirms the findings of ISS95 with respect to the scaling exponents which reads: τ ISS 95 E = 0.079a τ Bohm ρ 2.21 R ν 0.71 * 0.65 P * 0.04 0.59 β n 0.16 0.51 e a 0 * B ι 0.53 0.4 2 / 3 (2) More details, in particular on assessments of assumptions on the regression ansatz can be found in Ref. [2]. Figure 2 shows the ISS04 scaling results in conjunction with the configuration factors used for deriving the ISS04 scaling law. In addition to the

ISS04 subset of ISCDB (1721 from 3335 datasets) the remaining ISCDB data are plotted in Fig. 2 for the purpose of comparison. The configuration factor was derived as for the ISS04 subgroups. The dependence of f on the reference scaling is vanishing if the derived scaling is used iteratively. In addition, the normalization factor ansatz is independent from the choice of reference; here ISS95 [1], Lackner-Gottardi [5] and LHD [6] scaling laws were employed as initial reference. For all references the iterative determination converges within one standard deviation of the regression exponents results against τ ISS 04, E = f ( configuration) 0.137 a τ Bohm ρ 1.22 * ν * 0.11 β 0.06 a 0 * 2.40 R 0.64 P 0.57 n 0.55 e B ι 0.90 0.09 2 / 3 (3) The RMSE value of the scaling is 0.0262 (0.2361 if the original data were used as reference). This scaling law is close to ideal gyro-bohm scaling laws [7]. The dependence on the rotational transform differs substantially from ISS95, but nevertheless the ι dependence is significant regardless which reference is initially assumed. Due to technical settings and physical constraints the conditioning of a global confinement data-set from fusion devices is always restricted. This appears, e.g. in Ref. [2] as an apparent correlation of the scaling assumptions (weighting, normalization) and the scaling exponents. Such correlations were described in detail in Ref. [2] and the strategy to identify the role of correlations and data clustering is a detailed assessment of the data-sets both with respect to the physical properties (next section) as well as consistency of scaling models with the data as discussed in the over-next section. The collinearity between f and other engineering parameters is to be studied in future.

3. Assessment of ISCDB data The ISS04 scaling was derived from a subset of the International Stellarator Confinement Database. Partly, the subset for scaling was restricted to data which are expected to store most of their energy as thermal energy. Therefore, density and power scans from the HSX device were not included in scaling studies. The dependence on the heating scheme was investigated by applying the ISS04 scaling procedure on electron cyclotron heated plasma data. The resulting scaling exponents are consistent with the ISS04 scaling law except for the scaling in density and magnetic field (τ ~ n 0.63±0.02 B 0.63±0.03 instead of τ ~ n 0.54±0.01 B 0.84±0.02 ). The differences, however, may in part be attributed to a density magnetic field correlation in the data of ECRH plasmas (correlation coefficients r ECH = 0.81 between log(n) and log(b) vs. r ISS04 = 0.18 for the full data set). An example of such a subset is a density scan of the confinement time in the transition from normal confinement to the high density H-mode (HDH) in Wendelstein 7-AS. This transition is shown in the left frame of Fig. 3. This example shows that as the density rises, the energy confinement time jumps to a different confinement mode. At still higher densities, confinement then degrades with densitiy as the plasma becomes detached from the divertor targets in divertor operation [8,9]. A second example of regression parameter variations is the comparison of the ι dependence of a W7-AS subset [10] and TJ-II data [11] shown in the right plot of Fig. 3. Taking account for the configuration factor, the ι scaling of different configurations can be compared. This approach allows one to make use of the large parameter variations of an inter-machine database. Since ISS04 corresponds to the τ exp =fτ ISS04

line in the right plot of Fig. 3 the gross scaling of ι appears to be consistent from the inter-machine comparison. However, the variation of the data indicates possible uncertainties even in the gross scaling; for comparison the green line indicates what would be expected for vanishing ι dependence in ISS04. The fine structure of the ι dependence due to low order rationals is a salient feature of low-shear stellarator confinement and is extensively documented and discussed in [10,11]. Beyond scaling behavior the documentation of high confinement modes is also started to be addressed in ISCDB. Fig. 4 shows selected examples of the configuration factor under mode transitions, e.g. for HDH bifurcations and a single example for a L- to H-mode transition in W7-AS. Please note that the renormalization factors of the single L-H-mode transition are within the error margins of the renormalization factor of W7-AS limiter discharges (see Fig. 2) but show the confinement enhancement in H- modes. The configuration scan in LHD documents the enhancement of confinement due to inward shift of the magnetic configuration [12,13]. 4. Assessment of physical models A crucial issue in the assessment of scaling laws is the estimate of trustworthiness not only for reference purposes but ultimately for a predictive statement. One possibility to link the scaling relation with physics models is to derive dimensional constraints on the scaling exponents from invariance requirements of the basic equations describing plasma behavior connected with collisionality ν* or plasma β [14]. In order to compare these physics models on the basis of data from the ISCDB, Bayesian probability theory is employed. For model comparison an inclusion of uncertainties is necessary. Neglecting these uncertainties means assigning equal weight

to every entry in the database. However, the uncertainties of the data vary substantially being the result of different diagnostics measured in different machines during different experimental campaigns. The error statistics taking into consideration the uncertainty both in the regression target (τ E with error σ τ ) and the machine variables (so called errors in variables approach), can be written in terms of a probability distribution, i.e. likelihood, as [15,16]: 2 τ E α k xk p r r r const ( ) 1 = k x, τ E α, σ exp K K 2 2 2 σ + ( ) 2 2 τ α kσ k 2π σ k k k (4) where x r represents the K = 6 ISS regression variables (a,r,p,n,b,ι). α k and σ k are their respective regression exponents and errors.the denominator in the argument of the exponent marks the difference with ordinary least squares fitting. Since the present content of ISCDB covers different ranges of collisionality and plasma β, the validity of low/high collisionality and low/high models can be performed. The model probability is proportional to the posterior probability of the data, i.e. marginalizing the above likelihood over all α k and normalizing the result with respect to the available models. The sum of the model probabilities for the four models becomes unity. A methodological summary is to be found in Ref. [17,18]. A first study focuses on the low/high β and ν data from W7-AS. As can be seen in the right graph of Fig. 5, the low β data set is given by ν < 23 and β < 0.004 giving a

total of 223 shots, on the other hand the high β data set to ν* > 100 and β > 0.015 giving a total of 43 shots. As for ISS04, the rotational transform subdivides the W7-AS data and shots with ι<0.4 were taken into account only. The result of model comparison calculations is shown in the left graph of Fig. 5. The models under consideration are reviewed in Ref. [19]. It appears that the set designated as low-βdata can be described by a collisional low β Connor-Taylor model (electrostatic Fokker-Planck). It is to be noted that this subset result is not in agreement with the dimensionless formulation of ISS04. Since the model probabilities do not prefer a distinct solution, the model assessment fails for the high β data. The methods applied now allows one for a successive enhancement of significance of model assessment and is used as a guideline for future work, particularly for the β scaling of confinement in stellarator/heliotron devices. 5. Configuration factor assessment ISS04 is derived with renormalized data and it confirms the previous scaling law results. The variety of configuration factors derived, however, also allows one to explore correlations of the renormalization factor with configurations. The configuration factor should represent different physics/configurations in subgroups. Furthermore, a practical demand for a configuration factor is its accessibility. The LHD configuration scan data clearly shows the importance of the magnetic configuration for confinement (cf. Fig. 4). Therefore, in a first step, properties of the magnetic field structure are explored with respect to the derived configuration factor.

Tokamak confinement scaling laws indicate significant dependence on the elongation [20]. Moreover, the ITER design [21] also includes the triangularity, indicating again the relevance of magnetic geometry for confinement. The purely geometrical definition of elongation used in tokamaks fails for stellarators being insufficient to characterize drift orbits in the helical field. Instead, a definition employing magnetic (Boozer) coordinates is used with generalizes the elongation in 3-d geometries (see [22]): 2 ε t κ = (5) b 10 where b 10 describes the average toroidal curvature in these coordinates. This definition is equivalent to the geometrical definition in the limiting 2D (tokamak) case. A number of caveats must be raised if this factor is used as a configuration factor for stellarators: The factor b 10 increases linearly with minor radius in classical stellarator configurations and this dependence should cancel the minor radius factor in the inverse aspect ratio ε t = r/r. However, there may be a substantial deviation from linear behavior near the last closed flux surface, for example, in the case of shifted heliotron configurations,. Therefore, we chose the κ value at 2/3 of the minor radius just as we choose for ISS95 and ISS04 to evaluate ι at the same radius.

The use of elongation must fail for quasi-helically symmetric configurations since the corresponding b 10 component vanishes in these configurations. At high, b 10 may change significantly. Further accessible parameters from the magnetic geometry are the plateau factor and the effective helical ripple [22]; a more detailed discussion on these parameters is given in Ref. [2]. The plateau factor is closely related to elongation but also includes helical terms in the mod B structure for plateautype scaling laws such as Lackner-Gottardi scaling [5]. The effective helical ripple is a neoclassical measure for the long-mean-free-path transport of trapped particles. All three configuration describing parameters account for the departure of drift from flux surfaces, but weight the mod B structure differently. Plotting accessible elongation, effective helical ripples and plateau factors vs. the configuration factor shows correlations (see Fig. 6). A more detailed view shows that, e.g., the radial electric field (e.g. high configuration factor values of W7-A in the upper plot) is to be considered and cannot be described by the elongation. Hence, an inclusion of more relevant transport effects is required. Additional information on the profiles of temperature and densities should be incorporated in order to assess the capabilities of the configuration factor candidates discussed here. 6. Conclusion Continuing analysis of the International Stellarator Confinement Database has allowed us to refine the conclusions of initial studies by examining the systematics of coherent datasets from the devices W7-AS and LHD.

Bayesian analysis of the W7-AS database confirms the validity of the Connor- Taylor electrostatic Fokker-Planck model for low- discharges. The situation is unclear for high data. The search for a physical mechanism behind the configuration factors determined in the multi-device ISS04 analysis will require more detailed consideration of profile and neo-classical transport effects, as it is difficult to decide between several possibilities (e.g. helical ripple, elongation etc.) on the basis of the global data included in the database so far. The ISCDB is an ongoing stellarator community effort. An essential goal aiming at comparative studies is the documentation of global confinement data. Beyond scaling studies, confinement data of further interest are reviewed as well, such as data from operation limits or high and parameter scans. For the ISS04 subset of data, the gross parameter scaling of previous studies is confirmed. A clear result is the impact of the magnetic configuration on confinement. Although this paper did not attempt to extrapolate the existing data to predict future performance, well documented data of subgroups from LHD and Wendelstein 7- AS exist which are compatible with the ITER predictions and typical ELMy H-mode tokamak confinement data in figures of stellarator scaling laws. For a comparison of envisioned tokamak and stellarator reactors, however, further discussion requires the consideration of differences in the operational parameters for burning plasmas, e.g., the density [23]. Beyond documentation and a continuous refinement of data and their quality, the techniques implemented allow for the next steps of ISCDB, i.e. tests of physical models on data and the refinement of possible candidates for a physical configuration factor

substituting the empirical configuration factor. Extensions of ISCDB with high quality data for scans in the regression variables and configurations are the key for further comparative assessments of the stellarator/heliotron concept. Acknowledgments We would like to thank the International Global H-Mode Confinement Database Working Group for data provision and O. Kardaun for extensive discussion and help in accessing the standard datasets for L- and H-mode. The authors appreciate discussions with V. Dose, Y. Feng, Yu. Igitkhanov, H. Maaßberg, D. Mikkelsen, J. Nührenberg, F. Sardei, Yu. Turkin and F. Wagner.

References [1] U. STROTH et al, Nucl. Fusion 36, 1063 (1996). [2] H. YAMADA et al, Nucl. Fusion 45, 1684 (2005). [3] The H-mode database is available at http://efdasql.ipp.mpg.de/hmodepublic/. The standard subset of DB3V10 was chosen for this paper. The L-mode database ITERLDB2v10_May10 is available at ftp://ftp.ipp.mpg.de (username: iter) and the standard data selection as used in A. KUS, O. KARDAUN, A. DINKLAGE, R. PREUSS, INT. H-MODE CONFINEMENT DATABASE GROUP, INT. STELLARATOR CONFINEMENT DATABASE GROUP: Empirical Scalings of Energy Confinement Times for Fusion Experimental Devices presented at the 25 th EMS Conference (Oslo, 2005) (available at http://www.ipp.mpg.de/~otto.kardaun/rep.html). [4] O. KARDAUN, Plasma Phys. Contr. Fusion 41, 429 (1999) [5] K. LACKNER, N. GOTTARDI, Nucl. Fusion 30, 767 (1990) [6] S. SUDO et al., Nucl. Fusion 30, 11 (1990) [7] R.E. WALTZ et al., Phys. Rev. Lett. 65, 2390 (1990) [8] P. GRIGULL et al., J. Nucl. Mat. 313, 1287 (2003) [9] Y. FENG et al., J. Nucl. Mat. 313, 857 (2003) [10] R. BRAKEL, Nucl. Fusion 42, 903 (2002) [11] E. ASCASIBAR, Nucl. Fusion 45, 276 (2005) [12] S. MURAKAMI et al, Nucl. Fusion 42, L19 (2002) [13] H. YAMADA et al, Plasma Phys. Contr. Fusion 43, A55 (2003) [14] J.W. CONNOR, J.B. TAYLOR, Nucl. Fusion 17, 1047 (1977) [15] R. PREUSS, V. DOSE: Errors in all Variables in BAYESIAN INFERENCE AND MAXIMUM ENTROPY METHODS IN SCIENCE AND ENGINEERING: 25th International Workshop on Bayesian Inference and Maximum Entropy Methods in

Science and Engineering, Eds. K.H. Knuth, A.E Abbas, R. D. Morris, J. P. Castle, AIP Conf. Proc. Vol. 803, (2005) pp.448-455 [16] R. PREUSS et al.: Stellarator Scaling Condidering Uncertainties in Machine Variables in Proc. 32th EPS Conf. on Plasma Physics and Controlled Fusion, Eds. C. Hidalgo and B. van Milligen, European Physical Society, ECA Vol. 29C, P-1.115, (2005) [17] V. DOSE, R. PREUSS, W. v.d. LINDEN, Phys. Rev. Lett. 81, 3407 (1998) [18] R. PREUSS et al., Nucl. Fusion 39, 849 (1999) [19] J.W. CONNOR, Plasma Phys. Contr. Fusion 30, 619 (1988) [20] ITER PHYSICS EXPERT GROUP et al., Nucl. Fusion 39, 2175 (1999) [21] ITER final design report (summary, document number G A0 FDR 4 01-07-21 R0.4), available at www.iter.org [22] E. RODRIGUEZ-SOLANO RIBEIRO, K.C. SHAING, Phys. Fluids 30, 462 (1987) [23] H. WOBIG et al., Nucl. Fusion 43, 889 (2003)

Footnote (a) The ISCDB resources are jointly hosted by the National Institute for Fusion Science (NIFS) in Japan and the Max-Planck-Institut für Plamaphysik (IPP) in Germany, see http://iscdb.nifs.ac.jp and http://www.ipp.mpg.de/iss

log 10 [τ E (s)] 10 1 10 0 10 1 10 2 ATF CHS Heliotron E Heliotron J HSX LHD TJ II Wendelstein 7 A Wendelstein 7 AS 10 3 10 4 H mode (DB3v10) L mode (DB2v10) ITER prediction (PPCF41,429(1999)) 10 4 10 3 10 2 10 1 10 0 10 1 log 10 [τ ISS95 (s)] MS318905 - Fig. 1

log 10 [τ E (s)] 10 1 10 0 10 1 10 2 ATF CHS Heliotron E Heliotron J HSX LHD TJ II Wendelstein 7 A Wendelstein 7 AS 10 3 10 4 10 4 10 3 10 2 10 1 10 0 10 1 log 10 [f τ ISS04 (s)] MS318905 - Fig. 2 (left)

ATF CHS Heliotron E Heliotron J HSX LHD (inward, oblate) LHD (inward, prolate) LHD (outward shifted, R = 3.9 m) LHD (standard, R = 3.75 m) LHD (inward shifted, R = 3.6 m) TJ II W7 A W7 AS (high ) W7 AS (high iota) W7 AS (low iota) MS318905 - Fig. 2 (right) 0 0.5 1 1.5 2 f= (τ exp /τ ISS95 )/1.63

25 20 W7 AS HDH Density Scans τ exp normal confinement τ exp HDH (attached) τ exp HDH (detached) f τ ISS04 (f = 1.61) 1.9 f τ ISS04 (f = 1.61) 15 τ E (ms) 10 5 0 0 1 2 3 4 n (m 3 ) x 10 20 MS318905 - Fig. 3 (left)

TJ II W7 AS Iota Variations (scaling agreement) 2 W7 AS ι scan [Brakel, NF 42, 903 (2002)] (f=0.60) TJ II (ISS04 data) (f=0.30) 1.8 ISS04 without ι dependence 1.6 1.4 τ ISS04 ι 0.41 a 2.28 τ exp / f τ ISS04 1.2 1 0.8 0.6 0.4 0.2 0 0 0.5 1 1.5 2 2.5 ι / 2π MS318905 - Fig. 3 (right)

2 1.5 τ E exp /τe ISS04 1 0.5 0 W7 AS L mode W7 AS H mode W7 AS normal confinement W7 AS HDH attached W7 AS HDH part. detached LHD outward LHD standard LHD inward prolate LHD inward oblate LHD inward MS318905 - Fig. 4

0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 low coll. low high coll. high Model probability 0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 low coll. low high coll. high Model probability 0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 low coll. low high coll. high Model probability 0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 low coll. low high coll. high Model probability MS318905 - Fig. 5

τ E exp /τ E ISS04 2 1 LHD R ax =3.6 LHD R ax =3.75 LHD R ax =3.9 LHD κ=0.8 LHD κ=1.38 W7-AS low β W7-AS high β TJ-II W7-A ATF 0 0 1 2 3 κ 2/3 MS318905 - Fig. 6 (left)

MS318905 - Fig. 6 (center)

MS318905 - Fig. 6 (right)

MS328905 Figure Captions Fig. 1: Experimental confinement data of stellarators vs. ISS95 scaling. The tokamak data [3] were rephrased employing the plasma elongation and assuming a typical q profile according to Ref. [1]. For the 1998 design of ITER [20], a prediction with uncertainty analysis [4] is plotted against the ISS95 scaling energy confinement value. Fig. 2: Experimental confinement data of stellarators vs. ISS04 scaling. The right panel is a box-plot of the configuration factors as used for ISS04 [2]. Fig. 3: Examples for parameter scans/variations documented in ISCDB. Left: Energy confinement time in the high-density H mode in Wendelstein 7-AS. Right: Comparison of iota scaling in W7-AS and TJ-II in figures of renormalized confinement times. The f values refer to the ISCDB data subset displayed. Fig. 4: Configuration factor for different confinement modes and configuration scans in stellarator/heliotron devices. Fig. 5: Results of Bayesian model comparison of two W7-AS subsets for low ι data. The normalization of collisionality (ν*) is chosen to be the transition of plateau to Pfirsch-Schlüter regime (tokamak like). Fig. 6: Correlations of elongation κ (left, at 2/3 of effective radius), effective helical ripple ε eff (center) and plateau factor Γ stell /Γ tokamak (right) vs. configuration factor [2].