*$ stands for the convolution product), where $\theta$ is a radially symmetrical function, such that $\int_{\mathbb R ^3} \theta (\textbf{x})d\textbf{x}=1$ and $ \theta_\epsilon (\textbf{x})= \frac{1}{\epsilon^3}\theta (\textbf{x}/\epsilon)$[@RobVar08].

Numerical study
===============

We would like now to study the statistical properties of the velocity field defined by Eq. (\[eq:K41Field\]). Analytical formulas are difficult to obtain, thus we will focus on numerical simulations. To do so, one has to choose the short time scale $t=\tau$. It is easy to check that the variance of the matrix $\textbf{S}_0$ (defined as $\langle \mbox{tr } \textbf{S}_0^2\rangle$) goes to infinity as the small scale parameter $\epsilon$ goes to zero. So we take for $\tau$ the local normalizing value $\tau =  (\mbox{tr } \textbf{S}_0^2)^{-1/2}$, in the spirit of the RFD closures provided in Ref. [@CheA2].

The simulation is performed in a 1-periodic box with $N^3$ collocation points. The infinitesimal volume is given by $dV=dx^3$, with $dx=1/N$. We choose as a regularizing function and large-scale cut-off the isotropic normalized Gaussian function $\varphi_{L} (\textbf{x}) = \theta_{L}(\textbf{x}) = \left(\frac{6}{\pi L^2}\right)^{3/2}\exp (-6|\textbf{x}|^2/L^2)$. This allows to compute analytically the regularized norm of a vector $\textbf{x}$, useful for numerical purposes: $|\textbf{x}|_\epsilon = \frac{\epsilon}{\sqrt{6\pi}}e^{-6\frac{|\textbf{x}|^2}{\epsilon^2}} + \left(|\textbf{x}| + \frac{\epsilon^2}{12|\textbf{x}|}\right)\mbox{erf}\left( \frac{\sqrt{6}|\textbf{x}|}{\epsilon}\right)\mbox{,}$ with erf the error function. The small-scale cut-off is chosen as $\epsilon=2dx$ and the large one as $L=1/2$. The kernels
* $ stands for the convolution product), where $ \theta$ is a radially symmetrical affair, such that $ \int_{\mathbb radius ^3 } \theta (\textbf{x})d\textbf{x}=1 $ and $ \theta_\epsilon (\textbf{x})= \frac{1}{\epsilon^3}\theta (\textbf{x}/\epsilon)$[@RobVar08 ]. 

 Numerical study 
 = = = = = = = = = = = = = = = 

 We would like now to learn the statistical properties of the velocity playing field define by Eq. (\[eq: K41Field\ ]). Analytical formula are difficult to obtain, therefore we will focus on numerical simulations. To do therefore, one has to choose the short time scale $ t=\tau$. It is easy to see that the variance of the matrix $ \textbf{S}_0 $ (defined as $ \langle \mbox{tr } \textbf{S}_0 ^ 2\rangle$) goes to eternity as the small scale argument $ \epsilon$ goes to zero. then we take for $ \tau$ the local normalizing value $ \tau =   (\mbox{tr } \textbf{S}_0 ^ 2)^{-1/2}$, in the heart of the RFD closures provided in Ref. [ @CheA2 ]. 

 The simulation is perform in a 1 - periodic box with $ N^3 $ collocation points. The infinitesimal volume is given by $ dV = dx^3 $, with $ dx=1 / N$. We choose as a regularizing function and large - scale cut - off the isotropic normalized Gaussian function $ \varphi_{L } (\textbf{x }) = \theta_{L}(\textbf{x }) = \left(\frac{6}{\pi L^2}\right)^{3/2}\exp (-6|\textbf{x}|^2 / L^2)$. This allows to compute analytically the regularized norm of a vector $ \textbf{x}$, useful for numerical purpose: $ |\textbf{x}|_\epsilon = \frac{\epsilon}{\sqrt{6\pi}}e^{-6\frac{|\textbf{x}|^2}{\epsilon^2 } } + \left(|\textbf{x}| + \frac{\epsilon^2}{12|\textbf{x}|}\right)\mbox{erf}\left (\frac{\sqrt{6}|\textbf{x}|}{\epsilon}\right)\mbox{,}$ with erf the error function. The belittled - scale cut - off is chosen as $ \epsilon=2dx$ and the large one as $ L=1/2$. The kernel
*$ stwnds for the convolution product), where $\jhwta$ is a radjally syometrical function, such that $\ibt_{\matybb R ^3} \theta (\textbf{x})d\tdxtbf{x}=1$ anf $ \theta_\wpsioin (\textbf{x})= \frac{1}{\epsljon^3}\tgcta (\tzxvbf{x}/\epsilon)$[@RobVsr08].

Numerican study
===============

We woulg uine now to study the statistical pro[erties ov the velocity fiekq derpntd by Eq. (\[eq:K41Field\]). Analytical formhlas art difficult to obtsin, thus we will focus on jumegical simulations. Ho do so, ong haf to choose tfe short tpke scale $t=\jau$. It is easy to check that the varicnce of the mqtrlf $\textbf{S}_0$ (dxfined as $\langle \mnpx{tr } \dextbf{S}_0^2\tangle$) goes to inhiniry as the small scale parameter $\epsilon$ goes to sexo. So we take for $\tau$ tye lowal tormxoizkng velus $\tau =  (\muox{tr } \textgf{S}_0^2)^{-1/2}$, in the wpirit of the RFD ckofltes provided in Res. [@SheA2].

The simulation is performed in a 1-pegiodjc box with $N^3$ collocatiin points. The infinitgsimal voltme is given by $dV=dx^3$, with $dx=1/N$. We choose as a regunarizmne fbkgtiov ajd large-scale cut-off the isotropic normalized Faisxian function $\yarphi_{L} (\textbf{x}) = \tneha_{K}(\jextbf{x}) = \left(\ftac{6}{\pi L^2}\xjggt)^{3/2}\exp (-6|\textbf{x}|^2/L^2)$. Thid allowf to xompute agalyyically the regularized norn of a vectog $\teztbf{x}$, useful for nbmerical pur'oses: $|\jextbf{c}|_\epsilon = \frac{\epsilon}{\sqxt{6\pi}}e^{-6\fdac{|\textbf{x}|^2}{\eosilon^2}} + \lswt(|\textbf{x}| + \frac{\eosikot^2}{12|\textbf{x}|}\right)\mbox{erf}\left( \fwac{\sqrt{6}|\teetbf{x}|}{\zpsilon}\rkght)\kbox{,}$ wyth erf thf error function. The smapl-scape cut-off is chosen as $\epsilon=2dx$ and the lacje one as $L=1/2$. Tne kegnels
*$ stands for the convolution product), where a symmetrical function, that $\int_{\mathbb R \theta_\epsilon \frac{1}{\epsilon^3}\theta (\textbf{x}/\epsilon)$[@RobVar08]. Numerical =============== We would now to study the statistical properties the velocity field defined by Eq. (\[eq:K41Field\]). Analytical formulas are difficult to obtain, we will focus on numerical simulations. To do so, one has to choose short scale It easy to check that the variance of the matrix $\textbf{S}_0$ (defined as $\langle \mbox{tr } \textbf{S}_0^2\rangle$) to infinity as the small scale parameter $\epsilon$ to zero. So we for $\tau$ the local normalizing $\tau (\mbox{tr } in spirit the RFD closures in Ref. [@CheA2]. The simulation is performed in a 1-periodic box with $N^3$ collocation points. The infinitesimal is given with $dx=1/N$. choose a function and large-scale isotropic normalized Gaussian function $\varphi_{L} (\textbf{x}) \left(\frac{6}{\pi L^2}\right)^{3/2}\exp (-6|\textbf{x}|^2/L^2)$. This allows to compute analytically regularized norm a vector $\textbf{x}$, useful for numerical $|\textbf{x}|_\epsilon = \frac{\epsilon}{\sqrt{6\pi}}e^{-6\frac{|\textbf{x}|^2}{\epsilon^2}} + \left(|\textbf{x}| + \frac{\epsilon^2}{12|\textbf{x}|}\right)\mbox{erf}\left( \frac{\sqrt{6}|\textbf{x}|}{\epsilon}\right)\mbox{,}$ erf the error function. The small-scale cut-off is chosen as $\epsilon=2dx$ and the large one The kernels
*$ stands for the convolution prOduct), where $\Theta$ Is a RadIaLly sYmmeTrical function, SUch tHat $\int_{\mathbb R ^3} \theta (\textBf{x})d\tExTBf{x}=1$ aND $ \tHeta_\ePsilon (\tEXtBF{X})= \frAc{1}{\EpSilOn^3}\THeTa (\texTbf{X}/\epsiloN)$[@RobVar08].

NumEriCaL study
===============

We woulD LiKe now to stuDy tHe statisticaL prOpertiEs Of tHE veloCitY fielD definED by Eq. (\[eQ:K41Field\]). AnAlYTical fORmulas aRE DiFficUlt to obtain, thus we WIlL Focus on numericAl simuLaTIoNS. to dO so, One has to chOoSe the SHort timE ScALE $T=\taU$. it is easy to cheCk that the vaRIanCe of thE mAtrIX $\textbF{S}_0$ (defInED as $\Langle \mbox{tR } \texTbf{S}_0^2\ranglE$) goes tO InfinitY As the smAll scaLe pAraMeteR $\EpSiLon$ GoES to ZErO. So WE taKe for $\tau$ ThE lOcal nOrmaLIZINg vaLue $\Tau =  (\mBox{tr } \Textbf{S}_0^2)^{-1/2}$, in the sPirIt of THe RfD cloSures ProvIdEd in REf. [@CheA2].

the siMuLation is performEd in A 1-periodic Box WiTh $N^3$ CoLlocaTIon poiNts. the InfinitEsimal vOLumE iS GIVeN by $dV=dx^3$, with $dx=1/N$. We chOoSE As A regularIzing fUNcTiON and largE-sCalE cut-OFF the iSotrOPiC normaliZed GauSSiAn FunctioN $\vArphi_{L} (\TeXtbF{x}) = \tHeta_{L}(\TExtbF{x}) = \left(\Frac{6}{\pi L^2}\rIght)^{3/2}\eXP (-6|\textbf{x}|^2/L^2)$. This aLLows to compute ANaLYTiCAlly The Regularized Norm OF a veCtor $\TExTbf{X}$, UsefuL for nUmERiCAl purposes: $|\textbf{x}|_\epSiLon = \fraC{\epsiLon}{\sqrt{6\pi}}e^{-6\fraC{|\textbf{x}|^2}{\epSILOn^2}} + \left(|\teXtbf{X}| + \FrAC{\epsilon^2}{12|\textbf{X}|}\righT)\mbox{erf}\leFT( \frac{\sqrT{6}|\textBf{x}|}{\epsilOn}\right)\mbOX{,}$ With erf tHe eRroR fuNctION. THe small-scale cUT-Off iS cHosen as $\EpsIlon=2dx$ aNd tHe lArgE onE aS $L=1/2$. The kernEls
*$ stands for the convolut ion produc t), w her e $ \t heta $ is  a radially sy m metr ical function, such th at $\ in t _{\m a th bb R^3} \th e ta ( \te xt bf {x} )d \ te xtbf{ x}= 1$ and$ \theta_\ eps il on (\textbf{ x }) = \frac{1} {\e psilon^3}\th eta  (\tex tb f{x } /\eps ilo n)$[@ RobVar 0 8].

N umericalst u dy
=== = ======= = = == 

We  would like now t o  s t udy the statis ticalpr o pe r t ies  of  the veloc it y fie l d defin e db y Eq. (\[eq:K41Fiel d\]). Analy t ica l form ul asa re dif ficul tt o o btain, thus  wewill focu s on n u merical simulat ions.Todoso,o ne  h asto cho o se  th e  sh ort time  s ca le $t =\ta u $ . It i s e asyto ch eck that thevar ianc e  of  thematri x $\ te xtbf{ S}_0$(defi ne d as $\langle \ mbox {tr } \te xtb f{ S}_ 0^ 2\ran g le$) g oes  to  infini ty as t h e s ma l l sc ale parameter $\ep si l o n$  goes to  zero. So  w e  take fo r$\t au$t h e loc al n o rm alizingvalue$ \t au  =  (\m bo x{tr }  \ tex tbf {S}_0 ^ 2)^{ -1/2}$ , in the  spir i t of the RFD c l osures provid e di n  R e f. [ @Ch eA2].

Thesimu l atio n is pe rfo r med i n a 1 -p e ri o dic box with $N^3$co llocat ion p oints. The in finitesima l v olume is  giv e nb y $dV=dx^3$, w ith $ dx=1/N$. W e  chooseas aregulari zing func t i on and l arg e-s cal e c u t -o ff the isotro p i c no rm alizedGau ssian f unc tio n $ \va rp hi_{L} (\ textbf{x })  =  \ th eta _{L}( \ textbf{x })  =\l eft (\fra c {6}{\p i L^2 }\ri gh t) ^ {3/ 2}\exp( -6 | \ text bf {x }|^2 /L^ 2) $. Th is a l low s to co mpute ana lyt i call yth e regul arized norm o fa vector $ \t ext bf{x}$ , useful f or numerical purposes:$ |\textb f{x }|_\e psil on = \fra c{\ epsilo n}{ \ sqrt{6 \pi}}e ^{-6\ fr ac{ | \ textb f { x} |^2 }{ \epsilon^2 } }  +\left (| \tex tbf{x}|  + \frac{\epsilon^ 2 }{1 2|\textbf{x}| }\r ight ) \ mb ox{ e rf } \le ft (  \f r a c{\sqrt{6}|\tex tbf{x}|}{\ ep s il on}\right) \ mbo x{ ,}$ wit h erf t he er r or func tion. The  small-sc al e cu t - off  is chosen  as $\ep silon=2dx $  andt he  larg e o ne as$L =1/ 2$. T he ker n els 
*$ stands_for the_convolution product), where $\theta$_is a_radially_symmetrical function,_such_that $\int_{\mathbb R_^3} \theta (\textbf{x})d\textbf{x}=1$_and $ \theta_\epsilon (\textbf{x})=_\frac{1}{\epsilon^3}\theta (\textbf{x}/\epsilon)$[@RobVar08].

Numerical study
===============

We_would_like now to study the statistical properties of the velocity field defined by Eq._(\[eq:K41Field\])._Analytical formulas_are_difficult_to obtain, thus we will_focus on numerical simulations. To_do so,_one has to choose the short time scale_$t=\tau$._It is easy_to check that the variance of the matrix $\textbf{S}_0$_(defined as $\langle \mbox{tr } \textbf{S}_0^2\rangle$)_goes to infinity_as_the_small scale parameter $\epsilon$_goes to zero. So we take_for $\tau$ the local normalizing value_$\tau =  (\mbox{tr } \textbf{S}_0^2)^{-1/2}$, in_the spirit of the RFD closures_provided in Ref. [@CheA2].

The simulation_is performed_in a 1-periodic box with_$N^3$ collocation points._The infinitesimal_volume is given_by $dV=dx^3$, with $dx=1/N$. We choose_as a regularizing_function and large-scale cut-off the isotropic_normalized_Gaussian function $\varphi_{L}_(\textbf{x})_=_\theta_{L}(\textbf{x}) =_\left(\frac{6}{\pi L^2}\right)^{3/2}\exp (-6|\textbf{x}|^2/L^2)$._This_allows to_compute_analytically the regularized norm of a_vector_$\textbf{x}$, useful for numerical purposes: $|\textbf{x}|_\epsilon =_\frac{\epsilon}{\sqrt{6\pi}}e^{-6\frac{|\textbf{x}|^2}{\epsilon^2}} + \left(|\textbf{x}| +_\frac{\epsilon^2}{12|\textbf{x}|}\right)\mbox{erf}\left(_\frac{\sqrt{6}|\textbf{x}|}{\epsilon}\right)\mbox{,}$ with erf the_error function. The small-scale cut-off_is chosen as $\epsilon=2dx$ and the_large one_as $L=1/2$._The kernels