Nicholas Hall University of Oklahoma nicholashall@ou.edu REU Program - Summer 2006 University of Wisconsin - Madison Madison, WI 53706 Research projects of other REU students

# A Comparative Study of Velocity Statistics of Hydrodynamic and Magnetohydrodynamic Turbulence

## Abstract

Turbulence in an incompressible fluid with and without a magnetic field as well as moderately compressible MHD turbulence are compared. The results of three numerical simulation models in three dimensions of resolution up to $512^{3}$ are used for this purpose. The compatibility of the spectra of all three models with the Kolmogorov spectrum is confirmed. For the magnetohydrodynamic (MHD) models the probability distribution functions of the velocity components perpendicular to the external magnetic field are like the incompressible hydrodynamic (HD) model while those parallel to the field have a smaller range of velocities. The probability distribution functions of the transverse velocity increments for the MHD models decline slower than the incompressible HD model. The similarity of incompressible HD and both incompressible and compressible MHD turbulence persists over high order longitudinal structure function scaling exponents measured in the global reference frame as well as for motions perpendicular to the local mean field. In these two frames the longitudinal scaling exponents of both MHD models seem to follow theoretical incompressible HD dissipation structure predictions while the transverse scaling exponents of the incompressible MHD model seem to follow the predictions for incompressible MHD. In the local magnetic system the motions parallel to the local mean field for both MHD models are different from incompressible HD motions. The fields of the MHD simulations are decomposed into Alfvenic, fast, and slow modes. In the global reference frame and for motions perpendicular to the local mean field the Alfvenic mode is mostly responsible for the fact that the longitudinal components of both MHD models follow the incompressible HD model and the incompressible HD theoretical dissipation structure predictions. In the global reference frame for the incompressible MHD model both the Alfvenic and slow modes seem to contribute significantly to the fact that the transverse components follow the incompressible MHD theoretical predictions. For motions perpendicular to the magnetic field for the incompressible MHD model it appears that the slow mode contributes most to the fact that the transverse component seems to follow the incompressible MHD predictions.

## Introduction

Given that turbulent motions (velocities) are inherently stochastic, velocity field statistics are an essential tool for extracting its physical properties. In 1941 Kolmogorov published his theory (hereafter K41) on the local velocity field structure of turbulence for incompressible hydrodynamic (hereafter IHD) fluids with high Reynolds numbers (see kol41), in which he hypothesized the existence of homogeneity, isotropy, and an inertial range (hereafter K41). The inertial range in the kinetic energy spectrum is a range in wave number, $k\sim1/l$ where $l$ is the separation length between two points often regarded as the eddie size, which lies between the typical large scales and the dissipative scales and is therefore independent of external input such as turbulence drive or output such as dissipation (see bis03). In this range the spectrum is also independent of viscosity and exhibits a power-law dependence of the form E(k)=K$\bar{\epsilon}^{-2/3}$$k^{-5/3}$, where K is the universal Kolmogorov constant and $\bar{\epsilon}$ is the average rate of energy dissipation per unit mass.
It is known that the inclusion of an external magnetic field in the turbulence of an incompressible fluid brings about many substantial complications and can drastically modify the turbulent motions. The inclusion of compressibility also causes significant changes in the properties of turbulence.

got02 thoroughly and successfully investigated the velocity field statistics of high-resolution direct numerical simulations for IHD turbulence. The purpose of this paper is to address the question, are the velocity field statistics completely different from IHD when a magnetic field and compressibility are included? We extend the analysis of got02 to incompressible magnetohydrodynamic (hereafter IMHD) and compressible magnetohydrodynamic (hereafter CMHD) turbulence in an attempt to answer this question. This is considered mostly through the use of high order scaling exponents of the velocity structure functions.
Emphasis is also placed on the differences between the scaling exponents of longitudinal and transverse structure functions. For high order scaling exponents the findings of got02 are confirmed for structure functions measured in the global reference frame for IHD turbulence. Scaling exponents are then explored for IMHD and CMHD turbulence in the global frame and the local frame, defined by the local mean magnetic field. The scaling exponents are then decomposed into Alfvenic, fast, and slow modes and the contributions of individual modes are explored. The paper is organized as follows. Sec. II presents the numerical aspects of the three numerical simulation models considered (IHD, IMHD, and CMHD). Sec. III contains analysis of the kinetic energy spectra of the three models. In Sec. IV probability distribution functions are considered. Structure functions are explored in Sec. V. Sec. VI presents scaling exponents of the structure functions. In Sec. VII the scaling exponents are decomposed into the three MHD modes. Finally, Sec. VIII includes a discussion and summary.
The IHD data used in this study came from simulations with a grid size of $256^3$. For both MHD models, a high resolution grid size of $512^3$ was used and the external magnetic field is in the positive x-direction.

## Kinetic Energy Spectrum

Figure 1 - Kinetic energy spectra for the three models as well as Kolmogorov's power-law dependence in the inertial range for comparison.

The three-dimensional kinetic energy spectrum for each model is shown in figure 1. The inertial range of the spectrum, the flat part of the curve, is of most interest for this study. The shape of the curves at wave numbers less than those of the inertial range depends on how the turbulent motions are excited. The shape of the curves at wave numbers greater than those of the inertial range is affected significantly by viscosity, the bottleneck effect (fal94), and the resolution of the model. Despite the fact that the Kolmogorov phenomenology, E(k) proportional to k^{-5/3}, is derived for IHD, in the inertial range all three curves demonstrate very similar slopes, suggesting that the energy transfer process in the inertial range is consistent for all three models. This may be surprising considering the many changes resulting from the inclusion of a magnetic field, such as Alfvenic waves propagating parallel to the magnetic field and elongation of eddies in the velocity field. The inclusion of compressibility introduces significant changes as well, such as shocks. gol95 provide an explanation of this behavior. See cl05 also for a review. Deviations from the 5/3 power law have become an issue of debate (see bol06, ber06, and references therein). We do not address this issue.

## Probability Distribution Functions

Figure 2 - PDF's of the velocity field parallel (x) and perpendicular (y \& z) to $B_{ext}$ for the three models.

The probability distribution function (hereafter PDF) is used to obtain information about one-point velocity statistics. The single-point PDF's of different components of the velocity field for each model shown in figure 2 are histograms, normalized to the grid size, of the values of the velocity at each point. For the IHD model the x, y, and z components were averaged. For the two MHD models, the histograms were separated into components parallel (x direction) and perpendicular (y \& z directions) to the external magnetic field, ($B_{ext}$); therefore, the PDF's for the x components are presented separately and the y and z components are averaged.
The PDF for the IHD model is somewhat asymmetric. Any asymmetry or skewness in the PDF's should disappear when averaged over multiple time steps (see got02). The PDF's of the MHD simulations parallel to $B_{ext}$ are similar to each other, fairly symmetric, wider toward the top than a gaussian, and at large amplitudes decay faster than a gaussian. Their slopes are the steepest of the PDF's presented and given that the area under the curve must be the same, their ranges of velocities are smaller. For IMHD turbulence the divergence of the velocity vector is zero meaning that there are no sources, sinks, or shocks and that only eddies can exist. This is not the case for CMHD simulations in which shocks are produced at high velocities which dissipate the energy perhaps contributing to the PDF of the parallel CMHD model having the least high velocity points.
The PDF's of the IMHD and CMHD models perpendicular to $B_{ext}$ approximately follow each other and are even more similar to each other than the parallel components. They are very symmetric and gaussian-like, have the most shallow slopes of all, and therefore have the most high velocity points.
When an external magnetic field is introduced into the turbulent motions, an Alfvenic mode is produced with perturbations perpendicular to the direction of the field. This additional contribution to the velocity is responsible for the higher velocities found in the perpendicular PDF's. The additional kinetic energy given to the perpendicular direction explains the separation of these two components seen in the top of figure 2. This appears to be equally valid for both positive and negative values of the velocity.

Figure 3 - PDF's, normalized to match their maxima, of the velocity field parallel and perpendicular to $B_{ext}$ for the one time step of the CMHD simulation used throughout the paper and the average over 10 time steps of a 256 resolution simulation with error bars showing expected random fluctuations.

An average over 10 time steps of the velocity PDF's for a CMHD simulation with 256 resolution was performed. The 10 parallel time steps were averaged as well as the 20 perpendicular time steps, 10 for y and 10 for z. These PDF's are plotted in figure 3 along with the CMHD PDF's mentioned above. The PDF's are normalized to match their maxima in order to better compare their tails. The parallel and perpendicular PDF's for the single time step lie sufficiently close to the expected random errors obtained from the PDF's averaged over multiple time steps to justify using a single time step and averaging the perpendicular components.
 Figure 4 - PDF's of the transverse increments of the x component of the velocity field taken along y and z averaged together for the three models. Figure 5 - PDF's of the transverse increments of the y component of the velocity field taken along x and z presented separately for the three models.

Figures 4 and 5 show PDF's of the transverse increments of the velocity field which are histograms, normalized to the grid size, of the values of the difference between the velocity field at each point and the field at an adjacent point along the direction specified. Positive values suggests that $v$ is increasing in the positive direction while negative values suggest that it is decreasing. In figure 4 the increments of the x-component of the velocity field taken along the y and z directions, parallel to $B_{ext}$ for the MHD models, are averaged. Figure 5 shows the increments of the y-component along the x and z directions. The results are not averaged given that one is parallel to $B_{ext}$ while the other is perpendicular.
In both figures, all three models demonstrate near perfect symmetry. This is true despite the orientation with respect to $B_{ext}$. In figure 5 as expected for the IHD model, x and z are nearly identical because there is no magnetic field to bring about differences. On the other hand for the MHD models, x extends to higher values of the velocity increment than z.
Because of the fact that the difference in range between the x components for the IMHD and CMHD models is much greater than the difference in range between the z components for the IMHD and CMHD models, it is possible that most of the increase in range in going from IMHD models to CMHD models by introducing compressibility is coming from the components parallel to $B_{ext}$. This is supported by the fact that for the PDF's in figure 2 the perpendicular components follow each other closely while the parallel components are quite different from each other.

## Structure Functions

Structure functions are used to obtain information about two-point velocity statistics. Structure functions describe the distribution of velocity structures of size $l$ and are defined as $$S^{(p)}(l)=\langle | \delta v (\vec{l}) | ^p \rangle, \label{eqn:str_func}$$ where $\delta v$ is a velocity increment and $\langle ... \rangle$ represents ensemble averaging. Structure functions are defined with the absolute value of the velocity increment to ensure their positivity for all orders. Structure functions are separated into a longitudinal component naturally defined by the velocity increment $$\delta v_L(l) = [\vec{v}(\vec{x} + \vec{l}) - \vec{v} (\vec{x})]\cdot\hat{l}, \label{eqn:vel_incr}$$ and transverse components $\delta v_T (l)$ chosen randomly in the plane perpendicular to the longitudinal component. When higher order structure functions are used the PDF tails contribute more, where as you can see in figure 3 the errors are much larger.
For each separation length $l$ from 1 to 128, $\delta v$ is calculated at 100,000 points randomly distributed throughout the entire periodic velocity field. In the global reference frame the direction of $\hat{l}$, or the direction of the longitudinal structure function, is also random. The structure functions for the IHD model are calculated in the global reference frame. For the MHD models $\hat{l}$ can be chosen randomly (global frame), parallel to the local mean magnetic field, or perpendicular to the local mean field.

Figure 6 - For (a) points are chosen along the local mean field, hereafter referred to as local parallel. For (b) points are chosen in the plane perpendicular to the local mean field so that both the longitudinal and the transverse structure functions lie in this plane, hereafter referred to as local perpendicular.

According to cho00a and cho02, for MHD turbulence in a strongly magnetized medium eddies are not aligned along the direction of the total mean field but rather are aligned along the local mean field. This being the case, the local mean field is the physically relevant background for the eddie dynamics. The local frame is separated into two types which are represented in figure 6. In (a) the points from which the structure functions are calculated are chosen along the local mean field, this configuration will be referred to as local parallel. In (b) the points are chosen in the plane perpendicular to the local mean field, this configuration will be referred to as local perpendicular.

Figure 7 - Second order structure functions of the velocity field for the three models. The top row are longitudinal, the middle row are transverse, and the bottom row are total structure functions. The first column is in the global reference frame, the middle row is the local parallel frame, and the bottom row is the local perpendicular frame. In the local frames the HD structure functions from the global frame are shown for comparison.

Figure 7 shows the second order longitudinal, transverse, and total structure functions of the velocity field. The column on the left is the global frame, the middle column is the local parallel frame, and the column on the right is the local perpendicular frame. In the local frames the structure functions of the IHD model calculated in the global frame are shown for comparison.
The inertial range of the turbulent motions manifests itself in the structure functions as the range of $l$ over which the slope is constant in a log-log plot. The first maximum of a structure function is related to the size of the largest structures in the corresponding field, the radius of the largest eddies for example. It will be found at the separation length at which the largest $\delta v$ is observed.
We confirm the result of got02 that in the global frame for IHD turbulence the values of the transverse structure functions are greater than those of the longitudinal. In extending our analysis to IMHD and CMHD we find that the structure functions of both MHD models in the global and local perpendicular reference frames exhibit this behavior as well. For each MHD model in the local parallel reference frame the values of the longitudinal and transverse structure functions are very similar. In the global and local parallel frames the IMHD structure functions are greater than those of the CMHD model at lower values of the separation length where they appear more like the IHD structure functions in the global frame. In summary, the MHD models in the global frame and in the plane perpendicular to the local mean field are significantly more similar to the IHD structure functions than the structure functions of the MHD models in the local parallel frame. Also, the shapes of the structure functions in the global and local parallel frames for the IMHD model are more like that of the IHD model than the CMHD model is.

## Scaling Exponents

Self-similarity or scale-invariance implies that fluid turbulence at larger scales can be reproduced by the magnification of turbulent motions at smaller scales. At the dissipation scales self similarity fails as turbulence forms non-Gaussian dissipation structures (see bis03). Self-similarity is even an approximation in the inertial range where turbulent motions become increasingly sparse in space and time at smaller scales, a property known as intermittency. In 1962 Kolmogorov published a modification of his 1941 theory including the effects of intermittency or scale-variance (see kol62). One way to investigate intermittency is by studying the scaling exponents of velocity fluctuations in the inertial range, given that for a scale-invariant flow the scaling exponents are linear functions of the order of the structure functions.

Figure 8 - The sixth order longitudinal structure function as a function of the third order longitudinal structure function both calculated in the local parallel frame for the three models, along with their values of scaling exponents obtained from a minimum chi-square linear fit. Compare to figure 10 (a).

Despite the fact that the origin is not yet fully understood, the scaling behavior of structure functions is very useful and interesting (see laz06 for a review). Within the inertial range, a structure function of arbitrary order scales with the separation length $l$ as $S^{(p)}(l)\propto{l^{\zeta(p)}}$. From this scaling relation and a few transformations self-similarity between two structure functions of arbitrary orders, $p$ and $p'$, can be obtained: $$S^{(p)}(l) \propto{[S^{(p')}(l)]^{\zeta(p)/\zeta(p')}}, \label{eqn:sf_relation}$$ demonstrating how the structure functions $S^{(p)}(l)$ and $S^{(p')}(l)$ are related (see bis03). It also brings to light the possibility of using a structure function of an arbitrary order $p'$ as a normalization. In figure 8 the sixth order longitudinal structure function is plotted as a function of the third order longitudinal structure function where $\hat{l}$ is parallel to the local mean magnetic field, demonstrating how self similarity can be used to obtain scaling exponents. A minimum chi-square linear fit was used to determine $\zeta(6)$. The scaling relation $S^{(p)}(l) \propto{l^{\zeta(p)}}$ in the inertial range has been derived analytically in Kolmogorov phenomenology yielding $\zeta(p)=p/3$. For this reason the third moment is used as the normalization. If the effects of dissipation on the structure functions are independent of the order one may conclude that not only the third order moment but any moment $S^{(p')}(l) \propto{l^{\zeta(p')}}$ can be used (see bis03).
she94 proposed a scaling relation that contains three parameters (see pol95 and mul00): $g$ related to the scaling $\delta v\sim l^{\frac{1}{g}}$, $x$ related to the energy cascade rate $t^{-1}\sim l^{-x}$, and $C$ the codimension of the dissipative structures: $$\zeta(p) = \frac{p}{g}(1-x)+C[1-(1-\frac{x}{C})^{\frac{p}{g}}].$$ For IHD turbulence using $g=3$, $x=\frac{2}{3}$, and $C=2$ she94 obtained $$\zeta^{SL1}(p)=\frac{p}{9}+2[1-(\frac{2}{3})^{\frac{p}{3}}],$$ suggesting that dissipation occurs over one-dimensional structures, such as filaments or vortices (hereafter SL1). For three-dimensional MHD turbulence mul00 in extending the she94 model to IMHD turbulence and using $g=3$, $x=\frac{2}{3}$, and $C=1$ obtained $$\zeta^{SL2}(p)=\frac{p}{9}+1-(\frac{1}{3})^{\frac{p}{3}}$$ assuming the dissipative structures to be two-dimensional, such as sheets (hereafter SL2). Finally, for CMHD turbulence pad04 showed that the velocity field for subsonic MHD turbulence follows SL1, whiile for highly supersonic MHD turbulence it follows SL2.

 Figure 9 - Scaling exponents, normalized to the third order, as a function of the order of the structure functions in the global reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, (b) - transverse, and (c) - total). Figure 10 - Scaling exponents, normalized to the third order, as a function of the order of the structure functions in the local parallel reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, (b) - transverse, and (c) - total). For the IHD model the global scaling exponents are shown for comparison. Figure 11 - Scaling exponents, normalized to the third order, as a function of the order of the structure functions in the local perpendicular reference frame for the three models as well as K41, SL1, and SL2 ((a) - longitudinal, (b) - transverse, and (c) - total). For the IHD model the global scaling exponents are shown for comparison.

Figures 9 through 11 contain the scaling exponents as a function of the order of the longitudinal, transverse, and total structure functions normalized to the third order. For the MHD models the figures are organized as follows. In figure 9 the global reference frame is used, in figure 10 the local parallel reference frame is used, and in figure 11 the local perpendicular reference frame is used. For the IHD model the global scaling exponents are also presented in figures 10 and 11 for comparison, given that the local reference frame has no meaning without a magnetic field. Despite the fact that the error in the scaling exponents increases with order from the increasing fluctuation in the structure functions, emphasis is placed on high order structure functions ($p>4$) given that at low orders the theoretical predictions considered here differentiate themselves very little.
Beginning with figure 9 in the global reference frame the longitudinal scaling exponents of all three models follow SL1 quite well. In the transverse case although the slope of the IHD model initially follows SL1 it begins to flatten at around $p=7$. Its shape is in agreement with the results of got02. The IMHD model follows SL2 quite well while the CMHD model is between the two. The IHD results mentioned here agree with those obtained by got02 in that for $p>4$ the longitudinal scaling exponents fall slightly below SL1 while the transverse exponents are significantly lower. For all three models, despite the fact that the longitudinal scaling exponents follow SL1 much better than the transverse follow SL2, the scaling exponents calculated from the total structure functions are much more SL2-like than SL1-like for all three models. Their behavior is very similar to the transverse case, suggesting that the transverse scaling exponents have a significantly larger amplitude than the longitudinal, which can be seen to be true in figure 7. The same holds for figures 10 and 11.
In figure 10 in the local parallel reference frame for the longitudinal component, parallel to the local mean magnetic field in this case, the IMHD model is between SL1 and SL2 while the CMHD model mostly follows SL2. For the transverse component, perpendicular to the local mean magnetic field in this case, the IMHD model is between SL1 and SL2 while the CMHD model follows SL2.
In figure 11 in the local perpendicular reference frame for the longitudinal component, perpendicular to the local mean magnetic field in this case, the best agreement to SL1 is found for all three models at all values of $p$. For the transverse component, also pependicular to the local mean magnetic field in this case, we do not find a similar agreement; the IMHD model is more SL2-like while the CMHD model remains in between.
We confirm the results of got02 that for IHD turbulence the transverse scaling exponents, $\zeta_T(p)$ are less than the longitudinal $\zeta_L(p)$ for orders greater than four. In extending the analysis to incompressible and compressible MHD turbulence we find that in the global reference frame and in the local perpendicular reference frame this behavior continues. This is not necessarily the case in the local parallel reference frame. This difference is probably related to the fact that in the local parallel frame the longitudinal component necessarily probes the direction parallel to the local mean field while the transverse component probes the plane perpendicular to the local mean field, thus incorporating a consistent difference between the two components that does not exist in the other two frames.
The most conspicuous SL1 cases above are both longitudinal, namely $\zeta^{glo}_{L}(p)$ and $\zeta^{loc, \bot}_{L}(p)$ in which all models follow SL1 suggesting physically that the dissipation structures are one dimensional hydrodynamic-like vortices. The most conspicuous SL2 cases are both the transverse cases corresponding to the longitudinal SL1 cases above, namely $\zeta^{glo}_{T}(p)$ and $\zeta^{loc, \bot}_{T}(p)$ in which the IMHD model mostly follows SL2 suggesting physically that the dissipation structures are two dimensional sheets dominated by MHD processes. It is possible and perhaps more likely that the dissipation structures only appear one or two dimensional in slices of the three dimensional MHD turbulence (see cho03a). The cases in which the scaling exponents fall in between SL1 and SL2 may be some mixture of the two physical pictures. Also, in the global reference frame and in the local perpendicular reference frame the MHD models are similar. In these two cases it seems that the dissipation structures are similar for incompressible and compressible MHD turbulence.
HD and MHD turbulence are similar over high order longitudinal scaling exponents measured in the global reference frame as well as for motions perpendicular to the local mean field. In the local magnetic system the motions parallel to the local mean field for both MHD models are different from hydrodynamic motions.

## Fast, Slow, and Alfvenic Modes

Figure 12 - Separation method. We separate Alfven, slow, and fast modes in Fourier space by projecting the velocity Fourier component ${\bf v_k}$ onto bases ${\bf \xi}_A$, ${\bf \xi}_s$, and ${\bf \xi}_f$, respectively. Note that ${\bf \xi}_A = -\hat{\bf \varphi}$. Slow basis ${\bf \xi}_s$ and fast basis ${\bf \xi}_f$ lie in the plane defined by ${\bf B}_0$ and ${\bf k}$. Slow basis ${\bf \xi}_s$ lies between $-\hat{\bf \theta}$ and $\hat{\bf k}_{\|}$. Fast basis ${\bf \xi}_f$ lies between $\hat{\bf k}$ and $\hat{\bf k}_{\perp}$ (from Cho \& Lazarian 2003.

 Figure 13 - Decomposed scaling exponents, normalized to the third order, as a function of the order of the structure functions of the velocity field in the global reference frame for the IMHD and CMHD models as well as K41, SL1, and SL2 ((a) - Alfvenic, (b) - fast, and (c) - slow). Figure 14 - Decomposed scaling exponents, normalized to the third order, as a function of the order of the structure functions in the local parallel reference frame for the IMHD and CMHD models as well as K41, SL1, and SL2 ((a) - Alfvenic, (b) - fast, and (c) - slow). Figure 15 - Decomposed scaling exponents, normalized to the third order, as a function of the order of the structure functions in the local perpendicular reference frame for the IMHD and CMHD models as well as K41, SL1, and SL2 ((a) - Alfvenic, (b) - fast, and (c) - slow).

The compressible MHD equations support three types of linear waves or modes known as Alfvenic, fast, and slow. In the incompressible limit the fast mode does not exist. The velocity fields can be decomposed into these three components and scaling exponents for each component can be obtained, as explained in cho03b. Figure 12 gives a pictorial representation of the decomposition method (see cl03). Figures 13 through 15 contain decomposed scaling exponents and are organized as follows (similar to figures 9 through 11). In figure 13 the global reference frame is used, in figure 14 the local parallel reference frame is used, and in figure 15 the local perpendicular reference frame is used. In each figure the (a) is the Alfvenic mode, the (b) is the fast mode (containing only the CMHD model), and (c) is the slow mode. For the IMHD model the field of the slow mode was calculated by subtracting the Alfvenic field from the total velocity field. Each plot in figures 13 through 15 contains the scaling exponents derived from the longitudinal and transverse structure functions in order to expose the differences. For the Alfv\'en mode in figure 14, given that in the direction parallel to $B_{ext}$ there is no Alfv\'en mode, the longitudinal components have no physical meaning and are not included.
Beginning with figure 13 in the global reference frame, for the Alfv\'en mode the longitudinal components of both models follow SL1 while the transverse follow closer to SL2. For the the fast mode the CMHD longitudinal and transverse components are almost indistinguishable and follow SL2. For the slow mode the transverse components of both models are more SL2-like.
In figure 14 in the local parallel reference frame, for the Alfv\'en mode the CMHD transverse component follows SL2 quite well while the longitudinal component falls between SL1 and SL2. For the fast mode the CMHD transverse component falls well below SL2 while the longitudinal component falls between SL1 and SL2. For the slow mode the IMHD longitudinal component falls between SL1 and SL2 while the CMHD longitudinal component follows SL2 quite well. The transverse components of both models fall well below SL2 and for the IMHD model the curve actually becomes flat.
In figure 15 in the local perpendicular reference frame, for the Alfv\'en mode all cases follow SL1 except the IMHD transverse component which falls between SL1 and SL2. For the CMHD model the longitudinal and transverse components are indistinguishable. For the fast mode both components of the CMHD model follow closer to SL2, the transverse component more so. For the slow mode the IMHD longitudinal component follows closer to SL1 while the transverse component follows SL2 fairly well. For the CMHD model this is the only case where the transverse component is significantly greater than the longitudinal component. The transverse component follows SL2 while the longitudinal component falls well below SL2.
There are three cases in which scaling exponents fall well below SL2, namely the CMHD longitudinal component of $\zeta^{loc, \|}_{fast}(p)$, the transverse component for both models of $\zeta^{loc, \|}_{slow}(p)$, and the CMHD longitudinal component of $\zeta^{loc, \bot}_{slow}(p)$, suggesting that intermittency is greater in these cases.
Also, $\zeta_L(p)\geq\zeta_T(p)$ for all cases except for the CMHD model of $\zeta^{loc, \bot}_{slow}(p)$. This is not reflected in figure 11 where for the CMHD model $\zeta_L(p)>\zeta_T(p)$. This suggests that the amplitude of the slow mode for this model is small compared to the amplitude of the Alfvenic and fast modes. The two cases where $\zeta_L(p)=\zeta_T(p)$ are for the CMHD models of $\zeta^{glo}_{fast}(p)$ and $\zeta^{loc, \bot}_{alfv}(p)$.
In agreement with the above-mentioned property, $\zeta_L(p)\geq\zeta_T(p)$, it seems that the longitudinal components are more likely to be SL1-like while the transverse components are more likely to be SL2-like. This is the case for $\zeta^{glo}_{alfv}(p)$ and $\zeta^{loc, \bot}_{slow}(p)$, which is in direct agreement with the conspicuous cases mentioned in the discussion of figures 9 through 11. Combining the results of figures 9 and 13 it seems that in the global reference frame the Alfvenic mode is mostly responsible for the fact that the longitudinal components of both MHD models follow SL1. For the IMHD model both the Alfvenic and slow modes seem to contribute significantly to the fact that the transverse components appear to follow SL2. Combining the results of figures 11 and 15 it seems that in the local perpendicular reference frame once again the Alfvenic mode is mostly responsible for the fact that the longitudinal components of both MHD models follow SL1. For the IMHD model in this case the slow mode seems to contribute most to the fact that the transverse components appear to follow SL2.
It is also true that for the decomposed scaling exponents of the two MHD models the longitudinal components are more alike than the transverse components, in every case except for $\zeta^{loc, \bot}_{slow}(p)$. This makes sense in the local parallel reference frame where there is a bias in that the longitudinal components are chosen in the direction of the local mean field while the transverse components are chosen randomly to be perpendicular to the longitudinal components, but this observation also holds in the global reference frame where there is no bias because the longitudinal components are chosen randomly. If the dissipation structures are similar for compressible and incompressible MHD then it is possible that the longitudinal scaling exponents would do a better job of probing this.

## Discussion and Summary

In this paper turbulence in an incompressible fluid with and without a magnetic field as well as moderately compressible MHD turbulence are considered. We summarize our findings as follows:
• Substantial similarity was shown between the spectra for all three models.
• The PDF's of the velocity components perpendicular to the external magnetic field are IHD-like while those parallel to the field have a significantly smaller range of velocities. The PDF's of the transverse velocity increments for the MHD models decay slower than the IHD model.
• The structure functions of the MHD models in the global frame and in the plane perpendicular to the local mean field are IHD-like in that the values of the transverse structure functions are greater than those of the longitudinal. In these two frames the shapes of the structure functions of the IMHD model are more IHD-like than the CMHD model.
• The similarity of HD, IMHD, and CMHD turbulence persists over high order longitudinal structure function scaling exponents measured in the global reference frame as well as for motions perpendicular to the local mean field. In these two frames the longitudinal scaling exponents of both MHD models seem to follow the predictions for one dimensional IHD-like dissipations structures while the transverse scaling exponents of the IMHD model seems to follow the predictions for two dimensional IMHD-like dissipation structures. In the local magnetic system the motions parallel to the local mean field for both MHD models are different from hydrodynamic motions.
• In the global reference frame and for motions perpendicular to the local mean field the Alfvenic mode is mostly responsible for the fact that the longitudinal components of both MHD models follow the IHD model and the IHD theoretical dissipation structure predictions. In the global reference frame for the IMHD model both the Alfvenic and slow modes seem to contribute significantly to the fact that the transverse components follow the IMHD theoretical predictions. For motions perpendicular to the magnetic field for the IMHD model the slow mode seems to contribute most to the fact that the transverse component seems to follow the IMHD theoretical predictions. It is also true that for the decomposed scaling exponents of the two MHD models the longitudinal components are more alike than the transverse components, in every case except for the slow mode in the plane perpendicular to the local mean field.

## Bibliography

kol41 - Kolmogorov, A. N., "The local structure of turbulence in incompressible viscous fluid for very large reynolds numbers, " Dok. Akad. Nauk SSSR p30, 9 (1941).
bis03 - Biskamp, D., "Magnetohydrodynamic Turbulence," (Cambridge University Press, Cambridge, 2003).
got02 - Gotoh, T., Fukayama,D., & Nakano,T., "Velocity field statistics in homogeneous steady turbulence obtained using a high-resolution direct numerical simulation," Phys. Fluids p14, 1065 (2002).
bor96 - Borue, V. & Orszag, S.A., "Numerical study of three-dimensional Kolmogorov flow at high Reynolds numbers," J. Fluid Mech., 306, 293 (1996).
cl02 - Cho, J. & Lazarian, A., "Compressible sub-Alfvenic MHD turbulence in low-$\beta$ plasmas," Phys. Rev. Lett., p88, 245001 (2002).
cho02 - Cho, J., Lazarian, A., & Vishniac, E. T., "Simulations of MHD turbulence in a strongly magnetized medium," Astrophys. J., p564, 291 (2002).
cl03 - Cho, J. & Lazarian, A., "Compressible magnetohydrodynamic turbulence: mode coupling, scaling relations, anisotropy, new regime and astrophysical implications," Mon. Not. R. Astron. Soc., p345, 325 (2003).
fal94 - Falkovich, G., "Bottleneck phenomenon in developed turbulence," Phys. Fluids p6, 1411 (1994).
gol95 - Goldreich, P., Sridhar, S., "Toward a theory of interstellar turbulence. 2: Strong alfvenic turbulence," Astrophys. J., p438, 763 (1995).
cl05 - Cho, J., & Lazarian, A., "Generation of compressible modes in MHD turbulence," Theor. and Comp. Fl. Dynam., p19, 127 (2005).
bol06 - Boldyrev, Stanislav, "Spectrum of magnetohydrodynamic turbulence," Phys. Rev. Lett., p96, 115002 (2006) ber06 - Beresnyak, A., Lazarian, A., "Polarization intermittency and its influence on MHD turbulence," Astrophys. J., p640, 175 (2006).
cho00a - Cho, J. & Vishniac, E., "The anisotropy of MHD Alfvenic turbulence," Astrophys. J., p539, 273 (2000a).
kol62 - Kolmogorov, A. N., "A refinement of previous hypothesis concerning the local structure of turbulence in a viscous incompressible fluid at high Reynolds number," J. Fluid Mech. p13, 82 (1962).
laz06 - Lazarian, A., "Intermittency of magnetohydrodynamic turbulence: an astrophysical perspective," Int. J. Mod. Phys. D, p15(6), 1 (2006).
pol95 - Politano, H. & Pouquet, A., "Model of intermittency in magnetohydrodynamic turbulence," Phys. Rev. E, p52, 636 (1995).
she94 - She, Z. & Leveque, E., "Universal scaling laws in fully developed turbulence," Phys. Rev. Lett. p72, 336 (1994).
mul00 - M\"uller, W.-C. & Biskamp, D., "Scaling properties of three-dimensional magnetohydrodynamic turbulence," Phys. Rev. Lett., p84(3), 475 (2000).
pad04 - Padoan, P., Jimenez, R., Nordlund, A., & Boldyrev, S., "Structure function scaling in compressible super-Alfvenic MHD furbulence," Phys. Rev. Lett., p92, 191102 (2004).
cho03a - Cho, J., Lazarian, A., & Vishniac, E. T., "Ordinary and viscosity-damped MHD turbulence," Astrophys. J., p595, 812-23 (2003a).
cho03b - Cho, J., Lazarian, A., & Vishniac, E. T., "MHD turbulence: scaling laws and astrophysical implications," Lect. Notes Phys. p614, 56 (2003b).
cho00b - Cho, J. & Vishniac, E., "The generation of magnetic fields through driven turbulence," Astrophys. J., p538, 217 (2000b).