Linear Wave Modes Research Articles

Dynamics and modulation of cosmic ray modified magnetosonic waves in a galactic gaseous rotating plasma

The influence of the presence of cosmic fluid on the magnetosonic waves and modulation instabilities in the interstellar medium (ISM) of spiral galaxies is investigated. The fluid model is developed by modifying the pressure equation in such dissipative rotating magnetoplasmas incorporating thermal ionized gas and cosmic rays. Applying the normal mode analysis, a modified dispersion relation is derived to study linear magnetosonic wave modes and their instabilities. The cosmic rays influence the wave damping by accelerating the damping rate. The standard reductive perturbation method is employed in the fluid model leading to a Korteweg–de Vries–Burgers (KdVB) equation in the small-amplitude limit. Several nonlinear wave shapes are assessed by solving the KdVB equation, analytically and numerically. The cosmic ray diffusivity and magnetic resistivity are responsible for the generation of shock waves. The modulational instability (MI) and the rogue wave solutions of the magnetosonic waves are studied by deriving a nonlinear Schrödinger equation from the obtained KdVB equation under the assumption that the cosmic ray diffusion and magnetic resistivity are weak and the carrier wave frequency is considerably lower than the wave frequency. The influence of various plasma parameters on the growth rate of MI is examined. The modification of the pressure term due to cosmic fluid reduces the MI growth in the interstellar medium. In addition, a quantitative analysis of the characteristics of rogue wave solutions is presented. Our investigation's applicability to the interstellar medium of spiral galaxies is traced out.

Fast and slow magnetosonic shock waves in non-Maxwellian plasmas

Obliquely propagating nonlinear fast and slow magnetosonic wave modes in a hot non-Maxwellian dissipative plasma are investigated in the current work. Modified temperatures have been derived for the non-extensive Q- and (r,q)-distributions that correspond to the physical properties of such non-Maxwellian plasmas. The reductive perturbation technique has been employed to derive the linear dispersion relation (LDR) and Kadomstev-Petvashvilli-Burgers (KPB) equation for slow and fast magnetosonic wave modes in two dimensions. We then investigated the effect of non-extensive parameter Q, spectral indices (r,q) and kinematic viscosity ν on the LDR and nonlinear propagation of KPB shock profiles for both the slow and fast mode magnetosonic waves. We found that linear and nonlinear propagation of fast and slow magnetosonic wave modes have been considerably modified in such non-Maxwellian plasmas. The results presented here would depict a realistic picture of the propagation of linear and nonlinear fast and slow magnetosonic wave modes in non-Maxwellian plasmas.

Legolas 2.0: Improvements and extensions to an MHD spectroscopic framework

We report on recent extensions and improvements to the Legolas code, which is an open-source, finite element-based numerical framework to solve the linearised (magneto)hydrodynamic equations for a three-dimensional force- and thermally balanced state with a nontrivial one-dimensional variation. The standard Fourier modes imposed give rise to a complex, generalised non-Hermitian eigenvalue problem which is solved to quantify all linear wave modes of the given system in either Cartesian or cylindrical geometries. The framework now supports subsystems of the eight linearised MHD equations, allowing for pure hydrodynamic setups, only one-dimensional density/temperature/velocity variations, or the option to treat specific closure relations. We discuss optimisations to the internal datastructure and eigenvalue solvers, showing a considerable performance increase in both execution time and memory usage. Additionally the code now has the capability to fully visualise eigenfunctions associated with given wave modes in multiple dimensions, which we apply to standard Kelvin-Helmholtz and Rayleigh-Taylor instabilities in hydrodynamics, thereby providing convincing links between linear stability analysis and the onset of non-linear phenomena.

Theory of magnetic turbulence and shock formation induced by a collisionless plasma instability

Magnetic fields are ubiquitous in universe, space, and laboratory plasmas. Especially, self-generated magnetic fields are important to know the mind of nature. The formation of Weibel-mediated collisionless shock is studied theoretically as a structure formation by the linear plasma wave growth, nonlinear saturation, and mode–mode coupling. Following a series of computer simulations and experimental studies of the physics, a simple model equation is proposed here to describe the time evolution of magnetic turbulence. Weibel instability is saturated by magnetic pressure, and thicker filaments continue to be generated by current coalescence (magnetic reconnection) mechanism. The model equation concludes the fact that the filament spacing increases linearly in time, and the magnetic energy power spectrum is given as Bk2 ∝ k−2. The time evolution of the turbulence is characterized with the cascade toward smaller k. Such inverse cascade is well-known in 2D hydrodynamic turbulence such as a typhoon or hurricane formation and is known to have Kolmogorov spectrum k−5/3. Although only a small difference in power, the physics of inverse cascades is very different as shown in the present paper. With use of Alfvén current limit condition, the criteria of collisionless shock formation are evaluated. The present theory is compared to corresponding experiments done with Omega and NIF lasers and a variety of PIC simulations. The theory is also applied to evaluate the strength of magnetic field near the shock front of the supernova remnant SN1006. The enhancement of magnetic field of about 25 μG is concluded in the present theory. Finally, a universality of the model equation is shown by applying the theory to the turbulent mixing due to Rayleigh–Taylor instability at the contact surface of two fluids in a gravitational or inertial force, which is very important in compressing plasma such as inertial confinement fusion by implosion. It is shown that the well-known evolution physics, mixing layer of the two fluids grows in proportion to (time)2, can be explained with the same model equation.

Stratospheric vacillations in a minimal contour dynamics model

AbstractA low-dimensional dynamical system that describes dynamical variability of the stratospheric polar vortex is presented. The derivation is based on a linearized, contour-dynamics representation of quasigeostrophic shallow water flow on a polar f-plane. The model consists of a single linear wave mode propagating on a near-circular patch of constant potential vorticity (PV). The PV jump at the vortex edge serves as an additional degree of freedom. The wave is forced by surface topography, and interacts with the vortex through a simplified parameterization of diabatic wave/mean flow interaction. The approach can be generalized to other geometries.The resulting three-component system depends on four non-dimensional parameters, and the structure of the steady state solutions can be determined analytically in some detail. Despite its extreme simplification, the model exhibits variability that is closely analogous to the Holton-Mass model, a well-known and more complex dynamical model of stratospheric variability. The present model exhibits two stable steady solutions, one consisting of a strong vortex with a small amplitude wave and the second consisting of a weak vortex with a large amplitude wave. Periodic and aperiodic limit cycles are also identified, analogous to similar solutions in the Holton-Mass model. Model trajectories also exhibit a number of behaviors that have been identified in observations. A key insight is that the time-mean state of the vortex is predominantly controlled by the properties of the linear mode, while the strength of the topographic forcing plays a far weaker role away from bifurcations.

Legolas: A Modern Tool for Magnetohydrodynamic Spectroscopy

Abstract Magnetohydrodynamic (MHD) spectroscopy is central to many astrophysical disciplines, ranging from helio- to asteroseismology, over solar coronal (loop) seismology, to the study of waves and instabilities in jets, accretion disks, or solar/stellar atmospheres. MHD spectroscopy quantifies all linear (standing or traveling) wave modes, including overstable (i.e., growing) or damped modes, for a given configuration that achieves force and thermodynamic balance. Here, we present Legolas, a novel, open-source numerical code to calculate the full MHD spectrum of one-dimensional equilibria with flow, balancing pressure gradients, Lorentz forces, centrifugal effects, and gravity, and enriched with nonadiabatic aspects like radiative losses, thermal conduction, and resistivity. The governing equations use Fourier representations in the ignorable coordinates, and the set of linearized equations is discretized using finite elements in the important height or radial variation, handling Cartesian and cylindrical geometries using the same implementation. A weak Galerkin formulation results in a generalized (non-Hermitian) matrix eigenvalue problem, and linear algebraic algorithms calculate all eigenvalues and corresponding eigenvectors. We showcase a plethora of well-established results, ranging from p and g modes in magnetized, stratified atmospheres, over modes relevant for coronal loop seismology, thermal instabilities, and discrete overstable Alfvén modes related to solar prominences, to stability studies for astrophysical jet flows. We encounter (quasi-)Parker, (quasi-)interchange, current-driven, and Kelvin–Helmholtz instabilities, as well as nonideal quasi-modes, resistive tearing modes, up to magnetothermal instabilities. The use of high resolution sheds new light on previously calculated spectra, revealing interesting spectral regions that have yet to be investigated.

Mixed Rossby–Gravity Wave–Wave Interactions

AbstractMixed triad wave–wave interactions between Rossby and gravity waves are analytically derived using the kinetic equation for models of different complexity. Two examples are considered: initially vanishing linear gravity wave energy in the presence of a fully developed Rossby wave field and the reversed case of initially vanishing linear Rossby wave energy in the presence of a realistic gravity wave field. The kinetic equation in both cases is numerically evaluated, for which energy is conserved within numerical precision. The results are validated by a corresponding ensemble of numerical model simulations supporting the validity of the weak-interaction assumption necessary to derive the kinetic equation. Since they are generated by nonresonant interactions only, the energy transfers toward the respective linear wave mode with vanishing energy are small in both cases. The total generation of energy of the linear gravity wave mode in the first case scales to leading order as the square of the Rossby number in agreement with independent estimates from laboratory experiments, although a part of the linear gravity wave mode is slaved to the Rossby wave mode without wavelike temporal behavior.

Incompressible magnetohydrodynamic modes in the thin magnetically twisted flux tube

Context. Observations have shown that twisted magnetic fields naturally occur, and indeed are omnipresent in the Sun's atmosphere. It is therefore of great theoretical interest in solar atmospheric waves research to investigate the types of magnetohydrodynamic (MHD) wave modes that can propagate along twisted magnetic flux tubes. Aims. Within the framework of ideal MHD, the main aim of this work is to investigate small amplitude incompressible wave modes of twisted magnetic flux tubes with m more or equal 1. The axial magnetic field strength inside and outside the tube will be allowed to vary, to ensure the results will not be restricted to only cold plasma equilibria conditions. Methods. The dispersion equation for these incompressible linear MHD wave modes was derived analytically by implementing the long wavelength approximation. Results. It is shown, in the long wavelength limit, that both the frequency and radial velocity profile of the m = 1 kink mode are completely unaffected by the choice of internal background magnetic twist. However, fluting modes with m more or equal 2 are sensitive to the particular radial profile of magnetic twist chosen. Furthermore, due to background twist, a low frequency cut-off is introduced for fluting modes that is not present for kink modes. From an observational point of view, although magnetic twist does not affect the propagation of long wavelength kink modes, for fluting modes it will either work for or against the propagation, depending on the direction of wave travel relative to the sign of the background twist.

Study of turbulence and interacting inertial modes in a differentially rotating spherical shell experiment

We present a study of inertial modes in a differentially rotating spherical shell (spherical Couette flow) experiment with a radius ratio of $\eta = 1/3$. Inertial modes are Coriolis-restored linear wave modes which often arise in rapidly rotating fluids. Recent experimental work has shown that inertial modes exist in a spherical Couette flow for $\Omega_{i}<\Omega_{o}$, where $\Omega_i$ and $\Omega_o$ is the inner and outer sphere rotation rate. A finite number of particular inertial modes has previously been found. By scanning the Rossby number from $-2.5 < Ro = (\Omega_{i}-\Omega_{o})/\Omega_{o} < 0$ at two fixed $\Omega_{o}$, we report the existence of similar inertial modes. However, the behavior of the flow described here differs much from previous spherical Couette experiments. We show that the kinetic energy of the dominant inertial mode dramatically increases with decreasing Rossby number that eventually leads to a wave-breaking and an increase of small-scale structures at a critical Rossby number. Such a transition in a spherical Couette flow has not been described before. The critical Rossby number scales with the Ekman number as0 $E^{1/5}$. Additionally, the increase of small-scale features beyond the transition transfers energy to a massively enhanced mean flow around the tangent cylinder. In this context, we discuss an interaction between the dominant inertial modes with a geostrophic Rossby mode exciting secondary modes whose frequencies match the triadic resonance condition.

Solitons and conservation laws of coupled Ostrovsky equation for internal waves

We study the soliton like behavior and conservation laws of a model that describes weakly nonlinear oceanic internal waves when two distinct linear long wave modes have nearly coincident phase speeds. The equations governing the model is the systems equivalent of the Ostrovsky equation.

Coupled Ostrovsky equations for internal waves in a shear flow

In the context of fluid flows, the coupled Ostrovsky equations arise when two distinct linear long wave modes have nearly coincident phase speeds in the presence of background rotation. In this paper, nonlinear waves in a stratified fluid in the presence of shear flow are investigated both analytically, using techniques from asymptotic perturbation theory, and through numerical simulations. The dispersion relation of the system, based on a three-layer model of a stratified shear flow, reveals various dynamical behaviours, including the existence of unsteady and steady envelope wave packets.

VALIDITY OF THE TAYLOR HYPOTHESIS FOR LINEAR KINETIC WAVES IN THE WEAKLY COLLISIONAL SOLAR WIND

The interpretation of single-point spacecraft measurements of solar wind turbulence is complicated by the fact that the measurements are made in a frame of reference in relative motion with respect to the turbulent plasma. The Taylor hypothesis---that temporal fluctuations measured by a stationary probe in a rapidly flowing fluid are dominated by the advection of spatial structures in the fluid rest frame---is often assumed to simplify the analysis. But measurements of turbulence in upcoming missions, such as Solar Probe Plus, threaten to violate the Taylor hypothesis, either due to slow flow of the plasma with respect to the spacecraft or to the dispersive nature of the plasma fluctuations at small scales. Assuming that the frequency of the turbulent fluctuations is characterized by the frequency of the linear waves supported by the plasma, we evaluate the validity of the Taylor hypothesis for the linear kinetic wave modes in the weakly collisional solar wind. The analysis predicts that a dissipation range of solar wind turbulence supported by whistler waves is likely to violate the Taylor hypothesis, while one supported by kinetic Alfven waves is not.

Coupled Ostrovsky Equations for Internal Waves, with a Background Shear Flow

We study the behaviour of weakly nonlinear oceanic internal waves in the presence of background rotation and shear flow, when two distinct linear long wave modes have nearly coincident phase speeds. The waves are described by a system of coupled Ostrovsky equations, derived from the full set of Euler equations for incompressible density stratified fluid with a free surface and rigid bottom boundary conditions. We report here some preliminary results obtained when there is a background shear flow.

On strongly interacting internal waves in a rotating ocean and coupled Ostrovsky equations

In the weakly nonlinear limit, oceanic internal solitary waves for a single linear long wave mode are described by the KdV equation, extended to the Ostrovsky equation in the presence of background rotation. In this paper we consider the scenario when two different linear long wave modes have nearly coincident phase speeds and show that the appropriate model is a system of two coupled Ostrovsky equations. These are systematically derived for a density-stratified ocean. Some preliminary numerical simulations are reported which show that, in the generic case, initial solitary-like waves are destroyed and replaced by two coupled nonlinear wave packets, being the counterpart of the same phenomenon in the single Ostrovsky equation.

USING SYNTHETIC SPACECRAFT DATA TO INTERPRET COMPRESSIBLE FLUCTUATIONS IN SOLAR WIND TURBULENCE

Kinetic plasma theory is used to generate synthetic spacecraft data to analyze and interpret the compressible fluctuations in the inertial range of solar wind turbulence. The kinetic counterparts of the three familiar linear MHD wave modes---the fast, Alfven, and slow waves---are identified and the properties of the density-parallel magnetic field correlation for these kinetic wave modes is presented. The construction of synthetic spacecraft data, based on the quasi-linear premise---that some characteristics of magnetized plasma turbulence can be usefully modeled as a collection of randomly phased, linear wave modes---is described in detail. Theoretical predictions of the density-parallel magnetic field correlation based on MHD and Vlasov-Maxwell linear eigenfunctions are presented and compared to the observational determination of this correlation based on 10 years of Wind spacecraft data. It is demonstrated that MHD theory is inadequate to describe the compressible turbulent fluctuations and that the observed density-parallel magnetic field correlation is consistent with a statistically negligible kinetic fast wave energy contribution for the large sample used in this study. A model of the solar wind inertial range fluctuations is proposed comprised of a mixture of a critically balanced distribution of incompressible Alfvenic fluctuations and a critically balanced or more anisotropic than critical balance distribution of compressible slow wave fluctuations. These results imply that there is little or no transfer of large scale turbulent energy through the inertial range down to whistler waves at small scales.

Wave Instabilities of a Collisionless Plasma in Fluid Approximation

AbstractWave properties and instabilities in a magnetized, anisotropic, collisionless, rarefied hot plasma in fluid approx‐imation are studied, using the 16‐moments set of the transport equations obtained from the Vlasov equations. These equations differ from the CGL‐MHD fluid model (single fluid equations by Chew, Goldberger, and Low [5,9]) by including two anisotropic heat flux evolution equations, where the fluxes invalidate the double polytropic CGL laws. We derived the general dispersion relation for linear compressible wave modes. Besides the classic incompressible fire hose modes there appear four types of compressible wave modes: two fast and slow mirror modes – strongly modified compared to the CGL model – and two thermal modes. In the presence of initial heat fluxes along the magnetic field the wave properties become different for the waves running forward and backward with respect to the magnetic field. The well known discrepancies between the results of the CGL‐MHD fluid model and the kinetic theory are now removed: i) The mirror slow mode instability criterion is now the same as that in the kinetic theory. ii) Similarly, in kinetic studies there appear two kinds of fire hose instabilities ‐ incompressible and compressible ones. These two instabilities can arise for the same plasma parameters, and the instability of the new compressible oblique fire hose modes can become dominant. The compressible fire hose instability is the result of the resonance coupling of three retrograde modes ‐ two thermal modes and a fast mirror mode. The results can be applied to the theory of solar and stellar coronal and wind models (© 2011 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

Selection of inertial modes in spherical Couette flow

Spherical Couette flow involves fluid sheared between concentric coaxially rotating spheres. Its scientific relevance lies not only in the simplicity of the system but also in its applicability to astrophysical objects such as atmospheres, oceans, and planetary cores. One common behavior in all rotating flows, including spherical Couette flow, is the presence of inertial modes, which are linear wave modes restored by the Coriolis force. Building on a previous identification of inertial modes in a laboratory spherical Couette cell, here we propose selection mechanisms to explain the presence of the particular modes we have observed. Mode selection depends on both amplification and damping. Our experimental observations are consistent with amplification and selection by over-reflection at a shear layer, and we would expect other spherical Couette devices to behave similarly. Damping effects, due in part to the presence of an inner sphere, add further constraints which are likely to play a role in mode selection in planetary atmospheres and cores, including the core of earth.

X and Y waves in the spatiotemporal Kerr dynamics of a self-guided light beam

The nonlinear stage of development of the spatiotemporal instability of the monochromatic Townes beam in a medium with self-focusing nonlinearity and normal dispersion is studied by analytical and numerical means. Small perturbations to the self-guided light beam are found to grow into two giant, splitting Y pulses featuring shock fronts on opposite sides. Each shocking pulse amplifies a co-propagating X wave, or dispersion- and diffraction-free linear wave mode of the medium, with super-broad spectrum.

Unstable magnetohydrodynamical continuous spectrum of accretion disks

We present a detailed study of localised magnetohydrodynamical (MHD) instabilities occuring in two--dimensional magnetized accretion disks. We model axisymmetric MHD disk tori, and solve the equations governing a two--dimensional magnetized accretion disk equilibrium and linear wave modes about this equilibrium. We show the existence of novel MHD instabilities in these two--dimensional equilibria which do not occur in an accretion disk in the cylindrical limit. The disk equilibria are numerically computed by the FINESSE code. The stability of accretion disks is investigated analytically as well as numerically. We use the PHOENIX code to compute all the waves and instabilities accessible to the computed disk equilibrium. We concentrate on strongly magnetized disks and sub--Keplerian rotation in a large part of the disk. These disk equilibria show that the thermal pressure of the disk can only decrease outwards if there is a strong gravitational potential. Our theoretical stability analysis shows that convective continuum instabilities can only appear if the density contours coincide with the poloidal magnetic flux contours. Our numerical results confirm and complement this theoretical analysis. Furthermore, these results show that the influence of gravity can either be stabilizing or destabilizing on this new kind of MHD instability. In the likely case of a non--constant density, the height of the disk should exceed a threshold before this type of instability can play a role. This localised MHD instability provides an ideal, linear route to MHD turbulence in strongly magnetized accretion disk tori.

Green’s function of compressible Petschek-type magnetic reconnection

We present a method to analyze the wave and shock structures arising from Petschek-type magnetic reconnection. Based on a time-dependent analytical approach developed by Heyn and Semenov [Phys. Plasmas 3, 2725 (1996)] and Semenov et al. [Phys. Plasmas 11, 62 (2004)], we calculate the perturbations caused by a delta function-shaped reconnection electric field, which allows us to achieve a representation of the plasma variables in the form of Green’s functions. Different configurations for the initial conditions are considered. In the case of symmetric, antiparallel magnetic fields and symmetric plasma density, the well-known structure of an Alfvén discontinuity, a fast body wave, a slow shock, a slow wave, and a tube wave occurs. In the case of asymmetric, antiparallel magnetic fields, additionally surface waves are found. We also discuss the case of symmetric, antiparallel magnetic fields and asymmetric densities, which leads to a faster propagation in the lower half plane, causing side waves forming a Mach cone in the upper half plane. Complex effects like anisotropic propagation characteristics, intrinsic wave coupling, and the generation of different nonlinear and linear wave modes in a finite β plasma are retained. The temporal evolution of these wave and shock structures is shown.