Researcher profile

James Juno

James Juno contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

7 published item(s)

preprint2026arXiv

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species $f_s$ on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where $f_s$ is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6-60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

preprint2025arXiv

The effect of shear flow on the resistive tearing instability

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-$Ψ$ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-$Ψ$ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

preprint2022arXiv

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E$\times$B) flow relative to the unmagnetized ions can cause the Electron Cyclotron Drift Instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in collisionless shock ramps in space, and in $\rm{E \times B}$ discharge devices such as Hall thrusters. ECDI can induce an electron flow parallel to the background E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport may lead to particle thermalization at space shocks, and may cause unfavorable plasma flows towards the walls of E$\times$B devices. The development of ECDI and anomalous transport is often considered fully-kinetic. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a two-fluid framework without any collision. By treating electrons and ions on an equal footing, the free energy allows the growth of an anomalous electron flow parallel to the background E field. We first present linear analyses of the instability in the two-fluid 5- and 10-moment models, and compare them against the fully-kinetic theory. At lower temperatures, the two-fluid fastest-growing mode is in good agreement with the kinetic result. Also, by including more ($>=10$) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on various parameters is also explored. We then carry out direct numerical simulations of the cross-field setup using the 5-moment model. The growth of the instability and the anomalous transport is confirmed. Finally, 5-moment and Vlasov simulations using identical parameters in the lower-temperature regime are performed, showing reasonable agreement.

preprint2022arXiv

Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species interactions are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, satisfies the H-Theorem unlike analogous Bhatnagar-Gross-Krook operators.

preprint2020arXiv

A Deep Dive into the Distribution Function: Understanding Phase Space Dynamics with Continuum Vlasov-Maxwell Simulations

In collisionless and weakly collisional plasmas, the particle distribution function is a rich tapestry of the underlying physics. However, actually leveraging the particle distribution function to understand the dynamics of a weakly collisional plasma is challenging. The equation system of relevance, the Vlasov-Maxwell-Fokker-Planck (VM-FP) system of equations, is difficult to numerically integrate, and traditional methods such as the particle-in-cell method introduce counting noise into the distribution function. In this thesis, we present a new algorithm for the discretization of VM-FP system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin (DG) finite element method for the spatial discretization and a third order strong-stability preserving Runge-Kutta for the time discretization, we obtain an accurate solution for the plasma's distribution function in space and time. We both prove the numerical method retains key physical properties of the VM-FP system, such as the conservation of energy and the second law of thermodynamics, and demonstrate these properties numerically. These results are contextualized in the history of the DG method. We discuss the importance of the algorithm being alias-free, a necessary condition for deriving stable DG schemes of kinetic equations so as to retain the implicit conservation relations embedded in the particle distribution function, and the computational favorable implementation using a modal, orthonormal basis in comparison to traditional DG methods applied in computational fluid dynamics. Finally, we demonstrate how the high fidelity representation of the distribution function, combined with novel diagnostics, permits detailed analysis of the energization mechanisms in fundamental plasma processes such as collisionless shocks.

preprint2020arXiv

Alias-free, matrix-free, and quadrature-free discontinuous Galerkin algorithms for (plasma) kinetic equations

Understanding fundamental kinetic processes is important for many problems, from plasma physics to gas dynamics. A first-principles approach to these problems requires a statistical description via the Boltzmann equation, coupled to appropriate field equations. In this paper we present a novel version of the discontinuous Galerkin (DG) algorithm to solve such kinetic equations. Unlike Monte-Carlo methods we use a continuum scheme in which we directly discretize the 6D phase-space using discontinuous basis functions. Our DG scheme eliminates counting noise and aliasing errors that would otherwise contaminate the delicate field-particle interactions. We use modal basis functions with reduced degrees of freedom to improve efficiency while retaining a high formal order of convergence. Our implementation incorporates a number of software innovations: use of JIT compiled top-level language, automatically generated computational kernels and a sophisticated shared-memory MPI implementation to handle velocity space parallelization.

preprint2020arXiv

Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields

We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.