Researcher profile

A. Thyagaraja

A. Thyagaraja contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph physics.comp-ph physics.plasm-ph

Trust snapshot

Quick read

Trust 15 - Baseline

3works

0followers

3topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2012arXiv

CENTORI: a global toroidal electromagnetic two-fluid plasma turbulence code

A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of electron and ion fluid equations and Maxwell equations in toroidal configurations with axisymmetric equilibria. Uniquely, the equilibrium is co-evolved with the turbulence, and is thus modified by it. CENTORI is applicable to tokamaks of arbitrary aspect ratio and high plasma beta. A predictor-corrector, semi-implicit finite difference scheme is used to compute the time evolution of fluid quantities and fields. Vector operations and the evaluation of flux surface averages are speeded up by choosing the Jacobian of the transformation from laboratory to plasma coordinates to be a function of the equilibrium poloidal magnetic flux. A subroutine, GRASS, is used to co-evolve the plasma equilibrium by computing the steady-state solutions of a diffusion equation with a pseudo-time derivative. The code is written in Fortran 95 and is efficiently parallelized using Message Passing Interface (MPI). Illustrative examples of output from simulations of a tearing mode in a large aspect ratio tokamak plasma and of turbulence in an elongated conventional aspect ratio tokamak plasma are provided.

preprint2003arXiv

Field-guided proton acceleration at reconnecting X-points in flares

An explicitly energy-conserving full orbit code CUEBIT, developed originally to describe energetic particle effects in laboratory fusion experiments, has been applied to the problem of proton acceleration in solar flares. The model fields are obtained from solutions of the linearised MHD equations for reconnecting modes at an X-type neutral point, with the additional ingredient of a longitudinal magnetic field component. To accelerate protons to the highest observed energies on flare timescales, it is necessary to invoke anomalous resistivity in the MHD solution. It is shown that the addition of a longitudinal field component greatly increases the efficiency of ion acceleration, essentially because it greatly reduces the magnitude of drift motions away from the vicinity of the X-point, where the accelerating component of the electric field is largest. Using plasma parameters consistent with flare observations, we obtain proton distributions extending up to gamma-ray-emitting energies (>1MeV). In some cases the energy distributions exhibit a bump-on-tail in the MeV range. In general, the shape of the distribution is sensitive to the model parameters.

preprint2000arXiv

Azimuthally symmetric MHD and two-fluid equilibria with arbitrary flows

Magnetohydrodynamic (MHD) and two-fluid quasi-neutral equilibria with azimuthal symmetry, gravity and arbitrary ratios of (nonrelativistic) flow speed to acoustic and Alfven speeds are investigated. In the two-fluid case, the mass ratio of the two species is arbitrary, and the analysis is therefore applicable to electron-positron plasmas. The methods of derivation can be extended in an obvious manner to several charged species. Generalized Grad-Shafranov equations, describing the equilibrium magnetic field, are derived. Flux function equations and Bernoulli relations for each species, together with Poisson's equation for the gravitational potential, complete the set of equations required to determine the equilibrium. These are straightforward to solve numerically. The two-fluid system, unlike the MHD system, is shown to be free of singularities. It is demonstrated analytically that there exists a class of incompressible MHD equilibria with magnetic field-aligned flow. A special sub--class first identified by S. Chandrasekhar, in which the flow speed is everywhere equal to the local Alfven speed, is compatible with virtually any azimuthally symmetric magnetic configuration. Potential applications of this analysis include extragalactic and stellar jets, and accretion disks.