Engineering topics
Rognlien, T. D.
Publications and source records attributed to Rognlien, T. D..
2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries
UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.
UEDGE modeling of plasma detachment of CFETR with ITER ‐like divertor geometry by external impurity seeding
Abstract Efficient handling of high heat flux on the plasma‐facing components, particularly the divertor targets, poses a significant challenge for the Chinese Fusion Engineering Testing Reactor (CFETR) with fusion power of Gigawatt. This work investigates the divertor plasma detachment of CFETR with a standard ITER‐like divertor geometry by neon (Ne) or argon (Ar) impurity seeding using UEDGE code. The cross‐field drifts terms are switched off, and fluid neutral models and a “fixed‐fraction” impurity model are applied to enable efficient simulations for the study of CFETR detachment. In order to reduce the heat load on the divertor targets below the acceptable level (<10 MW/m 2 ), the impurity fraction ( f ), pumping speed ( S ), and upstream density are varied to identify the suitable operations window during Ne seeding. The effects of Ne and Ar impurities on the plasma detachment are compared. It is found that with the power across the core‐edge interface P SOL = 200 MW and separatrix density of 2.8 10 19 , Ne impurity fraction ≥1.7%, and Ar impurity fraction ≥0.24% can achieve the partial detachment. Achieving similar total radiation power (˜148 MW), the Ne fraction is 2.3% and the Ar fraction is 0.24%. Moreover, the simulation results indicate that Ar exhibits better power radiation efficiency and core compatibility compared with Ne.
Modeling of deuterium and carbon radiation transport in MAST-U tokamak advanced divertors
Modest effects of deuterium and carbon radiation opacity in the Super-X and snowflake divertor plasmas are predicted for MAST Upgrade tokamak with core plasma input power 2.5–5 MW and plasma current 1 MA. The radiation transport modeling is based on the SOLPS-EIRENE and UEDGE code divertor plasma predictions. Two radiation transport models are used: one is based on a full radiation transport equation implemented in the radiation transport and collisional-radiative code CRETIN (without feedback on the background plasma), and another is an internal self-consistent UEDGE model with ionization, recombination, and heating rates corrected for Ly α line trapping based on the escape probability model implemented in CRETIN. In MAST-U, the Super-X and snowflake divertor plasmas are predicted to reach detached regimes at lower upstream densities than the standard divertor, and the conclusion still holds with radiation transport effects included. At neutral densities m -3 , modest Ly α deuterium line trapping with optical depths 10–15 is predicted in the Super-X divertor. Divertor plasmas are optically thin to other Lyman and Balmer lines, as well as to strong C III and C IV lines that are responsible for most of divertor radiated power. Insignificant changes (within a few percent) to divertor deuterium ionization and recombination rates are found. Radiation fluxes on outer divertor target are modified within a factor of 2–3 when the radiation transport is accounted for, and a similar variation is found due to the line shape models that define the absorption and emission line profiles in the radiation transport modeling. The predicted Lyman and Balmer spectral intensities are significantly modified due to radiation trapping. A measurement of divertor radiation transport effects is discussed using the Ly β /Ba α line ratio. In the snowflake divertor configuration, divertor plasmas are found to be optically thin to Lyman series lines within a large range of parameter variations that include upstream density, divertor transport coefficients, and magnetic configurations. Finally, modest radiation transport effects are only found in a few cases with strongest divertor transport and magnetic configurations closest to the ideal snowflake configuration, however, the plasma background models that were used are yet to be validated with an experiment.
Modeling Snowflake Divertors in MAST-U Tokamak
We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.
FESS Design Simulations: methods, tools, & issues for the edge/scrape-off-layer region
This topic describes the plasma and neutral particles in the transition boundary region between the hot core plasma and the surrounding material walls. A key geometrical transition that occurs in this region is where the equilibrium magnetic field changes topology from a set of closed, nested magnetic flux surfaces inside the magnetic separatrix to flux surfaces, and therefore magnetic field lines that intersect material walls. Because the plasma exhaust heat flows very rapidly along the field lines, these intersection locations can have heat fluxes much higher than the walls can withstand. The most promising strategy pursued here to avoid this problem is injection of moderate-Z impurities that radiate the exhaust power over a much larger surface area on the walls, thus keeping the peak heat flux at or below the goal of 10 MW/m 2 . Further, the intrusion of injected and wall-sputtered impurities into the core region must be kept below certain limits to prevent degradation of the fusion power generated in the core.
Edge and scrape-off layer modeling for a Fusion Nuclear Science Facility with tungsten walls; a summary report for 2019-21
This report summarizes model development and simulations for the edge/scrape-off layer (SOL) region of a Fusion Nuclear Science Facility (FNSF) as part of the DOE Fusion Energy Systems Studies project. An overview of the FNSF device is given in Ref. 1. Our earlier related modeling of FNSF in the 2015-16 timeframe is reported in Ref. 2, and similar work on the ARIES ACT-1 tokamak device is described in Ref. 3. During 2017-18, we contributed to the analysis of a liquid lithium wall for FNSF [4].
Comparison of a collisional-radiative fluid model of H2 in UEDGE to the kinetic neutral code EIRENE
Not Available
Overview of the fusion nuclear science facility, a credible break-in step on the path to fusion energy
The Fusion Nuclear Science Facility (FNSF) is examined here as part of a two step program from ITER to commercial power plants. This first step is considered mandatory to establish the materials and component database in the real fusion in-service environment before proceeding to larger electricity producing facilities. The FNSF can be shown to make tremendous advances beyond ITER, toward a power plant, particularly in plasma duration and fusion nuclear environment. A moderate FNSF is studied in detail, which does not generate net electricity, but does reach the power plant blanket operating temperatures. The full poloidal Dual Coolant Lead Lithium (DCLL) blanket is chosen, with alternates being the Helium Cooled Lead Lithium (HCLL) and Helium Cooled Ceramic Breeder/Pebble Bed (HCCB/PB). Several power plant relevant choices are made in order to follow the philosophy of targeted technologies. Any fusion core component must be qualified by fusion relevant neutron testing and highly integrated non-nuclear testing before it can be installed on the FNSF in order to avoid the high probability of constant failures in a plasma-vacuum system. A range of missions for the FNSF, or any fusion nuclear facility on the path toward fusion power plants, are established and characterized by several metrics. A conservative physics strategy is pursued to accommodate the transition to ultra-long plasma pulses, and parameters are chosen to represent the power plant regime to the extent possible. An operating space is identified, and from this, one point is chosen for further detailed analysis, with R = 4.8 m, a = 1.2 m, IP = 7.9 MA, BT = 7.5 T, βN Gr = 0.9, fBS = 0.52, q95 = 6.0, H98 ∼1.0, and Q = 4.0. The operating space is shown to be robust to parameter variations. A program is established for the FNSF to show how the missions for the facility are met, with a He/H, a DD and 5 DT phases. The facility requires ∼25 years to complete its DT operation, including 7.8 years of neutron production, and the remaining spent on inspections and maintenance. The DD phase is critical to establish the ultra-long plasma pulse lengths. The blanket testing strategy is examined, and shows that many sectors have penetrations for heating and current drive (H/CD), diagnostics, or Test Blanket Modules (TBMs). The hot cell is a critical facility element in order for the FNSF to perform its function of developing the in-service material and component database. The pre-FNSF R&D is laid out in terms of priority topics, with the FNSF phases driving the time-lines for R&D completion. A series of detailed technical assessments of the FNSF operating point are reported in this issue, showing the credibility of such a step, and more detailed emphasis on R&D items to pursue. These include nuclear analysis, thermo-mechanics and thermal-hydraulics, liquid metal thermal hydraulics, transient thermo-mechanics, tritium analysis, maintenance assessment, magnet specification and analysis, materials assessments, core and scrape-off layer (SOL)/divertor plasma examinations.
Impact of ion temperature anisotropy on 2D edge-plasma transport
Not Available
Modeling snowflake divertors in MAST-U tokamak using UEDGE code
Not Available
Summary of 2020 edge-plasma modeling for FNSFand plans for 2021
This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.
Low-frequency flute instabilities of a bounded plasma column.
Derivation of exact solutions for unstable waves (called flute waves) which occur in a radially bounded plasma column at frequencies below the ion cyclotron frequency. Both analytical and numerical solutions are presented for the m = 1 and m = 2 azimuthal modes for a variety of radial electric field profiles. It is shown that the behavior of the flute waves can depend sensitively on the radial extent of the plasma column. Moreover, it is found that the m = 1 mode and the m = 2 mode do not respond in the same way to changes in the radial boundary position or in the electric field profile.
Low-frequency flute instabilities of a hollow cathode arc discharge - Theory and experiment.
The characteristics of two low-frequency electrostatic flute instabilities of a low-pressure hollow cathode arc discharge are reported. Mode I has azimuthal mode number m = 1, and occurs when the radial electric field is negative (directed inward), while mode II has m = - 1 and occurs when the field is positive. The radial electric field is controlled by varying the potential of a secondary anode cylinder located close to the outer discharge radius. A linear perturbation analysis, based on the two-fluid equations, is given for a low-beta, collisionless, cylindrical plasma column, immersed in a uniform axial magnetic field, having a Gaussian density profile and an arbitrary radial electric field profile. Reasonable correlation between theory and experiment is demonstrated for both modes.
Low-frequency macroscopic instabilities of fully ionized magnetoplasma
Studies are described of low-frequency quasi-static instabilities in a fully ionized plasma. The plasma is assumed to be immersed in a uniform magnetic field, and is either uniform or has a number density gradient perpendicular to the magnetic field. A moment equation description of the ion and electron dynamics is used; collisions are assumed to have a strong effect on electron motion along the magnetic field. Before considering specific modes, a stability analysis is developed which allows a classification of wave growth characteristics to be made for a bounded system from solutions to the dispersion relation for an infinite system. Also, a method is given for calculating the normal mode frequencies and wave profiles by using the reflection coefficients at the boundaries. For wave propagation perpendicular to the magnetic field, the flute wave is studied in cylindrical geometry. The destabilizing effect of a radial electric field is considered by solving a differential equation.
Interpretation of dispersion relations for bounded systems.
Treatment of the problem of constructing normal modes for an arbitrarily bounded system from roots of the linear dispersion relation D(omega, k) = 0 for the corresponding infinite or periodically bounded system. For a system described by continuous macroscopic variables, and of general cylindrical form, each transverse eigenmode gives rise to a set of axial normal modes constructed from a pair of dominant roots of D = 0 satisfying the boundary conditions which are characterized by complex reflection coefficients for the dominant waves. The implications of the results for the interpretation of experiments on plasma waves and instabilities on finite cylinders are discussed, with particular reference to the effects of end-plate damping and axial current on Q-machines.
Ion-acoustic instability of a two-temperature, collisional, fully ionized plasma
Long wavelength ion acoustic instability of two temperature collisional fully ionized plasma with heat transfer, noting additional destabilizing currents