Barbora Gulejová 1 of 10PSI 2008 abstract 5/12/2007

SOLPS modelling of Type I ELMing H-mode on JET

Barbora Gulejová, Richard Pitts, David Coster, Xavier Bonnin,

Marc Beurskens, Stefan Jachmich, Arne Kallenbach, Janos Marki

Rehearsal of abstract for PSI 2008

OUTLINE

* Type I ELMing H-mode on JET * Benchmarking code – experiment * Benchmarking code - code

Barbora Gulejová 2 of 10PSI 2008 abstract 5/12/2007

Type I ELMing H-mode on JET

# 70224

Wdia

PIN ~ 20 MW PRAD (below Xpoint)

PRAD (core)

Dalpha Parameters

Bt = 3 T Ip = 3 MA

ne=6x1019 m-3

P(SOL)= 19 MW

No GAS PUFF !!

ν*ped ~0.03 -0.08 (expected for ITER!)

ELM parametersfelm ~ 2 HzΔWELM ~ 1 MJ

Simulated for the first time ! … with SOLPS5

Core LIDAREdge LIDAR Li beam HRTS ECE

Langmuir probes

time

=baseline scenario for QDT=10 burning plasma operation on ITER !!! To avoid divertor damage => maximum ΔWELM (ITER) ~ 1MJ = this is achievable on JET

Bolometry (A.Huber)radiation between ELMs

Barbora Gulejová 3 of 10PSI 2008 abstract 5/12/2007

SScrape-crape-OOff ff LLayer ayer PPlasma lasma SSimulationimulation Suite of codes to simulate transport in edge plasma of tokamaks

B2B2 - solves 2D multi-species fluid equations on a grid given from magnetic equilibrium

EIRENE EIRENE - kinetic transport code for neutrals based on

Monte - Carlo algorithm

SOLPS 5SOLPS 5 – coupled EIRENE + B2.5

Main inputs: magnetic equilibrium Psol = Pheat – Prad

core upstream separatrix density ne

Free parameters: cross-field transport coefficients (D┴, ┴, v┴)

B2 plasma background =>recycling fluxes

EIRENE

Sources and sinks due to neutrals and molecules

measured

systematicallyadjusted

Barbora Gulejová 4 of 10PSI 2008 abstract 5/12/2007

Simulation vs. experiment - upstream

D [m2.s-1]Χ [m2.s-1]

ne [m-3]

Te [eV]

Ti [eV]

vperp [m.s-1]

inward pinch!

Core LIDAREdge LIDAR Li beam HRTS ECE

Barbora Gulejová 5 of 10PSI 2008 abstract 5/12/2007

Simulation vs. experiment - targets (LP)Outer targetInner target SOLPS LP LP

r-rsep r-rsep

r-rsep r-rsep

8050

0.5

50

1,2

30

Barbora Gulejová 6 of 10PSI 2008 abstract 5/12/2007

Benchmarking code-code

SOLPS 5

B2.5 + EIRENE

fluid (Braginskii) kinetic (Monte Carlo), neutrals

EDGE2D + NIMBUS

fluid (Braginskii) kinetic (Monte Carlo), neutrals

(less complex then EIRENE!)

vs.

Barbora Gulejová 7 of 10PSI 2008 abstract 5/12/2007

Type I ELMing H-mode on JET# 58569 Succesfully modelled by EDGE2D + NIMBUS by Arne Kallenbach (PPCF 46,2004)

Core LIDAREdge LIDAR Li beam ECE

Parameters

Bt = 2 T Ip = 2 MA

ne= 4x1019 m-3

P(SOL)= 12 MW

GAS PUFF from inner divertor !!

ELM parametersfelm ~ 30 HzΔWELM ~ 200 kJ

PIN ~ 14 MW PRAD (core)

Dalpha

Wdia

time

Barbora Gulejová 8 of 10PSI 2008 abstract 5/12/2007

Simulation vs. experiment-upstreamEDGE2D+NIMBUS

SOLPS5

-0.1 r-rsep [m]

0

22

0.04

Same AnsatzD,Chi,v

Same ne,Te,Tiupstreamprofiles

Barbora Gulejová 9 of 10PSI 2008 abstract 5/12/2007

Simulation vs. experiment-targets(LP)EDGE2D+NIMBUSSOLPS5

inner target outer target outer target

LP

250300

25 25

7

5

Barbora Gulejová 10 of 10PSI 2008 abstract 5/12/2007

Thank you for attention