Divertor Heat Flux Research Articles

EAST steady-state long pulse H-mode with core-edge integration for CFETR

A recent EAST experiment has successfully demonstrated long pulse steady-state high plasma performance scenario and core-edge integration since the last IAEA in 2018. A discharge with a duration over 60 s with β P ∼ 2.0, β N ∼ 1.6, H98y2 ∼ 1.3 and an internal transport barrier on the electron temperature channel is obtained with multi-RF power heating and current drive. A higher β N (β N ∼ 1.8, β p ∼ 2.0, H98y2 ∼ 1.3, n e/n GW ∼ 0.75) with a duration of 20 s is achieved by using the modulated neutral beam and multi-RF power, where several normalized parameters are close or even higher than the phase III 1 GW scenario of CFETR steady-state. High-Z impurity accumulation in the plasma core is well controlled in a low level by using the on-axis ECH. Modeling shows that the strong diffusion of TEM turbulence in the central region prevents tungsten impurity from accumulating. More recently, EAST has demonstrated compatible core-edge integration discharges in the high β p scenario: high confinement H98y2 > 1.2 with high β P ∼ 2.5/β N ∼ 2.0 and f bs ∼ 50% is sustained with reduced divertor heat flux at high density n e/n GW ∼ 0.7 and moderate q 95 ∼ 6.7. By combining active impurity seeding through radiative divertor feedback control and strike point splitting induced by resonant perturbation coil, the peak heat flux is reduced by 20–30% on the ITER-like tungsten divertor, here a mixture of 50% neon and 50% D2 is applied.

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of ‘kink’ phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor ‘wetted’ area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

Development and implementation of Divertor Fast Particles Injection for EAST tokamak

The divertor is the main device in the tokamak that is subjected to heat loads, excessive high heat flux could cause damage to the divertor target. Divertor heat flux control is very important for magnetically confined fusion devices. As a method for divertor heat flux control, the Divertor Fast Particles Injection (DFPI) is successfully developed and implemented on the Experimental Advanced Superconducting Tokamak (EAST). The components of DFPI are a high-pressure gas source (0.6–1 MPa), a pulse solenoid valve, a nozzle and the pipeline etc. Compared with the Supersonic Molecular Beam Injection (SMBI) particle injection from midplane, the particles injected by DFPI are more difficult to enter the plasma core and have less impact on the performance of the core plasma; Compared with the gas puffing (GP) particle injection from the divertor, the particles injected by DFPI can enter the divertor area faster, the delay time of DFPI is about 25 ms. The successful implementation of DFPI provides an important tool for divertor heat flux control on EAST.

Prediction of divertor heat flux width for ITER pre-fusion power operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second pre-fusion power operation (PFPO) phases specified in the new ITER research plan using BOUT++ transport code (Li et al 2018 Comput. Phys. Commun. 228 69–82). The initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer. An anomalous thermal diffusivity scan is performed with E × B and magnetic drifts. The results in two scenarios identifying two distinct regimes: a drift-dominant regime when diffusivity is smaller than the respective critical thermal diffusivity χ c and a turbulence-dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λ q . From transport simulations, we obtain the critical thermal diffusivity χ c = 0.5 m2 s−1 for the PFPO-1 scenario with toroidal magnetic field B = 1.77 T and plasma current I p = 5 MA, and χ c = 0.3 m2 s−1 for the PFPO-2 scenario with toroidal magnetic field B = 2.65 T and plasma current I p = 7.5 MA. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift-dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity using ITER scenarios with fixed safety factor q 95, major radius R, aspect ratio R/a, and the separatrix temperature T sep, and establishes the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q 95, R, R/a, and T sep fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χ c, possibly yielding a transition from turbulence- to drift-dominant regime.

Operating a full tungsten actively cooled tokamak: overview of WEST first phase of operation

WEST is an MA class superconducting, actively cooled, full tungsten (W) tokamak, designed to operate in long pulses up to 1000 s. In support of ITER operation and DEMO conceptual activities, key missions of WEST are: (i) qualification of high heat flux plasma-facing components in integrating both technological and physics aspects in relevant heat and particle exhaust conditions, particularly for the tungsten monoblocks foreseen in ITER divertor; (ii) integrated steady-state operation at high confinement, with a focus on power exhaust issues. During the phase 1 of operation (2017–2020), a set of actively cooled ITER-grade plasma facing unit prototypes was integrated into the inertially cooled W coated startup lower divertor. Up to 8.8 MW of RF power has been coupled to the plasma and divertor heat flux of up to 6 MW m−2 were reached. Long pulse operation was started, using the upper actively cooled divertor, with a discharge of about 1 min achieved. This paper gives an overview of the results achieved in phase 1. Perspectives for phase 2, operating with the full capability of the device with the complete ITER-grade actively cooled lower divertor, are also described.

Advances in the long-pulse steady-state high beta H-mode scenario with active controls of divertor heat and particle fluxes in EAST

Since the last IAEA-Fusion Energy Conference, the Experimental Advanced Superconducting Tokamak (EAST) research program has been, in support of ITER and CFETR, focused on development in terms of the long-pulse steady-state (fully noninductive) high beta H-mode scenario with active controls of the stationary and transient divertor heat and particle fluxes. The operational domain of the steady-state H-mode plasma scenario has been significantly extended with ITER-like tungsten mono-block divertor, plasma control and heating schemes. EAST has achieved several important milestones in the development of high β p H-mode scenario and its key physics and technologies. A 60 s-scale long-pulse steady-state high β p H-mode discharge with the major normalized plasma parameters similar to the designed performance of the CFETR 1 GW fusion power operation scenario has been successfully established and sustained by pure RF heating and current drive. Several feedback control schemes have been developed for a sustained detachment with good core confinement. This includes control of the total radiation power, target electron temperature, and particle flux measured using divertor Langmuir probes or a combination of the control of target electron temperature and AXUV radiation near the X point. The detachment feedback control schemes have been integrated with small-ELM regimes and high β p scenario via neon seeding, enabling a core and edge compatible integrated high-beta scenario applicable to long-pulse operations. ELM suppression has been achieved using various methods, including resonant magnetic perturbations and impurity seeding. Full suppression of ELMs by using n = 4 RMPs has been demonstrated for ITER for the first time in low input torque plasmas in EAST. EAST has been operated with helium to support the ITER research requirements for the first time. For a long-pulse, high bootstrap current fraction operation, a new lower tungsten divertor with active water-cooling has been installed, along with improvements in the heating and current drive capability.

Jitter correction for EAST infrared videos using edge-based phase correlation method

EAST has currently equipped with four wide-angle infrared (IR)/visible integrated endoscope diagnostics for routinely monitoring the status of plasma facing components (PFCs) and providing temperature information to ensure safe operation of the device. However, due to the existence of various sources of vibration force during the plasma experiments, image jitter can be clearly observed from the recorded IR videos, which increases the errors of temperature extraction from region of interest (ROI) within IR images and hence reduces the accuracy of heat flux calculation. In order to solve this problem, a method based on edge detection and phase correlation has been recently developed for correcting jitter of the IR videos on EAST, which has the advantages of high robustness, none accumulated error and high efficiency. Using this method, jitter of the IR video is successfully detected and corrected. Comparison between divertor heat flux results calculated from IR videos with and without jitter correction proves that the application of the edge-based phase correlation method for jitter correction can effectively improve the accuracy of the heat flux calculation.

Investigating the n = 1 and n = 2 error fields in W7-X using the newly accelerated FIELDLINES code

No fusion device can be created without some uncertainty; there is always a slight deviation from the geometric specification. These deviations can add up create a deviation of the magnetic field. This deviation is known as the (magnetic) error field. Correcting these error fields is desired as they cause asymmetries in the divertor loads and can thus cause damage to the device if they grow too large. These error fields can be defined by their toroidal (n) and poloidal number (m). The correction of the n = 1 and n = 2 fields in Wendelstein 7-X (W7-X) is investigated in this work. This investigation focuses on field line diffusion to the divertor, a proxy for divertor heat flux. Such work leverages the 25× speedup obtained through the implementation of a new particle-wall collision model. The n = 1 and n = 2 error fields of the as-built coils model of W7-X are corrected by scanning phase and amplitude of the trim and control coils. Reductions in the divertor load asymmetry by factors of four are demonstrated using error field correction. It is found that the as-built coils model has a significantly lower error field than found in experiments (Bozhenkov et al 2018 Nucl. Fusion 59 026004).

Impact of edge harmonic oscillations on the divertor heat flux in NSTX

Previously, edge-localized mode (ELM)-free and inter-ELM divertor peak heat flux reduction induced by an edge harmonic oscillation (EHO) was observed in NSTX [Gan et al., Nucl. Fusion 57, 126053 (2017)]. This paper introduces new analysis of the EHO impact on the divertor heat flux. It was found that enhanced edge turbulence significantly altered the divertor heat flux footprint in NSTX H-modes, relative to ELM-free discharges. When the background edge fluctuation level was low, the EHO significantly increased the heat flux width and decreased the divertor peak heat flux. When the background edge fluctuation level was high, the EHO actually increased the divertor peak heat flux. Furthermore, it was also found that the heat flux width increased with the increasing frequency of the EHO.

Fluid turbulence simulations of divertor heat load for ITER hybrid scenario using BOUT++

The BOUT++ six-field turbulence code is used to simulate the ITER 11.5 MA hybrid scenario and a brief comparison is made among ITER baseline, hybrid and steady-state operation (SSO) scenarios. Peeling–ballooning instabilities with different toroidal mode numbers dominate in different scenarios and consequently yield different types of ELMs. The energy loss fractions (ΔW ped/W ped) caused by unmitigated ELMs in the baseline and hybrid scenarios are large (∼2%) while the one in the SSO scenario is dramatically smaller (<1%), which are consistent with the features of type-I ELMs and grassy ELMs respectively. The intra ELM divertor heat flux width in the three scenarios given by the simulations is larger than the estimations for inter-ELM phase based on Goldston’s heuristic drift model. The toroidal gap edge melting limit of tungsten monoblocks of divertor targets imposes constraints on ELM energy loss, giving that the ELM energy loss fraction should be smaller than 0.4%, 1.0%, and 1.2% for ITER baseline, hybrid and SSO scenarios, correspondingly. The simulation shows that only the SSO scenario with grassy ELMs may satisfy the constraint.

Physics design of new lower tungsten divertor for long-pulse high-power operations in EAST

A new lower tungsten divertor has been developed and installed in the EAST superconducting tokamak to replace the previous graphite divertor with power handling capability increasing from <2 MW m−2 to ∼10 MW m−2, aiming at achieving long-pulse H-mode operations in a full metal wall environment with the steady-state divertor heat flux of ∼10 MW m−2. A new divertor concept, ‘corner slot’ (CS) divertor, has been employed. By using the ‘corner effect’, a strongly dissipative divertor with the local buildup of high neutral pressure near the corner can be achieved, so that stable detachment can be maintained across the entire outer target plate with a relatively lower impurity seeding rate, at a separatrix density compatible with advanced steady-state core scenarios. These are essential for achieving efficient current drive with low-hybrid waves, a low core impurity concentration and thus a low loop voltage for fully non-inductive long-pulse operations. Compared with the highly closed small-angle-slot divertor in DIII-D, the new divertor in EAST exhibits the following merits: (1) a much simpler geometry with integral cassette body structure, combining vertical and horizontal target plates, which are more suitable for actively water-cooled W/Cu plasma facing components, facilitating installation precision control for minimizing surface misalignment, achieving high engineering reliability and lowering the capital cost as well; (2) it has much greater flexibility in magnetic configurations, allowing for the position of the outer strike point on either vertical or horizontal target plates to accommodate a relatively wide triangularity range, δ l = 0.4–0.6, thus enabling to explore various advanced scenarios. A water-cooled copper in-vessel coil has been installed under the dome. Five supersonic molecular beam injection systems have been mounted in the divertor to achieve faster and more precise feedback control of the gas injection rate. Furthermore, this new divertor allows for double null divertor operation and slowly sweeping the outer strike point across the horizontal and vertical target plates to spread the heat flux for long-pulse operations. Preliminary experimental results demonstrate the ‘corner effect’ and are in good agreement with simulations using SOLPS-ITER code including drifts. The EAST new divertor provides a test-bed for the closed divertor concept to achieve steady-state detachment operation at high power. Next step, a more closed divertor, ‘sharp-cornered slot’ divertor, building upon the current CS divertor concept, has been proposed as a candidate for the EAST upper divertor upgrade.

Analyses of monoblocks divertor at one-sided high heat flux using RPI wall boiling model

The divertor is one of the core components in the Tokamak magnetic confinement nuclear fusion device, which is responsible for eliminating high-flux energy flow and particle flow. The divertor in China Fusion Engineering Test Reactor (CFETR) is designed to withstand a steady-state thermal load of 10 MW/m2 and a transient thermal load up to 20 MW/m2. The water-cooled scheme is a promising candidate scheme. The model under high heat flux divertor tube system was established by using RPI wall boiling model and the Eulerian multiphase model. On this basis, the fluid-thermal-structure coupling method is used to simulate the complete structure of divertor blocks. In the steady-state condition, three heat transfer patterns are found and compared in the circumferential and axial directions, and methods to enhance heat transfer (including flow velocity and tube structure) are further evaluated. The results reveal the superiority of the screw tube and the swirl tube in the enhancement of heat transfer with high heat flow. For unsteady conditions, the response of temperature and stress for different components are different. The change of flow boundary quickly affects the whole system, while the thermal shock only poses a great threat to external components, and its influence on internal components is lagging and slight.

I-mode operation in helium plasma with pure radio frequency wave heating and ITER-like tungsten divertor on EAST

The characteristics of I-mode operation in helium (He) plasma with radio-frequency (RF) wave heating and ITER-like tungsten divertor have been investigated for the first time on EAST, under the favorable B T configuration. Stationary I-mode in He plasma lasting for ∼5.4 s was achieved at I p = 0.5 MA, B T = 2.4 T with a total injected RF power of 3.5 MW, along with the appearance of the 30–90 kHz weakly coherent mode and the ∼13 kHz edge temperature ring oscillation. The pulse length was only limited by the heating power duration. He concentration was found to play a critical role in I-mode operation, as L–I transition threshold power was increased and the energy confinement was degraded with increasing He concentration. The experimental investigation result on energy confinement and divertor heat flux between I- and H-mode plasma was also given at a similar He concentration.

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

Simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

Comparing spontaneous and pellet-triggered ELMs via non-linear extended MHD simulations

Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. This paper presents non-linear simulations of spontaneous type-I ELMs and pellet-triggered ELMs in ASDEX Upgrade performed with the extended MHD code JOREK. A thorough comparison of the non-linear dynamics of these events is provided. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in A Cathey et al (2020 Nuclear Fusion 60 124007). Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to accurately capture the balance between stabilising and destabilising terms for the edge instabilities. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n = 1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and a narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluence is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.

Effects of a coherent mode on the redistribution of divertor particle flux in HL-2A H-mode plasmas

The effects of a coherent mode (CM) at a frequency of f = 13–16 kHz on the redistribution of divertor particle flux have been studied in HL-2A H-mode plasmas using multiple diagnostics. The experimental results demonstrate that a CM localized at the middle of the pedestal results in outward particle transport, limiting the growth of the pedestal density gradient in the recovery phase of the pedestal, leading to a substantial increase in the localized particle flux striking the outer divertor target, which clearly deviates from the primary peak at the strike point. These observations indicate that the CM might play a crucial role in redistributing the divertor particle flux, which is important for understanding and controlling the divertor heat flux in future fusion-relevant devices.

Theory-based scaling laws of near and far scrape-off layer widths in single-null L-mode discharges

Theory-based scaling laws of the near and far scrape-off layer (SOL) widths are analytically derived for L-mode diverted tokamak discharges by using a two-fluid model. The near SOL pressure and density decay lengths are obtained by leveraging a balance among the power source, perpendicular turbulent transport across the separatrix, and parallel losses at the vessel wall, while the far SOL pressure and density decay lengths are derived by using a model of intermittent transport mediated by filaments. The analytical estimates of the pressure decay length in the near SOL is then compared to the results of three-dimensional, flux-driven, global, two-fluid turbulence simulations of L-mode diverted tokamak plasmas, and validated against experimental measurements taken from an experimental multi-machine database of divertor heat flux profiles, showing in both cases a very good agreement. Analogously, the theoretical scaling law for the pressure decay length in the far SOL is compared to simulation results and to experimental measurements in TCV L-mode discharges, pointing out the need of a large multi-machine database for the far SOL decay lengths.

Real-time feedback system for divertor heat flux control at COMPASS tokamak

The reduction of the incident heat flux onto the divertor will be a necessity for the future thermonuclear reactors. Impurity seeding is recognized as an efficient way to achieve the partial detachment regime, which allows to dissipate a large fraction of power flux by radiation. This paper presents a heat flux real-time feedback system (RTFS) based on impurity seeding controlled by a combined ball-pen and Langmuir probe divertor array in the COMPASS tokamak. A number of features of the system have been studied, such as the type of impurity, seeding location, constants used in the real-time controller and the diagnostic selections. A detailed description of the designed RTFS and the results of the implementation are presented. The findings confirm the applicability of the RTFS for reduction and control of the divertor heat fluxes. Another important implication of this research is the ability of installing such systems in next-step devices.

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink. A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m2. It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

Progress of Divertor Heat and Particle Flux Control in EAST for Advanced Steady-State Operation in the Last 10 Years

Active control of the excessively high heat and particle fluxes on the divertor target plates is of fundamental importance to the steady state operation of tokamaks, especially for fusion reactors. A series of experiments have been carried out on this critical issue to relieve the plasma-wall interactions in the experimental advanced superconducting tokamak (EAST) in the last ten years, not only contributing to the long pulse operation of EAST itself, but also providing physical understandings and potential techniques to the next-generation devices like ITER. We have characterized the power deposition pattern and broadened the divertor footprint width effectively. The plasma-wetted area is actively handled using either 3-dimentional edge magnetic topology or advanced plasma equilibrium, thereby peak heat flux around the strike point is reduced. Active control of detachment or radiation compatible with core plasma performance has progressed significantly in very recent years, with a series of active feedback control modules developed and utilized successfully, based on the divertor physics advances with both experiments and simulation. The upper divertor of EAST was upgraded from graphite to active water-cooling ITER-like tungsten in 2014, exhibiting much enhanced heat removal capability. As for the particle exhaust including both fueling and impurity particles, in addition to wall conditioning and impurity source control, the efficiency of particle flux exhaust is optimized by making full use of the divertor closure and the plasma drifts in both scrape-off layer and divertor volume. These heat and particle exhaust advances contribute greatly to a series of EAST achievements like H-mode operation over 100 s. A brief near-term plan on the integrated control of divertor plasma-wall interactions in long-time scale will also be introduced, aiming to provide favorable divertor operation solution for ITER and CFETR.