Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ION CYCLOTRON”

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 325 records · Page 18

Impact of suprathermal ions on neutron yield in the pre-DT phase of ITER operation

An assessment of neutron production during the pre-fusion power-operation (PFPO) phase has been carried out for a representative set of plasma scenarios predicted by the ITER Research Plan. A range of heating systems, namely neutral beam injection (NBI) (hydrogen), electron cyclotron resonance heating (ECRH), and ion cyclotron resonance heating (ICRH) are planned to be used for PFPO studies in helium, hydrogen, and mixed hydrogen–helium plasmas. Fast ions (protons and 3He) originating from NBI and ICRH systems can increase neutron production in PFPO plasmas by directly interacting with intrinsic Be impurities or through secondary processes, as also evidenced at JET. The generation of fast ions in ITER PFPO scenarios has been modelled using the ASTRA-NBI and TORIC-SSFPQL codes. A significant impact of the synergy between hydrogen NBI and hydrogen-minority ICRH on neutron production in helium plasmas is reported. In addition, the stability of the toroidicity-induced Alfvén eigenmodes (TAE) is analyzed for PFPO plasmas with a high pressure of suprathermal ions and a weak reversed shear. The possible impact of sawtooth oscillations and TAEs on neutron production is discussed, based on a linear stability analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗

Harvesting 88 Zr from heavy-ion beam irradiated tungsten at the National Superconducting Cyclotron Laboratory

Tungsten is a commonly used material at many heavy-ion beam facilities, and it often becomes activated due to interactions with a beam. Many of the activation products are useful in basic and applied sciences if they can be recovered efficiently. In order to develop the radiochemistry for harvesting group (IV) elements from irradiated tungsten, a heavy-ion beam containing 88 Zr was embedded into a stack of tungsten foils at the National Superconducting Cyclotron Laboratory and a separation methodology was devised to recover the 88 Zr. The foils were dissolved in 30% hydrogen peroxide, and the 88 Zr was chemically purified from the tungsten matrix and from other co-implanted radionuclides (such as 85 Sr and 88 Y) using strong cation-exchange (AG MP-50) chromatographic resin in sulfuric acid media. The procedure provided 88 Zr in approximately 60 mL 0.5 M sulfuric acid with no detectable radio-impurities. The overall recovery yield for 88 Zr was (92.3 ± 1.2)%. In conclusion, this proof-of-concept experiment has facilitated the development of methodologies to harvest from tungsten and tungsten-alloy parts that are regularly irradiated at heavy-ion beam facilities.

07 ISOTOPE AND RADIATION SOURCES↗

Cyclotron resonance accelerators for industrial applications

Here, this paper describes novel configurations for cyclotron resonance acceleration of electrons and ions that have several attractive features including: a compact robust room-temperature single-cell RF cavity as the accelerator structure; and continuous high current accelerated un-bunched beam output with self-scanning, obviating need for a separate beam scanner. An electron accelerator version, the electron Cyclotron Resonance Accelerator (eCRA), is under development to be an efficient source for high power electron and x-ray beams for medical, research, sterilization, and National Security applications, so as to replace radioactive materials. An ion accelerator version, the ion Cyclotron Auto-Resonance Accelerator (iCARA) is described here, suggesting its potential to produce, as an example, a high-current multi-MeV beam of deuterons which could be highly competitive with that produced either with linacs or cyclotrons. Such a deuteron beam could produce a high flux of fast neutrons via deuteron stripping, for applications including the transmutation of used nuclear fuel, material studies relevant for a fusion reactor inner wall, tritium breeding and medical isotope production. For the high-current, high efficiency simulated performance for eCRA and iCARA as described in this paper, the particle beams produced may not exhibit the low emittance values that are important for most discovery research. Rather, the beams could be useful for industrial applications where higher emittance and some energy spread can be tolerated, in favor of high beam power.

43 PARTICLE ACCELERATORS↗

Time-dependent scenario modeling for the ST-E1 fusion power plant

ST-E1 is a low aspect ratio fusion power plant being designed by Tokamak Energy targeting 1.5 GW of fusion power. Characterization of the ST-E1 flat-top scenario is described elsewhere McNamara et al (2026 Nucl. Fusion 66 086008); here we focus on addressing the question of how to ramp-up the ST-E1 plasma from an initial state following breakdown and flux-surface formation to the target flat-top state. Being low-aspect ratio, the available solenoid flux of ST-E1 is limited. Therefore, particular consideration is placed on developing ramp-up scenarios that predominantly use inductive flux provided by external vertical field coils. Through time-dependent modeling with METIS, we show that this is possible when the ramp-up is performed at relatively high plasma density: although auxiliary current drive efficiency is reduced, this is significantly outweighed by (1) higher electron-ion collisional equilibration, (2) higher fusion power once ions become sufficiently hot, (3) higher poloidal beta for increased vertical-field flux, and (4) potentially favorable exhaust compatibilities. Ultimately, we show the target ST-E1 flat-top performance can be reached after a ramp-up period lasting 150 s using less than 40 Vs of solenoid flux (with vertical field providing ∼ 90 Vs of flux). The sensitivity to model assumptions are presented, with the general observation that deleterious effects can be mitigated through minor alterations of the auxiliary power temporal waveform and/or total auxiliary power level. The impact of a solenoid and the auxiliary power mix (electron cyclotron heating only versus electron and ion cyclotron heating) on the ST-E1 ramp-up success are also discussed in appendices. On this latter topic, we show that the effect of direct-ion heating during ramp-up is obscured by the uncertainty in the pedestal dynamics, identifying a clear line of future work required to make a definite decision on the ST-E1 auxiliary power mix.

current ramp-up↗

Hybrid Simulation of Proton Cyclotron Waves Upstream of Mars Generated by Pickup Ion Beam Distribution

Linear instability analysis as well as a corresponding two‐dimensional hybrid simulation are performed to examine the excitation of the proton cyclotron waves observed upstream of Mars. The waves are believed to be excited by the pickup ions produced from the ionization of the Martian hydrogen exosphere. And previous statistical analysis of wave observations suggested that the waves are mostly related to pickup ion beam velocity distributions. While earlier linear instability analysis of pickup ion beam distributions has mainly been focused on the parallel unstable modes, our analysis reveals that the maximum growth rate occurs at very oblique propagation. The corresponding hybrid simulation confirms the linear analysis results and further demonstrates that the pickup ions are scattered toward an isotropic shell velocity distribution by the waves excited. Interestingly, the waves at oblique propagation gradually damp out and the system is eventually dominated by waves of quasi‐parallel propagation.

79 ASTRONOMY AND ASTROPHYSICS↗

Isotope effects on transport in LHD

Abstract Isotope effects are one of the most important issues for predicting future reactor operations. Large helical device (LHD) is the presently working largest stellarator/helical device using super conducting helical coils. In LHD, deuterium experiments started in 2017. Extensive studies regarding isotope effects on transport have been carried out. In this paper, the results of isotope effect studies in LHD are reported. The systematic studies were performed adjusting operational parameters and nondimensional parameters. In L mode like normal confinement plasma, where internal and edge transport barriers are not formed, the scaling of global energy confinement time ( τ E ) with operational parameters shows positive mass dependence ( M 0.27 ; where M is effective ion mass) in electron cyclotron heating plasma and no mass dependence ( M 0.0 ) in neutral beam injection heating plasma. The non-negative ion mass dependence is anti-gyro-Bohm scaling. The role of the turbulence in isotope effects was also found by turbulence measurements and gyrokinetic simulation. Better accessibility to electron and ion internal transport barrier (ITB) plasma is found in deuterium (D) plasma than in hydrogen (H). Gyro kinetic non-linear simulation shows reduced ion heat flux due to the larger generation of zonal flow in deuterium plasma. Peaked carbon density profile plays a prominent role in reducing ion energy transport in ITB plasma. This is evident only in plasma with deuterium ions. New findings on the mixing and non-mixing states of D and H particle transports are reported. In the mixing state, ion particle diffusivities are higher than electron particle diffusivities and D and H ion density profiles are almost identical. In the non-mixing state, ion particle diffusivity is much lower than electron diffusivity. Deuterium and hydrogen ion profiles are clearly different. Different turbulence structures were found in the mixing and non-mixing states suggesting different turbulence modes play a role.

Physics↗

Ion concentration ratio measurements of ion beams generated by a commercial microwave electron cyclotron resonance plasma source

A commercially available electron cyclotron resonance (ECR) plasma source (GenII Plasma Source, tectra GmbH) is widely used for surface processing. This plasma source is compatible with ultrahigh vacuum systems, and its working pressure is relatively low, around 10 –6 –10 –4 Torr even without differential pumping. Here, we report ion flux concentration ratios for each ion species in an ion beam from this source, as measured by a mass/energy analyzer that is a combination of a quadrupole mass spectrometer, an electrostatic energy analyzer, and focusing ion optics. In this study, the examined beams were those arising from plasmas produced from feed gases of H 2 , D 2 , N 2 , O 2 , Ar, and dry air over a range of input power and working pressures. H 2 (D 2 ) plasmas are widely used for nuclear fusion applications and, hence, the ion concentration ratios of H + , H 2 + , and H 3 + reported here will be useful information for research that applies this plasma source to well-controlled plasma-material interaction studies. Ion energy distributions, stability of operation, and impurity concentrations were also assessed for each of the plasma species investigated.

47 OTHER INSTRUMENTATION↗

Polarization and properties of low-frequency waves in warm magnetized two-fluid plasma

This paper presents the derivation of a general wave dispersion relation for warm magnetized plasma under the two-fluid formalism. The discussion is quite general except for the assumption of low frequency and slow phase speed, for which the displacement current is negligible, under the implicit assumption that the plasma is sufficiently dense to satisfy the condition ωpe>ωce, where ωpe and ωce denote the plasma oscillation frequency and electron gyro frequency, respectively. The present discussion does not invoke charge neutrality associated with the fluctuations although it is implicitly satisfied. The resulting dispersion relation that includes the fluid thermal effects shows that there are three eigen modes, which include those corresponding to ideal MHD, namely, fast, slow, and kinetic Alfvén waves, as well as higher-frequency modes including the ion and electron cyclotron and lower-hybrid resonances. The fluid effects in the ideal MHD wave branches are influenced by the finite Larmor radius scales, and when the wave number in the cross field direction is comparable to these values, the fluid effects become significant. It is found that the Larmor radius should be interpreted in the sense as ion-acoustic gyro-radius instead of ion thermal gyro radius only. That is, it is found that the electrons also contribute to the non-ideal effect associated with the kinetic Alfvén wave. A comprehensive explanation of the polarization of each mode is also presented. The present findings indicate that the polarity may change its sign only for the kinetic Alfvén mode branch and that such a transition is based on the propagation angle. When such a change does take place, it is found that the kinetic Alfvén wave transits to an ion-acoustic mode. For each branch, it is also found that the electric field along the ambient magnetic field is purely transverse.

Physics↗

Compact, Achromatic Non-scaling FFAG Accelerator for HEP, Commercial and Medical Applications (CRADA FRA-2008-0005 Final Report)

Proton and light-ion accelerators have many research and medical applications, but the current state of technology in accelerators severely limits their widespread use: (1) synchrotrons have an intrinsic low duty cycle and become quite large for light ions; and (2) cyclotrons do not have energy variability and require complex superconducting magnets for high energy. A new concept in non-scaling Fixed-Field Alternating-Gradient (FFAGs) has been invented that incorporates both the energy variability of the synchrotron and the high duty cycle of the cyclotron. The concept uses normal-conducting, combined-function magnets that apply only constant (dipole) and linear-gradient (quadrupole) fields to stabilize the accelerator. This project will develop an optimized design for a new non-scaling and cost-effective innovation in FFAG accelerators. The accelerator design will have optics that can stably accelerate protons to 250 MeV (or higher) and carbon ions to =200 MeV, an energy considered ideal for radiation therapy. Phase I will optimize, fully simulate, and demonstrate the feasibility of this accelerator. Commercial Applications and other Benefits as described by the awardee: The fixed-field acceleration concept should eliminate some of the most pronounced technical difficulties, expense, maintenance, and required expertise faced in conventional proton and light-ion accelerators. In addition to the application for high energy physics, these accelerators should have application in cancer treatment, radiopharmaceuticals, and medical isotope production, and materials science.

43 PARTICLE ACCELERATORS↗

Ion temperature clamping in Wendelstein 7-X electron cyclotron heated plasmas

The neoclassical transport optimization of the Wendelstein 7-X stellarator has not resulted in the predicted high energy confinement of gas fueled electron-cyclotron-resonance-heated (ECRH) plasmas as modelled in (Turkin et al 2011 Phys. Plasmas 18 022505) due to high levels of turbulent heat transport observed in the experiments. The electron-turbulent-heat transport appears non-stiff and is of the electron temperature gradient (ETG)/ion temperature gradient (ITG) type (Weir et al 2021 Nucl. Fusion 61 056001). As a result, the electron temperature T e can be varied freely from 1 keV–10 keV within the range of P ECRH = 1–7 MW, with electron density n e values from 0.1–1.5 × 10 20 m –3 . By contrast, in combination with the broad electron-to-ion energy-exchange heating profile in ECRH plasmas, ion-turbulent-heat transport leads to clamping of the central ion temperature at T i ~ 1.5 keV ± 0.2 keV. In a dedicated ECRH power scan at a constant density of $\langle n_{e} \rangle$ = 7 × 10 19 m –3 , an apparent 'negative ion temperature profile stiffness' was found in the central plasma for (r/a < 0.5), in which the normalized gradient ∇T i /T i decreases with increasing ion heat flux. The experiment was conducted in helium, which has a higher radiative density limit compared to hydrogen, allowing a broader power scan. This 'negative stiffness' is due to a strong exacerbation of turbulent transport with an increasing ratio of T e /T i in this electron-heated plasma. This finding is consistent with electrostatic microinstabilities, such as ITG-driven turbulence. Theoretical calculations made by both linear and nonlinear gyro-kinetic simulations performed by the GENE code in the W7-X three-dimensional geometry show a strong enhancement of turbulence with an increasing ratio of T e /T i . The exacerbation of turbulence with increasing T e /T i is also found in tokamaks and inherently enhances ion heat transport in electron-heated plasmas. This finding strongly affects the prospects of future high-performance gas-fueled ECRH scenarios in W7-X and imposes a requirement for turbulence-suppression techniques.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Benchmark between antenna code TOPICA, RAPLICASOL and Petra-M for the ICRH ITER antenna

ITER will be equipped with three plasma heating systems: neutral beam (NB), electron cyclotron (EC), and ion cy-clotron resonance heating (ICRH). The latter consists of two identical ICRH antennas to deliver 20 MW to the plasma (baseline, upgradable to 40 MW). ICRH will play a crucial role in the ignition and sustainment of burning plasmas in ITER. A high fidelity and robust modeling effort to understand the interaction of the IC waves with the scrape-off-layer (SOL) plasma is a very important aspect. Among the main important research topics, we have the assessment of the antenna loading for different plasma scenarios, the role of the lower hybrid resonance in front of the antenna and how to include it in our models, and the RF sheath boundary conditions to evaluate the antenna impurity generation. In this work, we tackle the first of these by reporting on ICRF simulations employing the Petra-M code, which is an electromagnetic simulation tool for modeling RF wave propagation based on MFEM [http://mfem.org] for the ITER ICRH antenna. Moreover, a benchmark between the well tested antenna codes TOPICA, RAPLI-CASOL, which is based on COMSOL [www.comsol.com], and the Petra-M code is also presented. In conclusion, S- and Z-matrices and wave electric field are compared showing an excellent agreement among these codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗

Progress towards 28 GHz operations at the Facility for Rare Isotope Beams (FRIB)

This paper provides an overview on the design, operational parameters as well as recent development and recent performance, of the superconducting electron cyclotron resonance (ECR) ion source at the Facility for Rare Isotope Beams (FRIB). The ion source has been used in support of beam operations since the end of 2022 and was designed to operate at 28 GHz and be capable of delivering high intensity uranium beams to help support beam delivery of up to a final beam power of 400 kW. Currently operating at up to 20 kW many primary beams have been developed and routinely used with the superconducting ECR including uranium. Initially operated at 18 GHz, the ion source has now been coupled to a 28 GHz gyrotron and high intensity beams have been demonstrated ahead of the facility planned power ramp up. text with the article abstract. Replace this text with the article abstract.

Machicoane, G. A. [Michigan State University, East↗

Application of simulations and theory of sub-cyclotron frequency modes to DIII-D

Theoretical and numerical study of Alfvén eigenmodes (AEs) in the sub-cyclotron frequency range explains the observed frequency scaling with beam parameters in experiments across different devices. Global Alfvén (GAEs) and compressional Alfvén eigenmodes (CAEs) are frequently excited during neutral beam injection (NBI) in the National Spherical Torus Experiment (NSTX/NSTX-U), as well as other beam-heated devices such as MAST and DIII-D. These modes are driven unstable through the Doppler shifted cyclotron resonance with the NBI ions and can be excited in ITER due to super-Alfvénic velocities and anisotropy of the fast ion distribution. Numerical models and theory for sub-cyclotron frequency modes, previously developed and used to study the excitation of GAEs/CAEs in the NSTX(-U), have been successfully applied to explain the DIII-D observations. Simulations for DIII-D demonstrate that the modes with ω/ω ci ≈ 0.6, previously mis-identified as compressional Alfvén eigenmodes, have shear polarization and should be identified as the GAEs. Simulation results match the observed frequencies and estimated toroidal mode numbers in DIII-D experiments. Good agreement is found between scaling of the experimentally observed GAE frequencies with NBI parameters for NSTX, NSTX-U and DIII-D and scaling predicted by the theory and simulations. Here, the DIII-D simulation results also demonstrate that GAEs can be excited for larger aspect-ratio devices and smaller (sub- Alfvénic) beam injection velocities, therefore they can be unstable in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Materials Plasma Exposure eXperiment: Status of the Physics Basis Together With the Conceptual Design and Plans Forward

The Materials Plasma Exposure eXperiment (MPEX) is a linear plasma device that will address the plasma–material interaction (PMI) science for future fusion reactors and will enable testing of plasma-facing components (PFCs). It is designed as a steady-state device eventually delivering an ion fluence of up to 10 31 m -2 to the target. The device will be designed to handle neutron-activated materials. These capabilities should push the technical readiness level of PFCs up to six for some end-of-life aspects. In order to achieve the relevant plasma conditions, as they are expected in future fusion reactor divertors, MPEX will utilize a novel plasma source concept. This plasma source concept uses a high-power helicon (200 kW, 13.56 MHz), an electron cyclotron heating (ECH) system, which will heat electrons via electron Bernstein wave (EBW) heating (about 400 kW, 70 GHz), and an ion cyclotron heating (ICH) system (400 kW, 6–9 MHz). The physics basis for this plasma source concept including the heating physics and transport is based on experiments on Proto-MPEX. An overview of the experimental results and the physics basis will be given. The physics basis directly translates into functional requirements of MPEX and the conceptual design. The status of the conceptual design of MPEX will be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Low-energy H- beam injectors for particle accelerators: upgrade requirements, challenges and future plans [Slides]

High brightness, negative hydrogen ion sources are used extensively in large, accelerator-based, user facilities operating worldwide. Negative hydrogen beams have become the preferred means of filling circular accelerators and storage rings as well as enabling efficient extraction from cyclotrons. Much larger beams of H- ions are also utilized for neutral beams injection into fusion machines and will be discussed in a separate companion presentation. Some of the accelerator facilities include the US Spallation Neutron Source (SNS), Japan Proton Accelerator Research Complex (J-PARC), Rutherford Appleton Laboratory (RAL-ISIS), Los Alamos Neutron Science Center (LANSCE), Fermi National Accelerator Laboratory (FNAL), Brookhaven National Laboratory (BNL), the CERN LHC injector, the Chinese Spallation Neutron Source (CSNS) as well as numerous installations of D-Pace (licenced by TRIUMF) ion sources used mainly with cyclotrons. Many facilities are currently working on various improvement / upgrade projects, both in the near and long term, which are driving further development of their ion source and LEBT (Low Energy Beam Transport) and, in some cases, their overall front end injector system. This report will first provide a simple description of each facilities existing ion source and LEBT and summarize the operational parameters which are currently being routinely injected into their accelerators. Next, the parametric goals of each of the facilities upgrade projects are specified as well as development efforts currently underway to meet these requirement and improvement goals. It is hoped that that this work will not only capture the current state-of-the-art of worldwide H- beam injectors and clarify the research goals and efforts of the community in general but will also encourage further inter laboratory collaborations.

43 PARTICLE ACCELERATORS↗

Drift-cyclotron loss-cone instability in 3-D simulations of a sloshing-ion simple mirror

The kinetic stability of collisionless, sloshing beam-ion (45° pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over 1 – 10 μs, drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool ∼1keV ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗