Multiphysics model of thermomechanical and helium-induced damage of tungsten during plasma heat transients
Introduction
In nuclear fusion reactors, plasma-facing components (PFCs) (e.g. the divertor and the first wall), are subjected to high fluxes of energetic photons, neutrons, hydrogen, and helium ions. The impingement of these energetic atom fluxes leads to rapid and transient surface heating of the surface. Bombardment by helium isotopes leads to helium-induced damage accompanying micro-structural evolution (e.g. material swelling and the formation of blisters [1], [2], and dislocation loops and helium holes or bubbles [3], [4]). Several recent experiments (e.g., [5]) have shown that damage in the surface region and inside the material may degrade the thermophysical properties of tungsten. Other experiments have indicated that under helium and hydrogen ion bombardment conditions, some near-surface grains are ejected from the bulk to the surface region [6].
Divertors in fusion reactors are subjected to transient plasma events characterized by short durations of high thermal energy. One class of these events is the edge-localized mode (ELM), which is a highly nonlinear event accompanied by high thermal energy (3–10% of the core thermal energy). Typical values fall between 0.1 and 0.5 MJ/m2 for the Joint European Torus (JET) and between 1 and 5 MJ/m2 for the International Thermonuclear Experimental Reactor (ITER). The duration of these events is relatively short, usually between 0.1 and 1 ms.
Recently, tungsten has become a primary candidate material for PFCs because it possesses numerous favorable thermophysical properties [7]. Models that describe and study change and damage in the structure of tungsten are critical in determining the limits for tungsten’s operating conditions in extreme environments. The main objective of this paper is to develop a computational multiphysics model to investigate thermomechanical and helium-induced damage in tungsten under transient energetic ion bombardment conditions. The following sections present in Sections 2 Physical model description, 3 Finite element implementation, 4 Results and discussion, 5 Conclusions.
Section snippets
Thermal conduction
In fusion reactors, PFC materials are exposed to transient high heat fluxes that diffuse inside the material by heat conduction. This is expressed through the following heat equation, which, in the absence of a convective flow, takes the form:where T is the temperature and k is the thermal conductivity. To simulate the transient conditions of the plasma, the tungsten surface is subjected to a heat flux from one side in the form:where n is the normal vector to the
Finite element implementation
The multiphysics model in Section 2 was developed and solved within the finite element analysis framework to study the synergistic effects of transient plasma events and the existence of helium bubbles concentrated along grain boundaries on the thermomechanical damage of tungsten. The modeled geometry was a two dimensional (700 μm × 600 μm) block, as shown in Fig. 3. The method of Voronoi diagrams was used to divide the geometry into a random number of grains with different sizes and orientations.
Results and discussion
During the applied transient heat flux, the temperature rises from its initial value and peaks at the end of the period of application. When the heat flux is turned off, the temperature starts to decrease and diffuses into the bulk of the material until it eventually reaches the temperature of the infinite element and becomes uniform throughout the material. This is illustrated in Fig. 4, Fig. 5 for the two tested heat loads. A comparison of the evolution of the tungsten surface temperature in
Conclusions
High heat fluxes and helium bombardment are characteristics of nuclear reactions in fusion reactors. We investigated the combined effect of severe transient plasma heat loads and helium irradiation on tungsten damage. A multiphysics model that combines thermoelasticity, crack damage, and thermal conduction was developed using the finite element framework. The model examined the relation between temperature magnitudes and distribution, pressure in helium bubbles, and the onset of damage manifest
References (13)
- et al.
Journal of Nuclear Materials
(2007) - et al.
Journal of Nuclear Materials
(2009) - et al.
Journal of Nuclear Materials
(2007) Journal of Nuclear Materials
(2004)- H. Iwakiri, K. Yasunaga, K. Morishita, N. Yoshida. Microstructure evolution in tungsten during low-energy helium ion...
Journal of Nuclear Materials
(2004)