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 dynamics of the magnetic domain wall ( dw ) driven by electric current - induced spin - transfer torque ( stt ) is under extensive both experimental@xcite and theoretical@xcite investigations in recent years . as known to all , in general the dw motion can be driven by an external magnetic field@xcite and/or spin - polarized electric current.@xcite for the magnetic field - driven case , although the well - known walker s theory@xcite is usually used to understand the phenomena of dw movement induced by magnetic field , the origin of this motion is beyond the scope of walker s theory . however , x. r. wang _ et al_.@xcite recently found that the mechanism of dw propagation by external magnetic field in a nanowire can attribute to the energy dissipation that is owing to gilbert damping . on the other hand for the current - induced case , when an electron passes through a magnetic dw , it will be scattered by the noncollinear magnetic structure , which results in the phenomenon of magnetoresistance ( mr).@xcite meanwhile , the conduction electron can also transfer the spin angular momentum to the local spin when it flows through the dw . so a spin - transfer torque will exert on the local magnetic moment and then the electric current can be used to manipulate the magnetic structure of dw . this phenomenon of stt including current - induced dw motion was first predicted and observed in the experiments by l. berger and coworkers.@xcite moreover , this exotic phenomenon will also have some potential applications such as magnetic random - access memories ( mrams),@xcite racetrack memories,@xcite spin transfer nano oscillators ( stnos)@xcite etc . , which have attracted a great deal of attention recently .    for a dw in the adiabatic approximation , the spin of the incident electrons can be consistently aligned with the local spin moments . the spatial derivative of this adiabatic spin current yields the adiabatic stt which is written as @xmath0 , where @xmath1 is the local spin of dw and @xmath2 is coordinate in transport direction.@xcite here @xmath3 is parameter with unit [ m / s ] in proportion to current density . however , there are always some conduction electron spins that can not follow the local spin moments as the width of dw decreasing . the process that the spin of conduction electron relax toward the local spin moments will cause a nonadiabatic torque@xcite called the @xmath4 term , which is in the form of @xmath5 , where @xmath4 is the dimensionless parameter and @xmath6 is the unit vector of local spin moments . further , s. zhang and z. li@xcite have pointed out that the adiabatic stt only contributes the initial velocity of dw movement while the nonadiabatic term mentioned above determines the terminal velocity observed in experiments .    however , for the relative narrow dw a nonlocal oscillating torques have been predicted theoretically by several groups@xcite and its quantum origin is similar to the rkky oscillation.@xcite moreover , in general the ratio between the maximum of nonadiabatic and adiabatic stt , such as the coefficient @xmath4 mentioned above , represents the degree of nonadiabaticity which can determine the velocity of dw movement.@xcite so one of the purposes in this paper is to find out whether there is a nonlocal oscillating or localized torques in real ferromagnetic materials by calculating the stt of dw using the first - principles method . and the other purpose is to study the difference in magnitude of nonadiabatic stt among the traditional ferromagnetic dws ( fe , co and ni ) .    on the other hand , m. thorwart and r. egger@xcite have derived a torque in form of the second - order spatial derivative of the local spin by a gradient expansion scheme . hence , it is tempted to introduce a higher - order of nonadiabatic stt such as @xmath7 into the current - induced dw dynamics.@xcite so in this paper we will figure out whether there is higher - order torque in dw for real ferromagnetic materials . in addition , the effect of such higher - order stt on the dw movement is also needed for better understanding of current - induced dw dynamics .    in this paper , we will study the current - induced stt of defect - free dw in ballistic limit by the first - principles electronic structure calculations.@xcite our results show that the nonadiabatic stt of fe dw has localized form while the out - of - plane stts of co and ni are typical nonlocal oscillating torques . the degree of nonadiabaticity is also much greater for fe dw and increases exponentially with decreasing the width of dw . in addition , the results also show that the distribution of our calculated nonadiabatic stt for fe dw can well describe by the third - order spatial derivative term as @xmath7 . finally , the dynamics of local spin in dw pushed by this third - order stt is simulated using landau - lifshitz - gilbert ( llg ) equation . the obtained time evolution of dw movement demonstrates that with the same amplitude of the nonadiabatic stt the velocity owing to the third - order stt will be reduced to about half of the wall speed induced by the first - order term . our calculation of stt for dw is based on the scattering wave function matching ( wfm ) method with tight - binding linear muffine - tin orbital ( tb - lmto ) basis.@xcite for a typical nel or bloch dw , the left and right domains act as leads and the wall structure is regarded as the scattering region , which defines the completed scattering problem of layered system as shown in fig . [ f1 ] . for this layered system , we assumed that the magnetization of dw has lattice translation invariant in the plane perpendicular to transport direction , so the operator and scattering states can be characterized by a lateral @xmath8 wave vector in two - dimension ( 2d ) brillouin zone ( bz ) .    in this paper , the stt is calculated from the spin current , which is defined as @xmath9 . \label{js - define}\ ] ] where @xmath10 is pauli spin matrix and @xmath11 is velocity operator . in the mixed representation for a special @xmath8 with real space quasi one - dimension ( 1d ) tight - binding model , the spin current operator from * r*th to * * r**th site is@xcite @xmath12 . \label{js - opt}\ ] ] where @xmath13 is noncollinear tb - lmto hamiltonian matrix element in the global quantum axis representation . here @xmath14 , where @xmath15 and @xmath16 are the azimuthal and magnetic quantum numbers respectively . in order to obtain @xmath13 of dw , firstly the collinear hamiltonian from self - consistent one - electron effective potential of the magnetic materials is calculated by collinear electronic structure calculation in atom sphere approximation ( asa).@xcite next , by introducing the rigid potential approximation,@xcite we rotate the above hamiltonian , which is diagonal in @xmath17 spin space at local quantum axis representation , by unity rotation matrix @xmath18 in spin space to construct the noncollinear hamiltonian in the global quantum axis representation . here the rotation angle @xmath19 is determined by the local spin moments configuration of dw , which is written as @xmath20\label{lambda}\ ] ] where @xmath21 is the the characteristic length for dw and @xmath22 is the polar angle of local spin on the * * r**th site with position @xmath23 along transport direction . here , the injected current is along fcc(111 ) direction for co and ni dw and bcc(001 ) direction for fe dw . in our calculation , we choose the azimuthal angle @xmath24 for modeling nel - type dw as shown in fig . [ f1 ] . -@xmath25 plane denote the local spin moments whose polar angle is @xmath19 . note that the current ( @xmath26 ) flows from the left ( l ) lead to the right ( r ) lead and @xmath27 is layer index . , title="fig:",width=325 ] +    using the wfm method , we can obtain the scattering wave function corresponding to the noncollinear hamiltonian of dw . therefore , for scattering state with lateral wave vector @xmath8 , the expectation value of stt acting on local spin can be determined by the difference between the incoming and outgoing spin current on * * r**th site , i.e.,@xcite @xmath28 where @xmath27 is the index of principal layer in quasi-1d model and the * * r**th site belongs to the layer @xmath27 . here the superscript @xmath29 denotes that the scattering state for evaluating expectation value is induced by the injected electron in spin @xmath30 with respect to the local quantum axis of lead . then in the linear response regime , the total torque under a small bias @xmath31 is calculated by summing all the stt of @xmath8 states in 2d bz , which can be written as@xcite @xmath32 , \label{t - v}%\ ] ] where @xmath33 is the number of @xmath8 states in 2d bz at fermi level . here @xmath34 and @xmath35 denote the stt induced by the left and right incoming electrons from the lead region respectively . in our calculations , the stts of dw are performed with a @xmath36 @xmath37-mesh points in the lateral 2d bz , which can insure convergence of the results . [ ml ] at unit bias . ( a ) @xmath2 component of stt . ( b ) @xmath25 component of stt . the solid line denote the fitting adiabatic stt in form of @xmath0 . it can be seen that our calculated in - plane stts agree well with the term of @xmath0 . ( c ) the fitting results of @xmath3 for different dw width at current density @xmath38 $ ] . here the solid line denote the predicted values of @xmath39 . , title="fig:",width=325 ] +    .the parameter @xmath3 at current density @xmath40 [ a / m@xmath41 for @xmath42[ml ] , where @xmath43 is the calculated polarization @xmath44 , where @xmath45 is sharvin conductance of leads for spin up ( down ) channel at fermi level in our transport calculation . one can see that the values of our fitting @xmath3 are close to @xmath46 . [ cols="^,^,^,^,^,^,^",options="header " , ]     in this paper , we will investigate the nonadiabatic stt and corresponding dynamics of dw movement for the typical ferromagnetic materials fe , co and ni . firstly , fig . [ f2](a ) and ( b ) show the @xmath2 and @xmath25 component of stt for dw width @xmath42[ml ] . according to the configuration of nel dw in fig . [ f1 ] , here the @xmath2 and @xmath25 component of torque is so - called in - plane stt and the @xmath47 component is the out - of - plane stt . it can be seen that all these in - plane torques including fe , co and ni agree well with the prediction of adiabatic approximation,@xcite which has the form of @xmath0 . it is noted that the stt shown in fig . [ f2](a ) and ( b ) is at unit bias and the coefficient @xmath3 is proportional to the current density . so combined with the corresponding ballistic conductance @xmath48[@xmath49 of dw , we can use formula @xmath0 to fit the in - plane stt and obtain the coefficient @xmath3 at unit current density . here table . 1 shows the fitting results of parameter @xmath3 at current density @xmath38 $ ] for @xmath42[ml ] . one can observe that our fitting coefficients @xmath3 for fe , co and ni dws are close to the values of @xmath46 , which is the prediction of the semiclassical transport theory,@xcite where @xmath43 is the spin polarization of current , @xmath26 is current density and @xmath50 is saturated magnetization of ferromagnet . here note that we use the calculated polarization @xmath44 , where @xmath45 is sharvin conductance of leads for spin up ( down ) channel at fermi level in our transport calculation , instead of the values extracted from the experiments . the fitting results of @xmath3 for different dw width are also shown in fig . [ f2](c ) . it can be seen that the @xmath3 is almost independent on dw width even for the narrow wall with @xmath51[ml ] .    . one can observe that the out - of - plane stt of fe dw has localized form while the spatial distribution of stts for co and ni dw is nonlocal oscillating . , title="fig:",width=325 ] +    next , we will concentrate on the obtained out - of - plane nonadiabatic stt which is beyond the above in - plane adiabatic stt . fig . [ f3 ] shows the spatial distribution of out - of - plane stt for different width @xmath21 in fe , co and ni dws . one can observe that for co and ni dws in the region far away from the dw center , where the magnetization is uniform , there is nonlocal oscillation in the distribution of out - of - plane torque . this nonlocal oscillating torques have already been discovered by several groups@xcite and g. tatara _ et al_.@xcite pointed out that the oscillating torques can be summed up as a collective force on the dw . in our calculations , this oscillating out - of - plane stt is contributed by the propagating states which precess around the local spin moments and has the form of @xmath52 , where @xmath53 is wave vector in transport direction of propagating state.@xcite considering that the summation of propagating states over the 2d bz in the transport calculation , the cancellation effect@xcite in the different @xmath8 make the oscillating behavior of out - of - plane torque depend on the shape of fermi surface . therefore , as shown in fig . [ f3 ] the amplitude and diffusion length of torque oscillation is significant larger for ni dw than that for co dw , which is similar to the case of layered spin valve system.@xcite for the same reason as above , in fig . [ f3](a ) and ( d ) one can see that our calculated out - of - plane stt of fe dw has localized form without oscillation far away from dw center except a small oscillation near the center region . [ ml ] in the experimental current spin polarization . the black solid line denotes the total weighted stt , while the red and blue dash lines denote the stt induced only by incoming majority and minority electron respectively . , title="fig:",width=325 ] +    in our ballistic calculations , the spin polarization of the current is different with the experimental one . so in order to simulate the case of stt in experiments , firstly using the eq . [ t - v ] we can obtain the stt induced only by the incoming majority ( @xmath54 ) and minority ( @xmath55 ) electron respectively . then the total stt can be calculated by @xmath56 , where the weighting factors @xmath57 is determined by @xmath58 . here @xmath59 is the experimental current spin polarization . [ fig - suppl ] shows the distribution of the out - of - plane stts for dws with the experimental current spin polarization . one can observe that the results do not qualitatively change at all compared to the ballistic calculations , namely the out - of - plane stt of co and ni dws have nonlocal oscillating forms while that of fe dw has the localized form . is defined as the ratio between the out - of - plane and in - plane stt at dw center , namely @xmath60 . it is noted that for ni dw with @xmath61[ml ] the oscillation of out - of - plane stt is so larger that the maximum is no longer at the center of dw . for @xmath62[ml ] the nonadiabaticity of fe and co can be well fitted by @xmath63 , where @xmath64[ml ] . , title="fig:",width=325 ] +    for the out - of - plane torque in the form of @xmath5 , it is believed that the coefficient @xmath4 , which is just the ratio between the maximum of out - of - plane and in - plane stt , determine the final velocity @xmath65 of dw motion by the expression@xcite of @xmath66 . and coefficient @xmath4 is usually viewed as the degree of nonadiabaticity in stt . however , the form of out - of - plane nonadiabatic stt in our results is not exactly the same with the @xmath4 term , but we still calculated the ratio between the out - of - plane and in - plane stt at the center of dw to demonstrate the degree of nonadiabaticity for dw with different materials . as shown in fig . [ f4 ] , one can see that with decreasing the dw width the nonadiabaticity will increase exponentially , which agree well with the calculation of free - electron stoner model.@xcite further , the nonadiabaticity for co and fe dws have the same growth rate with decreasing @xmath21 . moreover , one can also observe that the nonadiabaticity of fe dw is about one order larger than that of co and ni dws . in the meantime , the value of nonadiabaticity for fe dw even can be reach at a very high value ( @xmath67 ) for the narrow width @xmath51[ml ] . in general , there are two contributions to the nonadiabatic stt in dw . one is the spin relaxation process of the conducting electron spin toward local spin , the other is the mistracking between the incoming conducting electron spin and local spin without any relaxation . in our calculations , it is noted that there is no any dynamical relaxation process taken into account , thus only the mistracking contribution exists , which is similar to the case of free - electron stoner model.@xcite from fig . [ f4 ] , one can see that this mistracking contributed stts for co and ni are too small to be measured in experiments with reasonable dw width . therefore , it implied that the relaxation contributed nonadiabatic stt dominates in co and ni dws .     term and the torques in the third - order of the spatial derivative respectively . it shows that the term of @xmath68 can well describe the calculated out - of - plane stt of fe dw . , title="fig:",width=325 ] +    fig . [ f5 ] shows the comparison between the calculated out - of - plane stt and the torque of spatial derivative of the local spin . here the torques in first- and third - order derivative are both depicted and their maximums are fitted to the calculated stt at center . at first glance , one can see that the configuration of the @xmath4 term is much different with the calculated results from fig . in particular , the width of stt peak for @xmath4 term at the dw center is much larger than that of the calculated stt . hence , it suggested that this deviation may have a significant effect on the velocity of dw movement . however , the study of the dw movement using the first - principles calculated stt directly will cost too much computing time so that we can hardly obtain the results for a long time intervals . so in order to investigate the time evolution of dw more efficiently , here we will model a simple analytical nonadiabatic stt term substituted into the equations of dw dynamical evolution . as shown in fig . [ f5 ] , we found that compared to the @xmath4 term , the nonadiabatic stt in the form of @xmath68 , which is proportional to the third - order spatial derivative of the local spin , can well describe the obtained out - of - plane stt in fe dw . here @xmath69 is the parameter in proportion to current density . the origin of this nonadiabatic stt may be similar with the @xmath4 term . according to the gradient expansion scheme for stt@xcite , the first - order expansion will result in both adiabatic stt and nonadiabatic @xmath4 term . meanwhile , the higher - order expansion will lead to the second - order stt.@xcite therefore , based on the gradient expansion scheme , it is reasonable to obtain the torques in the third - order spatial derivative in our calculations .    on the other hand , at two sides of fe dw center region there is still small difference between the calculated and the third - order stt , where the calculated stt also has small oscillating form . it implied that the reason of our calculated stt in fe dw being well fitted by the third - order term maybe is due to the dephasing effect , where the oscillating stt can be canceled each other by summing the @xmath8 states at fermi surface . for co and ni dw the majority fermi surface is similar to that of free - electron thus the dephasing effect is very weak , which makes the nonadiabatic stt still has oscillating form@xcite as shown in fig . [ fig - suppl](b ) and ( c ) . however , the fermi surface of fe is far beyond the free - electron - like one so the dephasing effect is much stronger in fe dw . meanwhile , the dephasing effect is not effective near the dw center and increases its effectiveness relatively rapidly with increasing distance from the dw center . therefore , the corresponding nonadiabatic stt has the form of the third - order term . in particular , when the distance from dw center is not very large , the dephasing has not taken fully effect in such regions , so the oscillating behavior of stt will still emerge as shown in fig . [ f5 ] . nevertheless , due to the dynamics of dw is mainly attribute to the stt in dw center , which is well fitted by the third - order stt as shown in fig . [ f5 ] , in this paper we will focus on how the nonadiabatic stt term in the form of @xmath68 will influence the velocity of dw movement .     for fe dw with different width at current density @xmath38 $ ] , where @xmath70 is the distance between the monolayers and the unit time is @xmath71 . the dependence of the third - order of spatial derivative at dw center ( @xmath72 ) on width @xmath21 is also depicted in this figure . it can be seen that although @xmath69 grows linearly with increasing @xmath21 , the amplitude of the third - order torque still decays due to the term of @xmath73 . , title="fig:",width=325 ] +    for fe dw with different width @xmath21 at current density @xmath38 $ ] , the fitting results of dimensionless coefficient @xmath74 is shown in fig . [ f6 ] , where @xmath69 is in unit of [ m@xmath75/s ] , @xmath76 is the distance between the monolayers and the unit time is @xmath77 $ ] . here , our obtained dimensionless coefficient @xmath74 can provide the magnitude of the third - order torque which maybe included in some micrmagnetic simulations such as oommf . in addition , it is noted that the factor of the third - order torque @xmath69 increases linearly with increasing @xmath21 . however , the exponential decay of @xmath73 term will lead to the obtained third - order torque still decreasing with increasing @xmath21 . [ ml ] at @xmath38 $ ] below the critical current density . here the unit time is @xmath78 . one can observe that the velocity driven by the third - order nonadiabatic stt is about half of the wall speed driven by the first - order @xmath4 term . the inset shows the configuration of out - of - plane component of the local spin at @xmath79 for the three cases . , title="fig:",width=325 ] +    in order to study the movement of fe dw driven by current , here we consider a one dimension spin chain model , where the local spin of dw is assumed to be uniform in the plane perpendicular to the transport direction . according to the llg equation and the nonadiabatic stt modeled above , the dynamical equation for the unit vector of local spin in layer @xmath80 is @xmath81 where the extra @xmath4 term @xmath82 is only for comparision . here @xmath80 is the site index of one dimension spin chain and the distance between two sites is @xmath70 , which is unit length used in our dynamical calculations . in the meanwhile , we also use the dimensionless time where the unit time is @xmath71 . the gilbert damping coefficient is chosen as @xmath83 and @xmath84 is the gyromagnetic ratio . @xmath85 is effective field for the local spin on site @xmath80 , which is written as@xcite @xmath86 where @xmath87 is dimensionless exchange constants and @xmath88 is dimensionless anisotropy field . for the equilibrium dw without current , the dw width is determined by @xmath89[ml ] . here we take @xmath90 and @xmath91 so that the width of fe dw in our simulation is 10[ml ] . in addition , in our calculations the current density is @xmath38 $ ] , so the coefficients in stt terms based on our first - principles calculated results are @xmath92 and @xmath93 , which are measured in unit time @xmath94 and unit length @xmath70 . further , in order to demonstrate the difference of dw movement between the case of @xmath4 term and the third - order stt , we also take @xmath95 , which is chosen to make the maximum of @xmath4 term be equal to that of the third - order stt . hence , given the initial configuration of dw , the above llg equations for each unit vector of local spin can be solved numerically by runge - kutta method .    the position of fe dw center as function of time with adiabatic and two kinds of nonadiabatic stts are shown in fig . as the prediction by the previous studies,@xcite our calculated result for the case of only adiabatic stt ( @xmath96 ) shows that the dw center start to move at first and then will stop inevitably at last . this phenomenon is attributed to the existence of out - of - plane component of local spin,@xcite i.e. , the local spin will be tilted in @xmath47 direction after the current have been injected into the dw . as shown in the inset of fig . [ f7 ] , one can observe that the out - of - plane component of dw local spin moments in this case will raise up when the dw center stopped . owing to the shape of dw is film - like , the out - of - plane component of local spin moments can produce the demagnetization field @xmath97 , which can produce a torque resisting the adiabatic in - plane stt . therefore , once the demagnetization field increase to be large enough , the velocity of dw center will decrease to zero at last . the dw movements driven by the nonadiabatic stt of the third- ( @xmath98 ) and first - order ( @xmath99 ) term are also shown in fig . it can be seen that in both cases the dw will no longer stop and the velocity is nearly constant as the dw center move away from its initial position . the dw velocity induced by the third - order stt is about half of wall speed driven by the @xmath4 term . from fig . [ f7 ] our obtained wall speed caused by @xmath4 term is about 5.0[m / s ] , which is exactly equal to the value of @xmath66 . note that the dw center move backward in contrast to the case of adiabatic stt . in the inset of fig . [ f7 ] , we also show the configuration of the out - of - plane component of local spin at @xmath79 for these two cases . compared with the case of only adiabatic stt , the out - of - plane component of local spin moments with third- and first - order stt are much greater . so the precession around the relative larger demagnetization field will not only resist the adiabatic in - plane stt but also can push the dw moving backward that is in the opposite direction for the case of only adiabatic stt . moreover , as shown in fig . [ f5 ] the area under the nonadiabatic stt distribution curve for the third - order term is about half as much as that for the @xmath4 term . therefore , one can see that the out - of - plane component of local spin at dw center for the third - order stt is also half of the value for the case of @xmath4 term . as consequence , it will lead to about 50% decrease in the values of both corresponding demagnetization field and velocity for third - order stt compared to the case of @xmath4 term . in summary , we have calculated the current - induced adiabatic and nonadiabatic stts of fe , co and ni dws by first - principles noncollinear scattering wave function matching method in the frame of tb - lmto with asa approximation . we found that in co and ni dws the spatial distribution of the nonadiabatic stt are in the form of typical nonlocal oscillating torques . however , in contrast to the case of co and ni dws , the out - of - plane component of nonadiabatic stt in fe dw has localized form . the calculated results also show that the degree of nonadiabaticity in stt of fe dw is much larger than that of co and ni dws . further , we found that the distribution of localized nonadiabatic stt in fe dw can be well fitted by the term of @xmath68 , which is in the third - order spatial derivative of the local spin instead of the well - known first - order nonadiabatic @xmath4 term . the coefficient @xmath69 is also obtained for different width @xmath21 in fe dw . finally , the dynamics of local spin in fe dw driven by the third - order nonadiabatic stt is calculated using llg equation . the results show that with the same amplitude of nonadiabatic stt , the velocity of dw center driven by the third - order stt is about half of wall speed caused by the first - order @xmath4 term . the authors acknowledge prof . ke xia for suggesting the problem and dr . shuai wang for useful discussion about the calculations . the authors are also grateful to dr . yuan xu , dr . yong wang and dr . rui wang for technical assistance . we are grateful to : ilja turek for his tb - lmto - sgf layer code ; anton starikov for the tb - mto code based upon sparse matrix techniques . the authors also acknowledge the financial support from the national natural science foundation of china ( grant no:11104247 ) , china postdoctoral science foundation ( grant no:2012m520666 ) , and the provincial natural science foundation of zhejiang ( grant no : y201121807 and y13a040032 ) . m. hayashi , l. thomas , y. b. bazaliy , c. rettner , r. moriya , x. jiang and s. s. p. parkin , _ phys . * 96 * ( 2006 ) 197207 . m. yamanouchi , d. chiba , f. matsukura , t. dietl and h. ohno , _ phys . * 96 * ( 2006 ) 096601 . c. burrowes , a. p. mihai , d. ravelosona , j. v. kim , c. chappert , l. vila , a. marty , y. samson , f. garcia - sanchez , l. d. buda - prejbeanu , i. tudosa , e. e. fullerton and j. p. attane , _ nat . phys . _ * 6 * ( 2010 ) 17 . d. ilgaz , j. nievendick , l. heyne , d. backes , j. rhensius , t. a. moore , m.  . ni@xmath100o , a. locatelli , t. o. mente , a. v. schmidsfeld , a. v. bieren , s. krzyk , l. j. heyderman and m. klui , _ phys . lett . _ * 105 * ( 2010 ) 076601 . l. san emeterio alvarez , k. y. wang , s. lepadatu , s. landi , s. j. bending and c. h. marrows , _ phys . * 104 * ( 2010 ) 137205 . i. m. miron , t. moore , h. szambolics , l. d. buda - prejbeanu , s. auffret , b. rodmacq , s. pizzini , j. vogel , m. bonfim , a. schuhl and g. gaudin , _ nat . mater . _ * 10 * ( 2011 ) 419 . a. j. schellekens , a. van den brink , j. h. franken , h. j. m. swagten and b. koopmans , _ nat . commun . _ * 3 * ( 2012 ) 847 . t. koyama , k. ueda , k .- j . kim , y. yoshimura , d. chiba , k. yamada , j .- jamet , a. mougin , a. thiaville , s. mizukami , s. fukami , n. ishiwata , y. nakatani , h. kohno , k. kobayashi and t. ono , _ nat . * 7 * ( 2012 ) 635 . u. bauer , s. emori and g. s. d. beach , _ appl . lett . _ * 101 * ( 2012 ) 172403 . p. yan , x. s. wang and x. r. wang , _ phys . lett . _ * 107 * ( 2011 ) 177207 . p. chureemart , r. f. l. evans and r. w. chantrell , _ phys . b _ * 83 * ( 2011 ) 184416 . a. v. khvalkovskiy , v. cros , d. apalkov , v. nikitin , m. krounbi , k. a. zvezdin , a. anane , j. grollier and a. fert , _ phys . b _ * 87 * ( 2013 ) 020402 . g. tatara , h. kohno and j. shibata , _ physics reports _ * 468 * ( 2008 ) 213 . i. garate , k. gilmore , m. d. stiles and a. h. macdonald , _ phys . b _ * 79 * ( 2009 ) 104416 . k. gilmore , i. garate , a. h. macdonald and m. d. stiles , _ phys . b _ * 84 * ( 2011 ) 224412 . t. taniguchi , j. sato and h. imamura , _ phys . b _ * 79 * ( 2009 ) 212410 .
Abstract:

Answer: current - induced spin - transfer torques ( stts ) have been studied in fe , co and ni domain walls ( dws ) by the method based on the first - principles noncollinear calculations of scattering wave functions expanded in the tight - binding linearized muffin - tin orbital ( tb - lmto ) basis . the results show that the out - of - plane component of nonadiabatic stt in fe dw has localized form , which is in contrast to the typical nonlocal oscillating nonadiabatic torques obtained in co and ni dws . meanwhile , the degree of nonadiabaticity in stt is also much greater for fe dw . further , our results demonstrate that compared to the well - known first - order nonadiabatic stt , the torque in the third - order spatial derivative of local spin can better describe the distribution of localized nonadiabatic stt in fe dw . the dynamics of local spin driven by this third - order torques in fe dw have been investigated by the landau - lifshitz - gilbert ( llg ) equation . the calculated results show that with the same amplitude of stts the dw velocity induced by this third - order term is about half of the wall speed for the case of the first - order nonadiabatic stt .


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: cylindrical floating bodies are ubiquitous in ocean engineering . in part , this is due to the fact that the most accessible cases in the linear theory of wave - structure interaction conform to the geometries for which the field equation  laplace s equation  is separable . in three spatial dimensions , the most straightforward such geometry is the circular cylinder . due to the accessibility , and wealth of prior studies ( a brief review of which is given later in the introduction ) , the heaving truncated cylinder is an exceedingly popular model for a wave energy converter ( wec ) , able to absorb energy from waves independent of direction . due to this very ubiquity , we feel the design parameters  in this case limited to cylinder size and damping  deserve a separate , careful investigation . furthermore , we aim to extend prior work by considering all three modes of motion available to an axisymmetric body . the intention of this study is to present a rather comprehensive account of the hydrodynamics of a wec based on relative motion of two cylindrical bodies , which are allowed to move in heave , as well as sway and roll . this entails an investigation of different wec sizes and damping coefficients and their performance in wind - generated sea surface waves given by a pierson - moskowitz ( pm ) spectrum . the spectral description of the sea surface allows one to derive values for the _ significant displacements _ of the wec , and thus add a measure of survivability to the design considerations . finally , an illustrative grading system can be devised to categorize the various performance metrics of the self - reacting wecs .    following the work of falco @xcite we classify wecs into three different types based on different working principles : oscillating water columns , overtopping devices , and oscillating body devices . for the most part , offshore devices fall into the last category . furthermore , since energy production by oscillating bodies is governed by the wave - induced motion , the dimensions of such wecs are constrained significantly by the incident wavelengths , and will likely be roughly similar under given conditions . subsequently , oscillating body devices may be further subdivided into those reacting against a fixed frame of reference ( either the sea bed or a bottom fixed structure ) , and those reacting against one or more floating bodies , called _ wave - activated bodies _ or _ self - reacting devices . _ these last devices may be installed in deep water , where the large distance between the sea - bed and the surface might otherwise be prohibitive . the mooring system for such devices has the sole role of counteracting drift and current forces , allowing the weight of moorings and anchors to be relatively small ( see @xcite and references therein ) . we consider vertically floating cylinders , for which a rich literature exists . the radiation problem in heave only was addressed by ursell in 1949 @xcite , and the scattering problem by dean and ursell @xcite . miles and gilbert @xcite later employed a variational approximation to provide the far field potential for scattering by a circular dock , along with the lateral forces on the dock . however , their results were subsequently found to contain several inaccuracies , in particular in their calculations of the radiation forces . this prompted garrett@xcite to take up the problem afresh , and establish the scattering forces for a circular dock . subsequently , black et al.@xcite revisited the application of variational methods to the radiation and scattering problem by several cylindrical geometries , employing haskind s theorem to give the wave forces . this latter , variational approach did not yield the added mass and damping coefficients . hence in 1981 , yeung@xcite studied the radiation problem of a vertical cylinder floating on the water surface and undergoing the combined motions of heave , sway and roll , and obtained these hydrodynamic coefficients . more recently , bhatta@xcite also gave the added mass and damping coefficient of a vertical cylinder undergoing heave motion , in terms of the two dimensionless ratios characterizing the problem ( depth to radius and draft to radius ) . while prior work had focused on the finite depth case , in 2013 , finnegan et al @xcite treated by means of an analytical approximation due to leppington the forces on a truncated vertical cylinder in water of infinite depth .    in the context of wave energy , the consideration of floating cylinders as models of wecs goes back at least to berggren & johansson @xcite , who approximated a device described by hagerman by two floating , axisymmetric cylinders oscillating in heave , albeit without any considerations of captured power . more recently , garnaud and mei @xcite revisited the single buoy with the intention of studying it in densely packed arrays , giving the captured power for buoys hanging from a large frame . this floating point absorber was also employed , e.g.  by child and venugopal @xcite in their discussion of optimization of wec arrays , or by borgarino et al @xcite as a generic model to investigate wave interaction effects . similarly , teillant et al @xcite employ an axisymmetric , heaving 2-body device for their study of wec economics , without detailed hydrodynamic considerations . meanwhile , a slightly different construction of the floating body was considered by engstrm et al @xcite , who added a sphere under the floating cylinder . this two - body configuration of floating cylinder and submerged sphere was then assumed connected to the sea bed by a generator , and its performance analyzed . zheng et al @xcite , in a generalization of berggren & johansson to three modes of motion , considered the hydrodynamics of two unconnected , coaxial floating cylinders , again without considering power capture . the power capture for a heave - only two - cylinder wec was recently obtained for attacking monochromatic incident waves by wu et al @xcite , albeit with a rather terse discussion of their results . we combine features of several previous authors , and consider the novel case of two floating cylinders with an idealized power take  off ( pto ) , represented by a linear damper of constant characteristics , whose optimization is part of the design procedure .. we undertake our parametric study with an eye towards applications , and thus consider irregular waves in the form of a pm spectrum ( see e.g.  recent work on optimizing a floating box - barge under irregular waves by  bdai & srinil @xcite ) while scatter diagrams may be available for some sites where an assessment of the wave resource has been carried out , where this is not the case estimates based on wind speed will need to be made . to this end , we present our data nondimensionalized on the basis of wind speed , which uniquely determines the pm spectrum . values of significant wave height and peak period may be readily derived therefrom , and the data recast in these terms if desired .    the paper is organized as follows : in section [ sec : physical preliminaries ] we present the physical set - up of the problem . this consists in presenting the twin cylinder wec and characterizing its geometry , and subsequently presenting the pm spectra for design and survivability considerations . in section [ sec : governing equations ] we present , very briefly , the basic mathematical formulation of the governing equations and sketch the solution procedure . subsequently , we employ the hydrodynamic coefficients and forces found from solving the equations of section [ sec : governing equations ] to characterizing wec design under monochromatic waves in section [ sec : design of the wec ] , and under irregular waves given by a pierson - moskowitz spectrum in section [ section : irregular design ] . a discussion of these results with a view to applications is given in section [ section : discussion ] , which is subdivided into section [ subsec : power capture ] on power capture , section [ subsec : survivability ] on survivability , and section [ subsec : grading ] which presents a synthesis of the preceding sections . finally , section [ sec : conclusions ] presents some concluding remarks and perspectives . the geometry and basic parameters of the twin - cylinder wec are depicted in fig [ fig : geometry ] . the @xmath0 plane is the still water surface and the @xmath1-axis points upwards . ( @xmath2 ) are polar coordinates in the horizontal plane , such that @xmath3 and @xmath4 the upper cylinder floats on the water surface with a draft @xmath5 . to provide for flotational stability , it is important to note that the mass of this cylinder is not uniformly distributed , but is divided into two parts with drafts @xmath6 and @xmath7 and densities @xmath8 and @xmath9 , respectively . the lower cylinder is entirely submerged with a draft @xmath10 , and like the upper cylinder is assumed divided into two parts with densities @xmath11 and @xmath12 and drafts @xmath13 and @xmath14 , respectively . the distance between the two cylinders is @xmath15 . both of them have the same radius @xmath16 , and the water depth @xmath17 is taken to be very large compared to the attacking wave length , with the intention of approximating deep - water conditions .    as shown in fig.[fig : geometry ] , the two cylinders are connected by a continuously distributed dashpot , which connects the upper edge of the lower cylinder with the lower edge of the upper cylinder at @xmath18 . the integrated dashpot coefficient is @xmath19 , which results in a dashpot coefficient per length @xmath20 . the dashpot is considered to represent a pto , which generates energy from both the relative heave and roll motion of the cylinders . twincylinders_new(0.2 )    since the two cylinders are axisymmetric , only these three modes are studied . the heave and sway motions will give rise to relative motions in @xmath1 and @xmath21 directions , respectively . for waves propagating in the @xmath22direction , the two cylinders roll around the @xmath23-axis in the mean free surface ( @xmath24 ) , yielding a relative angle about this axis . this formulation of the problem leaves us with thirteen parameters ( @xmath25 ) characterizing the wec . before proceeding , we will make several restrictions to ensure that the problem remains manageable ; nevertheless , we shall see that a wealth of interesting phenomena and properties of the wec are still accounted for .    for simplicity , we will take the drafts and distance between the cylinders identical to their radius , and denote the single size parameter by @xmath26 , i.e. @xmath27 for the density distribution of the cylinders , we shall assume @xmath28 where @xmath29 is the density of the water . thus , the design problem is reduced to two parameters , a size @xmath26 and dashpot coefficient @xmath19 , whose interplay with incoming waves of certain frequencies is the issue at hand . we shall see that suspending the lower cylinder at a depth @xmath30 below the still water surface has the desired effect of rendering its motion rather small , and thus creating a relatively stable point for the upper cylinder to react against . there are several reasonable criteria which may govern the design of a wec . evidently , the wec should capture as much of the incoming wave energy as possible . at the same time , as economic viability is the prime driver behind wave energy technology , the costs  which may be assumed to scale with device size  should be kept low ; in practical terms , this means that device size should be kept small . competing with this are concerns over the survivability of the converter , which dictate that displacements of the wec not be too large under severe conditions , favoring larger devices . we shall return to these issues in detail in later sections . one of the most common descriptions of a sea - state for engineering purposes is the unidirectional pierson moskowitz ( pm ) spectrum , here given as a function of wavenumber @xmath31 @xmath32 where @xmath33 is the mean wind speed at a height of @xmath34 m above the mean surface level , and @xmath35 is the gravitational acceleration . this empirically derived formula gives the energy distribution for wind waves in deep water , and differs from the jonswap spectrum only by the addition of a spectral  peak enhancement factor . this spectrum ( [ pm2 ] ) readily yields a number of important values associated with the sea - state : @xmath36 where @xmath37 is the significant wave height and @xmath38 is the wave number of the spectral peak for a given wind speed @xmath33 . this makes it easy to present subsequent results in an alternative form when desired . a monochromatic wave with wavenumber @xmath38 and the same wave energy density as the pm spectrum will have an amplitude @xmath39    for subsequent illustration it will be necessary to have some concrete , physical examples , which means specifying a sea - state via a wind speed value @xmath40 our design conditions ( denoted by a subscript @xmath41 ) will correspond to a wind speed @xmath42 m / s , while we will consider two  severe states " ( denoted by subscripts @xmath43 and @xmath44 ) with regard to the survivability , corresponding to wind speeds @xmath45 m / s and @xmath46 m / s . these are summarized in table [ swss ] . @ l c c c @ & @xmath47&@xmath48&@xmath49 +    @xmath50 & 10&15&20 + @xmath51&2.47 & 5.55&9.87 + @xmath52 & 0.065 & 0.029 & 0.016 + @xmath53 & 96.30 & 216.67 & 385.19 +    [ swss ] our approach to solving the wave - structure problem for the twin - cylinder wec relies on domain decomposition and eigenfunction expansion methods ( in the context of floating cylinders , see black & mei @xcite , who give a comprehensive description of the method , or more generally , linton & mciver @xcite , or zheng et al @xcite for a recent application to floating cylinders ) . as the full formulation is rather lengthy , we only indicate the most important equations , and refer the interested reader to work cited above .    the fluid is assumed to be incompressible and inviscid , and the flow irrotational . introducing a velocity potential @xmath54 and assuming periodic motion of frequency @xmath55 the potential is separated into the spatial and temporal parts , @xmath56 where @xmath57 satisfies the laplace equation : @xmath58 subject to the linearized boundary conditions on the free surface @xmath59 and on the bed @xmath60 : @xmath61 where @xmath62    at the interface between structure and fluid , the normal velocity of the structure must equal that of the adjacent fluid particles , written in terms of the potential : @xmath63 where @xmath64 is the component of the structure s velocity in the direction of the outward pointing normal vector @xmath65 which may be applied at the equilibrium surface under the assumptions of linearity . owing to this very linearity , we continue with a decomposition of the problem into two parts : one due to the waves ( @xmath66 ) scattered from the structure ( which is assumed fixed ) by the incident wave field , and one due to the waves ( @xmath67 ) radiated by the motion of the structure , such that @xmath68 @xmath66 is decomposed further into the potential due to the incident wave @xmath69 and that due to the waves diffracted from the fixed structure @xmath70 where @xmath71 the remaining radiated part of the potential @xmath67 must then satisfy , where the normal velocities are to be determined from the equations of motion . we shall consider an incident monochromatic wave with amplitude @xmath72 so that @xmath69 is known _ a priori . _    introducing the as - yet unknown displacements of the upper ( @xmath73 ) and lower ( @xmath74 ) cylinder for the three modes of motion @xmath75 where @xmath76 , @xmath77 and @xmath78 are the complex amplitudes of the corresponding displacements , we can write the boundary condition for the spatial part of the total potential @xmath79 in the following form @xmath80 where is posed on the bottom of the upper cylinder , on the top and bottom of the lower cylinder , on the sides of the upper cylinder , and on the sides of the lower cylinder . these conditions are supplemented by sommerfeld s radiation condition , requiring the diffracted and radiated waves to be outgoing as @xmath81    subdomains_new(0.2 )    due to the configuration of two axisymmetric floating cylinders , we must consider three fluid regions , one between the coaxial cylinders ( region ii ) , one between the lower cylinder and the bed ( region iii ) , and one outside the vertical extension of the cylinders ( region i ) , as depicted in figure [ fig : subdomains ] . subsequently the scattering problem is divided into three problems , one in each subdomain , and the radiation problem for each of the three modes of each of the two cylinders is divided into three problems . the reader may appreciate the effort involved in keeping track of , solving , and subsequently matching solutions , of 21 problems for the potentials involved . these potentials are then applied in calculating the forces on the two cylinders , in the form of pressures from the surrounding fluid . the full expressions for the exciting , hydrodynamic , and hydrostatic forces are lengthy and will not be given . we note only that we have found excellent agreement between our results and published work @xcite . the forces due to the fluid , as well as the forces due to the dampers are employed with newton s second law to yield the body motions . the first two - equate the vertical ( heave ) forces with the masses and accelerations of the upper and lower cylinder , respectively . the next - are those for the horizontal ( sway ) forces . the final pair - equate the torques about the @xmath82axis to the angular acceleration times moment of inertia of the upper and lower cylinder , respectively . @xmath83 the terms appearing in the above equations have the following interpretation :    @ l m9 cm @ @xmath84 & exciting forces / torques on cylinder @xmath85 in the @xmath86 and @xmath1 directions + @xmath87 & hydrodynamic force / torque of the @xmath88 motion of cylinder @xmath85 in the @xmath89 direction of cylinder @xmath90 + @xmath91 & hydrostatic forces in the @xmath23 and @xmath1 direction on cylinder @xmath92 + @xmath93 & damping forces in the @xmath23 and @xmath1 direction on cylinder @xmath92 + @xmath94 & displacement amplitudes of cylinder @xmath85 in sway ( @xmath95 ) , heave ( @xmath96 ) , and roll ( @xmath97 ) . + @xmath98 & mass of cylinder @xmath92 + @xmath99 & moment of inertia of cylinder @xmath85 about the @xmath23-axis . the masses and moments of inertia have the explicit form ( see , ) @xmath100 and the damping forces are @xmath101 after the displacements of the cylinders are obtained , the captured power can then be calculated as follows : @xmath102 where @xmath103 and @xmath104 are the complex conjugates of @xmath76 and @xmath78 , respectively . we now undertake to examine the design of the wec , based on the three parameters characterizing the environmental conditions @xmath105 and @xmath106 the two wec parameters @xmath26 and @xmath107 and the seven wec performance parameters calculated from the wave - structure interaction problem @xmath108 and @xmath109 an application of buckingham s @xmath110 theorem @xcite yields that there will be nine dimensionless quantities that characterize this problem : @xmath111 , @xmath112 , @xmath113 , @xmath114 , @xmath115 , and @xmath116 in the sequel , we will make use of a @xmath117 to denote nondimensional variables , i.e. , the nine dimensionless quantitites above will be @xmath118 and @xmath119      for simplicity of presentation and ease of understanding we initially consider only the heave mode , motivated by the fact that , while sway and roll are generally coupled , they are both independent of heave . the response of the wec under incoming monochromatic waves is first considered , where our physical test - case corresponds to a monochromatic wave of wavelength @xmath120 m equal to that at the peak of the design spectrum @xmath121 and an amplitude @xmath122 such that the total energy density of the wave is equal to that of @xmath121 see and table [ swss ] . we first choose the dashpot coefficient @xmath19 to be zero , which means that the two cylinders are freely floating . in this case , once the incident monochromatic wave is given , the only wec parameter to be determined is @xmath26 . dimensional analysis can then be applied to the problem of determining the quantity of interest @xmath26 for the motions of the upper cylinder @xmath123 and the lower cylinder @xmath124 separately . once again , buckingham s @xmath110 theorem yields that , for the variables @xmath125 and @xmath126 there exist exactly two nondimensional quantities , which must be related by a relation @xmath127 the maximum displacement of the cylinder @xmath128 as a function of size @xmath129 thus corresponds to the extrema of @xmath130 equation is plotted in figure [ fig : free - heave - z ] for the upper and lower cylinders .    ) under the design monochromatic wave . @xmath131 : upper cylinder , thick line ; @xmath132 : lower cylinder , thin line . ] as we are ultimately interested in relative displacements of the cylinders , the global maximum of @xmath133 and the local maximum of @xmath134 which occur at @xmath135 yield the chosen design size . the maximum displacement in fig .  [ fig : free - heave - z ] is related to the resonance between the cylinders and the incident monochromatic wave . the introduction of a damper , while changing the magnitude of the displacement , can be shown to have no effect on the location of the resonant maximum , which remains @xmath135 ( see fig .  [ fig : free - heave - z ] ) even for various values of @xmath19 . thus , the size of the wec determined from the freely floating case is used to specify the dashpot coefficient @xmath19 . given the unique relationship between @xmath26 and @xmath126 independent of @xmath107 described above , the dimensional analysis of the power absorption involves the quantities @xmath136 and @xmath137 which can be written in three dimensionless ratios : @xmath138 where @xmath139 is plotted in fig . [ fig : pald - heave ] for the wec size as determined in the last section ( @xmath135 ) .     and the dashpot coefficient @xmath140 for the heave motion induced by the design monochromatic wave , where @xmath135 . ] we elect to determine the dashpot coefficient @xmath19 from the maximum of captured power @xmath141 in figure [ fig : pald - heave ] , calculated from the heave terms only in ( [ pa ] ) . this results in @xmath142 and @xmath143 . thus , the wec design for a monochromatic wave has been determined . taking the design wave introduced in the beginning of section [ sec : design of the wec ] as a physical example , the wec has the dimensions @xmath144 and @xmath145 , and can capture @xmath146 from a monochromatic wave 96.3 m long and 0.87 m in amplitude . having treated the simpler case of heave - only motion , we now consider the general case in which the wec is additionally allowed to undergo sway and roll motions . akin to the previous section which only dealt with the heave motion , the design procedure of the wec in the general case is also divided into two steps , as illustrated in detail below . we start with the freely floating case , where the dashpot coefficient @xmath147 . using the equations of motion ( [ eq : mo1])-([eq : mo6 ] ) , we can obtain the displacements of the two cylinders in the @xmath21 and @xmath1 directions , and the angle around the @xmath23 axis . once the monochromatic wave is given , or equivalently , once the mean wind speed for the corresponding pm spectrum is given , the physical process of determining the size of the wec can be written in the following dimensionless form : @xmath148 where @xmath76 and @xmath77 denote the amplitudes of the vertical and horizontal displacements respectively , @xmath78 is the amplitude of the angle around the @xmath23 axis , and @xmath149 corresponds to the upper and lower cylinder , respectively . we now seek the maxima of the functions @xmath150 and @xmath151 presented in fig.[freefloat-3 m ] . due to the increase in number of modes , the picture of the displacements is more complex than in the preceding section . it may be observed that the heave mode is decoupled from the sway and roll modes , yielding again the global maximum at @xmath152 the sway and roll modes are coupled , and are observed to present a global maximum for relative displacement at @xmath153 resulting in an ambiguous situation for determining the size of the wec .      as in the preceding section , we now suppose that the size of the wec is given . the captured power @xmath141 then depends on the dashpot coefficient @xmath19 . the determination of optimal power absorption as a function of dashpot coefficient is described in dimensionless form by @xmath154 where , as we have seen , there is some flexibility in choice of @xmath155 the function @xmath156 is plotted in fig . [ fig : pald-3 m ] for both @xmath157 and @xmath158 . for the device operating optimally in heave ( @xmath159 thick line ) there is a unique maximum at @xmath160 with @xmath161 ( denoted case e ) , very close to the heave - only case discussed in section [ subsec : heave only ] . for the roll  sway optimized device ( @xmath162 thin line ) there are two local maxima @xmath163 and @xmath164 with corresponding @xmath165 and @xmath166 ( denoted case a1 and a2 ) respectively .     of the wec in the combined motions versus the dashpot coefficient @xmath19 . thin line : @xmath157 ; thick line : @xmath158 . ] the situation for monochromatic incident waves is summed up in table [ 3m : mono : ah ] , which shows the nondimensional size , optimal damping , captured power , and displacement amplitudes for the cases discussed above . as we have observed , introducing roll and sway motions leads to a two - fold branching in the design procedure . firstly , in free motion , one value of @xmath129 is found to yield the largest roll and sway displacements , while another value yields the largest heave displacements . while the heave - optimized case has a unique maximum @xmath167 as a function of damping , the roll / sway - optimized case admits two local maxima of @xmath168 one with relatively low , the other with relatively high damping @xmath169 compared to the heave case ( see figure [ fig : pald-3 m ] ) . this opens up the possibility that the overall optimal design may not coincide with a design optimized for roll / sway or heave alone , but occupying some middle ground . the performance of such intermediate devices ( cases b , c , and d ) , as well as devices somewhat larger than case e ( cases f , g and h ) is explored for the monochromatic design wave in table [ 3m : mono : ah ] . in each of cases a1 through h , a damping @xmath170 has been chosen to maximize the captured power . .the size , damping , displacement and captured power of 3-mode wecs in monochromatic waves . [ cols="<,<,<,<,<,<,<,<,<,<",options="header " , ]     we have developed an example framework for survivability for the twin - cylinder wec in the three spectral sea states considered , which is presented in table [ table : survivability ] . the maximum allowed relative vertical travel @xmath171 is limited to @xmath172 while the maximum allowable roll is @xmath173 those cases which exceed these values are marked red ( pdf only ) . a vertical travel of more than @xmath174 or a roll of more than @xmath175 is marked orange , while a vertical travel of more than @xmath176 or a roll more than @xmath177 is marked yellow ( pdf only ) . device motions smaller than these are marked green ( pdf only ) . recall that these nondimensional quantities depend only on the relations @xmath178 and @xmath179 as specified in section [ sec : physical preliminaries ] , and the concomitant changes in significant wave - height and peak wavenumber . for illustrative purposes , if the design spectrum @xmath47 is generated by a fresh breeze ( @xmath180 or 5 beafort , 2.47 m @xmath181 ) , then the first severe state @xmath48 may be thought generated by a high wind ( 7 beaufort , 5.5 m @xmath181 ) . the second severe state @xmath49 occurs under conditions between gale and severe gale ( 89 beaufort , 9.9 m @xmath181 ) . these extremely harsh conditions represent an energy density more than 16 times that of the design spectrum , and may be expected to challenge the device design . while the increase in significant wave - height between the design spectrum @xmath47 and the severe case @xmath49 may seem dramatic , there is no doubt that such conditions will be encountered within the operational life of a wec . for example , while deep water conditions for the eastern mediterranean off israel s coasts may see significant wave heights greater than 2 m only 6 % of the time , and wave heights in summer rarely exceed 11.5 m , nevertheless storms with @xmath181 in excess of 5 m occur almost yearly . the 10-year return period significant wave height is nearly 7 m , which clearly falls within the expected operational life of a converter .    from a pure survivability standpoint , it is immediate only that the smallest converter a1 is not viable . in particular , the very small damping of this configuration ( see table [ table : dimensional power capture ] ) , while allowing for efficient power capture from the roll mode , also leads to overly large displacements even for design conditions . with survivability as the central aim of design , larger structures will necessarily fare better , though the differences between devices d , e , and f are in practice rather small . while other authors ( e.g.  maisondieu @xcite or brown et al @xcite ) have investigated survivability of wecs , they have been forced to do so without reference to the hydrodynamics and actual displacements of a floating device , but rather purely based on estimations of the incident wave power . we shall now make a preliminary attempt to sum up the results of the preceding sections . the intricacies of wec economics , as well as the many factors which are outside the scope of the present study , such as moorings , specifics of the pto , control strategies , power conversion and transmission , and other environmental factors from seasonal variability to extreme events , will need to be taken into account for a fuller analysis . in addition , wec cost will not be considered , and is likely to impact significantly the ultimate design . ( solid line ) , the first severe spectrum @xmath48 ( dashed line ) , and the second severe spectrum @xmath49 ( dash - dotted line ) , where @xmath42 m / s . ]    figure [ fig : comparison abs power ] compares the dimensional absorbed power for the nine devices considered under different conditions . once again , we point out that the monochromatic case is essentially artificial , and included only for illustrative purposes . the lessons to be drawn from this comparison will likely change as wave - power technology matures . while current oscillating - body devices may be rather small , and situated in shallow water with the intention of keeping costs down , future developments will naturally lead to a move into the more powerful wave - regimes further offshore ( see stiassnie et al @xcite for a discussion ) .    as an example , while there is a 15 % reduction in absorbed power between case e and case h under the design spectrum , the corresponding increase in absorbed power for severe case @xmath49 is upwards of 60 % . the fact that , off the eastern mediterranean coast , some 45 % of average wave power comes during storm events that occur only 5 % of the time indicates the utility of the larger design @xcite . this is compounded by the increase in potential survivability of the larger devices as indicated in the previous section . depending on the variability of the wave - energy resource , more or less weight may ultimately be given to each of the considerations just outlined . the fact that the larger devices exhibit smaller relative motions may also be a benefit for their reliability , in terms of limiting loading during normal operation . ultimately , an effort will have to be made to weigh the additional cost of a larger device against the increase in survivability . both of these in turn will need to be weighed against the potential of continuing operation during high - energy events , while sustaining a slight performance decrease for low - energy sea states . we have investigated in detail the hydrodynamics of a model wec consisting of two floating , axisymmetric cylinders connected at their upper and lower perimeters by a continuously distributed damper  allowing power capture from heave and roll modes . while other authors have studied various aspects of the problem of floating cylinders , the present work has addresses for the first time the problem of a twin cylinder wec allowed to move in three degrees of freedom . the inclusion of a floating , submerged cylinder as a mechanical reference for power extraction makes this design viable in deep water . with further development of the wave energy industry , it may be anticipated that wecs will follow wind turbines in moving further offshore , making such self - reacting devices more and more relevant @xcite . our design procedure initially focused on optimizing device behavior for a damper of constant characteristic in monochromatic waves . at the outset , the heave - only case was considered , presenting a simple situation where a single device ( characterized by a size parameter @xmath26 and a damping parameter @xmath19 ) , coinciding with the resonant maximum of a freely floating body , outperformed all others . allowing the device also to sway and roll was seen to introduce additional complexity , and a differentiation was observed between devices operating preferentially in roll / sway and those operating preferentially in heave .    despite the multiplicity of possible designs when the device is allowed to undergo heave , sway , and roll motions , the monochromatic case presents a clear picture from the standpoint of power absorption : the device closest to heave resonance is found to perform best . this conclusion is an artifact of the idealization represented by the monochromatic theory  a fact established by the subsequent investigation of wec performance under an irregular sea .    for our design purposes , a pierson - moskowitz spectrum , characterized by wind speed , was chosen to evaluate the designs obtained from the monochromatic case . under this spectrum , the maxima of absorbed power were found to shift markedly with respect to the monochromatic case , reflecting the potential for detuning to increase captured power in wec design . larger values of absorbed power under the design spectrum were found for devices slightly larger and slightly smaller than the monochromatic optimum , raising the question of how to determine device sizing in light of other criteria .    to this end , we have devised some example metrics for grading the sizes of our twin - cylinder wec . we note that wave energy presents particular difficulties in many respects . while a fixed offshore structure may be designed for survival with very high safety factors , this is inappropriate for oscillating body wecs ; by their nature , they must undergo the largest possible motions in order to extract energy . at the same time , device loading should be minimized to avoid fatigue and failure . taking into account the fact that wecs may be expected to be operational for on the order of 25 years ( see starling @xcite ) , and it becomes clear that survival is a paramount issue . we have presented an example approach to quantitatively evaluate the competing aims of survivability and power extraction within the framework of our floating twin - cylinder device .    to a certain extent all renewable energy technologies , wecs more than most , can not control their operating conditions , but must work within their environment , subject to the resulting fluctuations of the resource . it must be expected that , like wind turbines , oscillating body wecs will be designed with a  survival mode " , when normal operation cease , and the device changes its characteristics in order to avoid extreme loads . ( we might note that overtopping wecs or owcs ( see section [ sec : introduction ] ) , due to a different working principle and resulting size , will likely have a very different survivability analysis than oscillating body designs . ) this may mean increasing the damping , altering the water plane area or mass ( see stallard et al @xcite ) , or other approaches ( see coe and neary @xcite ) . due to the nascent state of commercial wave - energy technology , it is difficult to offer concrete design recommendations based on the results for floating twin - cylinders . our discussion does bear out the fact that a slight over - engineering may be preferable , given the large relative contribution of infrequent , high - energy events to the annual energy budget at many sites , and the demands of survival and robustness . we believe these results to be applicable more broadly to oscillating - body converters , constrained in size as they are by the incident wavelength , indicated by the striking similarities in performance between our twin - cylinder configuration and a single bottom - referenced cylinder . while many topics of critical importance have not been touched upon  from moorings to electricity transmission to specifics of power take - off  we hope to have presented a rather comprehensive picture of the design of a floating twin - cylinder wec in deep water . there is considerable room for future work , e.g taking into account seasonal resource variability , extreme wave statistics , pto control , and device cost , as well as detailed studies of reliability under regular operation , where it is hoped that some of the gaps left by this work will be filled in . this research was supported by the israel science foundation ( grant 464/13 ) .                    adam brown , robert paasch , irem  y tumer , pukha lenee - bluhm , justin hovland , annette von jouanne , and ted brekken . . in _ asme 2010 4th int . energy sustain . _ , pages 917927 . american society of mechanical engineers , 2010 .
Abstract:

Answer: we discuss the hydrodynamics of a wave energy converter consisting of two vertically floating , coaxial cylinders connected by dampers and allowed to heave , sway and roll . this design , viable in deep water and able to extract energy independent of the incident wave direction , is examined for monochromatic waves as well as broad - banded seas described by a pierson moskowitz spectrum . several possible device sizes are considered , and their performance is investigated for a design spectrum , as well as for more severe sea states , with a view towards survivability of the converters . in terms of device motions and captured power , a quantitative assessment of converter design as it relates to survival and operation is provided . most results are given in dimensionless form to allow for a wide range of applications .