Science

A prize-winning story of metallic melting in tokamaks

30 Oct 2023 - Richard Pitts, Experiments & Plasma Operation Section Leader
FEC2023_award_ukaea_2_OK.jpg
Professor Svetlana Ratynskaia of Sweden's Royal Institute of Technology receiving the 2023 Nuclear Fusion prize from IAEA Director General Rafael Grossi during a special awards ceremony held on the first morning of the recent 29th IAEA Fusion Energy Conference. Photo: UKAEA

An obvious property of metals is that they melt, and metallic components directly facing the plasma are one thing ITER has a whole lot of. A total surface area of about 800 m² to be more precise. What happens when they do melt was in the spotlight at the recent 29th IAEA Fusion Energy Conference when a paper focusing on the physics of transient-heat-pulse-induced melting was awarded the prestigious Nuclear Fusion prize for 2023. 

In its 2016 baseline configuration, the ITER machine features beryllium-armoured first wall panels in the "main chamber" region (~650 m²), with another ~140 m² of tungsten targets in the divertor area at the bottom of the vacuum chamber, where more than 90 % of the plasma thermal and particle exhaust will be intercepted. Beryllium has two very important advantages from the point of view of ITER operation. It is an excellent getter for oxygen, reducing the level of impurities in the vacuum vessel, making plasma start-up easier. And even more importantly it is a low atomic number material, so that burning plasmas are much more tolerant to wall atoms eroded by plasma particle impact that make their way into the hot core. Unfortunately, this beneficial property comes with the price of a quite low melting temperature (1287 °C), as well as a low threshold for sputter erosion and a high propensity to trap fusion fuel (deuterium and tritium) in the deposited layers which will grow on plasma-facing surfaces as a result of this erosion. 

In contrast, tungsten is a remarkable material in terms of robustness and has the highest melting point (3422 °C) of all known elements, an elevated energy threshold for sputtering, and negligible fuel retention issues in comparison with beryllium. These are the principal reasons for deploying it in the divertor region, where steady heat and particle fluxes are extremely high. It does, though, have one serious drawback: as a very high atomic number material, only very small concentrations will be tolerable in ITER's fusion plasmas before the burn is compromised. Nevertheless, tungsten is the most widely favoured candidate in most magnetic confinement fusion reactor designs, and the long-term objective of ITER has always been to operate eventually with tungsten as sole plasma-facing material after starting first nuclear operations with beryllium main wall armour to ease the pathway to burning plasmas. In view of technical issues which have arisen during machine construction, this long-term plan is now being brought forward in a new proposal currently under development at the ITER Organization to go "full-tungsten" from the very beginning of operations. A review of this proposal was the sole objective of last month's STAC meeting (ITER Council Science and Technology Advisory Committee) held at ITER Headquarters.

But even the lofty tungsten melting temperature is no match for the reactor-scale plasmas ITER will explore. Thermal energies in ITER fusion plasmas (up to ~350 MJ) and energy stored in the magnetic field (up to ~600 MJ) when operating at high plasma current dwarf those accessible in smaller research tokamaks. Uncontrolled releases of this energy—for example during fast transient magnetohydrodynamic activity and plasma disruptions, or in cases of slow power flux excursions to the divertor targets—can melt even tungsten surfaces. This is a result of the specific character of magnetically confined plasmas: the energy is focused onto very local areas due to the geometry of the magnetic field so that energy densities can be extremely high. 

Pitts_1V.jpg
From the prize-winning paper comparing simulation (using the MEMOS-U code, a much upgraded version of the MEMOS-3D package originally developed at KIT) with two fast-transient-induced tungsten melting experiments on ASDEX Upgrade (upper) and JET (lower). The excellent quantitative reproduction of the key features of the empirical deformation profiles demonstrates that the model adequately accounts for all the key physical processes.
ITER will be equipped with a very sophisticated plasma control system and a barrage of mitigation systems to intercept high-energy transients, but it is not credible to imagine that all possible events can be avoided or mitigated over a decades-long machine lifetime. While not posing a threat to machine integrity, these off-normal heat fluxes can locally modify the surface topology of impacted plasma-facing components and, in the worst cases, require replacement—not an operation to be undertaken lightly or frequently in a device of ITER's scale and complexity.

In fact, tungsten melting in tokamaks is rather different than what would be found in most other situations. At melting temperature, tungsten is a very strong electron emitter (the phenomenon of thermionic emission). When a melt layer forms on a tungsten plasma-facing component, thermionic electrons emitted from the surface towards the plasma leave behind a positive charge, which is neutralized by an electric current (the "replacement current") flowing through the component and the molten layer. These currents can be enormous, easily in the range of millions of amperes per square metre. If they flow across the strong tokamak magnetic fields, Lorentz forces can then drive very rapid melt acceleration (hundreds of metres per second squared). Not only does the surface melt, but the melt layer can be swept over large surface areas. What then does the affected material surface look like topologically, how resilient will it be to subsequent power loading and how much damage can be tolerated before repair is required?

Answering such questions has been a priority at the ITER Organization in recent years. The tokamak community has been long been aware of the issue, but it has lacked a quantitative, experimentally benchmarked framework to make predictions at the reactor scale. Pioneering experiments on deliberate tungsten melting and melt motion under steady heat loading were performed in the TEXTOR tokamak nearly a decade and a half ago, but it was only in 2014 that the first experimental demonstration of "controlled" rapid transient tungsten melting in a tokamak was achieved at JET in a series of tests performed in collaboration with ITER. Numerical simulations were attempted of both the TEXTOR and JET experiments using a melt code (MEMOS-3D) developed at the Karlsruhe Institute of Technology (KIT). Results were reasonable, but something better was required to be confident of extrapolations to ITER. Enter a team led by 2023 Nuclear Fusion prize laureate Prof. Svetlana Ratynskaia from the Space and Plasma Physics Division of the Royal Institute of Technology, Stockholm (KTH). They were approached by the ITER Organization in January 2016 at an International Tokamak Physics Activity (ITPA) Topical Group meeting in Frascati, Italy, and asked to consider joining the effort to improve the framework for melt predictions.  Experts in many detailed theoretical aspects of plasma-material interactions, but not previously involved in tokamak melt simulations, they accepted the challenge and what followed has been an outstanding example of how the fusion research community can mobilize to address key issues for ITER.

Pitts_2V.jpg
From a new paper just submitted to Nuclear Fusion by the KTH group, this is an example of the new MEMENTO melt code applied to the simulation of a dedicated experiment performed at ASDEX Upgrade designed to study the flow of molten tungsten over castellation gaps in plasma-facing components (of which there will be hundreds of thousands in ITER).
Supported in part by the ITER Organization through contracts and co-sponsorship of graduate students, together with funding from EUROfusion and the Swedish Research Council, the KTH team began by refining and adding physics to the KIT simulation tool, guided by results from the older JET exposures and new similar, but better diagnosed, experiments on ASDEX Upgrade. The result was the MEMOS-U code, providing a unified description of metallic plasma-facing component melting in magnetic confinement devices, summarized in the Nuclear Fusion prize-winning paper, just one of the very many journal articles which have emerged from this multiyear effort.

And the work goes on. The physics models describing melting and melt motion have been subject to yet more experimental tests in the past two years on the ASDEX Upgrade, WEST and DIII-D tokamaks, now coordinated through a multi-machine ITPA activity. In parallel, the KTH team has used the improved understanding gained over the years since their first involvement to produce a completely new, modern and flexible modelling tool (the MEMENTO code) which is being applied successfully to the interpretation of the most recent melt experiments and will be applied to others still being planned. This new framework is also now being employed as the main predictive tool for the assessment of gross topological erosion due to transient-induced armour melting on ITER, closing the loop on this truly global effort in support of ITER and providing a key modelling framework for similar assessments at the fusion reactor scale.

Journal article reference: Ratynskaia, S., Thorén, E., Tolias, P., Pitts, R.A., Krieger, K., Vignitchouk, L., Iglesias D., the ASDEX-Upgrade team, the JET contributors, and the EUROfusion MST1 team (2020). Resolidification-controlled melt dynamics under fast transient tokamak plasma loads. Nuclear Fusion, 60. https://iopscience.iop.org/article/10.1088/1741-4326/abadac