Engineering Papers⌕ Search

Engineering topics

Maingi, Rajesh

Publications and source records attributed to Maingi, Rajesh.

At least 19 records

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U – A national User Facility to close the key fusion gaps

This presentation lays out how the NSTX-U Research Program will help close the science and technology gaps for moving towards a Fusion Pilot Plant (FPP) and how this work is aligned with the Bold Decadal Vision. The science gap focuses on determining whether the good confinement trends in Spherical Tokamaks extends to lower collisionalities that are comparable to those expected in FPPs. The technology gap is one of mitigating the power exhaust in FPPs, and NSTX-U is uniquely positioned to test Liquid Metal components in the highest heat flux environment of any existing tokamak, spherical or otherwise.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparative study of boron and neon injections on divertor heat fluxes using SOLPS-ITER simulations

Here, based on the EAST equilibrium, the effects of boron (B) and neon (Ne) injected at different locations on the target heat load, and the distributions of B and Ne particles were investigated by transport code SOLPS-ITER. It was found that the B injection was more sensitive to the injection location for heat flux control than impurity Ne. The high electron and ion densities near the inner target in the discharge with impurity B injected from over X-point (R 1 ) led to plasma detachment only at the inner target, and the localized B ions in the cases with injection from outer target location (R 2 ) and upstream location (R 3 ) led to far-SOL detachment at the outer target, but not at the inner target. In contrast, for Ne, the spatial distributions of Ne ions and electrons were found to be similar in all the cases at the three injection locations, and the detached plasma was achieved at the inner target and the electron temperature was reduced at the outer target. For locations R 2 and R 3 , impurity B showed a more pronounced effect on the heat flux at the far-SOL of the outer target. Further analysis indicated that Ne atoms came mainly from the recycling sources, whereas B atoms came mainly from injection, and that their distinct atomic distributions resulted from the difference in the ionization threshold and ionization mean free path. In addition, the radiation proportion of B in the divertor region was larger than that of Ne when the total radiation power was similar, which suggests that B has less influence on the core region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Coupled model for liquid lithium plasma facing components

Numerical analysis provides the design choice and operating window of liquid metal Plasma Facing Components (PFC) concepts. Coupled analysis of boundary plasma together with the surrounding boundary structures is required. Here, to achieve this goal, PPPL is developing a comprehensive multi-physics model for modeling of PFCs in fusion devices. The model includes the fluid-kinetic code SOLPS-ITER and the flow and heat transfer code CFX from ANSYS. SOLPS-ITER was augmented with a liquid metal boundary condition algorithm, allowing direct two-way coupling of the plasma analysis with the two-dimensional analytical slab flow model which includes heat convection in the liquid metal PFC. The target heat flux resulting from this coupled analysis is used as a boundary condition for detailed 3D Computational Fluid Dynamics (CFD) Magneto Hydro Dynamics (MHD) and heat transfer analysis. A new formulation of MHD equations is introduced in the numerical procedure ensuring current conservation of the discretized equations. Results of the 3D analysis are used for final validation of the coupled model. A PFC design where a porous wall is used to stabilize the liquid metal surface, while MHD drive is used to push the liquid metal flow inside the PFC, will be investigated in the regimes where vapor shielding is created for enhanced volumetric plasma heat dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NSTX-U National Research Program: White Paper in Response to Call from FESAC Sub-Committee

Both scientific and technical innovation is needed for the realization of an attractive engineering solution for a timely and cost-effective Pilot Plant, the design and construction of which is the overarching recommendation of the FESAC Long Range Plan, and the 2021 NASEM Pilot Plant reports, which underpin the Bold Decadal Vision. The two most significant plasma physics gaps to close for a Compact Pilot Plant (CPP) are core confinement improvement and heat flux mitigation, neither of which have been closed in an integrated fashion for any planned fusion power production device. High core confinement and stability are essential for producing majority self-driven plasmas in CPPs with reduced size and auxiliary heating power requirements, with an improvement in confinement being the major driver for cost reduction of a CPP. The National Spherical Tokamak Experiment - Upgrade (NSTX-U) is a unique low aspect ratio research facility that will address the fundamental challenge of developing the science and technology basis for a CPP design that integrates high core and edge confinement with the ability to mitigate very high incident heat fluxes. NSTX-U capabilities will enable the high performance, already achieved on NSTX, to extend into physics regimes much closer to those anticipated in Spherical Tokamak (ST)-based CPPs. These confinement and stability properties will be assessed by a full complement of diagnostics and analysis tools, which will also aid in the development of the underlying theory and predictive models needed for further optimization. Both conventional and transformative heat flux mitigation methods, such as liquid lithium plasma-facing components, will be developed and tested in-situ in NSTX-U at incident heat fluxes of ~100 MW/m 2 , and will inform plans and reduce risk for a subsequent major upgrade to the device to fully heated, high-Z wall and full liquid lithium divertor capability, a technology that potentially could then be implemented on any magnetic confinement device at any aspect ratio. NSTX-U research is fully complementary to programs performed on other STs, nationally and internationally. Furthermore, NSTX-U research has a direct connection to the private sector by informing design choices for future power production facilities being developed by these companies. The NSTX-U program will operate as a national User Facility, with collaborating researchers, engineers, and graduate students from 19 outside institutions, and open to participation and experiments led by researchers from both public and private entities. The research program will advance workforce development through training of young scientists, engineers, and technicians, and it will also serve for further diagnostic innovation, especially for high heat flux and high-Z wall environments, and implementation of advanced artificial intelligence (AI) for plasma and heat flux control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Report of the FESAC Facilities Construction Projects Subcommittee

On December 1, 2023, the Director of the of the Office of Science (SC) charged all the Department of Energy Office of Science Federal Advisory Committees to respond to look toward the scientific horizon and identify what new or upgraded facilities will best serve the SC community needs in the next ten years (2024-2034). This report by the Fusion Energy Sciences Advisory Committee (FESAC) assess facilities for the Fusion Energy Sciences (FES) program.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and Analysis of Liquid Lithium Plasma Facing Components

In this study, liquid metal (LM) on plasma-facing components (PFCs) creates a renewable protective cover enhancing power exhaust and protecting the solid substrate. Additionally, energy confinement improvement through particle pumping can be achieved in the case of liquid lithium (LL). A new variant of PFCs was recently introduced at Princeton Plasma Physics Laboratory (PPPL) where a porous wall is used to stabilize the LM surface, while using magnetohydrodynamic (MHD) drive to push the LM flow inside the component. This arrangement allows efficient heat exhaust, and enhanced control of the LM surface temperature, leading to spatial control of evaporation and condensation of LL on the plasma interface. This feature is particularly attractive when vapor shielding is introduced to allow heat flux redistribution. This system has the advantage that as the heat flux increases, the evaporation rate will increase, while decreasing heat flux will decrease the evaporation rate of lithium, ideally creating a feedback effect which could self-regulate the amount of lithium evaporated. Analytical and numerical models for LL PFCs were developed in PPPL. To calculate the target temperature distribution for the case of evaporation from the divertor, we apply an iterative process allowing two-way coupling between the fluid-kinetic analysis of plasma using SOLPS-ITER code and the flow and heat transfer analysis of the PFC using an analytical model. In the last stage, the results are validated using computational fluid dynamics (CFD) analysis with customized version of the CFX code from ANSYS. CFX was modified at PPPL to allow MHD analysis at the high magnetic field typical for fusion applications. Results of the analysis for NSTX-U tokamak conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis and design of fast flow liquid Li divertor for fusion nuclear science facility (FNSF) using coupled plasma boundary and LM MHD/heat transfer codes *

The SOLPS-ITER code is utilized to analyze the boundary plasma associated with a fast-flow lithium (Li) divertor configuration in the fusion nuclear science facility (FNSF) tokamak and identify operational regimes with acceptable divertor and core conditions. Plasma transport from the SOLPS-ITER code has been coupled with a liquid metal (LM) MHD/heat transfer code to model a Li open-surface divertor design and assess its impact on the scrape-off-layer (SOL) and core plasma performance. Simulations with only Neon (Ne) impurity seeding have been conducted to evaluate its impact on meeting FNSF design demands for the divertor and upstream plasma parameters. Simulation results indicate that Ne seeding significantly mitigates divertor heat flux but potentially reduces both upstream electron and main ion density due to fuel dilution. The combined application of Ne seeding and deuterium (D 2 ) puffing is required to satisfy the FNSF design requirements on upstream density ($n^{OMP}_{e,sep}$~1× 10 20 m -3 ) and divertor energy flux ($q^{Odiv}_{\bot,max}$<10 MW m -2 ). D 2 puffing plays a role in counteracting upstream density drops and augmenting energy and momentum losses through atomic and molecular processes. The inlet Li flow velocity is systematically varied across a wide range to identify acceptable flows and corresponding LM surface temperatures. This comprehensive analysis identifies the acceptable Li flow parameters, LM surface temperature, and emitted Li fluxes necessary to meet the major design constraints. The emitted Li fluxes exhibit minimal impact on the main plasma at surface temperatures up to approximately ~525 °C, corresponding emitted Li fluxes of up to φ Li ~2X10 23 atoms s -1 . Uncertainties in the Li emission processes from the surface are also investigated, primarily influencing Li loss in the lower surface temperature range (<525 °C), with simulation results indicating a minor impact on the divertor and upstream plasma. Conversely, evaporation predominantly drives the Li loss processes at higher surface temperature ranges (>525 °C), contaminating both the divertor and upstream plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conference Report on the 7th International Symposium on Liquid metals Applications for fusion (ISLA-7)

Supported by the world magnetic fusion research community, a series of International Symposia on Liquid metals Applications for fusion (ISLA) have been held biannually since 2010. The 7th edition (ISLA-7) was held for the period from 12 December through 16 December 2022, at Chubu University located in Kasugai, Aichi, Japan. For the first time in the history of this series of symposia, ISLA-7 was held in a hybrid fashion, due to the COVID-19 situation. The total number of the participants was 60, 34 out of whom attended the symposium in person, and the rest participated online. As to the presentation statistics, 29 papers were presented in person, whereas 21 presentations were delivered online but real-time by the presenters in China, Spain, the UK, and the USA. Both of the presentations delivered in person and online were recorded, and the video has been shared by all participants. These participants represent 11 countries: China, Czech, Italy, Japan, Latvia, Netherlands, Russia, Thailand, the UK, and the USA. All these numbers are among the largest in this series of symposia. Covered by these presentations are; in session-2, program overviews and liquid metal research review; in session-3, liquid metal flows, and MHD issues; in session-4, liquid metal facilities; in sessions-5 and 6, liquid metal experiments and modeling; in session-7, divertor physics and heat flux mitigation; in session-8, plasma and liquid metals interactions; in session-9 liquid metal plasma-facing components, erosion, and wettability. In addition, there were an opening session whereby several opening addresses were delivered and also a closing session where all technical session summaries were presented by the respective session chairs.

divertor↗

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Realization of thousand-second improved confinement plasma with Super I-mode in Tokamak EAST

Mastering nuclear fusion, which is an abundant, safe, and environmentally competitive energy, is a great challenge for humanity. Tokamak represents one of the most promising paths toward controlled fusion. Obtaining a high-performance, steady-state, and long-pulse plasma regime remains a critical issue. Recently, a big breakthrough in steady-state operation was made on the Experimental Advanced Superconducting Tokamak (EAST). A steady-state plasma with a world-record pulse length of 1056 s was obtained, where the density and the divertor peak heat flux were well controlled, with no core impurity accumulation, and a new high-confinement and self-organizing regime (Super I-mode = I-mode + e-ITB) was discovered and demonstrated. These achievements contribute to the integration of fusion plasma technology and physics, which is essential to operate next-step devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗

Plasma facing components with capillary porous system and liquid metal coolant flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 °C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently, a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [A. Khodak and R. Maingi, Nucl. Mater. Energy 26, 100935 (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In the case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a fusion nuclear science facility [Kessel et al., Fusion Sci. Technol. 75, 886 (2019)]: 10 T toroidal field and 10 MW/m 2 peak incident heat flux. Finally, computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma Facing Components with Capillary Porous System and Liquid Metal Coolant Flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [Khodak and Maingi (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile, as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a Fusion Nuclear Science Facility [Kessel et al. (2019)]: 10T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

High Heat Flux↗

Critical role of current-driven instabilities for ELMs in NSTX

The impact of different extended-magnetohydrodynamic (MHD) contributions on the stability thresholds of peeling-ballooning modes in ELMing and ELM-free plasmas in the spherical tokamak NSTX is investigated with the initial value code M3D-C1. We show that ELMing discharges in NSTX are limited by resistive current-driven peeling modes, whereas non-ELMing wide-pedestal H-mode discharges are located near the ideal pressure-driven ballooning threshold. It is demonstrated that extended-MHD can lead to more reliable edge stability predictions than existing ideal-MHD models. Resistive peeling-ballooning modes are found to exist well before the ideal stability threshold is met, and kink-peeling modes exhibit considerable sensitivity to plasma resistivity. Other effects not considered in ideal-MHD models affect PB modes in NSTX in a weaker way. Gyroviscous stress appears stabilizing such that the stability boundary lies closer to the experimental point. Equilibrium rotation can suppress ideal core modes and thus isolate edge modes. These results are important for the development of a predictive pedestal model for low-aspect ratio tokamaks.

pedestal↗

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

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 60s with β P ~2.0, β N ~1.6, H 98y2 ~1.3 and internal transport barrier on electron temperature channel is obtained with multi-RF power heating and current drive. A higher β N (β N ~1.8, β p ~2.0, H 98y2 ~1.3, n e /n GW ~0.75) with a duration of 20s 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 1GW 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. Modelling shows that the strong diffusion of TEM turbulence in the central region prevents tungsten impurity to accumulate. More recently, EAST has demonstrated compatible core-edge integration discharges in the high βp scenario: high confinement H 98y2 >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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of Liquid Lithium Divertor Geometry Using SOLPS-ITER

Here, plasma-facing component (PFC) geometries are evaluated for a Fusion Nuclear Science Facility (FNSF) design to find solutions that are compatible with flowing liquid lithium components while satisfying requirements from the perspective of plasma and neutral particle transport. Flowing liquid metal (LM) divertor systems offer important advantages due to continual regeneration of the surface material, but may require modifications from standard optimized tokamak divertor designs. Scrape Off Layer Plasma Simulator for the International Thermonuclear Experimental Reactor (SOLPS-ITER) simulations are used to compare standard vertical target, open, and balanced baffled geometries. It is found that the open geometry offers little control of neutral particles and has narrow operating windows where the divertor fluxes can be lowered to acceptable levels, while maintaining the required upstream density. The addition of neutral baffles allows greater control over the divertor conditions with less impact on the upstream density. The divertor neutral pressure in the baffled geometry is slightly higher than the open divertor, but not as large with a vertical target. A first-pass coupling of the SOLPS solution with sheath and implantation/erosion models to compute the lithium emission shows that the lithium is largely confined to the near-surface region due to the small ionization mean free path and large main ion flow toward the surface. The lithium source was varied by three orders of magnitude. At the highest rate, when the sourced lithium approaches the level of the main ion recycling flux, a high core lithium concentration (~10%) is observed with only moderate power dissipation to affect the divertor plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗