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: axis is antero - posterior : anterior ( head ) towards the left , posterior ( abdomen ) towards the right . the @xmath0 axis is medio - lateral . crosses : approximate positions of severing , along the midline in the scutellum ( blue ) and off - axis in the scutum ( black ) . + ( b ) epithelial cell apical junctions marked by e - cadherin : gfp just before the severing , @xmath1 , here in an old pupa ( see text for classification ) . blue circles : two concentric circles define the annular severed region ; the distance between circles corresponds to about one cell size . yellow arrows : macrochaete used as spatial references to position the severed region . + ( c ) first image after severing , @xmath2  s. yellow : fitted ellipse @xcite ( see methods ) . ( d ) time @xmath3  s after the severing , showing a larger opening along @xmath0 than @xmath4 . bars  @xmath5 mm ( a ) , @xmath6 m ( b - d ) . [ fig : metamorphosis ] [ fig : fit_ellipse ] ]    an epithelial tissue is a sheet of cells that acts as a barrier , separating for instance the outside and the inside of a multicellular organism . its biological function relies in part on the formation of a network of adherens junction belts , connected to the acto - myosin cytoskeleton , where cells adhere to each other @xcite , and which transmits mechanical information over several cell diameters  @xcite . a key issue is to understand and model the role of tissue mechanics ( forces , displacements , and their time evolution ) in the coordinated changes of cell shape and position which determine morphogenetic flows at the tissue level @xcite . several models describe tissues using _ continuum mechanics _ @xcite . one precondition is the existence of a mesoscopic scale defining a domain over which averages of cell properties are well - defined @xcite . this description further relies on the assumption that the tissue mechanical state can be quantified , at the same mesoscopic scale , by two variables @xcite : the stress characterizing in which directions , and to what extent , the domain is under tension or under compression ; and the strain characterizing how far the domain geometry is from that of a relaxed state . such description has the advantage of accounting for individual cell properties using only a few spatially averaged parameters , which determine for instance how fast and through which succession of states the domain reacts to external solicitations .    beyond simple observation , _ in vivo _ mechanical measurement techniques include elastography @xcite , photoelasticity @xcite , magnetic micromanipulation @xcite , tonometry @xcite or nanoindentation @xcite . a large literature ( for review , see _ e.g. _ @xcite ) has established laser ablation of individual cell junctions as a tool to measure the tensions within an epithelium , in particular during _ drosophila _ dorsal closure @xcite . this technique has allowed to measure the material relaxation time @xmath7 and contributed to a better characterization of morphogenetic processes in _ @xcite . the analysis of single cell junction ablation is usually based on reasonable assumptions @xcite : that the tissue is at mechanical equilibrium just before severing ; that the ablation is effective in removing at least part of the cell junction tension ; and that during relaxation the velocity remains small ( reynolds number much smaller than 1 ) . within these assumptions , the initial retraction velocity yields the value of the cell junction tension removed by the ablation , up to an unknown prefactor which depends on the dissipation @xcite . the tensor @xmath8 denotes the component of the stress removed by ablation . it arises as an average over several individual cell junctions within a region of the tissue . by ablating straight lines in either of two perpendicular directions @xcite or by performing statistics on single cell junction ablations in several directions @xcite , it is possible to measure step by step the anisotropy of the ratio @xmath9 ( where @xmath10 denotes the tissue viscosity ) . the stress to viscosity ratio @xmath9 is the initial strain rate after the severing : the rate of displacement that the epithelium would spontaneously undergo if it were free . the continuous description applies to other cellular materials such as foams  @xcite . in the present work , we investigate experimentally whether it also applies to epithelial tissues . more specifically , we ask whether the epithelium strain , stress and material properties can be directly measured on a mesoscopic scale . we address these questions in a model system , the dorsal thorax epithelium of _ drosophila _ pupa . the pupa is the life stage during which the larva starves , feeds upon its reserves and metamorphoses into an adult @xcite . the dorsal thorax is a single layer of cells whose apical surfaces , surrounded by the adherens junction belts , face the cuticle , which protects the pupa ; their basal surfaces face the hemolymph , which acts as a nutrient transporter . it is composed of a large anterior region , the scutum , and a posterior tip , the scutellum . the thorax , and especially the scutellum , undergoes extensive morphogenetic changes during the metamorphosis that shapes the adult fly ( figure  [ fig : metamorphosis ] ( a ) ) , making it a useful model to decipher the mechanisms that control tissue morphogenesis during development @xcite . it has a bilateral symmetry axis : the midline ( axis @xmath4 , dashed line in figure  [ fig : fit_ellipse ] ( a - b ) ) . it displays large cells , the macrochaete ( precursor cells of adult sensory hairs , figure  [ fig : metamorphosis ] ( b ) ) , whose positions are precisely reproducible @xcite . we introduce an original type of severing experiment . we sever by short pulse laser the adherens junctions in an annular region around a @xmath1130  @xmath12 m radius circular tissue domain ( see figure  [ fig : fit_ellipse ] ( b - d ) , methods and movies 1 - 3 ) . the inner and outer wound margins display a comparable speed , displacement and anisotropy . note that the same anisotropy in speed and displacement results in the outer ellipse boundary having minor and major axes orthogonal to the minor and major axes of the inner one . after severing , the wound heals within tens of minutes and the pupa develops to adult . experiments are done in the scutellum ( blue cross in figure  [ fig : metamorphosis ] ( a ) ) . we classify them into three groups , according to the developmental age at which the severing is performed : earlier than 15 - 16  h , approximately 18 - 20  h , and later than 24 - 26  h after pupa formation . for brevity we call them  young " ( @xmath13 ) ,  middle - aged " ( @xmath14 ) and  old " ( @xmath15 ) pupae , respectively . similar experiments can be performed in other locations . as a proof of principle , we include three experiments performed in the scutum ( black cross in figure  [ fig : metamorphosis ] ( a ) ) of middle - aged ( @xmath16 ) and old ( @xmath17 ) pupae . we analyse the retraction of the inner tissue domain and its margin , within the set of assumptions relevant for the retraction of a single severed cell - junction . the boundary of the retracting domain is fitted at each time point by an ellipse @xcite ( see methods and figure  [ fig : fit_ellipse ] ( c - d ) ) , yielding values of the @xmath4 and @xmath0 axes as a function of time ( figure  [ fig : tau_x_vs_v ] ) . our protocol offers several advantages . ( i ) the experiment directly yields measurements averaged over the severed tissue domain , comprising typically a hundred of cells . as discussed below , the domain size is thus an adequate mesoscopic scale , crossing - over from detailed cell - level to tissue - level continuous descriptions @xcite . ( ii ) the comparison between the initial and final states yields a direct measurement of the strain which existed before severing . ( iii ) several cell - cell junctions in various directions are severed at once ( rather than step by step ) and in the same pupa . hence one experiment yields a direct measurement of all stress to viscosity tensor components . it takes only a few tens of seconds , the typical value of the relaxation time . ( iv ) last but not least , by isolating in time and space an epithelium domain from its neighboring cells , we obtain a complete set of well - defined initial , boundary and final conditions for a spatio - temporal model of the retracting domain . we thus introduce a continuous model , in the spirit of mayer _ et al . _ @xcite and taking into account a coupling between space and time dependences . as explained below , fitting such a model to the displacements in the domain bulk tests the continuous description , suggests an interpretation of the relaxation time , and probes the competition between internal viscosity and external friction . model - independent measurement of strain @xmath18 , stress to viscosity ratio @xmath9 , and relaxation time @xmath7 . the wound margin position ( ellipse semi - axis ) is plotted _ versus _ time after severing ; data from figure  [ fig : fit_ellipse ] ( b - d ) , esm movie 1 , along the @xmath0 axis . the difference between the initial ( @xmath19 ) and final ( @xmath20 ) positions ( blue dashed lines ) directly yields the value of @xmath21 . the initial velocity @xmath22 ( slope of the orange dashed line ) yields an estimate of @xmath23 . inset : velocity , estimated by finite differences of successive positions , _ versus _ the position during the first @xmath24 s. an arrow indicates the direction of increasing time @xmath25 . the slope of a linear fit ( purple line ) yields the inverse of the relaxation time , @xmath26 ( see equation  [ eq : measure_tau ] ) . ] severing the epithelium reveals the displacement that the epithelium would spontaneously undergo if it were free . we measure the strain using a logarithmic definition @xcite , hencky s  true strain " @xcite : @xmath27 where @xmath19 is the initial radius , and @xmath28 is the ellipse semi - axis at the end of the relaxation . similarly @xmath29 . the advantages of this definition are that @xcite : it is valid at all amplitudes ; it is adapted for tensorial measurements ; and finally it increases the range of validity of the linear elasticity approximation which applies here even to the highest strain value we measure , @xmath30 . each severing is followed by a relaxation to a final domain strictly smaller than the initial disk ( figure  [ fig : fit_ellipse ] ( d ) ) , indicating that before the severing the tissue had a positive strain in all directions . since the midline is a symmetry axis , we expect that the strain axes are parallel and perpendicular to it . we check that this is the case for the fitted ellipse axes , and that accordingly the shear strain @xmath31 is indistinguishable from 0 . we thus plot @xmath32 and @xmath33 ( figure  [ fig : strain_anisotropy ] ( a ) ) , obtained with absolute precision better than @xmath34 . the values are clustered according to the three pupa age groups : @xmath18 is low and isotropic at young age ( green cluster ) , moderate and isotropic at middle age ( red ) , high and anisotropic ( @xmath35 ) at old age ( blue ) . similarly , in each experiment , the initial retraction of the wound margins ( figure  [ fig : fit_ellipse ] ( c ) ) indicates the sign of the stress : before the severing the tissue was under tensile stress in all directions . this is reminiscent of positive cell junction tensions observed in the wing @xcite and in the notum @xcite . the initial retraction velocity divided by @xmath19 yields the stress to viscosity ratio @xmath9 ( figure  [ fig : tau_x_vs_v ] ) , in a way which is likely to be independent of any tissue rheological model , as suggested by the following argument . when ablating a single cell - cell junction , force balance shows that the force @xmath36 exerted by the ablated junction on the neighbouring vertex is proportional to the initial recoil velocity @xmath37 : @xmath38 , where @xmath39 is a friction coefficient @xcite . coarse - graining this relation , and assuming that the main source of dissipation is identical in our experiments which sever many junctions at once , yields @xmath40 , where the factor @xmath19 is introduced by dimensional analysis . results for the stress to viscosity ratio @xmath9 ( figure  [ fig : sigma_anisotropy ] ( b ) ) are similar to those for @xmath18 : its amplitude and anisotropy are initially small and increase with age ; the values are clustered according to the three age groups ; no shear stress is detected . values of @xmath9 span almost two decades , reflecting the sensitivity and precision of the method . the relaxation time , @xmath7 , is obtained from the time evolution of the ellipse axes . instead of fitting an exponential to the data , we plot the velocity _ versus _ position ( inset of figure  [ fig : tau_x_vs_v ] ) . this method is robust and independent of any assumption . a linear regression of the velocity as a function of position yields a slope of @xmath41 , as seen , _ e.g. _ , in the equation : @xmath42 we find that the relaxation time is approximately isotropic ( esm figure  3 ( a ) ) and use @xmath43 as an estimate of @xmath7 . it is close to 10 s and only slightly varies with the pupa age group ( figure  [ tau_vs_logxi ] ( c ) , vertical axis ) . [ fig : sigma_anisotropy ] [ tau_vs_logxi ] mechanical state and material properties . color code according to pupa development ages : green , young ; red , middle - aged ; blue , old . the edges of the rectangular regions represent the mean values @xmath44 standard deviations for each group . the experiments in the scutum are in yellow . ( a ) strain anisotropy : @xmath33 _ versus _ @xmath32 . note the difference in horizontal and vertical scales ; the solid line is the first bisectrix @xmath45 , indicating the reference for isotropy . ( b ) same for the severed stress to viscosity ratio @xmath46 . ( c ) relaxation time @xmath7 and dimensionless friction to viscosity ratio @xmath47 ; values are the averages of the measurements along the @xmath4 and @xmath0 axes ( esm figure  3 ) . note the semi - log scale . the blue rectangle takes into account two very small values of @xmath48 , of order @xmath49 and @xmath50 ( below the plotted range ) . ]          to describe the bulk of the retracting domain , we track the displacements of around a hundred features in a band of tissue ( figure  [ fig : klt ] and methods ) . spatially averaging their positions by binning them into 8 equal - size bins improves the signal to noise ratio . we obtain the initial velocity profile _ versus _ position ( figure  [ fig : fit : xi ] ) . for old pupae , it is spatially linear ; for middle - aged pupae it is spatially nonlinear ; in the young pupae , immediately after severing only the boundaries move significantly .    to account for these observations , we formulate a spatio - temporal , viscoelastic , kelvin - voigt model in one dimension of space , where the coordinate @xmath51 denotes either @xmath4 or @xmath0 , and @xmath25 denotes time after severing ( see esm ) . thanks to the definition ( [ eq : def : hencky ] ) , the strain @xmath52 is related without approximation to the velocity @xmath53 through : @xmath54 we model the effect of friction of the epithelium against the hemolymph and the cuticle as an external fluid friction with coefficient @xmath55 . we find that when internal viscosity dominates , the strain remains uniform and all parts of the tissue relax exponentially with the same viscoelastic relaxation time @xmath56 where @xmath57 is the young modulus ( esm figure  4 ( a ) ) . when external friction dominates , strain diffuses from the boundaries with a diffusion coefficient @xmath58 and remains thus inhomogeneous at times @xmath25 smaller than @xmath59 ( esm figure  4 ( b ) ) . in the general case , the dynamical equation for the strain field @xmath60 reads : @xmath61 with the conditions @xmath62 initially , @xmath63 at boundaries , and @xmath64 at the end . the relative importance of external friction _ versus _ internal viscosity is quantified by the dimensionless parameter : @xmath65 solving equation  [ eq : general : d ] at short time and integrating over space equation  [ eq : velocity_strain ] at @xmath66 yields the initial velocity profile : @xmath67 whose curvature is controlled by @xmath48 . we thus estimate @xmath48 by fitting a hyperbolic sine to experimental initial velocity profiles ( figure  [ fig : fit : xi ] ) . as for the relaxation time , we check that the values @xmath68 and @xmath69 determined in the directions @xmath4 and @xmath0 are roughly consistent ( esm figure  3 ( b ) ) : we thus use @xmath70 as an estimate of @xmath48 . we find that @xmath48 varies more than @xmath7 ( figure  [ tau_vs_logxi ] ( c ) , horizontal axis ) . the middle aged group is characterized by intermediate values of @xmath48 ( figure  [ fig : fit : xi ] , red ) , of order of 10 . according to equation  [ eq : def : xibar ] , this indicates a possible effect of the external friction on the relaxation for scales larger than @xmath71 m . for both other age groups , the values of @xmath48 should be considered as bounds on the order of magnitude . for old pupae we find @xmath72 smaller than a few times unity ( figure  [ fig : fit : xi ] , blue ) . for young pupae , we find @xmath72 larger than a few tens or hundred ( figure  [ fig : fit : xi ] , green ) . this model depends on three parameters measured independently of each other : the initial strain @xmath73 , the viscoelastic time @xmath7 , and the friction to viscosity ratio @xmath48 . to validate the model , we use it to numerically compute the whole space - time map of strains and displacements we simulate each experiment using the values of @xmath7 and @xmath48 measured as described above . figure  [ boxtraj ] shows that we find a good agreement with experimental data . moreover , we plot the square deviation between experimental and numerical data as a function of the parameter values @xmath74 ( see methods ) . we check that the error landscapes we obtain are broadly consistent with our experimental estimates of @xmath48 and @xmath7 ( esm figure  5 ) . model - dependent measurement of the friction to viscosity ratio @xmath48 . initial velocity profiles immediately after severing are plotted _ versus _ initial position prior to severing , for three typical pupae : young ( green ) , middle - aged ( red ) and old ( blue ) . open squares : features , from figure  [ fig : klt ] ( b ) . closed triangles : same , spatially averaged in 8 bins . lines : fit by a sinh function ( equation  [ eq : vinit ] ) , yielding for @xmath48 a value either above , in , or below the measurable range : @xmath75 ( green ) , @xmath765 ( red ) , @xmath77 ( blue ) , respectively . ] validation of the model . blue triangles : positions of features from figure  [ fig : klt ] ( b ) , spatially averaged in bins _ versus _ time after severing ; bars : standard deviation . red lines : numerical resolution of equation  3 , using the measured values @xmath78 , @xmath79 . ] we report relaxation times @xmath7 in the range 5 to 20 s , and initial velocities in the range 0.1 to 0.6 @xmath12m.s@xmath80 . this agrees with values provided by single cell junction ablation : 10 - 100 s for the relaxation times , 0.1 to 1 @xmath12m.s@xmath80 for the initial velocities in other _ tissues @xcite . the model allows to interpret @xmath7 as a viscoelastic time . its order of magnitude is typical of the viscoelastic times corresponding to cell _ internal _ degrees of freedom @xcite , also measured in cell aggregates @xcite . conversely _ external _ degrees of freedom ( such as cell - cell rearrangements or cell divisions ) would have much larger viscoelastic times : in cell aggregates they are typically of order of several hours @xcite . the time scales of morphogenesis in this tissue are also of order of hours @xcite , as observed in other tissues as well , such as the _ drosophila _ pupal wing @xcite . indeed continuous models of morphogenetic processes presuppose the existence of such a separation of time scales .    the present experimental data are well described assuming that a single viscoelastic time is relevant . other experiments display a wide distribution of viscoelastic times @xcite : the model can easily be modified so as to take this feature into account . note that in the presence of friction the model already predicts that different parts of the severed domain relax at different rates ( esm figure  4 ( b ) ) . measurements of @xmath18 , @xmath9 , @xmath7 and @xmath48 are performed as averages over the severed disk . the requirements on the averaging scale are that it should be ( i ) large enough to obtain a sufficient signal to noise ratio yielding significant results , and ( ii ) small enough to be able to focus simultaneously on the whole severed ring despite the epithelium surface curvature which is observed in old pupae . the @xmath81 symmetry with respect to the midline observed in figures  [ fig : fit_ellipse ] ( b - d ) , [ fig : fit : xi ] , and [ boxtraj ] acts as a control of the signal to noise ratio . the experimentally available range of scales corresponds to 60 - 100 cell areas , in other words the semi - axis length @xmath19 corresponds to 4 - 5 cell diameters . for old pupae , scales smaller than @xmath19 are also probed , corresponding to the scale of bins ( at least @xmath82 or @xmath83 , see figs . [ fig : fit : xi ] and [ boxtraj ] ) . we experimentally determine the order of magnitude of the dimensionless friction to viscosity ratio , @xmath48 ( figure  [ tau_vs_logxi ] ( c ) , horizontal axis ) . obtaining more precise values of @xmath48 would depend on model ingredients , such as depth of ablation along the apico - basal axis , detailed cell geometry , tissue compressibility , or type of friction .    for simplicity we have modeled by a fluid friction @xcite the dissipation on either surfaces of the epithelium , namely the hemolymph at basal side and the cuticle at apical side . we do not observe any threshold effect due to solid friction , and the agreement between model and data justifies _ a posteriori _ our assumption .    for young pupae , the initial velocity field is not significantly higher than the noise level ( figure  [ fig : fit : xi ] , green ) . in this age group , it can not be used to validate the model , and the estimate of @xmath48 strongly depends on the motion of domain parts located close to the boundary ( equation  [ eq : vinit ] ) .      in the dorsal thorax , the severed domain fulfills the theoretical requirements for a continuous description : it contains a number of cells much larger than one ( large enough to allow for averaging ) , while being much smaller than the total number of cells in the whole tissue .    one could ask whether averages at scale @xmath19 are meaningful , due to variation of cell properties from cell to cell . we speculate that if heterogeneities within the domain dominated the average , the domain boundary could in principle adopt any arbitrary shape compatible with the @xmath81 symmetry . however , in practice , we observe that the domain boundary at the beginning , as well as at the end , of the relaxation also displays a @xmath84 symmetry ( figure  [ fig : fit_ellipse ] ( b - d ) and movies 1 - 3 ) and can be fitted by an ellipse : this is consistent with the measurements of @xmath9 and @xmath18 as tensors averaged over scale @xmath19 . limits of our model lie in its simplifications : one dimensional treatment , linear elasticity and viscosity , isotropy and homogeneity of material parameters including friction , identical properties of all cells contained in the severed region . however , the agreement between model and experiment is already good , and validates these approximations _ a posteriori_. we emphasize that one major goal of this model is to provide an order of magnitude of the friction to viscosity ratio @xmath48 . we do not claim that this model is unique , nor that it is general . in the future we hope to further test the model ingredients and relevance on this epithelium and other quasi - planar ones . since its ingredients are general , we expect our model to be of broad relevance when studying the relaxation of quasi - planar epithelia over similar time and length scales . however , the simplifications involved deserve further comments . in fact , completely different hypotheses could also be invoked to explain the spatial non - linearity in the velocity profile observed for young and middle - aged pupae . for instance , the viscous component of biological tissues could be a non - linear function of time , strain , and strain rate @xcite . a strain - dependent viscosity coefficient may lead to a non - linear velocity profile without the need to invoke any external friction . alternatively , a spatially heterogeneous friction might also explain the observed velocity profiles . since a linear velocity profile is observed in older pupae , the amplitude of these additional , more complex ingredients would have to decrease strongly during development . we believe that our model has the merit of simplicity . it successfully accounts for the description of the velocity profiles at any time . each of its hypotheses could be tested experimentally and , if experimental results require it , hypotheses could be relaxed to lead to more detailed explorations . we expect that additional ingredients could strongly affect the value of @xmath48 , but that other measurements presented here ( namely @xmath18 , @xmath9 and @xmath7 ) should be robust . the model thus provides a flexible and robust route to a mechanical description of quasi - planar epithelia . building upon classical single - junction laser ablations , our experiments enable to measure the time evolution of the mechanical state ( strain , and stress to viscosity ratio ) and material properties ( relaxation time , and friction to viscosity ratio ) during pupa metamorphosis , as well as their anisotropy . this provides relevant ingredients for the modeling of a tissue as a continuous , linear , viscoelastic material . moreover , our approach provides a method to explore the biological causality of the observed properties , by analysing mutant phenotypes or by micro - injecting drugs . we conjecture that @xmath7 would be affected by low doses of drugs known to modify cytoskeletal rheology such as cytochalasin d @xcite , ck-666 @xcite , y-27632 @xcite which respectively inhibit actin polymerization , arp2/3 complex activity or rho - kinase activity . we anticipate that modulating cell - cell adhesion by knock - down of e - cadherin or catenins @xcite could affect the viscosity . although we can not modulate the cell - cuticle adhesion , we could tune cell - matrix adhesion by knock - down of integrin or its linker to the cytoskeleton @xcite to affect the friction and thereby modify @xmath48 .    in principle , by establishing a complete map of @xmath8 at different positions , one could determine the length scale at which @xmath8 varies within the tissue , and check _ a posteriori _ whether it is larger than the scale @xmath19 used to measure @xmath8 . however , the present experiments measure the stress up to a dissipative prefactor , the tissue viscosity . similarly , single cell - junction ablation experiments measure force up to a prefactor ( a friction coefficient ) , rarely discussed or measured in the literature . to separate the stress and viscosity variations , alternative methods to directly measure forces and viscosities in live epithelia are thus called for . these include the measurement of cell junction tensions from movies observations @xcite , where the unknown prefactor , namely the average tension over the whole image , is by definition the same for all junctions . external mechanical manipulation @xcite yields direct determination of out - of - plane elastic and/or viscous moduli , but the in - plane moduli are only indirectly determined . another possible strategy could rely on microrheology @xcite : either passive , by tracking the brownian diffusion of beads within the cells ; or active , by moving ( magnetically or optically ) a bead back and forth within a cell at a specified frequency . such measurements could yield a direct access to the value of @xmath10 , but at the intracellular level . measuring @xmath10 at the scale of several cells either directly , or indirectly by measuring @xmath57 and @xmath7 and using the relation @xmath85 , could probably be possible , but we are not aware yet of a published method . since @xmath18 is a dimensionless geometrical quantity , its value is insensitive to cell size : the observed @xmath117-fold increase in @xmath18 can not be related to changes of cell size over time . it should involve other causes , that remain to be determined . part of the measured increase in @xmath9 values may be due to cell divisions which decrease the average cell area by a factor of @xmath86 between young and old pupae . the observed increase is a factor of @xmath87 , and thus should involve other causes as well . one candidate is a decrease of the tissue viscosity @xmath10 . however the near constancy of the viscoelastic time would then imply a corresponding increase in tissue elasticity , a rather unlikely coincidence . we therefore expect that stronger tensile forces are at work in the tissue in older pupae . from 18 hours apf onwards , the lateral part of the scutellum undergoes apical cell contractions that shape the lateral domain of the tissue @xcite . this observation led us to divide the ablation experiments into the three age groups . it may also explain the changes in mechanical properties measured in the tissue for middle - aged and old pupae . to precisely analyse the causes of the increases in @xmath18 and @xmath9 would therefore require to identify genes that specifically affect the lateral contraction of the tissue . our experimental method should be applicable to other tissues . performing experiments in an epithelium which is flat enough to significantly increase @xmath19 , and with a sufficient signal to noise ratio on the estimate of @xmath48 , would enable to check the scaling of @xmath48 _ versus _ @xmath19 given by equation  [ eq : def : xibar ] . _ drosophila melanogaster _ larvae were collected at the beginning of metamorphosis ( pupa formation ) @xcite . the pupae were kept at 25@xmath88c , then dissected and mounted as described in @xcite ( esm figure  1 ) . using the adherens junction protein e - cadherin fused to green fluorescent protein @xcite we imaged the apical cell junctions by fluorescent microscopy @xcite . the time - lapse laser - scanning microscope lsm710 nlo ( carl zeiss microimaging ) collected 512@xmath89512 pixel images in mono - photon mode at 488 nm excitation ( pixel size in the range 0.24 to 0.29 @xmath12 m ) through a plan - apochromat 63@xmath891.40 oil objective . before and after severing , the time interval between two consecutive frames was either 393 ms or 970 ms according to the scanning mode ( bi- or mono - directional ) , without any apparent effect on the results presented here . we focused on the adherens junctions to be severed ; due to the epithelial surface curvature especially in old pupae , adherens junctions of innermost and outermost regions were slightly out of focus ( figure  [ fig : fit_ellipse ] and movies 1 - 3 ) .    we defined a region of interest as an annular region between two concentric circles . we severed it at the 10@xmath90 frame ( figure  [ fig : fit_ellipse ] ) and defined @xmath19 as the radius of its inner circle . laser severing was performed using ti : sapphire laser ( mai tai deepsee , spectra physics ) in two - photon mode at 890  nm , @xmath91  100  fs pulses , 80  mhz repetition rate , @xmath11 0.2  w at the back focal plane , used at slightly less than full power to avoid cavitation ( esm figure  2 ) @xcite . severing itself had a duration ranging from 217 to 1300 ms ( according to the size of the region of interest and the scanning mode ) during which no image was acquired . we could check that the severing had been effective to remove their adherens junctions , rather than simply bleach them ( esm figure  2 ) : after retraction the cells at the wound margin moved apart more than a cell diameter ; moreover , there was no fluorescence recovery ( figure  [ fig : fit_ellipse ] ( d ) ) .      at the tissue scale , for each frame after the severing the inner severed tissue boundary was fitted with an ellipse ( figure  [ fig : fit_ellipse ] ( c - d ) ) using the ovuscule @xcite imagej plugin @xcite . we have modified it so that each fitting procedure started from the ellipse fitted on the preceding   are linearly fitted over a 5 frames sliding window in order to improve the signal to noise of position and velocity determinations . independently , at the cell scale , images were denoised with safir software @xcite . we took the logarithm of the grey levels , to obtain comparable intensity gradients in differently contrasted parts of the image ( figure  [ fig : klt ] ( a ) ) . this allowed us to use kanade - lucas - tomasi ( klt ) @xcite tracking algorithm @xcite by selecting @xmath92  features " of interest , corresponding to most cell vertices ( figure  [ fig : klt ] ( a ) ) . features were tracked from frame to frame : they followed the cell vertex movement towards the center , and we have observed no cell neighbour swapping . we checked that the features center of mass had only a small displacement , which we subtracted from individual feature displacements without loss of generality . to implement a quasi-1d analysis , we collected the features in a band along the @xmath0 axis and similarly along the @xmath4 axis ( figure  [ fig : klt ] ( b - c ) ) , and both axes were analysed separately . the band width ( here @xmath93 of the initial circle diameter ) was chosen to be large enough to perform large statistics on features , and small enough to treat the features displacement as one - dimensional . we numerically solved adimensionalized equations ( see esm ) using the matlab solver ` pdepde ` . for a given set of parameter values @xmath94 , we calculated @xmath73 from the strain of the feature of largest initial amplitude , and rescaled the solution of esm equation  21 so as to satisfy esm equation  10 . using the experimental initial positions of features @xmath95 and the calculated strain field @xmath52 , we simulated the trajectories @xmath96 , @xmath97 , where @xmath98 was the number of features . the mismatch @xmath99 with our computation was evaluated by comparison with the experimental trajectories , @xmath100 : @xmath101 where @xmath102 denoted the euclidean norm , and @xmath103 the number of movie images . we gratefully thank f. molino for modifying the ovuscule plugin , p. thvenaz and s. birchfield for advices regarding their softwares , h. oda for reagents , j .- allain for suggesting to use klt algorithm , lab members for discussions , and the pict - ibisa@bdd ( umr3215/u934 ) imaging facility of the institut curie . this work was supported by grants to y.b . from arc ( 4830 ) , anr ( blan07 - 3 - 207540 ) , erc starting grant ( cepodro 209718 ) , cnrs , inserm , and institut curie ; and by postdoc grants to i.b . by the frm ( spf20080512397 ) , to f.b . by the nwo ( 825.08.033 ) . 90    lecuit t , lenne pf 2007 cell surface mechanics and the control of cell shape , tissue patterns and morphogenesis . nat rev mol cell biol _ 8:633644 . chen cs , tan j , tien j 2004 mechanotransduction at cell - matrix and cell - cell contacts . _ annu rev biomed eng _ 6:275302 . ma x , lynch he , scully pc , hutson ms 2009 probing embryonic tissue mechanics with laser hole drilling . _ phys biol _ 6:036004 . lecuit t , le  goff l 2007 orchestrating size and shape during morphogenesis . _ nature _ 450:189192 . mammoto t , ingber de 2010 mechanical control of tissue and organ development . _ development _ 137:14071420 . fung yc 1993 _ biomechanics : mechanical properties of living tissues _ ( springer ) . hufnagel l , teleman aa , rouault h , cohen sm , shraiman bi 2007 on the mechanism of wing size determination in fly development . _ proc natl acad sci usa _ 104:3835 - 3840 . bittig t , wartlick o , kicheva a , gonzalez - gaitan m , jlicher f 2008 dynamics of anisotropic tissue growth . _ new j phys _ 10:063001 . blanchard g , _ et  al . _ 2009 tissue tectonics : morphogenetic strain rates , cell shape change and intercalation . _ nat meth _ 6:458 - 464 . butler l , _ _ 2009 cell shape changes indicate a role for extrinsic tensile forces in drosophila germ - band extension . _ nat cell biol _ 11:859 - 864 .    ranft j , _ _ 2010 fluidization of tissues by cell division and apoptosis . _ proc natl acad sci usa _ 107:2086320868 . aigouy b , _ _ 2010 cell flow reorients the axis of planar polarity in the wing epithelium of drosophila . _ cell _ 142:773786 . batchelor gk 2000 _ an introduction to fluid dynamics _ ( cambridge university press )    ophir j , cespedes i , garra b , ponnekanti h , huang y , maklad n , 1996 elastography : ultrasonic imaging of tissue strain and elastic modulus in vivo _ eur j ultrasound _ 3:49 - 70 . nienhaus u , aegerter - wilmsen t , aegerter cm 2009 determination of mechanical stress distribution in drosophila wing discs using photoelasticity . _ mech dev _ 126:942949 . desprat n , supatto w , pouille p , beaurepaire e , farge e 2008 tissue deformation modulates twist expression to determine anterior midgut differentiation in drosophila embryos . _ dev cell _ 15:470477 . fleury v , _ et  al . _ 2010 introducing the scanning air puff tonometer for biological studies . _ phys rev e _ 81:021920 . peaucelle a , braybrook s , le guillou l , bron e , kuhlemeier c , hfte h , 2011 pectin - induced changes in cell wall mechanics underlie organ initiation in _ arabidopsis_. _ curr biol _ 21:1720 - 1726 . rauzi m , lenne pf 2011 cortical forces in cell shape changes and tissue morphogenesis . _ curr top dev biol _ 95:93144 . hutson ms , _ _ 2003 forces for morphogenesis investigated with laser microsurgery and quantitative modeling . _ science _ 300:145149 . hutson ms , ma x 2007 plasma and cavitation dynamics during pulsed laser microsurgery in vivo . _ phys rev lett _ 99:158104 . hutson ms , _ et  al . _ 2009 combining laser microsurgery and finite element modeling to assess cell - level epithelial mechanics . _ biophys j _ 97:30753085 . kiehart dp , galbraith cg , edwards ka , rickoll wl , montague ra 2000 multiple forces contribute to cell sheet morphogenesis for dorsal closure in drosophila . _ j cell biol _ 149:471490 . peralta xg , _ et  al . _ 2007 upregulation of forces and morphogenic asymmetries in dorsal closure during drosophila development . _ biophys j _ 92:25832596 . farhadifar r , rper j , aigouy b , eaton s , jlicher f 2007 the influence of cell mechanics , cell - cell interactions , and proliferation on epithelial packing . _ curr biol _ 17:20952104 . rauzi m , verant p , lecuit t , lenne pf 2008 nature and anisotropy of cortical forces orienting drosophila tissue morphogenesis . _ nat cell biol _ 10:14011410 . landsberg kp , _ _ 2009 increased cell bond tension governs cell sorting at the drosophila anteroposterior compartment boundary . _ curr biol _ 19:19501955 . r , simoes sdm , rper j , eaton s , zallen ja 2009 myosin ii dynamics are regulated by tension in intercalating cells . _ dev cell _ 17:736743 . bosveld f , bonnet i , guirao b _ et al . _ 2012 mechanical control of morphogenesis by fat / dachsous / four - jointed planar cell polarity pathway . _ science _ , in press . mayer m , depken m , bois js , jlicher f , grill sw 2010 anisotropies in cortical tension reveal the physical basis of polarizing cortical flows . _ nature _ 467:617621 . martin ac , gelbart m , fernandez - gonzalez r , kaschube m , wieschaus ef 2010 integration of contractile forces during tissue invagination . _ j cell biol _ 188:735749 . janiaud e , graner f , 2005 foam in a two - dimensional couette shear : a local measurement of bubble deformation . _ j fluid mech _ 532:243267 . graner f , dollet b , raufaste c , marmottant p 2008 discrete rearranging disordered patterns , part i : robust statistical tools in two or three dimensions . _ eur phys j e _ 25:349369 . bainbridge sp , bownes m 1981 staging the metamorphosis of _ drosophila melanogaster_. _ j embryol exp morph _ 66:57 - 80    langevin j , _ et  al . _ 2005 drosophila exocyst components sec5 , sec6 , and sec15 regulate de - cadherin trafficking from recycling endosomes to the plasma membrane . _ dev cell _ 9:365376 . thvenaz p , delgado - gonzalo r , unser m 2011 the ovuscule . _ ieee trans pattern anal mach intell _ 33:382393 . hoger a , 1987 the stress conjugate to logarithmic strain . _ int j solids struct _ 23:1645 - 1656 . farahani k , naghdabadi r , 2000 conjugate stresses of the seth - hill strain tensors . _ int j solids struct _ 37:5247 - 5255 . tanner r , tanner e 2003 heinrich hencky : a rheological pioneer . _ rheologica acta _ 42:93101 . wottawah f , _ _ 2005 optical rheology of biological cells . _ phys rev lett _ 94:098103 . marmottant p , _ et  al . _ 2009 the role of fluctuations and stress on the effective viscosity of cell aggregates . _ proc natl acad sci usa _ 106:1727117275 . guevorkian k , gonzalez - rodriguez d , carlier c , dufour s , brochard - wyart f 2011 mechanosensitive shivering of model tissues under controlled aspiration . _ proc natl acad sci usa _ 108:1338713392 . guevorkian k , colbert m , durth m , dufour s , brochard - wyart f 2010 aspiration of biological viscoelastic drops . _ phys rev lett _ 104:218101 . flanagan md , lin s 1980 cytochalasins block actin filament elongation by binding to high affinity sites associated with f - actin . _ j biol chem _ 255:835 - 838 . nolen bj , tomasevic n , russell a , pierce dw , jia z , mccormick cd , hartman j , sakowicz r , pollard td 2009 characterization of two classes of small molecule inhibitors of arp2/3 complex . _ nature _ 460:1031 - 1034 . uehata m , ishizaki t , satoh h , ono t , kawahara t , morishita t , tamakawa h , yamagami k , inui j , maekawa m , narumiya s. 1997 calcium sensitization of smooth muscle mediated by a rho - associated protein kinase in hypertension . _ nature _ 389:990 - 994 . papusheva e , heisenberg cp 2010 spatial organization of adhesion : force - dependent regulation and function in tissue morphogenesis . _ embo j. _ 29:2753 - 2768 . brodland gw , _ et  al . _ 2010 video force microscopy reveals the mechanics of ventral furrow invagination in drosophila . _ proc natl acad sci usa _ 107:2211122116 . mackintosh f , schmidt c 1999 microrheology . _ curr opinion coll interface science _ 4:300 - 307 . panorchan p , lee jsh , daniels br , kole tp , tseng y , wirtz d , 2007 probing cellular mechanical responses to stimuli using ballistic intracellular nanorheology . _ methods cell biol _ 83:115 - 140    crocker jc , hoffman bd 2007 multiple - particle tracking and two - point microrheology in cells . _ methods cell biol _ 83:141178 . roberts db 1998 _ drosophila : a practical approach _ ( oxford university press , usa ) , 2 edition . sgalen m , _ et  al . the fz - dsh planar cell polarity pathway induces oriented cell division via mud / numa in drosophila and zebrafish . _ dev cell _ 19:740752 . oda h , tsukita s 2001 real - time imaging of cell - cell adherens junctions reveals that drosophila mesoderm invagination begins with two phases of apical constriction of cells . _ j cell sci _ 114:493501 . / thevenaz / ovuscule/.    kervrann c , boulanger j 2006 optimal spatial adaptation for patch - based image denoising . _ ieee trans image process _ 15:28662878 . boulanger j , _ et  al . _ 2010 patch - based nonlocal functional for denoising fluorescence microscopy image sequences . _ ieee trans med imaging _ 29:442454 . lucas b , kanade t 1981 an iterative image registration technique with an application to stereo vision . _ international joint conference on artificial intelligence _ pp 674679 . tomasi c , kanade t 1991 detection and tracking of point features . _ carnegie mellon university technical report _ cmu - cs-91 - 132 . feature tracker . http://www.ces.clemson.edu/~stb / klt/.    electronic supplementary material ( esm ) the external dissipation is not accessible to single cell junction ablation experiments and its estimate requires modelling @xcite . our model is analogous to models of viscoelastic relaxation proposed for local severing experiments of acto - myosin bundles in living cells @xcite . we use an eulerian description . since the severed domain area visibly decreases with time , we do not enforce any in - plane incompressibility condition . note also that in - plane incompressibility would fix the velocity profile , which is incompatible with the experimental observation of different profiles ( figure  5 ) . we rather assume for simplicity that the out - of - plane thickness is a passive variable : hence the 3d poisson ratio is @xmath104 , but the in - plane effective 2d poisson ratio vanishes . we test this hypothesis by checking _ a posteriori _ that the model predictions are compatible with experimental observations . within this simplification , the antero - posterior @xmath4 and medio - lateral @xmath0 axes ( figure  1 ) then decouple . finally , we define @xmath8 as the part of stress which has been removed by ablation . we thus consider only one spatial dimension , with the variable @xmath51 denoting position along either @xmath4 or @xmath0 . we model the retracting slab of tissue as a linear viscoelastic kelvin - voigt material ( that is , with a steady response at long time ) with a young modulus @xmath57 and a viscosity @xmath10 . we neglect the line tension around the tissue domain : first , because we observe that the ellipse does not remain round ( this suggests the line tension is weak , without being a proof , since we could also imagine that it has the same anisotropy as the strain ) ; second , because even if the line tension exists , it is small enough that the ellipse contour remains visibly rough , see movies 1 - 3 . two terms contribute explicitly to @xmath8 : an elastic one ( which implicitly includes the external pressure ) and a viscous one : @xmath105 taking into account a fluid external friction , with a friction coefficient @xmath55 of dimension pa.s.m@xmath106 , the force balance equation reads : @xmath107 combining equations  3 , [ eq : stress ] and [ eq : newton ] yields equation  4 . the initial length of the retracting slab is @xmath108 : @xmath109 $ ] . assuming that the initial state respects mechanical equilibrium ( over the relevant time scale of 10 to 100  s ) , the initial strain is uniform , and we have for all @xmath51 : @xmath110 at the boundary of the severed region , we have for all @xmath25 : @xmath111 since the final state is at mechanical equilibrium , we have for all @xmath51 : @xmath112 injecting equation  3 into equation  [ eq : stress ] , using conditions given by equations  [ eq : ic]-[eq : fin : sigma ] yields by integration over time : @xmath113 where the relaxation time is viscoelastic : @xmath114 . equations  3 , [ eq : ic ] , [ eq : bc : e ] yield an approximate analytical prediction for the initial velocity profile at time @xmath115 , as follows . an expansion of @xmath18 at first order in @xmath116 : @xmath117 defines the initial strain rate as : @xmath118 substituting equations  [ eq : e1:expand]-[eq : e1:def ] into equation  4 , we obtain an ordinary differential equation for @xmath119 : @xmath120 with a diffusion coefficient @xmath121 . with the boundary condition @xmath122 , the solution of equation  [ eq : e1:eqdiff ] reads : @xmath123 where @xmath124 , or @xmath125 since @xmath126 . spatial integration of equation  [ eq : e1:sol ] using equations  3 and [ eq : e1:def ] yields equation  6 .      esm figure  [ sfig : twolimits ] shows two limit regimes of equation  4 . @xmath127 when internal viscosity dominates :    external friction may be neglected , and equation  [ eq : newton ] reduces to @xmath128 so that @xmath8 is uniform and its value is fixed at the boundaries : @xmath129 for all @xmath130 . integration yields : @xmath131 the strain is uniform at all times ( esm figure  [ sfig : twolimits ] ( a ) ) . since the center of the slab does not move , @xmath132 for all @xmath25 , the velocity profile is linear : @xmath133 as expected , in the limit of small friction @xmath134 , equation  6 reduces to equation  [ eq : initv_modelfree ] when @xmath135 . in this regime , all parts of the tissue relax exponentially with the same time scale @xmath7 . @xmath127 when external friction dominates :    the internal viscosity of the tissue may be neglected . combining equations  3 , [ eq : stress ] and [ eq : newton ] yields a diffusion equation for the strain field : @xmath136 with conditions given by equations  [ eq : ic]-[eq : fin : sigma ] . strain diffuses from the boundaries on a time scale @xmath137 , and the strain field remains inhomogeneous until @xmath138 ( esm figure  [ sfig : twolimits ] ( b ) ) . one may note that the strain diffusion equation ( equation  [ eq : diff : e ] ) is identical to that obtained in the context of gel dynamics @xcite . taking into account interstitial flow ( _ i.e. _ the transport of water across the epithelium ) would only modify the diffusion constant : @xmath139 , where @xmath140 is the permeability coefficient .      in the general case , to solve equation  4 numerically we normalize variables by defining : @xmath141 , @xmath142 and @xmath143 . equation  4 then becomes : @xmath144 with the initial condition : @xmath145 the boundary conditions : @xmath146 and the final conditions : @xmath147    since @xmath148 , the limit regimes are recovered when @xmath149 ( internal viscosity dominates , @xmath150 ) , and @xmath151 ( external friction dominates , @xmath152 ) . 90    mayer m , depken m , bois js , jlicher f , grill sw 2010 anisotropies in cortical tension reveal the physical basis of polarizing cortical flows . _ nature _ 467:617621 . colombelli j , et  al . 2009 mechanosensing in actin stress fibers revealed by a close correlation between force and protein localization . _ j cell sci _ 122:16651679 . charras gt , mitchison tj , mahadevan l 2009 animal cell hydraulics . _ j cell sci _ 122:32333241 . doi m 2009 gel dynamics . _ j phys soc jpn _ 78:052001 . ( esm ) set - up for live imaging . ( a ) the pupal case ( light brown ) is partially removed , while the cuticle ( transparent ) is kept intact . blue cross : approximate position of severing , along the symmetry axis of the scutellum . ( b ) the pupa ( white ) is mounted between slide and coverslip with its back facing the objective ( dark grey ) . antero - posterior axis is horizontal ( anterior towards the left , posterior towards the right).,scaledwidth=70.0% ]     ( esm ) range of acceptable severing powers . images shown are just before the severing ( top ) and 24.3 s after the severing ( bottom ) as a function of the proportion ( indicated at top left ) of the full laser power , here at 890  nm . below 20% , no severing occurs . from 30 to 50% , severing is partial . from 60 to 80% , severing is complete and cell - cell junctions move significantly ; experiments presented in this work are performed in this range . above 90% , cavitation occurs . these ranges vary slightly from one experiment to the other , so that the quality of severing must be carefully assessed _ a posteriori_.,scaledwidth=70.0% ]     [ sfig : tau_anisotropy ] ( esm ) check of @xmath7 and @xmath48 isotropy . ( a ) relaxation time @xmath7 : for each experiment , represented by one point , we plot the value obtained along the @xmath0 direction ( estimated as in the inset of figure  2 ) _ versus _ the value obtained along the @xmath4 direction . age is color - coded as in figure  3 : green , young ; red , middle - aged ; blue , old . ( b ) same plot for the friction / viscosity ratio @xmath48 , fitted by equation  6 as in figure  5 . note the logarithmic scale . two values of @xmath48 of order @xmath49 and @xmath50 are below the plotted range . , scaledwidth=50.0% ]    , @xmath153 s , @xmath154 m . ( a ) uniform strain relaxation , when internal viscosity dominates ( low @xmath48 , initial velocity profile similar to the blue curve in figure  5 ) : @xmath155 , @xmath156 . ( b ) strain diffusion , when external friction dominates ( high @xmath48 , initial velocity profile similar to the green curve in figure  5 ) : @xmath157 , @xmath158 . strain value in each grid square is color - coded from 0 ( dark blue ) to 1 ( dark red ) . [ sfig : twolimits],scaledwidth=50.0% ]     ( esm ) error maps in the @xmath159 plane , semi - log axes . experiments are those of figure  5 and movies 1 - 3 ; error is estimated according to equation  7 . the minimal region of the error landscape ( its deepest valley ) is too wide to allow a direct measurement of the parameters . however , depending on the orientation of the minimal region of the error landscape , the computed error maps fall broadly into three categories , which agree with the three age groups . there is also a good agreement with hatched regions , which are the results of the independent measurements of the parameters @xmath7 and @xmath48 ( figure  3 ( c ) and esm figure  [ sfig : xi_anisotropy ] ) . ( a ) small @xmath48 , old pupa ( movie 3 ) . minimal valley parallel to the @xmath48 axis : this suggests that only an upper bound , of the order of @xmath160 ( white bar ) , can be given for the value of @xmath48 . friction is negligible . ( b ) finite @xmath48 , middle - aged pupa ( movie 2 ) . minimal valley intermediate , parallel to the @xmath7 axis : friction can not be neglected , and @xmath48 admits a finite estimate . ( c ) large @xmath48 , young pupa ( movie 1 ) . minimal valley parallel to the @xmath7 axis , at very large @xmath48 : this suggests that only a lower bound can be given for the values of @xmath7 and @xmath48 . ]
Abstract:

Answer: during development , epithelial tissues undergo extensive morphogenesis based on coordinated changes of cell shape and position over time . continuum mechanics describes tissue mechanical state and shape changes in terms of strain and stress . it accounts for individual cell properties using only a few spatially averaged material parameters . to determine the mechanical state and parameters in the _ drosophila _ pupa dorsal thorax epithelium , we sever _ in vivo _ the adherens junctions around a disk - shaped domain comprising typically hundred cells . this enables a direct measurement of the strain along different orientations at once . the amplitude and anisotropy of the strain increase during development . we also measure the stress to viscosity ratio and similarly find an increase in amplitude and anisotropy . the relaxation time is of order of ten seconds . we propose a space - time , continuous model of the relaxation . good agreement with experimental data validates the description of the epithelial domain as a continuous , linear , visco - elastic material . we discuss the relevant time and length scales . another material parameter , the ratio of external friction to internal viscosity , is estimated by fitting the initial velocity profile . together , our results contribute to quantify forces and displacements , and their time evolution during morphogenesis .


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: our understanding of stellar structure and stability has been built upon classical analytical studies of equilibrium configurations of incompressible self - gravitating systems by maclaurin , jacobi , dedekind , darwin , riemann , and especially @xcite . with the help of modern computers and numerical methods , such as the self - consistent - field ( scf ) technique @xcite and hydrodynamics techniques , these studies have been extended to compressible configurations by a variety of researchers @xcite . however , most three - dimensional ( 3d ) hydrodynamical studies have been confined to a limited parameter space , i.e. , starting from two - dimensional ( 2d ) axisymmetric equilibrium models ( compressible analogues of maclaurin spheroids ) because we have not been able to build 3d compressible equilibrium models with complicated flows . to the author s knowledge , techniques have only been developed to build compressible equilibrium models of nonaxisymmetric configurations for a few systems with simplified rotational profiles , e.g. , rigidly rotating systems @xcite , irrotational systems @xcite , and configurations that are stationary in the inertial frame @xcite . our aim is to find a method by which a much wider class of nonaxisymmetric configurations can be readily constructed . this will permit us to more clearly understand the viability of evolutionary sequences of triaxial configurations and will provide models whose stability properties can be readily analyzed . previous linear stability analyses @xcite have revealed a number of hydrodynamical instabilities that might arise in incompressible riemann ellipsoids , but no corresponding numerical work using state - of - the - art 3d hydrodynamic techniques has been carried out . in particular , riemann s - type ellipsoids have been found to be subject to a hydrodynamic strain instability associated with elliptical stream lines @xcite , which raises concerns about the stability of certain types of geophysical flows and leads to suspicions about the evolutionary path of stars that are driven by gravitational - radiation - reaction ( grr ) forces toward the dedekind - sequence . in a recent nonlinear study of the secular bar - mode instability induced by grr @xcite , a uniformly rotating , axisymmetric neutron star , which was secularly unstable and driven by artifically enhanced grr , evolved into a bar - like configuration with a very slow pattern speed . however , this special bar configuration became unstable to some kind of turbulence - like instability while evolving toward the dedekind - sequence , that is , toward a stationary ellipsoidal structure in the inertial frame maintained purely by the internal motion of the fluid . it is suspected that this turbulence - like instability may be the elliptical strain instability identified in the earlier linear stability analyses . but no definite conclusion can be drawn because the stability of only one dedekind - like configuration was studied . following @xcite , one might consider generating additional ellipsoidal structures hydrodynamically , but the generation of each additional ellipsoidal model would be expensive in the sense that one has to evolve an initially axisymmetric model for a very long time , even with artificially enhanced grr @xcite . this approach would make it impractical to undertake a full stability investigation through the entire parameter space of ellipsoidal models . here we present a new 3d scf technique that is capable of building the full range of incompressible riemann s - type ellipsoids with nontrivial internal flows ( i.e. , not just the uniformly rotating jacobi configurations , which are thought to be the end point of secularly unstable stars driven by viscosity , and stationary dedekind configurations ) as well as compressible triaxial models that share the same divergence - free velocity field as riemann s - type ellipsoids . our compressible models satisfy the steady - state euler s equation exactly , but only satisfy the steady - state continuity equation approximately , hence they are only in quasi - equilibrium . however , the violation of the steady - state continuity equation is very small for nearly incompressible models , hence these models can be used effectively as initial configurations for a hydrodynamical study of the stability of compressible analogues of riemann s - type ellipsoids . we note that , to date , hydrodynamical techniques have not been used to examine the stability of riemann ellipsoids because the sound speed is formally infinite for incompressible models , while the time step in hydrodynamic codes is usually inversely proportional to the sound speed . on the other hand , since the elliptical strain instability associated with elliptical stream lines would set in on a dynamical time scale , we believe that our new `` quasi - equilibrium '' models will prove to be good starting points for hydrodynamical studies of this instability . our numerical method inherits one central idea of the traditional scf method in that it makes use of the properties of the velocity field of riemann s - type ellipsoids to turn the vector form of the steady - state euler s equation into a scalar form , thus yielding a uniform bernoulli s function . this enables us to follow the procedures of the traditional scf method to use several boundary conditions to determine quantities elsewhere . to be concise , we will refer readers to @xcite and references therein for details of the traditional scf technique , and focus on the new aspects of our new method . figure [ triaxial ] reproduces key elements of a plot that ( * ? ? ? * see his chapter 7 , figure 15 ) used to discuss sequences of riemann s - type ellipsoids . in this ( @xmath0,@xmath1 ) parameter plot , @xmath2 , @xmath3 , and @xmath4 are the three principal semiaxes of an ellipsoid , and each riemann s - type ellipsoid with angular velocity @xmath5 and uniform vorticity @xmath6 in the rotating frame resides between the two self - adjoint sequences ( curves lp and lk ) with @xmath7 , where @xmath8 the velocity field of a riemann s - type ellipsoid as viewed from a frame rotating with angular velocity @xmath5 takes the following form : @xmath9 where @xmath10 is a constant that determines the magnitude of the internal motion of the fluid , and the origin of the x - y coordinate system is at the center of the ellipsoid . this velocity field @xmath11 is designed so that velocity vectors everywhere are always aligned with elliptical stream lines by demanding that they be tangent to the equipotential contours , which are concentric ellipses . now , let s turn to a general self - gravitating system with some compressibility . if we adopt a polytropic equation of state , @xmath12 ( @xmath13 is pressure , @xmath14 is the polytropic constant and @xmath15 is the polytropic index ) , euler s equation for steady - state flows in the rotating frame is , @xmath16          - 2 \vec{\omega } \times \vec{v } \ , ,   \label{euler}\ ] ] where h is enthalpy , and @xmath17 is the gravitational potential . with the velocity field specified by eq . ( [ vel ] ) , one can show that eq . ( [ euler ] ) is equivalent to the following equation , @xmath18 =       -\nabla [ h+\phi-\frac{1}{2}\omega^2 ( x^2+y^2 ) ] \nonumber\\       -\nabla [ \omega \lambda(\frac{b}{a}x^2 + \frac{a}{b}y^2 ) ] \,.\end{aligned}\ ] ] hence , within the configuration the following bernoulli s function must be uniform in space : @xmath19 where @xmath20 is a constant . it should be mentioned that in this bernoulli s function , @xmath5 and @xmath10 are interchangeable . this reflects the fact that , for any direct configuration ( @xmath21 ) in which fluid motion is dominated by figure rotation , the adjoint configuration ( @xmath22 ) in which fluid motion is dominated by internal motions can be obtained by simply interchanging @xmath5 and @xmath10 . we find it is useful to define an effective potential @xmath23 as : @xmath24 hence , the bernoulli s function can also be written as : @xmath25 in our later discussion of equi - potential ( iso - density ) surfaces , we will often use the phrase `` equi - potential '' referring to @xmath23 .    on the other hand , the steady - state continuity equation in the rotating frame requires , @xmath26 this is trivially satisfied for incompressible riemann s - type ellipsoids because @xmath27 is zero everywhere , and the velocity vectors in the rotating frame are always tangent to the concentric equi - potential contours and perpendicular to the unit vector that is normal to the ellipsoidal surface . in our compressible models , the divergence - free velocity field given by eq . ( [ vel ] ) will continue to be adopted , so if the equi - potential ( and iso - density ) contours continue to be concentric ellipses ( as assumed by * ? ? ? * ) , the velocity vectors will always be orthogonal to the density gradient everywhere ; this would satisfy the steady - state continuity equation exactly . however , as we demand that euler s equation be satisfied via the solution of eq . ( [ bernoulli ] ) , the iso - density contours for our compressible models are not guaranteed to be self - similar ellipses . therefore , in our current numerical method , the extent to which the steady - state continuity equation is satisfied will depend on how far the density shells deviate from self - similar ellipses . for the remainder of this paper , we will only focus on obtaining a solution of bernoulli s function given by eq . ( [ bernoulli ] ) . we will return to issues related to the steady - state continuity equation in section [ discussion ] .    for a given choice of axes ( @xmath2 , @xmath3 , @xmath4 ) and a mass distribution @xmath28 ( hence @xmath29 is determined by poisson s equation ) , there are three unknown constants in eq . ( [ bernoulli ] ) : @xmath20 , @xmath5 and @xmath10 . on the other hand , we have three boundary conditions on the three axes at the surface of the star , i.e. , points a , b and c as illustrated in figure [ boundary ] , where @xmath30 drops to zero . ( note that points a , b and c are located on @xmath31 , @xmath32 , and @xmath33 axes at @xmath34 , @xmath35 , and @xmath36 , respectively . ) these boundary conditions give the following equations : @xmath37 subtracting equation ( [ bc2 ] ) from ( [ bc1 ] ) , we obtain @xmath38 adding equation ( [ bc2 ] ) and ( [ bc1 ] ) together , we obtain @xmath39 \label{ch2 } \,.\ ] ] since the right - hand - side of this last expression is non - negative , it follows that @xmath40 . it can be shown that equations ( [ ch1 ] ) and ( [ ch2 ] ) are equivalent to equations ( 29 ) and ( 30 ) in chapter 7 of @xcite . the solutions of equations ( [ ch1 ] ) and ( [ ch2 ] ) for a direct configuration are : @xmath41 where @xmath42 for adjoint configurations , we only need to interchange @xmath43 and @xmath44 . the signs of @xmath5 and @xmath10 are determined by the sign of @xmath45 ( which is predetermined by the specified axis ratios ) , through eq . ( [ eqn : f ] ) and the relationship , @xmath46 with these constants and the bernoulli s function in hand , we can solve for @xmath30 throughout the interior of the configuration .    in the incompressible case ( @xmath47 ) , we know from the classical analytical results that the adopted velocity field generates equilibrium structures that are uniform density ellipsoids . hence , we can obtain riemann s - type ellipsoids directly by setting up a mass distribution that is confined by the ellipsoidal surface defined by three axes @xmath2 , @xmath3 and @xmath4 ; then we solve poisson s equation to get the potential everywhere ; finally , @xmath5 and @xmath10 can be determined by equations ( [ omega ] ) and ( [ lambda ] ) . compressible polytropic models ( @xmath48 ) having the same flow pattern as that of riemann ellipsoids ( given by eq . [ vel ] ) , can also be constructed by an iterative procedure that is very similar to hachisu s scf method @xcite . the main steps of this method include : ( i ) set up a trial ellipsoidal mass distribution defined by the choice of @xmath2 , @xmath3 , and @xmath4 , then solve poisson s equation to obtain the gravitational potential everywhere ( the details of our poisson solver are discussed in * ? ? ? * ) ; ( ii ) calculate @xmath20 , @xmath5 and @xmath10 using the three boundary conditions as discussed above ; ( iii ) calculate the enthalpy everywhere inside the configuration using eq . ( [ bernoulli ] ) ; ( iv ) calculate the new  trial " density distribution according to the relationship between the density and enthalpy for a polytrope . these steps are repeated until the model converges . in this section , we show our results for incompressible riemann s - type ellipsoids and their quasi - equilibrium , compressible counterparts . since the available parameter space occupies a large portion of figure [ triaxial ] , we choose to build models only along the vertical dashed lines , which largely represent most of the configurations . our models are constructed on a cylindrical grid with a resolution of @xmath49 in the @xmath50 ( cylindrical radius ) , @xmath51 ( vertical ) , and @xmath52 ( azimuthal ) directions , respectively . following @xcite , we adopt a set of polytropic units in which the gravitational constant @xmath53 , the radius of the entire grid @xmath54 , and maximum density @xmath55 are all set to 1 . for all of our models , we set the semiaxis @xmath2 to @xmath56 . we choose to present our results in such a way that readers can readily normalize our data to any unit system they are familiar with for the purpose of comparison . table [ n0 ] shows the results for direct configurations of incompressible ( @xmath47 ) riemann s - type ellipsoids . tables [ n05 ] , [ n10 ] , and [ n15 ] show the same results but for compressible counterparts of riemann ellipsoids with @xmath57 , @xmath58 , and @xmath59 , respectively . tables [ n0a ] , [ n05a ] , and [ n10a ] show results for adjoint configurations with @xmath60 , @xmath57 and @xmath58 , respectively .    in these tables , @xmath61 is the rotational energy of the system , @xmath62 is the total mass , @xmath63 is the total angular momentum , @xmath64 is the mean density , @xmath65 is the gravitational potential energy , and @xmath66 is the thermal energy . the quantity @xmath67 measures how well our models are in virial equilibrium ; ideally , this ratio should have a value of 0.5 . all of our models have values that are very close to 0.5 , especially for compressible configurations . in table [ n0 ] , the parameters @xmath68 and @xmath69 are analytically computed results for @xmath5 and @xmath10 , respectively . several properties of our models are worth discussing . first , for direct configurations with a fixed value of @xmath0 , as one proceeds from larger values of @xmath1 to smaller values , @xmath10 starts from a positive value , then continuously decreases ; it reaches zero when the rigidly rotating jacobi sequence is hit , and then becomes negative after passing the jacobi sequence . the parameter @xmath45 also switches sign when crossing the jacobi / dedekind sequences . for adjoint configurations , a similar behavior is observed for @xmath5 , while @xmath45 approaches @xmath70 around the dedekind sequence because @xmath71 . second , in the parameter space above the jacobi / dedekind sequence where @xmath72 , @xmath5 is positive for direct configurations , and negative for adjoint configurations , which reflects the fact that the jacobi mode in maclaurin spheroids is the `` forward moving '' wave but the dedekind mode is the `` backward moving '' wave . this also implies that in the adjoint configurations with @xmath72 , the ellipsoidal pattern is moving in a direction that is opposite to the overall rotation ! furthermore , when the critical jacobi / dedekind sequence is passed , @xmath45 switches sign ; so does @xmath10 in direct configurations and @xmath5 in adjoint configurations . in maclaurin spheroids , this is where the frequency of the backward moving wave becomes positive in the inertial frame ; hence , it is where the secular bar - mode instability sets in through gravitational radiation ( chandrasekhar - friedman - schutz instability ) or viscosity @xcite . finally , for softer equations of state , the scf technique failed to converge to models with relatively small values of @xmath1 because , in these cases , the centrifugal force at the boundary point a ( see figure [ boundary ] ) would exceed the gravitational force . in fact , for @xmath73 and @xmath59 , even the jacobi sequence is not reached . this is consistent with previous studies that have indicated that the jacobi sequence can not exist for @xmath74 @xcite . furthermore , as was pointed out by @xcite in their analysis of uniformly rotating polytropes , the mass - shedding limit never extends much beyond the bifurcation point ( where the jacobi / dedekind sequence branches off the maclaurin spheriod sequence ) even for very small @xmath75 : for @xmath76 , this mass - shedding limit corresponds to @xmath77 . our @xmath78 models are consistent with their results . therefore , the available parameter space for compressible counterparts of riemann ellipsoids becomes smaller as the equation of state becomes softer . figures [ direct05 ] and [ adjoint05 ] show equatorial - plane iso - density contours and velocity fields for direct and adjoint configurations of our @xmath57 model with @xmath79 and @xmath80 . both the rotating frame velocity field @xmath11 and inertial frame velocity field @xmath81 are plotted for comparison . the flow patterns in these two models are different from each other although they have identical geometrical shapes . it is also clear that the fluid motion in the adjoint configuration is dominated by internal motion rather than rotation . as another illustration , figures [ direct ] and [ adjoint ] illustrate the structure of our @xmath73 model with @xmath79 and @xmath82 . notice that @xmath5 , which is also the pattern speed of the ellipsoid , is negative in both fig . [ adjoint05 ] and fig . [ adjoint ] , so the ellipsoidal figures in both plots are spinning retrograde ( clockwise ) relative to the overall fluid rotation . comparisons between previous published results and our models yield very good agreement . for example , from table iv in chapter 6 of @xcite , a jacobi ellipsoid with @xmath83 and @xmath0=0.60 , @xmath84 , has @xmath85 and @xmath86 , where @xmath87 . the closest match among our models shown in table [ n0 ] , with @xmath79 , @xmath88 , has @xmath89 , @xmath90 , and an almost vanishing @xmath91 , which shows that this model from our table is very close to the jacobi sequence . it should be noticed that , due to the discrete nature of our grid , we can not construct models with axis ratios that will place the models precisely on the jacobi / dedekind sequence , but we can construct models that are very close to this sequence ( with very small values of @xmath10 or @xmath5 ) . as another comparison , the last four columns in table [ n0 ] give analytically computed results of @xmath92 and @xmath93 for riemann s - type ellipsoids having the specified axes ratios @xmath0 and @xmath1 , as well as the ratio between numerical and analytical results ( only models with @xmath94 have been computed ) . to obtain these analytical results , we used standard incomplete elliptic integrals to evaluate the potential field of a homogeneous ellipsoid ( * ? ? ? * see their table 2 - 2 ) . in virtually all cases , the error in @xmath5 is @xmath95 a few percent . on the other hand , for the direct configurations the fractional error in @xmath10 can be as large as 50% ; the error is largest for models whose @xmath10 almost vanishes . this is understandable because a small variation in the magnitude of @xmath10 will result in a large percentage difference when @xmath10 is approaching zero . also some mismatch between numerical and analytical results may arise because the discrete nature of our grids makes it difficult for us to define an exact axis ratio for our models .    for compressible equations of state , we can compare our results with those presented by @xcite , who used an energy variational method to determine approximate equilibrium configurations for ellipsoidal self - gravitating systems . @xcite assumed that the iso - density contours of each configuration are self - similar concentric ellipses , which as we have already emphasized is not the case . in their table 4 , for compressible analogues of jacobi ellipsoids , a configuration with @xmath96 and @xmath97 has @xmath98 ; the closest match among our @xmath57 models ( see table [ n05 ] ) with @xmath99 and @xmath100 has @xmath101 . other quantities also match well , after a proper normalization is applied as prescribed by their equation ( 3.27 ) . we omit more comparisons here and leave the rest to interested readers . other models that are of interest are the @xmath102 irrotational sequence @xcite . it turns out that the @xmath45 parameter is very sensitive to the axis ratios , so we could not build models that are very close to the irrotational sequence due to the discrete nature of our computational grids . our model that comes closest has @xmath103 ; it is a direct configuration with @xmath57 and its other parameters are listed in table [ n05 ] . we have presented a new method to construct 3d models for riemann s - type ellipsoids ( @xmath47 ) , and the method has been extended to construct compressible counterparts that share the same velocity field as that of riemann s - type ellipsoids . with this method , we have been able to build 3d models that cover almost the entire parameter space of riemann s - type ellipsoids . our results are in good agreement with previous studies . we expect that this method can also be straightforwardly extended to the construction of roche - riemann ellipsoids .    by design , our compressible models satisfy the steady - state euler s equation . however , it is not guaranteed that they simultaneously satisfy the steady - state continuity equation exactly , because this additional constraint was not built into our method . because the adopted velocity flow - field @xmath11 in our models is divergence - free , the violations of the steady - state continuity equation will only arise to the extent that the iso - density contours deviate from concentric ellipses . as we can see from figures [ direct05 ] - [ adjoint ] , the outer low - density shell in compressible models becomes somewhat pointed along the major axis . however , the deviation from concentric ellipses is small and almost negligible in the inner , high - density regions . furthermore , a comparison between fig . [ direct05 ] and fig . [ direct ] shows that this deviation becomes smaller in less compressible models . in order to quantitatively measure the violation of the steady - state continuity equation , we calculated @xmath104 for the models shown in figs [ direct05 ] , [ adjoint05 ] , and [ direct ] . as long as @xmath105 for a given configuration , we can consider that the model has a reasonablly good quasi - equilibrium structure . in figure [ divrhov ] , we have plotted @xmath106 in the equatorial plane as a function of cylindrical radius @xmath50 along different azimuthal angles , @xmath107 and @xmath108 . in fact , this quantity remains small ( @xmath109 for @xmath57 and @xmath110 for @xmath58 ) throughout the majority of the configuration , but not surprisingly near the surface it asymptotically approaches infinity because @xmath111 goes to zero . as the figure illustrates , @xmath106 is smaller for the stiffer equation of state . we also note that @xmath106 is nearly zero along the x and y axes because the velocity vectors in the rotating frame along the x ( y ) axis only have a nonzero y ( x ) component and are therefore perpendicular to the density gradient there . thus , we expect our compressible models to represent fairly good quasi - equilibrium configurations , especially in nearly incompressible cases . on the other hand , this point was brought to our attention by an anonymous referee : the relatively large violation of the steady - state continuity equation on the very surface layer can potentially change the shape of the surface in hydrodynamical evolutions . to understand how the shape of the surface might be changed by the flow field , it is important to know the sign of @xmath112 along the surface . since we are adopting a divergence free velocity field , the sign of @xmath112 is determined by the sign of @xmath113 , which also measures to what degree the velocity vectors are orthogonal to density gradient . if the velocity vector points slightly outward ( inward ) from the surface , @xmath113 will be negative ( positive ) , because the density gradient is negative outward . [ divrhov2 ] and [ divrhov3 ] show contour maps of @xmath112 in the equatorial planes of the direct and adjoint models shown in figs . [ direct05 ] and [ adjoint05 ] , respectively . ( we note that the contour pattern for direct configurations will switch from fig . [ divrhov2 ] to fig . [ divrhov3 ] if @xmath10 flips sign . ) a @xmath114 symmetry is observed for both patterns . we also note that , in each quadrant , @xmath112 flips sign , which implies that the direction of velocity vectors will change from inward ( outward ) to outward ( inward ) from the surface within one single quadrant . accordingly , fig . [ shape ] shows a combined schematic drawing of how the velocity vectors in the rotating frames are aligned with the surface for these two patterns ( the panel on the left is for the adjoint figure shown in fig . [ adjoint05 ] , the panel on the right is for the direct figure shown in fig . [ direct05 ] ) . in this figure , the black ellipse denotes the actual surface ( which should not be a perfect ellipse ) , red arrows denote velocity vectors on the surface , purple arrows denote density gradient vectors , which are drawn pointing inward because of their negative sign . for direct configurations , in which the figure rotation is larger than the internal motion , these mismatches ( nonvanishing @xmath113 ) on top of a large rotation are likely to generate spiral arms ; whereas , for adjoint configurations , the internal motion is larger than figure rotation , so these mismatches tend to generate waves along the surface .    using the 3d hydrodynamics code developed by the lsu group @xcite , we have carried out nonlinear 3d hydrodynamical evolutions for a subset of the nearly incompressible models discussed here in table [ n05 ] . despite the fact that the models do not fully satisfy the steady - state continuity equation , they are in very good force - balance initially and , as predicted by the above discussion , the models are fairly long - lived , holding their initial ellipsoidal shapes over the time that is followed by our simulations ( @xmath115 15 - 20 dynamical times ) . we also observed the formation of narrow spiral arms ( low amplitude surface waves ) containing low density materials in direct ( adjoint ) configurations . the results of these hydrodynamical evolutions will be presented in a separate paper in the context of our ongoing study of the stability of riemann s - type ellipsoids . since our ellipsoidal models are fairly long - lived and have exactly the same elliptical flow field as riemann s - type ellipsoids , they should provide good initial states for a nonlinear , hydrodynamical study of dynamical instabilities in riemann s - type ellipsoids . we especially have in mind the elliptical strain instability that is associated with elliptical stream lines , which should set in on a dynamical time scale . through such a study , we hope to provide insight into the stability of ellipsoidal configurations that will complement the work of @xcite . this should also help answer whether the evolutionary path of a grr - driven star toward the dedekind sequence is viable . i would like to thank joel tohline for valuable suggestions , discussions and carefully reading through this paper . i also thank juhan frank and patrick motl for helpful comments . it is a pleasure to thank norman lebovitz for drawing our attention to the particular analytical expression of the velocity field of riemann s - type ellipsoids used in this study . i would also thank howard cohl for developing the poisson solver that was an important ingredient in this work . i thank the anonymous referee for insightful suggestions to analysis our models . this work was partially supported by nsf grants ast-0407070 and phy-0326311 , and by the center for computation & technology at lsu . binney , j. , & tremaine , s. , 1987 , galactic dynamics , new york : academic press cazes , j. e. , & tohline , j. e. 2000 , , 532 , 1051 chandrasekhar , s. 1969 , equilibrium figures of equilibrium , new haven , ct : yale univ . press    cohl , h. s. , & tohline , j. e. 1999 , , 527 , 86 friedman , j. , & schutz , b. f. 1978 , , 222 , 281 hachisu , i. 1986 , , 61 , 479 hachisu , i. 1986 , , 62 , 461 hachisu , i. , eriguchi , y. , 1982 , prog . theor.phys . , 68 , 206 hachisu , i. , eriguchi , y. , 1984 , , 36 , 239 james , r. a. , 1964 , , 140 , 552      lebovitz , n. r. , & lifschitz , a. 1996 , , 458 , 699 lebovitz , n. r. , & saldanha , k. i. 1999 , physics of fluids , 11 , 3374 lifschitz , a. , & lebovitz , n. 1993 , , 408 , 603 motl , p. m. , tohline , j. e. , & frank , j. 2002 , , 138 , 121                    0.90 & 0.923 & 1.030 & 0.7550 & 6.6e-3 & 0.03183 & 0.8232 & -0.6989 & 0.4989 & -1.5 & & + 0.90 & 0.795 & 1.144 & 0.4310 & 0.046 & 0.07018 & 0.7086 & -0.5429 & 0.4992 & -0.76 & 1.14704 & 0.997 & 0.43181 & 0.998 + 0.90 & 0.641 & 1.119 & 0.1394 & 0.101 & 0.07764 & 0.5720 & -0.3759 & 0.4988 & -0.25 & 1.13137 & 0.989&0.15077 & 0.925 + 0.90 & 0.590 & 1.099 & 0.0569 & 0.122 & 0.07589 & 0.5261 & -0.3249 & 0.4987 & -0.10 & 1.10661 & 0.993 & 0.06406 & 0.888 + 0.90 & 0.564 & 1.090 & 0.0210 & 0.132 & 0.07452 & 0.5033 & -0.3007 & 0.4989 & -0.04 & 1.09034 & 0.999&0.02033&1.033 + 0.90 & 0.538 & 1.069 & -0.0250 & 0.144 & 0.07268 & 0.4801 & -0.2768 & 0.4991 & 0.05 & 1.07148 & 0.998&-0.02324 & 1.075 + 0.90 & 0.487 & 1.015 & -0.1205 & 0.167 & 0.06814 & 0.4341 & -0.2316 & 0.4992 & 0.24 & 1.02639&0.989 & -0.10880 & 1.108 + 0.90 & 0.333 & 0.797 & -0.3885 & 0.246 & 0.04860 & 0.2970 & -0.1167 & 0.4983 & 0.98 & 0.79257 & 1.006&-0.39224&0.990 + 0.74 & 0.769 & 1.099 & 0.5604 & 0.030 & 0.03807 & 0.5682 & -0.3744 & 0.4982 & -1.1 & & + 0.74 & 0.692 & 1.133 & 0.3840 & 0.058 & 0.04678 & 0.5116 & -0.3134 & 0.4986 & -0.71 & 1.13215 & 1.001&0.38599&0.995 + 0.74 & 0.538 & 1.094 & 0.0774 & 0.121 & 0.04836 & 0.3982 & -0.2034 & 0.4990 & -0.15 & 1.09508 & 0.999&0.08149 & 0.950 + 0.74 & 0.513 & 1.078 & 0.0303 & 0.133 & 0.04734 & 0.3791 & -0.1857 & 0.4991 & -0.06 & 1.07883 & 0.999&0.03474&0.872 + 0.74 & 0.487 & 1.054 & -0.0208 & 0.145 & 0.04594 & 0.3560 & -0.1700 & 0.4988 & 0.04 & 1.05883 & 0.995 & -0.01363 & 1.526 + 0.74 & 0.436 & 1.003 & -0.1184 & 0.171 & 0.04284 & 0.3223 & -0.1400 & 0.4987 & 0.25 & 1.00963 & 0.993&-0.10852 & 1.091 + 0.74 & 0.333 & 0.847 & -0.3230 & 0.228 & 0.03389 & 0.2463 & -0.0863 & 0.4983 & 0.80 & 0.85686 & 0.988&-0.31006&1.042 + 0.74 & 0.282 & 0.715 & -0.4523 & 0.259 & 0.02842 & 0.2081 & -0.0634 & 0.4968 & 1.3 & 0.72756 & 0.983&-0.43781 & 1.033 +    0.59 & 0.718 & 0.948 & 0.7677 & 0.014 & 0.01197 & 0.4208 & -0.2243 & 0.4924 & -1.9 & & + 0.59 & 0.538 & 1.087 & 0.2338 & 0.093 & 0.02876 & 0.3157 & -0.1374 & 0.4982 & -0.49 & 1.08732 & 0.999 & 0.23529 & 0.994 + 0.59 & 0.436 & 1.029 & 0.0108 & 0.145 & 0.02691 & 0.2554 & -0.0949 & 0.4985 & -0.02 & 1.03198&0.997 & 0.01472&0.734 + 0.59 & 0.410 & 1.008 & -0.0417 & 0.159 & 0.02592 & 0.2404 & -0.0853 & 0.4991 & 0.09 & 1.00774&1.000&-0.04012&1.039 + 0.59 & 0.359 & 0.943 & -0.1534 & 0.188 & 0.02342 & 0.2104 & -0.0673 & 0.4989 & 0.37 & 0.94638 & 0.996 & -0.14830&1.034 + 0.59 & 0.231 & 0.584 & -0.5426 & 0.267 & 0.01475 & 0.1351 & -0.0300 & 0.4910 & 2.1 & 0.60698 & 0.962&-0.51563&1.013 +    0.41 & 0.385 & 0.971 & 0.1369 & 0.132 & 0.01220 & 0.1563 & -0.0407 & 0.4978 & -0.40 & 0.97108 & 1.000&0.14159&0.967 + 0.41 & 0.333 & 0.930 & 0.0016 & 0.165 & 0.01112 & 0.1347 & -0.0313 & 0.4998 & -0.005 & 0.92963 & 1.000 & 0.00331&0.483 + 0.41 & 0.308 & 0.900 & -0.0685 & 0.183 & 0.01061 & 0.1254 & -0.0275 & 0.4991 & 0.22 & 0.89982&1.000 & -0.06224&1.101 + 0.41 & 0.256 & 0.808 & -0.2074 & 0.217 & 0.00896 & 0.1046 & -0.0199 & 0.4937 & 0.73 & 0.81275&0.994 & -0.20138&1.030 + 0.28 & 0.256 & 0.809 & 0.0246 & 0.183 & 0.00472 & 0.0717 & -0.0103 & 0.4982&-0.12 & 0.80944&0.999 & 0.03668&0.671 + 0.28 & 0.231 & 0.780 & -0.0572 & 0.204 & 0.00431 & 0.0646 & -0.0086 & 0.4988&0.28 & 0.77651&1.004&-0.04714&1.213 + 0.28 & 0.205 & 0.726 & -0.1505 & 0.224 & 0.00382 & 0.0575 & -0.0069 & 0.4967&0.80 & 0.72853&0.996&-0.13511&1.114 +      0.90 & 0.795 & 0.970 & 0.3660 & 0.045 & 0.0258 & 0.3812 & 0.5398 & -0.1751 & 0.5000 & -0.76 + 0.90 & 0.744 & 0.975 & 0.2831 & 0.061 & 0.0273 & 0.3541 & 0.5371 & -0.1544 & 0.5000 & -0.58 + 0.90 & 0.692 & 0.970 & 0.2050 & 0.079 & 0.0276 & 0.3266 & 0.5344 & -0.1344 & 0.5001&-0.43 + 0.90 & 0.641 & 0.958 & 0.1312 & 0.097 & 0.0227 & 0.2987 & 0.5309 & -0.1150 & 0.5001 & -0.28 + 0.90 & 0.590 & 0.940 & 0.0617 & 0.117 & 0.0256 & 0.2695 & 0.5263 & -0.0964 & 0.5001 & -0.13 + 0.90 & 0.564 & 0.929 & 0.0287 & 0.127 & 0.0246 & 0.2546 & 0.5240 & -0.0873 & 0.5001 & -0.06 + * 0.90 & 0.538 & 0.920 & 0.0001 & 0.137 & 0.0233 & 0.2391 & 0.5198 & -0.0783 & 0.5001 & -0.0002 + 0.90 & 0.513 & 0.908 & -0.0281 & 0.147 & 0.0218 & 0.2234 & 0.5169 & -0.0695 & 0.5001 & 0.06 + 0.90 & 0.487 & 0.898 & -0.0513 & 0.157 & 0.0201 & 0.2066 & 0.5115 & -0.0607 & 0.5001 & 0.12 + 0.74 & 0.692 & 0.958 & 0.3232 & 0.055 & 0.0166 & 0.2719 & 0.5380 & -0.0991 & 0.4997 & -0.71 + 0.74 & 0.641 & 0.959 & 0.2343 & 0.074 & 0.0171 & 0.2489 & 0.5345 & -0.0852 & 0.4999 & -0.51 + 0.74 & 0.590 & 0.949 & 0.1514 & 0.094 & 0.0167 & 0.2252 & 0.5314 & -0.0717 & 0.5001 & -0.33 + 0.74 & 0.538 & 0.932 & 0.0750 & 0.114 & 0.0156 & 0.2005 & 0.5263 & -0.0587 & 0.5001&-0.17 + * 0.74 & 0.487 & 0.911 & 0.0081 & 0.135 & 0.0137 & 0.1737 & 0.5192 & -0.0458 & 0.5002 & -0.02 + 0.74 & 0.462 & 0.904 & -0.0171 & 0.144 & 0.0123 & 0.1588 & 0.5156 & -0.0393 & 0.5002 & 0.04 + 0.59 & 0.718 & 0.792 & 0.6631 & 0.012 & 0.0037 & 0.2241 & 0.5418 & -0.0713 & 0.4943 & -1.9 + 0.59 & 0.538 & 0.922 & 0.2021 & 0.085 & 0.0090 & 0.1587 & 0.5304 & -0.0397 & 0.4997 & -0.50 + 0.59 & 0.487 & 0.911 & 0.1137 & 0.105 & 0.0081 & 0.1373 & 0.5256 & -0.0313 & 0.5001 & -0.29 + 0.59 & 0.462 & 0.905 & 0.0754 & 0.114 & 0.0073 & 0.1250 & 0.5230 & -0.0265 & 0.5002 & -0.19 + 0.59 & 0.436 & 0.871 & -0.0337 & 0.153 & 0.0105 & 0.1382 & 0.5453 & -0.0031 & 0.5088 & 0.09 +      0.90 & 0.795 & 0.803 & 0.3029 & 0.039 & 0.00893 & 0.2053 & 0.2921 & -0.0577 & 0.5001 & -0.76 + 0.90 & 0.744 & 0.804 & 0.2329 & 0.053 & 0.00910 & 0.1870 & 0.2867 & -0.0492 & 0.5001 & -0.58 + 0.90 & 0.692 & 0.796 & 0.1684 & 0.067 & 0.00872 & 0.1679 & 0.2799 & -0.0409 & 0.5002 & -0.43 + 0.90 & 0.641 & 0.782 & 0.1090 & 0.081 & 0.00788 & 0.1476 & 0.2711 & -0.0328 & 0.5002 & -0.28 + 0.90 & 0.590 & 0.760 & 0.0555 & 0.093 & 0.00657 & 0.1255 & 0.2611 & -0.0249 & 0.5002 & -0.14 + 0.90 & 0.564 & 0.747 & 0.0314 & 0.097 & 0.00572 & 0.1134 & 0.2563 & -0.0210 & 0.5002 & -0.08 + 0.74 & 0.641 & 0.784 & 0.1917 & 0.060 & 0.00513 & 0.1262 & 0.2805 & -0.0254 & 0.5003 & -0.51 + 0.74 & 0.590 & 0.769 & 0.1241 & 0.074 & 0.00450 & 0.1087 & 0.2731 & -0.0197 & 0.5003 & -0.34 + 0.74 & 0.538 & 0.743 & 0.0626 & 0.084 & 0.00349 & 0.0893 & 0.2680 & -0.0143 & 0.5003 & -0.18 + 0.74 & 0.512 & 0.722 & 0.0322 & 0.086 & 0.00289 & 0.0786 & 0.2676 & -0.0115 & 0.5003 & -0.09 + 0.59 & 0.590 & 0.748 & 0.244 & 0.046 & 0.00237 & 0.0870 & 0.2836 & -0.0136 & 0.5007 & -0.75 + 0.59 & 0.564 & 0.743 & 0.2011 & 0.052 & 0.00220 & 0.0797 & 0.2817 & -0.0118 & 0.5009 & -0.62 +      0.90 & 0.795 & 0.631 & 0.2351 & 0.030 & 0.00273 & 0.1068 & 0.1531 & -0.0181 & 0.5002&-0.75 + 0.90 & 0.744 & 0.622 & 0.1752 & 0.040 & 0.00260 & 0.0945 & 0.1467 & -0.0147 & 0.5003&-0.57 + 0.90 & 0.692 & 0.604 & 0.1187 & 0.048 & 0.00229 & 0.0817 & 0.1397 & -0.0115 & 0.5003&-0.40 + 0.90 & 0.641 & 0.574 & 0.0640 & 0.055 & 0.00183 & 0.0683 & 0.1324 & -0.0085 & 0.5003&-0.22 + 0.90 & 0.615 & 0.554 & 0.0365 & 0.057 & 0.00157 & 0.0614 & 0.1292 & -0.00711 & 0.5003&-0.13 + 0.74 & 0.692 & 0.609 & 0.2024 & 0.033 & 0.00148 & 0.0717 & 0.1491 & -0.00931 & 0.5006&-0.70 + 0.74 & 0.641 & 0.592 & 0.1378 & 0.041 & 0.00130 & 0.0612 & 0.1430 & -0.00713 & 0.5005&-0.49 + 0.64 & 0.670 & 0.587 & 0.2503 & 0.023 & 0.00086 & 0.0586 & 0.1524 & -0.00666 & 0.5014&-0.94 + 0.64 & 0.640 & 0.584 & 0.2116 & 0.027 & 0.00084 & 0.0542 & 0.1499 & -0.00584 & 0.5013&-0.80 +      0.90 & 0.923 & -0.755 & -1.030 & 0.0064 & 0.0307 & 0.8232 & 0.99999 & -0.6989 & 0.4987 & -2.744 + 0.90 & 0.641 & -0.139 & -1.119 & 0.100 & 0.0770 & 0.572 & 0.99999 & -0.3758 & 0.4980 & -16.1 + 0.90 & 0.590 & -0.0569 & -1.099 & 0.121 & 0.0754 & 0.5261 & 0.99999 & -0.3249 & 0.4979 & -38.8 + 0.90 & 0.564 & -0.0211 & -1.090 & 0.132 & 0.0741 & 0.5033 & 0.99999 & -0.3007 & 0.4981 & -104.2 + 0.90 & 0.538 & 0.0250 & -1.069 & 0.143 & 0.0723 & 0.4801 & 0.99999 & -0.2768 & 0.4982 & 86.0 + 0.90 & 0.487 & 0.121 & -1.015 & 0.166 & 0.0678 & 0.4342 & 0.99999 & -0.2314 & 0.4983 & 16.9 + 0.90 & 0.436 & 0.219 & -0.945 & 0.191 & 0.0623 & 0.3886 & 0.99999 & -0.1900 & 0.4979 & 8.689 + 0.90 & 0.333 & 0.389 & -0.797 & 0.246 & 0.0485 & 0.2970 & 0.99999 & -0.1167 & 0.4978 & 4.128 + 0.74 & 0.744 & -0.495 & -1.117 & 0.035 & 0.0376 & 0.5496 & 0.99999 & -0.3540 & 0.4943 & -4.7 + 0.74 & 0.590 & -0.171 & -1.116 & 0.093 & 0.0465 & 0.4357 & 0.99999 & -0.2377 & 0.4933 & -13.6 + 0.74 & 0.513 & -0.030 & -1.078 & 0.127 & 0.0452 & 0.3791 & 0.99999 & -0.1866 & 0.4932 & -74.3 + 0.74 & 0.487 & 0.021 & -1.054 & 0.139 & 0.0440 & 0.3595 & 0.99999 & -0.1699 & 0.4928 & 106.8 + 0.74 & 0.436 & 0.118 & -1.002 & 0.165 & 0.0414 & 0.3223 & 0.99999 & -0.1400 & 0.4929 & 17.7 + 0.74 & 0.385 & 0.210 & -0.944 & 0.193 & 0.0376 & 0.2841 & 0.99999 & -0.1117 & 0.4934 & 9.4 + 0.74 & 0.333 & 0.323 & -0.847 & 0.223 & 0.0332 & 0.2463 & 0.99999 & -0.0863 & 0.4940 & 5.5 + 0.74 & 0.282 & 0.453 & -0.714 & 0.257 & 0.0281 & 0.2081 & 0.99999 & -0.0634 & 0.4944 & 3.4 + 0.59 & 0.615 & -0.417 & -1.085 & 0.045 & 0.0199 & 0.3608 & 0.99999 & -0.1729 & 0.4834 & -5.5 + 0.59 & 0.487 & -0.119 & -1.063 & 0.100 & 0.0239 & 0.2853 & 0.99999 & -0.1152 & 0.4807 & -20.4 + 0.59 & 0.436 & -0.011 & -1.029 & 0.126 & 0.0235 & 0.2553 & 0.99999 & -0.0949 & 0.4801 & -217.5 + 0.59 & 0.410 & 0.042 & -1.008 & 0.140 & 0.0229 & 0.2404 & 0.99999 & -0.0853 & 0.4807 & 55.2 + 0.59 & 0.359 & 0.153 & -0.943 & 0.170 & 0.0213 & 0.2104 & 0.99999 & -0.0672 & 0.4811 & 14.1 + 0.59 & 0.282 & 0.271 & -0.856 & 0.204 & 0.0190 & 0.1803 & 0.99999 & -0.0509 & 0.4822 & 7.2 + 0.59 & 0.231 & 0.402 & -0.729 & 0.242 & 0.0161 & 0.1502 & 0.99999 & -0.0365 & 0.4843 & 4.1 +      0.90 & 0.923 & -0.642 & -0.8766 & 0.0063 & 0.0116 & 0.4476 & 0.5442 & -0.2297 & 0.4997 & -2.7 + 0.90 & 0.590 & -0.062 & -0.9401 & 0.116 & 0.0255 & 0.2695 & 0.5263 & -0.0964 & 0.4997 & -30.7 + 0.90 & 0.538 & -1.47e-4&-0.9195 & 0.137 & 0.0232 & 0.2391 & 0.5197 & -0.0782 & 0.4999 & -12606.1 + 0.90 & 0.513 & 0.028 & -0.9080 & 0.147 & 0.0218 & 0.2234 & 0.5169 & -0.0695 & 0.5000 & 65.2 + 0.74 & 0.795 & -0.537 & -0.9101 & 0.018 & 0.0102 & 0.3167 & 0.5421 & -0.1283 & 0.4966 & -3.5 + 0.74 & 0.692 & -0.323 & -0.9584 & 0.052 & 0.0154 & 0.2719 & 0.5379 & -0.0991 & 0.4963 & -6.2 + 0.74 & 0.538 & -0.075 & -0.9317 & 0.112 & 0.0150 & 0.2005 & 0.5263 & -0.0587 & 0.4975 & -25.9 + 0.74 & 0.487 & -8.2e-3 & -0.9114 & 0.134 & 0.0133 & 0.1737 & 0.5192 & -0.0458 & 0.4992 & -233.1 + 0.74 & 0.462 & 0.0171 & -0.9044 & 0.145 & 0.0121 & 0.1587 & 0.5157 & -0.0393 & 0.5009 & 110.2 + 0.59 & 0.615 & -0.355 & -0.9180 & 0.044 & 6.93e-3 & 0.1877 & 0.5368 & -0.0529 & 0.4894 & -5.9 + 0.59 & 0.538 & -0.202 & -0.9222 & 0.076 & 7.66e-3 & 0.1587 & 0.5304 & -0.0397 & 0.4910 & -10.4 + 0.59 & 0.487 & -0.114 & -0.9111 & 0.100 & 7.21e-3 & 0.1373 & 0.5256 & -0.0310 & 0.4948 & -18.3 + 0.59 & 0.462 & -0.075 & -0.9050 & 0.113 & 6.65e-3 & 0.1250 & 0.5230 & -0.0265 & 0.4992 & -27.4 +      0.90 & 0.795 & -0.303 & -0.8030 & 0.039 & 0.00887 & 0.2053 & 0.2921 & -0.0577 & 0.5000 & -5.33 + 0.90 & 0.692 & -0.168 & -0.7963 & 0.067 & 0.00869 & 0.1679 & 0.2799 & -0.0409 & 0.5003 & -9.51 + 0.90 & 0.590 & -0.056 & -0.7605 & 0.094 & 0.00657 & 0.1255 & 0.2611 & -0.0249 & 0.5009 & -27.59 + 0.90 & 0.564 & -0.031 & -0.7469 & 0.099 & 0.00573 & 0.1134 & 0.2563 & -0.0210 & 0.5014 & -47.78 + 0.74 & 0.641 & -0.192 & -0.7842 & 0.061 & 0.00499 & 0.1262 & 0.2805 & -0.0254 & 0.5013 & -8.54 + 0.74 & 0.538 & -0.062 & -0.7429 & 0.091 & 0.00354 & 0.0893 & 0.2680 & -0.0142 & 0.5071 & -24.80 + 0.74 & 0.512 & -0.032 & -0.7220 & 0.096 & 0.00295 & 0.0786 & 0.2676 & -0.0115 & 0.5097 & -46.74 + 0.59 & 0.590 & -0.244 & -0.7480 & 0.050 & 0.00217 & 0.0870 & 0.2836 & -0.0136 & 0.5049 & -7.01 + 0.59 & 0.564 & -0.201 & -0.7433 & 0.059 & 0.00210 & 0.0797 & 0.2817 & -0.0118 & 0.5083 & -8.47 +
Abstract:

Answer: we introduce a new technique for constructing three - dimensional ( 3d ) models of incompressible riemann s - type ellipsoids and compressible triaxial configurations that share the same velocity field as that of riemann s - type ellipsoids . our incompressible models are exact steady - state configurations ; our compressible models represent approximate steady - state equilibrium configurations . models built from this method can be used to study a variety of relevant astrophysical and geophysical problems .