## Abstract

The spatial structure of a steady state plasma flow is shaped by the standing modes with local phase velocity exactly opposite to the flow velocity. The general procedure of finding the wave vectors of all possible standing MHD modes in any given point of a stationary flow requires numerically solving an algebraic equation. We present the graphical procedure (already mentioned by some authors in the 1960’s) along with the exact solution for the Alfvén mode and approximate analytic solutions for both fast and slow modes. The technique can be used to identify MHD modes in space and laboratory plasmas as well as in numerical simulations.

This is a preview of subscription content, access via your institution.

## References

Krisko, P.H., Hill, T.W.: Two-dimensional model of a slow-mode expansion fan at Io. Geophys. Res. Lett.

**18**, 1947–1950 (1991). doi:10.1029/91GL02541Lacombe, C., Pantellini, F.G.E., Hubert, D., Harvey, C.C., Mangeney, A., Belmont, G., Russell, C.T.: Mirror and Alfvenic waves observed by ISEE 1-2 during crossings of the earth’s bow shock. Ann. Geophys.

**10**, 772–784 (1992)Pantellini, F., Griton, L., Varela, J.: Rarefaction and compressional standing slow mode structures in Mercury’s magnetosheath: 3D MHD simulations. Planet. Space Sci.

**112**, 1–9 (2015). doi:10.1016/j.pss.2015.04.007Siscoe, G.L., Sanchez, E.: An MHD model for the complete open magnetotail boundary. J. Geophys. Res.

**92**, 7405–7412 (1987). doi:10.1029/JA092iA07p07405Spreiter, J.R., Alksne, A.Y.: Solar-wind flow past objects in the solar system. Annu. Rev. Fluid Mech.

**2**, 313–354 (1970). doi:10.1146/annurev.fl.02.010170.001525Verigin, M., Slavin, J., Szabo, A., Kotova, G., Gombosi, T.: Planetary bow shocks: Asymptotic MHD Mach cones. Earth Planets Space

**55**, 33–38 (2003)

## Acknowledgements

The work of Léa Griton was supported by the Centre National d’Études Spatiales (CNES) and the Observatoire de Paris (contract ref. 5100016058). Part of the presented work has been completed by Laura Munteanu during a one month internship at LESIA in June 2015.

## Author information

### Affiliations

### Corresponding author

## Appendices

### Appendix A: Changing reference frame

The calculations presented in this paper have been established in a particular reference frame (hereafter the primed reference frame) defined by the flow velocity vector \(\vec{u}_{f}\) and the magnetic field vector \(\vec{B}\) as illustrated in Fig. 1. This particular frame is well suited for the theoretical treatment of mode propagation related issues but is rarely the most practical one as both velocity and magnetic field orientations do generally change as a function of space. In most applications, however, there is one unique frame for the whole system (for example the simulation frame as in Fig. 14 or the frame defined by the step geometry in the Krisko and Hill 1991 experiment of Fig. 13). We denote this frame as the unprimed frame. It is therefore useful to establish the matrix \(M\) which transforms the orthonormal right handed basis vectors of an arbitrary frame into the orthonormal basis vectors (also right handed) of the primed reference frame. The plasma flow being stationary in both frames there is no relative motion and the transformation matrix \(M\) may be viewed as the product of two rotations implying \(\det (M)=1\). \(M\) is easily obtained by writing the basis vectors \(\vec{e}_{x}\,'\), \(\vec{e}_{y}\,'\) and \(\vec{e}_{z}\,'\) of the primed reference frame in terms of the two unitary vectors \(\vec{\mu }_{f}\equiv \vec{u}_{f}/u _{f}\) and \(\vec{b}\equiv \vec{B}/B\) in the unprimed frame, i.e.

where \(s = \mathrm{sign} (\vec{\mu }_{f}\cdot \vec{b})\) and \(0 < \Delta _{f}\leq \pi /2\) is the angle between the directions of \(\vec{b}\) and \(\vec{\mu }_{f}\). In the primed frame the magnetic field is therefore either parallel or anti-parallel with respect to \(\vec{e}_{z}\,'\). If \(\Delta_{f}\neq 0\), the components of the three primed basis vectors are fully specified by the constraint that \(\vec{\mu }_{f}\) is in the \((\vec{e}_{x}\,',\vec{e}_{z}\,')\) plane with negative components \(\mu_{fx}\,'=\vec{\mu }_{f}\cdot \vec{e}_{x}\,'<0\) and \(\mu_{fz}\,'= \vec{\mu }_{f}\cdot \vec{e}_{z}\,'<0\). In the singular case \(\Delta_{f}=0\), \(\vec{e}_{y}\,'\) may be any unitary vector perpendicular to \(\vec{b}\) and \(\vec{e}_{x}\,' =\vec{e}_{y}\,'\times \vec{e}_{z}\,'\) as in the general case (A.1).

The transformation of the basis vectors \(\vec{e}_{x} =(1,0,0)\), \(\vec{e}_{y} =(0,1,0)\) and \(\vec{e}_{z} =(0,0,1)\) to the primed frame is thus given by

where the elements of the \(3\times 3\) transformation matrix \(M\) are merely the components of the primed basis vectors:

Accordingly, the components of an arbitrary vector projected onto the basis vectors of the primed frame are obtained by applying the transposed matrix \(M^{T}\) to the components of the vector in the unprimed frame, i.e.

where \(\vec{v}\,'=(\vec{v}\cdot \vec{e}_{x}\,',\vec{v}\cdot \vec{e} _{y}\,',\vec{v}\cdot \vec{e}_{z}\,')\).

### A.1 Special case: \(\mu_{fy} = b_{y} = 0\)

The transformation matrix \(M\) is particularly simple in the case where the \(y\) components of \(\vec{\mu }_{f}\) and \(\vec{b}\) are zero in the unprimed frame. In this case the basis vectors (A.1) reduce to:

where \(\sigma \equiv \mathrm{sign} (-\mu_{fx} b_{z} + \mu_{fz} b_{x})\). Accordingly, the transformation matrix \(M\) for this particular case is

### Appendix B: Parallel compressibility and Alfvén ratio

### B.1 The parallel compressibility

The parallel compressibility of a plane mode with wave vector \(\vec{k}\) is defined as

where \(\delta n\) and \(\delta B_{\parallel }\) are the variations of density and magnetic field along a path parallel to \(\vec{k}\). The subscript ∥ in \(\delta B_{\parallel }\) denotes the variation of the magnetic field parallel to itself so that (for example) \(\delta B_{\parallel }=0\) through a rotational discontinuity. Assuming a wave vector \(\vec{k}=(k,0,0)\):

and

In case of an arbitrary orientation \(\vec{k}\) the variations \(\delta B_{\parallel }\) and \(\delta n\) over an infinitesimal distance \(\vec{\delta x}=\epsilon \vec{k}\) are

and

respectively. The parallel compressibility can then be written as:

For the incompressible Alfvén mode the parallel compressibility is zero. For both the slow and the fast mode the parallel compressibility is given by

where \(\theta \) is the angle between \(\vec{k}\) and the magnetic field \(\vec{B}\), \(c_{A}\) is the Alfvén speed, \(c\) the adiabatic sound speed and \(u_{\phi }=\omega /k\) the phase velocity (3) of the corresponding mode. We note that for non-zero values of the Alfvén speed, the fast mode does always propagate faster than the sound speed. Thus, according to (B.7) the compressibility of the fast mode is always positive. On the contrary, slow modes do always propagate slower than the sound speed (except for the special case \(c=c_{A}\) and \(\theta =0\)). Thus, the denominator in (B.7) is always negative implying a negative compressibility for the slow mode. A typical example of parallel compressibility profiles for both compressible MHD modes is shown in Fig. 16.

### B.2 The Alfvén ratio

The Alfvén ratio is defined as

where the perpendicular direction is now to be considered with respect to \(\vec{k}\), i.e. \(\delta v_{\perp ,k}=|\vec{\delta v}\times \vec{k}/k|\).

For the three MHD modes the Alfvén ratio is given by

The phase velocity of the Alfvén mode being \(u_{\phi }^{2}=c_{A} ^{2} \cos^{2}\theta \) implies that \(R_{A}=1\), independently of the propagation angle \(\theta \). Sample profiles of the Alfvén ratio for the three MHD modes are shown in Fig. 16.

## Rights and permissions

## About this article

### Cite this article

Pantellini, F., Griton, L. Identification of standing fronts in steady state fluid flows: exact and approximate solutions for propagating MHD modes.
*Astrophys Space Sci* **361, **335 (2016). https://doi.org/10.1007/s10509-016-2921-y

Received:

Accepted:

Published:

### Keywords

- Magnetohydrodynamics (MHD)
- Waves
- Planets and satellites: magnetic fields
- Methods: analytical