Source author record

Amitava Bhattacharjee

Amitava Bhattacharjee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

27 published item(s)

preprint2022arXiv

Connection between quasi-symmetric magnetic fields and anisotropic pressure equilibria in fusion plasmas

The stellarator as a concept of magnetic confinement fusion requires careful design to confine particles effectively. A design possibility is to equip the magnetic field with a property known as quasisymmetry. Though it is generally believed that a steady-state quasisymmetric equilibrium can only be exact locally (unless the system has a direction of continuous symmetry such as the tokamak), we suggest in this work that a change in the equilibrium paradigm can ameliorate this limitation. We demonstrate that there exists a deep physical connection between quasisymmetry and magnetostatic equilibria with anisotropic pressure, extending beyond the isotropic pressure equilibria commonly considered.

preprint2022arXiv

Phases and phase-transitions in quasisymmetric configuration space

We explore the structure of the space of quasisymmetric configurations identifying them by their magnetic axes, described as 3D closed curves. We demonstrate that this topological perspective divides the space of all configurations into well-separated quasisymmetric phases. Each phase is characterized by the self-linking number (a topological invariant), defining different symmetry configurations (quasi-axisymmetry or quasi-helical symmetry). The phase-transition manifolds correspond to quasi-isodynamic configurations. By considering some models for closed curves (most notably torus unknots), general features associated with these phases are explored. Some general criteria are also built and leveraged to provide a simple way to describe existing quasisymmetric designs. This constitutes the first step in a program to identify quasisymmetric configurations with a reduced set of functions and parameters, to deepen understanding of configuration space, and offer an alternative approach to stellarator optimization that begins with the magnetic axis and builds outward.

preprint2022arXiv

Research Opportunities in Plasma Astrophysics

Major scientific questions and research opportunities are described on 10 unprioritized plasma astrophysics topics: (1) magnetic reconnection, (2) collisionless shocks and particle acceleration, (3) waves and turbulence, (4) magnetic dynamos, (5) interface and shear instabilities, (6) angular momentum transport, (7) dusty plasmas, (8) radiative hydrodynamics, (9) relativistic, pair-dominated and strongly magnetized plasmas, (10) jets and outflows. Note that this is a conference report from a Workshop on Opportunities in Plasma Astrophysics (WOPA, https://w3.pppl.gov/conferences/2010/WOPA/) in January 2010, that attracted broad representation from the community and was supported by the U.S. Department of Energy, National Aeronautics and Space Administration, National Science Foundation, American Physical Society's Topical Group for Plasma Astrophysics and Division of Plasma Physics, and Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. Although there has been much planning and many developments in both science and infrastructure since the report was written, most of the motivation, priorities, problems and technical challenges discussed therein remain unaddressed and are relevant at the time of posting.

preprint2022arXiv

Weakly Quasisymmetric Near-Axis Solutions to all Orders

We show that the equations satisfied by weakly quasisymmetric magnetic fields can be solved to arbitrarily high order in powers of the distance from the magnetic axis. This demonstration does not consider force balance. The existence of solutions requires an appropriate choice of parameters, most notably the toroidal current or rotational transform profiles. We do not prove that the expansion converges (it is likely divergent but asymptotic), and thus the demonstration here should not be taken as definitive proof of the existence of global solutions. Instead, we provide a systematic construction of solutions to arbitrarily high order.

preprint2021arXiv

Adjoint methods for quasisymmetry of vacuum fields on a surface

Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared to finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasisymmetry and a rotational transform value on a surface. To measure quasisymmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.

preprint2021arXiv

Measures of quasisymmetry for stellarators

Quasisymmetric stellarators are an attractive class of optimised magnetic confinement configurations. The property of quasisymmetry (QS) is in practice limited to be approximate, and thus the construction requires measures that quantify the deviation from the exact property. In this paper we study three measure candidates used in the literature, placing the focus on their origin and the comparison of their forms. The analysis shows clearly the lack of universality in these measures. As these metrics do not directly correspond to any physical property (except when exactly QS), optimisation should employ additional physical metrics for guidance. Especially close to QS minima, treating QS metrics through inequality constraints so that additional physics metrics dominate optimisation is suggested. The impact of different quasisymmetric measures on optimisation is presented through an example.

preprint2021arXiv

Small-Scale Dynamo in Stably Stratified Turbulence

We present numerical investigations into three principal properties of the small-scale dynamo in stably stratified turbulence: the onset criterion, the growth rate, and the nature of the magnetic field anisotropy in the kinematic regime. The results suggest that all three dynamo properties are controlled by the scale separation between the Ozmidov scale and the viscous or resistive scale. In addition to the critical magnetic Reynolds number, this allows for the definition of critical buoyancy and magnetic buoyancy Reynolds numbers for stratified small-scale dynamo onset in the high and low magnetic Prandtl number regimes, respectively. The presence of a small-scale dynamo in stellar radiative zones could affect dynamics through resulting Maxwell stresses and/or influence large-scale dynamo mechanisms in regions of differential rotation. Taking the solar radiative zone as a representative example and applying the onset criterion, we find that the stratification is strong enough to make the small-scale dynamo marginally active in the stably stratified turbulence of the solar tachocline.

preprint2020arXiv

Necessary and Sufficient Conditions for Quasisymmetry

A necessary and sufficient set of conditions for a quasisymmetric magnetic field in the form of constraint equations is derived from first principles. Without any assumption regarding the magnetohydrodynamic (MHD) equilibrium of the plasma, conditions for quasisymmetry are constructed starting from the single-particle Lagrangian to leading order. The conditions presented in the paper are less restrictive than the set recently obtained by Burby et al. [arXiv:1912.06468], and could facilitate ongoing efforts towards investigating the existence of global quasisymmetric MHD equilibria. It is also shown that quasisymmetry implies the existence of flux surfaces, regardless of whether the field corresponds to an MHD equilibrium.

preprint2020arXiv

Solving the problem of overdetermination of quasisymmetric equilbrium solutions by near-axis expansions: I. Generalised force balance

It is well known that the process of construction of quasisymmetric magnetic fields in magnetostatic equilibrium with isotropic pressure suffers from the problem of overdetermination. This has led to the widespread belief that global quasisymmetric solutions are likely not to exist. We develop a general near-axis expansion procedure that does not rely on the assumption of magnetostatic equilibria with isotropic pressure. We then demonstrate that in equilibria with anisotropic pressure, it is possible to circumvent the problem of overdetermination and carry out the power-series solutions to higher order. This suggests, contrary to current belief, that the existence of globally quasisymmetric fields is likely if one relaxes the assumption of magnetostatic equilibria with isotropic pressure.

preprint2020arXiv

Solving the problem of overdetermination of quasisymmetric equilbrium solutions by near-axis expansions. II. Circular axis stellarators

We apply the near-axis expansion method for quasisymmetric magnetic fields with anisotropic pressure (developed in the companion paper, Part I [E. Rodriguez, A. Bhattacharjee, arXiv:2008.04715 (2020)]) to construct numerical solutions to circular axis stellarators. The solutions are found to second order in the distance from the axis, not possible in the standard Garren-Boozer construction [D. A. Garren and A. H. Boozer, Physics of Fluids B: Plasma Physics, 2822 (1991)], which assumes magnetostatic equilibria with isotropic pressure. In the limit of zero anisotropy, it is shown that a subset of coefficients can be chosen to avoid the overdetermination problem.

preprint2020arXiv

Structure and overstability of resistive modes with runaway electrons

We investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfvén velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C1 linear simulation, and good agreement is found in both cases.

preprint2016arXiv

Adjoint Fokker-Planck equation and runaway electron dynamics

The adjoint Fokker-Planck equation method is applied to study the runaway probability function and the expected slowing-down time for highly relativistic runaway electrons, including the loss of energy due to synchrotron radiation. In direct correspondence to Monte Carlo simulation methods, the runaway probability function has a smooth transition across the runaway separatrix, which can be attributed to effect of the pitch angle scattering term in the kinetic equation. However, for the same numerical accuracy, the adjoint method is more efficient than the Monte Carlo method. The expected slowing-down time gives a novel method to estimate the runaway current decay time in experiments. A new result from this work is that the decay rate of high energyelectrons is very slow when E is close to the critical electric field. This effect contributes further to a hysteresis previously found in the runaway electron population.

preprint2016arXiv

Fluid moments of the Landau collision operator

One important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multi-dimensional Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The proposed methodology provides a systematic and rigorous means of extending the validity of fluid models that have an underlying inverse-square force particle dynamics to weakly collisional and strong flow regimes.

preprint2016arXiv

Hall current effects in mean-field dynamo theory

The role of the Hall term on large scale dynamo action is investigated by means of the First Order Smoothing Approximation. It is shown that the standard $α$ coefficient is altered, and is zero when a specific double Beltrami state is attained, in contrast to the Alfvénic state for MHD dynamos. The $β$ coefficient is no longer positive definite, and thereby enables dynamo action even if $α$-quenching were to operate. The similarities and differences with the (magnetic) shear-current effect are pointed out, and a mechanism that may be potentially responsible for $β< 0$ is advanced. The results are compared against previous studies, and their astrophysical relevance is also highlighted.

preprint2016arXiv

On the Value of the Reconnection Rate

Numerical simulations have consistently shown that the reconnection rate in certain collisionless regimes can be fast, on the order of $0.1 v_A B_u$, where $v_A$ and $B_u$ are the Alfvén speed and the reconnecting magnetic field upstream of the ion diffusion region. This particular value has been reported in myriad numerical simulations under disparate conditions. However, despite decades of research, the reasons underpinning this specific value remain mysterious. Here, we present an overview of this problem and discuss the conditions under which the "0.1 value" is attained. Furthermore, we explain why this problem should be interpreted in terms of the ion diffusion region length.

preprint2016arXiv

The magnetic shear-current effect: generation of large-scale magnetic fields by the small-scale dynamo

A novel large-scale dynamo mechanism, the magnetic shear-current effect, is discussed and explored. The effect relies on the interaction of magnetic fluctuations with a mean shear flow, meaning the saturated state of the small-scale dynamo can drive a large-scale dynamo -- in some sense the inverse of dynamo quenching. The dynamo is nonhelical, with the mean-field $α$ coefficient zero, and is caused by the interaction between an off-diagonal component of the turbulent resistivity and the stretching of the large-scale field by shear flow. Following up on previous numerical and analytic work, this paper presents further details of the numerical evidence for the effect, as well as an heuristic description of how magnetic fluctuations can interact with shear flow to produce the required electromotive force. The pressure response of the fluid is fundamental to this mechanism, which helps explain why the magnetic effect is stronger than its kinematic cousin, and the basic idea is related to the well-known lack of turbulent resistivity quenching by magnetic fluctuations. As well as being interesting for its applications to general high Reynolds number astrophysical turbulence, where strong small-scale magnetic fluctuations are expected to be prevalent, the magnetic shear-current effect is a likely candidate for large-scale dynamo in the unstratified regions of ionized accretion disks. Evidence for this is discussed, as well as future research directions and the challenges involved with understanding details of the effect in astrophysically relevant regimes.

preprint2015arXiv

Coherent nonhelical shear dynamos driven by magnetic fluctuations at low Reynolds numbers

Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. The primary results -- serving as a follow up to the results of Squire & Bhattacharjee [arXiv:1506.04109 (2015)] -- pertain to the "magnetic shear-current effect" as a viable mechanism to drive large-scale magnetic field generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, demonstrating the importance of rotation in determining the relative dominance of each mechanism.

preprint2015arXiv

Electromotive force due to magnetohydrodynamic fluctuations in sheared rotating turbulence

This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action -- the magnetic analogue of the "shear-current" effect. In addition, consideration of $α$ effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of $α$, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.

preprint2015arXiv

Generation of large-scale magnetic fields by small-scale dynamo in shear flows

We propose a new mechanism for turbulent mean-field dynamo in which the magnetic fluctuations resulting from a small-scale dynamo drive the generation of large-scale magnetic fields. This is in stark contrast to the common idea that small-scale magnetic fields should be harmful to large-scale dynamo action. These dynamos occur in the presence of large-scale velocity shear and do not require net helicity, resulting from off-diagonal components of the turbulent resistivity tensor as the magnetic analogue of the "shear-current" effect. Given the inevitable existence of non-helical small-scale magnetic fields in turbulent plasmas, as well as the generic nature of velocity shear, the suggested mechanism may help to explain generation of large-scale magnetic fields across a wide range of astrophysical objects.

preprint2015arXiv

Iris si iv line profiles: An indication for the plasmoid instability during small-scale magnetic reconnection on the sun

Our understanding of the process of fast reconnection has undergone a dramatic change in the last 10 years driven, in part, by the availability of high-resolution numerical simulations that have consistently demonstrated the break-up of current sheets into magnetic islands, with reconnection rates that become independent of Lundquist number, challenging the belief that fast magnetic reconnection in flares proceeds via the Petschek mechanism that invokes pairs of slow-mode shocks connected to a compact diffusion region. The reconnection sites are too small to be resolved with images but these reconnection mechanisms, Petschek and the plasmoid instability, have reconnection sites with very different density and velocity structures and so can be distinguished by high-resolution line-profiles observations. Using IRIS spectroscopic observations we obtain a survey of typical line profiles produced by small-scale events thought to be reconnection sites on the Sun. Slit-jaw images are used to investigate the plasma heating and re-configuration at the sites. A sample of 15 events from two active regions is presented. The line profiles are complex with bright cores and broad wings extending to over 300 km/s. The profiles can be reproduced with the multiple magnetic islands and acceleration sites that characterise the plasmoid instability but not by bi-directional jets that characterise the Petschek mechanism. This result suggests that if these small-scale events are reconnection sites, then fast reconnection proceeds via the plasmoid instability, rather than the Petschek mechanism during small-scale reconnection on the Sun.

preprint2015arXiv

The island coalescence problem: scaling of reconnection in extended fluid models including higher-order moments

As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of today's simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell (PIC) simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here we perform the complementary resistive MHD, Hall MHD and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.

preprint2015arXiv

The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support

Recent increases in supercomputing power, driven by the multi-core revolution and accelerators such as the IBM Cell processor, graphics processing units (GPUs) and Intel's Many Integrated Core (MIC) technology have enabled kinetic simulations of plasmas at unprecedented resolutions, but changing HPC architectures also come with challenges for writing efficient numerical codes. This paper describes the Plasma Simulation Code (PSC), an explicit, electromagnetic particle-in-cell code with support for different order particle shape functions. We focus on two distinguishing feature of the code: patch-based load balancing using space-filling curves, and support for Nvidia GPUs, which achieves substantial speed-up of up to more than 6x on the Cray XK7 architecture compared to a CPU-only implementation.

preprint2015arXiv

Variational formulation of relaxed and multi-region relaxed magnetohydrodynamics

Ideal magnetohydrodynamics (IMHD) is strongly constrained by an infinite number of microscopic constraints expressing mass, entropy and magnetic flux conservation in each infinitesimal fluid element, the latter preventing magnetic reconnection. By contrast, in the Taylor relaxation model for formation of macroscopically self-organized plasma equilibrium states, all these constraints are relaxed save for global magnetic fluxes and helicity. A Lagrangian variational principle is presented that leads to a new, fully dynamical, \emph{relaxed magnetohydrodynamics} (RxMHD), such that all static solutions are Taylor states but also allows flow. By postulating that some long-lived macroscopic current sheets can act as barriers to relaxation, separating the plasma into multiple relaxation regions, a further generalization, \emph{multi-region relaxed magnetohydrodynamics} (MRxMHD) is developed.

preprint2014arXiv

Magnetorotational instability: nonmodal growth and the relationship of global modes to the shearing box

We study the magnetorotational instability (MRI) using nonmodal stability techniques. Despite the spectral instability of many forms of the MRI, this proves to be a natural method of analysis that is well-suited to deal with the non-self-adjoint nature of the linear MRI equations. We find that the fastest growing linear MRI structures on both local and global domains can look very different to the eigenmodes, invariably resembling waves shearing with the background flow (shear waves). In addition, such structures can grow many times faster than the least stable eigenmode over long time periods, and be localized in a completely different region of space. These ideas lead -- for both axisymmetric and non-axisymmetric modes -- to a natural connection between the global MRI and the local shearing box approximation. By illustrating that the fastest growing global structure is well described by the ordinary differential equations (ODEs) governing a single shear wave, we find that the shearing box is a very sensible approximation for the linear MRI, contrary to many previous claims. Since the shear wave ODEs are most naturally understood using nonmodal analysis techniques, we conclude by analyzing local MRI growth over finite time-scales using these methods. The strong growth over a wide range of wave-numbers suggests that nonmodal linear physics could be of fundamental importance in MRI turbulence.

preprint2014arXiv

Nonmodal growth of the magnetorotational instability

We analyze the linear growth of the magnetorotational instability (MRI) in the short time limit using nonmodal methods. Our findings are quite different to standard results, illustrating that shearing wave energy can grow at the maximum MRI rate, $-dΩ/d \ln r,$ for any choice of azimuthal and vertical wavelengths. In addition, by comparing the growth of shearing waves with static structures, we show that over short time-scales shearing waves will always be dynamically more important than static structures in the ideal limit. By demonstrating that fast linear growth is possible at all wavelengths, these results suggest that nonmodal linear physics could play a fundamental role in MRI turbulence.

preprint2014arXiv

Observations of supra-arcade fans: instabilities at the head of reconnection jets

Supra-arcade fans are bright, irregular regions of emission that develop during eruptive flares, above flare arcades. The underlying flare arcades are thought to be a consequence of magnetic reconnection along a current sheet in the corona. At the same time, theory predicts plasma jets from the reconnection site which would be extremely difficult to observe directly because of their low density. It has been suggested that the dark supra-arcade downflows (SADs) seen falling through supra-arcade fans may be low density jet plasma. The head of a low density jet directed towards higher density plasma would be Rayleigh-Taylor unstable, and lead to the development of rapidly growing low and high density fingers along the interface. Using SDO/AIA 131A images, we show details of SADs seen from three different orientations with respect to the flare arcade and current sheet, and highlight features that have been previously unexplained, such as the splitting of SADs at their heads, but are a natural consequence of instabilities above the arcade. Comparison with 3-D magnetohydrodynamic simulations suggests that supra-arcade downflows are the result of secondary instabilities of the Rayleigh-Taylor type in the exhaust of reconnection jets.

preprint2014arXiv

Statistical simulation of the magnetorotational dynamo

Turbulence and dynamo induced by the magnetorotational instability (MRI) are analyzed using quasi-linear statistical simulation methods. It is found that homogenous turbulence is unstable to a large scale dynamo instability, which saturates to an inhomogenous equilibrium with a strong dependence on the magnetic Prandtl number (Pm). Despite its enormously reduced nonlinearity, the dependence of the angular momentum transport on Pm in the quasi-linear model is qualitatively similar to that of nonlinear MRI turbulence. This indicates that recent convergence problems may be related to large scale dynamo and suggests how dramatically simplified models may be used to gain insight into the astrophysically relevant regimes of very low or high Pm.