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: the existence of some uncertainties both in the experimental and in the theoretical information concerning the form factors in the amplitudes of various pion decays [ 1 ] stimulated the development of a plethora of models for these processes . the most popular are the quark - loop models , whose central problem is the definition of the quark - meson coupling . in some of these models the quark - meson coupling is introduced via effective lagrangians and pcac [ 1,2,3 ] . these models have a hybrid character due to the simultaneous presence of quarks and mesons . more consistent than these models are those based on the bethe salpeter equation [ 4 ] , which assume that a meson is a bound state of a quark and an antiquark ; however , they have some defficiences related to the form of the interaction kernel , the validity of the instantaneous approximation and the appearance of the relative time variable .    recently , without defining any special binding mechanism , simple models of mesons , like products of free quark and antiquark states with a suitable relativistic coupling of the spins and a particular distribution of the relative momentum have been proposed [ 5 ] . horbatsch and koniuk went a step further by adding a gluon to the @xmath4 state [ 6 ] in order to incorporate some confining effects and possible gluon exchange . some time ago we also proposed [ 7 ] a form of this kind , but , instead of the gluon , we introduced a scalar vacuum - like field @xmath5 intended to represent globally the fluctuations of the quark - gluonic field which can not be described in terms of a few elementary excitations . apppearing under different names like condensates in the shifman , vainshtein and zakharov approach [ 8 ] , or bag in the the bag models [ 9 ] , nonelementary excitations have been formely introduced in many quark models . in deep inelastic processes [ 10 ] they got the name of sea , because they looks like a sea    of highly virtual particles which becomes a real source of new particles at high energy . the main point is that it is almost always necessary to consider besides the valence quarks some nonperturbative effects which are essential in the description of the binding and of the generation of new particles .    in the following section we present the particular features of the model and derive the expressions of the pion decay amplitudes . the comparison of our results with the experimental data is performed in section 3 . in section 4 our conclusions are presented . the appendix contains some useful notations and integrals used throughout the paper , as well as the expression of the partial width of the weak radiative decay and the inadvertencies found by us here . the model used for calculating the decay amplitudes is based on the following picture we have now for hadrons . hadrons look like bags , cavities or bubbles containing a few bound valence quarks in equilibrium with the excitation of the surrounding quark - gluonic field . the equilibrium is stable and the internal distribution of momenta does not change as long as the hadrons remain free . during the interaction , the equilibrium breaks down , the fluctuations of the quark - gluonic field become very intense and new states can emerge . a long experience proves that perturbation theory is unable to describe both the binding effects and the hadron interaction at low energy . this fact shows clearly that the fluctuations of the quark - gluonic field generating these effects are far from being elementary . for this reason we avoid considering some elementary excitations like a small number of quarks or gluons besides the valence quarks and introduce a phenomenological field @xmath5 to describe globally the fluctuations of the quark - gluonic field inside the hadrons . the field @xmath5 has the quantum numbers of the vacuum and carries its own momentum which is not subject of any mass shell constraint because @xmath5 does not represent an elementary excitation . also , in order to bypass the perturbative treatment of the interaction , we do not start from the elementary quark gluon coupling but we prefer looking to the time translation operators , @xmath6 and @xmath7 , which are more appropiate in a phenomenological , nonperturbative approach . we describe by @xmath6 and @xmath7 the modifications produced in the structure of the free hadronic states by an external interaction .    in the cases studied in this paper , the pion is the single hadron involved and nothing can change its structure up to the emission of a weak or electromagnetic quanta . therefore , the time translation operator @xmath7 , where @xmath8 is the momentum of the first emission , acting on the single pion state leaves it unchanged since there was no external interaction to perturbe the internal equilibrium .    in the case of the radiative decays , the emission of the second quanta takes place in a state already perturbed by the recoil of the first emission . due to the smalness of the involved momenta , we assume however , that the perturbation of the equilibrium state can be neglected . therefore , the time translation operator @xmath9 , describing the effect of the quantum fluctuations in the time interval between the two emissions will also be approximated by unity .    according to the ideas outlined above , we propose the following expression for the @xmath0-meson state :    @xmath10 @xmath11 where @xmath12 and @xmath13 are the creation operators of the valence quark and antiquark whose colour , spin and flavour are denoted by the collective indices @xmath14 and @xmath15 . the quark creation and annihilation operators satisfy the free field canonical anticommutation relations , the only nonvanishing being @xmath16_+= \left [ b^+_\alpha({\bf p } ) , b_\beta({\bf q})\right]_+= ( 2\pi)^3 { e\over m}\delta^{(3)}({\bf p } - { \bf q } ) \delta_{\alpha\beta}. \eqno(2)\ ] ] @xmath17 denotes the creation of a nonelementary excitation having the internal quantum numbers of the vacuum and carrying the momentum @xmath18 . it represents the quantum fluctuations of the quark gluonic - field in equilibrium with the valence quarks.it is an independent component of the pion wave function and hence it is supposed to commute with the quark operators . the physical meaning of the phenomenological field @xmath19 allows us to assume that @xmath20    in the same time , keeping into account that in any process the conservation of the quark momentum is ensured by the commutation relations ( 2 ) and that @xmath21 is an independent component of the hadrons , one must impose from outside the separate conservation of its 4-momentum . we assume then @xmath22 where @xmath23 is a constant having the dimension of mass introduced for dimensional reasons . the binding effects in the expression ( 1 ) of the pion state are introduced by means of the internal wave function @xmath24 which must cut off the large relative momenta in order to allow the fulfilment of the orthogonality relation @xmath25 the factor @xmath26 in equation ( 5 ) is due to the equation ( 4 ) following from the nonelementary character of the phenomenological field @xmath5.it forces us to modify the phase space density as it follows @xmath27 where @xmath28 is the probability distribution of the pion mass around a central value @xmath29 . if @xmath28 has a small width , like in our case , the integral over m can be performed from the beginning and one can return to the old phase space density . we pass now to the calculation of the pion decay amplitudes . @xmath30    in the lowest order of the weak interaction , the @xmath31 matrix element of pion decay into leptons is @xmath32 where the weak hamiltonian @xmath33 @xmath34 is written in the limit of local quark - lepton interaction.in equation ( 8) @xmath35 is the fermi constant , @xmath36 is the cabibbo angle , @xmath37 , @xmath38 are the up and down quark fields and @xmath39 , @xmath40 are the lepton fields . the time translation operators @xmath41 , @xmath42 in equation @xmath43 represent the global effect of the quantum fluctuations due to some strong external interactions which are absent both in the single pion and in the final state where only leptons are present.then , as we have already mentioned above , it follows that the time translation operators @xmath41 , and @xmath44 in equation @xmath43 can be replaced by unity .    after introducing in equation ( 7 ) the plane wave decomposition of the free fields ( a1,a2 ) , the expression ( 7 ) of the @xmath31 matrix element can be written as it follows @xmath45 @xmath46 where @xmath47 and @xmath48 are the dirac spinors for the leptons . the matrix element of the quark current in equation ( 9 ) is calculated using the expression ( 1 ) of the pion state.one gets easily @xmath49 @xmath50 @xmath51 @xmath52 @xmath53 where @xmath54 .    in equation ( 10 ) , the condition ( 4 ) has been used.in this case it means that the annihilation of the valence quark and antiquark of the pion into leptons proceeds in a `` bare '' state , i.e. in the absence of other excitations of the quark - gluonic field described by the phenomenological field @xmath5 . using the definition of the pion decay constant @xmath55 , ( a6 ) , one gets from equation ( 10 ) @xmath56 where @xmath57 , and @xmath58~[1-{(m_u - m_d)^2\over m^2}]}$ ] . the factor 3 in equation ( 11 ) is due to the quantum number of colour and is present whenever one encounters a quark loop . it is worth noticing that equation ( 10 ) looks like those appearing in other quark loop models , with the difference that the spin projectors are replaced there by propagators . in our model the valence quark and antiquark are real , on mass shell particles and , because of this fact , the integral over their momenta in equation ( 10 ) leads to the finite result ( 11 ) . this is an essential difference between the present model and other quark loop models [ 2 ] , where @xmath59 is infinite and hence can not be calculated directly in the absence of a consistent renormalization procedure . proceeding in the same manner , we calculate the neutral pion constant @xmath60 defined by ( a5 ) and obtain @xmath61 where @xmath62    obviously , the values of @xmath63 and @xmath60 coincide in the limit of equal quark masses . we associate the diagram in fig . 1 to the expression of the matrix element ( 10 ) . it is important to notice that in our diagram the quark lines correspond to projectors on the states of positive and negative energy , not to propagators like in most quark - loop models [ 2 ] .    _ @xmath64    in the lowest order of perturbation with respect to the weak and electromagnetic interactions , the element of the @xmath31 matrix corresponding to the decay of the positive charged pion can be written @xmath65 @xmath66 @xmath67 where the electromagnetic hamiltonian is @xmath68 @xmath69~a^\lambda(x)\eqno(14)\ ] ] and @xmath70 denote the quark charges . then , just like in @xmath71 case , the time translation operators @xmath72 and @xmath73 can be replaced by unity.we also replace by unity the operators @xmath74 and @xmath75 because the perturbation produced in the equilibrium state by the emission of a quanta can be neglected , as outlined in the introduction , due to the smalness of the involved momenta . we shall comment more about this approximation in the last section . the calculation of the @xmath31 matrix element proceeds then as usually , leading to the following expression : @xmath76 @xmath77 @xmath78g^{\mu\nu } \eqno(15)\ ] ] where the first term is called the bremsstrahlung term and the second one -the structure dependent term.the first one describes the emission of the photon by the charged lepton , while the second one describes the emission of the photon by the meson itself .    after using wick s theorem to pass from time ordered to normal ordered products of quark operators , we can write the matrix elements involving two currents in equation ( 15 ) in the following manner @xmath79 @xmath80 @xmath81 @xmath82 @xmath83 @xmath84 @xmath85 @xmath86 @xmath87\ .\right . \eqno(17)\ ] ] the bremsstrahlung and the structure dependent terms can be represented diagramatically like in figures ( 2a ) and ( 2b ) .    after performing the integrals over the spatial coordinates and over the internal momenta ( a10,a11,a15,a20 ) in equations ( 15 ) , ( 16 ) and ( 17 ) , we can put the the amplitude of the weak radiative decay into its standard form [ 1 ] @xmath88 @xmath89 ( 1+\gamma_5)v(l)\right.\ ] ] @xmath90 \epsilon^\lambda\ \bar u_{\nu_{l}}(l')\gamma^{\mu } \left(1-\gamma_5\right)v(l)\right\}\eqno(18)\ ] ] where @xmath91[1-(m_u+m_d)^2/m^2]}$ ] , @xmath92 , @xmath93 , @xmath94 is given by eq.(11 ) and the form factors @xmath95 and @xmath96 express in terms of the usual adimensional variable @xmath97 and of the quark masses as it follows @xmath98 @xmath99\left[{2\over3 } { m_u\over m}\ln{e+p\over e - p}-{1\over 3 } { m_d\over m}\ln{\varepsilon+p\over \varepsilon - p}+{2 p(m_d - m_u)\over m^2}\right ] ; \eqno(19b)\ ] ] @xmath100 @xmath101\left [ { 2p\ ( m_u+m_d)\over m^2}-{2p\ ( m_d - m_u)\over 3m^2}\right.\ ] ] @xmath102\eqno(20b)\ ] ] @xmath103\ ( m_d - m_u ) \left[{p\over 3m^2}-{p(m_d^2-m_u^2)\over m^4}\right.\ ] ] @xmath104.\eqno(20c)\ ] ] it is worthwile mentioning here that treating in a unitary manner the bremsstrahlung and the structure dependent terms simply by using the weak and electromagnetic hamiltonians of the standard model , one obtains automatically the manifest gauge invariant form ( 18 ) of the radiative decay amplitude . this result is a direct consequence of the consistent treatment of the electromagnetic interaction . it is also important to notice that , contrary to the current opinion , our form factors ( 19 ) and ( 20 ) depend essentially on the photon energy . both of them have poles at @xmath105 , i.e. at vanishing photon energy , which lead to an infrared divergence in the axial contribution to the decay width . although this fact is not surprising since it appears in the bremsstrahlung part too , it requires a careful analysis of the soft photon limit , which is beyond our purpose now . @xmath106    the calculation of the neutral pion decay amplitude can be treated in the same manner we did in the case of weak radiative decay of the charged pion . then , just like in that case , one has @xmath107 @xmath108 @xmath109 \right\vert\pi^{(0)}(p)\right\rangle\ ] ] @xmath110 @xmath111 + \left[{x\to y}\atop { \mu\to\nu}\right]\right\vert\pi^0(p)\right\rangle\ ] ] @xmath112    @xmath113 \eqno(21)\ ] ] where @xmath114 and @xmath115 . writing the @xmath31 matrix element of the neutral pion decay into its usual form @xmath116 and performing the trace and the integration over the spatial coordinates and internal momenta , one gets from ( 21 ) and ( 22 ) @xmath117 \eqno(23)\ ] ] where @xmath118 , and @xmath119 . as it can be seen from the relations ( 11 ) , ( 12 ) , ( 19 ) , ( 20 ) and ( 23 ) , the expressions obtained for @xmath55,@xmath60,@xmath96,@xmath95 and @xmath120 depend on the quark masses and on the value of the internal wave function in a given point @xmath121.in order to reduce as much as possible the dependence on the particular form of the internal function , we shall consider the ratios @xmath122,@xmath123,@xmath124 , @xmath125 which are functions of the quark masses only , to be compared with the experimental values . we recall that , contrary to most quark loop models , in our case the valence quarks are real , on mass shell particles and their momenta must fulfill the condition @xmath126 , which is impossible to be satisfied by the constituent quarks . we try therefore to fit the experimental data with current quark masses and analyse our results in the range of masses @xmath127 , @xmath128 as quoted in particle data [ 10 ] .    first we consider the function    @xmath129 and plot its contour lines in this domain ( fig.3a ) . the plot shows a very quick variation of @xmath130 over the mass range . it also shows that there is a continuous sequence of masses with @xmath131 leading to a perfect agreement between our result and the experimental one .    as concerns the form factors occuring in the weak radiative decay , we notice that one can not perform a direct comparison with the experimental values because , as pointed above , in our case the form factors @xmath95 and @xmath96 depend significantly on the photon energy , while the experimental data [ 11,12 ] have been fitted with some constant values . however , in order to test the predictive power of the model , we treat these constant values like mean values of the form factors over the space region where the measurements have been made . accordingly , we have to compare @xmath132 with @xmath133 where @xmath134 is the mean value of @xmath135 over the phase space region . since there are two different kinematical weights @xmath136 and @xmath137 given by the equations ( a26 ) and ( a27 ) which correspond to the left and right photons , one can define two mean values @xmath138 and @xmath139 : @xmath140 where @xmath141 [ 1 ] , @xmath142 and @xmath143 are given in the appendix and @xmath144 is the volume in phase space . following v. n. bolotov et al . [ 12 ] , we choose @xmath145 with @xmath146 and find that @xmath147 . in fig.3b we ploted the contur lines of the function @xmath148    @xmath149 . just like in @xmath130 case one observes a very quick variation over the quark mass range . also , for @xmath150 one finds @xmath151 which represents a quite good agreement with the values quoted for @xmath95 [ 11 ] . we try the same procedure for the axial form factor but in this case the comparison is not so relevant because of the existence of two contributions to the axial form factor : one with a simple pole and the other with a double pole at @xmath105 . defining as above @xmath152 and @xmath153 we find @xmath154 for the volume @xmath144 . as it can be seen from figs.3c and 3d , the values of @xmath155 and @xmath156 are rather large for the quark masses in the above range , where the agreement is reasonable in the previous cases . it follows then that it is very hard to reach a small mean value for @xmath157 as quoted in particle data [ 11 ] , unless there is an almost perfect compensation of the two contributions which are of different signs . the compensation would roughly occur for @xmath158 , but , as it can be seen from fig.3e , the function    @xmath159 does not reach this value in the above range of current quark masses . then the difference between our prediction and the experimental value of @xmath96 can not be reduced under an order of magnitude , which is quite bad . the model presented in this paper has some attractive features when compared with other quark models for pion decays . it is essentially relativistic and allows the consistent treatment of the pion structure and of the electromagnetic interaction , leading to manifest gauge invariant amplitudes for radiative decays . it also leads to finite expressions for the pion decay constants @xmath55 and @xmath60 which could in principle be compared with the experimental values . however , taking into account their dependence on the internal wave function which could in principle be different in the two cases , a comparison between @xmath55 and @xmath60 is not relevant for the predictive power of the model . the dependence of the results on the particular internal function of the pion is significantly reduced by considering the ratios @xmath160 which are functions on the quark masses only . the first of these ratios is in agreement with the experimental value for current quark masses in the range quoted in particle data [ 11 ] . as concerns the other two ratios , a significant comparison would require some information about the eventual dependence of the form factors on the photon energy . in this sense , a new fit of the measurements with @xmath161 dependent form factors including also the interference term which was neglected up to now , would be most interesting and could possibly remove the uncertainty in the quoted value of @xmath96 [ 12 ] .    in the preceding section we pointed out a very quick variation of @xmath162 and of @xmath163 in the region of masses selected for the best fit due to the presence in the expressions ( 19 ) and ( 20 ) of logarithms of nearly zero arguments . it must be however stressed that this fit can not be used for the precise determination of current quark masses without a quantitative analysis of the recoil effects on the internal equilibrium in the quark system after the first emission of a quanta . the last comment concerns the infrared behaviour of the form factors . as it can be seen from the appendix , only the axial contribution to the structure dependent part leads to an infrared catastrophe , i.e. to a divergence in the expression of the decay width . according to the standard procedure used in quantum electrodynamics , the infinities of this kind in the cross section are eliminated by compensation with those appearing in the interference term between the amplitude with radiative corrections and without them [ 14,15 ] . in our case this compensation mechanism does not work because of the particular form of the decay amplitude . the problem deserves a further study following closely the arguments of bloch and nordsieck [ 16 ] which are more appropriate in this case . appendix          using the criterion of covariance we write the integrals @xmath179 as it follows @xmath180 we observe that @xmath181 @xmath182 and obtain @xmath183 in a similar way , using covariance , one has @xmath184 and , proceeding as above , one gets for @xmath185 @xmath186 @xmath187 @xmath188 we solve the above system and write the integrals @xmath189 as it follows @xmath190 @xmath191 \ { k_\mu\over\omega}{k_\nu\over\omega}.\eqno(a20)\ ] ] _ c _ ) . @xmath192 decay width    following the classic paper of brown and bludman [ 13 ] and ref.1 , the amplitude of the @xmath192 decay writes as follows @xmath193 where @xmath194 @xmath195 @xmath196 , \eqno(a23)\ ] ] whith @xmath197 and @xmath198 $ ] .    taking the square of the modulus of the above amplitude and summing over the polarisations one gets the expression of the partial decay rate [ 11,13 ] @xmath199 where @xmath200 @xmath201\ ] ] @xmath202 @xmath203.\eqno(a24)\ ] ] here @xmath204\eqno(a25)\ ] ] @xmath205\eqno(a26)\ ] ] @xmath206\eqno(a27)\ ] ] @xmath207\eqno(a28)\ ] ] @xmath208 g  =  1-y+rx(x+y-1-r)(a29)@xmath209 where @xmath210 , @xmath211 and @xmath212 . the first one concerns the sign of the last term in the square bracket in @xmath213 , which is different from that appearing in ref [ 13 ] and the forthcoming papers . this difference is negligible in the case of pion decay into light leptons , but it can be important in other decays involving heavier leptons . the second one concerns the interchange of @xmath96 with @xmath214when passing from @xmath215 decay ( ref . [ 11 ] ) to @xmath216 decay . ref.[13 ] does not mention this difference . a. bryman , p. depommier , c. leroy : phys . rep . 88(1982)151 2 m . moreno , j. pestieau : phys . rev . * d13*(1976)175 ; n. paver , m. d. scadron : nuovo cim . * 78a*(1983)159 ; l.i.ametler , c. y. ayala , a. bramon : phys . rev . * d29*(1984)916 ; ching - yieh lee : phys . rev . * d32 * ( 1985)658 3s . l. adler : phys . rev . * 177*(1969)2426 ; j. s. bell , r. jackiw : nuovo cim . * 60a*(1969(46 4d . gromes : z. phys . c * 57*(1993)631 ; i. guiasu , r. koniuk : phys . lett . * b314*(1993)408 5p . colangelo , g. nardulli , m. pietroni : phys . rev . * d43*(1991)3003 ; p. l. chung , f. coester : phys . rev . * d44*(1991)229 ; p. j. odonnell , h. k. k. tung : phys . rev . * d44*(1991)741 6 m . horbatsch , r. koniuk : phys . rev . * d47*(1993)210 7l . micu : czechoslovak j. phys . * 41*(1991)129 8 m . a. shifman , a. i. vainshtein , v. i. zakharov : nucl * b147*(1979)385 , 448 , 519 9a . chodos , r. l. jaffe , k. johnson , c. b. thorn , v. f. weisskopf : phys . rev . * d9*(1974)3471 10r.feynman in photon - hadron interactions w. a. benjamin , n. y. 1972 ; j. d. bjorken , e. paschos : phys . 186*(1969)1975 11particle data group : phys . rev . * d45*(1992)vii.1 , vii.4 12r . eichler et al . : . lett . * b175*(1986)101 ; a. bay et al . : phys . lett . * b174*(1986)445 ; v. n. bolotov et al . : * b243 * ( 1990)308 13s . g. brown , s. a. bludman : phys . rev . * 136*(1964)b1160 14j . m. jauch , f. r. rohrlich : helv . acta * 27*(1954)613 15 t . goldman , w. j. wilson : phys . rev . * d15*(1977)709 16f . bloch , a. nordsieck : phys . 52*(1937)54      fig . 1 quark diagram for the pion decay constant @xmath59 . the pion is represented by the quark and the antiquark lines . the bubble represents the pion internal wave function . 2 a ) the diagram of the bremsstrahlung term ;
Abstract:

Answer: the amplitudes of @xmath0-meson decays are calculated in a relativistic nonperturbative quark model which assumes that mesons are made of a quark antiquark pair and of a scalar neutral component representing the contribution of the nonelementary fluctuations of the quark gluonic field . the experimental data can be fitted in a satisfactory manner using current quark masses @xmath1 . -0.6 cm    _ ifa - ft-395 - 1994 _    * *    @xmath2    * and @xmath3 decays *    * *    institute of atomic physics    theoretical ph ysics department    bucharest pobox mg-6 , romania


Instruction: Write an abstract for the following acadamic paper.

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

Paper: in his famous work on statistical hydrodynamics @xcite , onsager has predicted that a system of vortices can have negative temperature states , which can be interpreted as a clustering of vortices of the same sign . more specifically , onsager considered a microcanonical ensemble of parallel point - vortices in an incompressible inviscid fluid confined to a finite region of physical space and calculated the entropy of the system @xmath0 as a function of its energy @xmath1 . taking into account that the coordinates of point vortices obey hamiltonian equations and thus are canonically conjugate , he concluded that the phase volume of the system @xmath2 is also finite , so that the density of states @xmath3 should have a maximum value at some ( critical ) energy @xmath4 . correspondingly , the statistical temperature @xmath5 , where @xmath6 is the entropy , should be positive at @xmath7 and negative at @xmath8 . using these arguments , onsager has concluded that `` the large compound vortices formed in this manner will remain as the only conspicuous features of the motion ; because the weaker vortices , free to roam practically at random , will yield rather erratic and disorganised contributions to the flow '' . this discovery of the negative temperature states has given rise to a number of theoretical considerations , where onsager s arguments were revisited , discussed and elaborated @xcite ( for review , see eyink and sreenivasan @xcite and campa et al . @xcite ) . in particular , it was shown that the key assumption of a finite region of physical space is not crucial because the accessible phase space is restricted by a constant of motion , which is an analog of enstrophy @xcite . recently , we have proposed a statistical model for a fluid flow , which considers the flow as a collection of spatially localized ( coherent ) structures and uses the principle of maximum entropy to determine the most probable size distribution of the structures @xcite . a principal difference of this model from the previous models in which a similar approach was used to describe inviscid flows ( for a review of such models , see , e.g. , refs . or ) is that the structures are assumed to be composed of elementary cells in which the behavior of the particles ( atoms or molecules ) is uncorrelated . the characteristic size of elementary phase space volumes corresponding to these cells is determined by the kinematic viscosity of the fluid , which makes it possible to introduce the reynolds number and thus to extend the statistical description of the flow to viscous fluids . the model successfully describes the transition to turbulence at large reynolds numbers and some other characteristic properties of turbulent flows @xcite . since the dependence of the energy and entropy of the flow on the reynolds number is available in this model , the statistical temperature of the flow state can be calculated as a function the reynolds number . in the present paper we show that at low reynolds numbers , associated with laminar motion , the temperature is positive , and at high reynolds numbers , associated with turbulent motion , it is negative . at intermediate reynolds numbers , representing the transition range , the temperature changes from positive to negative as the structure size exceeds some critical value , similar to what was discovered by onsager for the system of parallel point - vortices in a inviscid fluid . we also show that in the range of intermediate reynolds numbers the temperature exhibits a power - law divergence characteristic of second - order phase transitions .    the paper is organized as follows . section ii briefly describes the statistical model of the flow that was developed in ref . and serves as an introduction to sect . section iii introduces the statistical temperature of the flow states and presents the results of the study and their analysis , including the analogy with second - order phase transitions . section iv contains concluding remarks . the model is briefly described as follows ( for details , see ref . @xcite ) . let a fluid flow be represented by a system of @xmath9 identical particles ( atoms or molecules ) which are placed in a volume of physical space @xmath10 . assume that the particles can form spatially localized structures of @xmath11 particles ( @xmath12 ) , with the number of such ( @xmath13 ) structures being @xmath14 . these structures are viewed as coherent structures , i.e. , the structures in which the constituting particles execute a concerted motion on the structure scale . since such concerted motion should break down at the microscale level , each @xmath13 structure can be divided into @xmath15 elementary `` cells '' of volume @xmath16 ( each of @xmath17 particles ) , in which the behavior of the particles is uncorrelated . because the particles are identical , to each elementary cell there corresponds a 6-dimensional elementary volume of single - particle phase space , in which the state of the particles is uncertain . for a gas fluid , the velocity and linear scales of such a volume are , respectively , molecular thermal velocity @xmath18 and mean free path @xmath19 ( kusukawa @xcite ) . according to the kinetic theory of gases @xmath20 , where @xmath21 is the kinematic viscosity ( see , e.g. , ferziger and kaper @xcite ) , i.e. , the characteristic linear size of the elementary phase volume is determined by the kinematic viscosity @xcite . for a liquid fluid , the linear size of the elementary volume can be taken to be @xmath22 , where @xmath23 is the fluctuation of the ( molecular ) kinetic energy per unit mass , and @xmath24 is the mean residence time of the molecule in a settled state , so that a similar estimate is valid , @xmath25 ( e.g. , frenkel @xcite ) . alternatively , the space and velocity scales of the volume can be associated with the corresponding kolmogorov microscales @xcite , i.e. , @xmath26 and @xmath27 , respectively ( @xmath28 is the rate of dissipation per unit mass ) @xcite . the given system of coherent structures is characterized by two parameters , which are the total number of particles ( atoms or molecules ) in the system @xmath9 and the total energy @xmath1 , i.e. , the corresponding ensemble of the systems is a _ microcanonical _ ensemble . we assume that the total energy can be written as @xmath29 , where @xmath30 is the kinetic energy of the concerted motion of particles in @xmath13 structure . this definition is based on two assumptions : i ) the concerted motion of particles in each local region of the physical space is associated with one of the coherent structures and contributes to its kinetic energy , and ii ) the potential energy of interaction of coherent structures is either negligible in comparison to the kinetic energy of the system or , if not , can be considered as a part of its kinetic energy , as , e.g. , for the system of two - dimensional point vortices in inviscid fluid ( batchelor @xcite ) , which was studied by onsager @xcite . given the above definition of the total energy , the collection of the coherent structures can be treated as an ideal gas system , i.e. , a system in which the potential energy is negligible and the total energy is additive ( landau and lifshitz @xcite ) . it is important that the ideal gas approximation does not imply that the interaction in the system is absent at all . one example , which was used by ruelle @xcite to draw an analogy to the turbulence cascade , is the heat conductivity in a rarefied gas : while the gas is considered to be an ideal gas , the heat transfer exists due to the energy exchange between the molecules in the course of their short - range collisions . this example is also useful in another respect . although the process of the heat transfer is a non - equilibrium process , the molecular velocity distribution function , which corresponds to the the fourier law , is just slightly different from the equilibrium ( maxwell ) distribution , as a first - order perturbation of the latter ( e.g. , ferziger and kaper @xcite ) . therefore , if the statistical properties of the system is of interest , the equilibrium distribution function can be considered as a zero - order approximation . every distribution of particles ( atoms or molecules ) among the coherent structures that satisfies the conditions @xmath31 and @xmath29 presents a _ microstate _ of the system at given @xmath9 and @xmath1 , and the number of such microstates @xmath32 plays a role of the _ density _ of states . taking into account that the particles are identical , so that the permutations of the particles in the elementary cells , the permutations of the elementary cells in the coherent structures , and the permutations of the coherent structures themselves do not lead to new states , we have @xmath33^{m_i}m_{i}!}\ ] ] or , with the stirling approximation @xmath34 to be applicable @xmath35 using eq . ( [ eq1 ] ) , it is possible to calculate the most probable size distribution of the structures , @xmath36 , which maximizes the entropy @xmath37 ( the boltzmann constant is set to unity ) and determines the observational _ macrostate _ of the system . the variation of the entropy functional with respect to @xmath38 at two conservation conditions @xmath39 and @xmath40 yields @xmath41 where and @xmath42 and @xmath43 are the constants depending on @xmath9 and @xmath1 ( the lagrange multipliers ) . to specify the dependence of @xmath30 on @xmath44 , we followed kolmogorov @xcite , i.e. , we assumed that the energy of concerted motion of particles increased with distance as @xmath45 , where @xmath46 . since the coherent structures may be different in form , the law of dependence of @xmath30 on @xmath44 can also be different . to estimate the law , we calculated the kinetic energy in a parallelepiped with sides @xmath47 , @xmath48 and @xmath49 as @xmath50 , where @xmath51 , and related the obtained energy to the parallelepiped volume @xmath52 . the results of the calculations are presented in fig . [ aspect_rations ] . panel * a * shows typical coherent structures that were observed by moisy and jimnez in direct numerical simulations of isotropic turbulence @xcite , who characterized the structures by the aspect ratios @xmath53 and @xmath54 ( solid triangles ) . there are also shown the corresponding structures we chosen to calculate the dependence of the kinetic energy on the structure volume ( crosses ) . panel * b * of fig . [ aspect_rations ] depicts the obtained dependence of @xmath30 on @xmath55 for these structures . as is seen , it is well fitted by the function @xmath56 with @xmath57 , which is close to @xmath58 characteristic of the spherical structures @xcite . the same exponent @xmath59 was used by ruelle @xcite , who , instead of introducing the coherent structures , divided the physical space in cubes in which the kinetic energy of fluctuations was assumed to obey the kolmogorov theory @xcite . in the present paper , the calculations are made with @xmath59 . assuming the density of the fluid to be constant , we thus have @xmath60 .    strictly speaking , the kolmogorov relation assumes the inertial interval of scales @xmath61 , where @xmath62 and @xmath63 are the dissipation and external scales , respectively . however , for rough estimates , the range of its validity can be extended to the lower and upper bounds , i.e. , to @xmath64 @xcite . moreover , with the arguments given by landau and lifshitz @xcite , the lower bound can be associated with laminar motion @xcite . according to this , we associate the laminar state with the structures that are comparable in size with the elementary cell , i.e. , @xmath65 , and the turbulent state with those that contain many elementary cells , i.e. , @xmath66 .    after substituting the dependence @xmath67 on @xmath11 into eq . ( [ eq2 ] ) , we eventually have @xmath68 where @xmath69 is the number of the elementary cells in @xmath13 structure , and @xmath70 and @xmath71 are new constants that should be determined from the equations of conservation of the total number of particles and kinetic energy in the form @xmath72 and @xmath73 it is also possible , and is more convenient , to vary @xmath70 and @xmath71 to obtain @xmath74 and @xmath75 as functions of @xmath70 and @xmath71 @xcite .    the parameter @xmath74 can be associated with the reynolds number @xmath76 , where @xmath77 and @xmath63 are , respectively , the velocity and linear scales characterizing the system as a whole . as has been mentioned , the characteristic size of the elementary phase volume is determined by the kinematic viscosity , so that the number of particles in the volume @xmath78 . correspondingly , the total number of particles @xmath9 can be considered to be proportional to the total phase volume for the system , i.e. , @xmath79 . then @xmath80 , or assuming for simplicity that the coefficient of proportionality is equal to 1 @xmath81    the distribution given by eq . ( [ eq3 ] ) is drastically different from the conventional gibbs - boltzmann distribution for a many - particle system ( see , e.g. , landau and lifshitz @xcite ) . in particular , none of the lagrange multipliers @xmath70 and @xmath71 can be associated with temperature . as has been shown in our previous paper @xcite , these multipliers determine the range of the reynolds numbers where a turbulent state can be expected . to illustrate this , fig . [ q(re ) ] shows the average number of elementary cells in the coherent structures @xmath82 as a function of the reynolds number @xmath83 determined via eqs . ( [ eq3_1 ] ) and ( [ eq4 ] ) . the number of the elementary cells @xmath84 varied from 1 to @xmath85 ; a variation of @xmath86 does not change the overall picture , except that the ranges of variation of @xmath87 and @xmath83 extend to larger values of these quantities as @xmath86 increases @xcite . the values of @xmath70 and @xmath71 were randomly chosen from the uniform distributions within @xmath88 and @xmath89 , respectively ( for @xmath90 samples in each case ) . as was indicated in ref . , not every combination of @xmath70 and @xmath71 gives a physically reasonable bell - shaped distribution @xcite , i.e. , for some combinations , the fraction of which is typically within several percentages , @xmath91 does not vanish at @xmath92 . in fig . [ q(re ) ] , the points with such `` wrong '' @xmath93 combinations are excluded .    according to the present model , the points of fig . [ q(re ) ] for which @xmath94 should be associated with a laminar state , and those for which @xmath95 with a turbulent state . a boundary that seemingly separates these states lies at a reynolds number @xmath96 , which is between @xmath97 and @xmath98 . the positions of the lower and upper boundaries of the manifold of points representing @xmath87 as a function of @xmath83 are determined by the maximum values of @xmath71 and @xmath70 , respectively @xcite . the lower boundary shifts to lower values of @xmath99 as @xmath71 increases , and the upper boundary shifts to larger values of @xmath83 as @xmath70 increases , specifically as @xmath100 and @xmath101 , respectively . however , while the upper boundary shifts unlimitedly , the lower boundary changes its position until the value @xmath102 is reached , after that it `` freezes '' @xcite ( see also fig . [ q(re ) ] ) . it is significant that the dependence of @xmath103 on @xmath83 is not unique , i.e. , the state of the flow is not solely determined by the reynolds number ; rather , as is well - known , it can depend on the type of the flow , the inlet conditions , the flow environment , etc . @xcite . in this respect , it is important that for large values of @xmath70 , the states with @xmath104 are present at the reynolds numbers far above its `` critical '' value @xmath96 ; for example , at @xmath105 the values of @xmath106 are found up to @xmath107 . the model thus predicts that the flow should not only be laminar at @xmath108 but can also remain laminar at @xmath109 , as it was observed , e.g. , for pipe flows by reynolds @xcite and confirmed in many subsequent studies @xcite . we note that the specific properties of the coherent structures , i.e. , the structure shape and the velocity distribution inside the structure , do not seem to be of critical importance here . in particular , the structures in the form of the burgers vortices of finite length ( to mimic `` worms '' ) lead to the dependence of the structure size on the reynolds number that is qualitatively similar to that in fig . [ q(re ) ] @xcite . a highly nontrivial property of the present system , which distinguishes it from the conventional many - particle systems in statistical physics @xcite , is that the entropy of the system is non - additive . for example , if we divide the given system of @xmath9 particles into @xmath110 identical and noninteracting subsystems , the density states of the composite system will be @xmath111^{k}}% { \{\prod_{i}\left[(n!)^{n_{i}/n}(n_{i}/n)!\right]^{l_i}l_{i}!\}^{k}}\ ] ] where @xmath112 . comparison of this equation with eq . ( [ eq0 ] ) yields @xmath113 where @xmath114 and @xmath0 are the entropies of the composite and original systems , respectively . the non - additivity of entropy is consistent with the fact that the system of particles representing a fluid flow can not be divided into parts without a loss of its essential properties . according to eq . ( [ eq4 ] ) , the division of the system into subsystems will reduce the reynolds number characterizing a subsystem ( in the above example , in @xmath115 times ) , so that the resulting reynolds number for the subsystem will not correspond to the flow state in the subsystem , which is assumed to be the same as in the total system . in other words , the condition that the reynolds number is determined by the number of particles in the system [ eq . ( [ eq4 ] ) ] , originates a virtual interaction of the coherent structures that makes unfavorable the separation of the flow into domains at equilibrium . as can be seen from eq . ( [ eq0 ] ) , the entropy of the system is additive only at the level of the coherent structures , i.e. , each collection of the structures of a specific size gives an additive contribution to the system entropy .    according to eq . ( [ eq0 ] ) , the density of states @xmath116 and , correspondingly , the entropy @xmath117 , assume that the particles ( atoms or molecules ) are identical , and an elementary volume of the phase space exists in which the state of the particles is uncertain . the same assumptions play a key role in the derivation of the navier - stokes equation from the the liouville equation , and thus they are implicitly present in the navier - stokes equation @xcite ( for the procedure of the derivation of the navier - stokes equation for a gas fluid from the liouville equation see , e.g. , ferziger and kaper @xcite ) . more specifically , assuming the particles to be identical ( the first assumption of the present model ) , the liouville equation for the many - particle distribution function is reduced to a kinetic equation which contains single - particle and two - particle distribution functions ( the bogoliubov - born - green - kirkwood - yvon hierarchy ) . then , using the molecular chaos assumption ( which corresponds to the second assumption of the model ) , this equation is closed by expressing the two - particle distribution function through the single - particle one to yield the boltzmann equation . the navier - stokes equation is obtained from the boltzmann equation in the limit of small knudsen number @xmath118 , as the first - order correction in @xmath119 to the euler equation , which describes the gas motion in a state of local equilibrium ( the chapman - enskog expansion ) . therefore , although the navier - stokes equation is incomparably rich in the description of the laminar / turbulent transition , giving a dynamic picture of the process , it is reasonable to expect that the statistical properties of the flow calculated on the basis of the navier - stokes equation are , in a considerable degree , due to the above assumptions involved in the derivation of the navier - stokes equation ( they can be considered as a hidden statistical content of the navier - stokes equation ) . following onsager @xcite , let us consider local temperatures of flow states , which in our case will be related to the coherent structures of different size @xmath120 where @xmath121 and @xmath122 are the `` cumulative '' distributions of the energy and entropy , respectively , and the increment is calculated as @xmath123 . the quantities @xmath124 and @xmath125 play a role of the energy @xmath1 and entropy @xmath0 in the onsager theory @xcite . in particular , the energies are similar in two important respects . first , @xmath124 increases as the coherent structures of larger size are included into consideration , similar as @xmath1 increased because of clustering of vortices of the same size in the onsager theory . secondly , since the hamiltonian for the system of two - dimensional point vortices in inviscid fluid is the part of the kinetic energy of the system that depends on the relative positions of the vortices ( batchelor @xcite ) , the total energy @xmath1 considered by onsager represents the kinetic energy , i.e. , similar to what we assumed for the total energy of the system of coherent structures ( sect . the difference with the onsager theory is that in our case , because we consider a viscous fluid , the temperature distribution also depends on the reynolds number @xmath126 as on a parameter . taking into account eqs . ( [ eq1 ] ) , ( [ eq3 ] ) and ( [ eq3_2 ] ) , we have @xmath127 and @xmath128 where @xmath129 and @xmath130 . then , assuming @xmath9 and @xmath17 to be fixed , so that the reynolds number determined by eq . ( [ eq4 ] ) is fixed too , eq . ( [ eq5 ] ) gives @xmath131 to represent the temperature in the subsequent figs . [ t(q)]-[t_crit(re ) ] , the constant @xmath132 is set to unity . examination of the temperature distributions given by eq . ( [ eq8 ] ) shows that they are characteristically different in different ranges of variation of the reynolds number . specifically , the distributions of three types are observed , which are depicted in fig . [ t(q ) ] . at small @xmath83 the temperature is positive for all @xmath84 ( the red curve ) , at large @xmath83 it is negative for all @xmath84 ( the blue curve ) , and at intermediate values of @xmath83 , the temperature abruptly changes from positive to negative at some @xmath133 , which corresponds to the maximum of the entropy ( the black curve ) . in the latter case , the behavior of the temperature is qualitatively similar to what was predicted by onsager @xcite . the essential difference is that in the present case the distribution is dependent upon the reynolds number as on a parameter ; in particular , the value of @xmath134 changes with @xmath83 , although not regularly because different combinations of parameters @xmath70 and @xmath71 entering eq . ( [ eq8 ] ) can lead to very close values of @xmath99 ( as , e.g. , fig . [ q(re ) ] evidences ) .    to see how the temperature of state changes with the reynolds number , let us introduce average temperatures @xmath135 which are characteristic of different @xmath99 . if the temperature is positive or negative for all @xmath84 , as for @xmath136 and @xmath137 in fig . [ t(q ) ] , it was calculated as an overall average @xmath138 , where @xmath91 is the most probably size distribution of the structures , eq . ( [ eq3 ] ) . otherwise , if the sign of the temperature changes , as for @xmath139 in fig . [ t(q ) ] , it was calculated for the @xmath140 and @xmath141 segments of the dependence @xmath142 separately . the results of the calculations are presented in fig . [ t(re ) ] . it is seen that the overall average temperature drastically diverges , changing from positive to negative as the reynolds number increases . this divergence of the overall average temperature is accompanied by large fluctuations of the segment average temperatures in a relatively wide range of the reynolds numbers . the state of the flow that has the temperature of the same sign for all coherent structures , i.e. , the overall average temperature , should be associated with a certain state of the flow , which is laminar or turbulent . therefore , the laminar and turbulent states , which correspond to small and large reynolds numbers , are characterized by positive and negative temperatures , respectively . correspondingly , the range of the reynolds numbers where the positive and negative temperatures are in mixture should be associated with the laminar - turbulent transition range . figure [ t(re ) ] suggests that the laminar - turbulent transition is similar to a second - order phase transition @xcite in that at some critical reynolds number the temperature drastically diverges , which is accompanied by large fluctuations of the temperature . a similarity between these two phenomena has been noticed and discussed in refs . , where laminar - turbulent transition was studied experimentally . chabaud et al . @xcite investigated a jet of gaseous helium and found that in the transition range the behavior of the width of the velocity probability density function depending on the size of the spatial structures resembled the isotherms of a liquid - gas system near the critical point . tabeling and willaime @xcite studied the flatness of the velocity derivatives in the flow of gaseous helium between counter - rotating disks ( a von krmn swirling flow ) and showed that in the transition range the flatness experiences significant changes and has a power - law dependence on the reynolds number as @xmath143 below the critical number @xmath144 and @xmath145 above @xmath144 . more recently , cortet et al . @xcite investigated the susceptibility of the von krmn swirling liquid water flows to symmetry breaking initiated by different rotation frequencies of the disks and found that at a critical reynolds number @xmath144 the susceptibility diverges as @xmath146 . the quantity @xmath147 was interpreted as a temperature due to castaing @xcite , who introduced the turbulent flow temperature as a quantity that is conserved along the length scales in the energy cascade . the present model shows a similar qualitative behavior of the state temperature as a function of the reynolds number @xmath83 . to determine critical exponents , we averaged the data of fig . [ t(re ) ] by dividing the reynolds numbers into bins of length 0.1 ( all temperatures were taken into account , not only those shown in fig . [ t(re ) ] ) . in each bin , the mean values of the reynolds number ( @xmath148 ) and the temperature ( @xmath149 ) were calculated . the obtained dependence of @xmath149 on @xmath150 was fitted to the function @xmath151 separately for @xmath152 and @xmath153 , where @xmath154 and @xmath155 are the constants to be fitted , and @xmath156 is an expected critical reynolds number ( fig . [ t_crit(re ) ] ) . the constant @xmath155 plays a role of the critical exponent . to determine constants @xmath154 and @xmath155 for a selected value of @xmath156 , the functional @xmath157 ^ 2 $ ] was evaluated , where @xmath13 is the number of the bin , and @xmath158 is the standard deviation of @xmath159 in @xmath13 bin . the functional was minimized with respect to @xmath154 and @xmath155 in a close vicinity of @xmath156 ( @xmath160 ) . to find an optimal value of @xmath156 , the procedure was repeated for different values of @xmath156 that increased from 1.5 to 1.9 . the value of @xmath156 was considered to be optimal if the both final values of @xmath161 for @xmath162 and @xmath153 were smaller than the corresponding values of @xmath161 for the other values of @xmath156 . the best fit was achieved at @xmath163 , in which case the data are approximated by the functions @xmath164 at @xmath152 and @xmath165 at @xmath153 ( fig . [ t_crit(re ) ] ) . the fluctuations of the temperature sharply increase toward the critical reynolds number ( exponentially rather than by a power law ) and reach @xmath166 and @xmath167 at @xmath168 and @xmath153 , respectively . the observed behavior of the temperature with the reynolds number is similar to second - order phase transitions in several essential respects : the temperature exhibits a power - law divergence at the critical point ( the critical reynolds number ) , the critical exponents below and above the critical point are practically equal to each other and have the same order of magnitude as those in second - order phase transitions , and the fluctuations of the temperature in the vicinity of the critical point are as large or larger than the mean value of the temperature ( stanley @xcite ) . at the same time , drawing a direct analogy to second - order phase transitions is problematic , primarily because the flow can not be divided into parts without a loss of its essential properties ( sect . ii ) . consequently , the system of particles representing the flow does not have the property of self - similarity , which plays a key role in second - order phase transitions and critical phenomena ( kadanoff @xcite and wilson @xcite ) . the drastic change of the state temperature at @xmath139 offers a criterion to distinguish between laminar and turbulent states that is more definite than the dependence of @xmath87 upon @xmath83 ( fig . [ q(re ) ] ) . the unreasonably low value of @xmath156 ( @xmath169 ) we obtained is , in a considerable degree , a result of the flexibility in the definition of the size of the elementary volume . equation ( [ eq4 ] ) assumed that it was the product of the kolmogorov length and velocity dissipation microscales @xmath170 . however , the space scale must not be necessarily equal to @xmath62 . for example , the scale dividing the dissipation and inertial subranges @xmath171 ( pope @xcite ) can be taken as the space scale . correspondingly , because in the inertial subrange @xmath172 , @xmath173 can serve as the velocity scale . then , the value of the critical reynolds number @xmath156 increases @xmath174 times , which leads to a more reasonable value of @xmath175 . it would be of interest to determine the statistical temperature of the flow by experiment or numerical simulations and test the predicted behavior of the temperature with the reynolds number . for this , eq . ( [ eq8 ] ) can be rewritten as @xmath176 , where @xmath30 is the kinetic energy of @xmath13 structure , @xmath38 is the size distribution of the structures , and @xmath177 , where @xmath55 is the volume of @xmath13 structure , and @xmath62 is the kolmogorov microscale . according to this equation , all one needs to know to calculate the temperature is the volumes of the coherent structures , their kinetic energies , and the size distribution of the structures . such information can be obtained , in principle , by direct numerical simulations , similar to those of moisy and jimnez @xcite , in which the size distribution of the coherent structures was determined and the kinetic energies of the structures could be calculated . of the structures could be calculated . using the previously developed statistical model of laminar / turbulent states of a fluid flow , which considered the flow as a collection of coherent structures of various size @xcite , a statistical temperature of the flow states has been determined that depends on the reynolds number . it has been shown that at small reynolds numbers , associated with laminar states , the temperature is positive , while at large reynolds numbers , associated with turbulent states , it is negative . at intermediate reynolds numbers , the temperature abruptly changes from positive to negative as the structure size increases , similar to what onsager predicted for a system of parallel point - vortices in an inviscid fluid @xcite . it has also been shown that in the range of intermediate reynolds numbers , which is associated with the transition between laminar and turbulent states , the temperature drastically diverges , exhibiting a power - law behavior characteristic of second - order phase transitions . the critical exponents below and above critical reynolds number are close to each other and have the same order of magnitude as those in the second - order phase transitions . i thank r. khairulin for useful discussions of the critical phenomena .
Abstract:

Answer: using the previously developed model to describe laminar / turbulent states of a viscous fluid flow , which treats the flow as a collection of coherent structures of various size ( chekmarev , chaos , 2013 , 013144 ) , the statistical temperature of the flow state is determined as a function of the reynolds number . it is shown that at small reynolds numbers , associated with laminar states , the temperature is positive , while at large reynolds numbers , associated with turbulent states , it is negative . at intermediate reynolds numbers , the temperature changes from positive to negative as the size of the coherent structures increases , similar to what was predicted by onsager for a system of parallel point - vortices in an inviscid fluid . it is also shown that in the range of intermediate reynolds numbers the temperature exhibits a power - law divergence characteristic of second - order phase transitions .


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: experimental realization of bose - einstein condensation in a dilute gas of alkali atoms @xcite stimulated a number of theoretical studies of the quantum degenerate , weakly interacting many body system . while the lowest energy state of the celebrated gross - pitaevski equation @xcite describes remarkably well the properties of the condensate at zero temperature , the dynamics of the bose gas at non - zero temperatures remains a challenge . the most successful description of the temperature dependence of the condensate oscillation frequencies and their damping rates has been obtained using different versions of the two gas models @xcite . while useful , these models assume from the very beginning that the system at finite temperatures consists of two distinct factions : the condensate and the thermal cloud . a more fundamental approach would deduce the very existence of the condensed part directly from the dynamics of the many body system . in a series of papers , several groups @xcite proposed a classical field approximation to fulfill this goal . in this approximation , the bosonic field operator is replaced by a classical field . this dramatic simplification is rooted in quantum electrodynamics , where classical electric and magnetic fields arise naturally if the number of quanta is large in a given mode . the classical field approximation is a simple and convenient technique describing the condensate in dynamical equilibrium with the thermal cloud , at temperatures close to the critical one . we have shown that the whole isolated system may be viewed as a single classical field undergoing nonlinear dynamics leading to a steady state @xcite ( see also @xcite ) . the condensate is defined as the dominant term in the spectral decomposition of the time - averaged single - particle density matrix . two cases were discussed : the box with periodic boundary conditions @xcite and the realistic case of a spherically symmetric harmonic trap @xcite . in our procedure it is the observation process and the finite detection time that allow for splitting the system into the condensate and the thermal cloud . the aim of the present paper is to further corroborate the details and the applicability of the classical fields approximation .    in section [ method ] we briefly review the approximation based on the classical fields . we stress again the need for the time averaging of the single particle density matrix for unambiguous splitting of the system into the condensed and uncondensed phases . the classical field approximation is introduced with the high momentum cut - off . in section [ box ] we study in some detail the role of this cut - off in the model of a cubic box with the periodic boundary conditions . we also analyze the dependence of the cross - over energy on the value of the coupling constant . within our approximation , it is a growing function of the scattering length . in section [ vortices ] , for the first time , we apply the method to a dynamical dissipative process . we phase imprint the vortex on a finite temperature , partially condensed bose gas and track its decay . we conclude with the summary in section [ summary ] . we start with a brief description of our approach ( some other details can be found also in@xcite ) . first , we consider the simplest case of a system of @xmath0 bosonic atoms of mass @xmath1 interacting via two - body forces . the @xmath0-body hamiltonian written in the second quantization formalism has a form : @xmath2 where the first term is the kinetic energy and the second one is the energy of two - body interactions . in writing the hamiltonian ( [ hamiltonian ] ) we made a standard assumption : in the low energy limit a two - body interaction potential is approximated by a zero - range one . the interaction strength @xmath3 is determined by a single parameter  the s - wave scattering length @xmath4 , @xmath5 . @xmath6 is the field operator satisfying equal - time bosonic commutation relations :    @xmath7=\delta({\bf   r}-{\bf r}^{\prime}).\ ] ]    the field operator can be expanded into any complete set of single - particle wave functions . one of these basis sets is of particular importance  the one which corresponds to eigenfunctions of a single - particle density matrix@xcite . this is because these functions have direct physical interpretation as they define coherent modes of the system . these modes are related to , typically performed , single - particle measurements . note that a theoretical description of correlated measurements , which involve a simultaneous detection of a few particles , should be based on an appropriate multi - particle reduced density matrix . the single - particle mode of macroscopic occupation is a bose condensate . in the ideal gas case eigenfunctions of the single - particle density matrix are obviously one particle states of the external potential . the situation becomes more complex if particles interact . in general , the eigenmodes are unknown . the only exception is a stationary system of particles trapped in a box with periodic boundary conditions . here the symmetry of the problem imposes very strong constrains and eigenfunctions of the single - particle density matrix are simply plane waves  even in the presence of interactions . the field operator is therefore : @xmath8 where @xmath9 is a bosonic operator which annihilates a particle of momentum @xmath10 . a full operator solution of the heisenberg equation originating from the above hamiltonian is not available and some approximations are necessary . following bogoliubov , at zero temperature , when all particles occupy a single zero - momentum mode , the corresponding annihilation operator can be substituted by a c - number amplitude , @xmath11 . this is possible because the commutator @xmath12 = 1 \ll   n$ ] is much smaller than the number of particles occupying the zero - momentum state . the remaining terms of the expansion ( [ modeexpand ] ) represent quantum corrections which in the lowest order of approximation are neglected . let us note that the procedure leads to a substitution of the the full field operator by a c - number wave function :    @xmath13    the condensate wave function satisfies the gross - pitaevskii equation which is successfully used for description of low temperature behavior of the bose - einstein condensate : @xmath14 this equation preserves both the energy and the particle number , i.e. the normalization of the wave function @xmath15 . the gross - pitaevskii equation can also describe trapped condensates when a term describing the external potential is included . our approach is a simple generalization of the above procedure . let us notice first that at relatively high energies ( but below the critical one ) there exists a number of modes , say @xmath16 , whose occupation significantly exceeds unity , i.e. many different momentum states , up to a certain @xmath17 are highly populated . consequently , all corresponding annihilation operators in the expansion ( [ modeexpand ] ) can be substituted by complex amplitudes and the remaining terms may be neglected :    @xmath18    this approach evidently leads to the same gross - pitaevskii equation ( [ gp ] ) for a wave function describing the `` macroscopic part '' of the whole system , or , more precisely , to its finite lattice version . this kind of approach is frequently used in quantum electrodynamics . the only difference from the zero temperature case is that now this wave function does not correspond to the minimal energy solution . instead , it is a relatively high energy state . as it has been shown in@xcite , regardless of the particular choice of an initial state @xmath19 , different than the gross - pitaevskii eigenstate , the system evolves towards the same stationary state uniquely determined by the energy and the particle number . more precisely , occupations of all macroscopically populated eigenmodes @xmath20 fluctuate around their mean values@xcite . the self consistency of the c - number approximation requires that all populations have to be large , @xmath21 . to this end , the number of classical modes @xmath16 has to be carefully chosen for each energy and particle number@xcite . this fact can be checked only _ a posteriori _ when the stationary populations are obtained . the intimate link between the energy , the number of particles and the number of modes is the essential ingredient of the whole method . there is still one point which requires additional discussion . in fact , the `` quasi - stationary '' high energy solution of the gross - pitaevskii equation represents a pure state of the system with _ all particles _ in the same single state @xmath22 . therefore , it seems to be totally unjustified to interpret @xmath23 as relative occupations of different coherent modes of the system . the last concept is associated with a mixed state . detailed inspection shows , however , that the relative phases of complex amplitudes @xmath9 vary in time rapidly . on the other hand an observation process lasts for a finite period of time , typically of the order of hundreds of microseconds . it is sufficiently long for a `` decoherence '' of different eigenmodes :    @xmath24    in the stationary regime the time averaged quantities do not depend on time . on the contrary , a temporal wave function of the whole system varies very rapidly . therefore , what is being observed is related to a time averaged single - particle density matrix @xmath25 . it is the observation that reduces a pure state to a mixed one . for particles trapped in a periodic box we have :    @xmath26    eq.([time_box ] ) shows that eigenvalues of the time averaged single - particle matrix are normalized occupations of different coherent modes of the system . we want to stress that the averaging procedure is essential for the correct physical interpretation of the high energy solutions of the gross - pitaevskii equation . note that instead of averaging over time , spatial avergaing can be performed . their importance can be fully appreciated in a more realistic case of harmonically trapped bose condensates when eigenmodes of a single - particle matrix are not known . it is just diagonalization of the time averaged density matrix which gives simultaneously the eigenmodes and the populations . results of our calculations for a harmonically trapped bose condensate at finite energy are presented in@xcite . the time averaged density has a characteristic bimodal pattern : sharp peak of condensed atoms embedded in a broad thermal cloud . on the other hand , an instantaneous distribution consists of a number of irregular spiky structures@xcite . according to our numerical calculations , the average width of a single spike is about three times larger than the healing length , the distance over which the condensate wave function heals back when perturbed locally . the characteristic length scale of the irregular structure is slightly larger ( but of the same order ) than the grid spacing necessary for description of the system of a given energy and particle number . as we have already mentioned , the self - consistency of the approximation requires a very precise selection of a number of modes @xmath16 what in the present case directly translates into the grid spacing . in our approach @xmath16 becomes an important , physical parameter . in this section we explore our approach and present some important technical details . we consider here the system trapped in a periodic box because in this situation the coherent eigenmodes are uniquely determined by the symmetry of the potential . this greatly simplifies the problem . we numerically solve the gross - pitaevskii equation ( [ gp ] ) on a finite grid using the fast fourier transform split - operator method . the number of grid points is equal to the number of different momentum states . in order to explore the link between the number of macroscopically occupied modes , the number of particles and the energy we solve eq.([gp ] ) in a 3d geometry varying the number of eigenmodes @xmath16 only . in the calculations we keep the same initial state @xmath19 which guarantees that the energy per particle is constant . similarly , the value of the following product is fixed :    @xmath27    the results might seem surprising : the larger the number of modes @xmath16 ( with other parameters constant ) the higher the occupation of the condensate mode @xmath28 . this result requires a more detailed discussion . the total energy , the number of modes @xmath16 and the product of @xmath29 are the _ only _ control parameters in our method . the number of particles @xmath0 as well as the coupling strength @xmath3 do not enter the dynamical equation ( [ gp ] ) separately  it depends only on the product @xmath29 . nevertheless , the condition of validity of the classical field approximation allows to determine the above parameters quite accurately . an occupation of the highest eigenmode @xmath30 must be larger then one . we arbitrarily set this value to 5 . this way , having a stationary value of a relative occupation @xmath31 , we can determine both @xmath0 and @xmath3 . it occurs that the requirement of self - consistency indicates that in our calculations the number of particles grows with @xmath16 while the interaction strength @xmath3 decreases in order to satisfy the condition ( [ limit ] ) . only by examining the link between all parameters of the classical field approximation can we fully interprete our data . fraction of particles in the condensate @xmath32 as a function of particle number for _ gn = const_. for comparison we show the curve given by the equation @xmath33 . ]    in fig.[mean_vs_np ] a relative abundance of the condensate atoms is plotted as a function of the number of particles . the condensate occupation can be quite accurately approximated by the following analytic expression :    @xmath34    this behavior is in agreement with a recent finding@xcite . lieb and seiringer studied the properties of a true ground state @xmath0-body wavefunction of a system of @xmath0 bose particles interacting via two - body forces . they have shown that in the limit of @xmath35 with @xmath36 the relative population of the dominant eigenmode of the corresponding single particle density matrix aproaches @xmath37 while the eigenmode tends towards a ground state solution of the gross - pitaevskii equation . our method shows how this limit is approached .    in the next series of calculations we estimate the interaction - induced shift of the cross - over energy for bose - einstein condensation . we prefer to use the term `` cross - over energy '' instead of `` critical energy '' as , strictly speaking , the critical energy is defined in the thermodynamic limit only . for a finite system ( like the one we deal with here ) several measures of the corresponding `` cross - over '' energy can be defined as there is no sharp phase transition . one of them @xcite is based on the rapid change of the probability distribution of the condensate occupation . we study the occupation of the condensate mode @xmath38 . the total energy of the system is changed while @xmath0 and @xmath3 are kept constant . this requires an optimization of the number of classical fields @xmath16 to be taken into account : again , we set the occupation of the least occupied mode to 5 ( @xmath39 ) . the above constraint does not allow for a continuous change of the energy though . when a quasi - stationary state of a given energy is reached we trace a temporal population of the condensate mode @xmath40 for a long time ( about @xmath41 time steps ) . this way we obtain a probability distribution of the condensate population which is presented in fig.[histogram ] . in this figure we show two characteristic distributions for different energies of the system . one of them ( for the energy per particle @xmath42 ) is centered around some non - zero mean value of @xmath43 ( the grey histogram ) , while the other ( for the energy @xmath44 ) is peaked at @xmath45 . here @xmath46 is the energy of the first excited state . statistics of the condensate occupation in time . the grey histogram presents the case slightly below the cross - over energy ( energy per particle @xmath47 ) , while the transparent histogram corresponds to the situation above the cross - over energy ( energy per particle @xmath48 ) . note that the highest bar of the latter histogram is off the scale of the figure ( its value is 0.87 ) . here @xmath49 and @xmath50 . ]    despite the fact that the values of energy for which both histograms are plotted are very similar , the rapid change of the character of the statistics evidently signifies the onset of bose - einstein condensation . the two energies give the lower and upper limit on the critical energy for a given interaction strength . we can not determine the critical energy with better accuracy because of the previously discussed link between the parameters of our method . in the fig.[ekrit_vs_g ] we show the upper and lower limits for the critical energy as a function of the coupling strength @xmath3 . the critical energy grows with @xmath3 . the same must happen to the critical temperature . we do not estimate the critical temperature as the method based on fitting the thermal distribution to the outer wings of the time averaged density profiles introduces a large uncertaintity@xcite . there is still a controversy in the literature@xcite about the magnitude and the sign of the shift of critical temperature due to interactions . the accuracy of our calculations , at the present stage , does not allow for a precise determination of this quantity . nevertheless , they certainly indicate that the critical energy ( and temperature ) in the case of gas trapped in a box grows with interactions . the same conclusion for the classical field simulations in a box has been reached recently in @xcite . estimates of the critical energy of the phase transition based on the analysis of pairs of histograms like in fig.[histogram ] . total number of particles is @xmath49 . ] creation of quantum vortices has provided a long - sought clear evidence of supefluidity in bose - einstein condensates of dilute atomic gases @xcite . following the two pioneering experiments , various mechanisms of vortex nucleation have been investigated @xcite . more recently , experiments have focused on peculiar properties of vortex lattices @xcite .    in most experiments vortices have been created in pure , essentially zero temperature , condensates . very little has been said about finite temperature properties of these systems . @xcite , operating at temperatures below their detection limit , have measured finite vortex lifetimes of the order of 500 - 1000 ms . _ @xcite have seen no thermal damping on the time scale of 1 s. finite temperature decay of vortex lattices has been also investigated @xcite . similarly , few theoretical works approached the problem . fedichev and shlyapnikov @xcite show that it is the interaction of the thermal cloud with the vortex that provides a mechanism of energy dissipation . due to that , an off - center vortex spirals out to the condensate boundary where it decays through the creation of excitations . for typical experimental parameters , they estimate the vortex liftime to be in the range 0.1 - 10 s. temperature dependence of the critical rotation frequency for the creation of vortices has been calculated in @xcite . bogoliubov as well as hartree - fock - bogoliubov - popov theories have been used to study the stability of vortex states at finite temperature @xcite . the role of filling the vortex core by the thermal cloud as a stabilizing factor has been emphasized @xcite . trajectory of a vortex core . the intensity of lines and dots fades away with time . ] we have used the classical field method to study the dynamics of vortices in the presence of a thermal cloud . as an example , we consider @xmath41 @xmath51rb atoms in a spherically symmetric trap of frequency @xmath52 hz . first , a suitable phase pattern ( corresponding to a single quantum of circulation ) is imprinted on a finite - temperature equilibrium state of an interacting bose gas . experimentally , phase imprinting is realized by shining a far off - resonant light on the condensate through a glass plate whose opacity varies continuously with an azimuthal angle . as a result , the condensate atoms acquire a phase depending on their location @xcite . we situate the core of the vortex off the trap center . we see that the vortex created in this way starts moving around the trap center , slowly spiraling out to the border of the condensate . we follow the trajectory of the vortex core by time - averaging the gas density over short time . a sample trajectory of the vortex is shown in fig.[trajectory ] . the equilibrium state used here for the phase imprinting contained about 48% atoms in the condensate ( its energy being equal to 49.46 @xmath53 ) . the initial displacement of the vortex core from the trap center is @xmath54 . the thomas - fermi radius of a condensate consisting of 48000 @xmath51rb atoms in our trap is @xmath55 ( @xmath56 m being an oscillator unit of length )  note that the vortex disappears reaching roughly this distance from the trap center . time dependence of the distance of the vortex core from the trap center . time expressed in units of a trap period @xmath57 . ] fig.[radius ] shows the time dependence of the distance of the vortex core from the trap center . one can see a slow drift towards the condensate boundary superposed on top of small fluctuations coming from the interaction with the thermal part of the system . from fig.[radius ] one can conclude that the radial velocity of the vortex is roughly constant . note also the vortex lifetime @xmath58 ( here equal to about @xmath59 trap periods ) . angular velocity of the vortex line ( expressed in units of @xmath60 ) as a function of its distance from the trap center . for comparison , an estimate of fetter and svidzinsky for a pure condensate of 48000 @xmath51rb atoms in our trap is shown . ] fig.[phidot ] shows the angular velocity of the vortex core as a function of its ( varying with time ) distance from the trap center . one can see that the vortex moves faster and faster , in particular in the neighborhood of the condensate boundary . this result can be compared with theoretical predictions for pure , zero temperature , condensates  however , in this case there is no dissipation and off center singly quantized vortices in our case would follow _ closed _ circular trajectories . for pure condensates the angular velocity of the vortex line diverges at the condensate boundary @xcite . we also see this kind of behavior in our data for a partly condensed bose gas . to be specific , we have compared our data with a variational estimate provided by eq.(87 ) in @xcite ( here @xmath61 , as the trap is nonrotating ) taking 48000 @xmath51rb atoms . despite notable differences between the two situations ( presence or absence of dissipation due to the thermal cloud ) the agreement is rather good even at the quantitative level . 3-dimensional density of the gas cloud with a bent vortex line . the outer shell is a density isosurface near the condensate boundary . a plane with a surface density plot , perpendicular to the vortex line and through the trap center , is also shown . ] more features of the vortex structure can be seen in fig.[wir ] . in this figure a 3-dimensional density of the whole gas cloud , averaged over a short time , is depicted . the outer shell is an isosurface of density near the boundary of the condensate . the whole vortex line is visible in the upper part of the figure : one can see that it is bent at the edge of the condensate and it is much thicker there than in the center of the trap . the former observation agrees with recent experimental and theoretical results @xcite . a plane perpendicular to the vortex line and through the trap center is also shown  it contains a surface density plot with a distinct hole at the vortex line . so far we have kept the initial displacement of the vortex core fixed and simply studied the features of its trajectory . in the following , we look at the dependence of the vortex lifetime on its initial position in a trap . fedichev and shlyapnikov @xcite calculate that :    @xmath62    where @xmath63 is the spatial size of the condensate ( well approximated by the thomas - fermi radius ) , @xmath64 is the initial displacement of the vortex line from the trap center , @xmath57 is temperature and @xmath65 is the mass density of the thermal cloud ( for details of the derivation see @xcite ; however , note the misprint in eq.(12 ) of @xcite  the bracketed factor should be raised to the power of @xmath66 and not @xmath67 @xcite ) . dependence of vortex lifetime @xmath58 on its initial displacement from the trap center @xmath64 . points come from our simulations while the solid line corresponds to eq.([log ] ) . ]    in fig.[lifetime ] we plot the vortex lifetime as a function of its initial position ( calculated from our simulations ) as well as the values given by eq.([log ] ) ( note that temperature of our system , for a given condensate fraction , has been estimated with the help of the ideal - gas formula with the finite - size corrections taken into account ) . one can see that it does agree quite well with the analytical estimate . however , one should note that one of the assumptions leading to the derivation of eq.([log ] ) is that the condition @xmath68 holds . in our case @xmath69 which may be one of the reasons of the slight discrepances between our data and the analytically predicted approximate values . note that the total energy of the system is still fixed to 49.46 @xmath53 . we have also investigated the variation of the vortex lifetime with the total energy ( which amounts to its dependence on the condensate fraction ) . we have seen that this relation is much weaker than the one studied above . we have also noticed that even for a fixed initial displacement ( which does not specify its position uniquely ) the lifetime changes from shot to shot  in other words , the complex violent behavior of the thermal cloud in different regions of the trap makes the vortex lifetime a stochastic parameter . we have reviewed our version of the classical field approximation to the dynamics of a finite temperature bose gas . it provides an unambiguous splitting of the system into the condensed and uncondensed parts with the help of the time averaging of the single particle density matrix . in the case of a periodic box potential , we have investigated the role of the high momentum cut - off , essential in the method . in particular , we have shown that the cut - off going to infinity limit decribes the particle number going to infinity with the scattering length going to zero . in this weak interaction limit , the relative population of the condensate tends to unity . we have also shown that the cross - over energy , at which the probability distribution of the condensate occupation changes its character  one of several measures of the energies signifying the phase transition for the finite systems  grows with a growing scattering length . in the more physical case of the condensate in the harmonic trap we have investigated the dissipative dynamics of a vortex . we have compared the decay time and the velocities of the vortex with the available analytic estimates . the next step in going beyond the approximations employed in our model is the estimation of the influence of quantum corrections to the classical fields , in particular the study of the corresponding time scales . anderson , j.r . ensher , m.r . matthews , c.e . wieman , and e.a . cornell , science * 269 * , 198 ( 1995 ) ; k.b . davis , m .- o . mewes , m.r . andrews , n.j . van druten , d.s . durfee , d.m . kurn , and w. ketterle , phys . lett . * 75 * , 3969 ( 1995 ) ; c.c . bradley , c.a . sackett , j.j . tollett , and r.g . hulet , _ ibid . _ * 75 * , 1687 ( 1995 ) and erratum * 79 * , 1170(e ) ( 1997 ) . svistunov , j. moscow phys . soc . * 1 * , 373 ( 1991 ) ; k. damle , s.n . majumdar , and s. sachdev , phys . a * 54 * , 5037 ( 1996 ) ; yu . kagan and b.v . svistunov , phys . * 79 * , 3331 ( 1997 ) ; n.g . berloff and b.v . svistunov , phys . a * 66 * , 013603 ( 2002 ) ; k.staliunas , cond - mat/0110258 . o. penrose , phil . mag . * 42 * , 1373 ( 1951 ) ; o. penrose and l. onsager , phys . rev . * 104 * , 576 ( 1956 ) ; c.n . yang , rev . phys . * 34 * , 694 ( 1962 ) ; for a recent discussion , see also a.j . leggett , rev . phys . * 73 * , 307 ( 2001 ) . madison , f. chevy , v. bretin , and j. dalibard , phys . lett . * 86 * , 4443 ( 2001 ) ; c. raman , j.r . abo - shaeer , j.m . vogels , k. xu , and w. ketterle , _ ibid . _ * 87 * , 210402 ( 2001 ) ; p.c . haljan , i. coddington , p. engels , and e.a . cornell , _ ibid . _ * 87 * , 210403 ( 2001 ) . virtanen , t.p . simula , and m.m . salomaa , j. phys . c. : condensed matter * 13 * , l819 ( 2001 ) ; s.m.m . virtanen , t.p . simula , and m.m . salomaa , phys . 86 * , 2704 ( 2001 ) ; t. mizushima , t. isoshima , and k. machida , phys . rev . a * 64 * , 043610 ( 2001 ) .
Abstract:

Answer: we review our version of the classical field approximation to the dynamics of a finite temperature bose gas . in the case of a periodic box potential , we investigate the role of the high momentum cut - off , essential in the method . in particular , we show that the cut - off going to infinity limit decribes the particle number going to infinity with the scattering length going to zero . in this weak interaction limit , the relative population of the condensate tends to unity . we also show that the cross - over energy , at which the probability distribution of the condensate occupation changes its character , grows with a growing scattering length . in the more physical case of the condensate in the harmonic trap we investigate the dissipative dynamics of a vortex . we compare the decay time and the velocities of the vortex with the available analytic estimates .