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: extreme ultraviolet ( euv ) observations currently provide the most important source of information on the low solar corona . as new euv instruments are developed , the temporal , spatial and spectral resolution becomes ever finer , giving new insight on the coronal and chromospheric structure and dynamics . the _ atmospheric imaging assembly _ ( aia : ) aboard the _ solar dynamics observatory _ ( sdo : ) provides very fine temporal and spatial resolution of the sun at multiple wavelengths , and is having a large impact on the field . even as the community develop methods to digest the huge volume of data from aia / sdo , new instruments are planned and tested with even finer resolution ( e.g. _ high - resolution coronal imager _ ( hi - c ) , see ) .    despite the development of automated detection tools for euv observations ( e.g. ) , most scientific works begin from visual inspection of images . in particular , the volume of data provided by aia / sdo is so high that low - resolution images are often used as a starting point to find features of interest . the higher - resolution data are then used for further analysis . image processing is therefore an important step in analysing the data . the scientific return can be improved by the application of processing which better reveals features in the data , particularly in the early stages of analysis where visual inspection is most important . a common approach to process euv images is simply to display the square root ( or a gamma curve transformation ) , or alternatively the logarithm , of the original pixel values . this is a quick and easy way of reducing the dominance of the image contrast range by a few small bright regions . to reveal dynamic features , time - differencing is used , where a previous image is subtracted from the current image . such simple processes are commonly used not due to the quality of the output , but due to their simplicity . more advanced image processing methods are not so commonly used due to their complexity and computational expense . wavelet - based techniques have been used by and to greatly improve the visual information available from the _ extreme ultraviolet imaging telescope _ ( eit : ) aboard the _ solar and heliospheric observatory _ ( soho ) and the _ extreme ultraviolet imagers _ ( euvi : ) aboard the _ solar terrestial relations observatory _ ( stereo : ) . the technique involves the decomposition of images into different spatial scales , and the filtering / enhancements of features at the multiple resolutions . the technique is computationally expensive yet gives very good results . a less sophisticated yet very efficient technique to reveal features above the limb is based on techniques originally developed for coronagraphs @xcite . most relevant to this work is the recently - developed noise adaptive fuzzy equalization ( nafe ) method @xcite . the method is inspired by adaptive histogram equalization , where local statistics govern the output value of a pixel . the method uses a fuzzy membership function of a gaussian - weighted local set to enhance structural detail , preserve contextual detail , and to reduce noise . the nafe method results in very clear images without artifacts and with excellent noise reduction . its one downside is computational expense . this work summarizes the problems in visualising features in euv images ( section [ observation ] ) , introduces the new method ( section [ method ] ) , applies the method to several images from various instruments ( section [ results ] ) and closes with a brief summary ( section [ summary ] ) . an observation from 04 may 2005 00:00 ut by the 171   channel of aia is used as a working example throughout this section . the data is first pre - processed using the standard sdo solarsoft software ` aia_prep ' . this observation is shown without further image processing in the left panel of figure [ fig1 ] . the 171   channel of aia has a narrow bandwidth which is dominated by emission from highly - ionised iron which has a formation temperature which peaks at @xmath00.8 mk . the intensity of this image is therefore a function of the emitting plasma temperature and density . the line - of - sight and optical thickness of the plasma also has an effect on the measured intensity , as well as small contributions from other weaker spectral lines which share the same bandpass . the left panel of figure [ fig1 ] illustrates the main challenge of revealing information in euv images . for various physical reasons ( mostly due to density ) , the display range is dominated by the large difference between dark regions and the small bright regions near , or at , the base of active regions . most of the structural detail appears dark and is hidden . a quick and easy way of revealing some of this structure is by taking the square root of the image ( processing of this type is generally known as a gamma transformation ) . this is shown in the right panel of figure [ fig1 ] . although a large improvement on the unprocessed image , there is still much hidden structure , and bright areas become ` washed out ' . in particular , there is very little visible structure off the limb . the same is true of applying a logarithmic transformation , or any one - to - one mapping between input and output . the dominance of the image contrast range by small bright regions can be shown quantitatively by selecting regions and plotting histograms of pixel values . this is shown in figure [ fig2 ] for four areas - quiet sun , active region base , active region , and off - limb . note that these are pixel values normalized by the exposure time . the top @xmath075% of brightness values is taken by the small very brightest region of the active region ( cyan colour ) . this leaves the information from all other regions in the lower 25% . in fact , the off - limb and quiet sun regions have pixel values restricted to below @xmath0300 . furthermore , the bright active regions contain a large number of pixels which share the low values of the quiet sun . there are sharp boundaries between very dark and very bright regions . these characteristics become more extreme when flares occur . any image processing which attempts to reveal the information hidden in euv images must deal with the large range of brightness values , preserve the large spatial gradients in brightness and respect the fact that even the brightest regions contain dark patches . another challenge is the existence of structures at many spatial scales . we can not apply edge - detection or enhancement exclusively at small scales since this will enhance some features whilst other structures and the larger - scale context is lost . furthermore , with full resolution aia images of @xmath1 pixels taken every 11s ( and the anticipation of the resolution of future instruments ) , computational efficiency is very important . barring spurious calibration errors , the brightness values in the original euv image will be positive . let @xmath2 be the original image normalised by exposure time . @xmath3 is a 2d gaussian kernel of width @xmath4 pixels in the @xmath5 and @xmath6 image dimensions . throughout this work , @xmath4 denotes the one - sigma width of the gaussian . a normalized image @xmath7 can be computed by @xmath8    where    @xmath9\otimes k_w}.\ ] ]    for a given pixel , the numerator of equation ( [ eqn1 ] ) effectively subtracts the local ` mean ' ( weighted by the gaussian function centered on the pixel ) , and the denominator @xmath10 is the local ` standard deviation ' , also weighted by the gaussian kernel . thus , simply speaking , @xmath7 is locally normalized to a standard deviation of one and a mean of zero . an example of @xmath7 , computed with @xmath11 is shown in the left image of figure [ fig3 ] . this process of local normalisation is in essence similar to adaptive histogram equalization , where each pixel is scaled according to the statistics of the local group of pixels . the main difference is that histogram equalization enables precise control of the output pixel value range , whereas gaussian normalisation does not . this problem is solved by means of a further transformation applied to image @xmath7 :    @xmath12    since c is distributed across both negative and positive values , the @xmath13 function in formula ( 3 ) serves a purpose similar to a gamma transformation i.e. pixel values near the zero are amplified and extreme values far from zero are attenuated . this transform ensures the control over the output pixel value range and prevents saturation of output pixel values . an example of @xmath14 is shown in the right image of figure [ fig3 ] . the parameter @xmath15 controls the severity of the transformation . a default value of 0.7 gives good results for all the data shown in the results section . the pixel values of image @xmath7 have a standard deviation @xmath01 , and a value for @xmath15 of @xmath00.7 gives sensible results . a set of @xmath16 locally - normalized and @xmath13-transformed images , @xmath17 , is created for @xmath16 values of @xmath18 . values used throughout this work are @xmath19 . a global gamma - transformed image , @xmath20 is also created by    @xmath21    @xmath22 can be set at values of 2.5 to 4 . throughout this work , we use @xmath23 , and this gives good results for all instruments and bandpasses we have processed . @xmath24 and @xmath25 set the minimum and maximum input values respectively . for processing single images , they can be set to the minimum and maximum image values . for batch processing of a large number of observation , they can be set as the minimum and maximum values of the whole set of observations . for aia / sdo observations , these values can be easily found at the start of a batch process by reading the appropriate values stored in the data headers ( that is , without having to read in the full images ) . the final processed image , @xmath26 is calculated by a weighted average of the @xmath17 using weights @xmath27 , and addition of @xmath20 weighted by a global weight @xmath28 :    @xmath29    the @xmath27 can be calculated by a study of images consisting of only normally - distributed random noise . for narrow gaussian kernels , the mean of the local standard deviations across the whole image @xmath30 is an underestimate of the global standard deviation , and there is a larger variation in the local standard deviation . this is intuitive - the ratio of the mean of the local standard deviation to the global standard deviation approaches unity with a wider kernel , and there is less variation in the local standard deviation . this is shown in figure [ fig4 ] . therefore the weights @xmath27 are set at smaller values for small @xmath4 , and approach values near one as @xmath31 , as listed in figure [ fig4 ] . as described in equation ( [ eqn3 ] ) , the original ( normalized ) image @xmath20 is included to give contextual information of the largest scale structure . we use @xmath32 in this work . in practice , the weights @xmath27 and @xmath28 may be adjusted according to the desired output , and also according to the type of input image ( e.g. wavelength or channel ) . for most purposes , the @xmath27 can be set equal for all scales so that a straightforward mean of the locally - normalized @xmath17 is made .      1 . replace spurious negative pixels with zero or local mean / median . 2 .   create gaussian kernel of width @xmath18 . kernel elements should sum to unity . 3 .   convolve image with kernel to create local mean image @xmath33 . 4 .   calculate difference between image and the local mean image , square the difference , and convolve with kernel . square root the resulting image to give ` local standard deviation ' image @xmath10 ( equation ( [ eqn2 ] ) ) . calculate normalized image @xmath34 by subtracting the local mean image and dividing by the local standard deviation image ( equation ( [ eqn1 ] ) ) . store the result . 6 . apply @xmath13 transformation on @xmath34 to give @xmath17 . repeat 2 - 6 with the different kernel widths @xmath18 . take mean , or weighted mean if preferred , of the @xmath17 to give a weighted mean locally - normalised image . calculate a global gamma - transformed image @xmath20 by applying equation ( [ eqn2b ] ) . sum the weighted mean locally - normalised image with the global normalized image @xmath20 , with appropriate weight @xmath28 ( as equation ( [ eqn3 ] ) ) .    note that considerable efficiency is gained by creating a one - dimensional gaussian kernel for convolving first along the @xmath5-direction , then convolving the resulting image along the @xmath6-direction with a transposed kernel . this is far more efficient than convolving directly with a two - dimensional kernel . compared to the nafe @xcite , or wavelet - based routines , the mgn method is a simple one , and as gaussian smoothing is a standard practice employed by most programming languages , can readily be programmed in a few lines of code . without great rigour , we test the efficiency of the processing using the interactive data language ( idl ) on a macbook pro laptop with 8 gb ram and a 2.6ghz intel core i7 processor . the computing time increases approximately linearly from @xmath01s for a 500@xmath35500 image to 10s for a 2k@xmath352k image . a full aia / sdo 4096@xmath354096 image takes 40s . this is extremely fast compared to other methods such as the nafe ( at least an order of magnitude faster ) , and does not require a high - performance computer . figure [ fig5 ] shows the result of applying mgn to the example aia image of the previous sections . structure is enhanced down to fine spatial scales , although large - scale context is preserved ( that is , bright regions are still brighter than dark regions ) . off - limb structure , where present , is enhanced without being swamped by noise . structure is enhanced in the extremely bright region at the base of an active region ( the cyan - bordered region in figure [ fig2 ] ) , which originally appeared saturated . structure can also be seen in dark regions immediately neighbouring the bright region . it is also very useful to be able to trace structure from origins on the disk to off - limb heights . in this example , there are loops and fan - like bright features which extend from the disk and past the limb . movies created using the processing reveal dynamic features at very small scales such as flows within filament channels , flows along large coronal loops , small - scale activity within active regions , and movement in extended off - limb structures in response to dynamic events . an example aia movie is available online ( aia_movie.gif - an animated gif which can be viewed with a web browser ) which shows a large portion of the west disk and off - limb region observed in the 171   channel during 18 january 2013 . activity and small - scale flows are clearly seen , alongside larger - scale eruptions , and the movement / expansion of loops off the limb out to the edge of the field of view .    another online movie ( comparison_movie.gif ) compares the appearance of a sequence of aia observations processed by the mgn ( left panel ) , and by a simple gamma transformation ( right panel ) . this sequence shows a flaring active region and the eruption of a cme observed in the 171   channel during 08 march 2011 in the south - west . the disk and off - limb structure is seen far more clearly in the mgn images . individual loops are resolved , even in the complicated off - limb region to the north of the active region . in the original data , the sequence is dominated by the bright core of the active region , and when the flare peaks , by the small saturated regions centered on the flare . the mgn images show this activity clearly , and also succeed in maintaining a clear and stable view of surrounding structure throughout the movie , including the ` wobble ' of the active region loops in response to flaring activity . this enables a more comprehensive analysis of the event , showing connections between activity and movement of structure which can be hidden without mgn processing . figure [ fig6 ] shows a hi - c image processed with mgn . a movie showing the northernmost portion of this image is available online ( hic_movie.gif ) . various features such as magnetic braiding @xcite are clearly seen , as are plasma streams along the dark filament channels , and small - scale movement within the magnetic ` moss ' . such features can not be seen without appropriate processing . figure [ fig7 ] shows the application of mgn to an observation by the _ sun watcher using aps _ ( swap ) instrument aboard the _ project for onboard autonomy 2 _ ( proba2 ) satellite @xcite . this instrument has a coarser resolution than aia , but has an extended field of view . figure [ fig7 ] reveals an erupting filament out to the extremity of the field of view , and other quiescent structures to @xmath01.5@xmath36 . again , even low - signal structures are enhanced without too much amplification of noise . the mgn processing is not limited to euv observations . figure [ fig8 ] shows its application to a white light coronagraph image of the extended inner corona by the _ large - angle and spectrometric coronagraph _ ( lasco ) aboard soho . the left image shows an observation of 14 january 2011 , with a long - term background subtracted to remove instrumental stray light and f - corona . a point filter has also been applied to reduce spurious bright pixels . the right image shows a mgn - processed image , with parameters @xmath37 and @xmath38 . smaller - scale features are enhanced , in particular faint plumes over the poles which are difficult to see in the original image . in this respect , the mgn is an improvement over the nrgf @xcite , although the nrgf provides structural context closer to the true k - corona . it is difficult to make a qualitative comparison of images created using different processing . different processes may serve different analysis purposes - the important function for the process described here is to reveal information hidden in the original images . a non - rigorous comparison ` by eye ' of the same images processed using the mgn and of the nafe @xcite suggests that the results are similar , with the nafe giving slightly better clarity overall , and improved noise suppression - most obvious at large heights off the limb . for the mgn , noise reduction arises naturally from taking the ( weighted ) mean across many spatial scales . the noise suppression is not as effective as the nafe , and the mgn method lacks control over the degree of noise suppression . the nafe method is also better suited to reveal contrast in very bright regions . the mgn method normalizes an image by using the local mean and standard deviation calculated using a gaussian - weighted sample of local pixels . this normalised image is transformed by the @xmath13 function ( similar to a gamma transformation ) . this is applied over several spatial scales , and the final image is a weighted combination of the normalized components . the results compare well with multiresolution wavelet enhancement or the nafe procedure , but is far more computationally efficient . we hope that the mgn will become an established tool for researchers , offering a good compromise between computational time and clarity of the final images . the method is simple to implement , and the author is happy to provide the idl code by email request . we are grateful for comments by the anonymous referee which improved this work . huw morgan is grateful for funding from the coleg cymraeg cenedlaethol to prifysgol aberystwyth and support of shine grant 0962716 and nasa grant nnx08aj07 g to the institute for astronomy , university of hawaii . the work of miloslav druckmller was supported by grant agency of brno university of technology , project fsi - s-11 - 3 . we acknowledge the high resolution coronal imager instrument team for making the flight data publicly available . msfc / nasa led the mission and partners include the smithsonian astrophysical observatory in cambridge , mass . ; lockheed martin s solar astrophysical laboratory in palo alto , calif . ; the university of central lancashire in lancashire , england ; and the lebedev physical institute of the russian academy of sciences in moscow .                    , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , : 2011 , the atmospheric imaging assembly ( aia ) on the solar dynamics observatory ( sdo ) . _ _ solar phys . _ _ , 172 .
Abstract:

Answer: extreme ultraviolet images of the corona contain information over a large range of spatial scales , and different structures such as active regions , quiet sun and filament channels contain information at very different brightness regimes . processing of these images is important to reveal information , often hidden within the data , without introducing artifacts or bias . it is also important that any process be computationally efficient , particularly given the fine spatial and temporal resolution of _ atmospheric imaging assembly _ on the _ solar dynamics observatory _ ( aia / sdo ) , and consideration of future higher - resolution observations . a very efficient process is described here which is based on localized normalizing of the data at many different spatial scales . the method reveals information at the finest scales , whilst maintaining enough of the larger - scale information to provide context . it also intrinsically flattens noisy regions and can reveal structure in off - limb regions out to the edge of the field of view . the method is also successfully applied to a white light coronagraph observation .


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: an understanding of the origin , properties and evolution of our galaxy requires a careful census of its constituents , in particular its stellar members . of special importance are the intrinsic physical properties of these stars . the fundamental properties are mass , age and abundances , as these determine a star s history and future development . however , ages are not observable , and masses can only be directly obtained from some multiple systems . thus we must indirectly gain this information via the stellar spectrum , and a number of atmospheric parameters have been defined for this purpose . the main ones are the effective temperature , , the surface gravity , , and the metallicity , . to these can also be added the alpha abundances ,  ( which measure the devations away from the ` standard ' abundance ratios ) , the photospheric microturbulence velocity , , and the extinction by the interstellar medium , ( although not intrinsic to the star , it is necessary for determining its luminosity ) . masses and ages can then be determined from stellar structure and evolution models and with calibration via binary systems . it is important to realise that this modelling is complex , and a number of assumptions have to be made . there is , therefore , a limit to the precision with which we can determine physical properties . historically , spectroscopic parameters have been measured indirectly through the mk classification system ( morgan et al . @xcite ) or via colour - magnitude and colour - colour diagrams . in the mk system , the two parameters _ spectral type _ and _ luminosity class _ act as proxies for  and . originally a qualitative system relying on a visual match between observed spectra and a system of standards , much progress has been made in quantifying it with automated techniques ( e.g.  weaver & torres - dodgen @xcite ; bailer - jones et al.@xcite ) . the most commonly used classification techniques have been neural networks and @xmath4 matching to templates ( or more generally , minimum distance methods ) . a summary of recent progress in this area is given by von hippel & bailer - jones ( @xcite ) .    despite this focus on the mk system , it is not well suited to classifying data from the deep surveys which will be central to the future development of galactic astrophysics . this is for a number of reasons , but in particular because it lacks a measure of metallicity . although mk does make allowance for various ` peculiar ' stars , these are defined as exceptions , and the notation is not suited to a statistical , quantifiable analysis . this is problematic given the significance of metal poor halo stars in a deep survey . there is also now no good reason why we should not determine physical parameters directly from the observational data . some attempts have been made to determine the physical parameters of real spectra directly by training neural networks on synthetic spectra . gulati et al . ( @xcite ) used this approach to determine the effective temperatures of ten solar metallicity g and k dwarfs . taking the `` true '' effective temperature of these stars as those given by gray & corbally @xcite , they found a mean `` error '' in the network - assigned temperatures of 125k . bailer - jones et al . ( @xcite ) determined  for over 5000 dwarfs and giants in the range b5k5 , and also showed evidence of sensitivity of the parametrization models to metallicity . the accuracy with which physical parameters can be determined from a stellar spectrum depends upon , amongst other things , the wavelength coverage , spectral resolution and signal - to - noise ratio ( snr ) . from the point of view of designing a stellar survey project it is essential to know how well the stellar parameters can be determined for a given set of these observational parameters . moreover , given that there is always a limit to the collecting area and integration time available , there is always a trade - off between spectral resolution , sensitivity and sky coverage . the goal of this paper is to determine the accuracy with which physical stellar parameters can be determined from spectroscopic data at a range of snrs and resolutions which could realistically be achieved in a _ deep _ survey mission . this specification rules out high resolution spectra . the parametrization work has been carried out using neural networks ( sect . [ networks ] ) because they have been shown to be one of the best approaches for this kind of work . this is not to presuppose , however , that some other approach may not ultimately be better . the simulations have been made using a large database of synthetic spectra generated from kurucz atmospheric models ( sect . [ synspec ] ) . while these spectra do not show the full range of variation in real stellar spectra , they are adequate for a realistic demonstration of what is possible as a function of snr and resolution . the results are presented in sect . [ results ] and summarised and discussed in sect .  [ discussion ] . finally , the requirements for a complete survey - oriented classification system are given in sect . [ system ] . the simulations in this paper were partially inspired by the need to produce an optimal photometric / spectroscopic system for the gaia galactic survey mission . gaia is a candidate for the esa cornerstone 5 mission for launch in 2009 ( esa , in preparation ) . it is primarily an astrometric mission with a precision of a few microarcseconds , and will survey the entire sky down to v=20 , thus observing c.  10@xmath5 stars in our galaxy . radial velocities will be obtained on board down to v=17.5 , thus providing a 6d phase space survey ( three spatial and three velocity co - ordinates ) for stars brighter than this limit . a survey of this size will have a profound impact on galactic astrophysics , but to achieve this it is essential that the physical characteristics of the target objects are measured and correlated with their spatial and kinematic properties . as gaia is a continuously scanning satellite , a fixed total amount of integration time is available for each object , so there is a trade - off between resolution , signal - to - noise ratio and wavelength coverage . for various reasons , the current gaia design does not include a spectrograph ( other than a 1.5  resolution region between 8470 and 8700  intended for radial velocity measurements ) , but instead will image all objects in several medium and broad band filters ( table  [ filttab ] ) . three filter systems are shown : the system nominally selected for the mission plus two alternatives . the profiles of the two alternatives are represented as gaussians in this paper . the ability of these filter system to determine stellar parameters will be compared with that for spectra of various resolutions . .three multiband filter systems proposed for the gaia mission . all profiles are symmetric about the central wavelength , @xmath6 , and have a fwhm of @xmath7 . the profiles of the filters in the asiago and modified strmvil systems ( f.  favata 1999 , private communication ) are defined as gaussians ( although note that the former is only an approximation to the original asiago system in munari @xcite ) . the filters of the selected gaia system ( esa , in preparation ) have flatter tops and steeper sides than gaussians , and have defined relative peak transmissions , t. there is some ( intended ) redundancy within each filter system . [ cols="^,>,^,>,^,>,^ " , ]     the results of the parametrization process are shown in figs  [ results1][results3 ] and tabulated in tables  [ restab1][restab3 ] . before interpreting these results we should consider the limits which the data themselves place on the performance . first , the network will be unable to produce errors smaller than the smallest variations in the data set . if , to take a hypothetical example , the spectra did not change as the metallicity changed by 1.0 dex , we could not expect the network to determine  to much better than 0.5 dex . second , the grid of atmospheric models represents the physical parameters at a finite sampling , e.g.  a constant step size of 0.5 dex for . this sampling does not in itself limit the precision achievable ; it is perfectly possible for the network to legitimately give an error much smaller than the sampling because the network is minimising a continuous error function and not just obtaining the best match between a spectrum and a set of templates . nonetheless , the network input  output mapping is an _ interpolation _ of the training data , and the more coarsely sampled the parameter grid the harder it is for the network to get a reliable interpolation . consequently , while the network _ may _ be able to achieve sub - sampling accuracy , we should not be surprised if it can not . thus to avoid over - interpreting these results we should not compare two errors which are both smaller than half the sampling level . the _ average _ ` half - sampling ' values for  and  are 0.2 and 0.25 respectively , and for log  in the three temperature ranges ( cool , intermediate and hot ) are 0.01 , 0.01 and 0.03 respectively . the implication is that , if the network produces errors smaller than these half - sampling values ( as it does ) , we can not know whether the performance is limited by the network model or by the data themselves . a distinction will only be possible with a more sensitive and finely sampled grid of atmospheric models .    with the above caveat taken into account , i draw attention to some interesting features in figs  [ results1][results3 ] .    1 . good  determination is possible with all resolutions / filter systems and snrs . the larger error in  above 10000k may be an artifact of the larger half - sampling value in this region ( @xmath81000k ) . 2 .   only at high resolution can  be determined for the coolest stars and even then the determination is poor relative to the hotter stars . this is understandable , at least in part , because the spectral signature is primarily in the line widths which are only apparent at high resolution . although the three filter systems differ somewhat , they give essentially the same performance as each other . the filter systems ( each with 1015 input parameters ) have similar  and  as the @xmath9=400  spectra ( 35 inputs ) . 5 .   at low snr , the @xmath9=400  spectra and the filters give poor  and very poor  determination for all three temperature ranges . 6 .   at high snr ( 1000 ) all resolutions / filter systems appear to be equally good at determining any of the parameters . differences will probably become apparent with a more sensitive training grid . 7 .   at higher temperatures the accuracy is more sensitive to snr than at lower temperatures . 8 . metallicity determination is more difficult at higher temperatures , especially for the filters and low resolution spectra . this is understandable as at high temperature there are fewer and weaker metal lines which are only significant at high snr and/or resolution . 9 .   in most cases there is little difference between the performances of the @xmath9=25 , 50 and 100  spectra , at least for this data grid . the results demonstrate that a fully automated neural network can accurately determine the three principal physical parameters from spectroscopic or photometric stellar data , something which has not previously been demonstrated . moreover , this work has used spectra of considerably lower resolution than have been used before in automated classifiers . even at low resolution ( 50100  fwhm ) and snr ( 510 per resolution element ) , neural networks can yield good determinations of  and , and even for  for stars earlier than solar . still lower resolutions permit good results provided the snr is high enough ( @xmath3 ) . that good  can be achieved even at low resolution and snr is perhaps not surprising when we consider that the spectra have absolute fluxes , which will be the case with high precision parallax missions such as gaia . however , the more distant objects will have lower precision parallaxes and hence errors in the _ mean _ flux level . but even if we completely ignore distance information ( and flux normalise the spectra ) , the shape of the spectrum is still a strong indicator of : for example , bailer - jones et al .  ( @xcite ) obtained an mk spectral type precision of 0.8 subtypes ( @xmath10=0.0100.015 ) across a wide range of spectral types ( b2m7 ) using flux normalised spectra . this is similar to what can be achieved from broad band photometry , implying that  determination only requires very low resolution . the good performance of ` high ' resolution spectroscopy ( 25 ) at very low snr ( @xmath11 per pixel ) was not expected . it seems to imply that for a given amount of integration time it may be better to sacrifice snr for resolution . it is noteworthy that while the filters provide good , their ability to determine  and especially  is very limited at low snr . how do these results compare with classical parametrization methods ? gray ( @xcite ) compiles results showing that with photometric errors below 0.01 magnitudes , the b@xmath0v colour calibrates  to 23% ( 4% for hotter stars ) in the absence of reddening . slightly better precision can be obtained from the slope of the paschen continuum and size of the balmer discontinuity . the latter may also be used to measure  to @xmath2 dex . with spectra at a few   resolution over a similar wavelength range to that used here , cacciari et al . ( @xcite ) obtained uncertainties in log  and  of 0.01 and 0.04 respectively . sinnerstad ( @xcite ) made uvby,@xmath12 photometric measurements of b stars , and for uncertainties of 0.005 in @xmath12 and of 0.01 in u@xmath0b ( i.e.  snr @xmath13 200 ) , infers errors in log  and of 0.004 and 0.08 respectively . these are similar to or slightly better than the results for similar stars in tables  [ restab1][restab3 ] ( @xmath14 ) at the highest resolutions . high resolution ( @xmath15 )  spectra have generally been used to determine metallicity , and in a review , cayrel de strobel ( @xcite ) notes that metallicity can be determined to @xmath16 dex at snr=250 ( but only @xmath2 dex at snr=50 ) provided the effective temperature and gravity are approximately known . at lower snr ( 1020 ) , jones et al.(@xcite ) could determine [ fe / h ] to @xmath2 dex for g stars using a set of spectroscopic indices measured at 1  resolution in the range 40005000 , again using a known effective temperature . more recently , katz et al . ( @xcite ) have used a minimum distance method to parametrize spectra by finding the closest matching template spectrum . the template grid consisted of 211 flux calibrated spectra ( 39006800 , @xmath17 ) with 4000k@xmath18@xmath186300k , @xmath19[fe / h]@xmath20 , and  for dwarfs and giants . the _ internal accuracy _ of the method for log  ,  and  was 0.008 , 0.28 dex , and 0.16 dex respectively at snr=100 , and 0.009 , 0.29 dex and 0.17 dex at snr=10 . as expected , their results for  are much better than those in this paper at the similar temperature range ( @xmath21 in table  [ restab2 ] ) , presumably due to their much higher resolution . in contrast , their performance for  is similar and for  somewhat worse than that in this paper at 500 times lower resolution . their results also confirm that at high resolution a lower snr leads to very little loss in performance . snider et al.(@xcite ) trained and tested neural networks on a set of 182 real f , g and k spectra over the range 36304890  at intermediate resolution ( @xmath131 ) , and achieved 1@xmath22 errors in log  ,  and  of 0.015 , 0.41 dex and 0.22 dex respectively , based on training and testing a network with a set of 182 real f , g and k spectra .    when judging the relative values of the different resolution / snr combinations in this paper , we must also take account of their implementation ` costs ' , specifically the relative integration times required . usually for a survey , a fixed total amount of integration time is available for all filters / spectra . in the case of gaia  which is continuously rotating  a star moves across a focal plane covered with a mosaic of ccds which are clocked at the rotation rate . the different filters are fixed to different ccds , so that as a star moves across the mosaic it is recorded in different wavelength ranges . thus less numerous and/or broader filters would achieve a higher snr than more or narrower filters . some filters could be replaced with a slitless spectrograph ( e.g.  a prism or grism ) . this disperses every point on the sky and thus gives the full integration time for all wavelengths , but at the expense of increased sky noise and object confusion . these could be reduced by using one or more dichroics to redirect the light to two or more focal planes . ( confusion would be reduced further with gaia by the fact that each area of sky is observed at many different position angles over the mission life . ) an alternative approach is a set of many medium band filters ( @xmath23 for @xmath9=100  over the complete wavelength range , although omission of some filters could be achieved ) . while this avoids the two principal disadvantages of the slitless spectrograph , the integration time per wavelength interval is dramatically reduced . the development of a complete parametrization system will require further research , much of which needs to be directed at taking better account of the true nature of the observational data . directions and suggestions for the course of this work are now given . essentially all of the work in the literature on automated classification deals with preselected objects . in contrast , an unpointed survey will pick up a whole range of objects , necessitating a filtering system to select the stars . such a system could make use of both object morphology and spectral features , and systems based on neural networks ( e.g.  odewahn et al . @xcite ; miller & coe @xcite ; serra - ricart et al . @xcite ) and principal components analysis ( bailer - jones et al.@xcite ) have been proposed . such a system must be relatively robust and always allow for ` unknown ' objects which can be dealt with manually . it will be necessary to have a stellar database for training which takes better account of the larger range of variation present in the galactic stellar population . ideally , a large set of real spectra across a wide range of physical parameters should be obtained for this purpose . good atmospheric models and synthetic spectra are nonetheless still required for determining their physical parameters and thus for training the network . there are two possible approaches to training . the first is to train on synthetic spectra suitably preprocessed to be in the same form as the observed spectra ( e.g.  bailer - jones et al.@xcite ) . the alternative is to obtain a representative sample of real spectra with the survey system , calibrate them , and then use them to train a network . in theory the latter method gives a better sampling of the true cosmic variance in the spectra , but of course requires that a representative sample is selected from the survey data . this sample could be improved as the survey progressed . neural networks are fast to train and apply , so it is realistic to expect that even for a database of @xmath24 objects the network could be retrained and applied to the whole database in less than a day . more advanced model atmospheres are required for a number of reasons :    1 .   ,  and  do not uniquely describe a true spectrum . models sensitive to different abundance ratios and which include chromospheres ( for example ) are necessary . kurucz models assume lte which is known to break down in a number of regimes ( e.g.  for very hot stars ) . 3 .   both the atmospheric models and the line lists lack water opacity , known to be important for cool stars , thus setting the current lower  limit of about 4000k . yet more advanced models ( which include dust ) are required for very cool stars ( l and t dwarfs ) and brown dwarfs , of which many will be found by gaia .      of particular importance is insterstellar extinction ( reddening ) , especially in deep surveys . the extinction can , in theory , be determined by the network by training it on artificially reddened synthetic spectra and providing the network with a `` reddening '' output parameter ( or parameters ) . this has been demonstrated on limited data sets by weaver & torres - dodgen @xcite and gulati , gupta & singh @xcite , who determined e(b@xmath0v ) to within 0.05 and 0.08 magnitudes respectively . the latter made use of 6  resolution uv spectra ( 48505380 ) . the former used red spectra ( 58008900 ) at 15  resolution and found that the spectral type and luminosity class classifications did not degrade much as reddening was added . it is therefore to be expected that the parametrizations in this paper will be robust to reddening , particularly as the spectra have a much larger wavelength coverage . the filter systems proposed for gaia were of course designed with interstellar extinction in mind , and a study of its impact has been carried out ( esa , in preparation ) . this work shows that suitable q parameters ( non - linear combinations of the filter fluxes ) used to determine the physical parameters are largely insensitive to reddening . it also claims that narrow band filters are not necessary for overcoming reddening . in some parts of the parameter space , reddening is more problematic ( e.g.  for k stars ) , largely due to a degeneracy between it and  and . however , at intermediate and high galactic latitudes it is expected that e(b@xmath0v ) can be determined to within 0.002 magnitudes . munari ( @xcite ) similarly shows that reddening - free indices exist for the asiago filter system . as a neural network also forms non - linear combinations of the filter fluxes , it is reasonable to suppose that it too will be robust to redenning , although this will be the subject of future work . the parametrization model used in a real survey must confront the fact that most stars are in spatially unresolved multiple systems . independent measurement of the physical properties of each component is desirable and in principle achievable  when the brightness ratio is large enough  by training the network with composite spectra . in this case the network model would need to have multiple sets of outputs to deal with each component . an alternative approach is to use ` probabilistic outputs ' in which the single output for , say , , is replaced with a series of outputs in which each value of  ( 6000 , 6250 , 6500 etc . ) is represented separately . the network then evaluates the probability that each temperature is present in the input spectrum . this method is not recommended , however , as it eliminates the intrinsically continuous nature of the physical parameters . it would also greatly increase the number of outputs and hence the number of free parameters ( weights ) in the network . object confusion should not result in any overlapped spectrum being rejected entirely . rather , it would be better to have a parametrization model which is robust to missing data . this is a major challenge for the feedforward network models used in this and most other papers on automated classification , and will presumably require some transformation of the input spectrum . an analysis of the effect of wavelength coverage on the parameter determination accuracy is important because a smaller spectral coverage ( or _ coverages _  it need not be contiguous ) would also reduce this confusion . finally , the model should make use of all available data . in the case of gaia , this means including the data from the high resolution spectrograph ( 84708700  at 0.75 )  used to measure radial velocities . as the inputs to the network need not be homogenous , there should be no problem incorporating different types of data . i would like to thank fabio favata , gerry gilmore and michael perryman for useful discussions on this work , in particular within the context of the gaia mission . i am also grateful to robert kurucz for use of his model atmospheres and compiled line lists , and to richard gray for the use of his synthetic spectra generation program .
Abstract:

Answer: large scale , deep survey missions such as gaia will collect enormous amounts of data on a significant fraction of the stellar content of our galaxy . these missions will require a careful optimisation of their observational systems in order to maximise their scientific return , and will require reliable and automated techniques for parametrizing the very large number of stars detected . to address these two problems , i investigate the precision to which the three principal stellar parameters ( , , )  can be determined as a function of spectral resolution and signal - to - noise ( snr ) ratio , using a large grid of synthetic spectra . the parametrization technique is a neural network , which is shown to provide an accurate three - dimensional physical parametrization of stellar spectra across a wide range of parameters . it is found that even at low resolution ( 50100  fwhm ) and snr ( 510 per resolution element ) , and  can be determined to 1% and 0.2 dex respectively across a large range of temperatures ( 400030000k ) and metallicities ( @xmath03.0 to @xmath11.0 dex ) , and that  is measurable to @xmath2 dex for stars earlier than solar . the accuracy of the results is probably limited by the finite parameter sampling of the data grid . the ability of medium band filter systems ( with 1015 filters ) for determining stellar parameters is also investigated . although easier to implement in a unpointed survey , it is found that they are only competitive at higher snrs ( @xmath3 ) .


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: in classical phenomenology , turbulence arises after a sequence of symmetry breakings , successive instabilities or bifurcations , which however progressively restore the system symmetries in a statistical sense . @xcite the study of these instabilities traditionally proceeds from ( linear or non - linear ) perturbations of the so - called  basic state " , the stationary laminar solution of the navier - stokes equation at low reynolds number . at finite reynolds number there is no general well - defined criterion to discriminate between the turbulent or laminar nature of the flow : one can not define clearly a critical threshold beyond which the flow is turbulent and below which it is laminar . the consensus view is that the flow is turbulent when the reynolds number is large enough and when a well established spatio - temporal energy spectrum is observed with broad - band power laws . in such a case a possible statistical equivalent of the laminar `` basic state '' can be defined using the ( statistically or time ) averaged flow . this flow is stationary by construction , but differs from an usual basic state in the sense that it is solution of the ensemble averaged navier - stokes equation instead of a solution of the plain navier - stokes equation . a natural question then arises : what happens once the statistically stationary turbulent state is reached ? is this the end of the story or can new instabilities of the averaged turbulent flow develop ? experimentally , the answer seems to be positive , since spontaneous bifurcations and flow reversals in fully developed turbulence have already been observed : i ) in a wake flow , there is a mean pattern transition at a critical reynolds number . this corresponds to the so - called `` drag - crisis '' ( see ref .  ) . this leads to a dramatic decrease of the mean drag of a sphere or a cylinder in a turbulent flow for a critical value of the reynolds number @xmath7 . the wake becomes narrower as the mean flow pattern changes . ii ) spontaneous flow reversals have been observed in thermally driven ( rayleigh - bnard ) convection , both experimentally@xcite and numerically.@xcite in a cylindrical rayleigh - bnard geometry the reversals of large scale circulation was reported to be induced by reorientation along the azimuthal direction.@xcite however it is difficult to experimentally characterize a bifurcation from a non - stationary or periodic flow . iii ) in a magnetohydrodynamic turbulence , the large scale magnetic field spontaneously generated in a liquid metal at @xmath8 was shown@xcite to undergo spontaneous reversals , with a dynamics governed by a few magnetic modes despite the strongly turbulent background . this observation led to a simple model of geodynamos , with a few modes undergoing chaotic dynamics.@xcite an open problem common to all these situations is to find suitable tools to study and characterize these type of instabilities . the purpose of this paper is to present some answers to these questions in the specific case of a von krmn flow . this flow is generated by two counter - rotating impellers in a cylindrical vessel . at low reynolds number , the laminar flow is axi - symmetric and divided into two toric recirculation cells separated by an azimuthal shear layer . its transition to turbulence for various counter - rotations and aspect ratio has been extensively studied theoretically , numerically and experimentally , with smooth or rough ( fitted with blades ) impellers ( see e.g. refs .   and references therein ) . the study in ref .   was performed in the exact - counter - rotating case . it was reported that a first bifurcation from the stationary and axi - symmetric laminar state occurs at @xmath9 , yielding a stationary flow with an azimuthal modulation , _ i.e. _ for which the axisymmetry is broken . this stationary state persists up to typically @xmath10 where time - dependence arises . the transition to turbulence further proceeds through plain , modulated or chaotic traveling waves , until @xmath11 , where a  fully developed turbulent state " seems to have been reached . in that state , the energy spectrum is broad and the dimensionless dissipation does not depend on reynolds number anymore . @xcite in that turbulent state , the flow once time averaged regains the structure of the basic laminar flow , made of two shearing toric recirculation cells . in the present paper , we have worked using the same experimental setup . we demonstrate that the flow is organized into large scale coherent structures . these structures are subject to a sequence of eckhaus type instabilities leading to an increasing azimuthal wavenumber as the reynolds number varies from @xmath1 to @xmath2 , similar to the sequence of instabilities at much lower @xmath0 . since this sequence of instabilities occurs on coherent structures in a turbulent flow rather than a laminar  basic state " , we are facing a new paradigm that both requires new tools of investigations and opens interesting questions about turbulence parametrization . the tools are presented in section  [ sec : technical ] , the instability is described in section  [ sec : results ] and some theoretical consequences are discussed in section  [ sec : theocons ] . we have worked with a von krmn flow generated by two coaxial and counter - rotating impellers in a cylindrical vessel . more details about the experimental setup can be found in ref .  . the cylinder radius and height are @xmath12 mm and @xmath13 mm respectively . the impellers consist of @xmath14 mm diameter disks fitted with sixteen curved blades of @xmath15 mm height . the impellers are driven by two independent motors . the motors frequencies are respectively set to @xmath16 and @xmath17 . in the present work the impellers are rotating with the convex face of the blades going forwards , contrarily to ref .  . the working fluid is either pure water , a water - glycerol mixture ( 26%-74% in weight respectively ) or pure glycerol . the resulting accessible reynolds numbers , @xmath18 with @xmath19 the kinematic viscosity , vary from @xmath1 to @xmath2 .    , measured in a meridian plane using s - piv , at @xmath20 and @xmath21 . @xmath22 , @xmath23 $ ] . right , corresponding instantaneous fluctuating component of the velocity field @xmath24 . the color maps the azimuthal velocity @xmath25 ( normalized by @xmath26 ) whereas the arrows map the ( @xmath27,@xmath28 ) velocities.,title="fig:",scaledwidth=35.0% ] , measured in a meridian plane using s - piv , at @xmath20 and @xmath21 . @xmath22 , @xmath23 $ ] . right , corresponding instantaneous fluctuating component of the velocity field @xmath24 . the color maps the azimuthal velocity @xmath25 ( normalized by @xmath26 ) whereas the arrows map the ( @xmath27,@xmath28 ) velocities.,title="fig:",scaledwidth=35.0% ]    this setup is invariant under @xmath29 rotations around any radial axis passing through the center of the cylinder . measurements are done thanks to a stereoscopic particle image velocimetry ( s - piv ) system . the s - piv provides time series of the 3 components of the velocity ( radial @xmath30 , vertical @xmath31 and azimuthal @xmath32 ) on a @xmath33 points grid in a meridian plane @xmath34 $ ] composed of two planes dephased by @xmath35 in the cylindrical coordinates . @xmath36 , with the azimuthal angle @xmath23 $ ] and @xmath37 the radial distance from the impellers axis . @xmath38 is the vertical distance from the center of the cylinder , @xmath39 is the out - of - plane direction and @xmath40 is the time . the azimuthal component is then written in the cylindrical coordinates using @xmath41 . the size of the recorded images is @xmath42mm@xmath43 ( where @xmath44 mm and @xmath45 mm are the spatial resolutions of the s - piv ) . time series have sampling frequency in the range @xmath46 to 15hz and are composed of @xmath47 to 4200 samples . each experiment is started with the impellers and the fluid at rest . the impellers velocities are then suddenly increased to reach their target values . after a few tens of seconds , a statistically stationary state is reached and s - piv time series is finally acquired .      for low reynolds number @xmath48 , the flow is steady and laminar and only few fluctuations are present . when @xmath49 , the instantaneous flow is composed of two toric recirculation cells separated by an azimuthal shear layer located at @xmath50 reflecting the @xmath51-symmetry of the system . increasing the reynolds number , one expects to reach fully developed turbulence around @xmath52 as observed in ref .  , with increasing difference between the instantaneous flow and the time - averaged flow . above @xmath53 and at @xmath54 , the axisymmetry is broken by fluctuations of the instantaneous flow . however , this symmetry is restored for the time - averaged flow ( _ cf . _ fig.[fig : velocityfield_exp ] ) , at any reynolds number . by construction , this average flow is time - independent and will be the equivalent of a  basic state " . as the reynolds number is increased , this  basic state " is subject to increasingly stronger fluctuations . @xcite the question we want to address here is : are these fluctuations purely disorganized as in thermal noise or can we identify a pattern in them when varying the reynolds number,@xcite in a way similar to a sequence of bifurcations of ordinary instabilities ? to answer this question we need two steps : i ) to separate the fluctuations into  organized " motions and  thermal " noise ; ii ) to build a suitable tool to study the temporal behavior of these organized motions . these steps are described below . for step i ) , we use comparison with and projection onto beltrami modes . for step ii ) , we use spatio - temporal spectra , as used in wave turbulence.@xcite    [ cols="^,^,^,^ " , ]      finally the existence of organized motions in the turbulent flow means @xmath55 can be written as the sum of an in - phase term and a fluctuation term : @xmath56 , with @xmath57 the fluctuations and @xmath58 the coherent structure velocity field . a visualization of the spatial structure of one azimuthal mode @xmath59 can then be obtained using conditional average of the velocity field @xmath55 . the condition is given by suitable projection onto the corresponding beltrami mode , see eq.([projection ] ) . due to the geometry of our piv device , we however only have access to the velocity field in a meridional plane , i.e. for @xmath23 $ ] . since the azimuthal modes are rotating we may however replace the azimuthal integration in the scalar product by a suitable time integration . this is made possible for high enough reynolds number where the rotation frequency of the coherent structure is resolved . we used conditional average as follows : we first compute for each instantaneous experimental field @xmath55 the quantity : @xmath60 keeping the azimuthal wavenumber @xmath59 and the angle @xmath61 constant . a result is shown in fig.[fig : am](a ) obtained using the same experimental data shown in fig.[fig : azimuthalmodes](d ) at @xmath20 and the synthetic velocity field as in fig.[fig : velocityfield](d ) at @xmath62 . one sees quasi - sinusoidal fluctuations in the signal @xmath63 . we propose that these fluctuations correspond to the coherent structures passing through the piv plane . this is corroborated by looking at @xmath64 , the temporal power spectrum of @xmath63 in fig.[fig : am](b ) . the rotation frequency measured using the 2d - spectrum of the velocity field and the main peak frequency of @xmath65 ( @xmath66hz in this experiment ) are found identical . secondary peak is a sub - harmonic at @xmath67hz . for lower reynolds number , below @xmath68 , both amplitudes of these peaks decrease and are shifted toward lower frequencies . simultaneously the energy becomes roughly equally distributed along a plateau , which amplitude is increasing and width is decreasing with the reynolds number . down to @xmath69 a peak at small amplitude can still be distinguished at @xmath70hz . for the lower reynolds number , the plateau vanishes and the energy is mainly concentrated in the vicinity of @xmath71 . this confirms the behavior described in section[subsec : tempbehav ] .     with @xmath62 and @xmath20 at a fixed angle @xmath61 ( corresponding to the angle of fig.[fig : velocityfield](c ) ) as a function of time . @xmath63 was computed using @xmath72 from the time - series used in fig.[fig : azimuthalmodes](d ) and was finally low - pass filtered at 4hz . horizontal lines , see fig.[coherentstructure ] . ( b ) power spectrum @xmath64 of @xmath63 ( @xmath73 solid line ) for the data shown in ( a ) and for other reynolds numbers : @xmath74 ( thin red / gray line ) , 21 ( thin black line ) , 12 ( thick black line ) , 7 ( @xmath75 ) . the arrow marks the main peak evolution with decreasing re.,title="fig:",scaledwidth=34.5% ] with @xmath62 and @xmath20 at a fixed angle @xmath61 ( corresponding to the angle of fig.[fig : velocityfield](c ) ) as a function of time . @xmath63 was computed using @xmath72 from the time - series used in fig.[fig : azimuthalmodes](d ) and was finally low - pass filtered at 4hz . horizontal lines , see fig.[coherentstructure ] . ( b ) power spectrum @xmath64 of @xmath63 ( @xmath73 solid line ) for the data shown in ( a ) and for other reynolds numbers : @xmath74 ( thin red / gray line ) , 21 ( thin black line ) , 12 ( thick black line ) , 7 ( @xmath75 ) . the arrow marks the main peak evolution with decreasing re.,title="fig:",scaledwidth=34.0% ]    the velocity field @xmath55 was then averaged over times for which @xmath76 lies in a given range of values @xmath77 $ ] , representative of a given phase of the velocity field . an example is provided in fig.[coherentstructure ] for @xmath78 and @xmath62 from @xmath63 shown in fig.[fig : am ] . one sees the presumable @xmath62 coherent structure in two planes , dephased by @xmath35 . one indeed sees large scale structures looking like pairs of azimuthal vortices . by reproducing this conditional average using different phase @xmath79 for @xmath80 at a fixed @xmath59 , we obtained the complete movie ( see fig.[coherentstructure ] ) of the rotation of this pair of vortices ( enhanced online ) . the physical mechanism leading to this vortices can be the kelvin - helmoltz instability laying in the shear layer , a secondary or related instability or the centrifugal instability . it is interesting to compare this result with the @xmath62 velocity field shown in fig.[fig : velocityfield](d ) using the same phase @xmath79 as in fig.[coherentstructure](a ) . both the anti - symmetric nature of the velocity field with respect to the @xmath81 vertical axis and @xmath82 are recovered .    ) $ ] in the plane @xmath83 $ ] and ( right ) @xmath84)$ ] in the orthogonal plane @xmath85 $ ] . both are computed using the signal @xmath63 and the horizontal lines represented in fig.[fig : am](a ) for @xmath20 , @xmath62 and @xmath78 . we clearly distinguish two pairs of azimuthal vortices , rotating around the azimuthal axis . the observed discontinuity at @xmath86 is due to a slight shift of our s - piv laser sheet from the center of the tank . ( enhanced online),title="fig:",scaledwidth=35.0% ] ) $ ] in the plane @xmath83 $ ] and ( right ) @xmath84)$ ] in the orthogonal plane @xmath85 $ ] . both are computed using the signal @xmath63 and the horizontal lines represented in fig.[fig : am](a ) for @xmath20 , @xmath62 and @xmath78 . we clearly distinguish two pairs of azimuthal vortices , rotating around the azimuthal axis . the observed discontinuity at @xmath86 is due to a slight shift of our s - piv laser sheet from the center of the tank . ( enhanced online),title="fig:",scaledwidth=35.0% ] in 1970 , kraichnan wrote a paper entitled  instability in fully developed turbulence " ,  @xcite in which he focuses on the propagation of errors from one range of scale to another . this problem is related to predictability of turbulence . in the present work , we consider another aspect of the same problem related to the stability of the unsteady coherent structures of a turbulent flow ( in other words , to the dynamic of the perturbation of this given coherent structures ) . we have shown that using spatio - temporal spectra , and comparison or projection onto suitable beltrami modes , it is possible to identify coherent perturbations to the time - averaged flow . this is liable to an eckhaus instability analog to what is observed in systems in the laminar limit , with a few degrees of freedom . in classical phenomenology of turbulence the number of degrees of freedom scales like @xmath87 . given our range of reynolds number ( from @xmath1 to @xmath2 ) , we thus can expect our system to be described by a number of degrees of freedom in the range @xmath88 to @xmath89 . our results thus prove that a huge number of degrees of freedom are irrelevant to describe the instability of our turbulent flow . this sets a number of interesting questions , that we list below .    in the present paper , we have characterized experimentally the instability of the basic state . is there any analytical or numerical way to explain our observations ? classical instability analysis for a given velocity field @xmath90 starts from stationary solutions of the navier - stokes equation : @xmath91 in our case , we consider the instability of a time - average basic state @xmath92 , that obeys the equation : @xmath93 where @xmath94 is the reynolds stress tensor . this reynolds stress represents the influence of all the degrees of freedom of the flow onto its average , and can , in general , only be computed via full solution of the ns equation . therefore , the problem of instability of a mean turbulent flow can not be tackled analytically or is too demanding numerically , unless a prescription ( parametrization ) of the reynolds stress is provided . in the case of the plane couette turbulent flow , for example , this was attempted by tuckerman _ et al . _ @xcite via the k - @xmath95 closure model . they calculate steady 1d solution profiles of the k - @xmath95 model and their linear stability to 3d perturbations , but find no correspondence between this analysis and the onset of turbulent - laminar bands in experiments and simulations . in the same way , legras and villone@xcite address the problem of kolmogorov flow instability when molecular viscosity is replaced by a smagorinsky parametrization for small - scale turbulence . such a parametrization represents the motion at scales smaller than the large scale kolmogorov flow . they claim that it may provide hints on large - scale instabilities at large @xmath0 and , hopefully , on the character of such instabilities , but it has never been checked . in our case , it was observed that the sequence of eckhaus instability resembles the sequence of instability observed at much lower reynolds number . @xcite this suggest that the replacement of the reynolds stress by an eddy - viscosity might be a suitable parametrization by increasing the dissipation , and therefore decreasing the effective reynolds number ( and the effective number of degrees of freedom ) . theoretically , we note that the special geometry of our experiment and its parity symmetry properties make it plausible that the first term in the expansion of the reynolds stress as a function of the velocity gradient is indeed quadratic , resulting in a non - isotropic eddy - viscosity tensor ( in other words , there is probably no anisotropic kinetic alpha ( aka ) term in the reynolds stress expansion ) . @xcite it is however not clear that the correct parameters of the instability ( threshold , wavenumbers ) can be captured with a simple eddy - viscosity model , since this procedure fails in the case of a plane couette flow . @xcite an interesting alternative would be to derive directly the mean state as the critical points of a suitable arnold functional , making the problem liable to traditional tolls of bifurcation theory . such strategy was followed in ref .  , resulting in a full description of a spontaneous parity - breaking turbulent bifurcation . the generalization of this model to non - axisymmetric mean state is however an open problem .    summarizing , we have shown that the fluctuations of a turbulent von krmn flow obey a sequence of eckhaus - like instabilities with varying reynolds number that is similar to the sequence of echaus instability observed in the laminar flow , at much lower reynolds number and that may be explained through a suitable parametrization by increasing the dissipation , and therefore decreasing the effective reynolds number ( and the effective number of degrees of freedom ) . if this is indeed true , it may have interesting implications for other fields , especially astrophysics and geophysics where reynolds number are huge and our interest is mainly in the dynamics of the large - scale mean flow . for example , current climate models that have presently very low effective reynolds number may be able to capture efficiently certain types of instabilities of the mean ( climatic ) turbulent state . they may therefore be more predictable that kraichnan thought , in a very different meaning . m.  berhanu , r.  monchaux , s.  fauve , n.  mordant , f.  ptrlis , a.  chiffaudel , f.  daviaud , b.  dubrulle , l.  mari , l.  ravelet , m.  bourgoin , ph . odier , j .- f . pinton , r.  volk , magnetic field reversals in an experimental turbulent dynamo , , 77 , 59001 , 2007                      p .- cortet , e.  herbert , a.  chiffaudel , f.  daviaud , b.  dubrulle , and v.  padilla . susceptibility divergence , phase transition and multistability of a highly turbulent closed flow . , 2011(07):p07012 , 2011 . cortet , p.  diribarne , r.  monchaux , a.  chiffaudel , f.  daviaud , and b.  dubrulle . normalized kinetic energy as a hydrodynamical global quantity for inhomogeneous anisotropic turbulence . , 21(2):02510402510411 , 2009 . b.  legras and b.  villone . stability of turbulent kolmogorov flow . in p. d.  j.  peinke , p.  d.  h.  a.  kittel , d.  s.  barth , and p.  d .- i . m. oberlack , editors , _ progress in turbulence _ , number 101 in springer proceedings in physics , pages 99102 . springer berlin heidelberg , 2005 .
Abstract:

Answer: the notion of instability of a turbulent flow is introduced in the case of a von krmn flow thanks to the monitoring of the spatio - temporal spectrum of the velocity fluctuations , combined with projection onto suitable beltrami modes . it is shown that the large scale coherent fluctuations of the flow obeys a sequence of eckhaus instabilities when the reynolds number @xmath0 is varied from @xmath1 to @xmath2 . this sequence results in modulations of increasing azimuthal wavenumber . the basic state is the laminar or time - averaged flow at an arbitrary @xmath0 , which is axi - symmetric , _ i.e. _ with a @xmath3 azimuthal wavenumber . increasing @xmath0 leads to non - axisymmetric modulations with increasing azimuthal wavenumber from @xmath4 to @xmath5 . these modulations are found to rotate in the azimuthal direction . however no clear rotation frequency can be established until @xmath6 . above , they become periodic with an increasing frequency . we finally show that these modulations are connected with the coherent structures of the mixing shear layer . the implication of these findings for the turbulence parametrization is discussed . especially , they may explain why simple eddy viscosity models are able to capture complex turbulent flow dynamics .