 for $%
k^{\prime }>2$. We take $B_0^{+-}=0$. The Eqs. (\[rg1\],\[rg2\]) are solved iteratively. We continue iterating the equations till $%
B_{n+1}^{*}(k^{\prime })\approx B_n^{*}(k^{\prime })$, that is, till the solution converges. For $B_0^{initial}=1.0,$ the $B_n^{\prime }$s for various $n$ ranging from $0 \ldots 3$ is shown in Figure 1. Here the convergence is very fast, and after $n=3-4$ iterations $B_n^{*}(k)$ converges to an universal function $$f(k^{\prime })=1.24*k^{\prime -0.32}.$$ From the above arguments, we have shown that $B_n^{*}(k^{\prime })\ $is approximately proportional to $k^{\prime -1/3}$. The other parameter $%
B_n^{*+-}(k^{\prime })$ remains close to zero.

We infer from the above analysis that the mean magnetic field scales as $%
k^{-1/3}$, and the energy spectra scales as $k^{-5/3}$. Essentially, the scaling of $B_0$ leads to $k^{-5/3}$ energy (Kolmogorov-like) spectra in our scheme. We have calculated $B_n^{*}(k^{\prime })$ for $B_0^{initial}=1,2,10$ and found that for large $n$, $B_n^{*}(k^{\prime })\approx
1.25B_0^{initial}k^{\prime -1/3}$ or $$B_n(k)=1.25B_0^{initial}K^{1/2}\Pi ^{1/3}k^{-1/3}.$$

Calculation of $K$
------------------

We can calculate the Kolmogorov’s constant for MHD turbulence $K$ by calculating the cascade rate $\Pi $ [@Lesl]. In MHD the cascade rates are $$\Pi ^{+}(k)=\Pi ^{-}(k)=-\int_0^kdk^{\prime }T(k^{\prime })$$ The numerical solution of the cascade rate integral yields [@Lesl] $$\label{alpha}\frac{1.24B_0^{initial}}{K^{3/2}}=3.85$$ From the above
for $% 
 k^{\prime } > 2$. We take $ B_0^{+-}=0$. The Eqs. (\[rg1\],\[rg2\ ]) are solved iteratively. We continue iterating the equation till $% 
 B_{n+1}^{*}(k^{\prime }) \approx B_n^{*}(k^{\prime }) $, that is, till the solution converges. For $ B_0^{initial}=1.0,$ the $ B_n^{\prime } $ south for various $ n$ ranging from $ 0 \ldots 3 $ is shown in Figure 1. Here the overlap is very fast, and after $ n=3 - 4 $ iterations $ B_n^{*}(k)$ converges to an universal affair $ $ f(k^{\prime }) = 1.24*k^{\prime -0.32}.$$ From the above arguments, we have shown that $ B_n^{*}(k^{\prime }) \ $ is approximately proportional to $ k^{\prime -1/3}$. The early argument $% 
 B_n^{*+-}(k^{\prime }) $ remains close to zero. 

 We guess from the above analysis that the mean magnetic playing field scales as $% 
 k^{-1/3}$, and the department of energy spectra scale as $ k^{-5/3}$. basically, the scaling of $ B_0 $ leads to $ k^{-5/3}$ energy (Kolmogorov - like) spectrum in our scheme. We have calculated $ B_n^{*}(k^{\prime }) $ for $ B_0^{initial}=1,2,10 $ and find that for big $ n$, $ B_n^{*}(k^{\prime }) \approx 
 1.25B_0^{initial}k^{\prime -1/3}$ or $ $ B_n(k)=1.25B_0^{initial}K^{1/2}\Pi ^{1/3}k^{-1/3}.$$ 

 Calculation of $ K$ 
 ------------------ 

 We can calculate the Kolmogorov ’s constant for MHD turbulence $ K$ by forecast the cascade rate $ \Pi $ [ @Lesl ]. In MHD the cascade pace are $ $ \Pi ^{+}(k)=\Pi ^{-}(k)=-\int_0^kdk^{\prime } T(k^{\prime }) $ $ The numerical solution of the cascade rate integral yields [ @Lesl ] $ $ \label{alpha}\frac{1.24B_0^{initial}}{K^{3/2}}=3.85$$ From the above
 fog $%
k^{\prime }>2$. We take $B_0^{+-}=0$. The Tqs. (\[rg1\],\[rg2\]) are solvgd iteravively. Se contivue iterating the equations vill $%
B_{n+1}^{*}(k^{\peime })\approx B_n^{*}(k^{\prime })$, ghat is, tpll the silutmon converges. Foc $B_0^{initial}=1.0,$ the $N_n^{\prikx }$s for various $n$ ranging from $0 \ldots 3$ hs skown in Figure 1. Here the convergence is veru vast, and after $n=3-4$ iuerwtiohs $B_n^{*}(k)$ converges to an universal fhnction $$f(k^{\prime })=1.24*k^{\prike -0.32}.$$ From the above argumenhs, wf have shown that $H_n^{*}(k^{\prime })\ $iw ap[eoximately pfoportional to $k^{\prime -1/3}$. The other parameter $%
B_n^{*+-}(k^{\prime })$ femaius close to zwro.

Ag infer from the wbove analyslx that the mesn magnetic ficld sralew as $%
k^{-1/3}$, and the energy spectra scales as $k^{-5/3}$. Essenthamly, the scaling od $V_0$ leags tm $k^{-5/3}$ dberey (Iokmkgorov-pikx) spectra ih our schemw. We have calculatec $F_b^{*}(k^{\prime })$ for $G_0^{initiwl}=1,2,10$ and found that for large $n$, $B_n^{*}(k^{\prime })\apkrox
1.25B_0^{jnitial}k^{\prime -1/3}$ or $$B_n(k)=1.25B_0^{ibitial}K^{1/2}\Pi ^{1/3}k^{-1/3}.$$

Calculatioj of $K$
------------------

We san calculate the Kolmogorov’s constant for MHD tusbuleicd $K$ bn cauxupating the cascade rate $\Pi $ [@Lesl]. In MHD the cwacsdv rates are $$\Pi ^{+}(k)=\Pl ^{-}(k)=-\int_0^kdk^{\prime }T(k^{\ptile })$$ The numericau solucjoh of the cascade rwte intggral tields [@Lefl] $$\lsbel{alpha}\frac{1.24B_0^{initial}}{K^{3/2}}=3.85$$ From the above
for $% k^{\prime }>2$. We take $B_0^{+-}=0$. (\[rg1\],\[rg2\]) solved iteratively. continue iterating the B_n^{*}(k^{\prime that is, till solution converges. For the $B_n^{\prime }$s for various $n$ from $0 \ldots 3$ is shown in Figure 1. Here the convergence is fast, and after $n=3-4$ iterations $B_n^{*}(k)$ converges to an universal function $$f(k^{\prime })=1.24*k^{\prime From above we shown that $B_n^{*}(k^{\prime })\ $is approximately proportional to $k^{\prime -1/3}$. The other parameter $% B_n^{*+-}(k^{\prime })$ close to zero. We infer from the above that the mean magnetic scales as $% k^{-1/3}$, and energy scales as Essentially, scaling $B_0$ leads to energy (Kolmogorov-like) spectra in our scheme. We have calculated $B_n^{*}(k^{\prime })$ for $B_0^{initial}=1,2,10$ and found that for $n$, $B_n^{*}(k^{\prime -1/3}$ or ^{1/3}k^{-1/3}.$$ of ------------------ We can Kolmogorov’s constant for MHD turbulence $K$ cascade rate $\Pi $ [@Lesl]. In MHD the rates are ^{+}(k)=\Pi ^{-}(k)=-\int_0^kdk^{\prime }T(k^{\prime })$$ The numerical of the cascade rate integral yields [@Lesl] $$\label{alpha}\frac{1.24B_0^{initial}}{K^{3/2}}=3.85$$ the above
 for $%
k^{\prime }>2$. We take $B_0^{+-}=0$. The Eqs. (\[rg1\],\[rG2\]) are solved IteraTivEly. we ContInue Iterating the eqUAtioNs till $%
B_{n+1}^{*}(k^{\prime })\approx B_n^{*}(K^{\primE })$, tHAt is, TIlL the sOlution COnVERgeS. FOr $b_0^{inItIAl}=1.0,$ The $B_n^{\PriMe }$s for vArious $n$ ranGinG fRom $0 \ldots 3$ is shOWn In Figure 1. HeRe tHe convergencE is Very faSt, And AFter $n=3-4$ IteRatioNs $B_n^{*}(k)$ cONvergeS to an univErSAl funcTIon $$f(k^{\prIME })=1.24*k^{\PrimE -0.32}.$$ From the above arguMEnTS, we have shown thAt $B_n^{*}(k^{\pRiME })\ $iS APprOxiMately propOrTionaL To $k^{\primE -1/3}$. thE OTHer PArameter $%
B_n^{*+-}(k^{\prIme })$ remains cLOse To zero.

we InfER from tHe aboVe ANalYsis that the Mean Magnetic fIeld scALes as $%
k^{-1/3}$, aND the eneRgy speCtrA scAles AS $k^{-5/3}$. esSenTiALly, THe ScaLIng Of $B_0$ leads To $K^{-5/3}$ eNergy (kolmOGOROv-liKe) sPectRa in oUr scheme. We havE caLculATed $b_n^{*}(k^{\prIme })$ foR $B_0^{inItIal}=1,2,10$ anD found That fOr Large $n$, $B_n^{*}(k^{\prime })\aPproX
1.25B_0^{initial}K^{\prImE -1/3}$ or $$b_n(K)=1.25B_0^{iniTIal}K^{1/2}\Pi ^{1/3}K^{-1/3}.$$

CaLcuLation oF $K$
------------------

We can CAlcUlATE ThE Kolmogorov’s constaNt FOR MhD turbulEnce $K$ bY CaLcULating thE cAscAde rATE $\Pi $ [@LeSl]. In mhD The cascaDe rateS ArE $$\PI ^{+}(k)=\Pi ^{-}(k)=-\inT_0^kDk^{\primE }T(K^{\prIme })$$ the nuMEricAl soluTion of thE cascADe rate integral YIelds [@Lesl] $$\labeL{AlPHA}\fRAc{1.24B_0^{iNitIal}}{K^{3/2}}=3.85$$ From the AbovE
 for $%
k^{\prime }>2$. We  take $B_0 ^{+-} =0$ . T he  Eqs . (\ [rg1\],\[rg2\] )  are  solved iteratively. W e con ti n ue i t er ating  the eq u at i o nsti ll  $% 
B _ {n +1}^{ *}( k^{\pri me })\appr oxB_ n^{*}(k^{\pr i me  })$, that  is , till the s olu tion c on ver g es. F or$B_0^ {initi a l}=1.0 ,$ the $B _n ^ {\prim e  }$s fo r va riou s $n$ ranging fro m  $ 0  \ldots 3$ isshownin Fi g u re1.Here the c on verge n ce is v e ry f a st, and after $n= 3-4$ iterat i ons  $B_n^ {* }(k ) $ conv erges  t o  an  universalfunc tion $$f( k^{\pr i me })=1 . 24*k^{\ prime-0. 32} .$$F ro mthe  a b ove ar gum e nts , we hav esh own t hat$ B _ n ^{*} (k^ {\pr ime } )\ $is approx ima tely pro porti onalto $ k^ {\pri me -1/ 3}$.Th e other paramet er $ %
B_n^{*+ -}( k^ {\p ri me }) $  remai nsclo se to z ero.

W e  in fe r f ro m the above analys is t ha t the me an mag n et ic field sc al esas $ % k^{-1 /3}$ ,  a nd the e nergys pe ct ra scal es  as $k ^{ -5/ 3}$ . Ess e ntia lly, t he scali ng of $B_0$ leads to $k^{-5/3}$ en e rg y (K o lmog oro v-like) spe ctra in o ur s c he me. We ha ve ca lc u la t ed $B_n^{*}(k^{\pri me  })$ f or $B _0^{initial}= 1,2,10$ an d f ound tha t fo r  l a rge $n$, $B_n^ {*}(k ^{\prime } ) \approx1.25B _0^{init ial}k^{\p r i me -1/3} $ o r $ $B_ n(k ) = 1. 25B_0^{initia l } K^{1 /2 }\Pi ^{ 1/3 }k^{-1/ 3}. $$
Ca lcu la tion of $ K$
----- -- -- -- -- --- --

W e  can cal cu lat ethe  Kolm o gorov’ s con stan tfo r  MH D turbu l en c e  $K$  b ycalc ula ti ng th e ca s cad e rate$\Pi $ [@ Les l ]. I nMH D the c ascade ratesar e $$\Pi ^{ +} (k) =\Pi ^ { - }(k)=-\i nt_0^kdk^{\prime }T(k^{ \ prime } )$$  Thenume rical sol uti on ofthe cascad e rate  inte gr aly i elds[ @ Le sl]  $ $\label{al p h a}\ frac{ 1. 24B_ 0^{init ial}}{K^{3/2}}=3.8 5 $$From the abov e
 for_$%
k^{\prime }>2$._We take $B_0^{+-}=0$. The_Eqs. (\[rg1\],\[rg2\])_are_solved iteratively._We_continue iterating the_equations till $%
B_{n+1}^{*}(k^{\prime_})\approx B_n^{*}(k^{\prime })$, that_is, till the_solution_converges. For $B_0^{initial}=1.0,$ the $B_n^{\prime }$s for various $n$ ranging from $0 \ldots 3$_is_shown in_Figure_1._Here the convergence is very_fast, and after $n=3-4$ iterations_$B_n^{*}(k)$ converges_to an universal function $$f(k^{\prime })=1.24*k^{\prime -0.32}.$$ From_the_above arguments, we_have shown that $B_n^{*}(k^{\prime })\ $is approximately proportional to_$k^{\prime -1/3}$. The other parameter $%
B_n^{*+-}(k^{\prime_})$ remains close_to_zero.

We_infer from the above_analysis that the mean magnetic field_scales as $%
k^{-1/3}$, and the energy_spectra scales as $k^{-5/3}$. Essentially, the scaling_of $B_0$ leads to $k^{-5/3}$ energy_(Kolmogorov-like) spectra in our scheme._We have_calculated $B_n^{*}(k^{\prime })$ for $B_0^{initial}=1,2,10$_and found that_for large_$n$, $B_n^{*}(k^{\prime })\approx
1.25B_0^{initial}k^{\prime_-1/3}$ or $$B_n(k)=1.25B_0^{initial}K^{1/2}\Pi ^{1/3}k^{-1/3}.$$

Calculation of $K$
------------------

We_can calculate the_Kolmogorov’s constant for MHD turbulence $K$_by_calculating the cascade_rate_$\Pi_$ [@Lesl]._In MHD the_cascade_rates are_$$\Pi_^{+}(k)=\Pi ^{-}(k)=-\int_0^kdk^{\prime }T(k^{\prime })$$ The numerical_solution_of the cascade rate integral yields [@Lesl]_$$\label{alpha}\frac{1.24B_0^{initial}}{K^{3/2}}=3.85$$ From the above