Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test reactor”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integrated simulation of U-10Mo monolithic fuel swelling behavior

Here, a separate computational branch has been implemented within the DART (Dispersion Analysis Research Tool) computational code to simulate the swelling behavior of U-10Mo monolithic fuel under the operating conditions of high-power research and test reactors (RTRs). The monolithic branch of the DART code implements a mechanistic rate-theory-based fission-gas-behavior model for the calculation of fission gas swelling, as well as a suite of thermal, physical, and mechanical models to take into account various processes occurring in RTR fuels during irradiation. In order to accurately simulate and eventually predict U-10Mo monolithic fuel irradiation behavior, the code uses materials properties calculated with lower length-scale computational methods, such as gas atom diffusivity and U-Mo surface energy from atomic simulations and grain-morphology-specific recrystallization kinetics (recrystallized fuel volume fractions vs. fission density) predicted using the phase-field method. The remainder of fission gas behavior parameters used in the model were calibrated with measured intergranular bubble size distributions. With this integrated simulation approach, the swelling behavior of U-10Mo monolithic fuel was simulated for various initial grain sizes at different operating conditions and compared with measured data. Furthermore, because limited experimental data exist for parameter calibration detailed sensitivity studies for the important parameters used in the fission gas behavior model were performed in order to examine their impact on both intergranular gas bubble morphology at low fission density, and on total porosity at high fission density.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Synchrotron micro-computed tomography analysis of neutron-irradiated U-Mo fuel

The three-dimensional (3D) microstructure of neutron-irradiated uranium-10 wt.% molybdenum (U-10Mo) fuel with a burn-up of 9.8 × 10 21 fissions/cm 3 was characterized using a novel, multi-modal synchrotron micro-computed tomography approach combining propagation-based phase-contrast enhanced and absorption contrast techniques. The porosity development, porosity interconnectedness, swelling, composition, local thickness of the zirconium (Zr) diffusion barrier, and the influence of the fuel–cladding interaction on the local composition and pore morphology, were uniquely determined in 3D. Two cuboids were produced using a focused ion beam-scanning electron microscope at the Zr diffusion barrier–fuel interface and in the bulk fuel. The bulk fuel sample swelled by 53.3 [+9.7/−3.1]%, while the fuel near the Zr–fuel interface swelled by 63.3 [+14.7/−7.2]%. The average local thickness of the Zr diffusion barrier decreased by 53 %, compared to the expected pre-irradiated thickness. Four pore morphology regions were identified initiating parallel to the fuel–Zr interaction region: (1) an interaction layer of suppressed porosity, (2) a layer of elongated and interconnected porosity, (3) a transition zone of low porosity, and (4) a layer of unoriented porosity representative of the bulk fuel behavior. The increase in porosity near the diffusion barrier corresponded to a higher U concentration compared to that in the bulk fuel. The interconnected porosity in the fuel near the diffusion barrier was extensive and oriented parallel to the diffusion barrier, while the bulk fuel had more compact and isolated pore networks. The interaction layer, despite having suppressed porosity, was nearly 100 wt.% U. Porosity suppression at the diffusion barrier corresponds to the expected reduction in radiation-driven diffusion of Xe at the interface despite the anticipated increase in fission product nucleation originating from a higher U concentration. In conclusion, the novel 3D insights of the porosity, swelling, and compositional variations characterized herein can improve the fidelity of fuel performance codes for proliferation-resistant fuels for research and test reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

A critical review of the history of fabricating monolithic U-Mo fuel plates

The fabrication of monolithic uranium-molybdenum alloy fuels, specifically those developed for high-performance research and test reactors, began in the early 2000s. The primary fuel form consists of uranium alloyed with 10-wt% molybdenum in a thin foil coated with zirconium and encapsulated in aluminum-6061 cladding. Over the years, the process has evolved with different types of casting, heat treatments, rolling schedules, and cladding applications. This review examines the history of these fabrication processes and their impact on microstructure and fuel-swelling performance. Even though various fabrication methods were used, we found little correlation between fabrication variation and fuel swelling. This insensitive relationship between fabrication variation and fuel swelling is primarily due to inhomogeneous microstructures that formed during casting and grain refinement that occurred during rolling. We conclude that the fabrication processes we examined produced similar microstructures, indicating that the fuel microstructure is somewhat insensitive to the fabrication parameters evaluated. However, the relatively small amount of historic data, such as those for grain sizes, limited this analysis. More recently fabricated materials, such as those from ongoing irradiation experiment, Mini-Plate-1 and Mini-Plate-2, were also excluded from this analysis and are intended to be reviewed separately. The findings, that fuel microstructure is somewhat insensitive to the fabrication parameters, do not imply that any fabrication method is acceptable, given the uncertainties in data and fuel-swelling observations. For example, only arc melting and vacuum induction melting casting processes were previously explored in the historic fabrication efforts. In conclusion, the findings should not be extrapolated to other casting processes.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Micro X-ray computed tomography examination of mini plate fuel with hot isostatic pressed aluminum cladding

In order to minimize proliferation risks and improve security of nuclear material, the United States high performance research and test reactors (USHPRR) program is tasked with converting nuclear reactors that are fueled with highly enriched uranium (HEU) fuels to operate with low-enriched uranium (LEU) nuclear fuels. One favorable LEU fuel configuration is plate fuel with a metallic uranium-molybdenum foil clad within an aluminum alloy (AA 6061). In this fuel, the aluminum cladding is bonded with a hot isostatic pressing (HIP) method to seal the cladding around the fuel meat. However, the HIP process parameters influence the cladding performance, as a defective or incomplete bond can cause a pathway for corrosion. Micro X-ray computed tomography (XCT), a nondestructive technique that provides volumetric imaging, can be applied to inspect fuel plate cladding at the engineering scale. In this work, XCT methodology was developed and successfully utilized to not only observe the bond line of unirradiated mini fuel plates, but to also identify subsurface abnormalities in the plates’ cladding. Importantly, in future work, this technique could be applied to fuel plates, pre- and post-irradiation, to quantify irradiation effects on cladding defects and bond line integrity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling and simulation study for the design of the Fuel Interrogation and Examination using Submersible Tomography Analysis Mk II instrument

Here, the design of a submersible, gamma-ray tomography system for imaging irradiated nuclear fuel is described. The system—named Fuel Interrogation and Examination using Submersible Tomography Analysis (FIESTA) Mk. II—is a variation on a previous Mk. I I design, which was developed to non-destructively image fuel capsules irradiated in the Advanced Test Reactor at the Idaho National Laboratory. The FIESTA system uses a combination of transmission computed tomography and emission computed tomography to image the restructuring and fission product migration at different points of burnup. Changes made to FIESTA Mk. I reflect the revised design requirements and a need to reduce background noise, largely originating from downscattered photons from fuel and transmission source. The computational design and radiation transport simulations are described.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗
DOE OSTI2022

Sensitivity study of coupled chemical-CFD simulations for analyzing aluminum-clad spent nuclear fuel storage in sealed canisters

We report the United States Department of Energy (DOE) manages over 50 Metric Tons Heavy Metal (MTHM) of aluminum-clad spent nuclear fuel. One main source for DOE’s Aluminum-clad spent nuclear fuel (ASNF) inventory is the advanced test reactor (ATR) at the INL site, which makes this fuel of particular interest for storage scenarios. Road-ready and final disposition packaging configurations for the ATR fuel dictates storage within helium-backfilled, sealed DOE standard canisters. The conditions within these sealed canisters for extended (greater than50 year dry) storage is of interest. To further this goal, a three-dimensional (3D) multi-physics computational fluid dynamics (CFD) model is developed of the sealed DOE standard canisters. This 3D CFD model is one-way coupled with bulk gas radiolysis reactions considering sealed canisters with inert gas and possible trace amounts of air and water vapor. This study looks at the evolution of the thermal history of the canisters over a 50 year time period with a coupling to the chemical reactions occurring from radiolytic breakdown of residual water. A sensitivity study is then carried out over the parameters of the model including the fuel decay heat, residual water content, sealed pressure, canister external temperature, and canister emissivity. In pure helium, hydrogen generation rates are low, under 10 ppm, but hydrogen generation rates are affected greatly by the presence of even 1% residual air, increasing by 50-plus-fold, and nitric acid generation with residual air also occurs ranging from 500 to 4000 ppm after 50 years. The fuel decay heat and the residual water content show the most importance in the generation of hydrogen gas in pure air, and for nitric acid with a residual air condition. External temperature, canister emissivity and sealed pressure all show minor sensitivity effects to the generation of potentially harmful species.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Spark plasma sintering of fuel meats for U 3 O 8 based dispersion fuels

Research and test reactors often use dispersion-type fuel due to its increased thermal conductivity and burn-up capabilities compared to conventional fuel. Al-U 3 O 8 (aluminumtriuranium octaoxide) dispersion fuels have several advantages over their competitors, such as higher service temperature and better stability of oxygen stoichiometry. However, the two-step fabrication of dispersion fuel causes undesirable porosity in cold-pressed fuel meats that is preserved in co-extruded fuel plates. To combat this, spark plasma sintering (SPS) was used for the fabrication of Al-15, 20, and 30 vol% U 3 O 8 and 8 and 12 vol% Mo-U 3 O 8 fuel meats for the Al-U 3 O 8 time. The in situ SPS data was used to construct and validate Master Sintering Curves (MSCs) with accuracies in Al fuels at 0.02 g/cm 3 , and Mo fuels at 0.07 and 0.17 g/cm 3 . The as-sintered fuel meats were characterised using x-ray diffraction (XRD) and scanning electron microscopy (SEM) to understand chemical and physical changes following the SPS process. The pellets exhibited very high relative densities, the U 3 O 8 was observed to undergo reduction to UO 2 .

Aluminium↗
DOE OSTI2023

Mechanistic Insights into Nonoxidative Ethanol Dehydrogenation on NiCu Single-Atom Alloys

Ethanol dehydrogenation presents a promising pathway towards the production of acetaldehyde, a valuable building block in chemicals production. Under non-oxidative conditions, the reaction is facilitated by supported Cu nanoparticles which afford reasonable activity and high selectivity. The stability issues associated with Cu nanoparticle sintering can be addressed by the addition of small amounts of Ni, which further boost reactivity while retaining selectivity. Despite the promise of NiCu single-atom alloys for non-oxidative ethanol dehydrogenation, little is known about the role of each component and the pathway of this mechanistically complex process. Herein, kinetic investigations from reactor tests identify C-H bond scission as the rate limiting step, while 1-hydroxyethyl is detected as the intermediate via IR spectroscopy. Temperature program desorption studies are employed to examine the effect of Ni coverage and to demonstrate that Ni atoms activate ethanol selectively at lower temperatures, resulting in higher acetaldehyde yield than pure Cu. Temperature program desorption experiments also reveal the spillover of intermediates from the Ni atom to neighboring Cu sites as a relevant step in the reaction pathway. Density functional theory calculations are used to investigate the 2 reaction energetics and to confirm that C-H bond scission is the initial reaction step, while a clear effect of H 2 partial pressure on the reaction pathway is realized. Further, counter to the expected behavior that all reaction steps take place on the Ni atoms, our degree of rate control analysis reveals that a mechanism involving spillover of the 1-hydroxyethyl intermediate from the Ni atom to the Cu surface, where it will dehydrogenate further, is more likely. Furthermore, our combined kinetic, spectroscopic, and theoretical approach sheds light on this complex reaction mechanism and represents a promising method for the understanding and designing of highly active, selective, and stable single-atom alloys for other multistep catalytic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Selective Methane Oxidation to Methanol on ZnO/Cu 2 O/Cu(111) Catalysts: Multiple Site-Dependent Behaviors

Due to the abundance of natural gas in our planet, a major goal is to achieve a direct methane to methanol conversion at medium to low temperatures using mixtures of methane and oxygen. Here, we report an efficient catalyst, ZnO/Cu 2 O/Cu(111), for this process investigated using a combination of reactor testing, scanning tunneling microscopy, ambient-pressure x-ray photoemission spectroscopy, density functional calculations, and kinetic Monte Carlo simulations. The catalyst is capable of methane activation at room temperature and transforms mixtures of methane and oxygen to methanol at 450 K with a selectivity of ~30%. This performance is not seen for other heterogenous catalysts which usually require the addition of water to enable a significant conversion of methane to methanol. The unique coarse structure of the ZnO islands supported on a Cu 2 O/Cu(111) substrate provides a collection of multiple centers that display different catalytic activity during the reaction. ZnO-Cu 2 O step sites are active centers for methanol synthesis when exposed to CH 4 and O 2 due to an effective O-O bond dissociation, which enables a methane-to-methanol conversion with a reasonable selectivity. Upon addition of water, the defected O-rich ZnO sites, introduced by Zn vacancies, show superior behavior towards methane conversion and enhance the overall methanol selectivity to over 80%. Thus, in this case, the surface sites involved in a direct CH 4 → to CH 3 OH conversion are different from those engaged in methanol formation without water. Furthermore, the identification of the site-dependent behavior of ZnO/Cu 2 O/Cu(111) opens a design strategy for guiding efficient methane reforming with high methanol selectivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The effect of triangularity on predicted pedestals for the CFETR

In this work, the dominant unstable magnetohydrodynamic mode and its stable region in the pedestal of the China Fusion Engineering Test Reactor are studied by numerical modeling with the peeling-ballooning theory over a wide range of triangularities (δ) and collisionalities (υ*). A new accessible stable region is found at δ<−0.1 for large βN,ped and υ*. This new stable region expands toward lower βN,ped and υ* with decreasing δ and is totally covered by the peeling unstable region with δ=−0.5 because of the increasing trapped particle fraction (ft,eff). The sensitivity of this new stable region to the kinetic ballooning mode constraint and elongation is studied. For negative and low δ<0.1, the boundary of the first stable region is determined from the ballooning mode. For δ>0.1, the peeling mode becomes dominant as the boundary approaches low s with low υ*, while the ballooning mode is still dominant at the boundary with high υ*. When δ increases beyond 0.46, the first stable region is expanded, and access to the second stable region of the ballooning mode opens up. The entire newly connected region of the first and second stable regions keeps expanding with further increases in δ until δ=0.6. Beyond this value, the ballooning mode becomes more unstable with increasing δ, while the peeling mode is approximately unchanged with increasing δ in this range. The change in the dominant mode and the stable region with increasing δ can be explained by the trade-off between the stabilization effect from the deeper poloidal magnetic well and destabilization due to the enlarged drive term.

Physics↗
DOE OSTI2021

Parametric Study of Used Nuclear Oxide Fuel Constituent Dissolution in Molten LiCl-KCl-UCl 3

Prior work identified dissolution of used nuclear oxide fuel constituents from a uranium oxide matrix into molten LiCl-KCl-UCl 3 at 500°C, prompting a subsequent series of three progressive studies (including an initial scoping study, an electrolytic dissolution study, and a chemical-seeded dissolution study) to further investigate associated parameters and mechanisms. Thermodynamic calculations were performed to identify possible reaction mechanisms and their propensities in used oxide fuel constituent dissolution. Used nuclear oxide fuels with varying preconditions from fast and thermal test reactors were separately immersed in the subject salt system to assess fuel constituent migration from the bulk fuel matrix to the salt phase in an initial scoping study. Dissolution of expected fuel constituents, including alkali, alkaline earth, lanthanide, and transuranium oxides, into the chloride salt phase varied widely, ranging from 12% to 99% in the initial study. Uranium isotope blending between the salt phase and bulk fuel matrix was also observed, which was attributed to reducing conditions in the fuel matrix. Electrolytic and chemical-seeded dissolution studies were subsequently performed to effect reducing conditions in the fuel. Other parameters, including temperature (at 500°C, 650°C, 725°C, and 800°C) and uranium trichloride concentrations (at 6, 9, and 19 wt% uranium), were investigated in the latter two studies, resulting in fuel constituent dissolution above 90%. Extents of dissolution were based on initial and final fuel constituent concentrations in the oxide fuels following operations in the salt and subsequent removal of the salt via distillation. Finally, in this series of progressive studies, oxide fuel preconditioning and in situ reducing conditions, along with elevated temperature and uranium trichloride concentrations, were the primary parameters promoting used nuclear oxide fuel constituent dissolution in accordance with identified reaction mechanisms.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Low- n global ideal MHD instabilities in the CFETR baseline scenario

This article reports an evaluation on the linear ideal magnetohydrodynamic (MHD) stability of the China Fusion Engineering Test Reactor (CFETR) baseline scenario for various first-wall locations. The initial-value code NIMROD and eigen-value code AEGIS are employed in this analysis. Despite the distinctly different approaches in modeling the scrape off layer region, the dominant growth in each of the low-n (n=1–10) modes are consistent between the two codes. We report the higher-n modes are dominated by ballooning modes and localized in the pedestal region, while the lower-n modes have more prominent external kink components and broader mode profiles. The influence of the plasma-vacuum profile and wall shape are also examined using NIMROD. In the presence of a resistive wall, the low-n ideal MHD instabilities are further studied using AEGIS. For the designed first-wall location, the n=1 resistive wall mode is found to be unstable, which could be fully stabilized by uniform toroidal rotation above the 2.9% core Alfven speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Grassy ELM regime at low pedestal collisionality in high-power tokamak plasma

Spontaneous mitigation of giant ELMs and appearance of grassy ELMs have been observed repeatedly at low pedestal collisionality (${\nu }_{\mathrm{e},\mathrm{p}\mathrm{e}\mathrm{d}}^{{\ast}}\sim 0.15$) in the high-power (P inj > 13 MW) hybrid scenario in DIII-D. Higher β p and higher q 95 appear to be beneficial to achieving the grassy ELM regime. The grassy ELM H-mode plasma shows high energy confinement performance (H 98y2 up to over 1.5) under the conditions of high neutral beam torque and high core rotation. The pedestal width appears to exceed the EPED1.0 model prediction by more than 50%. Pedestal stability analysis performed with the ELITE code indicates that the stability against low-n kink/peeling modes is improved with increased auxiliary heating power and the operational point in the grassy ELM regime is located near the ballooning boundary. The pedestal stability characteristics during the grassy ELM crashes have been investigated in comparison with the giant ELM crashes based on plasma profiles experimentally measured with high resolution and accuracy. It has been found that the underlying mechanism for the observed small-amplitude ELM crashes is mainly the expansion of the ballooning stability boundary induced by an initial radially localized collapse in the pedestal, which helps to stop the growth of instabilities and further collapse of the pedestal. The effect of electron density pedestal on mitigating edge localized instabilities has been analyzed by numerical simulation, suggesting that the electron density pedestal characterized by high n e,sep /n e,ped and low density gradient helps to stabilize peeling-ballooning modes because of a low pressure pedestal gradient and to lower the ballooning boundary mainly because of a low ion diamagnetic frequency in the pedestal region, thus triggering ballooning instabilities and producing the intrinsic grassy ELMs. Numerical simulation of the Chinese fusion engineering test reactor (CFETR) with the SOLPS code indicates that the separatrix density might be insensitive to the electron diffusivity in the pedestal region and increase with the power flowing from the core region to the edge region. Furthermore, pedestal stability analysis suggests that the flat density pedestal with high separatrix density obtained in the high-power plasma in CFETR would make the operational point close to the ballooning boundary, which is considered to help destabilize ballooning instabilities and facilitate the access to the grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Scoping study of lower hybrid current drive for CFETR

The paper assesses the applicability of lower hybrid current drive (LHCD) for two potential operating scenarios for the China Fusion Engineering Test Reactor (CFETR): the “hybrid” scenario in which some of the plasma current is sustained by the Ohmic transformer, and the fully non-inductive “steady state” scenario. Here, the πScope workflow engine was used to set up a large number of ray tracing/Fokker-Planck simulations (> 10 4 ) with parametric scans in the antenna poloidal position and launched parallel refractive index (n || ) for both the hybrid and steady state scenarios. Modeling predicts efficient off-axis current drive (1.3 MA for 20 MW launched power) with a peak near ρ of 0.6-0.65 for waves launched from the high field side (HFS). Waves launched from the low field side (LFS) damp at larger radius (ρ > 0.73) with similar efficiency to HFS launch. Stability analysis of the CFETR scenarios favors current drive profiles peaked near the mid-radius, suggesting that HFS launch is preferable due to the current drive location. The effect of wave scattering from density blobs in the edge/scrape-off-layer region was assessed through rotation of the perpendicular wavenumber at the ray origin. Simulations show that this effect can be quite large both in efficiency and damping location, however by adjusting the launched n|| much of the unperturbed performance can be recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Theoretical study of the Alfven eigenmode stability in CFETR steady state discharges

The aim of this study is to analyze the stability of Alfven eigenmodes (AE) in the China Fusion Engineering Test Reactor (CFETR) plasma for steady state operations. The analysis is done using the gyro-fluid code FAR3d including the effect of the acoustic modes, EP finite Larmor radius damping effects and multiple energetic particle populations. Two high poloidal β scenarios are studied with respect to the location of the internal transport barrier (ITB) at r/a ≈ 0.45 (case A) and r/a ≈ 0.6 (case B). Both operation scenarios show a narrow TAE gap between the inner-middle plasma region and a wide EAE gap all along the plasma radius. In this work, the AE stability of CFETR plasmas improves if the ITB is located inwards, case A, showing AEs with lower growth rates with respect to the case B. The AEs growth rate is smaller in the case A because the modes are located in the inner-middle plasma region where the stabilizing effect of the magnetic shear is stronger with respect to the case B. Multiple EP populations effects (NBI driven EP + alpha articles) are negligible for the case A, although the simulations for the case B show a stabilizing effect of the NBI EP on the n = 1 BAE caused by α particles during the thermalization process. If the FLR damping effects are included in the simulations, the growth rate of the EAE/NAE decreases up to 70%, particularly for n > 3 toroidal families. Low n AEs (n < 6) show the largest growth rates. On the other hand, high n modes (n = 6 to 15) are triggered in the frequency range of the NAE, strongly damped by the FLR effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗