Source author record

P. H. Diamond

P. H. Diamond 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

22works
18topics
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

22 published item(s)

preprint2025arXiv

Disentangling Core and Edge Mechanisms of the Density Limit in DIII-D Negative Triangularity Plasmas

The density limit is investigated in the DIII-D negative triangularity (NT) plasmas which lack a standard H-mode edge. We find the limit may not be a singular disruptive boundary but a multifaceted density saturation phenomenon governed by distinct core and edge transport mechanisms. Sustained, non-disruptive operation is achieved at densities up to 1.8 times the Greenwald limit ($n_\mathrm{G}$) until the termination of auxiliary heating. Systematic power scans reveal distinct power scalings for the core ($n_e \propto P_\mathrm{SOL}^{0.27\pm0.03}$) and edge ($n_e \propto P_\mathrm{SOL}^{0.42\pm0.04}$) density limits. The edge density saturation is triggered abruptly by the onset of a non-disruptive, high-field side radiative instability that clamps the edge density below $n_\mathrm{G}$. In contrast, the core density continues to rise until it saturates, a state characterized by substantially enhanced core turbulence. Core transport evolves from a diffusive to an intermittent, avalanche-like state, as indicated by heavy-tailed probability density functions (kurtosis $\approx 6$), elevated Hurst exponents, and a $1/f$-type power spectrum. These findings suggest that the density limit in the low-confinement regime is determined by a combination of edge radiative instabilities and core turbulent transport. This distinction provides separate targets for control strategies aimed at extending the operational space of future fusion devices.

preprint2020arXiv

Turbulence model reduction by deep learning

A central problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce a new, data-driven method for parameterizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local \emph{gradient} of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

preprint2016arXiv

Cascades and Spectra of a Turbulent Spinodal Decomposition in 2D Symmetric Binary Liquid Mixture

We study the fundamental physics of cascades and spectra in 2D Cahn-Hilliard-Navier-Stokes (CHNS) turbulence, and compare and contrast this system with 2D MagnetoHydroDynamic (MHD) turbulence. The important similarities include basic equations, ideal quadratic invariants, cascades and the role of linear elastic waves. Surface tension induces elasticity, and the balance between surface tension energy and turbulent kinetic energy determines a length scale (Hinze scale) of the system. The Hinze scale may be thought of as the scale of emergent critical balance between fluid straining and elastic restoring forces. The scales between the Hinze scale and dissipation scale constitute the elastic range of the 2D CHNS system. By direct numerical simulation, we find that in the elastic range, the mean square concentration spectrum $H^ψ_k$ of the 2D CHNS system exhibits the same power law ($-7/3$) as the mean square magnetic potential spectrum $H^A_k$ in the inverse cascade regime of 2D MHD. This power law is consistent with an inverse cascade of $H^ψ$, which is observed. The kinetic energy spectrum of the 2D CHNS system is $E^K_k\sim k^{-3}$ if forced at large scale, suggestive of the direct enstrophy cascade power law of 2D Navier-Stokes (NS) turbulence. The difference from the energy spectra of 2D MHD turbulence implies that the back reaction of the concentration field to fluid motion is limited. We suggest this is because the surface tension back reaction is significant only in the interfacial regions. The interfacial regions fill only a small portion of the 2D CHNS system, and their interface packing fraction is much smaller than that for 2D MHD.

preprint2016arXiv

Intrinsic rotation drive by collisionless trapped electron mode turbulence

Both the parallel residual stress and parallel turbulent acceleration driven by electrostatic collisionsless trapped electron mode (CTEM) turbulence are calculated analytically using gyrokinetic theory. Quasilinear results show that the parallel residual stress contributes an outward flux of co-current rotation for normal magnetic shear and turbulence intensity profile increasing outward. This may induce intrinsic counter-current rotation or flattening of the co-current rotation profile. The parallel turbulent acceleration driven by CTEM turbulence vanishes, due to the absence of a phase shift between density fluctuation and ion pressure fluctuation. This is different from the case of ion temperature gradient (ITG) turbulence, for which the turbulent acceleration can provide co-current drive for normal magnetic shear and turbulence intensity profile increasing outward. Its order of magnitude is predicted to be the same as that of the divergence of the residual stress [Lu Wang and P.H. Diamond, Phys. Rev. Lett. {\bf 110}, 265006 (2013)]. A possible connection of these theoretical results to experimental observations of electron cyclotron heating effects on toroidal rotation is discussed.

preprint2016arXiv

Plus Charge Prevalence in Cosmic Rays: Room for Dark Matter in the Positron Spectrum

The unexpected energy spectrum of the positron/electron ratio is interpreted astrophysically, with a possible exception of the 100-300 GeV range. The data indicate that this ratio, after a decline between $0.5-8$ GeV, rises steadily with a trend towards saturation at 200-400GeV. These observations (except for the trend) appear to be in conflict with the diffusive shock acceleration (DSA) mechanism, operating in a \emph{single} supernova remnant (SNR) shock. We argue that $e^{+}/e^{-}$ ratio can still be explained by the DSA if positrons are accelerated in a \emph{subset} of SNR shocks which: (i) propagate in clumpy gas media, and (ii) are modified by accelerated CR \emph{protons}. The protons penetrate into the dense gas clumps upstream to produce positrons and, \emph{charge the clumps positively}. The induced electric field expels positrons into the upstream plasma where they are shock-accelerated. Since the shock is modified, these positrons develop a harder spectrum than that of the CR electrons accelerated in other SNRs. Mixing these populations explains the increase in the $e^{+}/e^{-}$ ratio at $E>8$ GeV. It decreases at $E<8$ GeV because of a subshock weakening which also results from the shock modification. Contrary to the expelled positrons, most of the antiprotons, electrons, and heavier nuclei, are left unaccelerated inside the clumps. Scenarios for the 100-300 GeV AMS-02 fraction exceeding the model prediction, including, but not limited to, possible dark matter contribution, are also discussed.

preprint2016arXiv

Poloidal rotation driven by nonlinear momentum transport in strong electrostatic turbulence

Virtually, all existing theoretical works on turbulent poloidal momentum transport are based on quasilinear theory. Nonlinear poloidal momentum flux - $\langle \tilde{v}_r \tilde{n} \tilde{v}_θ \rangle$ is universally neglected. However, in the strong turbulence regime where relative fluctuation amplitude is no longer small, quasilinear theory is invalid. This is true at the all-important plasma edge. In this work, nonlinear poloidal momentum flux $ \langle \tilde{v}_r \tilde{n} \tilde{v}_θ \rangle $ in strong electrostatic turbulence is calculated using Hasegawa-Mima equation, and is compared with quasilinear poloidal Reynolds stress. A novel property is that symmetry breaking in fluctuation spectrum is not necessary for a nonlinear poloidal momentum flux. This is fundamentally different from the quasilinear Reynold stress. Furthermore, the comparison implies that the poloidal rotation drive from the radial gradient of nonlinear momentum flux is comparable to that from the quasilinear Reynolds force. Nonlinear poloidal momentum transport in strong electrostatic turbulence is thus not negligible for poloidal rotation drive, and so may be significant to transport barrier formation.

preprint2016arXiv

Zonal Flow Patterns: How Toroidal Coupling Induces Phase Jumps and Shear Layers

A new, frequency modulation mechanism for zonal flow pattern formation is presented. The model predicts the probability distribution function of the flow strength as well as the evolution of the characteristic spatial scale. Magnetic toroidicity-induced global phase dynamics is shown to determine the spatial structure of the flow. A key result is the observation that global phase patterning can lead to zonal flow formation in the absence of turbulence inhomogeneity.

preprint2015arXiv

Ion-acoustic Shocks with Self-Regulated Ion Reflection and Acceleration

An analytic solution describing an ion-acoustic collisionless shock, self-consistently with the evolution of shock-reflected ions, is obtained. The solution extends the classic soliton solution beyond a critical Mach number, where the soliton ceases to exist because of the upstream ion reflection. The reflection transforms the soliton into a shock with a trailing wave and a foot populated by the reflected ions. The solution relates parameters of the entire shock structure, such as the maximum and minimum of the potential in the trailing wave, the height of the foot, as well as the shock Mach number, to the number of reflected ions. This relation is resolvable for any given distribution of the upstream ions. In this paper, we have resolved it for a simple "box" distribution. Two separate models of electron interaction with the shock are considered. The first model corresponds to the standard Boltzmannian electron distribution in which case the critical shock Mach number only insignificantly increases from M=1.6 (no ion reflection) to M=1.8 (substantial reflection). The second model corresponds to adiabatically trapped electrons. They produce a stronger increase, from M=3.1 to M=4.5. The shock foot that is supported by the reflected ions also accelerates them somewhat further. A self-similar foot expansion into the upstream medium is also described analytically.

preprint2015arXiv

Nonlinear parallel momentum transport in strong turbulence

Most existing theoretical studies of momentum transport focus on calculating the Reynolds stress based on quasilinear theory, without considering the \emph{nonlinear} momentum flux-$<\tilde{v}_r \tilde{n} \tilde{u}_{\|} >$. However, a recent experiment on TORPEX found that the nonlinear toroidal momentum flux induced by blobs makes a significant contribution as compared to the Reynolds stress [Labit et al., Phys. Plasmas {\bf 18}, 032308 (2011)]. In this work, the nonlinear parallel momentum flux in strong turbulence is calculated by using three dimensional Hasegawa-Mima equation. It is shown that nonlinear diffusivity is smaller than quasilinear diffusivity from Reynolds stress. However, the leading order nonlinear residual stress can be comparable to the quasilinear residual stress, and so could be important to intrinsic rotation in tokamak edge plasmas. A key difference from the quasilinear residual stress is that parallel fluctuation spectrum asymmetry is not required for nonlinear residual stress.

preprint2014arXiv

Anomalous Viscosity of the Quark-Gluon Plasma

The shear viscosity of the quark-gluon plasma is predicted to be lower than the collisional viscosity for weak coupling. The estimated ratio of the shear viscosity to entropy density is rather close to the ratio calculated by N = 4 super Yang-Mills theory for strong coupling, which indicates that the quark-gluon plasma might be strongly coupled. However, in presence of momentum anisotropy, the Weibel instability can arise and drive the turbulent transport. Shear viscosity can be lowered by enhanced collisionality due to turbulence, but the decorrelation time and its relation to underlying dynamics and color-magnetic fields have not been calculated self-consistently. In this paper, we use resonance broadening theory for strong turbulence to calculate the anomalous viscosity of the quark-gluon plasma for nonequilibrium. For saturated Weibel instability, we estimate the scalings of the decorrelation rate and viscosity and compare these with collisional transport. This calculation yields an explicit connection between the underlying momentum space anisotropy and the viscosity anomaly.

preprint2013arXiv

Analytic Solution for Self-regulated Collective Escape of Cosmic Rays from their Acceleration Sites

Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CR), are believed to maintain an average CR spectrum by diffusive shock acceleration (DSA) regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular paradigm of the CR injection into the ISM, according to which they passively and isotropically diffuse in the prescribed magnetic fluctuations as test particles. First, we treat the escaping CR and the Alfven waves excited by them on an equal footing. Second, we adopt field aligned CR escape outside the source, where the waves become weak. An exact analytic self-similar solution for a CR "cloud" released by a dimmed accelerator strongly deviates from the test-particle result. The CR diffusion coefficient $D_{NL}$ is strongly suppressed compared to its background ISM value $D_{ISM}$: $D_{NL}\sim D_{ISM}\exp(-Π)<< D_{ISM}$ for sufficiently high field-line-integrated CR partial pressure, $Π$. When $Π>>1$, the CRs drive Alfven waves efficiently enough to build a transport barrier that strongly reduces the leakage. The solution has a spectral break at $p=p_{br}$, where $p_{br}$ satisfies the following equation $D_{NL}(p_{br})\simeq z^{2}/t$.

preprint2013arXiv

Impact of Resonant Magnetic Perturbations on Zonal Modes, Drift-Wave Turbulence and the L-H Transition Threshold

We study the effects of Resonant Magnetic Perturbations (RMPs) on turbulence, flows and confinement in the framework of resistive drift-wave turbulence. This work was motivated, in parts, by experiments reported at the IAEA 2010 conference [Y. Xu {\it et al}, Nucl. Fusion \textbf{51}, 062030] which showed a decrease of long-range correlations during the application of RMPs. We derive and apply a zero-dimensional predator-prey model coupling the Drift-Wave Zonal Mode system [M. Leconte and P.H. Diamond, Phys. Plasmas \textbf{19}, 055903] to the evolution of mean quantities. This model has both density gradient drive and RMP amplitude as control parameters and predicts a novel type of transport bifurcation in the presence of RMPs. This model allows a description of the full L-H transition evolution with RMPs, including the mean sheared flow evolution. The key results are: i) The L-I and I-H power thresholds \emph{both} increase with RMP amplitude $|\bx|$, the relative increase of the L-I threshold scales as $ΔP_{\rm LI} \propto |\bx|^2 ν_*^{-2} \gyro^{-2}$, where $ν_*$ is edge collisionality and $\gyro$ is the sound gyroradius. ii) RMPs are predicted to \emph{decrease} the hysteresis between the forward and back-transition. iii) Taking into account the mean density evolution, the density profile - sustained by the particle source - has an increased turbulent diffusion compared with the reference case without RMPs which provides one possible explanation for the \emph{density pump-out} effect.

preprint2012arXiv

On the mechanism for breaks in the cosmic ray spectrum

The proof of cosmic ray (CR) origin in supernova remnants (SNR) must hinge on full consistency of the CR acceleration theory with the observations; direct proof is impossible because of the orbit stochasticity of CR particles. Recent observations of a number of galactic SNR strongly support the SNR-CR connection in general and the Fermi mechanism of CR acceleration, in particular. However, many SNR expand into weakly ionized dense gases, and so a significant revision of the mechanism is required to fit the data. We argue that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by a partial evanescence of Alfven waves that confine particles to the accelerator. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated. Using the recent Fermi spacecraft observation of the SNR W44 as an example, we demonstrate that the parent proton spectrum is a classical test particle power law $\propto E^{-2}$, steepening to $E^{-3}$ at $E_{br}\approx7GeV$.

preprint2012arXiv

Proton-Helium Spectral Anomaly as a Signature of Cosmic Ray Accelerator

The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants (SNR) must hinge on full consistency of acceleration theory with the observations; direct proof is impossible because of the orbit scrambling of CR particles. The PAMELA orbital telescope revealed deviation between helium and proton CR spectra deemed inconsistent with the theory, since the latter does not differentiate between elements of ultrarelativistic rigidity. By considering an initial (injection-) phase of the diffusive shock acceleration (DSA), where elemental similarity does not apply, we demonstrate that the spectral difference is, in fact, a unique signature of the DSA. Collisionless plasma SNR shocks inject more He2+ relative to protons when they are stronger and so produce harder helium spectra. The injection bias is due to Alfven waves driven by the more abundant protons, so the He2+ ions are harder to trap by these waves because of the larger gyroradii. By fitting the p/He ratio to the PAMELA data, we bolster the DSA-case for resolving the century-old mystery of CR origin.

preprint2011arXiv

Magnetic and density spikes in cosmic ray shock precursors

In shock precursors populated by accelerated cosmic rays (CR), the CR return current instability is believed to significantly enhance the pre-shock perturbations of magnetic field. We have obtained fully-nonlinear exact ideal MHD solutions supported by the CR return current. The solutions occur as localized spikes of circularly polarized Alfven envelopes (solitons, or breathers). As the conventional (undriven) solitons, the obtained magnetic spikes propagate at a speed $C$ proportional to their amplitude, $C=C_{A}B_{\rm max}/\sqrt{2}B_{0}$. The sufficiently strong solitons run thus ahead of the main shock and stand in the precursor, being supported by the return current. This property of the nonlinear solutions is strikingly different from the linear theory that predicts non-propagating (that is, convected downstream) circularly polarized waves. The nonlinear solutions may come either in isolated pulses (solitons) or in soliton-trains (cnoidal waves). The morphological similarity of such quasi-periodic soliton chains with recently observed X-ray stripes in Tycho supernova remnant (SNR) is briefly discussed. The magnetic field amplification determined by the suggested saturation process is obtained as a function of decreasing SNR blast wave velocity during its evolution from the ejecta-dominated to the Sedov-Taylor stage.

preprint2011arXiv

Mechanism for spectral break in cosmic ray proton spectrum from Supernova remnant W44

Recent observations of the supernova remnant W44 by the \emph{Fermi} spacecraft observatory strongly support the idea that the bulk of galactic cosmic rays is accelerated in such remnants by a Fermi mechanism, also known as diffusive shock acceleration. However, the W44 expands into weakly ionized dense gas, and so a significant revision of the mechanism is required. In this paper we provide the necessary modifications and demonstrate that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by Alfven wave evanescence leading to the fractional particle losses. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated and successfully fitted to the Fermi Observatory data. The parent proton spectrum is best represented by a classical test particle power law E^-2, steepening to E^-3 at E_br~7GeV due to deteriorated particle confinement.

preprint2010arXiv

On the Structure and Scale of Cosmic Ray Modified Shocks

Strong astrophysical shocks, diffusively accelerating cosmic rays (CR) ought to develop CR precursors. The length of such precursor $L_{p}$ is believed to be set by the ratio of the CR mean free path $λ$ to the shock speed, i.e., $L_{p}\sim cλ/V_{sh}\sim cr_{g}/V_{sh}$, which is formally independent of the CR pressure $P_{c}$. However, the X-ray observations of supernova remnant shocks suggest that the precursor scale may be significantly shorter than $L_{p}$ which would question the above estimate unless the magnetic field is strongly amplified and the gyroradius $r_{g}$ is strongly reduced over a short (unresolved) spatial scale. We argue that while the CR pressure builds up ahead of the shock, the acceleration enters into a strongly nonlinear phase in which an acoustic instability, driven by the CR pressure gradient, dominates other instabilities (at least in the case of low $β$ plasma). In this regime the precursor steepens into a strongly nonlinear front whose size scales with \emph{the CR pressure}as $L_{f}\sim L_{p}\cdot(L_{s}/L_{p})^{2}(P_{c}/P_{g})^{2}$, where $L_{s}$ is the scale of the developed acoustic turbulence, and $P_{c}/P_{g}$ is the ratio of CR to gas pressure. Since $L_{s}\ll L_{p}$, the precursor scale reduction may be strong in the case of even a moderate gas heating by the CRs through the acoustic and (possibly also) the other instabilities driven by the CRs.

preprint2010arXiv

Probing Nearby CR Accelerators and ISM Turbulence with Milagro Hot Spots

Both the acceleration of cosmic rays (CR) in supernova remnant shocks and their subsequent propagation through the random magnetic field of the Galaxy deem to result in an almost isotropic CR spectrum. Yet the MILAGRO TeV observatory discovered a sharp ($\sim10^{\circ})$ arrival anisotropy of CR nuclei. We suggest a mechanism for producing a weak and narrow CR beam which operates en route to the observer. The key assumption is that CRs are scattered by a strongly anisotropic Alfven wave spectrum formed by the turbulent cascade across the local field direction. The strongest pitch-angle scattering occurs for particles moving almost precisely along the field line. Partly because this direction is also the direction of minimum of the large scale CR angular distribution, the enhanced scattering results in a weak but narrow particle excess. The width, the fractional excess and the maximum momentum of the beam are calculated from a systematic transport theory depending on a single scale $l$ which can be associated with the longest Alfven wave, efficiently scattering the beam. The best match to all the three characteristics of the beam is achieved at $l\sim1$pc. The distance to a possible source of the beam is estimated to be within a few 100pc. Possible approaches to determination of the scale $l$ from the characteristics of the source are discussed. Alternative scenarios of drawing the beam from the galactic CR background are considered. The beam related large scale anisotropic CR component is found to be energy independent which is also consistent with the observations.

preprint1997arXiv

Dissipative Dynamics of Collisionless Nonlinear Alfven Wave Trains

The nonlinear dynamics of collisionless Alfven trains, including resonant particle effects is studied using the kinetic nonlinear Schroedinger (KNLS) equation model. Numerical solutions of the KNLS reveal the dynamics of Alfven waves to be sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. The combined effects of both wave nonlinearity and Landau damping result in the evolutionary formation of stationaryOA S- and arc-polarized directional and rotational discontinuities. These waveforms are freqently observed in the interplanetary plasma.

preprint1996arXiv

Fluid Models for Kinetic Effects on Coherent Nonlinear Alfven Waves. II. Numerical Solutions

The influence of various kinetic effects (e.g. Landau damping, diffusive and collisional dissipation, and finite Larmor radius terms) on the nonlinear evolution of finite amplitude Alfvenic wave trains in a finite-beta environment is systematically investigated using a novel, kinetic nonlinear Schrodinger (KNLS) equation. The dynamics of Alfven waves is sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. Numerical solution for the case with Landau damping reveals the formation of dissipative structures, which are quasi-stationary, S-polarized directional (and rotational) discontinuities which self-organize from parallel propagating, linearly polarized waves. Parallel propagating circularly polarized packets evolve to a few circularly polarized Alfven harmonics on large scales. Stationary arc-polarized rotational discontinuities form from obliquely propagating waves. Collisional dissipation, even if weak, introduces enhanced wave damping when beta is very close to unity. Cyclotron motion effects on resonant particle interactions introduce cyclotron resonance into the nonlinear Alfven wave dynamics.

preprint1995arXiv

Suppression of Cross-Field Transport of a Passive Scalar in Two-Dimensional Magnetohydrodynamic Turbulence

The theory of passive scalar transport in two dimensional turbulent fluids is generalized to the case of 2D MHD. Invariance of the cross correlation of scalar concentration and magnetic potential produces a novel contribution to the concentration flux. This pinch effect is proportional to the mean potential gradient, and is shown to drastically reduce transport of the passive scalar across the mean magnetic field when . Transport parallel to the mean magnetic field is unchanged. Implications for models of transport in turbulent magnetofluids are discussed. PAC NOS. 47.25.Jn, 47.65.+a