Report on IOS-JA4byJ.M. Park
Main Contributors:
M. Murakami, G. Giruzzi, J. Garcia, N. Hayasi, Y-S. Na, A. Polevoi,P. Bonoli, A. Fukuyama

Presented at the
7th ITPA-IOS Topical Group Meeting
Kyoto Univ, Kyoto, JapanOct 18–21, 2011

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 2011

ITER Steady-State Scenario Code-to-Code Benchmark• Published in Nuclear Fusion (M. Murakami et al, 2011 Nucl. Fusion 51103006)• Participating codes trough ITPA-IOS JA-4– FASTRAN/ONETWO, TRANSP, CRONOS, TOPICS-IB, ASTRAi, ASTRAk(ACCOME, TASK, …)• Generally, excellent agreement among the codes• But, also identified a number of small(?) differences, especially in– He transport– Rotation– H&CD modules (NB, EC, FW)

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 20112

Weak Magnetic Shear Scenario

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 20113

Plan for 2011 WAS to Resolve Differences We Found• He transport– Previous guideline: τp*/τE = 5⇒ benchmark with more specific He transport model (D, V, recycling, …)• Rotation– Most codes fixed the rotation profile (Ωφ=0)	

– No boundary condition provided yet⇒ a) Check momentum balance equation, GLF23 implementationb) Check beam-driven torque• EC– Benchmkark for ECH upgrade (wider range of steering angle)– Impact of momentum conservation, especially for ρEC > 0.5• FW– Most codes imported the guideline profile ⇒ need self-consistent modeling No much progress yet!!
J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 20114

Equation for toroidal rotation in variouscodes:
TRANSP [R J Goldston]:m =<R2>ωΣjnj*Mj, ω(√Φ) = Vϕ/R (sum over thermal ion species)

ONETWO and FASTRAN:

CRONOS [J F Artaud et al, NF 2010]:

ASTRA [G Pereverzev, P Yushmanov, IPP-2002]:
where F is specified by user. Torqueand losses should correspond to thechoice of F

JETTO solves the equation for Vtor
I. Voitsekhovitch, J. Garcia, J. Ferreira, D. Kalupin, X. Litaudo, J.M. Park	5

GLF23 equations for rotation[R. E. Waltz et al, Phys. Plasmas 4 (1997), 2482]

Γ is the ion particle flux

- torque from TRANSP to be recalculated to rotation source- χϕ  (dρ/dr)niηφeff- scenario with time evolving ion density ni(t,ρ) requires the changes ofequation for momentum implemented in transport codes	

I. Voitsekhovitch, J. Garcia, J. Ferreira, D. Kalupin, X. Litaudo, J.M. Park	

6

Case 1: computed momentum diffusivity with fixed
profiles	

ASTRA/CRONOS/FASTRAN(dashed)/JETTO/ONETWO(solid)ρI. Voitsekhovitch, J. Garcia, J. Ferreira, D. Kalupin, X. Litaudo, J.M. Park	7Benchmark between Full-wave AORSA and Ray-tracing
CURRAY codes
Table 3. Benchmark between full-wave AORSA and ray-tracing
CURRAY codes. CURRAY-1 includes beam and alpha absorption,
while CURRAY-2 does not.
P (electron) (%)
P (D-thermal) (%)
P (T-thermal) (%)
P (D-beam) (%)
P (Alpha) (%)
P (He4 ) (%)
P (impurity) (%)
I (FWCD) (kA)AORSA

CURRAY-1

CURRAY-2

64.3
0.4
32.8
1.0
1.5
0.0
0.0603

71.6
1.3
19.5
7.4
0.2
0.0
0645

79.3
0.2
20.30.2000654

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 2011

8

ITB Scenario

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 20119

Plan for 2011 WAS to Resolve Differences We Found• More participating codes/modelers (ASTRA, TRANSP,…)• Simplify the guideline for quick convergence• Employ theory-based transport model No much progress yet!!

J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 201110

Next Step?
• Close out JA-4?
– Mission completed??
– Documented well (Masanori’s NF paper)– Need to archive simulation datao ITER Profile database format
o Simply upload data files to the ITER ITPA homepage• Further utilize the guideline scenario to explore steady-stateoptimization?
– Now we have a good starting pointo Well benchmarked in terms of numerics and modelingprotocols– Sensitivities to the transport and pedestal models, density peaking(particle transport), …
– Optimization of SS scenario towards the steady-state Q=5 missiono Heating/CD upgrade study
J.M. Park, 7th ITPA-ISO Topical Group Meeting, Kyoto Univ., Japan, 201111