BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

C. S. Ng

C. S. Ng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.plasm-phastro-ph.SRphysics.space-phphysics.comp-phastro-ph.EPastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
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

12 published item(s)

preprint2021arXiv

A Turbulent Heating Model Combining Diffusion and Advection Effects for Giant Planet Magnetospheres

The ion temperature of the magnetospheres of Jupiter and Saturn was observed to increase substantially from about 10 to 30 planet radii. Different heating mechanisms have been proposed to explain such observations, including a heating model for Jupiter based on MHD turbulence with flux-tube diffusion. More recently, an MHD turbulent heating model based on advection was shown to also explain the temperature increase at Jupiter and Saturn. We further develop this turbulent heating model by combining effects from both diffusion and advection. The combined model resolves the physical consistency requirement that diffusion should dominate over advection when the radial flow velocity is small and vice versa when it is large. Comparisons with observations show that previous agreements, using the advection only model, are still valid for larger radial distance. Moreover, the additional heating by diffusion results in a better agreement with the temperature observations for smaller radial distance.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Kinetic Flux Ropes: Bernstein-Greene-Kruskal Modes for the Vlasov-Poisson-Ampère System

Electrostatic structures have been observed in many regions of space plasmas, including the solar wind, the magnetosphere, the auroral acceleration region. One possible theoretical description of some of these structures is the concept of Bernstein-Greene-Kruskal (BGK) modes, which are exact nonlinear steady-state solutions of the Vlasov-Poisson system of equations in collisionless kinetic theory. We generalize exact solutions of two-dimensional BGK modes in a magnetized plasma with finite magnetic field strength [Ng, Bhattacharjee, and Skiff, Phys. Plasmas {\bf13}, 055903 (2006)] to cases with azimuthal magnetic fields so that these structures carry electric current as well as steady electric and magnetic fields. Such nonlinear solutions now satisfy exactly the Vlasov-Poisson-Ampère system of equations. Explicit examples with either positive or negative electric potential structure are provided.

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
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

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

High Lundquist Number Resistive MHD Simulations of Magnetic Reconnection: Searching for Secondary Island Formation

Recently, secondary island formation due to the tearing instability of the Sweet-Parker current sheet was identified as a possible mechanism that can lead to fast reconnection (less sensitive dependence on Lundquist number $S$) both in numerical simulations using Particle-in-Cell (PIC) method [Daughton et al. 2009], as well as using resistive magnetohydrodynamics (MHD) [Lapenta 2008; Bhattacharjee et al. 2009]. This instability is thought to appear when $S$ is greater than a certain threshold. These recent results prompt us to perform more resistive MHD simulations of a basic reconnection configuration based on the island coalescence instability, using much higher resolutions and larger $S$. Our simulations are based on a fairly standard pseudo spectral code, which has been tested for accuracy, convergence, and compared well with codes using other methods [Ng et al. 2008]. In our simulations, formation of plasmoids were not found, except when insufficient resolution was used, or when a small amount of noise was added externally. The reconnection rate is found to follow the Sweet-Parker scaling when no noise is added, but increases to a level independent of $S$ with noise, when plasmoids form. Latest results with $S$ up to $2 \times 10^5$ will be presented.

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

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
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