No description available.

The study of low-frequency turbulence in magnetized plasmas is a difficult problem due to both the enormous range of scales involved and the variety of physics encompassed over this range. Much of the progress that has been made in turbulence theory is based upon a result from incompressible magnetohydrodynamics (MHD), in which energy is only transferred from large scales to small via the collision of Alfv ́n waves propagating oppositely along the mean magnetic field. Improvements in laboratory devices and satellite measurements have demonstrated that, while theories based on this premise are useful over inertial ranges, describing turbulence at scales that approach particle gyroscales requires new theory. In this thesis, we examine the limits of incompressible MHD theory in describing collisions between pairs of Alfvén waves. This interaction represents the fundamental unit of plasma turbulence. To study this interaction, we develop an analytic theory describing the nonlinear evolution of interacting Alfv ́n waves and compare this theory to simulations performed using the gyrokinetic code AstroGK. Gyrokinetics captures a much richer set of physics than that described by incompressible MHD, and is well-suited to describing Alfvénic turbulence around the ion gyroscale. We demonstrate that AstroGK is well suited to the study of physical Alfvén waves by reproducing laboratory Alfvén dispersion data collected using the LAPD. Additionally, we have developed an initialization alogrithm for use with AstroGK that allows exact Alfvén eigenmodes to be initialized with user specified amplitudes and phases. We demonstrate that our analytic theory based upon incompressible MHD gives excellent agreement with gyrokinetic simulations for weakly turbulent collisions in the limit that k⊥ ρi << 1. In this limit, agreement is observed in the time evolution of nonlinear products, and in the strength of nonlinear interaction with respect to polarization and scale. We also examine the effect of wave amplitude upon the validity of our analytic solution, exploring the nature of strong turbulence. In the kinetic limit where k⊥ ρi ≥ 1 where incompressible MHD is no longer a valid description, we illustrate how the nonlinear evolution departs from our analytic expression. The analytic theory we develop provides a framework from which more sophisticated of weak and strong inertial-range turbulence theories may be developed. Characterization of the limits of this theory may provide guidance in the development of kinetic Alfvén wave turbulence.

No description available.

The exact mechanisms responsible for heating the solar atmosphere in regions such as the chromosphere (partially ionized) and the corona (fully ionized) remain quantitatively unknown. This thesis demonstrates that the ions can be heated by Alfven waves with low frequencies in fully and partially ionized low beta plasmas, which is contrary to the customary expectation. For the partially ionized case, we find the heating process to be less efficient than the scenario with no ion-neutral collisions, and that the heating efficiency depends on the ratio of ion-neutral collision frequency to the ion gyrofrequency. For Alfven waves propagating obliquely to the background magnetic field in fully ionized plasmas, we find the heating process to be more efficient than the situation with Alfven waves propagating along the background magnetic field. Furthermore, the simulation results show the parallel kinetic temperature can become even larger than the perpendicular component for the case of obliquely propagating Alfven waves.

Possibilities of generating plasma turbulence to provide control of space weather processes have been of particular interest in recent years. Such turbulence can be created by chemical released into a magnetized background plasma. The released plasma clouds are heavy ions which have ring velocity distribution and large free energy to drive the turbulence. An electromagnetic hybrid (fluid electrons and particle ions) model incorporating electron inertia is developed to study the generation and nonlinear evolution of this turbulence. Fourier pseudo-spectral methods are combined with finite difference methods to solve the electron momentum equations. Time integration is accomplished by a 4th-order Runge-Kutta scheme or predicator-corrector method. The numerical results show good agreement with theoretical prediction as well as provide further insights on the nonlinear turbulence evolution. Initially the turbulence lies near harmonics of the ring plasma ion cyclotron frequency and propagates nearly perpendicular to the background magnetic field as predicted by the linear theory. If the amplitude of the turbulence is sufficiently large, the quasi-electrostatic short wavelength ion cyclotron waves evolve nonlinearly into electromagnetic obliquely propagating shear Alfven waves with much longer wavelength. The results indicate that ring densities above a few percent of the background plasma density may produce wave amplitudes large enough for such an evolution to occur. The extraction of energy from the ring plasma may be in the range of 10-15% with a generally slight decrease in the magnitude as the ring density is increased from a few percent to several 10's of percent of the background plasma density. Possibilities to model the effects of nonlinear processes on energy extraction by introducing electron anomalous resistivity are also addressed. Suitability of the nonlinearly generated shear Alfven waves for applications to scattering radiation belt particles is discussed. / Master of Science

The solar wind, magnetosphere and ionosphere are coupled by magnetohydrodynamic waves, and this gives rise to new and often unexpected behaviours that cannot be produced by a single, isolated part of the system. This thesis examines two broad instances of coupling: field-line resonance (FLR) which couples fast and Alfvén waves, and magnetosphere-ionosphere (MI-) coupling via Alfvén waves. The first part of this thesis investigates field-line resonance for equilibria that vary in two dimensions perpendicular to the background magnetic field. This research confirms that our intuitive understanding of FLR from 1D is a good guide to events in 2D, and places 2D FLR onto a firm mathematical basis by systematic solution of the governing equations. It also reveals the new concept of ‘imprinting’ of spatial forms: spatial variations of the resonant Alfvén wave correlate strongly with the spatial form of the fast wave that drives the resonance. MI-coupling gives rise to ionosphere-magnetosphere (IM-) waves, and we have made a detailed analysis of these waves for a 1D sheet E-region. IM-waves are characterised by two quantities: a speed v_{IM} and an angular frequency ω_{IM} , for which we have obtained analytic expressions. For an ideal magnetosphere, IM-waves are advective and move in the direction of the electric field with speed v_{IM}. The advection speed is a non-linear expression that decreases with height-integrated E-region plasma-density, hence, wavepackets steepen on their trailing edge, rapidly accessing small length-scales through wavebreaking. Inclusion of electron inertial effects in the magnetosphere introduces dispersion to IM-waves. In the strongly inertial limit (wavelength λ << λ_{e} , where λ_{e} is the electron inertial length at the base of the magnetosphere), the group velocity of linear waves goes to zero, and the waves oscillate at ω_{IM} which is an upper limit on the angular frequency of IM-waves for any wavelength. Estimates of v_{IM} show that this speed can be a significant fraction (perhaps half) of the E_{⊥} × B_{0} drift in the E-region, producing speeds of up to several hundred metres per second. The upper limit on angular frequency, ωIM , is estimated to give periods from a few hundredths of a second to several minutes. IM-waves are damped by recombination and background ionisation, giving an e-folding decay time that can vary from tens of seconds to tens of minutes. We have also investigated the dynamics and steady-states that occur when the magnetosphere-ionosphere system is driven by large-scale Alfvénic field-aligned currents. Steady-states are dominated by two approximate solutions: an ‘upper’ solution that is valid in places where the E-region is a near perfect conductor, and a ‘lower’ solution that is valid where E-region depletion makes recombination negligible. These analytic solutions are extremely useful tools and the global steady-state can be constructed by matching these solutions across suitable boundary-layers. Furthermore, the upper solution reveals that E-region density cavities form and widen (with associated broadening of the magnetospheric downward current channel) if the downward current density exceeds the maximum current density that can be supplied by background E-region ionisation. We also supply expressions for the minimum E-region plasma-density and shortest length-scale in the steady-state. IM-waves and steady-states are extremely powerful tools for interpreting MI-dynamics. When an E-region density cavity widens through coupling to an ideal, single-fluid MHD magnetosphere, it does so by forming a discontinuity that steps between the upper and lower steady-states. This discontinuity acts as part of an ideal IM-wave and moves in the direction of the electric field at a speed U = \sqrt{v_{IM} {+} v_{IM} {-}}, which is the geometric mean of v_{IM} evaluated immediately to the left and right of the discontinuity. This widening speed is typically several hundreds of metres per second. If electron inertial effects are included in the magnetosphere, then the discontinuity is smoothed, and a series of undershoots and overshoots develops behind it. These undershoots and overshoots evolve as inertial IM-waves. Initially they are weakly inertial, with a wavelength of about λ_{e}, however, strong gradients of ω_{IM} cause IM-waves to phase-mix, making their wavelength inversely proportional to time. Therefore, the waves rapidly become strongly inertial and oscillate at ω_{IM}. The inertial IM-waves drive upgoing Alfvén waves in the magnetosphere, which populate a region over the downward current channel, close to its edge. In this manner, the E-region depletion mechanism, that we have detailed, creates small-scale Alfvén waves in large-scale current systems, with properties determined by MI-coupling.

Rotating convection in the Earth's core produces columnar vortices of radius ~10 km or less near the inner core boundary. Small-scale motions in the core can travel as Alfvén waves in the face of Ohmic diffusion, provided the ratio of the magnetic diffusion time th to the Alfvén wave travel time tA (measured by the Lundquist number S0) is much greater than unity. These motions transfer angular momentum from the core to the mantle, a process that can help explain variations in length of day. Vortices subject to the combined influence of a magnetic field and background rotation give rise to fast and slow Magneto-Coriolis (MC) waves whose damping is not well understood. This thesis investigates the long-time evolution of magneto hydrodynamic (MHD) waves generated by an isolated, small-scale motion in an otherwise quiescent, electrically conducting fluid. The first part of the study focuses on the damping of small-scale Alfvén waves, which is independent of rotation. For a plausible magnetic field strength in the Earth's core, it is shown that flows of lengthscale ~ 5 km or larger can propagate across the core as damped Alfvén waves on sub-decadal timescales. The second part of the study looks at MC waves generated from an isolated blob under rotation and a uniform axial magnetic field. The decay laws for these waves are obtained by considering the decay of fast and slow waves individually. While the fast waves are subject to strongly anisotropic magnetic diffusion, the slow waves diffuse isotopically. New timescales are derived for the onset of damping and the transition from the wave-dominated to the diffusion-dominated (quasi-static) phase of decay. This study shows for the first time that MC waves originating from small-scale vortices of magnetic Reynolds number Rm ~ 1 can be long-lived. The results of this study are extendible to small-scale MHD turbulence under rotation, whose damped wave phase has not been adequately addressed in the literature. Furthermore, it is thought that this study would help place a lower bound on the poloidal magnetic field strength in the Earth’s core.

We consider a simplified model based on the Navier-Stokes-Fourier system coupled to a transport equation and the Maxwell system, proposed to describe radiative flows in stars. We establish global- in-time existence for the associated initial-boundary value problem in the framework of weak solutions. Next, we study a hydrodynamical model describing the motion of internal stellar layers based on compressible Navier-Stokes-Fourier-Poisson system. We suppose that the medium is electrically charged, we include energy exchanges through radiative transfer and we assume that the system is steadily rotating. We analyze the singular limit of this system when the Mach number, the Alfven number, the Peclet number and the Froude number go to zero in a certain way and prove convergence to a 3D incompressible MHD system with a stationary linear transport equation for transport of radiation intensity. Finally, we show that the energy equation reduces to a steady equation for the temperature corrector.