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: asteroseismology is the only systematic technique to study the insides of a star . pulsators are fortunately found all over the h - r diagram , and allow us the opportunity to look inside different stars in various evolutionary stages . white dwarf stars are the stellar remains of 9899% of stars in the sky @xcite , and contain an archival record of their main sequence lifetime . pulsating white dwarf stars serve as effective instruments to harness this archival record . white dwarf spectra reveal that 80% of them have atmospheres dominated by hydrogen ( das ; @xcite ) . hydrogen atmosphere da white dwarfs are observed to pulsate in an instability strip located in the temperature range 1080012500k for @xmath0 @xcite . recent studies determine that the da instability strip is only about 10001200k in width @xcite . these da variables ( davs ) are known as the zz  ceti stars , after the proto type zz  ceti ( r548 ) of the class .    in this paper we gather the pulsation spectra of all known non - interacting zz  ceti stars with the purpose of studying their ensemble pulsation characteristics , and how these properties change across the instability strip . such a systematic study of the zz  ceti pulsators was first undertaken by pioneers in the field such as @xcite , @xcite , and @xcite . @xcite was the first to demonstrate the distinct behavior of pulsation periods and amplitudes as a function of temperature systematically for a significant sample of davs . the sample size of known davs is now three times larger than the sample used in the last such characterization of the zz  ceti pulsators by @xcite . also , we now have better and internally consistent temperatures for these pulsators , calculated using the ml2/@xmath1 convection description ( see * ? ? ? * ; * ? ? ? this prescription for convection gives the best internal consistency between optical and ultra - violet ( uv ) effective temperatures , trigonometric parallaxes , v magnitudes , and gravitational redshifts @xcite . these two reasons compel us to re - examine the pulsation characteristics of the instability strip .      of the 35 new zz  ceti variables published in @xcite , hs0951 + 1312 , hs0952 + 1816 , wd1443 + 0134 , wd1524 - 0030 , and wd2350 - 0054 do not have reliable spectroscopic @xmath2 and @xmath3 fits ( see @xcite ) . our sample of new zz  ceti variables discovered using spectra from the sloan digital sky survey ( sdss ) then consists of 30 davs from @xcite and 11 davs from @xcite . for these 41 davs , d. eisenstein derived a homogeneous and internally consistent set of @xmath2 and @xmath3 values from the sdss using model atmospheres from d. koester ( see * ? ? ? * ) . we carefully and consistently re - analyzed all of our time - series photometry data on these pulsators , acquired using the prime focus ccd camera argos @xcite on the 2.1 m telescope at mcdonald observatory . we shall henceforth refer to this set of 41 davs as the sdss zz  ceti sample with homogeneous spectroscopic fits and homogeneous time - series photometry . @xcite and @xcite have published internally consistent temperatures and @xmath3 fits for 39 davs using their latest model atmospheres ; we will henceforth refer to this second set of 39 davs as the bg04 zz  ceti sample . we compiled a corresponding set of 39 pulsation spectra from the literature , and via private communication with our colleagues . the seismic data of the bg04 zz  ceti sample was acquired by different observers using different instruments and telescopes . however a substantial amount of time - series data exists on most davs in this sample , and we utilized practically all published pulsation spectra to carefully derive well - averaged values of weighted mean period and pulsation amplitudes which we present in this paper . @xcite show evidence of a relative shift of about 200k between the sdss and bg04 zz  ceti instability strips , which also differ in shape and width . the spectra of zz  ceti stars in these samples were analyzed using different techniques and with different model atmospheres . as homogeneity is imperative , we can not merge these two samples for analyses that involve the spectroscopic temperature . however we can merge these two samples when formulating plots based solely on seismic data such as figure 3 .    during the writing of this manuscript , 25 new zz  ceti stars were submitted for publication in separate papers : silvotti et al . 2005a , castanheira et al . 2005 , kepler et al . 2005b , gianninas 2005 , and silvotti et al . this brings the total number of non - interacting zz  ceti stars known to 107 . twenty of the new davs in these papers come from the sdss . their spectroscopic @xmath2 and @xmath3 values were derived using the same technique as the 41 sdss davs discovered previously , but with a different algorithm ( version auto23 vs. auto21 used earlier ) . hence we can not include the new davs in our sdss zz  ceti sample for plots based on effective temperature , such as figures 1 and 2 . also , we find that the main purpose of the pulsation spectra presented in these papers is to demonstrate discovery of variability . such pulsation spectra may be incomplete and are not well suited for direct inclusion in most plots of this paper where ensemble homogeneity plays a big role . however we will include these variables in figure 3 ( based on seismic data ) for completeness , and show the net effect of the inclusion on the original plot . @xcite , @xcite , @xcite , and @xcite show that the non - linear pulse shapes observed in some variable white dwarfs may arise as a result of relatively thick convection zones . many frequencies are evident in these stars because the non - linearities appear as harmonics and linear combinations in our fourier transforms . we are interested in studying the pulsation characteristics of self - excited real modes in the context of this paper . hence we will disregard all harmonics and linear combination frequencies in the observed pulsation spectra throughout this paper . when our program detects a linear combination or harmonic in an observed pulsation spectrum , we typically hand - pick the lowest amplitude mode as the linear combination mode . however if we find a mode involved in two or more linear combinations , we consider it to be a linear combination mode even if it is not the lowest amplitude mode . our simplistic approach could lead us to incorrect choices ; we could be misinterpreting resonances as linear combination modes in a few cases . however any small discrepancies in values evaluated for a single star should have a minimal impact on the conclusions we draw from the ensemble . short pulsation periods ( 100300s ) in the zz  ceti stars typically have pulsation amplitudes of a few percent or smaller . long pulsation periods ( 6001000s ) typically have large amplitudes that can be as high as 10% . this correlation between the pulsation periods and amplitudes of the zz  ceti stars has been established for a long time @xcite . however the correlation of these properties with temperature had to wait more than a decade for sufficiently accurate determinations of zz  ceti temperatures . subsequently the distinct behavior of pulsation periods and amplitudes as a function of temperature was systematically demonstrated for a significant sample of davs by @xcite and more recently by @xcite . we now demonstrate and discuss these trends for the new as well as previously known davs of our two homogeneous samples . we will begin with showing how the number of observed independent modes changes across the instability strip . we expect a larger ( than intrinsic ) scatter in the number of independent modes we measure for an ensemble of zz  ceti stars due to varied detection efficiency . this partially arises because we use different instruments on different telescopes , or observe for different durations of time , and also due to the variety of weather conditions and extinction values . a fraction of the scatter in our detection efficiency is caused because zz  ceti stars have different magnitudes and different pulsation amplitudes , and the pulsation amplitudes depend on the effective wavelength of observation @xcite .    by using the sdss zz  ceti sample alone , we reduce some of this expected scatter ( see subsection 1.1 ) . to reduce the scatter additionally , we choose to limit ourselves to the magnitude range @xmath4 and to a stellar mass @xmath5 to exclude intrinsic low amplitude davs ( see section 2.3 ) . reducing our magnitude range any more will take us further into the realm of small number statistics . for the subset of 23 sdss davs that fall in the chosen range , we find an average of 2.5modes / star for the hotter zz  ceti stars ( 1130011800k ) , and a slightly larger value of 3.1modes / star for the cooler zz  ceti stars ( 1080011300k ) . if we consider all the sdss davs in our sample , excluding the three massive variables ( wd0923 + 0120 , wd1711 + 6541 , & wd2159 + 1322 ) , then we find an average of 2.9modes / star for the hotter zz  ceti stars ( 1130011800k ) and 3.3modes / star for the cooler zz  ceti stars ( 1080011300k ) up to an amplitude limit of order a few mma . this may seem contrary to published literature , but we are attempting to restrict ourselves to independent self - excited modes only , rather than plot the total number of observed periodicities in the star .    the period spacing for nonradial g - modes with @xmath63 is larger than the asymptotic limit . this increase in the density of available eigenmodes as a function of period is consistent with the increase we find in the number of observed independent modes , comparing the hotter zz  ceti stars to the cooler zz  ceti stars . pulsations in zz  ceti models ( and in stars in general ) are self - driven oscillations . stochastic noise frequencies coincident with the eigenfrequencies are amplified by the driving mechanism to observable amplitudes . the blue edge of the zz  ceti instability strip occurs when the star is cool enough to have a hydrogen partial ionization zone sufficiently deep to excite global pulsations . our models suggest that the driving frequency is governed by the thermal timescale at the base of the hydrogen partial ionization zone ( see @xcite ) . as the model star cools , the base of the partial ionization zone moves deeper in the model , and the thermal timescale increases @xcite . numerical calculations by montgomery ( 2005 ) show that the thermal timescale at the base of the ionization zone @xmath7 is proportional to @xmath8 . however recent investigations by @xcite suggest that the thermal timescale is not consistent with the mean pulsation period , although both quantities increase in the models as we traverse from the blue to the red edge of the instability strip . figure 1 shows the mean pulsation period as a function of effective temperature for davs in both the sdss and bg04 zz  ceti samples . we determine the weighted mean period ( wmp ) for each dav by weighting each period with the corresponding observed amplitude ( see equation 1 ) . @xmath9    we can now appreciate the significance of figure 1 ; the weighted mean pulsation period provides us with a method of determining the temperature of a zz  ceti star independent of spectroscopy . the respective equations of the best fit lines through the sdss and the bg04 zz  ceti samples shown in figure 1 are : @xmath10 = -0.835@xmath11 @xmath12 + 10060@xmath13@xmath14 = -0.830@xmath15 @xmath16 + 10240@xmath17    where wmp is in units of seconds and @xmath2 is in units of kelvin . it is reassuring to see an agreement in the values of the slope for both the sdss and bg04 zz  ceti samples , and gives us more confidence in using this period - temperature correlation as a means to measure the effective temperature of a dav . we can now determine a temperature for the five variables with unreliable spectroscopic fits we excluded earlier , hs0951 + 1312 , hs0952 + 1816 , wd1443 + 0134 , wd1524 - 0030 , and wd2350 - 0054 ( see table 1 in the appendix ) . we compute a sum of the squares of the * horizontal * offsets of the points in figure 1 respectively from the above best fit lines , and arrive at a seismic estimate of the average uncertainty in our spectroscopic temperature . we determine @xmath18=188.3k for the sdss zz  ceti sample , and @xmath18=185.4k for the bg04 zz  ceti sample . it is not surprising that the average @xmath2 uncertainties in both these samples turned out to be similar because zz  ceti stars in both samples define an instability strip of the same width @xmath1910001200k@xcite . we can expect that if the average @xmath2 uncertainty in the sdss sample had been significantly larger than that of the bg04 sample , the empirical width of the sdss instability strip would also have been correspondingly larger . our current temperature determinations are getting better and may be internally consistent to about 200k at present . our empirical measure of the width of the instability strip appears to have reduced from 1500k @xcite to about 10001200k @xcite . but even with our current atmospheric models we do find exceptions to this period-@xmath2 correlation , where the spectroscopic temperature does not match the pulsation properties of the zz  ceti star . for example , @xcite report the zz  ceti star wd2350 - 0054 with a dominant period close to 300s as an unusual pulsator because its sdss spectroscopic determination places it below the red edge of the instability strip at @xmath20k . @xcite report a known complex pulsator g29 - 38 @xcite as a hot zz  ceti star with @xmath21k . better s / n spectra of these and other such zz  ceti stars may reveal temperatures consistent with their observed periods . unless the pulsation characteristics themselves reveal curious features , it seems premature to conclude that the zz  ceti star in question ( with a spectroscopic temperature inconsistent with its pulsation spectrum ) is strange or abnormal in any way . measuring the pulsation periods of a zz  ceti star is a model - independent and straight forward process as opposed to a spectroscopic determination of its temperature . the instability strip spans 10001200k in temperature and about 1300s in period . the internal uncertainty in measuring temperature using spectra is typically about 200k , which is 1720% of the width of the instability strip . while even low quality photometric data will yield a period precise to at least a few seconds . for the hotter zz  ceti stars , with little or no amplitude modulation , the uncertainty in wmp constitutes less than 1% of the width of the instability strip in period . cooler zz  ceti stars with a significant amplitude modulation can exhibit pulsation spectra with a wmp different by 3060s and in a few cases even as much as 100200s . for most of the cooler zz  ceti stars , the uncertainty in wmp represents less than 5% of the width of the instability strip . we therefore expect that this method of using the wmp to effectively measure the location of a dav within the instability strip is in general more accurate and reliable than spectroscopy . any dispersion in period due to differences in stellar mass , core composition , h / he layer masses , etc . , can increase the uncertainty of determining the temperature of a given dav using this relation . we expect this may be significantly smaller than our present typical spectroscopic @xmath2 uncertainties of a few 100k . note that we do not claim that the relationship between the wmp and the temperature is linear ; a straight line is merely the simplest model fit possible to the observations shown in figure 1 , considering the large amount of scatter . we do not necessarily require a linear relation between wmp and @xmath2 for this method to work as a temperature indicator . our collaborators at the university of texas are presently investigating the interpretation of these data in terms of the theoretical models ( kim et al . 2006 ; in preparation ) . the plot of wmp vs. @xmath2 is similar to comparing two independent temperature scales , each with its own independent source of uncertainties . the uncertainty in spectroscopic @xmath2 depends on the quality of the spectrum , the accuracy and completeness of the model atmosphere , and the details of the algorithm used to fit the observed spectrum with the template spectra . large amplitude variables have a corresponding higher uncertainty in temperature @xcite . the uncertainty in wmp comes from the quality of the photometric observations , and the amplitude modulation of the star . as long as we restrict our relative @xmath2 parameter in units of seconds in the wmp temperature scale , the uncertainty in our measurements is not related to spectroscopic temperature values at all . it is only when we attempt to translate wmp into a temperature in degree kelvin , that we have to use the relation between wmp and spectroscopic temperature . even in this case , we are still better off than the typical 200k uncertainty in spectroscopic temperature because the slopes of the best fit lines in figure 1 do not depend on a few stars , but on 3540 stars . using wmp directly as a temperature scale is non - intuitive at the present time . but this maybe worth thinking about as an alternative scale in the long run , once we improve our understanding of the relation between wmp and @xmath2 , both observationally and theoretically . we presently classify the zz  ceti stars into two groups : the hot davs and the cool davs . however a clear dividing point ( in temperature ) between these two classes does not exist to date in published literature . our present classification of these stars is thus vague in this context , although it has certainly proved to be a useful guide . the primary reason such a temperature based classification can not be well defined is because the location and width of the instability strip are model - dependent features . model atmospheres of dav stars treat convection with some parameterization , the choice of which can shift the instability strip in temperature by a few thousand k(@xcite ) . the pulsation characteristics of the zz  ceti stars helped us divide them into simple and complex pulsators in the late seventies . a decade later , we recognized that the simple pulsators that exhibit a few modes with short periods ( 200300s ) , small amplitudes ( few mma ) , sinusoidal or saw - toothed pulse shapes , and continued to show the same pulsation spectra over a few decades , were hot zz  ceti stars . we also realized that the complex pulsators that exhibit several long periods ( 6001200s ) with large amplitudes ( 40110mma ) , non - sinusoidal pulse shapes with fast rise times and slow decadence , and amplitude modulation over timescales of a few days to a few years , were cool zz  ceti stars . although we now use spectroscopic temperature as a means to classify these stars , the classification scheme came about only because we initially used the pulsation characteristics to separate these stars into two groups . we suggest a new zz  ceti classification scheme based on the wmp of these variables . we intend to retain the fundamental aspect of the previous scheme in using a temperature based classification . we have shown that wmp is also a temperature scale ( figure 1 ) . we expect that wmp serves as a more accurate and reliable @xmath2 scale than spectroscopic temperature because measuring the pulsation periods of a zz  ceti star constitutes a relatively simple , less uncertain , and model - independent exercise compared to measuring its spectroscopic temperature . we redefine the class of hot davs ( hdavs ) as zz  ceti stars with wmp@xmath22350s . we redefine the class of cool davs ( cdavs ) as zz  ceti stars with wmp@xmath23650s . we suggest introducing a new class of davs , to be called the intermediate davs ( idavs ) , as zz  ceti stars with 350@xmath24wmp@xmath24650 . this class merely forms the evolutionary subclass adjoining the hot and cool zz  ceti stars and typically encompasses those zz  ceti pulsators that show a large range of pulsation periods , e.g. hs0507 + 0435b .    for borderline cases such as g29 - 38 , when one season of observations place the dav in one class , and a second season places it in another class , then we suggest choosing the coolest class of the two possibilities . we find several cases where a cdav or an idav exhibit modes typical of the hotter pulsators , but we have not seen an instance of a hot dav exhibiting a mode typical of the cool zz  ceti stars . we show our suggested classification for most of the known non - interacting zz  ceti stars in tables 1 and 2 . note that the boundaries of 350s and 650s are arbitrary in some sense ; we merely wished to divide the instability strip into three parts and used a period histogram to fine - tune our choice . for simplicity and better readability , tables 1 and 2 in the appendix do not show all of the pulsation spectra we accrued for each star . we show multiple seasons of observations only for those davs that exhibit different frequencies and amplitudes at different times . this also helps us understand the intrinsic changes in weighted mean period and mean pulsation amplitude for such variables . to use such stars in the ensemble , we computed the weighted mean period and mean pulsation amplitude for each season of observations individually , and then included the average value in our analysis . we refer the reader to @xcite , @xcite , @xcite , @xcite , and @xcite for the pulsation spectra of the new davs . the physical quantity of interest in this subsection is the intrinsic pulsation amplitude of modes excited in the zz  ceti star . but the measured amplitudes of observed modes are most likely lower than the intrinsic amplitudes due to geometric cancellation . this effect has three independent causes : disk averaging , inclination angle , and limb darkening . we are unable to resolve the disk of the star from earth . hence the observed amplitude of each pulsation mode is lower due to a disk - averaging effect . this explains why the probability of detecting @xmath25=1 modes is higher than the detection probability for @xmath25=2 and @xmath25=3 modes . the inclination angle dictates the distribution of the bright and dark zones in our view for a given mode . eigenmodes with different @xmath26 values exhibit different cancellation patterns ( see @xcite ) . limb darkening effectively reduces the area of the stellar disk in our view , and this reduction in area depends on wavelength . at uv wavelengths , the increased limb darkening decreases the contribution of zones near the limb . as a result , modes of higher @xmath27 are canceled less effectively in the uv compared to the low @xmath27 modes @xcite . of these three independent causes of geometric cancellation , limb darkening is the only one that works in our favor , and provides us with a mode identification technique @xcite . the intrinsic pulsation amplitude depends on the mass of the star . nonradial g - modes have a non - negligible radial component , the amplitude of which scales with stellar mass and plays a role in dictating the amplitude of the nonradial component . the massive pulsators bpm37093 ( @xmath3=8.81 ) , wd0923 + 0120 ( @xmath3=8.74 ) , wd1711 + 6541 ( @xmath3=8.64 ) , and wd2159 + 1322 ( @xmath3=8.61 ) exhibit low amplitudes as a result of their high gravity , and hence we exclude them from section 2.3 .    figure 2 shows the square root of total power @xmath28 for 38 stars from the sdss sample ( top ) and 38 stars from the bg04 zz  ceti sample ( bottom ) , plotted as a function of their spectroscopic temperature . all of the pulsation amplitudes we report in this paper come from optical whole disk observations . we expect these are mostly low @xmath25 modes . for these reasons , we expect them to be lower than the corresponding intrinsic amplitudes due to geometric cancellation . the points that form the upper envelopes in both panels of figure 2 are then better indicators of the intrinsic amplitude at that temperature . note that in the few cases of under - resolved data , beating of closely spaced periodicities can lead us to determine a relatively smaller or larger amplitude than the observable amplitude . we expect that zz  ceti stars with a high pulsation amplitude have a corresponding higher uncertainty in temperature . for large amplitude variables , the surface temperature changes substantially during a pulsation cycle , by as much as a thousand degrees @xcite . depending on which phase of the pulsation cycle ( typical periods of 6001200s for cool davs ) we acquire the spectra and for how long , our measure of their effective temperature can be incorrect by a few hundred degrees . we can attempt to estimate this uncertainty by obtaining time - series photometry on the star simultaneous with the spectroscopic data .    at any given temperature , the apparently low amplitude variables could be suffering from extensive geometric cancellation . at the same time , the high amplitude variables with the expectedly least geometric cancellation have a proportionally high uncertainty in temperature due to the large intrinsic temperature fluctuation during a pulsation cycle . this makes the task of interpreting figure 2 difficult . there is fortunately a silver lining to this bleak cloud : zz  ceti stars are multi - mode pulsators . if a zz  ceti with 3 independent modes still exhibits a small amplitude , then it is unlikely that we are dealing with an unfavorable inclination angle in all three cases and @xmath26 . however we do not fully understand the mode selection mechanism for different @xmath26 values . also , they exhibit different cancellation patterns @xcite . ] .    instead of traditional histograms , we adopt a gaussian binning technique where each bin serves as a gaussian function instead of a box ( top - hat ) function . this is a better noise averaging technique because all points contribute to the value of each bin . although the contributions of distant points to the value of an individual bin are small , this method may not be appropriate to study sharp local trends . we are interested in slow trends across the width of the instability strip and do not hesitate to use this technique . we choose gaussian bins with @xmath29 that are an infinitesimal 1k apart from each other . we show the result as a solid line in both panels of figure 2 . both samples show an initial increase in pulsation amplitude , and are also suggestive of a decline prior to the red edge .    we established in section 2.2 how the weighted mean pulsation period correlates directly with the effective temperature of the davs for both samples . we now merge both the sdss and the bg04 samples , also including the five davs with unreliable spectroscopic fits ( see section 1.1 ) to plot the pulsation amplitude as a function of the wmp . we no longer have to worry about the internal consistency in their spectroscopic temperatures , and we can use their wmp as a @xmath2 scale . we show these 81 davs as filled circles in figure 3 . we also include the 19 new davs from @xcite , @xcite , and @xcite as filled squares . we have excluded the massive dav wd1337 + 0104 ( @xmath3=8.55 ) from @xcite and the low mass davs he0031 - 5525 ( @xmath3=7.65 ) & wd2135 - 0743 ( @xmath3=7.67 ) from @xcite in figure 3 . figure 3 shows a plot with better statistics than figure 2 due to the larger sample size of 100 davs . we used gaussian bins of width 75s , separated from adjacent bins by an infinitesimal amount of 1s , to produce both the curves shown in figure 3 . the solid line is the histogram determined from the 81 average mass davs of the bg04 and sdss samples , while the dotted line shows the effect of including the new 19 average mass davs . the minor difference between the two lines assures us that our conclusions are robust . we clearly see an initial increase in pulsation amplitude near the blue edge , followed by a gentle rise and then a decline prior to the red edge . we find several reasons why figure 3 looks much more convincing of a decline in amplitude before the red edge compared to figure 2 . changing from a x - axis based on spectroscopic temperature to weighted mean period allows us to include four red edge pulsators with low amplitudes , namely hs0952 + 1816 , wd1443 + 0134 , wd0249 - 0100 , and wd2307 - 0847 . secondly , we get better statistics by combining the two samples that are individually suggestive of a decline in amplitude near the red edge . lastly , we also find that the change in x - axis between figures 2 and 3 moves pulsators left and right , apparently leaving behind cleaner evidence of a decline prior to the red edge . for example , we find that the low amplitude cool dav gd154 moves from 11180k in figure 2 to 1154s in figure 3 , which corresponds to @xmath30k using the slope in the bottom panel of figure 1 as a conversion factor . we have eight pulsators with weighted mean periods of order or greater than 1000s with amplitudes near or less than 25mma . there are no high amplitude pulsators in this range , and the scatter in amplitude near the red edge is much smaller than the center of the zz  ceti instability strip . it is unlikely that unfavorable inclination angles can explain the low amplitudes of all the independent modes of all eight pulsators . @xcite established the observed red edge at 11000k , but did not see any decline in pulsation amplitude near the red edge within the instability strip , although they were searching for it . it is for the first time that we now have clear evidence of a decline in pulsation amplitude within the instability strip , just prior to the zz  ceti red edge . in other words , the star loses pulsation energy before pulsations shut down at the red edge . this work will have implications and new constraints for our pulsation models . we find that the present large sample of non - interacting davs conforms to the well established trend of increasing pulsation periods with decreasing temperature across the instability strip . investigations of the number of observed independent modes show that hot davs have one more mode on average compared to the cool davs . we also confirm the increase in pulsation amplitude near the blue edge and find strong evidence of a decline in amplitude prior to the red edge . @xcite established the observed zz  ceti red edge at 11000k , but did not see any decline in pulsation amplitude near the red edge within the instability strip . we present the first observational evidence that shows the red edge is not an abrupt feature in zz  ceti evolution , and that the star loses pulsation energy before pulsations shut down at the red edge . this work poses new constraints on our pulsation models . support for this work was provided by nasa through the hubble fellowship grant hst - hf-01175.01-a awarded by the space telescope science institute , which is operated by the association of universities for research in astronomy , inc . , for nasa , under contract nas 5 - 26555 . we thank the referee for helping us make this manuscript a better paper . llllllllllll wd0018 + 0031 & hdav & 11700 & 7.93 & 257.9 5.8 & & & & & & + wd0048 + 1521 & idav & 11290 & 8.23 & 615.3 24.8 & & & & & & + wd0102@xmath310032 & cdav & 11050 & 8.24 & 1043.4 15.3 & 926.3 34.7 & 830.3 35.1 & 752.2 19.4 & & & + wd0111 + 0018 & hdav & 11510 & 8.26 & 292.3 21.9 & 255.3 15.6 & & & & & + wd0214@xmath310823 & hdav & 11570 & 7.92 & 348.1 8.4 & 347.1 8.2 & 297.5 16.0 & 263.5 7.1 & & & + wd0318 + 0030 & cdav & 11040 & 8.07 & 844.9 15.3 & 826.4 27.3 & 695.0 8.9 & 587.1 10.6 & 536.1 11.1 & & + wd0332@xmath310049 & cdav & 11040 & 8.25 & 1143.7 7.4 & 938.4 6.7 & 765.0 15.2 & 402.0 4.1 & & & + wd0332@xmath310049 & cdav & 11040 & 8.25 & 770 23 & 910 10 & & & & & + wd0756 + 2020 & hdav & 11710 & 8.01 & 199.5 6.8 & & & & & & + wd0815 + 4437 & idav & 11620 & 7.93 & 787.5 7.2 & 311.7 18.9 & 309.0 10.2 & 258.3 6.8 & & & + wd0818 + 3131 & hdav & 11800 & 8.07 & 253.3 2.9 & 202.3 3.3 & & & & & + wd0825 + 4119 & idav & 11820 & 8.49 & 653.4 17.1 & 611.0 11.2 & & & & + wd0842 + 3707 & hdav & 11720 & 7.73 & 321.1 4.4 & 309.3 18.0 & 212.3 5.2 & & & & & + wd0847 + 4510 & hdav & 11680 & 8.00 & 200.5 7.0 & 123.4 3.0 & & & & & + wd0906@xmath310024 & idav & 11520 & 8.00 & 769.4 26.1 & 618.8 9.3 & 574.5 23.7 & 457.9 9.8 & 266.6 7.6 & & + wd0913 + 4036 & hdav & 11680 & 7.87 & 320.5 14.7 & 288.7 12.4 & 260.3 16.5 & 203.9 3.8 & & & + wd0923 + 0120 & cdav & 11150 & 8.74 & 655.7 4.4 & & & & & & & + wd0923 + 0120 & cdav & 11150 & 8.74 & 668.9 3.5 & & & & & & & + wd0939 + 5609 & hdav & 11790 & 8.22 & 249.9 7.2 & 48.5 5.9 & & & & & & + wd0942 + 5733 & idav & 11260 & 8.27 & 909.4 7.7 & 694.7 37.7 & 550.5 12.3 & 273.0 9.0 & & & & + wd0949@xmath310000 & idav & 11180 & 8.22 & 711.6 6.0 & 634.2 5.1 & 516.6 16.2 & 365.2 17.7 & 364.1 7.3 & 363.2 12.5 & & + wd0958 + 0130 & hdav & 11680 & 7.99 & 264.4 4.7 & 203.7 2.5 & 121.2 1.6 & & & & & + wd1002 + 5818 & hdav & 11710 & 7.92 & 304.6 5.3 & 268.2 6.8 & & & & & + wd1007 + 5245 & hdav & 11430 & 8.08 & 323.1 10.4 & 290.1 7.7 & 258.8 11.0 & & & & + wd1015 + 0306 & hdav & 11580 & 8.14 & 270.0 8.4 & 255.7 7.3 & 194.7 5.8 & & & & & + wd1015 + 5954 & idav & 11630 & 8.02 & 769.9 5.7 & 453.8 15.5 & 401.7 21.4 & 292.4 8.6 & & & & + wd1015 + 5954 & idav & 11630 & 8.02 & 768.4 7.5 & 455.3 17.7 & 401.4 19.8 & 294.0 9.1 & & & & + wd1015 + 5954 & idav & 11630 & 8.02 & 456.1 19.6 & 399.7 19.2 & 145.5 3.1 & 139.2 4.6 & & & & + wd1054 + 5307 & cdav & 11120 & 8.01 & 869.1 37.4 & & & & & & & + wd1056@xmath310006 & cdav & 11020 & 7.86 & 1024.9 31.6 & 925.4 60.3 & 670.6 12.0 & 603.0 11.5 & & & & + wd1122 + 0358 & cdav & 11070 & 8.06 & 996.1 17.3 & 859.1 34.6 & 740.1 10.0 & & & & & + wd1125 + 0345 & hdav & 11600 & 7.99 & 335.1 2.8 & 265.8 3.3 & 265.5 7.1 & 208.6 2.8 & & & & + wd1157 + 0553 & cdav & 11050 & 8.15 & 1056.2 5.8 & 918.9 15.9 & 826.2 8.1 & 748.5 5.6 & & & + wd1345@xmath310055 & hdav & 11800 & 8.04 & 254.4 2.4 & 195.5 3.9 & 195.2 5.5 & & & & & + wd1354 + 0108 & hdav & 11700 & 8.00 & 322.9 1.9 & 291.6 2.2 & 198.3 6.0 & 173.3 1.1 & 127.8 1.5 & & + wd1355 + 5454 & hdav & 11580 & 7.95 & 324.0 21.8 & & & & & & + wd1417 + 0058 & cdav & 11300 & 8.04 & 980.0 11.3 & 894.6 42.8 & 812.5 32.1 & 749.4 17.9 & 522.0 14.9 & & + wd1502@xmath310001 & idav & 11200 & 8.00 & 687.5 12.0 & 629.5 32.6 & 581.9 11.1 & 418.2 14.9 & & & & + wd1617 + 4324 & cdav & 11190 & 8.03 & 889.7 36.6 & 661.7 21.2 & 626.3 15.4 & & & & & + wd1700 + 3549 & cdav & 11160 & 8.04 & 1164.4 11.4 & 955.3 20.3 & 893.4 54.3 & 552.6 9.3 & & & + wd1711 + 6541 & cdav & 11310 & 8.64 & 1248.2 3.2 & 690.2 3.3 & 606.3 5.2 & 234.0 1.2 & & & & + wd1711 + 6541 & cdav & 11310 & 8.64 & 1186.6 3.3 & 934.8 2.9 & 612.6 5.7 & 561.5 3.0 & 214.3 1.7 & & & + wd1724 + 5835 & hdav & 11540 & 7.89 & 337.9 5.9 & 279.5 8.3 & 189.2 3.2 & & & & & + wd1732 + 5905 & cdav & 10860 & 7.99 & 1248.4 4.5 & 1122.4 8.0 & & & & & & + wd1732 + 5905 & cdav & 10860 & 7.99 & 1336 7.8 & 1090 8.0 & & & & & & + wd2159 + 1322 & cdav & 11710 & 8.61 & 801.0 15.1 & 683.7 11.7 & & & & & & + wd2214@xmath310025 & hdav & 11440 & 8.33 & 255.2 13.1 & 195.2 6.1 & & & & & & + wd1443 + 0134 & cdav & 10830 & & 1085.0 5.2 & 969.0 7.5 & & & & & & + wd1524@xmath310030 & cdav & 11160 & & 873.3 110.8 & 717.5 28.3 & 498.6 21.6 & 255.2 17.9 & & & + wd2350@xmath310054 & hdav & 11690 & & 391.1 3.1 & 304.1 16.3 & 272.8 16.2 & & & & & + wd2350@xmath310054 & hdav & 11710 & & 304.5 13.8 & 272.7 14.8 & 212.6 3.0 & & & & & + wd2350@xmath310054 & hdav & 11680 & & 391.1 7.5 & 304.3 17.0 & 273.3 6.3 & 206.7 3.2 & & & & + hs0951 + 1312 & hdav & 11740 & & 311.7 2.7 & 282.2 9.0 & 258.0 3.6 & 208.0 9.4 & & & & + hs0952 + 1816 & cdav & 10800 & & 1160.9 7.9 & 945.9 10.4 & & & & & & + hs0952 + 1816 & cdav & 10960 & & 1150 4.8 & 883 3.6 & 790 2.9 & 674.7 3.0 & & & & +    llllllllllllllll bpm30551 & cdav & 11260 & 8.23 & 920.5 18.5 & 741.4 21.7 & 655.4 17.4 & 442.8 6.5 & & & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 993.0 6.5 & 936.2 13.0 & 885.6 15.2 & 819.2 19.5 & 738.0 6.5 & 609.6 6.5 & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 963.8 13.0 & 844.5 19.5 & 799.2 10.9 & 682.7 13.0 & 606.8 13.0 & & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 920.4 8.7 & 751.6 13.0 & 682.7 10.9 & 606.8 15.2 & 546.1 5.4 & 496.5 4.3 & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 910.2 8.7 & 862.3 6.5 & 731.4 6.5 & 682.7 13.0 & 607.9 17.4 & & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 1129.9 10.9 & 1057.0 8.7 & 744.7 18.5 & 606.8 16.3 & & & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 1137.8 8.7 & 958.1 7.6 & 862.3 16.3 & 744.7 7.6 & 606.8 14.1 & & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 1092.3 5.4 & 993.0 8.7 & 936.2 9.8 & 744.7 7.6 & 712.3 6.5 & 606.8 11.9 & & & & & & @xcite + bpm30551 & cdav & 11260 & 8.23 & 862.3 6.5 & 799.2 8.7 & 744.7 7.6 & 668.7 7.6 & 612.5 5.4 & & & & & & & @xcite + r548 & hdav & 11990 & 7.97 & 333.6 0.6 & 318.0 0.9 & 274.8 3.8 & 274.3 4.8 & 213.1 7.4 & 212.8 4.7 & 187.3 0.9 & & & & & @xcite + mct0145 - 2211 & idav & 11550 & 8.14 & 823.2 17 & 727.9 19 & 462.2 25 & & & & & & & & & @xcite + bpm31594 & idav & 11540 & 8.11 & 617.9 48.0 & 401.6 16 & 416.1 5 & & & & & & & & & @xcite + hltau76 & idav & 11450 & 7.89 & 933 25.2 & 796 9.7 & 781 9.9 & 541 40.5 & 494 30.2 & 383 21.8 & & & & & & @xcite + g38 - 29 & cdav & 11180 & 7.91 & 1024.0 26.1 & 938.0 26.5 & & & & & & & & & & @xcite + g38 - 29 & cdav & 11180 & 7.91 & 1019.8 12.2 & 910.3 28.3 & & & & & & & & & & @xcite + g191 - 16 & cdav & 11420 & 8.05 & 892.9 100 & & & & & & & & & & & @xcite + hs0507 + 0435b & idav & 11630 & 8.17 & 743.0 13.9 & 557.7 18.7 & 446.2 11.0 & 444.8 13.6 & 355.8 22.7 & 354.9 6.8 & 286.1 3.6 & & & & & @xcite + hs0507 + 0435b & idav & 11630 & 8.17 & 743.4 8.3 & 588.7 3.9 & 583.8 1.6 & 559.6 2.7 & 557.6 17.4 & 557.2 3.1 & 556.5 7.8 & 555.3 18.0 & 446.1 15.1 & 445.3 3.0 & 444.6 12.9 & @xcite + & & & & 355.8 26.1 & 355.4 4.7 & 354.9 11.1 & & & & & & & & & + gd66 & hdav & 11980 & 8.05 & 649.4 2.2 & 441.9 1.6 & 301.7 6.7 & 271.7 8.4 & 196.5 3.6 & & & & & & & @xcite + gd66 & hdav & 11980 & 8.05 & 788.6 5.5 & 461.0 3.6 & 301.7 6.3 & 271.7 30 & 197.2 2.2 & & & & & & & @xcite + gd66 & hdav & 11980 & 8.05 & 304.5 8.8 & 271.1 14.8 & 256.5 9 & 197.8 7 & 123.1 2.3 & & & & & & & @xcite + he0532 - 5605 & idav & 11560 & 8.49 & 688.8 8.3 & 586.4 7.9 & & & & & & & & & & @xcite + kuv08368 + 4026 & idav & 11490 & 8.05 & 618.0 17.4 & 494.5 6.0 & & & & & & & & & & @xcite + gd99 & cdav & 11820 & 8.08 & 1151.0 1.9 & 1088.0 4.3 & 1058.0 8.3 & 1007.0 6.5 & 976.0 2.1 & 924.7 1.7 & 853.2 2.4 & 633.1 2.0 & 228.9 4.5 & 223.6 2.9 & 105.2 2.0 & chynoweth&thompson2005;prv.comm . + g117-b15a & hdav & 11630 & 7.97 & 304.4 8.2 & 271.0 7.3 & 215.2 23.9 & & & & & & & & & @xcite + kuv11370 + 4222 & hdav & 11890 & 8.06 & 462.9 3.5 & 292.2 2.7 & 257.2 5.8 & & & & & & & & & @xcite + g255 - 2 & cdav & 11440 & 8.17 & 898.5 9.4 & 855.4 11.2 & 819.7 11.3 & 775.2 15.2 & 681.2 24.9 & 607.9 13.1 & 568.5 6.6 & & & & & vauclair2005;prv.comm . + g255 - 2 & cdav & 11440 & 8.17 & 985.2 4.8 & 773.4 12.7 & 681.2 27.7 & 568.5 16.5 & & & & & & & & vauclair2005;prv.comm . + bpm37093 & idav & 11730 & 8.81 & 636.7 1.7 & 633.2 1.1 & 613.5 1.1 & 600.7 0.9 & 582.0 1.0 & 565.5 1.2 & 562.6 0.9 & 548.4 1.1 & & & & @xcite + bpm37093 & idav & 11730 & 8.81 & 660.8 0.5 & 637.2 0.7 & 633.5 1.3 & 565.9 0.5 & 549.2 0.8 & 531.1 1.2 & 511.7 0.7 & & & & & @xcite + he1258 + 0123 & cdav & 11410 & 8.04 & 1092.1 14 & 744.6 23 & 528.5 9 & & & & & & & & & @xcite + gd154 & cdav & 11180 & 8.15 & 1190.5 6.3 & 1186.5 16.7 & 1183.5 4.6 & 1092.1 3.0 & 1088.6 5.0 & 1084.0 5.6 & & & & & & @xcite + lp133 - 144 & hdav & 11800 & 7.87 & 327.3 4 & 306.9 4 & 304.5 4 & 209.2 10 & & & & & & & & @xcite + g238 - 53 & hdav & 11890 & 7.91 & 206.2 8 & & & & & & & & & & & @xcite + ec14012 - 1446 & idav & 11900 & 8.16 & 937 11 & 610 57 & 724 21 & 530 15 & 399 13 & & & & & & & @xcite + gd165 & hdav & 11980 & 8.06 & 249.7 0.7 & 192.8 0.9 & 192.7 2.4 & 192.6 1.9 & 166.2 0.4 & 146.4 0.5 & 120.4 1.8 & 120.36 4.8 & 120.3 1.4 & 114.3 0.6 & 107.7 0.4 & @xcite + l19 - 2 & hdav & 12100 & 8.21 & 350.1 1.1 & 348.7 0.5 & 192.6 6.5 & 193.1 0.9 & 192.1 0.8 & 143.4 0.6 & 143.0 0.3 & 118.9 0.3 & 118.7 1.2 & 118.5 2.0 & & @xcite + & & & & 114.2 0.3 & 113.8 2.4 & 113.3 0.6 & & & & & & & & + pg1541 + 651 & cdav & 11600 & 8.10 & 757 14 & 689 45 & 564 15 & 467 3 & & & & & & & & @xcite + r808 & cdav & 11160 & 8.04 & 833 81 & & & & & & & & & & & @xcite + g226 - 29 & hdav & 12270 & 8.28 & 109.5 2.8 & 109.3 1.1 & 109.1 2.5 & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 1176 22.6 & & & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 1050 9.1 & & & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 1086 13.2 & & & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 978 7.7 & & & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 1098 6.1 & & & & & & & & & & & @xcite + bpm24754 & cdav & 11070 & 8.03 & 1122 6.7 & & & & & & & & & & & @xcite + g207 - 9 & idav & 11950 & 8.35 & 557.4 6.3 & 318.0 6.4 & 292.0 5.0 & 259.1 1.7 & & & & & & & & @xcite + g185 - 32 & hdav & 12130 & 8.05 & 651.7 0.7 & 537.6 0.6 & 454.6 0.4 & 370.2 1.6 & 301.4 1.1 & 299.8 1.0 & 264.2 0.5 & 215.7 1.9 & 141.9 1.4 & 72.9 0.4 & & @xcite + g185 - 32 & hdav & 12130 & 8.05 & 370.2 1.3 & 301.6 1.5 & 215.7 1.9 & 141.9 1.5 & 72.6 0.7 & & & & & & & @xcite + gd385 & hdav & 11710 & 8.04 & 256.3 10.9 & 256.1 11.4 & & & & & & & & & & @xcite + gd244 & hdav & 11680 & 8.08 & 307.0 14 & 294.6 5.5 & 256.3 5 & 203.3 10.5 & & & & & & & & @xcite + pg2303 + 243 & cdav & 11480 & 8.09 & 900.5 16 & 794.5 56 & 675.4 8 & 623.4 15 & 570.7 8 & & & & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 1015.5 14.5 & 930.9 25.7 & 824.7 20.2 & 677.0 17.6 & 612.9 20.0 & & & & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 934.5 20.5 & 813.8 23.5 & 671.5 23.0 & 623.8 11.8 & & & & & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 859.6 24.6 & 648.7 7.8 & 614.3 31.3 & 498.3 5.8 & 401.3 4.4 & 400.4 7.0 & 399.6 8.8 & 283.9 3.5 & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 915.4 5.9 & 615.1 58.0 & 500.4 8.0 & 474.9 4.8 & 401.3 5.6 & 400.4 5.8 & 399.6 10.0 & 354.9 2.9 & 333.9 2.7 & 283.9 3.5 & 110.1 0.9 & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 770.7 1.5 & 503.5 8.6 & 495.0 11.8 & 401.2 9.7 & 400.5 1.3 & 399.7 4.9 & 283.9 6.4 & 177.1 0.8 & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 809.4 30.1 & 610.3 10.6 & 401.2 11.2 & 399.7 4.5 & 283.9 6.9 & & & & & & & @xcite + g29 - 38 & cdav & 11820 & 8.14 & 894.0 14.0 & 770.8 8.7 & 678.4 9.7 & 612.4 31.6 & 610.7 8.2 & 608.9 8.5 & 551.9 4.4 & 498.3 6.1 & 401.2 6.0 & 399.7 5.7 & 237.0 1.8 & @xcite + & & & & 236.4 1.8 & & & & & & & & + g30 - 20 & cdav & 11070 & 7.95 & 1068 13.8 & & & & & & & & & & & @xcite + ec23487 - 2424 & cdav & 11520 & 8.10 & 993.0 37.7 & 989.3 11.0 & 868.2 12.8 & 804.5 19.3 & & & & & & & & @xcite +
Abstract:

Answer: we present the observed pulsation spectra of all known non - interacting zz  ceti stars ( hydrogen atmosphere white dwarf variables ; davs ) and examine changes in their pulsation properties across the instability strip . we confirm the well established trend of increasing pulsation period with decreasing effective temperature across the zz  ceti instability strip . we do not find a dramatic order of magnitude increase in the number of observed independent modes in zz  ceti stars , traversing from the hot to the cool edge of the instability strip ; we find that the cool davs have one more mode on average compared to the hot davs . we confirm the initial increase in pulsation amplitude at the blue edge , and find strong evidence of a decline in amplitude prior to the red edge . we present the first observational evidence that zz  ceti stars lose pulsation energy just before pulsations shut down at the empirical red edge of the instability strip .


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: it is an important task to visualize , understand and then extract useful information in complex continuous dynamical systems from many science and engineering fields including flight simulation @xcite , wave propagation @xcite , bio - inspired fluid flow motions @xcite , ocean current model @xcite , and etc . it is , therefore , natural to see many existing tools to study dynamical systems . for example , a lot of work in dynamical systems focuses on understanding different types of behaviors in a model , such as elliptical zones , hyperbolic trajectories , chaotic attractors , mixing regions , to name just a few . one interesting tool is to partition the space - time domain into subregions based on certain quantity measured along with the passive tracer advected according to the associated dynamical system . because of such a lagrangian property in the definition , the corresponding partition is named the lagrangian coherent structure ( lcs ) . one commonly used quantity is the so - called finite time lyapunov exponent ( ftle ) @xcite which measures the rate of separation of a passive tracer with an infinitesimal perturbation in the initial condition , over a finite period of time . another popular diagnostic of trajectory separation in dynamical systems is the finite - size lyapunov exponent ( fsle ) @xcite which is especially popular in the applications from the oceanography . instead of introducing an infinitesimal perturbation in the initial condition of the fluid particles , fsle considers finite perturbation . it does not measure the rate of separation but computes the time required to separate two adjacent tracers up to certain distance . since all these quantities are long treated as a lagrangian property of a continuous dynamical system , most , if not all , numerical methods are developed based on the traditional lagrangian ray tracing method by solving the ode system using any well - developed numerical integrator . these approaches , however , require the velocity field defined at arbitrary locations in the whole space depending on the location of each individual particle . this implies that one has to in general implement some interpolation routines in the numerical code . unfortunately , it could be numerically challenging to develop an interpolation approach which is computationally cheap , monotone , and high order accurate . an eulerian approach to the flow map computation for the ftle has been first proposed in @xcite which has incorporated the level set method @xcite . based on the phase flow method @xcite , we have developed a backward phase flow method for the eulerian ftle computations in @xcite . based on these partial differential equation ( pde ) based algorithms , more recently , we have developed in @xcite an efficient eulerian numerical approach for extracting invariant sets in a continuous dynamical system in the extended phase space ( the @xmath1 space ) . in @xcite , we have proposed a simple algorithm called vials which determines the average growth in the surface area of a family of level surfaces . this eulerian tool relates closely to the ftle and provides an alternative for understanding complicated dynamical systems .    in this paper , we are going to further develop these eulerian tools for continuous dynamical systems . there are several main contributions of this paper . we will first improve the numerical algorithm for _ forward _ flow map construction . the papers @xcite have proposed to construct the _ forward _ flow map defined on a fixed cartesian mesh by solving the level set equations or the liouville equations in the _ backward _ direction . in particular , to compute the forward flow map from the initial time @xmath2 to the final time @xmath3 , one needs to solve the liouville equations _ backward _ in time from @xmath3 to @xmath2 with the initial condition given at @xmath3 . this implementation is numerically inconvenient , especially when incorporating with some computational fluid dynamic ( cfd ) solvers , since the velocity field is loaded from the current time @xmath3 backward in time to the initial time . this implies that the whole field at all time steps has to be stored in the desk which might not be practical at all . in this paper , we first propose an eulerian approach to compute the _ forward _ flow map _ on the fly _ so that the pde is solved _ forward _ in time . developed using this improved algorithm , we propose a simple eulerian approach to compute the isle . for each of the separation factor @xmath0 in the definition of the isle , typical lagrangian method requires to shoot an individual set of rays and to monitor the variations in these trajectories . if the scale @xmath0 is changed , these methods will require to restart the whole computations all over again . in this work , we are going to develop an efficient algorithm which provides the isle field for any arbitrary separation factor @xmath0 . for an individual value of @xmath0 , the method requires to extract only an isosurface of a volumetric data . this can be done easily by any well - developed interpolation algorithm or simply the function isosurface in matlab . another main contribution of this work is to provide a theoretical link between the ftle and the isle . even though these two quantities are measuring different properties of the flow , it has been reported widely , such as @xcite , that they give visually very similar solutions in many examples . we are going to show in particular that their ridges can be identified through the ridge of the largest eigenvalue of the associated cauchy - green deformation tensor of the flow . to the best of our knowledge , this theoretical result is the first one quantitatively revealing the relationship between the ftle ridges and the isle ridges . indeed , comparisons between ftle and fsle have been made in various studies . the paper @xcite has pointed out that the ftle might be unable to recognize the boundaries between the chaotic and the large - scale mixing regime . the article @xcite has argued that the fsle has several limitations in lagrangian coherence detection including aspects from local ill - posedness , spurious ridges and intrinsic jump - discontinuities . nevertheless , the work @xcite has stated that these two concepts might yield similar results , if properly calibrated , and could be even interchangeable in many flow visualizatons such as oceanography and chemical kinetics . in this paper , however , we would like to emphasize on the computational aspect of using these lyapunov exponents , rather than on comparing the advantages and the disadvantages of these different methods . the preference for one or the other could be based on the tradition in various fields , or the availability of the computational resources . this paper is organized as follows . in section [ sec : background ] we will give a summary of several important concepts including the ftle , fsle and isle , and also our original eulerian formulations for computing the flow maps and the ftle . then , our proposed eulerian algorithms will be given in section [ sec : proposal ] . we will then point out some theoretical properties and relationships between the ftle and the isle in section [ sec : relation ] . finally , some numerical examples will be given in section [ sec : examples ] . in this section , we will summarize several useful concepts and methods , which will be useful for the developments that we are proposing . we first introduce the definition of various closely related lyapunov exponents including the finite time lyapunov exponent ( ftle ) @xcite , the finite size lyapunov exponent ( fsle ) @xcite and also the infinitesimal size lyapunov exponent ( isle ) @xcite . then we discuss typical lagrangian approaches and summarize the original eulerian approaches as in @xcite . the discussions here , however , will definitely not be a complete survey . we refer any interested readers to the above references and thereafter .      in this paper , we consider a continuous dynamical system governed by the following ordinary differential equation ( ode ) @xmath4 with the initial condition @xmath5 . the velocity field @xmath6 is a time dependent lipschitz function where @xmath7 is a bounded domain in the @xmath8-dimensional space . to simplify the notation in the later sections , we collect the solutions to this ode for all initial conditions in @xmath9 at all time @xmath10 and introduce the flow map @xmath11 such that @xmath12 represents the arrival location @xmath13 of the particle trajectory satisfying the ode ( [ eqn : ode ] ) with the initial condition @xmath14 at the initial time @xmath15 . this implies that the mapping will take a point from @xmath16 at @xmath15 to another point @xmath13 at @xmath17 .    the finite time lyapunov exponent ( ftle ) @xcite measures the rate of separation between adjacent particles over a finite time interval with an infinitesimal perturbation in the initial location . mathematically , consider the initial time to be 0 and the final time to be @xmath18 , we have the change in the initial infinitesimal perturbation given by @xmath19 the leading order term of the magnitude of this perturbation is given by @xmath20^ * \nabla \phi_0^t(\mathbf{x } ) \delta \mathbf{x}(0 ) \right > } \ , .\ ] ] with the cauchy - green strain tensor @xmath21^ * \nabla \phi_0^t(\mathbf{x})$ ] , we obtain the maximum deformation @xmath22 } } \| \mathbf{e}(0)\| = e^ { \sigma_0^t({\mathbf{x } } ) |t| } \|\mathbf{e}(0 ) \| \ , , \ ] ] where @xmath23 is the eigenvector associated with the largest eigenvalue of the deformation tensor . using this quantity , the finite - time lyapunov exponent ( ftle ) @xmath24 is defined as @xmath25 } = \frac{1}{|t| } \ln \sqrt{\lambda_0^t({\mathbf{x } } ) } \ , . \label{eqn : ftle}\ ] ] where @xmath26 denotes the largest eigenvalue of the cauchy - green tensor . the absolute value of @xmath18 in the expression reflects the fact that we can trace the particles either _ forward _ or _ backward _ in time . in the case when @xmath27 , we are measuring the maximum stretch _ backward _ in time and this corresponds to the maximum compression _ forward _ in time . to distinguish different measures , we call @xmath24 the _ forward _ ftle if @xmath28 and the _ backward _ ftle if @xmath27 . a related concept is a widely used concept in oceanography called the finite size lyapunov exponent ( fsle ) which tries to measure the time it takes to separate two adjacent particles up to certain distance @xcite . there are two length scales in the original definition of the fsle . one is the initial distance between two particles which is denoted by @xmath29 . the other one is the so - called separation factor , denoted by @xmath30 . for any trajectory with initial location @xmath31 with @xmath32 , we first determine the shortest time @xmath33 so that @xmath34 the fsle at the point @xmath35 is then defined as @xmath36    numerically , on the other hand , the maximization over the constraint @xmath32 might require special considerations and the corresponding optimization problem might not be easily solved by typical optimization algorithms . one simple approximation is to consider only @xmath37-neighbors of @xmath35 with each point @xmath38 obtained by perturbing one coordinate of @xmath35 with distance @xmath29 .    to avoid introducing two separate length scales ( @xmath29 and @xmath0 ) and the approximation by @xmath37 neighbors in the final optimization step , one theoretically takes the limit as the length scale @xmath29 tends to 0 . the fsle can then be reduced to the infinitesimal size lyapunov exponent ( isle ) as defined in @xcite given by @xmath39 where @xmath40 is the shortest time for which @xmath41 . even with tremendous theoretical development , there are not many discussions on the numerical implementations . because the quantity is defined using lagrangian particle trajectories , most numerical algorithms for computing the ftle or the fsle are developed based on the ray tracing method which tries to solve the flow system ( [ eqn : ode ] ) using a numerical integrator . in the two dimensional case , for example , the flow map is first determined by solving the ode system with the initial condition imposed on a uniform mesh @xmath42 . then one can apply simple finite difference to determine the deformation tensor and , therefore , its eigenvalues . the computation of the ftle is a little easier . for a given final time , one simply solves the ode system up to that particular time level . then the ftle at each grid location can be computed directly from these eigenvalues . the isle , on the other hand , is a little more involved since for a particular separation factor , @xmath0 , one has to continuously monitor the change in these eigenvalues and record the precise moment when it reaches the thresholding value . computational complexity could be one major concern of these numerical tools . let @xmath43 and @xmath44 be the number of grid points in each physical direction and in the temporal direction , respectively . the total number of operations required to compute the flow map @xmath45 is , therefore , @xmath46 . note however that such complexity is needed for determining the ftle or the isle on one single time level . for all time levels , the computational complexity for these usual lagrangian implementation is @xmath47 . there are several nice ideas to improve the computational time . one natural idea is to concentrate the computations near the locations where the solutions have large derivatives by implementing an adaptive mesh strategy @xcite . another interesting approach is to avoid the computation of the overall map @xmath48 but to numerically decompose it into sub - maps @xmath49 @xcite . then , the extra effort from @xmath50 to @xmath51 , given by @xmath52 are simply some interpolations . all these proposed improvements are only on the flow map computations , and therefore the ftle . for the fsle or the isle computations on one single time level and one particular @xmath0 , say @xmath53 , one has to go through all @xmath54 mesh points in the temporal direction . for a different time level or even a different separation factor @xmath0 , we have to go though the same procedure all over again . we do not aware of many efficient algorithms for the fsle or the isle implementation . ( a ) @xcite . ( a ) lagrangian ray tracing from a given grid location @xmath55 at @xmath2 . note that @xmath56 might be a non - grid point . ( b ) eulerian values of @xmath57 at a given grid location @xmath56 at @xmath3 gives the corresponding take - off location at @xmath2 . note the take - off location might not be a mesh point.,title="fig:",width=217 ] ( b ) @xcite . ( a ) lagrangian ray tracing from a given grid location @xmath55 at @xmath2 . note that @xmath56 might be a non - grid point . ( b ) eulerian values of @xmath57 at a given grid location @xmath56 at @xmath3 gives the corresponding take - off location at @xmath2 . note the take - off location might not be a mesh point.,title="fig:",width=217 ]    we briefly summarize the eulerian approach based on the level set method and the louville equation . we refer interested readers to @xcite and thereafter . we define a vector - valued function @xmath58 . at @xmath2 , we initialize these functions by @xmath59 these functions provide a labeling for any particle in the phase space at @xmath2 . in particular , any particle initially located at @xmath60 in the extended phase space can be * implicitly * represented by the intersection of @xmath8 codimension-1 surfaces represented by @xmath61 in @xmath62 . following the particle trajectory with @xmath63 as the initial condition in a given velocity field , any particle identity should be preserved in the lagrangian framework and this implies that the material derivative of these level set functions is zero , i.e. @xmath64 this implies the following level set equations , or the liouville equations , @xmath65 with the initial condition ( [ eqn : levelsetequationic ] ) . the above * implicit * representation embeds all path lines in the extended phase space . for instance , the trajectory of a particle initially located at @xmath66 can be found by determining the intersection of @xmath8 codimension-1 surfaces represented by @xmath67 in the extended phase space . furthermore , the forward flow map at a grid location @xmath63 from @xmath2 to @xmath3 is given by @xmath68 where @xmath56 satisfies @xmath69 . note that , in general , @xmath56 is a non - mesh location . the typical two dimensional scenario is illustrated in figure [ fig : forwardbackward ] ( a ) . the solution to ( [ eqn : levelsetequation ] ) contains much more information than what was referred to above . consider a given mesh location @xmath56 in the phase space at the time @xmath3 , as shown in figure [ fig : forwardbackward ] ( b ) , i.e. @xmath70 in the extended phase space . as discussed in our previous work , these level set functions @xmath71 defined on a uniform cartesian mesh in fact give the backward flow map from @xmath3 to @xmath2 , i.e. @xmath72 . moreover , the solution to the level set equations ( [ eqn : levelsetequation ] ) for @xmath73 provides also backward flow maps for all intermediate times , i.e. @xmath74 .    to compute the forward flow map , on the other hand , @xcite has proposed to simply reverse the above process by initializing the level set functions at @xmath3 by @xmath75 and solving the corresponding level set equations ( [ eqn : levelsetequation ] ) backward in time . a typical algorithm of this type is given in algorithm 1 . note that in the original lagrangian formulation , the attracting lcs or the unstable manifold is obtained by backward time tracing , while the repelling lcs or the stable manifold is computed by forward time integration . in this current eulerian formulation , on the other hand , forward time integration of the liouville equations gives the attracting lcs and backward time marching provides the repelling lcs . ( 4,0 ) ( -0.5,0)(1,0)12     +    algorithm 1 : computing the forward flow map @xmath76 :    1 . discretize the computational domain to get @xmath77 . 2 .   initialize the level set functions on the last time level @xmath78 by @xmath79 3 . solve the liouville equations for each individual level set function @xmath80 @xmath81 from @xmath82 down to @xmath2 using any well - developed high order numerical methods like weno5-tvdrk2 @xcite with the boundary conditions @xmath83 where @xmath84 is the outward normal of the boundary . 4 .   assign @xmath85 . ( 4,0 ) ( -0.5,0)(1,0)12 in this section , we introduce two new eulerian algorithms for flow visualizations . the first algorithm is to improve the eulerian numerical method for computing the forward flow map . instead of solving the pde _ backward _ in time as in section [ subsec : oldeulerian ] and @xcite , we develop an eulerian pde algorithm to construct the forward flow map from @xmath2 to the final time @xmath3 so that we do not need to load the time - dependent velocity field from the terminal time level backward to the initial time . in other words , we construct the forward flow map _ on the fly_. this is computationally more natural and convenient . because of this simple algorithm , we propose to develop an efficient and systematic pde based method for computing the isle field . typical challenge in obtaining such field is that the choice of the separation factor @xmath0 is in general flow dependent . one has to vary @xmath0 to extract the necessary information . usual lagrangian methods unfortunately require shooting different sets of initial rays for each individual separation factor @xmath0 . this is computationally very inefficient . the proposed eulerian formulation requires either some simple interpolations or only one single level surface extraction which can be easily implemented using any well - developed contour extraction function such as in . one disadvantage about the previous eulerian approach for forward flow map computation as discussed in section [ subsec : oldeulerian ] and @xcite is that the level set equation @xmath86 has to be first solved _ backward _ in time from @xmath3 to @xmath2 . once we have the final solution at the time level @xmath2 , one can then identity the flow map @xmath76 by @xmath87 . this could be inconvenient especially when we need to access the intermediate forward flow map . in this section , we propose a new algorithm to construct the forward flow map _ on the fly_.    consider the forward flow map from @xmath2 to @xmath3 . suppose the time domain @xmath88 $ ] is discretized by @xmath89 discrete points @xmath90 , where @xmath91 and @xmath92 . then on each subinterval @xmath93 $ ] for @xmath94 , we use the same method as described in section [ subsec : oldeulerian ] to construct the forward flow map @xmath49 . in particular , we solve the level set equation ( [ eqn : levelset ] ) _ backward _ in time from @xmath95 to @xmath96 with the terminal condition @xmath97 imposed on the time level @xmath95 . then the forward flow map is given by @xmath98 . once we have obtained this one step _ forward _ flow map @xmath49 , the _ forward _ flow map from @xmath2 to @xmath3 can be obtained using the composition @xmath99 . and this can be easily done by typical numerical interpolation methods . to prevent extrapolation , we could enforce @xmath100 \times [ y_{\min},y_{\max}]$ ] .    here we emphasize that even though we solve equation ( [ eqn : levelset ] ) _ backward _ in time for each subinterval , we access the velocity field data @xmath101 prior to @xmath102 and therefore the forward flow map is indeed obtained _ on the fly_. as a simple example , we consider the first order euler method in the temporal direction . the forward flow map @xmath49 can be obtained by first solving @xmath103 with the terminal condition @xmath104 , and then assigning @xmath105 . higher order generalization is straight - forward . for example , if we use the tvd - rk2 method in the temporal direction @xcite , we can compute the forward flow map @xmath49 by @xmath106 with the terminal condition @xmath104 , and the functions @xmath107 and @xmath108 are two intermediate predicted solutions at the time levels @xmath109 and @xmath90 , respectively . then the flow map from the time level @xmath96 to @xmath95 is given by @xmath110 . it is rather natural to decompose the flow map @xmath111 into a composition of maps @xmath112 . for example , such idea has been used recently in @xcite to improve the computational efficiency of the lagrangian ftle construction between various time levels . define the interpolation operator by @xmath113 , which returns the interpolated mesh values from the non - mesh values . then the flow map @xmath48 can be decomposed into @xmath114 . if the maps @xmath115 for @xmath116 are all stored once they are computed , one can form @xmath48 by determining only @xmath117 . the idea in this work shares some similarities with what we are implementing . our approach can also re - use all intermediate flow maps , the work @xcite however has concentrated only on improving the computational efficiency of the lagrangian approach . we are using the idea of flow map decomposition to propose a forward computational strategy for the eulerian formulation .    to obtain a stable evolution in the flow map constructions , theoretically we require the interpolation scheme between two local flow maps to be monotone , i.e. the interpolation scheme should preserve the monotonicity of the given data points @xcite . this is especially important in the context of the phase flow maps for autonomous flows @xcite since the overshooting due to the interpolation could be significantly amplified in the time doubling strategy @xmath118 ^ 2 \ , .\ ] ]      in this section , we propose an efficient eulerian approach to compute the isle function for an arbitrary separation factor @xmath0 based on the techniques developed in the previous section .    to compute the isle at @xmath119 , one has to determine the minimum time @xmath120 at each location @xmath119 for which @xmath121 . it is , therefore , necessary to keep track of all intermediate values @xmath122 for all @xmath123 s during the flow map construction . at each grid point @xmath124 and each intermediate time step @xmath90 , we construct a quantity @xmath125 by @xmath126\ ] ] with @xmath127 . once we determine the flow map @xmath128 , we have also constructed an increasing sequence at each grid point @xmath124 in time @xmath129 , i.e. @xmath130 for each @xmath124 . now , we interpret this quantity as a spatial - temporal volumetric data in the @xmath1 space . due to such monotonicity property in the temporal direction , the isosurface @xmath131 forms a graph of @xmath119 , i.e. for each @xmath119 we have a unique @xmath132 so that @xmath133 . and more importantly , the value @xmath134 gives the shortest time so that @xmath135 . this implies that @xmath136 and so the level surface in the extended phase space , i.e. @xmath137 can be used to define @xmath138 . finally , the isle function can then be computed using @xmath139 the computational algorithm is summarized in algorithm 2 . ( 4,0 ) ( -0.5,0)(1,0)12     +    algorithm 2 ( our proposed eulerian method for computing the isle ) :    1 . discretize the computational domain @xmath140 2 . set @xmath141 . 3 .   for @xmath94 , 1 .   solve the level set equation ( [ eqn : levelset ] ) for @xmath142 with the condition @xmath143 . 2 .   define @xmath144 . 3 .   interpolate the solutions to obtain @xmath145 . 4 .   compute @xmath146 and determine @xmath147 \ , .\ ] ] 4 . given an arbitrary separation factor @xmath0 . for each @xmath124 , search for a time step @xmath90 such that @xmath148 . then determine @xmath138 by applying the piecewise linear interpolation using the values @xmath125 and @xmath149 . + the isle function is computed as @xmath150    ( 4,0 ) ( -0.5,0)(1,0)12     + there are several major advantages of the proposed eulerian approach . it provides an efficient yet simple numerical implementation of @xmath151 for multiple separation factors @xmath0 . the construction of the quantity @xmath152 can be first done in the background . for each particular given value of @xmath0 , one has to perform only one single isosurface extraction in the extended phase space , i.e. the @xmath1 space . unlike the typical lagrangian implementation where one has to repeatedly keep track of _ all _ particle trajectories for the constraint @xmath153 , the proposed eulerian algorithm provides a systematic way based on a simple thresholding strategy . note that in our implementation , if it happens that @xmath154 for all @xmath90 , we simply set @xmath155 . we conclude this section by discussing the computational complexity of our proposed algorithm for a two dimensional flow , i.e. @xmath156 . let @xmath157 and @xmath54 be the discretization size of one spatial dimension and time dimension respectively . since the liouville equation is a hyperbolic equation , we have @xmath158 by the cfl condition . at each time step @xmath90 , a short time flow map @xmath49 is obtained by solving the liouville equations from @xmath159 to @xmath160 , the computational effort is @xmath161 . then in computing the long time flow map @xmath162 , an interpolation @xmath163 takes @xmath161 operations . therefore the construction of flow map at each time step requires @xmath164 operations . computing @xmath165 and @xmath149 also needs @xmath161 operations , the complexity order is kept at @xmath161 . summing up this procedure in all time steps , the overall computational complexity is @xmath166 . this new eulerian method has significantly improved the overall computational complexity in the application to the isle computations . since the calculations involve the forward flow map from @xmath167 to @xmath168 for all @xmath169 , the original eulerian approach as discussed in section [ subsec : oldeulerian ] requires to solve the liouville equation from each individual time level @xmath168 backward in time to @xmath167 . the overall computational complexity in obtaining all flow maps is therefore given by @xmath170 which is one order magnitude larger than the newly developed eulerian method . with algorithms proposed in section [ sec : proposal ] , we can finally develop an efficient numerical approach to compute the _ forward _ ftle field by solving corresponding pdes _ forward _ in time and the isle field for arbitrary separation factor @xmath0 . the ftle and the isle ( or the fsle ) are indeed measuring different properties of a given flow , even though both quantities depend on the same function @xmath171 which relates the growth of an infinitesimal perturbation . in particular , ftle measures how much @xmath171 grows over a finite time span @xmath172 $ ] , while the isle measures the time required for the quantity @xmath171 to reach a threshold value prescribed by the separation factor . they are two different tools to study how rapid @xmath171 grows . a careful comparison of these two quantities can be found in , for example , @xcite .    in many numerical experiences , on the other hand , the ftle fields can show striking visual resemblance with the isle fields of certain separation factors . however , there has not been much studies on the theoretical relationship between these two quantities . in this section , we are going to demonstrate how these two quantities relate to each other by considering their corresponding _ ridges_. of course , the way to define a _ ridge _ is not unique and one can pick a convenient definition for a particular application @xcite . in this work , we adapt the following simple definition .    for a given scalar function @xmath173 , we define a @xmath174-ridge to be a codimension - one compact and @xmath175 surface , denoted by @xmath176 , such that for every point @xmath177 , @xmath173 is a local maximum along the normal direction to @xmath176 at @xmath119 .    at these locations , @xmath173 is also called a generalized maximum in some literatures . next , we define a particular neighbourhood of a ridge for later discussion . suppose @xmath178 is a codimension - one , compact and smooth surface . we define @xmath179 , the tubular neighbourhood of @xmath178 with radius @xmath180 , as the collection of any point @xmath181 that can be uniquely expressed as @xmath182 , for some point @xmath183 , a real number @xmath184 , and @xmath185 is an unit normal vector of @xmath178 at @xmath119 .    to ensure that any point @xmath181 on @xmath179 can be uniquely expressed as such , one can equivalently require that the boundary of @xmath179 does not intersect itself . for a smooth surface , this can always be guaranteed by choosing a small enough positive number @xmath186 . the definitions of both the ftle and the isle involve @xmath187 . it is therefore natural to link them through this particular quantity . the following simple theorem which identifies the @xmath188-ridges with the @xmath189-ridges , will be the very first step in building up our analysis . [ prop : ftle&lambda ] let @xmath188 and @xmath189 be the ftle function and the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 , respectively . then , a @xmath188-ridge is also a @xmath189-ridge . since @xmath191 is a just a scaling factor and the natural logarithm function is strictly increasing , all generalized maxima of @xmath188 as defined in ( [ eqn : ftle ] ) and @xmath189 coincide .    for convenience , in the time window @xmath172 $ ] , we impose @xmath192 if @xmath193 for all @xmath194 $ ] . if an isle function @xmath195 on an open set @xmath9 of the spatial domain , there is obviously no visible isle feature on @xmath9 , corresponding to this separation factor @xmath0 . the following theorem establishes a fundamental result that for certain choices of @xmath0 , the isle function @xmath196 can be positive on a tubular neighbourhood of a @xmath189-ridge . this explains the non - trivial isle features around the location of a ftle ridge , when one chooses a suitable separation factor @xmath0 . [ prop : funprop1 ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 . suppose @xmath189 is continuous everywhere in the spatial domain , and @xmath176 is a @xmath189-ridge with @xmath197 . then for any positive number @xmath198 , there is a tubular neighbourhood @xmath199 such that the isle function is positive on @xmath199 .    for each @xmath177 , using the definition of @xmath200 , we have @xmath201 . then by the continuity of @xmath189 at @xmath119 , we can find an open ball @xmath202 with center @xmath119 and radius @xmath203 such that @xmath204 for all @xmath205 . collect all these open balls at every point on @xmath176 into @xmath206 , which is clearly an open cover of @xmath176 . by the compactness of @xmath176 , we can find a finite subcover @xmath207 from @xmath206 . take @xmath208 . without loss of generality , assume @xmath186 is small enough such that the tubular neighbourhood @xmath199 is contained in the union of sets in @xmath209 . now for any @xmath210 , it is clear that @xmath211 . because the function @xmath212 is a continuous function of @xmath213 , we apply the intermediate value theorem and conclude that it must have reached @xmath0 at some @xmath214 $ ] . therefore @xmath215 is positive . the next theorem is to identify a @xmath189-ridge with a @xmath196-ridge of a certain isle function @xmath196 . to justify the condition , the points on the location of a @xmath189-ridge , are expected to exhibit highly chaotic behavior . it is reasonable to assume that we can find a tubular neighbourhood of a @xmath189-ridge , where the time monotonicity of @xmath216 holds for @xmath217 $ ] , and that no other points reach the minimum value of @xmath218 , except those on @xmath176 . [ prop : funprop2 ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 . suppose @xmath176 is a @xmath189-ridge with @xmath219 . if there exists a tubular neighbourhood @xmath199 of @xmath176 such that    1 .   for any fixed @xmath210 , the quantity @xmath220 is increasing in @xmath217 $ ] , and 2 . @xmath221 ,    then @xmath176 is also a @xmath222-ridge for the isle function @xmath222 . let @xmath177 , and @xmath223 be a point in the normal direction of @xmath176 at @xmath119 . by the second condition , we have @xmath224 . and since @xmath220 is increasing with @xmath213 , it follows that @xmath225 for @xmath217 $ ] . therefore @xmath226 by our convention . we can conclude that no other points on the normal direction of @xmath176 at @xmath119 has positive value of @xmath222 , except @xmath119 itself . this makes @xmath119 a generalized maximum of @xmath222 . theorem [ prop : funprop2 ] is a restricted result , in the sense that the whole @xmath189-ridge is preserved for a particular @xmath196 function . but it is straightforward to extend it to the following corollary , which looks for a certain portion of a @xmath189-ridge that can be conserved in a larger class of isle functions . [ prop : funcor ] let @xmath189 be the largest eigenvalue of the cauchy - green deformation tensor from time @xmath190 to time @xmath18 . suppose @xmath176 is a @xmath189-ridge with @xmath219 , and @xmath227 is a connected subset of @xmath176 with @xmath228 . for any positive number @xmath229 $ ] , if there exists a tubular neighbourhood @xmath230 of @xmath227 such that    1 . for any fixed @xmath231 , the quantity @xmath220 is increasing in @xmath217 $ ] , and 2 . @xmath232 ,    then @xmath227 is a @xmath196-ridge for the isle function @xmath196 . the proof is parallel to that of theorem [ prop : funprop2 ] , with @xmath176 and @xmath200 being substituted by @xmath227 and @xmath0 , respectively .    through theorem [ prop : funprop1 ] to corollary [ prop : funcor ] , we have obtained a close relationship between the @xmath189-ridges and the isle ridges . by the equivalence of the ftle ridges and the @xmath189-ridges as proved in theorem [ prop : ftle&lambda ] , we can then link the ftle ridges and the isle ridges . in practice , the estimated location of the isle ridges can be inferred by the location of the ftle ridges , while @xmath189 plays an important role for estimating a suitable separation factor defined in the isle ridge . to demonstrate this , suppose that there is a ftle ridge @xmath176 satisfying all conditions in theorem [ prop : funprop2 ] . then @xmath176 is an isle ridge with the separation factor @xmath233 . although the value of @xmath200 is actually an unknown , we can determine a number @xmath234 that is slightly less than @xmath235 manually . from the definition of the ftle , we have @xmath236 > e^{lt } \ , .\ ] ]    now , with any separation factor @xmath237 , theorem [ prop : funprop1 ] implies that there is a tubular neighborhood around @xmath176 with positive isle values . we can then increase the separation factor from @xmath238 approaching to @xmath200 , so that the tubular neighborhood containing @xmath176 becomes narrower , and thus leads to a better approximation of @xmath176 . in practice , the value of @xmath234 has to be chosen in a trial - by - error fashion and has to be determined case - by - case . based on the eulerian algorithm we developed in section [ subsec : neweulerian ] and section [ subsec : islecomp ] , we are now able to efficiently construct the isle field of each individual separation factor using one single isocontour extraction . in this section , we present numerical results to show the feasibility of the proposed eulerian approach . we will compare the solution from the new approach with the previous eulerian approach discussed in section [ subsec : oldeulerian ] . it is noteworthy that for a dynamical system in a fixed time window , we can now construct the isle function with an arbitrary choice of the separation factor in only one single numerical flow simulation . therefore , we can easily verify the theoretical connection between ftle and isle , developed in section [ sec : relation ] . this first example is a simple flow taken from @xcite to describe a periodically varying double - gyre . the flow is modeled by the following stream - function @xmath239 \sin(\pi y ) \ , , \ ] ] where @xmath240 in this example , we follow @xcite and use @xmath241 , @xmath242 . we discretize the domain @xmath243\times[0,1]$ ] using 513 grid points in the @xmath244-direction and 257 grid points in the @xmath245-direction . this gives @xmath246 .    in figure [ fig : ftle ] , we compare the ftle field @xmath247 of the double - gyre flow with @xmath248 computed using the lagrangian approach and our proposed eulerian approach , respectively . these two solutions match extremely well . we have also checked the @xmath249-error of the forward flow map @xmath250 using different @xmath251 s ranging from @xmath252 to @xmath253 . since we do not have the exact solution for this flow , the _ exact _ solutions are computed using the lagrangian ray tracing with a very small time step . figure [ fig : convergence ] shows the errors in the flow map @xmath254 with respect to different @xmath251 s . we find that the flow map computed using the proposed eulerian approach is approximately second order accurate . figures [ fig : isle1 ] shows the isle fields @xmath255 with several separation factors @xmath0 . in the implementations , the flow map is computed from @xmath2 up to @xmath256 . recall that we impose zero isle value at points whose local separation rate is less than the given separation factor . with the increasing separation factor , the regions with non - zero isle values , get smaller . nevertheless , the regions with non - zero or high isle values , always concentrate near the location of the ftle ridges shown in figure [ fig : ftle ] . we denote by @xmath176 the most prominent ridge originating from the central bottom of the domain . note that a rough estimate of the minimum ftle value on @xmath176 is approximately 0.3 . and so we have the estimate @xmath257 . using theorem [ prop : funprop1 ] , we can always find a tubular neighborhood near @xmath176 of non - zero isle using a separation factor @xmath0 smaller than approximately 20 . in figure [ fig : isle1 ] , we have computed the isle using various separation factors from @xmath258 to @xmath259 . as @xmath0 increases , we can see that there are indeed non - zero isle regions around the major ridge @xmath176 . and evidently , the isle ridges lie within these regions .    to better compare the ftle field and the isle field , we plot these two solutions together in figure [ fig : isosurface ] . the @xmath244- and the @xmath245- axis represent the computational domain @xmath9 while the @xmath260 axis denotes the temporal direction . in ( a ) , we plot the ftle field @xmath261 at @xmath262 and @xmath263 at @xmath264 . in ( b - d ) , we plot the isosurfaces of the isle field @xmath265 with @xmath266 , respectively . the @xmath260-values denote the time required for the local separation ratio to achieve @xmath0 for the first time , i.e. the @xmath260-values are the @xmath267 we defined before . near the major ftle ridge , we can see the values of @xmath120 are significantly lower than the remaining part of the domain . therefore the corresponding isle values are generalized maxima . this shows us the isle ridges match with the ftle ridges very well . we consider a simple analytical field @xcite given by @xmath268 and @xmath269 in this example . the computational domain is @xmath270 ^ 2 $ ] and we discretize the domain with @xmath271 . we impose the fixed inflow boundary condition and non - reflective outflow boundary condition as in @xcite . figure [ fig : ftle2 ] shows the ftle field @xmath272 computed using our original eulerian approach and the proposed eulerian approach . it is clear that the new approach can successfully handle the inflow and outflow of particles on the boundaries , just like the previous approach does . next , we are going to apply the results in section [ sec : relation ] to analyze the isle ridges . figure [ fig : ftle2 ] suggests that there is only one ftle ridge . and the minimum value of @xmath273 on the ridge is close to @xmath274 . by the definition of ftle , we have @xmath275 \ , .\ ] ] therefore we can conclude that @xmath276 is a reasonable estimate of the minimum value of @xmath277 on the ridge . we have tested with several separation factors that are near @xmath278 , and the results can be plotted in figure [ fig : isle2a ] . to better visualize the isle field , we plot in figure [ fig : isle2 ] the time @xmath267 required for different separation factors @xmath0 . the shorter the time , the bigger the corresponding isle .      [ ! htb ]    we consider the dynamical system governed by a duffing and a van der pol oscillator , as in @xcite . the system is given by @xmath279    the computational domain is @xmath280\times [ -1.5,1.5]$ ] , with mesh size @xmath281 . we simulate the flow from @xmath2 to @xmath256 . the discretization size is rather coarse for such a complicated dynamics , but our new approach can still demonstrate the fine details of the transport barrier as shown in figure [ fig : ftle3 ] . there are several major ftle ridges that are easily spotted . the minimum value of @xmath282 on those ftle ridges are very close to @xmath283 . so @xmath284 . theorem [ prop : funprop1 ] implies that the separation factor @xmath285 can approximately single out the isle ridges . we pick several separation factors that are close to @xmath286 , and we find that the isle ridges corresponding to @xmath287 , as shown in figure [ fig : isle3](c ) , have the greatest resemblance to the ftle ridges in figure [ fig : ftle3 ] .      [ ! htb ]    to further demonstrate the effectiveness and robustness of our proposed approach , we test our algorithm on a real life dataset obtained from the ocean surface current analyses real - time ( oscar ) . the oscar data were obtained from jpl physical oceanography daac developed by esr . it measures ocean surface flow in the region from @xmath288 to @xmath289 in latitude and from @xmath290 to @xmath291 longitude . the time resolution is about 5 days and the spatial resolution is @xmath292 in each direction .     in this numerical example , we consider the area near the line islands in the region from @xmath293 to @xmath294 in latitude and from @xmath295 to @xmath296 longitude . to have a better visualization of both the ftle and isle field , we interpolate the velocity data to obtain a finer resolution of @xmath297 days in the temporal direction and @xmath298 in each spatial direction . then we look at the ocean surface current within the first 50 days in the year 2014 to obtain the forward ftle field @xmath299 and the isle field @xmath300 . as demonstrated in figure [ fig : oscar_ftle ] , our new approach works well in extracting fine details of the transport barrier even for real data . there are a few major ftle ridges spotted in the figure . according to theorem [ prop : funprop1 ] , we found that the separation factor @xmath0 roughly equals to @xmath301 to @xmath302 can single out the isle ridges . we have tested several separation factors that are close to @xmath303 , as shown in figure [ fig : oscar_isle ] . we find that the isle ridges corresponding to @xmath259 , in figure [ fig : oscar_isle](b ) , have the greatest resemblance to the ftle ridges in figure [ fig : oscar_ftle ] . as a comparison , we have also plotted the isle field @xmath300 for @xmath304 and @xmath305 in figure [ fig : oscar_isle ] ( a ) and ( c ) , respectively . these results verify the theory developed in section [ sec : relation ] . to summarize , we emphasize the novelty and importance of both the two eulerian algorithms we developed in this paper . based on the first algorithm , we are now able to determine the _ forward _ flow map _ on the fly _ so that the pde is solved _ forward _ in time . unlike the original eulerian algorithm as developed in @xcite , there is no need to store the whole velocity field at all time steps at the beginning of the computations . this makes the forward flow map computations practical . more importantly , because of the first algorithm , we are now able to access all intermediate forward flow maps which is necessary in computing the isle field computations . this eulerian approach is extremely efficient in computing the isle field . to the best of our knowledge , this is the first efficient numerical approach for the isle computations . theoretically , this paper provides a rigorous mathematical analysis to explain the similarity between the ftle ridges and the isle ridges based on @xcite , which has been verified by the proposed eulerian approach . the work of you was supported by the talents introduction project of nanjing audit university , and the natural science foundation of jiangsu higher education institutions of china ( no.16kjb110012 ) . the work of leung was supported in part by the hong kong rgc grants 16303114 and 16309316 .                                                                          c.  w. shu . essentially non - oscillatory and weighted essentially non - oscillatory schemes for hyperbolic conservation laws . in b.  cockburn , c.  johnson , c.w . shu , and e.  tadmor , editors , _ advanced numerical approximation of nonlinear hyperbolic equations _ , volume 1697 , pages 325432 . springer , 1998 . lecture notes in mathematics . w.  tang , p.w . chan , and g.  haller . accurate extraction of lagrangian coherent structures over finite domains with application to flight data analysis over hong kong international airport . , 20:017502 , 2010 .
Abstract:

Answer: we propose a new eulerian numerical approach for constructing the forward flow maps in continuous dynamical systems . the new algorithm improves the original formulation developed in @xcite so that the associated partial differential equations ( pdes ) are solved forward in time and , therefore , the _ forward _ flow map can now be determined _ on the fly_. due to the simplicity in the implementations , we are now able to efficiently compute the unstable coherent structures in the flow based on quantities like the finite time lyapunov exponent ( ftle ) , the finite size lyapunov exponent ( fsle ) and also a related infinitesimal size lyapunov exponent ( isle ) . when applied to the isle computations , the eulerian method is in particularly computational efficient . for each separation factor @xmath0 in the definition of the isle , typical lagrangian methods require to shoot and monitor an individual set of ray trajectories . if the scale factor in the definition is changed , these methods have to restart the whole computations all over again . the proposed eulerian method , however , requires to extract only an isosurface of a volumetric data for an individual value of @xmath0 which can be easily done using any well - developed efficient interpolation method or simply an isosurface extraction algorithm . moreover , we provide a theoretical link between the ftle and the isle fields which explains the similarity in these solutions observed in various applications .