1. Introduction
Two design options are currently under investigation for the electron cyclotron heating and current drive (EC H&CD) system required in the European DEMO to heat up the plasma and to control the magneto-hydrodynamic (MHD) instabilities, especially the neoclassical tearing modes (NTMs), and also the thermal instability of the plasma [
1]. In the first option, the angular steering concept, consisting of antenna launchers with mirrors steered by dedicated actuators, is used to deploy the electron cyclotron (EC) power at different positions in the plasma. On the contrary, the second option features fixed-angle antenna launchers connected to frequency step-tunable gyrotrons [
2]. The frequency steering concept is used for a different radial EC power deposition by tuning the gyrotron frequency in the range of a few seconds. The target is the tuning in steps of 2–3 GHz over a range of ±10 GHz around the main two frequencies to be selected among 136, 170 and 204 GHz [
3].
Broadband diamond window solutions are thus required both on the torus and gyrotron sides of the EC H&CD system. The chemical vapor deposition (CVD) diamond Brewster-angle window is the primary choice, with a minimum diameter of 180 mm of the disk, required for a 63.5 mm aperture allowing 2 MW continuous wave (CW) propagation. Until 2020, several diamond growth experiments were carried out by Diamond Materials [
4] aiming at a target of an optical grade disk with 180 mm diameter and 2 mm thickness [
5,
6]. Excellent results were obtained with reference to the diamond deposition, polishing and quality (represented by loss tangent measurements). In addition, the design of such a window has been advanced with reference to cooling layout and characterized by a consistent and parametric set of analyses in view of the window prototyping activity foreseen in the European DEMO conceptual design phase (2021–2027) [
6].
However, in case of showstoppers, the double-disk diamond window is the backup broadband solution. This window concept was used to design and manufacture a torus window in the context of the EC H&CD system at ASDEX Upgrade for injection up to 1 MW microwaves at four selected frequencies between 105 and140 GHz [
7,
8]. In this paper, first, the possibility of using the double-disk diamond window within the boundary conditions of the EC H&CD system in DEMO was investigated by computational fluid dynamics (CFD) conjugated heat transfer analyses in order to check the fluid behavior and the thermal performance of the window with respect to the reference microwave beam scenario of 2 MW at 204 GHz (worst scenario). Subsequently, the sensitivity of the design with respect to different mass flow rate, loss tangent (tanδ), beam radius and frequency was checked. Last, the resulting temperature distribution was transferred to a FEM structural analysis to validate the design against the allowable limits.
Temperature and thermal stress results showed that the double disk window is a feasible window solution for DEMO, but safety margins against limits shall be increased by introducing design features able to make the fluid more turbulent. A first design change is proposed, showing that, in combination with a higher inlet flow rate, the maximum temperature in the disks can be reduced from 238 to 186 °C, leading, therefore, to lower thermal gradients and stresses in the window.
2. Geometry of the Window
The double disk window designed and manufactured for ASDEX Upgrade is shown in
Figure 1a. It consists of two parallel diamond disks of 106 mm diameter facing each other with a tunable gap of about 5 mm and having the resonant thickness of 1.80 mm for the frequency range limits of 105 and 140 GHz. The window has an aperture of 80 mm. Setting the gap to values in the range of 4.1 to 6 mm, EC power can be transmitted at other two intermediate frequencies with minimum reflected power (around 115 and 125 GHz). Each disk is integrated in a subsystem, forming a housing with the counterpart subsystem. Fine tuning of the disk interspace is performed by dedicated helical springs.
With reference to each subsystem, the disk is brazed to two copper cuffs of 1 mm thickness, allowing cooling only on one side by direct contact with water. These cuffs, named inner and outer cuffs, are brazed to a steel flange by a copper ring and spacer. An intermediate cuff is also used to force the water in the vicinity of the disk. The flange contains inlet and outlet pipes with 6 mm inner diameter and connections to the housing.
Thanks to the symmetry regarding the cooling path in the two subsystems of the window, only half of a subsystem was considered for the DEMO window analyses. The geometry for these analyses is shown in
Figure 1b. The thickness of the disk was increased to 1.85 mm, being in resonant condition with the main frequencies 136, 170 and 204 GHz currently of interest for DEMO. Small geometrical simplifications, not affecting the overall results, were introduced for better quality of the mesh to generate.
Figure 1.
(a) Geometry of the double-disk window designed and manufactured for ASDEX Upgrade; (b) geometry of the window used in the analyses. Reference system, nomenclature of the main parts and materials are shown.
Figure 1.
(a) Geometry of the double-disk window designed and manufactured for ASDEX Upgrade; (b) geometry of the window used in the analyses. Reference system, nomenclature of the main parts and materials are shown.
3. Method: CFD-Conjugated Heat Transfer Analysis
Steady-state CFD conjugated heat transfer analysis was first carried out as a reference case to investigate the cooling and the thermal performance of the window shown in
Figure 1b with respect to the foreseen worst beam scenario, i.e., 2 MW at 204 GHz. The code ANSYS CFX 2021 R1 was used for the analyses. Temperature-dependent properties for pure copper, CVD diamond and steel were taken from [
9,
10,
11], while the properties directly from CFX’s library were used for the water coolant. A fine mesh was generated for the window geometry with 13.3 × 10
6 elements. In particular, a mesh size of 0.5 mm was applied to the disk and cuffs, while, in the fluid model, a size of 0.4 mm was adopted and a very fine mesh (inflation layer with first element size of 10 µm) was generated at the boundary layer to model properly the near wall interactions (heat transfer). For the latter reason, in the analysis settings, the k-omega shear stress transport (SST) model was selected as the turbulence model.
Symmetry was applied to the xy plane (
Figure 1b). A mass flow rate of 0.167 kg s
−1 (10 L/min) was assumed for the inlet (0.0833 kg s
−1 in the symmetric model), while a reference pressure of 0 Pa was applied to the outlet. The inlet temperature of the water was set to 20 °C. A Gaussian mm-wave beam was considered as the analysis refers to a gyrotron output window. First, the absorbed power P
abs in one disk was calculated by [
12]:
where P
abs is the absorbed power, P
beam = 2 MW is the beam power, f = 204 GHz is the beam frequency, t = 1.85 mm is the disk thickness, tanδ = 3.5 × 10
−5 is the reasonable assumed loss tangent, ε
r = 5.67 is the dielectric constant of diamond and c
0 is the speed of light in vacuum. It turned out an absorbed power of 1847 W. As a benchmark, it is interesting to observe that, in the Brewster-angle disk configuration, with the same geometrical features and properties, using the appropriate formula also provided in [
12], the absorbed power in the disk only would have been 1430 W. It means that, due to the resonance thickness, there is higher absorption in the window planar configuration of about 1.3 times. The heat load was then applied to the disk in terms of volumetric power density along the radial coordinate q
(r) in W m
−3 by:
where P
abs = 1847 W is the absorbed power in the disk, r is the radial coordinate and w
0 = 20 mm is the assumed beam radius. The Gaussian distribution was normalized in order to obtain in the disk the calculated absorbed power of 1847 W.
The average y+ (the non-dimensional distance) at the cooling interface lower than 1 and average residuals of the solved equations lower than default target of 1 × 10−4, used in many engineering applications, provided confidence in the accuracy of the performed analysis. In addition, an analysis run with a much finer mesh than the reference case showed the grid independency of the results. In fact, with a mesh of 25 × 106 elements, i.e., having almost a double number of elements, no changes occurred in the results.
Then, additional CFD conjugated heat transfer analyses were carried out to investigate the sensitivity of this double-disk window design with respect to different mass flow rate, loss tangent, beam radius and frequency.
Last, a first conceptual design change was introduced to increase the fluid turbulence along the cooling path in order to improve the heat exchange and obtain lower temperatures in the window. The design change was checked by CFD conjugated heat transfer analysis to evaluate the impact on the safety margin with respect to the reference case. As shown in
Figure 2, the intermediate cuff was simply modified by increasing its thickness between the inner and outer cuffs and introducing a series of holes having a diameter of 2 mm. In this way, for the same flow rate, the velocity of the fluid along the cooling path is increased with the aim to make the heat removal more effective.
4. Method: Structural Analysis
A structural analysis was performed by the code ANSYS Workbench 2021 R1 to check, in the reference case, the thermal stresses generated in the window by the beam power absorption. The residual stresses due to the brazing process are not considered. The geometry in
Figure 1b (without the coolant) was used with the boundary conditions shown in
Figure 3a. Again, symmetry was applied in the xy plane, and the structure was made isostatic by applying a condition of zero displacement to all three coordinates at point B and zero displacement only to the x coordinate at the point A. The temperature distribution was applied as load and the stresses were taken in terms of first principal stresses for diamond and equivalent von Mises stresses for the metallic parts.
An accurate mesh was generated, especially at the critical regions where high stresses are expected, i.e., at the interface copper–diamond (
Figure 3b). The inner and outer cuffs at the interface with the disk experience plastic deformation as the stresses become close to the yield strength of copper for the temperature range of interest (minimum yield strength of copper for plate products is 57 MPa at 80 °C [
9]). A plastic steady-state structural analysis was, therefore, carried out with the temperature-dependent properties for pure copper, CVD diamond and steel taken from [
9,
10,
11]. The multilinear isotropic hardening was used as plasticity material model for copper with stress–strain curves given for plate products at different temperatures (up to 250 °C).
Figure 3.
(a) Boundary conditions applied in the structural analysis for the reference case; (b) details of the mesh generated at the disk interface with inner and outer cuffs.
Figure 3.
(a) Boundary conditions applied in the structural analysis for the reference case; (b) details of the mesh generated at the disk interface with inner and outer cuffs.
6. Conclusions
The double-disk CVD diamond window is the broadband backup window solution for DEMO. It was characterized by CFD-conjugated heat transfer and structural analyses for the 2 MW at 204 GHz worst beam scenario, also considering the impact of changing inlet mass flow rate, the frequency and radius of the beam and the tanδ of diamond. The temperature and consequent stress results showed that it is a feasible window solution for DEMO, but safety margins against limits shall be increased by introducing features aiming to make the fluid more turbulent. This would allow counteracting factors such as potential higher tanδ degradation during disk brazing.
A first conceptual design change was introduced and showed that, e.g., in combination with a higher flow rate, a maximum temperature of 186 °C might be achieved in the diamond disk, leading to lower thermal gradients and thus stresses in the window. In the near future, the double disk window shall be characterized for DEMO from the perspective of radio frequency beam transmission. This is important, for instance, to determine, depending on the possible distance to set between the two parallel diamond disks, the intermediate frequencies at which minimum reflection conditions are met (beyond the frequencies that already fulfil the resonant condition with the thickness of the two disks).