Transcript ITER sim

Comprehensive ITER Approach to Burn

L. P. Ku, S. Jardin, C. Kessel, D. McCune Princeton Plasma Physics Laboratory SWIM Project Meeting Oct. 15-17, 2007 Oak Ridge National Laboratory

A comprehensive simulation of an ITER discharge has been completed. Plasma states obtained are being used for testing various aspects of IPS.

Using IPS/EPA we have carried out an ITER discharge simulation (a hybrid scenario) using NUBEAM and TORIC for the beam and RF calculations, respectively, via the alternate internal coupling method.

We use this comprehensive ITER approach to burn to

test coupling to sophisticated source models in TRANSP

– –

compare with previous calculations with analytic models provide plasma states for IPS to run in re-play mode for

component development and testing,

 

comparative studies for resolution requirement, and testing porting to new computer platforms of other components.

LPK-101507 2

Outline of Presentation

Approach to ITER discharge simulation

– – –

EPA in IPS EPA internal coupling to TRANSP (client-server) TSC & EPA

Physics models and results

– – –

discharge scenario source and current drive models summary of results

Improvements, targets and plans

LPK-101507 3

In normal mode of operation the equilibrium and profile advance (EPA) component communicates with others via the plasma state and the controller.

Plasma State Dist.

Function RF/Beam Source Electric Fields & Others Equilibrium & Profile Advance MHD Stability Controller

LPK-101507 4

However, there is an option in EPA which uses models in TRANSP to calculate the source and current drive in conjunction with the equilibrium and profile evolution.

Plasma State equilibrium profiles control advance free-bndy equilibrium sources EPA -- TSC advance profiles equilibrium profiles signal file passing Init, step, save, kill request sources TRANSP control receive ready, error evolve sources heating & currents NB, RF, fusion products

LPK-101507 5

The engine driving EPA is the Tokamak Simulation Code TSC.

TSC performs free-boundary self-consistent transport evolution in any part of a tokamak discharge.

TSC models PF coils and 2-D passive structures with circuit equations. Arbitrary power supply models can be used.

TSC runs with feedback control systems

radial position, vertical position, Ip, and shape

– –

sensors and reconstructions stored energy, current profile, etc.

LPK-101507 6

Using internal coupling to TRANSP initially facilitates the development of other components in ISP.

TRANSP has several high fidelity source models

– – – – –

NB NUBEAM (orbit following Monte Carlo) ICRF TORIC /SPRUCE (full wave), CURRAY (ray-tracing) LH - LSC (ray-tracing 1D Fokker Planck) EC - TORAY (ray-tracing, relativistic)

-particles (orbit following Monte Carlo)

TRANSP runs in “interpretive” mode for experimental analysis, but this same approach has been modified to run in “predictive” simulations by passing T(

), n(

), q(

), and equilibrium geometry.

LPK-101507 7

The case study is a “hybrid” scenario

Discharge scenario – off-axis beam and ICRF heating:

– –

3 s ≤ t ≤ 250 s pre-programmed: z eff , n e , I p ; flattop @ t~240 s

n e20 =0.924 @ t=240 s

Zeff=1.2 (t<50 s), 1.8 (t<150 s) , 2.2 (t>150 s)

Ip=12 MA @ t=240 s

– – –

total internal energy clamped at 450 MJ particle confinement time=25 s energy confinement

• •

t t< 50 s Coppi/Tang profile consistent model

50 s relaxed GLF23 and neo-classical

LPK-101507

Flux contours for the free-boundary equilibrium @250 s B=5.3 T, R=6.2 m, a=2.0 m,

e

=1.8,

d

=0.45

b=1.77%

,

b

p =0.74, q(0)=1.17, q(1)=4.30

8

Models for heating and current drive

– – –

Neutral beams

off-axis (beam tangency 5.295 m), 1 st beam on @ 10s, with partial power and @ 50s with full power, 2 nd beam on @ 110 s. Co, 1 MeV D, full power 16.5 MW.

NUBEAM, sample size=1000.

ICRF

1 antenna, 20 MW, 52.5 MHz on @ 170 s, energy.

3 He minority (0.001), power level regulated by the total internal stored

TORIC, 31 poloidal modes, 128 radial, 64 angular grids Alpha particles

Monte Carlo fast ion slowing down, sample size=1000.

Source integration time – 25 ms.

LPK-101507 9

LPK-101507

Time Evolution of Power Balance Beam 1 Total ICH+Beam Beam 1+2

 10

t=150 s

150000 100000 50000 0.0

0.2

0.4

0.6

0.8

1.0

150000 100000 50000 0.0

LPK-101507 0.2

0.4

rho

0.6

0.8

1.0

t=250 s

500000 400000 300000 200000 100000 0.0

0.2

0.4

0.6

0.8

1.0

150000 100000 50000 0.0

0.2

0.4

rho

0.6

0.8

1.0

11

LPK-101507

Time Evolution of Central Temperatures 1 st beam (8 MW) 1 st beam (5 MW) 1 st beam (16.5 MW) 2 nd beam (16.5 MW) ICH (20 MW) electron Decrease due to increased Zeff Ion

12

LPK-101507

Time Evolution of Current Drives Pre-programmed plasma current

13

Temperature, Density and q Profiles at t=250 s T i T e

LPK-101507 14

LPK-101507

n i n e

15

Some simulation results during the first 50 s of an ITER discharge with beam turned on at t=10 s. Note the loss of power for t<15 s before the plasma is fully developed.

Requested power

LPK-101507 16

Requests of Improvement, Targets and Plans

Simulation under full IPS control

Ease of use

readily accessible to users (Kessel, Budney, and others)

– •

having templates setup for machines (ITER, KSTAR, JET, NSTX, CMOD, DIII) Options for NBI and

-particle slowing down

NUBEAM, Bounce-averaged FP (ACCOME?), ASTRA package

Options for ICRF

TORIC, AORSA (full wave), GENRAY (ray-tracing)

Options for LHCD

LSC/ACCOME, GENRAY (1D FP), GENRAY/CQL3D (2D FP)

Options for EC

TORAY, GENRAY, GENRAY/CQL3D

Routine use of stability codes

PEST, BALLON, DCON

NOVA-K

Implement redistribution of fast ions

from sawtooth

from other MHD modes

LPK-101507 17

Improvements, Targets and Plans (cont.)

EPA enhancement

Add prediction of plasma rotation profile

requires sources from all sources

requires options for momentum transport

Improve density prediction

3 hydrogen species and impurity

pellet fueling model

– –

Include fast ion pressure in equilibrium evolution Implement new GLF23 (TGLF23, available in Dec) and standardized transport interface

needs to be in parallel (10-100 processors)

– – –

Implement options for edge pedestal model (including ELMs) Porting and benchmarking on Jaguar Implement interface with machine description files

Application of IPS/EPA to additional ITER discharge simulations

ELMy H-mode

Q(0) ≤ 1, useful for testing sawtooth handling in IPS

Steady-State

High

b

n, stability

Study sensitivity to resolution, calling frequency, etc.

LPK-101507 18