BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

A. Bhattacharjee

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

ResearcherUnclaimed source record
physics.plasm-phastro-ph.SRphysics.space-phphysics.comp-phphysics.flu-dynastro-ph.EPastro-ph.HEastro-ph.IMastro-ph.GAphysics.geo-ph

Catalog footprint

What is connected

41works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

41 published item(s)

preprint2022arXiv

Laser-Driven, Ion-Scale Magnetospheres in Laboratory Plasmas. I. Experimental Platform and First Results

Magnetospheres are a ubiquitous feature of magnetized bodies embedded in a plasma flow. While large planetary magnetospheres have been studied for decades by spacecraft, ion-scale "mini" magnetospheres can provide a unique environment to study kinetic-scale, collisionless plasma physics in the laboratory to help validate models of larger systems. In this work, we present preliminary experiments of ion-scale magnetospheres performed on a unique high-repetition-rate platform developed for the Large Plasma Device (LAPD) at UCLA. The experiments utilize a high-repetition-rate laser to drive a fast plasma flow into a pulsed dipole magnetic field embedded in a uniform magnetized background plasma. 2D maps of magnetic field with high spatial and temporal resolution are measured with magnetic flux probes to examine the evolution of magnetosphere and current density structures for a range of dipole and upstream parameters. The results are further compared to 2D PIC simulations to identify key observational signatures of the kinetic-scale structures and dynamics of the laser-driven plasma. We find that distinct 2D kinetic-scale magnetopause and diamagnetic current structures are formed at higher dipole moments, and their locations are consistent with predictions based on pressure balances and energy conservation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Radio Dichotomy in Quasars with H$β$ FWHM greater than $15,000$ km\,s$^{-1}$

It has been inferred from large unbiased samples that $10\%$-$15\%$ of all quasars are radio-loud (RL). Using the quasar catalog from the Sloan Digital Sky Survey, we show that the radio-loud fraction (RLF) for high broad line (HBL) quasars, containing H$β$ FWHM greater than $15,000$ km s$^{-1}$, is $\sim 57 \%$. While there is no significant difference between the RL and radio-quiet (RQ) populations in our sample in terms of their black hole mass, Eddington ratio, and covering fraction (CF), optical continuum luminosity of the RL quasars are higher. The similarity in the distribution of their CF indicates that our analysis is unbiased in terms of the viewing angle of the HBL RL and RQ quasars. Hence, we conclude that the accretion disc luminosity of the RL quasars in our HBL sample is higher, which indicates a connection between a brighter disc and a more prominent jet. By comparing them with the non-HBL H$β$ broad emission line quasars, we find that the HBL sources have the lowest Eddington ratios in addition to having a very high RLF. That is consistent with the theories of jet formation, in which jets are launched from low Eddington ratio accreting systems. We find that the [O III] narrow emission line is stronger in the RL compared to RQ quasars in our HBL sample, which is consistent with previous findings in the literature, and may be caused by the interaction of the narrow line gas with the jet.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Weak Alfvénic turbulence in relativistic plasmas. Part 1. Dynamical equations and basic dynamics of interacting resonant triads

Alfvén wave collisions are the primary building blocks of the non-relativistic turbulence that permeates the heliosphere and low-to-moderate energy astrophysical systems. However, many astrophysical systems such as gamma-ray bursts, pulsar and magnetar magnetospheres, and active galactic nuclei have relativistic flows or energy densities. To better understand these high energy systems, we derive reduced relativistic MHD equations and employ them to examine weak Alfvénic turbulence, dominated by three-wave interactions, in reduced relativistic magnetohydrodynamics, including the force-free, infinitely magnetized limit. We compare both numerical and analytical solutions to demonstrate that many of the findings from non-relativistic weak turbulence are retained in the relativistic system. But, an important distinction in the relativistic limit is the inapplicability of a formally incompressible limit, i.e, there exists finite coupling to the compressible fast mode regardless of the strength of the magnetic field. Since fast modes can propagate across field lines, this mechanism provides a route for energy to escape strongly magnetized systems, e.g., magnetar magnetospheres. However, we find that the fast-Alfvén coupling is diminished in the limit of oblique propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast magnetic reconnection in highly-extended current sheets at the National Ignition Facility

Fast magnetic reconnection was observed between magnetized laser-produced plasmas at the National Ignition Facility. Two highly-elongated plasma plumes were produced by tiling two rows of lasers, with magnetic fields generated in each plume by the Biermann battery effect. Detailed magnetic field observations, obtained from proton radiography using a D$^3$He capsule implosion, reveal reconnection occurring in an extended, quasi-1D current sheet with large aspect ratio $\sim 100$. The 1-D geometry allowed a rigorous and unique reconstruction of the magnetic field, which showed a reconnection current sheet that thinned down to a half-width close to the electron gyro-scale. Despite the large aspect ratio, a large fraction of the magnetic flux reconnected, suggesting fast reconnection supported by the non-gyrotropic electron pressure tensor.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinetic simulations of electron pre-energization by magnetized collisionless shocks in expanding laboratory plasmas

Collisionless shocks are common features in space and astrophysical systems where supersonic plasma flows interact, such as in the solar wind, the heliopause, and supernova remnants. Recent experimental capabilities and diagnostics allow detailed laboratory investigations of high-Mach-number shocks, which therefore can become a valuable way to understand shock dynamics in various astrophysical environments. Using 2D particle-in-cell simulations with a Coulomb binary collision operator, we demonstrate the mechanism for generation of energetic electrons and experimental requirements for detecting this process in the laboratory high-Mach-number collisionless shocks. We show through a parameter study that electron acceleration by magnetized collisionless shocks is feasible in laboratory experiments with laser-driven expanding plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Major Scientific Challenges and Opportunities in Understanding Magnetic Reconnection and Related Explosive Phenomena in Solar and Heliospheric Plasmas

Magnetic reconnection underlies many explosive phenomena in the heliosphere and in laboratory plasmas. The new research capabilities in theory/simulations, observations, and laboratory experiments provide the opportunity to solve the grand scientific challenges summarized in this whitepaper. Success will require enhanced and sustained investments from relevant funding agencies, increased interagency/international partnerships, and close collaborations of the solar, heliospheric, and laboratory plasma communities. These investments will deliver transformative progress in understanding magnetic reconnection and related explosive phenomena including space weather events.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An extended MHD study of the 16 October 2015 MMS diffusion region crossing

The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A heuristic model for MRI turbulent stresses in Hall MHD

Although the Shakura-Sunyaev $α$ viscosity prescription has been highly successful in characterizing myriad astrophysical environments, it has proven to be partly inadequate in modelling turbulent stresses driven by the MRI. Hence, we adopt the approach employed by \citet{GIO03}, but in the context of Hall magnetohydrodynamics (MHD), to study MRI turbulence. We utilize the exact evolution equations for the stresses, and the non-linear terms are closed through the invocation of dimensional analysis and physical considerations. We demonstrate that the inclusion of the Hall term leads to non-trivial results, including the modification of the Reynolds and Maxwell stresses, as well as the (asymptotic) non-equipartition between the kinetic and magnetic energies; the latter issue is also addressed via the analysis of non-linear waves. The asymptotic ratio of the kinetic and magnetic energies is shown to be \emph{independent} of the choice of initial conditions, but it is governed by the \emph{Hall parameter}. We contrast our model with the Kazantsev prescription from small-scale dynamo theory, and the Hall term does not contribute in the latter approach, illustrating the limitations of the Kazantsev formalism. We indicate potential astrophysical applications of our model, including the solar wind where a lack of equipartition has been observed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of current singularity in a topologically constrained plasma

Recently a variational integrator for ideal magnetohydrodynamics in Lagrangian labeling has been developed. Its built-in frozen-in equation makes it optimal for studying current sheet formation. We use this scheme to study the Hahm-Kulsrud-Taylor problem, which considers the response of a 2D plasma magnetized by a sheared field under sinusoidal boundary forcing. We obtain an equilibrium solution that preserves the magnetic topology of the initial field exactly, with a fluid mapping that is non-differentiable. Unlike previous studies that examine the current density output, we identify a singular current sheet from the fluid mapping. These results are benchmarked with a constrained Grad-Shafranov solver. The same signature of current singularity can be found in other cases with more complex magnetic topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Magnetic reconnection in plasma under inertial confinement fusion conditions driven by heat flux effects in Ohm's law

In the interaction of high-power laser beams with solid density plasma there are a number of mechanisms that generate strong magnetic fields. Such fields subsequently inhibit or redirect electron flows, but can themselves be advected by heat fluxes, resulting in complex interplay between thermal transport and magnetic fields.We show that for heating by multiple laser spots reconnection of magnetic field lines can occur, mediated by these heat fluxes, using a fully implicit 2D Vlasov-Fokker-Planck code. Under such conditions, the reconnection rate is dictated by heat flows rather than Alfvènic flows. We find that this mechanism is only relevant in a high $β$ plasma. However, the Hall parameter $ω_c τ_{ei}$ can be large so that thermal transport is strongly modified by these magnetic fields, which can impact longer time scale temperature homogeneity and ion dynamics in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Turbulent Magnetohydrodynamic Reconnection Mediated by the Plasmoid Instability

It has been established that the Sweet-Parker current layer in high Lundquist number reconnection is unstable to the super-Alfvénic plasmoid instability. Past two-dimensional magnetohydrodynamic simulations have demonstrated that the plasmoid instability leads to a new regime where the Sweet-Parker current layer changes into a chain of plasmoids connected by secondary current sheets, and the averaged reconnection rate becomes nearly independent of the Lundquist number. In this work, three-dimensional simulation with a guide field shows that the additional degree of freedom allows plasmoid instabilities to grow at oblique angles, which interact and lead to self-generated turbulent reconnection. The averaged reconnection rate in the self-generated turbulent state is of the order of a hundredth of the characteristic Alfvén speed, which is similar to the two-dimensional result but is an order of magnitude lower than the fastest reconnection rate reported in recent studies of externally driven three-dimensional turbulent reconnection. Kinematic and magnetic energy fluctuations both form elongated eddies along the direction of local magnetic field, which is a signature of anisotropic magnetohydrodynamic turbulence. Both energy fluctuations satisfy power-law spectra in the inertial range, where the magnetic energy spectral index is in the range from $-2.3$ to $-2.1$, while the kinetic energy spectral index is slightly steeper, in the range from $-2.5$ to $-2.3$. The anisotropy of turbulence eddies is found to be nearly scale-independent, in contrast with the prediction of the Goldreich-Sridhar theory for anisotropic turbulence in a homogeneous plasma permeated by a uniform magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comparison of multi-fluid moment models with Particle-in-Cell simulations of collisionless magnetic reconnection

We introduce an extensible multi-fluid moment model in the context of collisionless magnetic reconnection. This model evolves full Maxwell equations, and simultaneously moments of the Vlasov-Maxwell equation for each species in the plasma. Effects like electron inertia and pressure gradient are self-consistently embedded in the resulting multi-fluid moment equations, without the need to explicitly solving a generalized Ohms's law. Two limits of the multi-fluid moment model are discussed, namely, the five-moment limit that evolves a scalar pressures for each species, and the ten-moment limit that evolves the full anisotropic, non-gyrotropic pressure tensor for each species. We first demonstrate, analytically and numerically, that the five-moment model reduces to the widely used Hall Magnetohydrodynamics (Hall MHD) model under the assumptions of vanishing electron inertia, infinite speed of light, and quasi-neutrality. Then, we compare ten-moment and fully kinetic Particle-In-Cell (PIC) simulations of a large scale Harris sheet reconnection problem, where the ten-moment equations are closed with a local linear collisionless approximation for the heat flux. The ten-moment simulation gives reasonable agreement with the PIC results regarding the structures and magnitudes of the electron flows, the polarities and magnitudes of elements of the electron pressure tensor, and the decomposition of the generalized Ohm's law. Possible ways to improve the simple local closure towards a nonlocal fully three-dimensional closure are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Helicity-Flux-Driven Alpha Effect in Laboratory and Astrophysical Plasmas

The constraint imposed by magnetic helicity conservation on the alpha effect is considered for both magnetically and flow dominated self-organizing plasmas. Direct numerical simulations are presented for a dominant contribution to the alpha effect, which can be cast in the functional form of a total divergence of an averaged helicity flux, called the helicity-flux-driven alpha ( H$α$) effect. Direct numerical simulations of the H$α$ effect are prese nted for two examples---the magnetically dominated toroidal plasma unstable to tearing modes, and the flow-dominated accretion disk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rapid Change of Field Line Connectivity and Reconnection in Stochastic Magnetic Fields

Magnetic fields without a direction of continuous symmetry have the generic feature that neighboring field lines exponentiate away from each other and become stochastic, hence the ideal constraint of preserving magnetic field line connectivity becomes exponentially sensitive to small deviations from ideal Ohm's law. The idea of breaking field line connectivity by stochasticity as a mechanism for fast reconnection is tested with numerical simulations based on reduced magnetohydrodynamics equations with a strong guide field line-tied to two perfectly conducting end plates. Starting from an ideally stable force-free equilibrium, the system is allowed to undergo resistive relaxation. Two distinct phases are found in the process of resistive relaxation. During the quasi-static phase, rapid change of field line connectivity and strong induced flow are found in regions of high field line exponentiation. However, although the field line connectivity of individual field lines can change rapidly, the overall pattern of field line mapping appears to deform gradually. From this perspective, field line exponentiation appears to cause enhanced diffusion rather than reconnection. In some cases, resistive quasi-static evolution can cause the ideally stable initial equilibrium to cross a stability threshold, leading to formation of intense current filaments and rapid change of field line mapping into a qualitatively different pattern. It is in this onset phase that the change of field line connectivity is more appropriately designated as magnetic reconnection. Our results show that rapid change of field line connectivity appears to be a necessary, but not a sufficient condition for fast reconnection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Secondary Rayleigh-Taylor type Instabilities in the Reconnection Exhaust Jet as a Mechanism for Supra-Arcade Downflows

Supra-arcade downflows (hereafter referred to as SADs) are low-emission, elongated, finger-like features usually observed in active-region coronae above post-eruption flare arcades. Observations exhibit downward moving SADs intertwined with bright upward moving spikes. Whereas SADs are dark voids, spikes are brighter, denser structures. Although SADs have been observed for decades, the mechanism of formation of SADs remains an open issue. In our three-dimensional resistive magnetohydrodynamic simulations, we demonstrate that secondary Rayleigh-Taylor type instabilities develop in the downstream region of a reconnecting current sheet. The instability results in the formation of low-density coherent structures that resemble SADs, and high-density structures that appear to be spike-like. Comparison between the simulation results and observations suggests that secondary Rayleigh-Taylor type instabilities in the exhaust of reconnecting current sheets provide a plausible mechanism for observed SADs and spikes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational integration for ideal magnetohydrodynamics with built-in advection equations

Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Current Sheet Statistics in Three-Dimensional Simulations of Coronal Heating

In a recent numerical study [Ng et al., Astrophys. J. {\bf 747}, 109, 2012], with a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of \al time in order to obtain converged statistical values. The heating rate found in these simulations saturate to a level that is independent of the Lundquist number. This scaling result was also supported by an analysis with the assumption of the Sweet-Parker scaling of the current sheets, as well as how the width, length and number of current sheets scale with Lundquist number. In order to test these assumptions, we have implemented an automated routine to analyze thousands of current sheets in these simulations and return statistical scalings for these quantities. It is found that the Sweet-Parker scaling is justified. However, some discrepancies are also found and require further study.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Filamentation instability of counter-streaming laser-driven plasmas

Filamentation due to the growth of a Weibel-type instability was observed in the interaction of a pair of counter-streaming, ablatively-driven plasma flows, in a supersonic, collisionless regime relevant to astrophysical collisionless shocks. The flows were created by irradiating a pair of opposing plastic (CH) foils with 1.8 kJ, 2-ns laser pulses on the OMEGA EP laser system. Ultrafast laser-driven proton radiography was used to image the Weibel-generated electromagnetic fields. The experimental observations are in good agreement with the analytical theory of the Weibel instability and with particle-in-cell simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Magnetic reconnection mediated by hyper-resistive plasmoid instability

Magnetic reconnection mediated by the hyper-resistive plasmoid instability is studied with both linear analysis and nonlinear simulations. The linear growth rate is found to scale as $S_{H}^{1/6}$ with respect to the hyper-resistive Lundquist number $S_{H}\equiv L^{3}V_{A}/η_{H}$, where $L$ is the system size, $V_{A}$ is the Alfvén velocity, and $η_{H}$ is the hyper-resistivity. In the nonlinear regime, reconnection rate becomes nearly independent of $S_{H}$, the number of plasmoids scales as $S_{H}^{1/2}$, and the secondary current sheet length and width both scale as $S_{H}^{-1/2}$. These scalings are consistent with a heuristic argument assuming secondary current sheets are close to marginal stability. The distribution of plasmoids as a function of the enclosed flux $ψ$ is found to obey a $ψ^{-1}$ power law over an extended range, followed by a rapid fall off for large plasmoids. These results are compared with those from resistive magnetohydrodynamic studies.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Plasmoid Instability in High-Lundquist-Number Magnetic Reconnection

Our understanding of magnetic reconnection in resistive magnetohydrodynamics has gone through a fundamental change in recent years. The conventional wisdom is that magnetic reconnection mediated by resistivity is slow in laminar high Lundquist ($S$) plasmas, constrained by the scaling of the reconnection rate predicted by Sweet-Parker theory. However, recent studies have shown that when $S$ exceeds a critical value $\sim10^{4}$, the Sweet-Parker current sheet is unstable to a super-Alfvénic plasmoid instability, with a linear growth rate that scales as $S^{1/4}$. In the fully developed statistical steady state of two-dimensional resistive magnetohydrodynamic simulations, the normalized average reconnection rate is approximately 0.01, nearly independent of $S$, and the distribution function $f(ψ)$ of plasmoid magnetic flux $ψ$ follows a power law $f(ψ)\simψ^{-1}$. When Hall effects are included, the plasmoid instability may trigger onset of Hall reconnection even when the conventional criterion for onset is not satisfied. The rich variety of possible reconnection dynamics is organized in the framework of a phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Plasmoid solutions of the Hahm--Kulsrud--Taylor equilibrium model

The Hahm--Kulsrud (HK) [T. S. Hahm and R. M. Kulsrud, Phys. Fluids {\bf 28}, 2412 (1985)] solutions for a magnetically sheared plasma slab driven by a resonant periodic boundary perturbation illustrate fully shielded (current sheet) and fully reconnected (magnetic island) responses. On the global scale, reconnection involves solving a magnetohydrodynamic (MHD) equilibrium problem. In systems with a continuous symmetry such MHD equilibria are typically found by solving the Grad--Shafranov equation, and in slab geometry the elliptic operator in this equation is the 2-D Laplacian. Thus, assuming appropriate pressure and poloidal current profiles, a conformal mapping method can be used to transform one solution into another with different boundary conditions, giving a continuous sequence of solutions in the form of partially reconnected magnetic islands (plasmoids) separated by Syrovatsky current sheets. The two HK solutions appear as special cases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-Lundquist Number Scaling in Three-Dimensional Simulations of Parker's Model of Coronal Heating

Parker's model is one of the most discussed mechanisms for coronal heating and has generated much debate. We have recently obtained new scaling results in a two-dimensional (2D) version of this problem suggesting that the heating rate becomes independent of resistivity in a statistical steady state [Ng and Bhattacharjee, Astrophys. J., 675, 899 (2008)]. Our numerical work has now been extended to 3D by means of large-scale numerical simulations. Random photospheric footpoint motion is applied for a time much longer than the correlation time of the motion to obtain converged average coronal heating rates. Simulations are done for different values of the Lundquist number to determine scaling. In the high-Lundquist number limit, the coronal heating rate obtained so far is consistent with a trend that is independent of the Lundquist number, as predicted by previous analysis as well as 2D simulations. In the same limit the average magnetic energy built up by the random footpoint motion tends to have a much weaker dependence on the Lundquist number than that in the 2D simulations, due to the formation of strong current layers and subsequent disruption when the equilibrium becomes unstable. We will present scaling analysis showing that when the dissipation time is comparable or larger than the correlation time of the random footpoint motion, the heating rate tends to become independent of Lundquist number, and that the magnetic energy production is also reduced significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Distribution of Plasmoids In High-Lundquist-Number Magnetic Reconnection

The distribution function $f(ψ)$ of magnetic flux $ψ$ in plasmoids formed in high-Lundquist-number current sheets is studied by means of an analytic phenomenological model and direct numerical simulations. The distribution function is shown to follow a power law $f(ψ)\simψ^{-1}$, which differs from other recent theoretical predictions. Physical explanations are given for the discrepant predictions of other theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Constrained Tectonics Model for Coronal Heating

An analytical and numerical treatment is given of a constrained version of the tectonics model developed by Priest, Heyvaerts, & Title [2002]. We begin with an initial uniform magnetic field ${\bf B} = B_0 \hat{\bf z}$ that is line-tied at the surfaces $z = 0$ and $z = L$. This initial configuration is twisted by photospheric footpoint motion that is assumed to depend on only one coordinate ($x$) transverse to the initial magnetic field. The geometric constraints imposed by our assumption precludes the occurrence of reconnection and secondary instabilities, but enables us to follow for long times the dissipation of energy due to the effects of resistivity and viscosity. In this limit, we demonstrate that when the coherence time of random photospheric footpoint motion is much smaller by several orders of magnitude compared with the resistive diffusion time, the heating due to Ohmic and viscous dissipation becomes independent of the resistivity of the plasma. Furthermore, we obtain scaling relations that suggest that even if reconnection and/or secondary instabilities were to limit the build-up of magnetic energy in such a model, the overall heating rate will still be independent of the resistivity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anisotropic MHD Turbulence

The solar wind and the interstellar medium are permeated by large-scale magnetic fields that render magnetohydrodynamic (MHD) turbulence anisotropic. In the weak-turbulence limit in which three-wave interactions dominate, analytical and high-resolution numerical results based on random scattering of shear-Alfvén waves propagating parallel to a large-scale magnetic field, as well as direct simulations demonstrate rigorously an anisotropic energy spectrum that scales as $k^{-2}_\perp$, instead of the famous Iroshnikov-Kraichnan (IK) spectrum of $k^{-3/2}$ for the isotropic case. Even in the absence of a background magnetic field, anisotropy is found to develop with respect to the local magnetic field, although the energy spectrum is globally isotropic and is found to be consistent with a $k^{-3/2}$ scaling. It is also found in direct numerical simulations that the energy cascade rate is much closer to IK scaling than a Kolmogorov scaling. Recent observations in the solar wind on cascade rates (as functions of the proton temperature and solar wind speed at 1 AU) seem to support this result [Vasquez et al. 2007].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast magnetic reconnection in laser-produced plasma bubbles

Recent experiments have observed magnetic reconnection in high-energy-density, laser-produced plasma bubbles, with reconnection rates observed to be much higher than can be explained by classical theory. Based on fully kinetic particle simulations we find that fast reconnection in these strongly driven systems can be explained by magnetic flux pile-up at the shoulder of the current sheet and subsequent fast reconnection via two-fluid, collisionless mechanisms. In the strong drive regime with two-fluid effects, we find that the ultimate reconnection time is insensitive to the nominal system Alfven time.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generalizations of a Nonlinear Fluid Model for Void Formation in Dusty Plasmas

New developments in the theory and numerical simulation of a recently proposed one-dimensional nonlinear time-dependent fluid model [K. Avinash, A. Bhattacharjee, and S. Hu, Phys. Rev. Lett. 90, 075001 (2003)] for void formation in dusty plasmas are presented. The model describes an initial instability caused by the ion drag, rapid nonlinear growth, and a nonlinear saturation mechanism that realizes a quasi-steady state containing a void. The earlier one-dimensional model has been extended to two and three dimensions (the latter, assuming spherical symmetry), using a more complete set of dynamical equations than was used in the earlier one-dimensional formulation. The present set of equations includes an ion continuity equation and a nonlinear ion drag operator. Qualitative features of void formation are shown to be robust with respect to different functional forms of the ion drag operator.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hall magnetohydrodynamic reconnection in the plasmoid unstable regime

A set of reduced Hall magnetohydrodynamic (MHD) equations are used to evaluate the stability of large aspect ratio current sheets to the formation of plasmoids (secondary islands). Reconnection is driven by resistivity in this analysis, which occurs at the resistive skin depth $d_η\equiv S_L^{-1/2} \sqrt{L v_A/γ}$, where $S_L$ is the Lundquist number, $L$ the length of the current sheet, $v_A$ the Alfvén speed, and $γ$ the growth rate. Modifications to a recent resistive MHD analysis [N.\ F.\ Loureiro, A.\ A.\ Schekochihin, and S.\ C. Cowley, Phys.\ Plasmas {\bf 14}, 100703 (2007)] arise when collisions are sufficiently weak that $d_η$ is shorter than the ion skin depth $d_i \equiv c/ω_{pi}$. Secondary islands grow faster in this Hall MHD regime: the maximum growth rate scales as $(d_i/L)^{6/13} S_L^{7/13} v_A/L$ and the number of plasmoids as $(d_i/L)^{1/13} S_L^{11/26}$, compared to $S_L^{1/4} v_A/L$ and $S^{3/8}$, respectively, in resistive MHD.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kolmogorov versus Iroshnikov-Kraichnan spectra: Consequences for ion heating in the solar wind

Whether the phenomenology governing MHD turbulence is Kolmogorov or Iroshnikov-Kraichnan (IK) remains an open question, theoretically as well as observationally. The ion heating profile observed in the solar wind provides a quantitative, if indirect, observational constraint on the relevant phenomenology. Recently, a solar wind heating model based on Kolmogorov spectral scaling has produced reasonably good agreement with observations, provided the effect of turbulence generation due to pickup ions is included in the model. Without including the pickup ion contributions, the Kolmogorov scaling predicts a proton temperature profile that decays too rapidly beyond a radial distance of 15 AU. In the present study, we alter the heating model by applying an energy cascade rate based on IK scaling, and show that the model yields higher proton temperatures, within the range of observations, with or without the inclusion of the effect due to pickup ions. Furthermore, the turbulence correlation length based on IK scaling seems to follow the trend of observations better.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large-Scale High-Lundquist Number Reduced MHD Simulations of the Solar Corona Using GPU Accelerated Machines

We have recently carried out a computational campaign to investigate a model of coronal heating in three-dimensions using reduced magnetohydrodynamics (RMHD). Our code is built on a conventional scheme using the pseudo-spectral method, and is parallelized using MPI. The current investigation requires very long time integrations using high Lundquist numbers, where the formation of very fine current layers challenge the resolutions achievable even on massively parallel machines. We present here results of a port to Nvidia CUDA (Compute Unified Device Architecture) for hardware acceleration using graphics processing units (GPUs). In addition to a brief discussion of our general strategy, we will report code performance on several machines which span a variety of hardware configurations and capabilities. These include a desktop workstation with commodity hardware, a dedicated research workstation equipped with four Nvidia C2050 GPUs, as well as several large-scale GPU accelerated distributed memory machines: Lincoln/NCSA, Dirac/NERSC, and Keeneland/NICS.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Onset of Fast Reconnection in Hall Magnetohydrodynamics Mediated by the Plasmoid Instability

The role of a super-Alfvénic plasmoid instability in the onset of fast reconnection is studied by means of the largest Hall magnetohydrodynamics simulations to date, with system sizes up to $10^{4}$ ion skin depths ($d_{i}$). It is demonstrated that the plasmoid instability can facilitate the onset of rapid Hall reconnection, in a regime where the onset would otherwise be inaccessible because the Sweet-Parker width is significantly above $d_{i}$. However, the topology of Hall reconnection is not inevitably a single stable X-point. There exists an intermediate regime where the single X-point topology itself exhibits instability, causing the system to alternate between a single X-point geometry and an extended current sheet with multiple X-points produced by the plasmoid instability. Through a series of simulations with various system sizes relative to $d_{i}$, it is shown that system size affects the accessibility of the intermediate regime. The larger the system size is, the easier it is to realize the intermediate regime. Although our Hall MHD model lacks many important physical effects included in fully kinetic models, the fact that a single X-point geometry is not inevitable raises the interesting possibility for the first time that Hall MHD simulations may have the potential to realize reconnection with geometrical features similar to those seen in fully kinetic simulations, namely, extended current sheets and plasmoid formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Reduced magnetohydrodynamic theory of oblique plasmoid instabilities

The three-dimensional nature of plasmoid instabilities is studied using the reduced magnetohydrodynamic equations. For a Harris equilibrium with guide field, represented by $\vc{B}_o = B_{po} \tanh (x/λ) \hat{y} + B_{zo} \hat{z}$, a spectrum of modes are unstable at multiple resonant surfaces in the current sheet, rather than just the null surface of the polodial field $B_{yo} (x) = B_{po} \tanh (x/λ)$, which is the only resonant surface in 2D or in the absence of a guide field. Here $B_{po}$ is the asymptotic value of the equilibrium poloidal field, $B_{zo}$ is the constant equilibrium guide field, and $λ$ is the current sheet width. Plasmoids on each resonant surface have a unique angle of obliquity $θ\equiv \arctan(k_z/k_y)$. The resonant surface location for angle $θ$ is $x_s = - λ\arctanh (\tan θB_{zo}/B_{po})$, and the existence of a resonant surface requires $|θ| < \arctan (B_{po} / B_{zo})$. The most unstable angle is oblique, i.e. $θ\neq 0$ and $x_s \neq 0$, in the constant-$ψ$ regime, but parallel, i.e. $θ= 0$ and $x_s = 0$, in the nonconstant-$ψ$ regime. For a fixed angle of obliquity, the most unstable wavenumber lies at the intersection of the constant-$ψ$ and nonconstant-$ψ$ regimes. The growth rate of this mode is $γ_{\textrm{max}}/Γ_o \simeq S_L^{1/4} (1-μ^4)^{1/2}$, in which $Γ_o = V_A/L$, $V_A$ is the Alfvén speed, $L$ is the current sheet length, and $S_L$ is the Lundquist number. The number of plasmoids scales as $N \sim S_L^{3/8} (1-μ^2)^{-1/4} (1 + μ^2)^{3/4}$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Effect of Magnetic Turbulence Energy Spectra and Pickup Ions on the Heating of the Solar Wind

In recent years, a phenomenological solar wind heating model based on a turbulent energy cascade prescribed by the Kolmogorov theory has produced reasonably good agreement with observations on proton temperatures out to distances around 70 AU, provided the effect of turbulence generation due to pickup ions is included in the model. In a recent study [Ng et al., J. Geophys. Res., 115, A02101 (2010)], we have incorporated in the heating model the energy cascade rate based on Iroshnikov-Kraichnan (IK) scaling. We showed that the IK cascade rate can also produce good agreement with observations, with or without the inclusion of pickup ions. This effect was confirmed both by integrating the model using average boundary conditions at 1 AU, and by applying a method [Smith et al., Astrophys. J., 638, 508 (2006)] that uses directly observed values as boundary conditions. The effects due to pickup ions is found to be less important for the IK spectrum, which is shallower than the Kolmogorov spectrum. In this paper, we will present calculations of the pickup ions effect in more details, and discuss the physical reason why a shallower spectrum generates less waves and turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Weakly collisional Landau damping and three-dimensional Bernstein-Greene-Kruskal modes: New results on old problems

Landau damping and Bernstein-Greene-Kruskal (BGK) modes are among the most fundamental concepts in plasma physics. While the former describes the surprising damping of linear plasma waves in a collisionless plasma, the latter describes exact undamped nonlinear solutions of the Vlasov equation. There does exist a relationship between the two: Landau damping can be described as the phase-mixing of undamped eigenmodes, the so-called Case-Van Kampen modes, which can be viewed as BGK modes in the linear limit. While these concepts have been around for a long time, unexpected new results are still being discovered. For Landau damping, we show that the textbook picture of phase-mixing is altered profoundly in the presence of collision. In particular, the continuous spectrum of Case-Van Kampen modes is eliminated and replaced by a discrete spectrum, even in the limit of zero collision. Furthermore, we show that these discrete eigenmodes form a complete set of solutions. Landau-damped solutions are then recovered as true eigenmodes (which they are not in the collisionless theory). For BGK modes, our interest is motivated by recent discoveries of electrostatic solitary waves in magnetospheric plasmas. While one-dimensional BGK theory is quite mature, there appear to be no exact three-dimensional solutions in the literature (except for the limiting case when the magnetic field is sufficiently strong so that one can apply the guiding-center approximation). We show, in fact, that two- and three-dimensional solutions that depend only on energy do not exist. However, if solutions depend on both energy and angular momentum, we can construct exact three-dimensional solutions for the unmagnetized case, and two-dimensional solutions for the case with a finite magnetic field. The latter are shown to be exact, fully electromagnetic solutions of the steady-state Vlasov-Poisson-Ampère system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effects of Line-tying on Magnetohydrodynamic Instabilities and Current Sheet Formation

An overview of some recent progress on magnetohydrodynamic stability and current sheet formation in a line-tied system is given. Key results on the linear stability of the ideal internal kink mode and resistive tearing mode are summarized. For nonlinear problems, a counterexample to the recent demonstration of current sheet formation by Low \emph{et al}. [B. C. Low and Å. M. Janse, Astrophys. J. \textbf{696}, 821 (2009)] is presented, and the governing equations for quasi-static evolution of a boundary driven, line-tied magnetic field are derived. Some open questions and possible strategies to resolve them are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Scaling laws of resistive magnetohydrodynamic reconnection in the high-Lundquist-number, plasmoid-unstable regime

The Sweet-Parker layer in a system that exceeds a critical value of the Lundquist number ($S$) is unstable to the plasmoid instability. In this paper, a numerical scaling study has been done with an island coalescing system driven by a low level of random noise. In the early stage, a primary Sweet-Parker layer forms between the two coalescing islands. The primary Sweet-Parker layer breaks into multiple plasmoids and even thinner current sheets through multiple levels of cascading if the Lundquist number is greater than a critical value $S_{c}\simeq4\times10^{4}$. As a result of the plasmoid instability, the system realizes a fast nonlinear reconnection rate that is nearly independent of $S$, and is only weakly dependent on the level of noise. The number of plasmoids in the linear regime is found to scales as $S^{3/8}$, as predicted by an earlier asymptotic analysis (Loureiro \emph{et al.}, Phys. Plasmas \textbf{14}, 100703 (2007)). In the nonlinear regime, the number of plasmoids follows a steeper scaling, and is proportional to $S$. The thickness and length of current sheets are found to scale as $S^{-1}$, and the local current densities of current sheets scale as $S^{-1}$. Heuristic arguments are given in support of theses scaling relations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theory of incompressible MHD turbulence with scale-dependent alignment and cross-helicity

(Abridged) An anisotropic theory of MHD turbulence with nonvanishing cross-helicity is constructed based on Boldyrev's alignment hypothesis and probabilities p and q for fluctuations v and b to be positively or negatively aligned. Guided by observations suggesting that the normalized cross-helicity and the probability p are approximately constant in the inertial range, a generalization of Boldyrev's theory is derived that depends on the three ratios w+/w-, epsilon+/epsilon-, and p/q. The theory reduces to Boldyrev's original theory when w+ = w-, epsilon+ = epsilon-, and p = q.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ultraviolet Coronagraph Spectroscopy: A Key Capability for Understanding the Physics of Solar Wind Acceleration

Understanding the physical processes responsible for accelerating the solar wind requires detailed measurements of the collisionless plasma in the extended solar corona. Some key clues about these processes have come from instruments that combine the power of an ultraviolet (UV) spectrometer with an occulted telescope. This combination enables measurements of ion emission lines far from the bright solar disk, where most of the solar wind acceleration occurs. Although the UVCS instrument on SOHO made several key discoveries, many questions remain unanswered because its capabilities were limited. This white paper summarizes these past achievements and also describes what can be accomplished with next-generation instrumentation of this kind.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Extension of the Electron Dissipation Region in Collisionless Hall MHD Reconnection

This paper presents Sweet-Parker type scaling arguments in the context of hyper-resistive Hall magnetohyrdodynamics (MHD). The predicted steady state scalings are consistent with those found by Chacón et al. [PRL 99, 235001 (2007)], though as with that study, no prediction of electron dissipation region \emph{length} is made. Numerical experiments confirm that both cusp-like and modestly more extended geometries are realizable. However, importantly, the length of the electron dissipation region, which is taken as a parameter by several recent studies, is found to depend explicitly on the level of hyper-resistivity. Furthermore, although hyper-resistivity can produce more extended electron dissipation regions, the length of the region remains smaller than one ion skin depth for the largest values of hyper-resistivity considered here. These electron dissipation regions are significantly shorter than those seen in many recent kinetic studies. The length of the electron dissipation region is found to depend on electron inertia as well, scaling like $(m_e/m_i)^{3/8}$. However, the thickness of the region appears to scale similarly, so that the aspect ratio is at most very weakly dependent on $(m_e/m_i)$. The limitations of scaling theories which do not predict the length of the electron dissipation region are emphasized.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fast Reconnection in High-Lundquist-Number Plasmas Due to Secondary Tearing Instabilities

Thin current sheets in systems of large size that exceed a critical value of the Lundquist number are unstable to a super-Alfvenic tearing instability. The scaling of the growth rate of the fastest growing instability with respect to the Lundquist number is shown to follow from the classical dispersion relation for tearing modes. As a result of this instability, the system realizes a nonlinear reconnection rate that appears to be weakly dependent on the Lundquist number, and larger than the Sweet-Parker rate by an order of magnitude (for the range of Lundquist numbers considered). This regime of fast reconnection appears to be realizable in a dynamic and highly unstable thin current sheet, without requiring the current sheet to be turbulent.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A comparison of spectral element and finite difference methods using statically refined nonconforming grids for the MHD island coalescence instability problem

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

0 citations0 reviewsNo code linkedOpen access