Toki 11-14 Dec. 2001

ITER CONFINEMENT AND STABILITY MODELING
Presented by A.R.POLEVOI, ITER JCT, Naka JWSOUTLINEBACKGROUND

MODEL DESCRIPTION (PRETOR & ASTRA)
HEAT, MOMENTUM AND PARTICLE TRANSPORTSAW TOOTH OSCILLATIONSHELIUM TRANSPORTBOUNDARY CONDITIONSMODEL VALIDATION:
ITER MODEL VS. EXPERIMENTAL DATAITER RELEVANT JET EXPERIMENTS

1.5D SENSITIVITY ANALYSES:INDUCTIVE SCENARIO:DIFFUSIVITY PROFILE
BOUNDARY CONDITIONSSAW TOOTH MIXING AND BALLOONING LIMITSTEADY STATE SCENARIO:MHD MODE STABILISATION BY CONDUCTING WALL

CONCLUSIONS

1

Toki 11-14 Dec. 2001

BACKGROUND
Comparative analysis of the different transport models available atpresent reveals [1] that no single model agrees well with theequilibrium temperature profiles measured in experiment. It was
found [2,3] that the models predict significantly different transportlevels for the same instabilities governing the radial transport.
3
2.5

Theory based
Empirical
Theory Based (Gyro-Fluid)

2
1.5
1
0.5
0
IFSPPPL IFSPPPL GLF 23
no ExBExB

CulhamMixed- Multi-Mode CDBM
Shear

RLW

Fusion power predicted by various models for ITER-98 [3].
So, at present it is premature to recommend any single 1D modelfor the ITER predictive analysis. And we use semi empiricalapproach based on energy confinement scaling.
Benchmarking of the ITER transport model vs. experimental data
and robustness of chosen scenarios vs. model assumptions andphysical limitations are subjects of our consideration.
[1] DeBoo J C et al 1999 Nucl. Fusion 39 1935
[2] Dimits A M et al 2000 Phys. Plasmas 7 969[3] ITER Physics Basis 1999 Nucl. Fusion 39 2178

2

Toki 11-14 Dec. 2001

ITERH-98P(y,2) Global Confinement Scaling
The Confinement Database and Modelling Expert Group recommended for
ITER design the ITERH-98P(y,2) confinement scaling,

H98(y,2)tE= 0.0562I

0.93 0.15 0.41 -0.69 1.97 0.78 0.58 0.19
B n19 PR ka e M

The point prediction for the thermal energy confinement time in ITER is tE =3.6 s.
The 2s log-linear interval was determined as ±20%.
By recent analysing the enlarged ITERH.DB3 (‘final’) dataset, the practicalreliability of the ITERH-98(y,2) scaling was confirmed and 2s log-linear interval
was reduced1
to ±14%.

10

1

t exp
E ,s 0.1

ASDEX 431
AUG
468
CMOD
45
COMPAS 16
DIII-D 300
JET1114
JFT-2M 59
JT-60U 54
PBX-M
59
PDX
97
TCV
11
TFTR
13
TdeV
3
START
8

ITER­

0.01
r.m.s.e.: 16%

0.001
0.001

0.01

0.1

t

1

1

10

H98(y,2),s

O. Kardaun, “Interval estimate of the global confinement time during ELMy H-mode
in ITER FEAT, based on the international multi-tokamak ITERH.DB3 dataset,” IPPIR 2001/5 1.1 http://www.ipp.mpg.de/ipp/netreports, in preparation; O. Kardaun, “Onestimating the epistemic probability to realise Q=Pfus/Paux larger than specified lower
bound in ITER,” Nucl. Fusion (2001), in press.
3

Toki 11-14 Dec. 2001

COMPARISON OF PREDICTIONS BY PHYSICS BASEDMODELS WITH ITER DESIGN POINT

[1]
[2]
[3]

A.H.Kritz, J.Kinsey,T.Onjun, I.Voitsekhovich, G.Bateman, R.Waltz, G.Staebler, "Burning
Plasma Projections with Internal Transport Barriers", ITPA Meeting on Burning Plasma
Transport, 10-12 September 2001, NIFS, Toki, Japan.
J.Weiland, "Predictive Simulations of ITER-FEAT Performance," 28th EPS Conference,
Madeira, 2001, P2.039.
G. Bateman, A. H. Kritz, T.Onjun and A. Pankin, Private communication, 7 Dec., 2001.

4

Toki 11-14 Dec. 2001

GLF23 MODEL PREDICTION FOR ITER

Temperature at the top of pedestal required for obtaining Q=10
drops with reducing Paux

5

Toki 11-14 Dec. 2001

MODEL DESCRIPTION(PRETOR & ASTRA)

HEAT, MOMENTUM AND PARTICLE TRANSPORT
1D transport modelling for Te, Ti, ne, nHe, y, Vtor evolution with
self-consistent 2D equilibrium (1.5D modelling)
Heat, toroidal momentum and particle diffusivities:
ci,e,f = ci,e,f(0)f (r) h (r)+ (1-h(r)) cneo,

cI,f (0)= 2 ce(0),D = D(0) f (r) h (r)+ (1-h(r)) cneo

D(0) = ce(0)Neoclassical edge pedestal transport cneo:
h (r) = 1

(r < 0.9),

h (r) = 0(r ≥ 0.9),

where r is the square root of the normalised toroidal flux.
Profile dependence, used for ITER simulations:
ASTRA:
PRETOR:f (r) = 1 + 3r2
f (r) by Rebut-Lallia-Watkins-Boucher [1]Semi empirical approach:ce(0) is fitted to provide the scaling dependence at the proper
phase of the discharge [2]:HH98 = tE/tE,H98(y,2) = 1

tE,H98(y,2) = 0.0562 I0.93B0.15n0.41R1.39P-0.69k0.78a0.58M0.19
Plasma heating and current drive, plasma fuelling by gas puffing,pellets and neutral beam are also simulated. Impurities (other thanHe) are prescribed as nzk = fkne, the fuel densities nD, nT are
calculated from the quasineutrality conditions:
ne=nD+nT+2 nHe+SkZknzk.
[1] P. H. Rebut, et al., in Proc. 12th Plasma Physics and Controlled Nuclear Fusion, Nice, 1988
(IAEA, Vienna, 1989) p.191[2] ITER PHYSICS BASIS, Nucl. Fusion, 39, 1999

6

Toki 11-14 Dec. 2001

SAW-TOOTH MIXINGFor ITER simulations two approaches were used:
ASTRA:

Semi-empirical approach is chosen:
ST Trigger: q = 1 at any radial position somewhere;
ST mixing width: r < 1.4r (q=1);

PRETOR: Complete reconnection model by F. Porcelli et. al [1]ST Trigger: dWmag > Wthr (perturbed magnetic energy > threshold)ST mixing width: rST is calculated from flux continuity ;In both approaches particles and temperatures are flattened over
the ST zone taking account of particle and energy conservation.
Pressure profile recovers faster than plasma current profile:tST >> tE. So, the details of the ST modelling have minor effect on
plasma performance.
T(0)

tE

tST
t
[1] F. Porcelli, D. Boucher and M. N. Rosenbluth, Plasma Phys. Control. Fusion 38 (1996) 2163

7

Toki 11-14 Dec. 2001

HELIUM TRANSPORT
Two approaches are tested:1.

No neutral He influx at the core boundary, He pumping at theboundary (t*He/tE) is controlled independently. The reference
operational point is chosen to be t*He/tE = 5, wheret*He = ònHedV/GHe,fus, GHe,fus = òSHedV is the fusion He source.
2.

By parameterisation of B2Eirene [1] calculation for SOL/DIV
we calculate boundary He density nHe(ra) and He atomic influx
GHe,atom self-consistently with core/SOL/DIV parameters.
Operational point t*He/tE is calculated.

nHe,2
nHe,1

GHe,fus0GHe,atom
1

X

Qualitative behaviour of He density profiles nHe,1,2 with
GHe,1 = GHe,fus,

GHe,2 = GHe,fus + GHe,atomHigher fuel dilution by He is expected for type 2 approach[1] A. S. Kukushkin, et al., “Basic Divertor Operation in ITER-FEAT”, 18th IAEA Fusion
Conference, Sorrento, Italy, Oct. 2000

8

Toki 11-14 Dec. 2001

BOUNDARY CONDITIONS
We consider separatrix as a boundary of plasma core 1.5 transport
analysis.
Two approaches are used:
1.

PRETOR Simplified SOL/DIV model, which gives
relatively high boundary density ne ~ 6 1019m-3 andtemperature T ~ 1 keV.

2.

Analytic interpolation of self-consistent B2-Eirene
SOL/Div simulations [1] for core boundary conditions.
This interpolation calculates the boundary conditions asfunctions of the pumping speed and particle circulationand heat/particle loss to the SOL/divertor region.
For the reference Pfus = 400 MW inductive operation it
gives lower boundary density ne(a) ~ 3 1019m-3 andtemperature T(a) ~ 200 eV for loss power < 100 MW.

9

Toki 11-14 Dec. 2001

MODEL VALIDATION:
MODEL PREDICTIONS VS. EXPERIMENTAL PROFILE DB1 [1]MODEL

DTe,std,%

DTi, std,%

Weiland

18

23

Multi-mode

13

15

GLF23

25

24

IPS/PPPL

24

16

CDBM

35

36

RLWB

20

21

Culham24

22

Mixed-shear

18

33

T11/SET

14

18

ITER2

14

12

Mean standard deviation DT,std = ( S (Ts–Tx)2)/( S Tx2))1/2, where Ts is1simulation, Tx is an experiment.
2

ITER model used experimental boundary conditions, other models

start from the top of the edge pedestal. ITER model is applied toreduced set of data (HH98 ~ 1, high density n/nGW > 0.5 with flat
density profile). So, direct comparison with other models is not
appropriate. It is presented just for scale to conclude, that
Semi-empirical model used for ITER predictions satisfactorily
reproduces experimental profiles from the profile database.
[1] ITER Physics Basis 1999 Nucl. Fusion 39 2178

10

Toki 11-14 Dec. 2001

ITER MODEL BENCHMARKING VS. PROFILE DATABASE.

Temperature Ti,Te and density ne profiles simulated by semiempirical model. Experimental profiles are shown by crosses. Input
power profiles and boundary conditions are taken from theexperimental DB.
JET#35156 H-mode:

n/nGW = 0.84

QNBI = 9 MW

Semi-empirical model used for ITER predictions satisfactorily
reproduces experimental profiles from the profile database.

11

Toki 11-14 Dec. 2001

ITER OPERATIONAL SPACE FROM 1.5D SIMULATIONS
Operational Domain for IP = 15 MA and Q = 10

Fusion Power (MW)

1000
800

1.5D simulation (PRETOR)
IP=15 MA, Q=10<ne>/nG = 1.0<ne>/nG = 0.85
<ne>/nG = 0.70600

bN = 2.5
bN = 2.0

400

PLOSS/PLH = 1.3

200
00.85 0.90 0.95 1.00 1.05 1.10 1.15
HH98(y,2)Operation boundaries (shaded area) are given by<ne>/nG = 1.0, bN = 2.5 and PLOSS/PL-H = 1.
REFERENCE OPERATIONAL POINT:
HH98(y,2) = 1.0 and <ne>/nG = 0.85, Pfus = 400 MW.
Calculations [1] were carried out by PRETOR code [2]bN is the normalised beta,nG is the Greenwald density limit,
PLOSS is the power loss,
PLH is the power required for the H-mode transition [3].
Good confinement is supposed for PLOSS > 1.3 PL-H.
[1] Y. Murakami, et al., J. Plasma and Fusion Res. 77 (2001) 712[2] D. Boucher, et. al., in Proc. 16h IAEA Fusion Energy Conference,Montreal, 1996 (IAEA,Vienna, 1997) 945.
[3] ITER Physics Basis, Nucl. Fusion 39 (1999) 2137.

12

Toki 11-14 Dec. 2001

ITER RELEVANT JET EXPERIMENTS

Experiments on JET tokamak demonstrated possibility of highconfinement HH98=1 in the configuration, similar to ITER withplasma density close to Greenwald limit n/nGW ~ 1.
[1] J.Ongena et. al. EPS-28, Madeira, 2001

13

Toki 11-14 Dec. 2001

ITER RELEVANT JET EXPERIMENTS

Experiments on JET tokamak demonstrated possibility of highconfinement HH98=1 in the configuration, similar to ITER withplasma density close to Greenwald limit n/nGW ~ 1.
[1] J.Ongena et. al. EPS-28, Madeira, 2001

14

Toki 11-14 Dec. 2001

2.0

0.0
25

5
0
2.0

PRETOR
ASTRA

1.5
j OH

1.0
0.5

j NB

PRETOR
ASTRA

20

Te (keV)

10

15
10
5
0
2.0

0.0
0.5

PRETOR
ASTRA

1.5
jT
1.0
0.5

j BS

0.0
2.0
PRETOR
ASTRA

0.4
0.3
Pa

0.2
0.1

3.0

1.0

j T, jBS (MA/m2 )15

PRETOR
ASTRA

4.0

PADD

0.0
0.0 0.2 0.4 0.6 0.8 1.0

ci (m 2/s)

T i (keV)j OH, jNB (MA/m2)

5.0

PRETOR
ASTRA

20

P a, PADD (MW/m3 )DIFFUSIVITY PROFILEqY

ne, nDT (1019 /m3 )1.5D SENSITIVITY ANALYSES:
12
ne
10
8
nDT
6
4
PRETOR
2
ASTRA
0
25

1.5

PRETOR
ASTRA

1.0
0.5
0.0
0.0 0.2 0.4 0.6 0.8 1.0

rr

ASTRA and PRETOR give similar profiles for the same boundary conditions
Central zone r < 0.5 transport seems essential cPRETOR ~ cASTRA

15

Toki 11-14 Dec. 2001

DIFFUSIVITY PROFILE

1.5D SENSITIVITY ANALYSES:

R/a (m/m)
BT (T)Ip(MA)
k95 / d95<ne> (1019m-3)
n/nG<Ti > (keV)<Te> (keV)
bT
(%)
bN
PFUS (MW)PNB (MW)PRF (MW)Q = PFUS / (PNB + PRF)
Wth (MJ)
PLOSS / P L-HtE (s)
fHe,ave (%)
Zeff, ave
P RAD (MW)li (3)ICD / IP (%)
IBS/ IP (%)
ci/ce
HH98 (y,2)
tHe*/ tE

ASTRA
6.2 / 2.0
5.3
15.0
1.7 / 0.33
10.1
0.85
8.3
8.9
2.53
1.78
394
33
7
9.84
327
1.8
3.73
3.1
1.66
43
0.84
7.9
16
2.0
1.0
5.0

PRETOR
¬
¬
¬
¬
¬
¬
8.0
8.8
2.49
1.76
400
¬
¬
10
320
1.8
3.71
3.2
1.66
47
0.85
7.6
15
¬
¬
¬

Semi empirical approach predicts for ITER similar results for similar boundary
conditions but different diffusivity profiles.

16

Toki 11-14 Dec. 2001

qY

ne, nDT (1019 /m3)

1.5D SENSITIVITY ANALYSES: BOUNDARY CONDITIONS
5.0
12
ne
PRETOR
10
4.0
ASTRA
8
nDT
3.0
6
2.0
4
PRETOR
1.0
ASTRA
2
0
0.0
40
40
PRETORASTRA

Te (keV)

100
2.0

PRETOR
ASTRA

1.5
j OH
1.0
0.5
0.0
0.6
0.5
0.4
0.3
0.2
0.1

j NB

PADD

20
10

PRETOR
ASTRA

jT
1.5
1.0
0.5

j BS

0.0
2.0

PRETORASTRA

Pa

PRETOR
ASTRA

30

0
2.0
j T, jBS (MA/m2)

20

ci (m 2/s)

P a, PADD (MW/m3 )j OH, jNB (MA/m2)

T i (keV)

30

0.0
0.0 0.2 0.4 0.6 0.8 1.0

1.5

PRETOR
ASTRA

1.0
0.5
0.0
0.0 0.2 0.4 0.6 0.8 1.0rr

Normalization HH98 =1 gives higher central temperatures (ASTRA) for lower
boundary values. Central zone r < 0.5 transport seems essential cPRETOR >cASTRA

17

Toki 11-14 Dec. 2001

1.5D SENSITIVITY ANALYSES:

R/a (m/m)
BT (T)Ip(MA)
kX / dX<ne> (1019m-3)
n/nG<Ti > (keV)<Te> (keV)
bT
(%)
bN
PFUS (MW)PNB (MW)PRF (MW)Q = PFUS / (PNB + PRF)
Wth (MJ)
PLOSS / P L-HtE (s)
fHe,ave (%)
Zeff, ave
P RAD (MW)li (3)ICD / IP (%)
IBS/ IP (%)
ci/ce
HH98 (y,2)
tHe*/ tE

BOUNDARY CONDITIONS

ASTRA
6.2 / 2.0
5.3
15.0
1.7 / 0.33
10.1
0.84
8.9
9.7
2.78
1.97
471
33
7
11.8
348
1.9
3.45
4.7
1.69
46
0.81
8.9
22
2.0
1.0
6.7

PRETOR
¬
¬
¬
¬
¬
0.85
8.0
8.8
2.49
1.76
400
¬
¬
10
320
1.8
3.71
3.2
1.66
47
0.85
7.6
15
¬
¬
5.0

Semi empirical approach predicts enhanced performance for B2-Eirine
boundary conditions.

18

Toki 11-14 Dec. 2001

1.5D SENSITIVITY ANALYSES:SAW-TOOTH MIXING ANDBALLOONING LIMITBACKGROUND FOR ANALYSIS
Saw Tooth mixing zone is large r/ra ~ 0.5 with high pressure
gradient p¢ before the saw-tooth and low magnetic shear s =rq¢/q << 1
Mercier criterion gives the limit for ballooning stable pressure
gradient s2 > - 8 rm0p¢ (1-q2)/B2
ANALYSIS OF IDEAL MHD MODE STABILITY
Is carried out by KINX code [1] coupled with ASTRA

Pressure gradient is close to ballooning/Mercier stability limit p¢ ~p¢lim in the saw-tooth mixing zone q ~ 1 before the saw-tooth in the
inductive scenario.
[1] Degtyarev L, Martynov A, Medvedev S, Troyon F, Villard L, Gruber R 1997 Comput. Phys.
Comm. 103 10

19

Toki 11-14 Dec. 2001

STEADY-STATE HIGH-Q OPRATION
Steady state operation with high Q > 5 requires high beta bN > 4 li
operation where the Ideal kink modes become unstable.
High bootstrap current fraction in the SS operation produces
reversed shear configurations q(0) > qmin.
a

b
2 -1

19

D, m sT, keVn, 10 m

2

40
35

T

30
25

T

-3

q

10

D

10

n

8

8

e

i

6

6

e

n

1

20

D+T

4

15
10

4

q

2

2

5
0

0
0

0.2

0.4

0.6

0.8

0

1

x

a)
b)

0

0.2

0.4

x

0.6

0.8

1

0

Te, TI and diffusivity D radial distributions
safety factor q, electron and fuel dencities ne, nD+T.

SS is possible for the same global parameters: geometry, B, Ip, newith different safety factor profiles (different qmin) and multiplication
factor Q (which decreases when qmin increases).

20

Toki 11-14 Dec. 2001

STEADY-STATE HIGH-Q OPRATION
SS is possible for the same global parameters: geometry, B, Ip, newith different safety factor profiles (different qmin) and multiplication
factor Q (which decreases when qmin increases).

Maximal distance to the conducting wall, which can provide thekink mode stabilisation, increases with qmin. So, the ITER design
wall position implies lower limit on qmin (upper limit on Q) for
chosen global parameters.
ITER plasma parameters for the SS WNS scenario different q
profiles.
Parameter

Value

Parameter

Value

R/a, m

6.35/1.85<ne>,1019m-3

6.74

d95/k95

0.41/1.84<Te>n/<Ti>n, keV 11/12-10.5/11

q95

5.16-5.13

Wth/Wfast , MJ

273/60-255/50

qmin

2.1-2.4

HH98

1.41-1.3*bN

2.8-2.56

Q

5.7-5

li

0.72-0.63

PNB/PLH, MW

34/29-33.7<Zeff>

2.2-2.17

Pfus/Ps, MW

361/93-338/97

Ip, MA

9tE, s

2.54-2.32

BR, Tm

32.86

n/nG

0.83

*

In SS simulations we used the neoclassical ion heat diffusivity inthe reversed shear zone.

21

Toki 11-14 Dec. 2001

IDEAL n=1 KINK MODE STABILISATION BY CONDUCTING WALL

STABILISING WALL POSITION aw/a VS. NORMALISED BETA bN FOR
DIFFERENT SCENARIO
No-wall limit is shown by dashed lines
Pressure scan is carried out for fixed q profile with different qmin
Maximal distance to the conducting wall that can provide the kinkmode stabilisation increases with qmin.
The ITER design wall position, aw/a ≈ 1.4, implies lower limit on qmin
(upper limit on Q) for chosen global parameters.

22

Toki 11-14 Dec. 2001

CONCLUSIONS
1.

Semi-empirical

model

used

for

ITER

predictions

satisfactorily

reproduces experimental profiles from the profile database.
2.

Semi-empirical approach predicts for ITER weak sensitivity of plasma
performance to diffusivity profiles for similar boundary conditions.
Simulation predicts enhanced performance for B2-Eirene compatible
boundary conditions.

3.

The details of saw-tooth modelling do not affect plasma performanceprovided the size of the mixing zone is similar and ST period is higher
than the pressure recovery time. The pressure gradient in the mixing
zone is marginally stable vs. ballooning modes in the reference
inductive scenario.

4.

High Q > 5 steady state operation would require stabilisation of low-n
ideal kink modes. There is an operational window for the stabilising wall
position compatible with the ITER design aw/a > 1.4.
SO, THE REFERENCE ITER SCENARIOS ARE ROBUST AGAINST THE
CONSIDERED EFFECTS.

23