Ales Necas

The challenge of developing a new energy source based on fusion power has led to dramatic growth in the private fusion industry over the past several years, with over 50 companies worldwide now pur...

A verification study of a semi-implicit energy-conserving electrostatic particle-in-cell algorithm is presented. The algorithm relaxes the time-step and mesh-size constraints that require resolution of the plasma period and Debye length associated with traditional explicit momentum-conserving particle-in-cell algorithms. Physical implications and applicability of using the semi-implicit scheme for modeling high-density plasmas are discussed. Where possible, numerical results are compared against analytical solutions. The simulation results indicate that the algorithm is stable at time steps larger than twice the inverse plasma frequency and cell sizes larger than the Debye length. It is found that the algorithm gives adequate results, provided that the distribution function and the spatiotemporal scales dictating the physics of the problem are resolved. As such, the algorithm may provide a robust method for kinetic modeling of high-density plasmas at scale.

Proton-boron-11 (p- ¹¹ B) fusion is an attractive concept for advanced fusion energy because of its nearly aneutronic characteristics. However, efficient use of beam-injected fast ions may be limited by their slowing down through Coulomb collisions, which lead to lower fusion reactivity. In this study we investigate whether ion cyclotron range of frequencies (ICRF) heating can ameliorate this effect by resonantly accelerating fast ions in the Large Helical Device (LHD) of a stellarator. Using the tangential neutral beam (NB) injection to supply energetic protons at ∼160-170 keV, together with boron seeding, we conducted a series of experiments with superimposed hydrogen minority ICRF heating for beam proton acceleration. Three complementary experimental methods, ICRF power modulation, NB termination, and NB energy scan, consistently showed that ICRF waves accelerate beam ions and enhance the p- ¹¹ B fusion yield. This led to measurable enhancements in the p- ¹¹ B fusion alpha-particle yield, as observed with a dedicated fusion-product detector. These findings provide the first direct experimental evidence that ICRF heating can improve the efficiency of beam-driven p- ¹¹ B fusion reactions.

We report on THz pulses with energy flux exceeding 2.5 J/sr and high directionality, generated from picosecond, joule laser irradiation of microwires. The impact of wire diameter, length, and laser irradiation orientation is explored with the irradiation of 160 μm diameter, 750 μm, sharpened wires at the dull base found to produce the highest energy and shortest wavelength THz pulses among the targets tested. A traveling wave antenna combined with coherent transition radiation model is developed to explain the emission pattern, pulse energy, and spectrum seen.

A one-year experimental program was performed using the Multi-TeraWatt (MTW) laser at the Laboratory for Laser Energetics (LLE) to investigate the efficiency and physical process of laser-driven—wire THz generation. This process was evaluated as a potential high-repetition-rate source of powerful THz pulses for use in THz pulsed polarimetry of a moderate-density magnetic-confinement fusion concept developed by TAE. The goal was to develop a source for short pulses in the THz regime (150 μm wavelength) containing 10’s of mJ. The campaign tested the impact of Cu wire diameter (25, 75, 160 μm), length (250 – 750 μm), and laser irradiation orientation on THz production efficiency in comparison to 20-μm thick Cu foils. Of these targets, the 160 µm diameter, 750 µm long wires with sharpened tips and irradiated at the base were found to produce the highest energy and shortest wavelength THz pulses, exceeding 2.5 Joule/sr and meeting the requirements for the polarimetry diagnostic. The brightest THz signal was radiated in the direction of the sharpened wire tip. A combination of two models – a traveling wave antenna combined with coherent transition radiation – describes the emission pattern, pulse energy, and spectrum observed in the experiments.

Plasma flow and acceleration in a magnetic mirror configuration are studied using a drift-kinetic particles-in-cell model in the paraxial approximation, with an emphasis on finite temperature effects and energy transport. Energy conversion between electrons and ions, overall energy balance, and axial energy losses are investigated. The simulations of plasma flow, acceleration, and energy transport in the magnetic mirror are extended into the high-density regimes with implicit particle-in-cell simulations. It is shown that profiles of the anisotropic ion temperatures and heat fluxes obtained with the full drift-kinetic model compare favorably with the results of a fluid model, which includes collisionless ion heat fluxes beyond the two-pressure adiabatic equations. The effects of collisions on trapped electrons and the resulting impacts on electron temperature and electric field profiles are investigated using a model collision operator.

The field reversed configuration (FRC), such as studied in the C-2W experiment at TAE Technologies, is an attractive candidate for realizing a nuclear fusion reactor. In an FRC, kinetic ion effects play the majority role in macroscopic stability, which allows global stability studies to make use of fluid-kinetic hybrid (also referred to as Ohm's law) models wherein ions are treated kinetically while electrons are treated as a fluid. The development and validation of such a hybrid particle-in-cell algorithm in the Exascale Computing Project code WarpX are reported here. Implementation of this model in the WarpX framework benefits from the numerical efficiency of WarpX as well as its scalability on large HPC systems and portability to different architectures. Performance benchmarks of the new algorithm for large, 3-dimensional, full device simulations from the Perlmutter supercomputer are presented. Results of a series of FRC simulations are discussed in which the impact of two-fluid effects on the tilt-mode growth rate was studied. It was observed that, in agreement with previous Hall-MHD studies, two-fluid effects have a stabilizing impact on the tilt mode.