1 / 22

Progress in ITER relevant ELM exhaust physics at JET

Progress in ITER relevant ELM exhaust physics at JET. Presented by W.Fundamenski on behalf of JET Task Force E and JET EFDA Contributors Part of an invited talk by R. A. Pitts at the 49 th Annual Meeting of the APS-DPP, Orlando, Florida, US, 12-16 November 2007.

cybil
Download Presentation

Progress in ITER relevant ELM exhaust physics at JET

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Progress in ITER relevant ELM exhaust physicsat JET Presented by W.Fundamenski on behalf of JET Task Force E and JET EFDA Contributors Part of an invited talk by R. A. Pitts at the49th Annual Meeting of the APS-DPP, Orlando, Florida, US, 12-16 November 2007

  2. with thanks to many co-authors A. Alonso1, P. Andrew2, G. Arnoux3, S. Brezinsek4, M. Beurskens5, J. P. Coad5, T. Eich6, G. Esser4, W. Fundamenski5, A. Huber4, S. Grünhagen7, B. Gulejova8, S. Jachmich9, M. Jakubowski10, A. Kirschner4, S. Knipe5, A. Kreter4, T. Loarer3, J. Likonen11, A. Loarte12, E. de la Luna1, J. Marki8, M. Maslov8, G. F. Matthews5, V. Philipps4, M. Rubel13, E. Solano1, M. F. Stamp5, J. D. Strachan14, D. Tskhakaya15, A. Widdowson5 and JET EFDA Contributors* 1Associacion Euratom/CIEMAT para Fusion, Madrid, Spain 2ITER Organization, Cadarache, France, 3Association EURATOM-CEA, DSM-DRFC, CEA Cadarache, 13108 Saint Paul lez Durance, France 4Institut für Plasmaphysik, Forschungszentrum Jülich GmbH, EURATOM Association, Trilateral Euregio Cluster, D-52425 Jülich, Germany 5Euratom/UKAEA Fusion Association, Culham Science Centre, Abingdon, OX14 3DB, UK 6Max-Planck-Institut für Plasmaphysik, IPP-EURATOM Association, D-85748 Garching, Germany 7FZ Karlsruhe, Postfach 3640, D-76021 Karlsruhe, Germany 8CRPP-EPFL, Switzerland, Association EURATOM-Swiss Confederation 9LPP, ERM/KMS, Association Euratom-Belgian State, B-1000, Brussels, Belgium 10Max-Planck-Institut für Plasmaphysik, Teilinstitut Greifswald, Germany 11VTT Technical research Centre of Finland, Association EURATOM-Tekes, Finland 12EFDA-Close Support Unit, Garching, Boltzmannstrasse 2, D-85748 Garching bei München, Germany 13Association EURATOM-VR, Fusion Plasma Physics, Stockholm, Sweden 14PPPL Princeton University, Princeton, NJ 0854, USA 15University of Innsbruck, Institute for Theoretical Physics, Association EURATOM-ÖAW, A-6020 Innsbruck, Austria *See appendix of M. Watkins et al., Fusion Energy 2006 (Proc. 21st Int. Conf. Chengdu, 2006) IAEA Vienna (2006)

  3. The problem with ELMs Important also in preparation for JET ITER-like wall and improved understanding of ELM SOL physics Material damage poses a limit on the maximum ELM size tolerable on ITER Current estimates indicate that ELM power fluxes (for CFC or W) must remain below ~0.5 MJm-2 at the ITER divertor targets JET Type I ELMs can approach 1 MJ  study the effects on first wall surfaces and edge plasma This implies an ELM energy loss, DWELM~ 1 MJ ~0.3% of stored energy in ITER QDT = 10 burning plasma! This is lower than any ELM energy so far achieved  mitigation strategies required. BUT …

  4. Large ELMs with low fueling #70226 – no gas fuelling Da (inner) PTOT(MW) WDIA(MJ) Te,ped (keV) ne,ped(1019m-3) H98Y Zeff(Brems) Time (s) Mostly NBI Vertical targets, MarkIIHD div.Specific JET sessionIp = 3.0MA, Bj = 3.0T, gas scanq95 ~ 3.1, d95 ~ 0.25Input energy ~195 MJEnergy Tile 3,7: 24.6, 70.1 MJ R. A. Pitts et al., ITPA, Garching, 2007

  5. Large ELMs with low fueling #70226 – no gas fuelling Da (inner) PTOT(MW) WDIA(MJ) Te,ped (keV) ne,ped(1019m-3) H98Y Zeff(Brems) Time (s) ITER Lowest fuelling cases at ITER relevant n*ped WELM/Wped ~ 0.2 for largest ELMs R. A. Pitts et al., ITPA, Garching, 2007

  6. Target surface temperatures ~600ºC ~200ºC #70228 – no gas fuelling Target surface temperatures from tangential view. Time resolution insufficient for power flux analysis Da (inner) #70228 Inner Outer Total wetted area ~1.0 m2 (cf. ITER ~ 3.5 m2) Tmax inner (ºC) Inter-ELM power loads higher at outer than inner as usualClear affect of surface layers on inner target (none on outer)Large ELMs: DTsurf (inner) ~ 600ºCDTsurf (outer) ~ 200ºCTsurf far from bulk sublimation Tmax outer (ºC) Time (s) J. Marki, T. Eich

  7. Radiation during large ELMs 0.85 MJ 1.29 MJ 1.08 MJ 0.58 MJ #70225, low fuelling A. Huber et al., EPS 2007 Da(inner) WDIA (MJ) PRAD (MW) Erad (MJ) Strong in-out asymmetry in ELM induced radiation for high DWELM  probably due to layers on inner targets and preferential inboard deposition of ELM energy Time (s)

  8. In-out ELM radiation asymmetry DERAD/DWELM ~ 0.5 if DWELM 0.6 MJ Evidence for a break at largerDWELM > ~ For DWELM 0.6 MJ radiation “spills over” separatrix – in-out radiation asymmetry reduced < ~ First ELM spikeonly Up to 70% DWELM radiated WELM = 0.45 MJ WELM = 0.85 MJ R. A. Pitts, ITPA 2007, A. Huber et al., EPS 2007

  9. ELMs can move carbon Non-linear dependence of carbon erosion on ELM energy  thermal decomposition of surface layers and favourable geometry rapidly increases QMB deposition Explains high deposition rates on water-cooled louvres during 1997 JET DT experiments  high T-retention A. Kreter, H. G. Esser et al., submitted to PRL

  10. Main wall ELM filaments ELM exposure superimposed on ambient background Difference frame: ELM – previous ELM-free frames #66515 DWELM ~ 200 kJt = 7.6 sExp. time 300 msFrame time 7.8 ms New wide angle IR camera diagnostic (E. Gauthier et al., CEA) using ITER-like front mirrors. 640x512 pixel FPA, max. full frame rate 100 Hz W. Fundamenski, M. Jakubowski, ITPA Garching May 2007, P. Andrew et al., EPS 2007

  11. Filament footprint field aligned Df = 22o,35o 68193, 57 s Filament IR footprint in main chamber closely aligned to pre-ELM field linesMode number (in this case) n ~360/Dj = 11-16. More cases  n = 10 - 50 Field aligned filaments also seen at upper dump plates: crude mode analysis gives n ~ 5 - 20

  12. How much ELM energy to walls? Main chamber IR camera too slow to follow single ELMs and filaments very asymmetric toroidally and poloidally Make energy balance for a single outboard poloidal limiter during H-mode phase, assume:Only ELMs can deposit energy on limitersNo energy to upper dump platesNo energy deposited in compound phasesSame energy on 16 limiters 68193, 57 s

  13. How much ELM energy to walls? 17.405 s 20.016 s #70226 Da (inner) Temp. (ºC) 13 12 11 Energy per tile (kJ) ∑Etile (15 tiles) Time (s) Main chamber IR camera too slow to follow single ELMs and filaments very asymmetric toroidally and poloidally Make energy balance for a single outboard poloidal limiter during H-mode phase, assume:Only ELMs can deposit energy on limitersNo energy to upper dump platesNo energy deposited in compound phasesSame energy on 16 limiters 68193, 57 s

  14. Wall loading and ELM size Ip = 3.0 MA, Bj = 3.0 T, gas scan. Separatrix-midplane outer wall gap fixed at ~5.0 cm. DWELM estimated for first ELM peak only 68193, 57 s Larger ELMs deposit more energy on outboard main chamber surfaces. How does this compare with theory?

  15. Compare with filament model > ~  W/W0 = 9.4% at limiter radius, cf. Experiment = 8.8%Excellent agreement given inherent approximations ELMs with <DWELM> 500 kJ deposit ~10% of their energy on the main chamber limiters (for separatrix-wall gap ~ 5 cm) Assume mid-pedestal paramsTe,0 = Ti,0 ~ 800 eVne,0 ~ 3.01019 m-3Dped ~ 4 cmvELM = 600 ms-1 Filament parallel energy loss model(W. Fundamenski, R. A. Pitts, PPCF 48 (2006) 109)

  16. Conclusions • JET can access ELM conditions which match new ITER specifications (DWELM ~ 1 MJ) in pulses at Ip = 3.0 MA with upstream lp ~ 5 mm and divertor wetted area ~1.0 m2 • Strong in-out divertor radiation asymmetry – up to 70% of the ELM energy drop can be radiated, mostly in the divertor volume. • Evidence that thermal decomposition of inner divertor surface layers increases radiation but DTsurf provoked by largest ELMs relatively modest (~ few 100 ºC) • ELM filaments seen clearly at main chamber limiters but only carry ~10% of DWELM for largest ELMs (<DWELM> > 0.5 MJ with fixed wall gap (~5 cm). • Fraction of energy to the wall in good agreement with a parallel loss model of ELM filament evolution

  17. Reserve slides

  18. Large ELMs with low fueling #70226 – no gas fuelling Da (inner) PTOT(MW) WDIA(MJ) Te,ped (keV) ne,ped(1019m-3) H98Y Zeff(Brems) Time (s) Large ELMs have large drop in Te,ped New data populate scaling beyond DTe,ELM/Te,ped = 0.4 R. A. Pitts et al., ITPA, Garching, 2007

  19. Filaments in fast visible light 1 2 3 6 5 4 Courtesy of J. A. Alonso, CIEMAT #70228 Da (inner)divertor WDIA(MJ) Time (s) Frame time 33 ms, main chamber view – filament-wall interaction seen during divertor Da rise. DWELM = 804 kJDERAD = 537 kJ

  20. Parallel ELM transport T, n tELM = 200 ms Post ELM 150 ms dR 2L|| = 80 m Nparticles = 0.8 – 5.0  106, Ncells = 6000High resolution, low noiseTped = 0.5 – 5 keVnped = 0.15 – 15  1019 m-3DWELM = 0.025 – 2.5 MJ Significant progress being made in realistic parallel transport modelling of ELM pulse with the BIT1 PIC code Treat ELM as a square wave pulse launched upstream over time tELM with specified Tped, nped DWELM ~ tELM3npedTped2pLpolRdR Plasma expelled into 1D SOL with cosine distribution centred on midpoint between targets.B = const., inclined targets (~5º) D. Tskhakaya et al., EPS 2007

  21. Test case: PIC vs. expt. PIC ONLY D. Tskhakaya Example: DWELM = 400 kJ Tped = 1.5 keVnped = 51019 m-3 tELM Clear separation of electron (~2 ms) and ion (~100 ms) transit times Assumed “ELM duration” 200 ms te ti

  22. Test case: PIC vs. expt. PIC + EXPT. IR data obtained at outer target (no layers) from coherent average of 20 similar ELMs with <DWELM> ~ 310 ± 66 kJ Time resolution artifically enhanced to 50 ms Good agreement in shape of pulse rise Width a question of time and shape of ELM pedestal loss PIC overestimates expt. by ~ factor 5 D. Tskhakaya, T. Eich, R. A. Pitts Factor ~2 due to known in-out ELM loading asymmetryFactor ~2 due to 1D nature of PICReasonable agreement given how DWELM specified in the code

More Related