No Arabic abstract
The favourable properties of tungsten borides for shielding the central High Temperature Superconductor (HTS) core of a spherical tokamak fusion power plant are modelled using the MCNP code. The objectives are to minimize the power deposition into the cooled HTS core, and to keep HTS radiation damage to acceptable levels by limiting the neutron and gamma fluxes. The shield materials compared are W2B, WB, W2B5 and WB4 along with a reactively sintered boride B0.329C0.074Cr0.024Fe0.274W0.299, monolithic W and WC. Of all these W2B5 gave the most favourable results with a factor of ~10 or greater reduction in neutron flux and gamma energy deposition as compared to monolithic W. These results are compared with layered water-cooled shields, giving the result that the monolithic shields, with moderating boron, gave comparable neutron flux and power deposition, and (in the case of W2B5) even better performance. Good performance without water-coolant has advantages from a reactor safety perspective due to the risks associated with radio-activation of oxygen. 10B isotope concentrations between 0 and 100% are considered for the boride shields. The naturally occurring 20% fraction gave much lower energy depositions than the 0% fraction, but the improvement largely saturated beyond 40%. Thermophysical properties of the candidate materials are discussed, in particular the thermal strain. To our knowledge, the performance of W2B5 is unrivalled by other monolithic shielding materials. This is partly as its trigonal crystal structure gives it higher atomic density compared with other borides. It is also suggested that its high performance depends on it having just high enough 10B content to maintain a constant neutron energy spectrum across the shield.
Here, we propose a two-dimensional tungsten boride (WB4) lattice, with the Gibbs free energy for the adsorption of atomic hydrogen, tending to be the ideal value of 0 eV at 3% strained state, to host a better hydrogen evolution reaction activity. Based on first-principles calculations, it is demonstrated that the multiple d-p-pi and d-p-sigma Dirac conjugations of WB4 lattice ensures its excellent electronic transport characteristics. Meanwhile, coupling with the d-orbitals of W, the p-orbitals of borophene subunits in WB4 lattice can modulate the d band center to get a good HER performance. Our results not only provide a versatile platform for hosting multiple Dirac semimetal states with a sandwich configuration, but also offer a guiding principle for discovering the relationship between intrinsic properties of the active centre and the catalytic activity of metal layer from the emerging field of low-dimensional noble-metal-free lattices.
Turbulence is a major factor limiting the achievement of better tokamak performance as it enhances the transport of particles, momentum and heat which hinders the foremost objective of tokamaks. Hence, understanding and possibly being able to control turbulence in tokamaks is of paramount importance, not to mention our intellectual curiosity of it.
The 3D quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam-plasma interactions. In this algorithm, the plasma response to a non-evolving laser or particle beam is calculated using Maxwells equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. At the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r-z and a gridless description in $phi$. We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations are derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code, a full 3D quasi-static code, and the new quasi-static PIC code with azimuthal decomposition are also presented.
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the Mega-Ampere Spherical Tokamak (MAST), is studied as a high Lundquist number and low plasma-beta magnetic reconnection experiment. In this paper, 2D fluid simulations are presented of this merging process in order to understand the underlying physics, and better interpret the experimental data. These simulations examine the individual and combined effects of tight-aspect ratio geometry and two-fluid physics on the merging. The ideal self-driven flux-rope dynamics are coupled to the diffusion layer physics, resulting in a large range of phenomena. For resistive MHD simulations, the flux-ropes enter the sloshing regime for normalised resistivity eta < 1E-5. In Hall-MHD three regimes are found for the qualitative behaviour of the current sheet, depending on the ratio of the current sheet width to the ion-sound radius. These are a stable collisional regime, an open X-point regime, and an intermediate regime that is highly unstable to tearing-type instabilities. In toroidal axisymmetric geometry, the final state after merging is a MAST-like spherical tokamak with nested flux-surfaces. It is also shown that the evolution of simulated 1D radial density profiles closely resembles the Thomson scattering electron density measurements in MAST. An intuitive explanation for the origin of the measured density structures is proposed, based upon the results of the toroidal Hall-MHD simulations.
A modular, maintainable and extensible particle beam simulation architecture is presented. Design considerations for single particle, multi particle, and rms envelope simulations (in two and three dimensions) are outlined. Envelope simulation results have been validated against Trace3D. Hybridization with a physics-centric contol-system abstraction provides a convenient environment for rapid deployment of applications employing model-reference control strategies.