Magnetohydrodynamic turbulence concerns the chaotic regimes of magnetofluid flow at high Reynolds number. Magnetohydrodynamics (MHD) deals with what is a quasi-neutral fluid with very high conductivity. The fluid approximation implies that the focus is on macro length-and-time scales which are much larger than the collision length and collision time respectively.
The incompressible MHD equations for constant mass density, , are
where
The third equation is the incompressibility condition. In the above equation, the magnetic field is in Alfvén units (same as velocity units).
The total magnetic field can be split into two parts: (mean + fluctuations).
The above equations in terms of Elsässer variables () are
where . Nonlinear interactions occur between the Alfvénic fluctuations .
The important nondimensional parameters for MHD are
The magnetic Prandtl number is an important property of the fluid. Liquid metals have small magnetic Prandtl numbers, for example, liquid sodium's is around . But plasmas have large .
The Reynolds number is the ratio of the nonlinear term of the Navier–Stokes equation to the viscous term. While the magnetic Reynolds number is the ratio of the nonlinear term and the diffusive term of the induction equation.
In many practical situations, the Reynolds number of the flow is quite large. For such flows typically the velocity and the magnetic fields are random. Such flows are called to exhibit MHD turbulence. Note that need not be large for MHD turbulence. plays an important role in dynamo (magnetic field generation) problem.
The mean magnetic field plays an important role in MHD turbulence, for example it can make the turbulence anisotropic; suppress the turbulence by decreasing energy cascade etc. The earlier MHD turbulence models assumed isotropy of turbulence, while the later models have studied anisotropic aspects. In the following discussions will summarize these models. More discussions on MHD turbulence can be found in Biskamp, [1] Verma. [2] and Galtier.
Iroshnikov [3] and Kraichnan [4] formulated the first phenomenological theory of MHD turbulence. They argued that in the presence of a strong mean magnetic field, and wavepackets travel in opposite directions with the phase velocity of , and interact weakly. The relevant time scale is Alfven time . As a results the energy spectra is
where is the energy cascade rate.
Later Dobrowolny et al. [5] derived the following generalized formulas for the cascade rates of variables:
where are the interaction time scales of variables.
Iroshnikov and Kraichnan's phenomenology follows once we choose .
Marsch [6] chose the nonlinear time scale as the interaction time scale for the eddies and derived Kolmogorov-like energy spectrum for the Elsasser variables:
where and are the energy cascade rates of and respectively, and are constants.
Matthaeus and Zhou [7] attempted to combine the above two time scales by postulating the interaction time to be the harmonic mean of Alfven time and nonlinear time.
The main difference between the two competing phenomenologies (−3/2 and −5/3) is the chosen time scales for the interaction time. The main underlying assumption in that Iroshnikov and Kraichnan's phenomenology should work for strong mean magnetic field, whereas Marsh's phenomenology should work when the fluctuations dominate the mean magnetic field (strong turbulence).
However, as we will discuss below, the solar wind observations and numerical simulations tend to favour −5/3 energy spectrum even when the mean magnetic field is stronger compared to the fluctuations. This issue was resolved by Verma [8] using renormalization group analysis by showing that the Alfvénic fluctuations are affected by scale-dependent "local mean magnetic field". The local mean magnetic field scales as , substitution of which in Dobrowolny's equation yields Kolmogorov's energy spectrum for MHD turbulence.
Renormalization group analysis have been also performed for computing the renormalized viscosity and resistivity. It was shown that these diffusive quantities scale as that again yields energy spectra consistent with Kolmogorov-like model for MHD turbulence. The above renormalization group calculation has been performed for both zero and nonzero cross helicity.
The above phenomenologies assume isotropic turbulence that is not the case in the presence of a mean magnetic field. The mean magnetic field typically suppresses the energy cascade along the direction of the mean magnetic field. [9]
Mean magnetic field makes turbulence anisotropic. This aspect has been studied in last two decades. In the limit , Galtier et al. [10] showed using kinetic equations that
where and are components of the wavenumber parallel and perpendicular to mean magnetic field. The above limit is called the weak turbulence limit.
Under the strong turbulence limit, , Goldereich and Sridhar [11] argue that ("critical balanced state") which implies that
The above anisotropic turbulence phenomenology has been extended for large cross helicity MHD.
Solar wind plasma is in a turbulent state. Researchers have calculated the energy spectra of the solar wind plasma from the data collected from the spacecraft. The kinetic and magnetic energy spectra, as well as are closer to compared to , thus favoring Kolmogorov-like phenomenology for MHD turbulence. [12] [13] The interplanetary and interstellar electron density fluctuations also provide a window for investigating MHD turbulence.
The theoretical models discussed above are tested using the high resolution direct numerical simulation (DNS). Number of recent simulations report the spectral indices to be closer to 5/3. [14] There are others that report the spectral indices near 3/2. The regime of power law is typically less than a decade. Since 5/3 and 3/2 are quite close numerically, it is quite difficult to ascertain the validity of MHD turbulence models from the energy spectra.
Energy fluxes can be more reliable quantities to validate MHD turbulence models. When (high cross helicity fluid or imbalanced MHD) the energy flux predictions of Kraichnan and Iroshnikov model is very different from that of Kolmogorov-like model. It has been shown using DNS that the fluxes computed from the numerical simulations are in better agreement with Kolmogorov-like model compared to Kraichnan and Iroshnikov model. [15]
Anisotropic aspects of MHD turbulence have also been studied using numerical simulations. The predictions of Goldreich and Sridhar [11] () have been verified in many simulations.
Energy transfer among various scales between the velocity and magnetic field is an important problem in MHD turbulence. These quantities have been computed both theoretically and numerically. [2] These calculations show a significant energy transfer from the large scale velocity field to the large scale magnetic field. Also, the cascade of magnetic energy is typically forward. These results have critical bearing on dynamo problem.
There are many open challenges in this field that hopefully will be resolved in near future with the help of numerical simulations, theoretical modelling, experiments, and observations (e.g., solar wind).
In physics and engineering, magnetohydrodynamics is a model of electrically conducting fluids that treats all interpenetrating particle species together as a single continuous medium. It is primarily concerned with the low-frequency, large-scale, magnetic behavior in plasmas and liquid metals and has applications in multiple fields including space physics, geophysics, astrophysics, and engineering.
The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).
In physics, the screened Poisson equation is a Poisson equation, which arises in the Klein–Gordon equation, electric field screening in plasmas, and nonlocal granular fluidity in granular flow.
Poisson's equation is an elliptic partial differential equation of broad utility in theoretical physics. For example, the solution to Poisson's equation is the potential field caused by a given electric charge or mass density distribution; with the potential field known, one can then calculate the corresponding electrostatic or gravitational (force) field. It is a generalization of Laplace's equation, which is also frequently seen in physics. The equation is named after French mathematician and physicist Siméon Denis Poisson who published it in 1823.
In physics, the dynamo theory proposes a mechanism by which a celestial body such as Earth or a star generates a magnetic field. The dynamo theory describes the process through which a rotating, convecting, and electrically conducting fluid can maintain a magnetic field over astronomical time scales. A dynamo is thought to be the source of the Earth's magnetic field and the magnetic fields of Mercury and the Jovian planets.
Quantum turbulence is the name given to the turbulent flow – the chaotic motion of a fluid at high flow rates – of quantum fluids, such as superfluids. The idea that a form of turbulence might be possible in a superfluid via the quantized vortex lines was first suggested by Richard Feynman. The dynamics of quantum fluids are governed by quantum mechanics, rather than classical physics which govern classical (ordinary) fluids. Some examples of quantum fluids include superfluid helium, Bose–Einstein condensates (BECs), polariton condensates, and nuclear pasta theorized to exist inside neutron stars. Quantum fluids exist at temperatures below the critical temperature at which Bose-Einstein condensation takes place.
In physics and mathematics, the Helmholtz decomposition theorem or the fundamental theorem of vector calculus states that certain differentiable vector fields can be resolved into the sum of an irrotational (curl-free) vector field and a solenoidal (divergence-free) vector field. In physics, often only the decomposition of sufficiently smooth, rapidly decaying vector fields in three dimensions is discussed. It is named after Hermann von Helmholtz.
In plasma physics, magnetic helicity is a measure of the linkage, twist, and writhe of a magnetic field.
In classical electromagnetism, magnetic vector potential is the vector quantity defined so that its curl is equal to the magnetic field: . Together with the electric potential φ, the magnetic vector potential can be used to specify the electric field E as well. Therefore, many equations of electromagnetism can be written either in terms of the fields E and B, or equivalently in terms of the potentials φ and A. In more advanced theories such as quantum mechanics, most equations use potentials rather than fields.
In physics, magnetosonic waves, also known as magnetoacoustic waves, are low-frequency compressive waves driven by mutual interaction between an electrically conducting fluid and a magnetic field. They are associated with compression and rarefaction of both the fluid and the magnetic field, as well as with an effective tension that acts to straighten bent magnetic field lines. The properties of magnetosonic waves are highly dependent on the angle between the wavevector and the equilibrium magnetic field and on the relative importance of fluid and magnetic processes in the medium. They only propagate with frequencies much smaller than the ion cyclotron or ion plasma frequencies of the medium, and they are nondispersive at small amplitudes.
In plasma physics, the Vlasov equation is a differential equation describing time evolution of the distribution function of collisionless plasma consisting of charged particles with long-range interaction, such as the Coulomb interaction. The equation was first suggested for the description of plasma by Anatoly Vlasov in 1938 and later discussed by him in detail in a monograph. The Vlasov equation, combined with Landau kinetic equation describe collisional plasma.
The Grad–Shafranov equation is the equilibrium equation in ideal magnetohydrodynamics (MHD) for a two dimensional plasma, for example the axisymmetric toroidal plasma in a tokamak. This equation takes the same form as the Hicks equation from fluid dynamics. This equation is a two-dimensional, nonlinear, elliptic partial differential equation obtained from the reduction of the ideal MHD equations to two dimensions, often for the case of toroidal axisymmetry. Taking as the cylindrical coordinates, the flux function is governed by the equation,
In relativistic physics, the electromagnetic stress–energy tensor is the contribution to the stress–energy tensor due to the electromagnetic field. The stress–energy tensor describes the flow of energy and momentum in spacetime. The electromagnetic stress–energy tensor contains the negative of the classical Maxwell stress tensor that governs the electromagnetic interactions.
In plasma physics, the Hasegawa–Mima equation, named after Akira Hasegawa and Kunioki Mima, is an equation that describes a certain regime of plasma, where the time scales are very fast, and the distance scale in the direction of the magnetic field is long. In particular the equation is useful for describing turbulence in some tokamaks. The equation was introduced in Hasegawa and Mima's paper submitted in 1977 to Physics of Fluids, where they compared it to the results of the ATC tokamak.
In quantum mechanics, the Pauli equation or Schrödinger–Pauli equation is the formulation of the Schrödinger equation for spin-1/2 particles, which takes into account the interaction of the particle's spin with an external electromagnetic field. It is the non-relativistic limit of the Dirac equation and can be used where particles are moving at speeds much less than the speed of light, so that relativistic effects can be neglected. It was formulated by Wolfgang Pauli in 1927. In its linearized form it is known as Lévy-Leblond equation.
The Chandrasekhar number is a dimensionless quantity used in magnetic convection to represent ratio of the Lorentz force to the viscosity. It is named after the Indian astrophysicist Subrahmanyan Chandrasekhar.
In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.
The distortion free energy density is a quantity that describes the increase in the free energy density of a liquid crystal caused by distortions from its uniformly aligned configuration. It also commonly goes by the name Frank free energy density named after Frederick Charles Frank.
The spin angular momentum of light (SAM) is the component of angular momentum of light that is associated with the quantum spin and the rotation between the polarization degrees of freedom of the photon.
Chandrasekhar–Kendall functions are the eigenfunctions of the curl operator derived by Subrahmanyan Chandrasekhar and P. C. Kendall in 1957 while attempting to solve the force-free magnetic fields. The functions were independently derived by both, and the two decided to publish their findings in the same paper.