Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: nonlinearities prevail in fluid dynamics when the reynolds number , @xmath4 , is large @xcite . for geophysical flows , the reynolds number is often larger than @xmath5 and for some astrophysical flows values of @xmath6 is not unreasonable . the number of degrees of freedom ( _ dof _ ) in the flow increases as @xmath7 for @xmath8 in the kolmogorov framework @xcite ( hereafter k41 ) . such a huge number of _ dof _ makes direct numerical simulations ( dns ) of turbulence at high @xmath4 infeasible on any existing or projected computer for decades to come . because of this intractability , simulations of turbulence are always carried out in regions of parameter space far from the observed values , either with : ( a ) an unphysical lack of scale separation between the energy - containing , inertial , and dissipative ranges while parameterizing the missing physics , or ( b ) a study of the processes at much smaller length scales , often with periodic boundaries ( unphysical at large scales but used under the hypothesis of homogeneity of turbulent flows ) . clearly , modeling of unresolved small scales is necessary . two main approaches have been developed over the years to model the effects of the unresolvable small scales in turbulence on the scales resolved in the simulations . the first approach is large eddy simulations ( les , see @xcite ) . les is widely used in engineering , in atmospheric sciences , and to a lesser extent in astrophysics . however , in the les approach , the reynolds number is not known . instead , one attempts modeling the behavior of the flow in the limit of very large @xmath4 . another approach models the effects of turbulence at higher reynolds numbers than are possible with a dns on a given grid , by using a variety of techniques that can be viewed as filtering of the small scales    a rather novel approach to modeling of turbulent flows employs _ regularization modeling _ as a sgs model @xcite . unlike closures which employ eddy - viscosity concepts ( modifying the dissipative processes ) , the approach of regularization modeling modifies the spectral distribution of energy . existence and uniqueness of smooth solutions can be rigorously proven , unlike many les models ( e.g. , eddy - viscosity ) , as well as the fact that the subgrid model recovers the navier - stokes equations in the limit of the filter width going to zero . their robust analytical properties ensure computability of solutions . proved the existence ( but not smoothness , or uniqueness ) of solutions to the navier - stokes equations in @xmath9 , n=2,3 using the leray model @xcite .    geurts and holm @xcite began using the leray model with ( a three - point invertible approximation of ) an inverse - helmholtz - operator filter of width @xmath2 . later it was dubbed leray@xmath0 and an upper bound for the dimension of the global attractor was established @xcite . the global existence and uniqueness of strong solutions for the leray model is a classical result @xcite . leray@xmath0 has been compared to dns simulations on a grid of @xmath10 in a doubly - periodic compressible channel flow domain @xcite . its performance was found to be superior to a dynamic mixed ( similarity plus eddy - viscosity ) model ( with an even greater reduction in computational cost ) . however , it possessed a systematic error of a slight over - prediction of the large scales accompanied by a slight under - prediction of the small scales . it possessed both forward- and back - scatter , but exhibited too little dissipation . the leonard tensor - diffusivity model @xcite ( sometimes known as the clark model @xcite ) is the first term of the reconstruction series for the turbulent sub - filter stress for all symmetric filters that possess a finite nonzero second moment . this leading order approximation of the subgrid stress is thus generic @xcite . in _ a priori _ testing , it reconstructs a significant fraction ( @xmath11 , but not all ) of the subgrid stress , provides for local backscatter along the stretching directions while remaining globally dissipative , and possesses a better reconstruction of the subgrid stress than the scale - similarity model . used as a les , in _ a posteriori _ testing the leonard tensor - diffusivity model required additional dissipation ( a dynamic smagorinsky term ) to achieve reasonable gains in computation speed for 3d periodic flows and for channel flows @xcite . the leonard tensor - diffusivity model does not conserve energy in the non - viscous limit . cao , holm , and titi @xcite developed a related ( conservative ) subgrid model which they dubbed clark@xmath0 . the clark@xmath0 model applies an additional inverse - helmholtz filter operation to the reynolds stress tensor of the clark model . the global well - posedness of the clark@xmath0 model and the existence and uniqueness of its solutions were demonstrated , and upper bounds for the hausdorff ( @xmath12 ) and fractal ( @xmath13 ) dimensions of the global attractor were found @xcite .    the third regularization model we will consider is the incompressible lagrangian - averaged navier - stokes ( lans@xmath0 , @xmath14model , also known as the viscous camassa - holm equation @xcite ) . it can be derived , for instance , by applying temporal averaging to hamilton s principle , where taylor s frozen - in turbulence hypothesis ( the only approximation in the derivation ) is applied as the closure for the eulerian fluctuation velocity in the reynold decomposition , at linear order in the generalized lagrangian mean description @xcite . in this derivation , the momentum - conservation structure of the equations is retained . for scales smaller than the filter width , lans@xmath0 reduces the steepness in gradients of the lagrangian mean velocity and thereby limits how thin the vortex tubes may become as they are transported , while the effect on larger length scales is negligible @xcite . lans@xmath0 may also be derived by smoothing the transport velocity of a material loop in kelvin s circulation theorem @xcite . consequently , there is no attenuation of resolved circulation , which is important for many engineering and geophysical flows where accurate prediction of circulation is highly desirable . an alternative interpretation of the @xmath14model is that it neglects fluctuations in the smoothed velocity field , while preserving them in the source term , the vorticity @xcite . lans@xmath0 has previously been compared to direct numerical simulations ( dns ) of the navier - stokes equations at modest taylor reynolds numbers ( @xmath15 @xcite , @xmath16 @xcite , and @xmath17 @xcite ) . lans@xmath0 was compared to a dynamic eddy - viscosity les in 3d isotropic turbulence under two different forcing functions ( for @xmath18 and @xmath19 ) and for decaying turbulence with initial conditions peaked at a low wavenumber ( with @xmath20 ) as well as at a moderate wavenumber ( with @xmath21 ) @xcite . lans@xmath0 was preferable in these comparisons because it demonstrated correct alignment between eigenvectors of the subgrid stress tensor and the eigenvectors of the resolved stress tensor and vorticity vector . the lans@xmath0 and the leray@xmath0 dynamic mixed model ( similarity plus eddy - viscosity ) in a turbulent mixing shear layer ( with @xmath22 ) @xcite . lans@xmath0 was found to be the most accurate of these three les candidates at proper subgrid resolution , but the effects of numerical contamination can be strong enough to lose most of this potential . while lans@xmath0 has the greatest grid - independent accuracy of the three models , it also requires the greatest resolution . from the les perspective , this could pose some limitations on the practical use and application of lans@xmath0 for high @xmath4 cases . indeed , recent high - resolution simulations of lans@xmath0 showed that energy artificially accumulates in the sub - filter - scales , giving as a result only a modest computational gain at very high reynolds number @xcite . we propose to pursue these previous studies of leray@xmath0 and lans@xmath0 further at high reynolds number , one goal is to as the three regularizations are related via truncation of sub - filter stresses , for lans@xmath0 , the predicted -scale spectra is @xmath23 @xcite . this scaling has been observed to be subdominant to an energy spectrum @xmath3 which corresponds to `` enslaved rigid body '' or `` polymerized '' portions of the fluid @xcite . the sub - filter scaling observed in the third - order structure function corresponded to the predicted @xmath24 scaling of the energy spectrum . however , regions were observed in the flow where no stretching was acting in the sub - filter scales . these regions , which give no contribution to the energy cascade , and hence do not affect the third order structure functions , are responsible for the @xmath3 scaling in the lans@xmath0 energy spectrum . for clark@xmath0 , the correct time scale for vortex stretching is difficult to determine and its spectrum is found to range between @xmath23 and @xmath25 @xcite . leray@xmath0 has the same difficulty and the spectrum can range between @xmath26 and @xmath27 @xcite . the determination of these scaling laws is needed to quantify the computational gain if each model is to be used as a sgs model . as a result , to determine empirically the sub - filter scale spectra . we present the three models and describe how they are related , derive the krmn - howarth equation for clark@xmath0 and leray@xmath0 ( from which exact scaling laws for third order quantities follow ) , and review theoretical predictions of inertial range scaling in section [ sec : models ] . in section [ sec : sfsstudy ] . we first compute a fully resolved dns of the navier - stokes equations at a resolution of @xmath28 ( @xmath29 , @xmath30 , and @xmath31 ) . we then perform model runs with the exact same conditions at a resolution of @xmath32 . we take @xmath2 to be @xmath33th the box size , which was found in an earlier study to be large enough to exhibit both navier - stokes and sub - filter - scale lans@xmath0 dynamics @xcite . this large filter case is important because it gives insight into the behavior of the models at scales much smaller than the filter width without requiring higher resolution than is feasible . we compare the three regularizations as subgrid models guided by a previous study of lans@xmath0 @xcite , we take @xmath2 to be @xmath34th the box size . this choice was found to produce an optimal @xmath14les ( in the sense of being optimal for the class of lans@xmath0 models , with respect to the value of @xmath2 ) . finally , we review bounds on the size of the attractors and use these bounds to comment on the computational savings of the three regularizations viewed as sgs models . the incompressible navier - stokes equations are given in cartesian coordinates by @xmath35 filtering these equations with a convolution filter , @xmath36 in which @xmath37 ( resp . @xmath38 ) denotes the filtered ( resp . unfiltered ) field , yields @xmath39 in which by convention we denote @xmath40 and the reynolds turbulence stress tensor , @xmath41 , represents the closure problem . it has been shown @xcite that for all symmetric filters possessing a finite nonzero second moment , the first term of the reconstruction series for the turbulent sub - filter stress is @xmath42 where @xmath43 is the fourier transform of the convolution kernel ( @xmath44 is the convolution kernel where @xmath45({\mathbf}{r})=\int g({\mathbf}{r}-{\mathbf}{r { ' } } ) z({\mathbf}{r { ' } } ) d^3{\mathbf}{r{'}}$ ] ) . this approximation of the subgrid stress is then generic and is known as the leonard tensor - diffusivity model @xcite ( or , often , the clark model @xcite ) .    related to this model , the clark@xmath0 model is @xcite @xmath46 and its subgrid stress for @xmath47 is given by @xmath48 here the filter is the inverse of a helmholtz operator , @xmath49 . the clark@xmath0 model conserves energy in the @xmath50 norm instead of the @xmath51 norm , @xmath52 in which the clark@xmath0 energy @xmath53 is expressed as @xmath54 and the clark@xmath0 energy dissipation rate is given by @xmath55 where @xmath56 and @xmath57 . for @xmath58 , we note that @xmath59 which implies that the turbulent sub - filter stress tensor for the tensor - diffusivity model given by eq . ( [ eq : carati ] ) is @xmath60 which is proportional to the clark@xmath0 stress tensor to second order in @xmath2 . hence , the _ a priori _ tests of @xcite should apply to clark@xmath0 , at least at this order .      in 1938 , krmn and howarth @xcite introduced the invariant theory of isotropic hydrodynamic turbulence , and derived from the navier - stokes equations the exact law relating the time derivative of the two - point velocity correlation to the divergence of the third - order correlation function . the corresponding krmn - howarth theorem for lans@xmath0 in the fluid case was derived in @xcite . the relevance of the krmn - howarth theorem for the study of turbulence can not be underestimated . as a corollary , rigorous scaling laws in the inertial range can be deduced . in this section , we derive these results for the clark@xmath0 case .    for the sake of simplicity , we consider the case @xmath61 , since the dissipative terms may be added at any point in the derivation . we denote @xmath62 and begin our investigation of the correlation dynamics by computing the ingredients of the partial derivative @xmath63 . the clark@xmath0 motion equation ( [ eq : clark ] ) may be rewritten as @xmath64 combining eqs . ( [ eq : subgrid ] ) and ( [ eq : sgs_clark ] ) , we arrive at the fluctuation - velocity equation , @xmath65 where @xmath66 ( @xmath67 ) . multiplying eq . ( [ eq : clarkv2 ] ) by @xmath68 and eq . ( [ eq : clarku2 ] ) by @xmath69 , then adding the result yields @xmath70 where we have used statistical homogeneity @xmath71 we symmetrize eq . ( [ fluct - eqn ] ) in the indices @xmath72 by adding the corresponding equation for @xmath73 . we then use homogeneity again as @xmath74 and define the tensors @xmath75 we can drop @xmath76 because the terms with the pressures @xmath77 and @xmath78 vanish everywhere , as follows from the arguments of isotropy @xcite . finally , we obtain @xmath79 this is the krmn - howarth equation for clark@xmath0 ( compare to eq . ( 3.8 ) in @xcite for lans@xmath0 ) .    by dimensional analysis , the energy dissipation rate in clark@xmath0 is @xmath80 and eq . ( [ eq : kh1 ] ) implies @xmath81 for large scales ( @xmath82 ) , we recover the navier - stokes scaling known as the four - fifths law , @xmath83 @xcite . here , @xmath84\cdot{\mathbf}{l}/l$ ] is the longitudinal increment of @xmath85 . strictly speaking , the four - fifths law expresses that the third - order longitudinal structure function of @xmath85 , @xmath86 , is given in the inertial range in terms of the mean energy dissipation per unit mass @xmath87 by @xmath88 or , equivalently , that the flux of energy across scales in the inertial range is constant . we also recover the kolmogorov 1941 @xcite ( hereafter , k41 ) energy spectrum , @xmath89 or , equivalently , @xmath90 for sub - filter scales ( @xmath91 ) , we have @xmath92 and the first and third right - hand terms in eq . ( [ eq : full_clark ] ) are equivalent . in this case , we are left with two different possible scalings depending on the prefactors in eq . ( [ eq : full_clark ] ) . if the first ( or third ) right - hand term is dominant , our scaling law becomes @xmath93 for our sub - filter scale energy spectrum we would then have @xmath94 , or , equivalently , @xmath95 this result is the same as for the @xmath14model @xcite . if , however , the second right - hand term in eq . ( [ eq : full_clark ] ) is dominant , then the k41 results are recovered , with @xmath96 substituted for @xmath85 . in that case , one finds the alternative clark@xmath0 sub - filter scale spectral energy scaling , @xmath97      we review here the derivation by dimensional analysis of the spectrum which follows the scaling ideas originally due to kraichnan @xcite and which is developed more fully in ref . @xcite . in examining the nonlinear terms in eq . ( [ eq : clark ] ) , it is not entirely clear which of three possible scales for the average velocity for an eddy of size @xmath1 , @xmath98 @xmath99 or @xmath100 should result . therefore , three corresponding `` turnover times '' , @xmath101 , for such an eddy may be proposed @xmath102 the term `` turnover time '' is used advisedly here , since only the velocity @xmath103 is composed of the fluid transport velocity . we define the ( omnidirectional ) spectral energy density , @xmath104 , from the relation @xmath105 since , @xmath106 , we have @xmath107 then , the total energy dissipation rate , @xmath108 is related to the spectral energy density by @xmath109 which yields , finally , the predicted energy spectra for clark@xmath0 , @xmath110 , @xmath111 for scales much larger than @xmath2 ( @xmath112 ) the kolmogorov scaling for navier - stokes is recovered , @xmath113 whereas for scales much smaller ( @xmath114 ) , the spectrum becomes @xmath115 these arguments constrain the clark@xmath0 sub - filter scale spectrum to lie between @xmath1 and @xmath116 . the leray model in cartesian coordinates is @xmath117 where the flow is advected by a smoothed velocity , @xmath96 . by comparison with eq . ( [ eq : subgrid ] ) we see that the leray model approximates the subgrid stress as @xmath118 , or , with @xmath119 , @xmath120 as has been noted previously @xcite , the subgrid stress of clark@xmath0 in eq . ( [ eq : sgs_clark ] ) is a truncation of the subgrid stress of leray@xmath0 , in eq . ( [ eq : sgs_leray ] ) . for leray@xmath0 , the @xmath121 norm is the quadratic invariant that is identified with energy , @xmath122 where @xmath123 and @xmath124 as was pointed out in @xcite , the incompressibility of the velocity field @xmath85 only implies a divergenceless filtered velocity @xmath96 under certain boundary conditions for leray@xmath0 . when @xmath125 , the energy @xmath126 is no longer conserved ( helicity and kelvin s theorem are not conserved for leray@xmath0 ) . in our numerical study , we employ periodic boundary conditions , for which @xmath127 implies @xmath128 and leray@xmath0 conserves energy in the usual sense of @xmath51 .      in this section we derive the krmn - howarth equation for the leray@xmath0 case . following section [ sec : kh - clark ] , we begin our investigation of the correlation dynamics by computing the ingredients of the partial derivative @xmath129 . ( [ eq : leray ] ) may be rewritten as : @xmath130 multiplying eq . ( [ eq : lerayv2 ] ) by @xmath131 yields @xmath132 we can make this equation symmetric in the indices @xmath72 , by adding the equation for @xmath133 . we define the tensors @xmath134 again , we may drop @xmath135 because the terms with the pressure @xmath136 vanish everywhere and thereby obtain @xmath137 this is the krmn - howarth equation for leray@xmath0 .    the energy dissipation rate for leray@xmath0 is denoted by @xmath138 and it satisfies @xmath139 . by dimensional analysis , eq . ( [ eq : kh1leray ] ) implies @xmath140 for large scales ( @xmath82 ) , we recover the navier - stokes scaling , eqs . ( [ eq : fourfifths ] ) and ( [ eq : k41 ] ) . for sub - filter scales ( @xmath91 ) our scaling law becomes @xmath141 for our small scale energy spectrum we would then have @xmath142 ( where we employed @xmath92 ) , or , equivalently , cf . ( [ eq : clarka_spectrum2 ] ) , @xmath143      we review here the derivation by dimensional analysis of the spectrum for leray@xmath0 as we did for clark@xmath0 in section [ sec : pheno_clark ] . this analysis is developed more fully in ref . we argue again that there are three possible scales for the average velocity for an eddy of size @xmath1 , eqs . ( [ eq : uk0 ] ) , ( [ eq : uk1 ] ) , and ( [ eq : uk2 ] ) , with the turn - over time , @xmath144 given by eq . ( [ eq : turnover ] ) . since , @xmath145 , we have @xmath146 then , the total energy dissipation rate , @xmath138 is related to the spectral energy density by @xmath147 which yields , finally , the predicted energy spectra for leray@xmath0 , @xmath148 , @xmath149 for scales much larger than @xmath2 ( @xmath112 ) the k41 spectrum is recovered , eq . ( [ eq : k41 ] ) , and for scales much smaller ( @xmath114 ) the spectrum is @xmath150 these arguments constrain the leray@xmath0 sub - filter scale spectrum to lie between @xmath151 and @xmath152 . lans@xmath0 is given by @xmath153 for lans@xmath0 , the usual choice of filter is again @xmath154 with this filter , the subgrid stress tensor is given by @xmath155 as has been previously noted @xcite , the subgrid stress of leray@xmath0 , eq . ( [ eq : sgs_leray ] ) , is a truncation of the subgrid stress of lans@xmath0 eq . ( [ eq : sgs_lans ] ) . like clark@xmath0 , energy is conserved in the @xmath156 norm instead of the @xmath51 norm . additionally , lans@xmath0 is the only model of the three examined here that conserves the helicity ( and kelvin s circulation theorem ) .    for lans@xmath0 in the fluid case the krmn - howarth theorem was derived in @xcite . we summarize here the dimensional analysis argument for the lans@xmath0 inertial range scaling that follows from this theorem , beginning from equation ( 3.8 ) in ref . @xcite . in the statistically isotropic and homogeneous case , without external forces and with @xmath157 , taking the dot product of eq . ( [ eq : lans ] ) with @xmath158 yields the equation    @xmath159    the trace of this equation is the fourier transform of the detailed energy balance for lans@xmath0 ; @xmath160 is the second - order correlation tensor while @xmath161 and @xmath162 are the third - order correlation tensors for lans@xmath0 and @xmath163 is the sub - filter scale stress tensor . for @xmath164 this reduces to the well - known relation derived by krmn and howarth . the energy dissipation rate for lans@xmath0 , @xmath165 , satisfies @xmath166 . by dimensional analysis in eq . ( [ eq : holm38 ] ) we arrive at @xmath167 for large scales ( @xmath82 ) , we recover the navier - stokes scaling eqs . ( [ eq : fourfifths ] ) and ( [ eq : k41 ] ) . for sub - filter scales ( @xmath91 ) our scaling law becomes eq . ( [ eq : lcube ] ) and our sub - filter scale spectra is given by @xmath168 in this case , by the phenomenological arguments , we know that eddies of size @xmath1 are advected by the smoothed velocity , eq . ( [ eq : uk2 ] ) . this scaling is confirmed in ref . but it coexists with a @xmath169 energy spectrum corresponding to `` enslaved rigid bodies '' or `` polymerized '' portions of fluid which do not contribute to the turbulent energy cascade . to this end , we compute numerical solutions to eqs . ( [ eq : ns ] ) , ( [ eq : clark ] ) , ( [ eq : leray ] ) , and ( [ eq : lans ] ) in a three - dimensional ( 3d ) cube with periodic boundary conditions using a parallel pseudospectral code @xcite . we employ a taylor - green forcing @xcite , @xmath170 \label{eq : tg}\ ] ] ( with @xmath171 ) , and employ dynamic control @xcite to maintain a nearly constant energy with time . the taylor - green forcing , eq . ( [ eq : tg ] ) , is not a solution of the euler s equations , and as a result small scales are generated rapidly . the resulting flow models the fluid between counter - rotating cylinders @xcite and it has been widely used to study turbulence , including studies in the context of the generation of magnetic fields through dynamo instability @xcite . we define the taylor microscale as @xmath172 and the mean velocity fluctuation as @xmath173 . the taylor microscale reynolds number is defined by @xmath174 and the reynolds number based on a unit length is @xmath175 .    to compare the three regularizations ( clark@xmath0 , leray@xmath0 , and lans@xmath0 ) we compute a fully resolved dns of the navier - stokes equations at a resolution of @xmath28 ( @xmath29 , @xmath176 ) and model runs with the exact same conditions at a resolution of @xmath32 . the details of the flow dynamics of the dns have already been given @xcite . in particular , the reynolds number based on the integral scale @xmath177 ( where @xmath178 is the total energy ) is @xmath179 , where @xmath180 is the _ r.m.s . _ velocity and the reynolds number based on the taylor scale is @xmath31 . the dns was run for nine turnover times @xmath181 ( in the following results , time @xmath182 is in units of the turnover time ) . we employ a filter width of @xmath183 for which lans@xmath0 exhibits both navier - stokes and lans@xmath0 inertial ranges in the third - order structure function @xcite . from these we hope to obtain the behavior of the models for scales much smaller than @xmath2     for leray@xmath0 and in the dns , @xmath184 for lans@xmath0 and clark@xmath0 ) . dns ( @xmath176 ) is shown as solid black lines , lans@xmath0 as dotted red , clark@xmath0 as dashed green , and leray@xmath0 as blue dash - dotted . the cyan long - dashed line represents a @xmath185 dns ( @xmath186 , @xmath187 ) . here each run is calculated only until it reaches a statistical steady state . leray@xmath0 reduces the dissipation , @xmath188 , and increases the time scale to reach a statistical turbulent steady - state . both effects are greater as @xmath2 is increased . by comparison with the @xmath189 run , we see that these two effects are consistent with a reduced effective reynolds number . a smaller reduction in flux ( but not an increase in time to steady - state ) is also observed for lans@xmath0 and is likely related to its rigid bodies.,width=338 ]    in fig . [ fig : enstrophy ] we present the time evolution of the enstrophy ( @xmath190 for leray@xmath0 and dns , @xmath184 for lans@xmath0 and clark@xmath0 ) , which is proportional to the dissipation ( @xmath191 or @xmath192 depending on the case ) . also shown is a well - resolved @xmath185 dns of a less turbulent flow ( @xmath193 , @xmath194 , @xmath187 , ( cyan online ) long - dashed line ) . here each run is calculated only until it reaches a statistically steady state . we see that the dissipation is reduced and the time - scale is increased for leray@xmath0 . we see by comparison with the @xmath194 dns that a reduced effective reynolds number in the leray@xmath0 run . for lans@xmath0 , the dissipation is decreased although the time to reach steady state is not increased . this is probably related to the enslavement of its rigid body regions which would have no internal dissipation . of the three models , clark@xmath0 the total dissipation for a large range of @xmath2 . indeed , as it is the order @xmath195 approximation of navier - stokes , this dissipation behavior for clark@xmath0 may continue to hold until @xmath2 becomes quite large . figure [ fig : eavg ] shows a comparison of the energy spectrum at the turbulent steady state , for all the runs in fig . [ fig : enstrophy ] . the isotropic energy spectra are calculated as follows , @xmath196 where @xmath197 ( the @xmath198 norm is employed for clark@xmath0 and lans@xmath0 ) . the length scale @xmath2 is indicated by a vertical dashed line and the plotted energy spectra are compensated by @xmath199 ( i.e. , leading to a flat k41 @xmath152 spectrum ) . from fig . [ fig : eavg ] it is clear that the clark@xmath0 spectral behavior is close to the predicted @xmath1 spectrum , rather than the @xmath200 from eq . ( [ eq : clarka_spectrum2 ] ) , or the other possible spectrum from eq . ( [ eq : clarka_spectrumn ] ) . likewise , the @xcite . leray@xmath0 , on the other hand , possesses a very steep sub - filter scale spectrum as well as enhanced large - scale energy as has been previously observed @xcite .    the spectrum of leray@xmath0 in fig . [ fig : eavg ] gives a good approximation to the @xmath189 dns in the range @xmath201 $ ] ( i.e. , to @xmath193 rather than to @xmath202 which was employed ) . this result , leray@xmath0 s increased characteristic time scales and its reduced dissipation imply that the leray@xmath0 model is operating at a much lower effective reynolds number . this is also clear from the rapid drop in the spectrum at small scales , shown in fig . [ fig : eavg ] . indeed , we can build an effective reynolds number in the large scales as @xmath203 . since @xmath204 is controlled in this simulation by the forcing scale , the drop in the dissipation rate implies a reduced nonlinearity in leray@xmath0 . this is also consistent with a direct comparison of the nonlinear terms in leray@xmath0 with , for instance , lans@xmath0 . the nonlinear terms in lans@xmath0 , eq . ( [ eq : lans ] ) , may be written as @xmath205 ( where the suffix @xmath206 denotes a transposition ) , while the nonlinear term in leray@xmath0 , eq . ( [ eq : leray ] ) , is only @xmath207 . both nonlinear terms in lans@xmath0 are of order @xmath208 ; so the absence of one of the nonlinear terms in leray@xmath0 be understood as a reduction in the nonlinearity . validation of the krmn - howarth equation scalings , eqs . ( [ eq : lcube ] ) and ( [ eq : lcubeleray ] ) , enables us to measure scaling laws in the inertial range and , thus , compare the intermittency properties of the models . the third - order correlations involved in the theorems , namely @xmath209 for clark@xmath0 and lans@xmath0 , and @xmath210 for leray@xmath0 , and @xmath211 for navier - stokes , are plotted versus @xmath212 in fig . [ fig : theorem ] . in fig . [ fig : theorem ] we can see validation of the krmn - howarth scaling for scales smaller than @xmath2 for both lans@xmath0 and clark@xmath0 . in particular , we note the observed scaling for clark@xmath0 verifies the @xmath213 scaling and not the ( theoretically possible ) @xmath214 ( @xmath215 ) scaling . the predicted scaling is not observed in leray due to its reduced effective reynolds number . with these scalings in hand , we may proceed to observe the scaling of the longitudinal structure functions , @xmath216 where we again replace the @xmath217 norm , @xmath218 , for the @xmath219 norm , @xmath220 , in the case of clark@xmath0 and lans@xmath0 . we utilize the extended self - similarity ( ess ) hypothesis @xcite which proposes the scaling @xmath221 and normalize the results by @xmath222 to better visualize the deviation from linearity ( which serves as a measure of intermittency ) . the results are displayed in fig . [ fig : exponents ] .    in fig . [ fig : exponents ] , we may observe the intermittency properties of the models at sub - filter scales . we note a reduced intermittency for both leray@xmath0 and the @xmath194 dns . this is consistent with the smoother , more laminar fields ( due to the reduction of the effective @xmath4 ) possessed by both . interestingly , though lans@xmath0 and clark@xmath0 both possess the same cascade scaling ( eq . ( [ eq : lcube ] ) , as confirmed in fig . [ fig : theorem ] ) , the clark@xmath0 model is markedly more intermittent than lans@xmath0 . if artificially truncated local interactions ( in spectral space ) is taken as a cause of enhanced intermittency @xcite , then the increased intermittency observed in clark@xmath0 is the expected result of truncation of the higher - order terms in the sub - filter stress tensor . moreover , if the lans@xmath0 s @xmath3 spectrum is indeed associated with rigid bodies , these would serve to decrease the intermittency ( no internal degrees of freedom being available in a rigid body ) which is consistent with the results shown here . due to this effect , lans@xmath0 of the three regularization models the high - order intermittency of @xmath223 corresponding to an optimal @xmath14les @xcite . note that the value of @xmath2 has been optimized for neither clark@xmath0 nor leray@xmath0 ; as a consequence , these models might perform better in other parameter regimes than the results indicate in this study .     for leray@xmath0 and in the dns , @xmath184 for lans@xmath0 and clark@xmath0 ) . dns ( @xmath176 ) is shown as solid black lines , lans@xmath0 as dotted red , clark@xmath0 as dashed green , and leray@xmath0 as blue dash - dotted . an under - resolved ( @xmath32 ) navier - stokes run is shown as a pink dash - triple - dotted line.,width=338 ]    fig . [ fig : enstrophy2 ] along with that of an under - resolved navier - stokes solution at a resolution of @xmath32 ( @xmath202 , ( pink online ) dash - triple - dotted line ) . we see that both lans@xmath0 and clark@xmath0 reproduce the proper amount of dissipation and are within @xmath224 of the time required by the dns to reach a statistical turbulent steady state . as has been observed before , leray@xmath0 is under - dissipative @xcite . we also note that it takes longer than the    compensated spectra averaged over several eddy turn - over times are shown for the sgs case ( i.e. , @xmath225 ) in fig . [ fig : eavg2 ] . note that as the subgrid models are averaged over a different time interval , no meaningful comparison to the dns is possible for @xmath226 . even without an optimal choice for the value of @xmath2 , clark@xmath0 best reproduces the dns spectrum for scales larger than @xmath2 . we compute root - mean - square spectral errors where @xmath227 is the wavenumber for the forcing scale , @xmath228 is the dns spectrum ( in the @xmath121 norm ) , and @xmath229 is the subgrid model spectrum ( in the appropriate norm ) . we see that only clark@xmath0 improves the estimate of the power spectrum at this resolution for leray@xmath0 performs the poorest of the three regularization models , but it is also not optimized . as previously argued , its effective reynolds number is too low to accurately model the dns flow . either a decrease in the viscosity @xmath230 , or a decrease in the filter size @xmath2 ( and , hence , an increase in the nonlinearity ) , or both would likely improve the accuracy of leray@xmath0 as an sgs model . due to its frozen - in ( or enslaved ) rigid - body regions , the lans@xmath0 model can not reproduce the dns spectrum unless @xmath2 is only a few times larger than the dissipation scale @xcite . [ cols="^,^,^ " , ]     another measure of the success of a subgrid model is the reproduction of structures in the flow . in figure [ fig : vapor ] we have 3d volume rendering of the enstrophy density @xmath231 ( @xmath232 for lans@xmath0 and clark@xmath0 ) for the dns , the three sgs - model simulations ( @xmath233 ) , the @xmath32 under - resolved navier - stokes solution , all at a reynolds number of @xmath234 , and the @xmath194 dns . due to the late times depicted ( longer than a lyapunov time ) there can be no point - by - point comparison between the simulations . instead , we note that there are four horizontal bands where the forcing causes a maximum shear . this large - scale feature of the flow is missing only from leray@xmath0 and the @xmath194 run . the three other runs reproduce this feature well ( note that the apparently thicker tubes present in clark@xmath0 are vortex tube mergers ) . the results lead again to the conclusion that the under - resolved navier - stokes , the clark@xmath0 , and the lans@xmath0 models are better subgrid models than leray@xmath0 due to its reduced effective @xmath4 .    for the sgs models , the predicted @xmath235 from the krmn - howarth theorem for navier - stokes is well - reproduced by all models at super - filter scales . in fig . [ fig : exponents2 ] the sgs model intermittency results . all models reproduce the intermittency up to the tenth - order moment within the error bars ( although there is a small decrease in intermittency for leray@xmath0 ) . thus , we conclude that with adequately chosen values of @xmath2 ( and of @xmath230 for leray@xmath0 ) , all three models can reproduce the intermittency of the dns ( to within the error bars ) . the rationale behind using is that it leads to adequate solutions at a reduced computational cost , because it computes fewer _ dof _ ; indeed , for an , the ratio of navier - stokes s _ dof _ to the model s _ dof _ , a prediction for memory savings and hence computation time savings for numerical simulation , is a crucial factor . consequently , analytical bounds on the sizes of the attractors for the three regularization subgrid models may be useful indicators of their computational savings . the _ dof _ for lans@xmath0 is derived in @xcite and confirmed in @xcite , @xmath236 where @xmath237 is the integral scale ( or domain size ) . we may compare this to the _ dof _ for navier stokes , @xmath238 which immediately yields @xmath239 it was found , however , that to reproduce the energy spectrum of an equivalent dns , the filter - width @xmath2 must be no larger than a few times the dissipation scale , @xmath240 @xcite . this is the result of the `` polymerization '' of the flow in lans@xmath0 , and the associated @xmath241 scaling at sub - filter scales . with this added caveat , it follows that the reduction in _ dof _ is independent of @xmath4 ( and a net factor of about @xmath242 ) . our study here illustrates that the high - order structure functions may be reproduced for much larger values of @xmath2 . therefore , in applications where the spectrum is not of great concern , much greater reduction in numerical resolution would be feasible . for clark@xmath0 there is an upper bound on the hausdorff , @xmath12 , and fractal , @xmath13 , dimensions of the attractor , @xmath243 where @xmath244 is the kolmogorov dissipation length scale corresponding to the clark@xmath0 model @xcite . from its observed @xmath1 spectrum , we may estimate @xmath244 or , equivalently , @xmath245 . for dissipation the large wavenumbers dominate and , therefore , combining the clark@xmath0 energy balance eq . ( [ eq : lansa_balance ] ) with its sub - filter scale energy spectrum eq . ( [ eq : clarka_spectrum ] ) allows us to implicitly specify its dissipation wavenumber , @xmath246 , by @xmath247 then we have , @xmath248 it follows that @xmath249 this is similar to the prediction for lans@xmath0 , but as energy spectra are more easily reproduced for larger values of @xmath2 than with lans@xmath0 ( but not the intermittency properties ) , it may be the case that @xmath2 is not tied to the kolmogorov dissipation scale @xmath240 . if so , then the computational saving might increase as @xmath250 which is promising for use of clark@xmath0as an les model . this conclusion is bolstered to the extent that the results in section [ sec : super ] for @xmath233 ( @xmath251 ) are acceptable . if , there is still a greater reduction in _ dof _ ( a factor of 20 ) for clark@xmath0 than for lans@xmath0 .    for leray@xmath0 , we have the following upper bounds on the hausdorff dimension , @xmath12 , and fractal dimension , @xmath13 , of the global attractor , @xmath252 where @xmath253 is the dissipation length scale for leray@xmath0 @xcite . again , we estimate the dissipation wavenumber for leray@xmath0 @xmath254 . then , from eqs . ( [ eq : leray_balance ] ) and ( [ eq : leray_spectrum ] ) , that is , assuming the @xmath151 spectrum resulting from the krmn - howarth equation , we find @xmath255 consequently we have , @xmath256 it follows that @xmath257 our results suggest that for an effective les the viscosity @xmath258 must be chosen to be _ smaller _ than @xmath230 . this leads to an upper bound on the computational savings for leray@xmath0 , @xmath259 if we further assume that @xmath2 is directly proportional to the dissipative scale @xmath240 , we arrive at @xmath260 which is not exceedingly promising for use as a les . all such estimates are , however , purely conjectural until the proper choices of @xmath2 and @xmath258 are determined . we derived the krmn - howarth equations for the leray@xmath0 and clark@xmath0 models . these two models may be viewed as successive truncations of the sub - filter scale stress of the lagrangian - averaged navier - stokes @xmath14model ( lans@xmath0 ) . in the case of clark@xmath0 two different inertial range scalings follow from the dimensional analysis of this equation . the case of leray@xmath0 is simpler as a single scaling is predicted . this is the case for navier - stokes and lans@xmath0 as well . solutions for a @xmath28 dns under periodic boundary conditions ( @xmath261 ) using a @xmath32 resolution under the same exact conditions for lans@xmath0 , leray@xmath0 , clark@xmath0 , and an under - resolved @xmath32 solution of the navier - stokes equations . we employed two different filter widths @xmath2 . the first choice @xmath262 was used to and the second choice @xmath263 was employed to test the sgs potential of the models . in comparing these two choices , we found for leray@xmath0 that an increase in @xmath2 substantially decreases the nonlinearity ( and hence decreases the effective reynolds number @xmath4 ) . for this reason , we were unable to confirm either the inertial range scaling from its krmn - howarth equation or its sub - filter scale energy spectrum . for clark@xmath0 we were able to determine the dominant krmn - howarth inertial range scaling to be @xmath264 associated @xmath1 energy spectrum ,    the performance of the three regularizations as sgs models ( for a resolution of @xmath32 and @xmath225 ) was comparable to that of the under - resolved navier - stokes solution in reproducing the dns energy spectrum . only clark@xmath0 showed a clear improvement in approximating the spectrum . from 3d volume rendering of enstrophy density we found that clark@xmath0 and lans@xmath0 were comparable to the under - resolved solution . even at @xmath265 , leray@xmath0 s 3d spatial structures are consistent with a significantly reduced @xmath4 flow ( e.g. , comparable to a @xmath266 dns ) . we note that the value of @xmath2 was chosen optimally for lans@xmath0 at the resolution of @xmath32 , and that for clark@xmath0 ( and especially for leray@xmath0 ) smaller resolutions may have comparable results for this value . such a comparison is beyond the scope of the present work . although lans@xmath0 and clark@xmath0 exhibit the same inertial range scaling arising from similarities in their krmn - howarth equations , clark@xmath0 is decidedly more intermittent than at the same time , lans@xmath0 is only slightly more intermittent than navier - stokes . these results are consistent with the artificial truncation of local nonlinear interactions ( in spectral space ) in the sgs stress tensor of each model . this effect is reduced for lans@xmath0 by the `` rigid - body regions '' enslaved in its larger scale flow which possess no internal degrees of freedom . the reduced intermittency observed for leray@xmath0 is related to its smoother , more laminar fields as a result of its reduced effective @xmath4 . finally , we analyzed the reduction in the number of _ dof _ in the models , as compared to navier - stokes ( and , hence , their les potential based on their computational savings ) . we noted that as lans@xmath0 reproduces the intermittency properties of a dns quite well even for larger values of @xmath2 , some further reduction in numerical saving might be achieved provided the contamination due to its @xmath169 rigid - body energy spectrum were not important in a given application . as clark@xmath0 possesses a similar reduction in _ dof _ to lans@xmath0 , its les potential is tied to the optimal value of @xmath2 for les . our study indicates that clark@xmath0 may be applicable ( especially with regards to the energy spectrum ) for larger values of @xmath2 than lans@xmath0 . in fact , if its optimal value is not a function of @xmath4 , the computational resolution savings increases as @xmath250 for clark@xmath0 . for the case of leray@xmath0 , the prediction is complicated by the effective reduction in @xmath4 as @xmath2 increases . prediction of optimized values of @xmath2 and of effective dissipation @xmath258 are required to assess its les potential . future work should include such a study for both leray@xmath0 and clark@xmath0 . all three regularizations were shown to be successful , in that their control of the flow gradient reduces the degrees of freedom and saves computation while preserving a properly defined reynolds number ( albeit for leray@xmath0 that definition is not yet ) . clark@xmath0 accurately reproduces the total dissipation , the time scale to obtain a turbulent statistical steady - state and the large - scale energy spectrum of a dns . these results seem to result from clark@xmath0 being an order @xmath195 approximation of navier - stokes . we have shown that leray@xmath0 reduces the effective reynolds number of the flow . the last of the three models , lans@xmath0 restores kelvin s circulation theorem ( advected by a smoothed velocity ) and the conservation of a form of helicity . using spectra as a measure of the success of a subgrid model , lans@xmath0 is less than optimal , due to its however , other measures of the success of a subgrid model are possible : for example , in regard to intermittency , lans@xmath0 may be considered a superior model . for clark@xmath0 , intermittency may be a function of filter width while for lans@xmath0 , intermittency does not vary much with @xmath2 . we remark that the computational gain thus far achieved by any of these regularizations is insufficient for applications at very high reynolds numbers , and the three subgrid stress tensors discussed here may need to be supplemented with an enhanced effective viscosity . this is a common practice when implementing the clark model ( see e.g. , ref . @xcite for a study of this model with an extra smagorinsky term ) . in this light , the present study may be useful as an analysis of the properties of the sgs tensors of the regularizations , and to pick best candidates , before the addition of enhanced dissipation . studies similar to that in ref . @xcite will also need to be done for the cases of clark@xmath0 and leray@xmath0 to quantify their computational savings before the addition of such dissipative terms . computer time provided at ncar is gratefully acknowledged , as well as nsf - cmg grant 0327888 . the work of ddh was partially supported by the us department of energy , office of science , advanced scientific computing research , and by the royal society of london wolfson award for meritorious research . a.  leonard . . in _ turbulent diffusion in environmental pollution ; proceedings of the second symposium , charlottesville , va . , april 8 - 14 , 1973 . volume a. ( a75 - 30951 13 - 47 ) new york , academic press , inc . , 1974 , p. 237 - 248 . nasa - supported research . _ , pages 237248 , 1974 . j.  pietarila graham , d.  holm , p.  mininni , and a.  pouquet . highly turbulent solutions of the lagrangian - averaged navier - stokes alpha model and their large - eddy - simulation potential . , 76:056310+ , 2007 .
Abstract:

Answer: and to a direct numerical simulation ( dns ) of the navier - stokes equations . we derive the krmn - howarth equation for both the clark@xmath0 and leray@xmath0 models . we confirm one of two possible scalings resulting from this equation for clark@xmath0 as well as its associated @xmath1 energy spectrum . clark@xmath0 reproduces the large - scale energy spectrum and intermittency properties of the dns . for the leray@xmath0 model , increasing the filter width , @xmath2 , decreases the nonlinearity and , hence , the effective reynolds number is substantially decreased . therefore even for the smallest value of @xmath2 studied leray@xmath0 was inadequate as a sgs model . the lans@xmath0 energy spectrum @xmath3 , consistent with its so - called `` rigid bodies , '' precludes a reproduction of the large - scale energy spectrum of the dns . we find , however , that this same feature reduces its intermittency compared to clark@xmath0 ( which shares a similar krmn - howarth equation ) . clark@xmath0 is found to be the best approximation for reproducing the total dissipation rate and the energy spectrum at scales larger than @xmath2 , whereas high - order intermittency properties for larger values of @xmath2 are best reproduced by lans@xmath0 .