Search

This paper presents the Hammir tandem mirror design based on Realta Fusion's first-of-a-kind model for axisymmetric magnetic mirror fusion performance. This model uses an integrated end plug simulation model including, heating, equilibrium, and transport combined with a new formulation of the plasma operation contours (POPCONs) technique for the tandem mirror central cell. Using this model, it is shown that an end plug utilizing high temperature superconducting magnets and modern neutral beams enables a classical tandem mirror pilot plant producing a fusion gain Q > 5. The approach here represents an important advance in tandem mirror design. The high fidelity end plug model enables calculations of heating and transport in the highly non-Maxwellian end plug to be made more accurately and the central cell POPCON technique allows consideration of a wide range of parameters in the relatively simple near-Maxwellian central cell, facilitating the selection of more optimal central cell plasmas. These advances make it possible to find more conservative classical tandem mirror fusion pilot plant operating points with lower $\beta$, temperatures, neutral beam energies, and end plug performance than designs in the literature. Despite being more conservative, it is shown that these operating points can still form the basis of a viable fusion pilot plant.

We demonstrate through high-fidelity particle-in-cell simulations a simple approach for efficiently generating 20+ GeV electron beams with the necessary charge, energy spread, and emittance for use as the injector for an electron arm of a future linear collider or a next generation XFEL. The self-focusing of an unmatched, relatively low quality, drive beam results in self-injection by elongating the wakefield excited in the nonlinear blowout regime. Over pump depletion distances, the drive beam dynamics and self-loading from the injected beam leads to extremely high quality and high energy output beams. For plasma densities of $10^{18} \ \text{cm}^{-3}$, PIC simulation results indicate that self-injected beams with $0.52 \ \text{nC}$ of charge can be accelerated to $\sim 20$ GeV energies with projected energy spreads, $\lesssim 1\%$ within the beam core, slice normalized emittances as low as $110 \ \text{nm}$, a peak normalized brightness $\gtrsim 10^{19} \ \text{A}/\text{m}^2/\text{rad}^2$, and energy transfer efficiencies $\gtrsim 54\%$.

The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively.

This paper introduces a dynamic logic extension of separation logic. The assertion language of separation logic is extended with modalities for the five types of the basic instructions of separation logic: simple assignment, look-up, mutation, allocation, and de-allocation. The main novelty of the resulting dynamic logic is that it allows to combine different approaches to resolving these modalities. One such approach is based on the standard weakest precondition calculus of separation logic. The other approach introduced in this paper provides a novel alternative formalization in the proposed dynamic logic extension of separation logic. The soundness and completeness of this axiomatization has been formalized in the Coq theorem prover.

In a laser wakefield accelerator (LWFA), an intense laser pulse excites a plasma wave that traps and accelerates electrons to relativistic energies. When the pulse overlaps the accelerated electrons, it can enhance the energy gain through direct laser acceleration (DLA) by resonantly driving the betatron oscillations of the electrons in the plasma wave. The particle-in-cell (PIC) algorithm, although often the tool of choice to study DLA, contains inherent errors due to numerical dispersion and the time staggering of the electric and magnetic fields. Further, conventional PIC implementations cannot reliably disentangle the fields of the plasma wave and laser pulse, which obscures interpretation of the dominant acceleration mechanism. Here, a customized field solver that reduces errors from both numerical dispersion and time staggering is used in conjunction with a field decomposition into azimuthal modes to perform PIC simulations of DLA in an LWFA. Comparisons with traditional PIC methods, model equations, and experimental data show improved accuracy with the customized solver and convergence with an order-of-magnitude fewer cells. The azimuthal-mode decomposition reveals that the most energetic electrons receive comparable energy from DLA and LWFA., Comment: 10 pages, 5 figures, to submit to Physics of Plasmas

Accelerator-based X-ray free-electron lasers (XFELs) are the latest addition to the revolutionary tools of discovery for the 21st century. The two major components of an XFEL are an accelerator-produced electron beam and a magnetic undulator which tend to be kilometer-scale long and expensive. Here, we present an ultra-compact scheme to produce 10s of attosecond X-ray pulses with several GW peak power utilizing a novel aspect of the FEL instability using a highly chirped, pre-bunched and ultra-bright electron beam from a plasma-based accelerator interacting with an optical undulator. The self-selection of electrons from the combination of a highly chirped and pre-bunched beam leads to the stable generation of attosecond X-ray pulses. Furthermore, two-color attosecond pulses with sub-femtosecond separation can be produced by adjusting the energy distribution of the electron beam so that multiple FEL resonances occur at different locations within the beam. Such a tunable coherent attosecond X-ray sources may open up a new area of attosecond science enabled by X-ray attosecond pump/probe techniques., Comment: 8 pages, 4 figures

The overall problem addressed in this paper is the long-standing problem of program correctness, and in particular programs that describe systems of parallel executing processes. We propose a new method for proving correctness of parallel implementations of high-level transition system specifications. The implementation language underlying the method is based on the model of active (or concurrent) objects. The method defines correctness in terms of a simulation relation between the transition system which specifies the program semantics and the transition system that is described by the correctness specification. The simulation relation itself abstracts from the fine-grained interleaving of parallel processes by exploiting a global confluence property of the particular model of active objects considered in this paper. As a proof-of-concept we apply our method to the correctness of a parallel simulator of multicore memory systems.

Treatments of plasma waves usually assume homogeneity, but the parallel gradients ubiquitous in plasmas can modify wave propagation and absorption. We derive a quasilocal inhomogeneous correction to the plasma dielectric for arbitrary distributions by expanding the phase correlation integral and develop a novel integration technique that allows our correction to be applied in many situations and has greater accuracy than other inhomogeneous dielectric formulas found in the literature. We apply this dielectric tensor to the lower-hybrid current drive problem and demonstrate that inhomogeneous wave damping does not affect the lower-hybrid wave's linear damping condition, and in the non-Maxwellian problem damping and propagation should remain unchanged except in the case of waves with very large phase velocities., Comment: 24th Topical Conference on Radio-frequency Power in Plasmas

Raytracing/Fokker-Planck (FP) simulations used to model lower-hybrid current drive (LHCD) often fail to reproduce experimental results, particularly when LHCD is weakly damped. A proposed reason for this discrepancy is the lack of "full-wave" effects, such as diffraction and interference, in raytracing simulations and the breakdown of raytracing approximation. Previous studies of LHCD using non-Maxwellian full-wave/FP simulations have been performed, but these simulations were not self-consistent and enforced power conservation between the FP and full-wave code using a numerical rescaling factor. Here we have created a fully-self consistent full-wave/FP model for LHCD that is automatically power conserving. This was accomplished by coupling an overhauled version of the non-Maxwellian TORLH full-wave solver and the CQL3D FP code using the Integrated Plasma Simulator. We performed converged full-wave/FP simulations of Alcator C-Mod discharges and compared them to raytracing. We found that excellent agreement in the power deposition profiles from raytracing and TORLH could be obtained, however, TORLH had somewhat lower current drive efficiency and broader power deposition profiles in some cases. This discrepancy appears to be a result of numerical limitations present in the TORLH model and a small amount of diffractional broadening of the TORLH wave spectrum. Our results suggest full-wave simulation of LHCD is likely not necessary as diffraction and interference represented only a small correction that could not account for the differences between simulations and experiment.

A new ARC-class, highly-radiative, pulsed, L-mode, burning plasma scenario is developed and evaluated as a candidate for future tokamak reactors. Pulsed inductive operation alleviates the stringent current drive requirements of steady-state reactors, and operation in L-mode affords ELM-free access to $\sim90\%$ core radiation fractions, significantly reducing the divertor power handling requirements. In this configuration the fusion power density can be maximized despite L-mode confinement by utilizing high-field to increase plasma densities and current. This allows us to obtain high gain in robust scenarios in compact devices with $P_\mathrm{fus} > 1000\,$MW despite low confinement. We demonstrate the feasibility of such scenarios here; first by showing that they avoid violating 0-D tokamak limits, and then by performing self-consistent integrated simulations of flattop operation including neoclassical and turbulent transport, magnetic equilibrium, and RF current drive models. Finally we examine the potential effect of introducing negative triangularity with a 0-D model. Our results show high-field radiative pulsed L-mode scenarios are a promising alternative to the typical steady state advanced tokamak scenarios which have dominated tokamak reactor development.

Lower-hybrid current drive (LHCD) actuators are important components of modern day fusion experiments as well as proposed fusion reactors. However, simulations of LHCD often differ substantially from experimental results, and from each other, especially in the inferred power deposition profile shape. Here we investigate some possible causes of this discrepancy; "full-wave" effects such as interference and diffraction, which are omitted from standard raytracing simulations and the breakdown of the raytracing near reflections and caustics. We compare raytracing simulations to state-of-the-art full-wave simulations using matched hot-plasma dielectric tensors in realistic tokamak scenarios for the first time. We show that differences between full-wave simulations and raytracing in previous work were primarily due to numerical and physical inconsistencies in the simulations, and we demonstrate that good agreement between raytracing and converged full-wave simulations can be obtained in reactor relevant-scenarios with large ray caustics and in situations with weak damping.

High-fidelity modeling of plasma-based acceleration (PBA) requires the use of 3D fully nonlinear and kinetic descriptions based on the particle-in-cell (PIC) method. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the beam-plasma interaction. In a QSA-PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA-PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA-PIC algorithm that combines another speedup technique known as azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can significantly reduce the algorithmic complexity. Modeling the laser-plasma interaction in a full 3D PIC algorithm is very computationally expensive due to the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible with the QSA PIC algorithms based on the azimuthal mode expansion. This algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described.