Researcher profile

Justin Ball

Justin Ball contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
1topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2020arXiv

Eliminating turbulent self-interaction through the parallel boundary condition in local gyrokinetic simulations

In this work, we highlight an issue that may reduce the accuracy of many local nonlinear gyrokinetic simulations - turbulent self-interaction through the parallel boundary condition. Given a sufficiently long parallel correlation length, individual turbulent eddies can span the full domain and "bite their own tails," thereby altering their statistical properties. Such self-interaction is only modeled accurately when the simulation domain corresponds to a full flux surface, otherwise it is artificially strong. For Cyclone Base Case parameters and typical domain sizes, we find that this mechanism modifies the heat flux by roughly 40% and it can be even more important. The effect is largest when using kinetic electrons, low magnetic shear, and strong turbulence drive (i.e. steep background gradients). It is found that parallel self-interaction can be eliminated by increasing the parallel length and/or the binormal width of the simulation domain until convergence is achieved.

preprint2016arXiv

Effect of the Shafranov shift and the gradient of $β$ on intrinsic momentum transport in up-down asymmetric tokamaks

Tokamaks with up-down asymmetric poloidal cross-sections spontaneously rotate due to turbulent transport of momentum. In this work, we investigate the effect of the Shafranov shift on this intrinsic rotation, primarily by analyzing tokamaks with tilted elliptical flux surfaces. By expanding the Grad-Shafranov equation in the large aspect ratio limit we calculate the magnitude and direction of the Shafranov shift in tilted elliptical tokamaks. The results show that, while the Shafranov shift becomes up-down asymmetric and depends strongly on the tilt angle of the flux surfaces, it is insensitive to the shape of the current and pressure profiles (when the geometry, total plasma current, and average pressure gradient are kept fixed). Next, local nonlinear gyrokinetic simulations of these MHD equilibria are performed with GS2, which reveal that the Shafranov shift can significantly enhance the momentum transport. However, to be consistent, the effect of $β'$ (i.e. the radial gradient of $β$) on the magnetic equilibrium was also included, which was found to significantly reduce momentum transport. Including these two competing effects broadens the rotation profile, but leaves the on-axis value of the rotation roughly unchanged. Consequently, the shape of the $β$ profile has a significant effect on the rotation profile of an up-down asymmetric tokamak.

preprint2016arXiv

Poloidal tilting symmetry of high order tokamak flux surface shaping in gyrokinetics

A poloidal tilting symmetry of the local nonlinear $δf$ gyrokinetic model is demonstrated analytically and verified numerically. This symmetry shows that poloidally rotating all the flux surface shaping effects with large poloidal mode number by a single tilt angle has an exponentially small effect on the transport properties of a tokamak. This is shown using a generalization of the Miller local equilibrium model to specify an arbitrary flux surface geometry. With this geometry specification we find that, when performing an expansion in large flux surface shaping mode number, the governing equations of gyrokinetics are symmetric in the poloidal tilt of the high order shaping effects. This allows us to take the fluxes from a single configuration and calculate the fluxes in any configuration that can be produced by tilting the large mode number shaping effects. This creates a distinction between tokamaks with mirror symmetric flux surfaces and tokamaks without mirror symmetry, which is expected to have important consequences for generating toroidal rotation using up-down asymmetry.

preprint2016arXiv

Scaling of up-down asymmetric turbulent momentum flux with poloidal shaping mode number in tokamaks

Breaking the up-down symmetry of tokamaks removes a constraint limiting intrinsic momentum transport, and hence toroidal rotation, to be small. Using gyrokinetic theory, we study the effect of different up-down asymmetric flux surface shapes on the turbulent transport of momentum. This is done by perturbatively expanding the gyrokinetic equation in large flux surface shaping mode number. It is found that the momentum flux generated by shaping that lacks mirror symmetry (which is necessarily up-down asymmetric) has a power law scaling with the shaping mode number. However, the momentum flux generated by mirror symmetric flux surface shaping (even if it is up-down asymmetric) decays exponentially with large shaping mode number. These scalings are consistent with nonlinear local gyrokinetic simulations and indicate that low mode number shaping effects (e.g. elongation, triangularity) are optimal for creating rotation. Additionally it suggests that breaking the mirror symmetry of flux surfaces may generate significantly more toroidal rotation.

preprint2016arXiv

Up-down asymmetric tokamaks

Bulk toroidal rotation has proven capable of stabilising both dangerous MHD modes and turbulence. In this thesis, we explore a method to drive rotation in large tokamaks: up-down asymmetry in the magnetic equilibrium. We seek to maximise this rotation by finding optimal up-down asymmetric flux surface shapes. First, we use the ideal MHD model to show that low order external shaping (e.g. elongation) is best for creating up-down asymmetric flux surfaces throughout the device. Then, we calculate realistic up-down asymmetric equilibria for input into nonlinear gyrokinetic turbulence analysis. Analytic gyrokinetics shows that, in the limit of fast shaping effects, a poloidal tilt of the flux surface shaping has little effect on turbulent transport. Since up-down symmetric surfaces do not transport momentum, this invariance to tilt implies that devices with mirror symmetry about any line in the poloidal plane will drive minimal rotation. Accordingly, further analytic investigation suggests that non-mirror symmetric flux surfaces with envelopes created by the beating of fast shaping effects may create significantly stronger momentum transport. Guided by these analytic results, we carry out local nonlinear gyrokinetic simulations of non-mirror symmetric flux surfaces created with the lowest possible shaping effects. First, we consider tilted elliptical flux surfaces with a Shafranov shift and find little increase in the momentum transport when the effect of the pressure profile on the equilibrium is included. We then simulate flux surfaces with independently-tilted elongation and triangularity. These two-mode configurations show a $60\%$ increase over configurations with just elongation or triangularity. A rough analytic estimate indicates that the optimal two-mode configuration can drive rotation with an on-axis Alfven Mach number of $1.5 \%$ in an ITER-like machine.

preprint2015arXiv

Radial penetration of flux surface shaping in tokamaks

Using analytic calculations, the effects of the edge flux surface shape and the toroidal current profile on the penetration of flux surface shaping are investigated in a tokamak. It is shown that the penetration of shaping is determined by the poloidal variation of the poloidal magnetic field on the surface. This fact is used to investigate how different flux surface shapes penetrate from the edge. Then, a technique to separate the effects of magnetic pressure and tension in the Grad-Shafranov equation is presented and used to calculate radial profiles of strong elongation for nearly constant current profiles. Lastly, it is shown that more hollow toroidal current profiles are significantly better at conveying shaping from the edge to the core.

preprint2014arXiv

Intrinsic momentum transport in up-down asymmetric tokamaks

Recent work demonstrated that breaking the up-down symmetry of tokamak flux surfaces removes a constraint that limits intrinsic momentum transport, and hence toroidal rotation, to be small. We show, through MHD analysis, that ellipticity is most effective at introducing up-down asymmetry throughout the plasma. We detail an extension to GS2, a local $δf$ gyrokinetic code that self-consistently calculates momentum transport, to permit up-down asymmetric configurations. Tokamaks with tilted elliptical poloidal cross-sections were simulated to determine nonlinear momentum transport. The results, which are consistent with experiment in magnitude, suggest that a toroidal velocity gradient, $(\partial u_{ζi} / \partial ρ) / v_{th i}$, of 5% of the temperature gradient, $(\partial T_{i} / \partial ρ) / T_{i}$, is sustainable. Here $v_{th i}$ is the ion thermal speed, $u_{ζi}$ is the ion toroidal mean flow, $ρ$ is the minor radial coordinate normalized to the tokamak minor radius, and $T_{i}$ is the ion temperature. Since other intrinsic momentum transport mechanisms scale poorly to larger machines, these results indicate that up-down asymmetry is the most feasible method to generate the current experimentally-measured rotation levels in reactor-sized devices.

preprint2013arXiv

Conditions for up-down asymmetry in the core of tokamak equilibria

A local magnetic equilibrium solution is sought around the magnetic axis in order to identify the key parameters defining the magnetic-surface's up-down asymmetry in the core of tokamak plasmas. The asymmetry is found to be determined essentially by the ratio of the toroidal current density flowing on axis to the fraction of the external field's odd perturbation that manages to propagate from the plasma boundary into the core. The predictions are tested and illustrated first with an analytical Solovev equilibrium and then using experimentally relevant numerical equilibria. Hollow current-density distributions, and hence reverse magnetic shear, are seen to be crucial to bring into the core asymmetry values that are usually found only near the plasma edge.