Researcher profile

J. O. Thurgood

J. O. Thurgood contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

8 published item(s)

preprint2016arXiv

Particle-in-cell simulations of the relaxation of electron beams in inhomogeneous solar wind plasmas

Previous theoretical considerations of electron beam relaxation in inhomogeneous plasmas have indicated that the effects of the irregular solar wind may account for the poor agreement of homogeneous modelling with the observations. Quasi-linear theory and Hamiltonian models based on Zakharov's equations have indicated that when a level of density fluctuations is above a given threshold, density irregularities act to de-resonate the beam-plasma interaction, restricting Langmuir wave growth on the expense of beam energy. This work presents the first fully kinetic particle-in-cell (PIC) simulations of beam relaxation under the influence of density irregularities. We aim to independently determine the influence of background inhomogeneity on the beam-plasma system, and to test theoretical predictions and alternative models using a fully kinetic treatment. We carry out 1D PIC simulations of a bump-on-tail unstable electron beam in the presence of increasing levels of background inhomogeneity using the fully electromagnetic, relativistic EPOCH PIC code. We find that in the case of homogeneous background plasma density, Langmuir wave packets are generated at the resonant condition and then quasi-liear relaxation leads to a dynamic increase of wavenumbers generated. No electron acceleration is seen - unlike in the inhomogeneous experiments, all of which produce high-energy electrons. For the inhomogeneous experiments we also observe the generation of backwards propagating Langmuir waves, which is shown directly to be due to the refraction of the packets off the density gradients. Our fully kinetic PIC simulations broadly confirm the findings of quasi-linear theory and the Hamiltonian model based on Zakharov's equations. Strong density fluctuations modify properties of excited Langmuir waves altering their dispersion properties.

preprint2015arXiv

Self-consistent particle-in-cell simulations of fundamental and harmonic plasma radio emission mechanisms

Aims. The simulation of three-wave interaction based plasma emission, thought to be the underlying mechanism for Type III solar radio bursts, is a challenging task requiring fully-kinetic, multi-dimensional models. This paper aims to resolve a contradiction in past attempts, whereby some studies indicate that no such processes occur. Methods. We self-consistently simulate three-waved based plasma emission through all stages by using 2D, fully kinetic, electromagnetic particle-in-cell simulations of relaxing electron beams using the EPOCH2D code. Results. Here we present the results of two simulations; Run 1 (nb/n0 = 0.0057, vb/Δvb = vb/Ve = 16) and Run 2 (nb/n0 = 0.05, vb/Δvb = vb/Ve = 8), which we find to permit and prohibit plasma emission respectively. We show that the possibility of plasma emission is contingent upon the frequency of the initial electrostatic waves generated by the bump-in-tail instability, and that these waves may be prohibited from participating in the necessary three-wave interactions due to frequency conservation requirements. In resolving this apparent contradiction through a comprehensive analysis, in this paper we present the first self-consistent demonstration of fundamental and harmonic plasma emission from a single-beam system via fully kinetic numerical simulation. We caution against simulating astrophysical radio bursts using unrealistically dense beams (a common approach which reduces run time), as the resulting non-Langmiur characteristics of the initial wave modes significantly suppresses emission. Comparison of our results also indicates that, contrary to the suggestions of previous authors, an alternative plasma emission mechanism based on two counter-propagating beams is unnecessary in an astrophysical context.

preprint2014arXiv

First direct measurements of transverse waves in solar polar plumes using SDO/AIA

There is intense interest in determining the precise contribution of Alfvénic waves propagating along solar structures to the problems of coronal heating and solar wind acceleration. Since the launch of SDO/AIA, it has been possible to resolve transverse oscillations in off-limb solar polar plumes and recently McIntosh et al. (2011, Nature, 475, 477) concluded that such waves are energetic enough to play a role in heating the corona and accelerating the fast solar wind. However, this result is based on comparisons to Monte Carlo simulations and confirmation via direct measurements is still outstanding. Thus, this letter reports on the first direct measurements of transverse wave motions in solar polar plumes. Over a 4 hour period, we measure the transverse displacements, periods and velocity amplitudes of 596 distinct oscillations observed in the 171 Åchannel of SDO/AIA. We find a broad range of non-uniformly distributed parameter values which are well described by log-normal distributions with peaks at $234$ km, $121$ s and $8$ km s$^{-1}$, and mean and standard deviations of $407\pm297$ km, $173\pm118$ s and $14\pm10$ km s$^{-1}$. Within standard deviations, our direct measurements are broadly consistent with previous results. However, accounting for the whole of our observed non-uniform parameter distribution we calculate an energy flux of $9-24$ W m$^{-2}$, which is $4-10$ times below the energy requirement for solar wind acceleration. Hence, our results indicate that transverse MHD waves as resolved by SDO/AIA cannot be the dominant energy source for fast solar wind acceleration in the open-field corona.

preprint2013arXiv

3D Alfven wave behaviour around proper and improper magnetic null points

Context: MHD waves and magnetic null points are both prevalent in many astrophysical plasmas, including the solar atmosphere. Interaction between waves and null points has been implicated as a possible mechanism for localised heating events. Aims: Here we investigate the transient behaviour of the Alfven wave about fully 3D proper and improper 3D magnetic null points. Previously, the behaviour of fast magnetoacoustic waves at null points in 3D, cold MHD was considered by Thurgood & McLaughlin (Astronomy & Astrophysics, 2012, 545, A9). Methods: We introduce an Alfven wave into the vicinity of both proper and improper null points by numerically solving the ideal, $β=0$ MHD equations using the LARE3D code. A magnetic fieldline and flux-based coordinate system permits the isolation of resulting wave-modes and the analysis of their interaction. Results: We find that the Alfven wave propagates throughout the region and accumulates near the fan-plane, causing current build up. For different values of null point eccentricity, the qualitative behaviour changes only by the imposition of anisotropic pulse dilation, due to the differing rates at which fieldlines diverge from the spine. For all eccentricities, we find that the fast and Alfven waves are linearly decoupled. During the driving phase, an independently propagating fast wave is nonlinearly generated due to the ponderomotive force. Subsequently, no further excitation of fast waves occurs. Conclusions: We find that the key aspects of the theory of Alfvén waves about 2D null points extends intuitively to the fully 3D case; i.e. the wave propagates along fieldlines and thus accumulates at predictable parts of the topology. We also highlight that unlike in the 2D case, in 3D Alfven-wave pulses are always toroidal and thus any aspects of 2D Alfven-wave-null models that are pulse-geometry specific must be reconsidered in 3D.

preprint2013arXiv

Nonlinear Alfvén wave dynamics at a 2D magnetic null point: ponderomotive force

Context : In the linear, β=0 MHD regime, the transient properties of MHD waves in the vicinity of 2D null points are well known. The waves are decoupled and accumulate at predictable parts of the magnetic topology: fast waves accumulate at the null point; whereas Alfvén waves cannot cross the separatricies. However, in nonlinear MHD mode conversion can occur at regions of inhomogeneous Alfvén speed, suggesting that the decoupled nature of waves may not extend to the nonlinear regime. Aims: We investigate the behaviour of low-amplitude Alfvén waves about a 2D magnetic null point in nonlinear, β= 0 MHD. Methods: We numerically simulate the introduction of low-amplitude Alfvén waves into the vicinity of a magnetic null point using the nonlinear LARE2D code. Results: Unlike in the linear regime, we find that the Alfvén wave sustains cospatial daughter disturbances, manifest in the transverse and longitudinal fluid velocity, owing to the action of nonlinear magnetic pressure gradients (viz. the ponderomotive force). These disturbances are dependent on the Alfvén wave and do not interact with the medium to excite magnetoacoustic waves, although the transverse daughter becomes focused at the null point. Additionally, an independently propagating fast magnetoacoustic wave is generated during the early stages, which transports some of the initial Alfvén wave energy towards the null point. Subsequently, despite undergoing dispersion and phase-mixing due to gradients in the Alfvén-speed profile (Grad c_A \neq 0) there is no further nonlinear generation of fast waves. Conclusions: We find that Alfvén waves at 2D cold null points behave largely as in the linear regime, however they sustain transverse and longitudinal disturbances - effects absent in the linear regime - due to nonlinear magnetic pressure gradients.

preprint2013arXiv

On ponderomotive effects induced by Alfven waves in inhomogeneous 2.5D MHD plasmas

Where spatial gradients in the amplitude of an Alfven wave are non-zero, a nonlinear magnetic-pressure gradient acts upon the medium (the ponderomotive force). We investigate the nature of such a force in inhomogeneous 2.5D MHD plasmas by analysing source terms in the nonlinear wave equations for the general case of inhomogeneous B and ρ, and consider supporting nonlinear numerical simulations. Our equations indicate there are two distinct classes of ponderomotive effect induced by Alfven waves in general 2.5D MHD, each with both a longitudinal and transverse manifestation; i) Geometric Effects: Gradients in the pulse geometry relative to the background magnetic field cause the wave to sustain cospatial disturbances, the longitudinal and transverse daughter disturbances - where we report on the transverse disturbance for the first time. ii) Grad (c_A) Effects: Where a pulse propagates through an inhomogeneous region (where the non-zero gradients in the Alfven-speed profile are non-zero), the nonlinear magnetic-pressure gradient acts to accelerate the plasma. Transverse gradients (phase mixing regions) excite independently propagating fast magnetoacoustic waves (generalising the result of Nakariakov et al. 1997, Solar Physics, 175, 93) and longitudinal gradients perturb along the field (thus creating static disturbances in β = 0, and slow waves in β \= 0). We additionally demonstrate that mode conversion due the nonlinear Lorentz force is a one-way process, and does not act as a mechanism to nonlinearly generate Alfven waves due to propagating magnetoacoustic waves. We conclude that these ponderomotive effects are induced by an Alfven wave propagating in any MHD medium, and have the potential to have significant consequences on the dynamics of energy transport and aspects of dissipation provided the system is sufficiently nonlinear and inhomogeneous.

preprint2012arXiv

Linear and nonlinear MHD mode coupling of the fast magnetoacoustic wave about a 3D magnetic null point

Context: Coronal magnetic null points have been implicated as possible locations for localised heating events in 2D models. We investigate this possibility about fully 3D null points. Aims: We investigate the nature of the fast magnetoacoustic wave about a fully 3D magnetic null point, with a specific interest in its propagation, and we look for evidence of MHD mode coupling and/or conversion to the Alfvén mode. Methods: A special fieldline and flux-based coordinate system was constructed to permit the introduction of a pure fast magnetoacoustic wave in the vicinity of proper and improper 3D null points. We considered the ideal, β = 0, MHD equations, which are solved using the LARE3D numerical code. The constituent modes of the resulting wave were isolated and identified using the special coordinate system. Numerical results were supported by analytical work derived from perturbation theory and a linear implementation of the WKB method. Results: An initially pure fast wave is found to be permanently decoupled from the Alfvén mode both linearly and nonlinearly for both proper and improper 3D null points. The pure fast mode also generates and sustains a nonlinear disturbance aligned along the equilibrium magnetic field. The resulting pure fast magnetoacoustic pulse has transient behaviour, which is found to be governed by the (equilibrium) Alfvén-speed profile, and a refraction effect focuses all the wave energy towards the null point. Conclusions: Thus, the main results from previous 2D work do indeed carry over to the fully 3D magnetic null points and

preprint2012arXiv

On the Periodicity of Oscillatory Reconnection

Oscillatory reconnection is a time-dependent magnetic reconnection mechanism that naturally produces periodic outputs from aperiodic drivers. This paper aims to quantify and measure the periodic nature of oscillatory reconnection for the first time. We solve the compressible, resistive, nonlinear MHD equations using 2.5D numerical simulations. We identify two distinct periodic regimes: the impulsive and stationary phases. In the impulsive phase, we find the greater the amplitude of the initial velocity driver, the longer the resultant current sheet and the earlier its formation. In the stationary phase, we find that the oscillations are exponentially decaying and for driving amplitudes 6.3 - 126.2 km/s, we measure stationary-phase periods in the range 56.3 - 78.9 s, i.e. these are high frequency (0.01 - 0.02 Hz) oscillations. In both phases, we find that the greater the amplitude of the initial velocity driver, the shorter the resultant period, but note that different physical processes and periods are associated with both phases. We conclude that the oscillatory reconnection mechanism behaves akin to a damped harmonic oscillator.