Instruction: Write an abstract for the following acadamic paper.

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

Paper: in this paper we study the effects of rotation on the radiation driven winds of early - type supergiants . we will focus on the explanation for the occurrence of disks around fast rotating b[e ] supergiants . b[e ] supergiants , also designated sgb[e ] stars ( lamers et al . @xcite ) , are b type supergiants that exhibit forbidden emission lines in their optical spectra . the observations of hybrid spectra , i.e. spectra with _ broad _ uv p cygni features and _ narrow _ emission lines and dust emission of b[e ] supergiants in the magellanic clouds led to the proposal of a disk wind model by zickgraf et al.(@xcite ) . this model postulates a dense disk of outflowing material in a fast line driven wind to explain the observed characteristics of the spectra of these stars ( see zickgraf  @xcite ) . the stellar wind at the equator is about ten times slower than that at the pole . also , the wind at the equator is about a hundred times denser than at the pole . additional evidence for a two - component outflow has been obtained from polarimetric measurements , e.g. by zickgraf & schulte - ladbeck  ( @xcite ) . the precise mechanism behind the formation of these disks is still a mystery ( cassinelli 1998 ) . it is however clear that the origin of an axisymmetric wind structure such as a disk may well be connected to the fast rotation of a star . two theories using rotation in a different way have been considered for the formation of these disk winds :     ( 1 ) the wind compressed disk ( wcd ) model of bjorkman & cassinelli  ( @xcite ) , that invokes the kinematics of the winds from rotating stars . the streamlines of the gas in the wind from both hemispheres of a rapidly rotating star cross in the equatorial plane . the concentration of the gas and the shock in the equatorial plane produce an outflowing equatorial disk , with a thickness on the order of a few degrees . owocki & cranmer ( 1994 ) have argued that the motion to the equatorial plane may be counteracted by the radiation force perpendicular from the plane due to lines . however , this `` wind - inhibition effect '' may not be effective in winds with a strong density gradient in the equatorial direction . we return to this in sect . ( 2 ) the rotationally induced bi - stability model ( rib ) of lamers & pauldrach   ( @xcite ) invokes an increase in the mass flux from the equator and a decrease in the equatorial wind velocity compared to the poles from the bi - stability jump . this jump in mass flux and in wind velocity is due to the temperature difference between the pole and the equator of a fast rotating b[e ] supergiant with gravity darkening . the jump will occur for stars with effective temperatures between 20  000 and 30  000 k.    in reality , both effects , i.e. the wind compression and the rotation induced bi - stability , may be operating together and amplifying one another in the line driven winds of rapidly rotating early - b stars . ( for a detailed explanation of both models , see lamers & cassinelli ( 1999 ) ( hereafter @xmath5 ) , chapter 11 . ) the theory of radiation driven winds for non - rotating stars was developed by castor et al . ( @xcite ) ( hereafter cak ) and predicts the mass loss rate and @xmath6 for spherical winds . the influence of rotation on line driven winds was investigated by friend & abbott ( @xcite ) . they found an increase in the mass - loss rate at the equator and a corresponding decrease in @xmath6 , but these authors did not consider the effects of gravity darkening and the oblate shape of the star . cranmer & owocki  ( @xcite ) described the finite disk correction factor in case of an oblate star . in this paper we will modify the line driven wind theory , including these effects of oblateness and gravity darkening as well as the rotational terms in the equation of motion . this will then be applied to rapidly rotating b[e ] supergiants to investigate whether these effects can explain the occurrence of disks . we will also apply recent bi - stability calculations to study the effect of the rib - model to explain the occurrence of outflowing disks around rapidly rotating b[e ] supergiants .    in sect . 2 we describe the theoretical background of the radiation driven wind theory and the bi - stability effect . in sect . 3 we derive the equations for the winds from oblate rotating stars with a temperature gradient between poles and equator due to gravity darkening . in sect . 4 we describe the method for solving the equations and for calculating the winds from rotating stars . these calculations will be applied to b[e ] supergiants in sect .  [ s_b[e ] ] . in sect . [ s_bs ] we will investigate the effect of bi - stability on rotating stars and sect . [ s_concl ] concludes with a summary of this work and a discussion on the formation of disks of b[e ] supergiants . the computation of the dynamics of line driven stellar winds is a complicated problem in radiation hydrodynamics . it involves the simultaneous solution of the equations of motion , the rate equations and the radiative transfer equations in order to calculate the radiative acceleration .    for hot luminous stars the line forces are the most important driving forces in the wind . cak have shown that the line acceleration can be parameterized in terms of the optical depth parameter @xmath7 as @xmath8 , where @xmath0 and @xmath1 are parameters that depend on composition and temperature of the wind . in this expression @xmath0 is a measure of the number of lines and @xmath1 is a measure of the distribution of the line strengths with @xmath9 or 1 for a pure mix of optically thin or thick lines respectively .    in this paper we adopt the cak formalism and simplify the equations of motion by assuming a stationary , radial flow of a viscousless fluid . the possible influence of magnetic fields will be ignored . these simplifications are subject to the following restrictions ( see e.g. abbott 1980 ) :    * the cak theory treats the line acceleration in a simplified way in terms of the force multipliers , described above . detailed calculations of the line driving , including millions of lines , show that the total line acceleration can be well described by this simple representation ( pauldrach et al . both observations and calculations show that line driven stellar winds are not stationary . even a wind in a stationary solution will develop shocks ( see e.g. lucy  @xcite ) . however , the time - averaged structure follows the stationary state quite well ( see owocki et al . @xcite , and feldmeier 1999 ) . therefore we will restrict ourselves to stationary models . * a sectorial model will be adopted . this means that for every latitude a one dimensional problem will be solved . wind compression as in the wcd model will thus be neglected , but the main effect of wind compression is expected to _ redistribute _ the mass loss and not to change the _ total _ mass that is lost from the rotating star . * the absence of viscosity is a good approximation at the high density of line driven winds ( see cak ) . the representation of the line force by a simple power law may seem to be a gross simplification of the underlying physics of myriad line absorption processes , but it can be shown to hold for a homogeneously distributed mix of optically thick and thin lines ( abbott  @xcite , gayley  @xcite ) . calculations of cak and more recently in nlte by vink et al . ( @xcite ) of realistic model atmospheres , based on calculations of unified atmosphere / wind models with the isa - wind code ( de koter et al . 1997 ) , confirm this to be generally true to good accuracy for the important part of the wind ( from about the sonic point up to a few stellar radii ) . therefore , the complex physics of ion populations can be ignored in the investigation of the various effects of rotation on the stellar wind . however , there is one notable exception which is called the bi - stability jump . this refers to a jump in @xmath6 around the temperature of 21000 k where @xmath10 climbs from 1.3 ( lower temperature ) to 2.6 ( higher temperature ) . this jump was observed by lamers et al . ( @xcite ) in b supergiants . it is linked to a shift in ionization states of fe in the lower part of the wind . the line driving in the lower part of the wind is dominated by iron . below about @xmath11 k , fe iv recombines to fe iii and since fe iii is a more efficient line driver than fe iv , the wind structure changes dramatically ( vink et al . @xcite ) . the bi - stability jump may also occur in the temperature difference between the pole and the equator of a fast rotating star with gravity darkening . therefore two sets of force multiplier parameters for the wind will be adopted to reflect the sudden change in ionization states ( a high - temperature set for the pole and a low - temperature set for the equator ) . intuitively it is clear that rotation has an effect on the shape of a star and on the motion of the gas in the wind . an additional effect that the rotation of a star can have is the darkening of the equatorial regions of a star via the von zeipel effect . these effects will modify the wind of a star considerably and we will incorporate them in the line driven wind theory for rotating stars . the shape of a uniformly rotating star with all its mass concentrated in the core is determined by the equipotential surfaces of the potential in a rotating frame ( roche model ) : @xmath12    where @xmath13 , @xmath14 is the co - latitude ( @xmath14 = 0 at the pole ) and @xmath15 . @xmath16 is of course independent of the longitude @xmath17 . note that @xmath18 . so the maximum oblateness is @xmath19 . the critical ( break - up ) velocity is defined in this paper as : @xmath20 where @xmath21 is the effective mass of the star ,    @xmath22    which is the reduced mass due to radiation pressure by electron scattering ( see below ) .    the surface of a rotating star is implicitly given by eq . [ eq_firot ] . solving this expression for @xmath23 is equivalent to the solution of a 3rd degree polynomial , which results in    @xmath24    because of the difference in radius    @xmath25    and the difference in rotational velocity between pole and equator , the radial escape velocity at the equatorial region is smaller than at the pole    @xmath26    we can then roughly estimate the effect of rotation on @xmath6 as this value has an almost linear dependence on @xmath27 ( see later in  [ s_ssol ] ) . the von zeipel theorem  ( @xcite ) for distorted stars states that the radiative flux from a point on the star is proportional to the local effective gravity : @xmath28 with @xmath29 the boltzmann constant . as @xmath30 , we can write down the flux of a rotating star as a function of the co - latitude @xmath14 . this ( lengthy ) expression for the flux as a function of @xmath14 is plotted in fig . [ fg_f ] for a few different values of @xmath31 . the rotation of the star can induce a @xmath14 component of the flow of the matter in the wind . in this study a possible @xmath14 component is neglected as was discussed in  [ s_physics ] , but the motion in the longitudinal @xmath17 direction must still be taken into account .    in the absence of forces that can exert a torque on the wind , the equation of motion for @xmath32 is given by the conservation of angular momentum . note however that there could be a torque from the line forces themselves , which is dependent on the velocity gradient of the gas . we assume that there are no external torques , so conservation of angular momentum gives :    @xmath33    the equation of motion for the radial direction becomes :    @xmath34    where @xmath35 is the radial velocity and @xmath36 is the stellar mass . the conservation of mass for a non - spherical sectorial wind can be written as    @xmath37    where @xmath38 is the  local mass loss rate ,  i.e. the total mass loss rate if the solution for this latitude where valid for a spherical star . the total mass loss rate from the star is    @xmath39    the equation of motion together with the conservation of mass governs the dynamics of the wind given the equation of state , @xmath40 ( @xmath41 is the isothermal sound speed ) , and the radiative acceleration @xmath42 . using the conservation of mass , the pressure term can be rewritten as : @xmath43    as in the cak - theory the temperature structure is specified a priori as a function of @xmath44 . ( the results depend only very weakly on the chosen temperature structure ; see pauldrach et al . @xcite ) . the radiative acceleration consists of two components . the continuum component due to the electron scattering , @xmath45 , and a second component , @xmath46 , due to line scattering and absorption processes . the continuum acceleration is given in terms of @xmath47 : @xmath48 where @xmath49 is the electron opacity , @xmath50 is the stellar luminosity and @xmath51 is the radiation flux . the contributions by other continuum opacities are small . nevertheless , they were included in the calculation of the force multipliers by vink et al . ( 1999 ) used in our bi - stable wind models ( sect . for a homogeneous spherical star , @xmath47 is given by as the ratio between the continuum force and gravity and not as the ratio between continuum force and the critical radiation force for a rotating star ] :    @xmath52    in the more general case of a non - homogeneous , non - spherical star the continuum acceleration can be defined as a correction to the classical continuum acceleration of eq . [ eq_gam ] : @xmath53 where @xmath54 is the continuum correction factor which is given by :    @xmath55    the line acceleration @xmath56 has a more complicated form , since it is a function of distance @xmath44 as well as the velocity @xmath35 and velocity gradient @xmath57 . combining the equation of motion ( eq .  [ eq_bewvgl ] ) with the rewritten pressure and continuum acceleration terms and multiplying by @xmath58 , one finds :    @xmath59    this equation is valid for each co - latitude @xmath14 . following cak , the critical point of this equation is found by imposing the singularity condition ( where the subscript @xmath60 indicates values at the critical point ) : @xmath61 and the regularity condition : @xmath62 note that the critical point is not the sonic point @xmath63 ( @xmath64 ) due to the fact that there is an additional dependence on @xmath65 in the line acceleration @xmath46 ( abbott  1980 ; @xmath5 chapt 3.3 ) . the line forces are described within the framework of cak in the sobolev approximation . this means that the intristic line absorption profiles are considered to be infinitely sharp . then the line force becomes a function of local properties of the wind only . the acceleration due to an ensemble of lines is given by the summation of the contributions of all individual lines with rest frequencies @xmath66 . @xmath46 is given by ( castor  @xcite ) : @xmath67    where @xmath68 is the absorption coefficient per gram for the @xmath69th line , @xmath70 , @xmath71 the intensity at the rest frequency @xmath72 and @xmath73 is the sobolev optical depth of @xmath69th line , defined as @xmath74 . the integration is performed over the complete visible disk of the star . following cak , the line acceleration can be rewritten in terms of the force multiplier as a function of the optical depth parameter @xmath75 . where @xmath75 is defined as : @xmath76 with @xmath77 being the thermal velocity of the protons and @xmath78 is some reference value for the electron scattering @xmath79 @xmath80  g@xmath81 ( see _ isw _ , chapt . the line acceleration in terms of the continuum acceleration is give by :    @xmath82    if the star is assumed to be a point source , @xmath83 can be approximated by a simple power law parameterization ( cak , abbott 1982 ) : @xmath84 where the @xmath85 accounts for the effect of the electron density on the ionization balance in the wind . here @xmath86 is the electron density in units of @xmath87 @xmath88 and @xmath89 is the geometric dilution factor . friend & abbott ( 1986 ) and pauldrach et al . ( 1986 ) showed the importance of the finite disk correction on the line force in case of an extended source . this finite disk correction factor @xmath90 is given by ( cak ) : @xmath91    where    @xmath92    this expression is completely general : the actual shape of the star and its intensity distribution enter through the integration domain and the @xmath93 dependence . however , here we assume the intensity @xmath94 to be locally isotropic , i.e. we neglect limb darkening . we have tested the effect of limbdarkening . the inclusion of limbdarkening results in a small @xmath14-independent correction on @xmath90 . the line acceleration @xmath46 including all the neccesary correction factors is now given by :    @xmath95    combining eqs . [ eq_bewvglf ] and  [ eq_gl ] gives the full equation of motion for the sectorial wind of a rotating star    @xmath96 for line driven winds of non - rotating stars solutions of the equation of motion have been retrieved by e.g. cak and pauldrach et al . ( 1986 ) . these will be presented to serve as an illustration of the solution of the equation of motion and as a basis to interpret the more complicated results from the full equation of motion ( [ eq_bewvglrot ] ) .        if the star is considered to be a point soure the integral of the line acceleration in eq . [ eq_gmeerlijn ] collapses to one point . in this case @xmath90 disappears from eq . [ eq_gl ] and the line acceleration becomes a simple function @xmath97 . the solution of eq . [ eq_bewvglf ] was found by cak and can easily be explained by neglecting the gas pressure terms ( @xmath98 ) . in this case [ eq_bewvglf ] reduces to    @xmath99    where @xmath100 is a constant containing the mass loss rate . the equation is solved by imposing uniqueness of the solution ( kudritzki et al . 1989 ) which then fixes the value of the constant and thus @xmath101 . the solution for the mass loss and velocity law is given by    @xmath102    @xmath103    with    @xmath104    this simplified solution is equal to the full cak solution , in the limit of small sound speed , @xmath105 ( see also @xmath5 chapt 8) .      in case of a homogeneous spherical star the finite disk correction factor @xmath90 ( eq .  [ eq_corr ] ) can be calculated analytically ( cak ) :    @xmath106    including the finite disk correction results in an increase in @xmath6 , a decrease in @xmath107 and a modification of the simple scaling laws that were found in the original cak approach ( see friend & abbott  @xcite ; pauldrach et al . the decrease of the mass loss rate compared to the point source case is due to the decrease in @xmath46 close to the star where @xmath101 is determined . close to the stellar photosphere the finite disk correction is smaller than one , viz . @xmath108 . the accompanying increase in @xmath6 of typically a factor of two is due to two effects : ( 1 ) a reduction of @xmath107 results in a smaller amount of material to be accelerated and ( 2 ) far from the photosphere , the correction factor @xmath90 becomes larger than 1 . the resulting dependence between @xmath6 and @xmath27 is approximated by friend & abbott ( 1986 ) :    @xmath109    they also found a relation for the modified mass loss rate :    @xmath110      the solution of eq . [ eq_bewvglrot ] gives the mass loss rate , or rather the local mass loss rate @xmath38 , and velocity structure of the wind . the solution is complicated however by the presence of @xmath90 . this is an integral of a velocity - dependent function times the surface intensity over the visible section of the star . the appearance of the star varies throughout the wind in intensity distribution as well as in shape , as the star is no longer spherical and @xmath111 is a function of latitude . since the analytic solution would be cumbersome , a numerical approach of the solution is chosen . this is done with the analytic solutions of simple , i.e. the non - rotating , models in mind for comparison . the numerical solution of eq . [ eq_bewvglrot ] is relatively straightforward in case the function @xmath90 is a given function of @xmath44 . therefore we solve eq . [ eq_bewvglrot ] with @xmath112 in an iterative way as follows :    1 . a @xmath113 law for @xmath114 is assumed , viz . @xmath112 is calculated using this velocity law . note that @xmath90 is not dependent on the actual value of @xmath114 but on the velocity gradient . 3 .   eq . [ eq_bewvglrot ] is solved using this correction factor @xmath112 , obtaining a new velocity law @xmath114 . step 2 - 3 are repeated until convergence is reached . typically three iterations are sufficient , since further iterations changed the obtained values by less than 0.5 % . we checked that the solutions converged to the same value , independent of the starting value . the calculation of the correction factor @xmath90 from the velocity law ( step 1 . in the scheme described above ) is performed by numerical evaluation of eq . [ eq_corr ] . this is a non - trivial task , since both the integrand ( eq .  [ eq_zeipel ] ) and the shape of the visible `` disk '' ( eq . [ eq_rtheta ] ) are dependent on @xmath44 . the main parameters needed to obtain @xmath90 for a fixed co - latitude @xmath116 are : the stellar surface @xmath117 and the surface temperature @xmath118 , both given as a function of the co - latitude @xmath14 , as well as the velocity law @xmath114 . it is convenient to change the integration variables in eq . [ eq_corr ] from the wind centered coordinates @xmath14 and @xmath17 to the star centered coordinates @xmath119 and @xmath17 . for the definition of these coordinates , see fig . [ fg_ofr ] . note that here @xmath17 is not the stellar longitude anymore , but the angle of rotation around the line from a point @xmath120 in the wind to the stellar center . @xmath121    where @xmath94 is the frequency - integrated intensity    @xmath122    the angles are related to one another via    @xmath123    with the stellar surface given by @xmath124 , and @xmath125 equation  [ eq_coint ] can be integrated , considering that the limb angle @xmath126 can be determined numerically , yielding the finite disk correction for the line acceleration of an arbitrary shaped star .    the correction factor @xmath54 for the electron acceleration ( defined in eq . [ eq_dc ] ) similarly becomes :    @xmath127    this correction is also included in the solution of the equation of motion . [ tb_res1 ]    [ cols="^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^ " , ]     ( 1 ) : the adopted stellar parameters are : @xmath128 ; @xmath129 ; @xmath130 k.      to solve the equation of motion it is neccesary to determine the conditions at the critical point : @xmath131 , @xmath132 and @xmath133 , or , equivalently , @xmath131 , @xmath132 and @xmath134 , as well as the local mass loss rate @xmath135 . apart from eqs . [ eq_sing ] and  [ eq_regu ] and the equation of motion ( [ eq_bewvglrot ] ) itself the additional constraint for the continuum optical depth @xmath136 at the stellar radius is neccesary ( pauldrach et al . 1986 ) to uniquely determine these quantities :    @xmath137    if the condition for @xmath138 is used one may adopt any small value , say @xmath139 km / s , since even a large error in @xmath138 will only yield a relatively small error in @xmath6 and @xmath135 . this is due to the fact that at the surface of the star @xmath114 and thus also @xmath140 vary exponentially with scale height @xmath141 . this means that the radius where eq . [ eq_tau2/3 ] is fulfilled is very close to @xmath142 . this has been checked a posteriori . in our calculations we have adopted this boundary condition for @xmath138 in the solution of the full momentum equation .    a first guess for the local mass loss rate obtained from the cak solution is used to integrate the equation of motion from the stellar surface outward . this gives an approximate value for the critical point where @xmath143 if @xmath144 . this approximate value for the critical point can then be used to consecutively improve the guessed value of @xmath135 using eqs .  [ eq_sing ] and  [ eq_regu ] . this iterative procedure is terminated when the local mass loss rate converges . the calculation is then completed by integrating past the critical point on the solution which extends to infinity ( see cak or abbott 1980 for a description of the topology of the solutions ) . to deal with a numerical singularity at the sonic point the equations are solved for @xmath145 rather than for @xmath114 . we have calculated models for rotating b[e ] stars using the above described method . for the first models we adopted the following stellar parameters , which are typical for b[e ] stars : @xmath146 , @xmath147 , @xmath148 , @xmath149 and solar abundances . the effective temperature is defined as @xmath150 where @xmath151 is the total surface of the distorted star . the following force multiplier parameters from pauldrach et al . ( @xcite ) were adopted to describe the line acceleration : @xmath152 , @xmath153 , @xmath154 . an overview of the results for this generic b[e ] star rotating at @xmath31 = 0 , 0.3 , 0.4 , 0.5 and 0.6 times the critical rotation speed is given in table 1 . the calculation of models with higher values of @xmath31 runs into numerical problems , because in the iteration , the approximate critical point can be below the sonic point . the table shows for different co - latitudes @xmath14 and different rotation velocities @xmath31 the values of @xmath6 , the local mass loss rate @xmath38 and the value of the velocity law parameter @xmath113 , obtained by fitting the calculated velocity structure from @xmath155 to @xmath156 to a @xmath113-law . we see that the value of @xmath113 does not change much . [ fg_vr ] shows the velocity structure of the wind and the location of the critical point for @xmath31 = 0 , 0.3 and 0.6 . note that the critical point is very close the star . this is true for all finite disk corrected solutions . figure  [ fg_vmp ] shows the behaviour of @xmath6 versus @xmath135 for various rotational speeds and latitudes . as expected , @xmath6 at the equator is smaller than at the pole because @xmath6 scales roughly with @xmath27 and due to the smaller value of @xmath27 at the equator . figure [ fg_mp ] shows the latitude dependence of the local mass - loss rate and of the terminal velocity from equator to pole , for different rotational speeds , indicated by @xmath31 . at the pole ( @xmath157 ) @xmath6 is about constant for various rotational speeds the mass - loss rate at the equator decreases slightly with increasing rotation . the increase of the mass loss rate at the equator as calculated by pauldrach et al . ( @xcite ) and friend & abbott ( @xcite ) is offset by a decrease in mass loss rate due to the smaller radiative flux at the equator due to the von zeipel effect . the mass loss rate at the pole increases strongly with increasing rotation rate due to the increase in luminosity at the pole . this is because the total luminosity of the star is fixed and the smaller luminosity at the equator must therefore be compensated by a higher luminosity at the pole . the overall effect of rotation is to increase the mass loss rate of a star . this can be seen in fig . [ fg_mp ] where the overall mass - loss rate ( i.e. @xmath135 integrated over the stellar surface ) is plotted vs. the rotational velocity . the mass - loss rate of the models with @xmath158 vary with @xmath31 as @xmath159 . figure  [ fg_ldens ] shows the resulting density contrast ( far from the star ) between pole and equator ( @xmath160 ) . the densities at pole and equator are essentially the same : the smaller @xmath6 at the equator is offset by the larger mass loss rate @xmath135 at the pole . [ fg_dc ]    we have also calculated models for different luminosities ( see table 2 ) , and different stellar masses ( see table 3 ) . the mass - loss rate for the different luminosities is also plotted in fig .  [ fg_mp ] . the density contrast for various luminosities is plotted in fig . [ fg_ldens ] . we see that the density at the pole compared to the equator increases for increasing luminosity and the effect on the mass - loss rate is only visible for the highest luminosity . this is an effect of the continuum acceleration @xmath54 . 7 illustrates the effect of the continuum acceleration on the forces in the wind for a star with high luminosity ( @xmath162 ) . the net effect of the continuum radiation pressure and the gravity , i.e. @xmath163 , is an inward force in the equatorial region and an outward force in the polar region where the flux is higher . for small luminosity the variation in @xmath164 through @xmath54 is not important because @xmath47 itself is small . for high luminosity @xmath165 increases and thus the effect of @xmath166 on @xmath164 becomes noticeable ( this is with all other properties of the star fixed ) . in eq . [ eq_vinfsimple ] we see that a larger value of @xmath166 at the pole causes a relatively larger decrease ( through @xmath164 ) in @xmath27 and thus in @xmath6 , whereas the mass loss rate at the equator will increase less because of the smaller @xmath54 . this is only a qualitative explanation and the actual dependence of @xmath167 does not follow the relation suggested by eqs . 29 and 31 . we have found that the actual values of @xmath167 could be approximated quite well ( within 10 percent ) with the following relation :    @xmath168    where the @xmath169 is the radius where @xmath54 reaches its minimum value . so we see that although rotation alone modifies the wind structure considerably , it can not be responsible for the density contrast between the equatorial and polar wind observed in b[e ] stars , in case only radial effects are considered . we conclude that at least one other physical effect must be responsible for the formation of disks of rotating b[e ] stars . the two most likely additional effects are the bi - stability and the flow of the wind material towards the equator . in the section above we have shown that rotation alone can not explain the observed large density contrast between the poles and the equator of b[e ] stars . one of the possible additional mechanisms to enhance this contrast is the bi - stability jump .    to investigate whether the bi - stability jump observed in normal b - supergiants can explain the formation of disks around b[e ] stars , we have calculated bi - stable wind models for a fast rotating typical b[e ] supergiant . spherical , nlte wind models for normal b supergiants have been calculated by vink et al . ( 1999 ) . from these models cak - like force multiplier parameters have been derived using a monte - carlo method to simulate photon - gas interactions . these models show the occurence of a bi - stability jump around @xmath11 k.    for two models on either side of the bi - stability jump ( @xmath170 = 17 500 and 30 000 k ) , the cak parameters were determined by fitting a power law from about the sonic point to about @xmath171 ( from @xmath172 to @xmath173 in optical depth , @xmath75 is defined in eq . [ eq_t ] ) to the calculated force multiplier values . the resulting values for @xmath0 , @xmath1 and @xmath2 are listed in table  [ tb_ka ] . note that the @xmath2 parameter is taken equal to zero as it was not possible to extract it explicitly from the models . its effect is hidden in the constant @xmath0 . as expected , the values for the force multiplier parameters are quite different for the two cases , because the ionization has changed dramatically over the bi - stability jump . cccc & @xmath0 & @xmath1 & @xmath2 + 17500 k & 0.57 & 0.45 & 0.0 + 30000 k & 0.06 & 0.65 & 0.0 +    [ tb_ka ]    the predicted bi - stability jump occurs around @xmath174 kk . therefore we calculated a model for a rotating b[e ] star with the following properties : @xmath175 k , @xmath176 , @xmath177 and solar abundances . this implies that the pole of the rapidly rotating star will be hotter than 25 kk and the equator will be cooler than 25 kk .    for the pole , the force multipliers @xmath0 and @xmath1 from the hot monte carlo model of 30 000 k were used , whereas for the equator , @xmath0 and @xmath1 from the cool model of 17 500 k were used . the resulting @xmath178 and mass loss rates @xmath38 are plotted in figs . [ fg_vtheta ] and [ fg_mdtheta ] for @xmath31 = 0.6 . clearly visible is the drastic decrease in @xmath6 and the drastic increase of @xmath135 towards the equator . this occurs around a co - latitude @xmath14 ( angle from the pole ) of about 60 degrees , where @xmath179 k. ( the precise location of the jump depends on the effective temperature of the star ) . since both @xmath135 and @xmath6 have a equator / pole contrast of about a factor of three , the resulting density contrast is about a factor of 10 in a disklike region with a half opening angle of 30 degrees . our calculations show that the density contrast in the wind between pole and equator due to the bi - stability jump is significantly larger than without the bi - stability jump . yet , the calculated value of about a factor 10 is not sufficient to explain the observed density contrast of a factor @xmath180 . so another mechanism is needed to enhance the density contrast even further . this is most likely the wind compression mechanism ( bjorkman &   cassinelli 1993 ; bjorkman 1998 ) . we have modified the line driven wind theory for rotating stars by including the effects of oblateness and gravity darkening as well as rotational terms in the equation of motion . we considered a sectorial wind , i.e. we neglected the effect of the motion of the wind towards the equator or towards the pole . this assumption is justified close to the stellar surface below the critical point of the momentum equation . this implies that our method will predict about the correct distibution of mass - loss rates from the star as a function of stellar latitude , but it may not be accurate enough to predict the velocity and density distribution further away from the star , if motions in the @xmath14 direction become important . the equation of motion was solved using an iterative numerical scheme , that includes the conservation of angular momentum and the correction factors to the radiative acceleration by lines and by the continuum , due to the non - spherical shape of the star and due to the latitude dependence of the radiative flux . this method was applied to study its possible effects on the formation of disks around fast rotating b[e ] supergiants . the models with constant force multiplier parameters @xmath0 , @xmath1 and @xmath2 show a decrease of _ both _ the mass loss rate and the terminal velocity from pole to equator . this is mainly due to two effects : the reduction of the escape velocity from pole to equator , resulting in a higher terminal velocity at the pole , and the reduction of the radiative flux from the pole to the equator due to gravity darkening , which results in a decrease in mass loss rate at the equator .    for a star with a fixed luminosity and a fixed polar radius , the temperature at the pole increases and the temperature at the equator decreases with increasing rotation rate . the terminal velocity of the wind from the poles is almost independent of @xmath181 but at the equator @xmath6 decreases with increasing rotation rate . the mass - loss rate at the pole increases with increasing @xmath181 , due to the increase in @xmath170 at the poles , but the mass loss rate from the equator is almost independent of the rotation . the combination of these effects alone produce a density contrast between the polar and the equatorial wind of a factor @xmath182 , except for fast rotating stars with high luminosity ( in which case @xmath183 ) . our results agree quantitatively with those obtained by maeder ( @xcite ) for the latitudinal dependence of @xmath135 . our predictions of the global mass loss rates are also in general agreement with those of maeder for models with luminosities far from the eddington limit . for high luminosities we find a strong polar outflow , whereas maeder s models predict an enhanced equatorial outflow on the basis of the original cak scaling laws . this difference is due to the fact that we have included an improved description of the radiation pressure in a rotating star . the difference between the polar wind and the equatorial wind is strongly enhanced when the bi - stability of radiation driven winds between @xmath149 and 30  000 k is taken into account . in this case the force multipliers @xmath0 and @xmath1 of a rapidly rotating star change drastically with stellar latitude if the pole is hotter than 25  000 k and the equator is cooler . in   [ s_bs ] we have applied the newly calculated force multipliers from vink et al . ( 1999 ) above and below the bi - stability jump to the models of rotating b[e ] supergiants . one might argue that it is not allowed to apply the force multipliers of ( bi - stable ) spherical wind models to aspherical winds of rotating stars . however , we have shown that a disk formed through bi - stability will be sufficiently thick ( typically about 30 degrees for a star with an effective temperature equal to about 25  000 ) to make these spherical models reasonable approximations for the conditions in the wind close to the star , where the mass loss rate is determined . we find that rotationally induced bi - stability models of b[e ] stars reach a density contrast of about a factor 10 between the dense equatorial wind and the less dense polar wind . this is less than the factor @xmath184 that is observed ( see e.g. bjorkman  @xcite ) . the extra density increase is most likely due to the wind compression . the gas that leaves the photosphere from a rotating wind , will follow an orbit in a tilted plane defined by the local rotation vector and the center of the star . if the rotational velocity is large or the wind velocity is small , this orbit will cross the equatorial plane where the streamlines of the wind from different stellar latitudes cross . the resulting shock will compress the gas into a thin outflowing disk with an opening angle of only a few degrees ( bjorkman & cassinelli 1993 ; @xmath5 chapt 11 ) . we have not considered this wind compression in our model . however , it is likely to occur in rotationally induced bi - stable winds , because the wind velocity in the equatorial plane is about a factor three smaller than from the poles , thus facillitating the wind compression . owocki & cranmer ( 1994 ) have argued that the flow of the wind towards the equatorial plane , predicted in the wcd - theory , may be offset by a @xmath14-component of the line acceleration towards the polar regions . this `` disk inhibition '' mechanism operates in wind models with constant force multipliers @xmath0 and @xmath1 . however , it is not clear that this mechanism is sufficiently strong to overcome the wind compression in the rib - model . this is because @xmath0 and @xmath1 change with latitude in the rib - model and there is a strong density gradient in the wind from the equator to the poles ( owocki et al . 1998 , puls et al . 1999 ) . the combination of the rib and the wcd mechanism may offer the best possibility for explaining the disks of b[e]-supergiants . whether the compression is strong enough to explain the observed high density contrast between the polar and the equatorial wind , remains to be seen . it can be calculated by combining the solutions of the wind momentum equation of rotating oblate stars with gravity darkening ( derived in this paper ) , and the calculation of the resulting trajectories of the wind . the combination of the rotation induced bi - stability model and the wind compressed disk model is promising for explaining the disks of b[e ] stars because the rib - mechanism explains the increased mass loss and the small velocity from the equatorial regions and the wcd - mechanism explains the strong compression of the disk . abbott d.c . , 1980 , apj 242 , 1183 abbott d.c . , 1982 , apj 259 , 282 bjorkman j.e . , 1998 , in _ b[e ] stars _ , eds a.m. hubert , c. jaschek kluwer : dordrecht , assl 233 , p 189 bjorkman j.e . , cassinelli j.p . , 1993 , apj 409 , 429 cassinelli 1998 in _ b[e ] stars _ , eds a.m. hubert , c. jaschek kluwer : dordrecht , assl 233 , p 177 castor j.i . , 1974 , mnras 169 , 279 castor j.i . , abbott d.c . , klein r.i . , 1975 , apj 195 , 157 cranmer s.r . , owocki s.p . , 1994 , apj 440 , 308 feldmeier a. , 1999 in _ variable and non - spherical stellar winds in luminous hot stars _ , lecture notes in physics , eds . b. wolf et al . springer : heidelberg , p 285 . friend d.b . , abbott d.c . , 1986 , apj 311 , 701 gayley k.g . , 1995 , apj 454 , 410 kudritzki r.p . , pauldrach a. , puls j. , abbott d.c . , 1989 , a&a 219 , 205 lamers h.j.g.l.m . , cassinelli j.p . , 1999 , introduction to stellar winds , cambridge univ . press . ( @xmath5 ) lamers h.j.g.l.m . , pauldrach a. , 1991 , a&a 244 , l5 lamers h.j.g.l.m . , , lindholm d.m . , 1995 , apj 455 , 269 lamers h.j.g.l.m . , zickgraf f .- j . , de winter d. , houziaux l. , zorec j. , 1998 , a&a 340 , 117 lucy l.b . , 1982 , apj 255 , 286 maeder a. , 1999 , a&a 347 , 185 owocki s.p . , castor j.i . , rybicki g.b . , 1988 , apj 335 , 914 owocki s.p . , cranmer r.c . , 1994 , apj 424 , 887 owocki s.p . , cranmer r.c . , gayley k.g . , 1998 , ap&ss 260 , 149 pauldrach a. , puls j. , kudritzki r.p . , 1986 , a&a 164 , 86 puls j. , petrenz p. , owocki s.p . , 1999 in _ variable and non - spherical stellar winds in luminous hot stars _ , lecture notes in physics , eds . b. wolf et al . springer : heidelberg , p 131 . , de koter a. , lamers h.j.g.l.m . , 1999 , a&a 350 , 181 von zeipel h. , 1924 , mnras 84 , 665 zickgraf f .- j . , 1992 , in _ astronomical ccd observing and reduction techniques _ howell s.b . , asp conference series vol . 22 , 75 zickgraf f .- j . , schulte - ladbeck r.e . , 1989 , a&a 214 , 274 zickgraf f .- j . , wolf b. , stahl o. , leitherer c. , klare g. , 1985 , a&a 143 , 421
Abstract:

Answer: we have formulated the momentum equation for sectorial line driven winds from rotating stars including : ( a ) the oblateness of the star , ( b ) gravity darkening ( von zeipel effect ) , ( c ) conservation of angular momentum , ( d ) line driving specified by the force multiplier parameters ( @xmath0 , @xmath1 , @xmath2 ) , ( e ) finite disk correction factors for an oblate star with gravity darkening for both the continuum and the line driving . the equations are solved numerically . we calculated the distribution of the mass flux and the wind velocity from the pole to the equator for the winds of b[e]-supergiants . rotation decreases the terminal velocity in the equatorial region but hardly affects the wind velocity from the poles ; it enhances the mass flux from the poles while the mass flux from the equator remains nearly the same . these effects increase with increasing rotation rates . we also calculated models with a bi - stability jump around 25  000 k , using force multipliers recently calculated with a monte carlo technique . in this case the mass flux increases drastically from the pole to the equator and the terminal velocity decreases drastically from pole to equator . this produces a density contrast in the wind @xmath3 of about a factor 10 independent of the rotation rate of the star . we suggest that the observed density contrast of a factor @xmath4 of the disks of b[e ] stars may be reached by taking into account the wind compression due to the trajectories of the gas above the critical point towards the equatorial plane .


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: we consider the problem of reconstruction of refractive indices or dielectric constants of unknown targets placed in a homogeneous domain from blind backscattered experimental data . we work with time resolved backscattering experimental data of wave propagation for a 3-d hyperbolic coefficient inverse problem ( cip ) . our data are generated by a single location of the point source . the backscattering signal is measured on a part of a plane . we present a combination of the approximately globally convergent method of @xcite with a finite element method ( fem ) for the numerical solution of this cip . given a certain function computed by the technique of @xcite , the fem reconstructs the unknown coefficient in an explicit form . as a result , we can reconstruct refractive indices and locations of targets . in addition , we estimate their sizes . we believe that these results can be used as initial guesses for locally convergent methods in order to obtain better shapes , see , e.g. section 5.9 in @xcite , where the image obtained by the globally convergent method for transmitted experimental data was refined via a locally convergent adaptivity technique . experimental data were collected by the device which was recently built at university of north carolina at charlotte . in our experiments we image targets standing in the air . a potential application of our experiments is in imaging of explosives . note that explosives can be located in the air @xcite , e.g. improvised explosive devices ( ieds ) . the work on real data for the case when targets are hidden in a soil is ongoing . we have collected backscattering time resolved experimental data of electrical waves propagation in a non - attenuating medium . as it was pointed out in @xcite , the main difficulty of working with such data is caused by a huge mismatch between these data and ones produced by computational simulations . conventional data denoising techniques do not help in this case . therefore , it is unlikely that any numerical method would successfully invert the raw data . to get the data , which would look somewhat similar with ones obtained in computational simulations , a heuristic data pre - processing procedure should be applied . the pre - processed data are used as the input for the globally convergent method . previously our research group has applied the method of @xcite to the simpler case of transmitted experimental data which were produced by a similar device ( chapter 5 of @xcite ) . the backscattering real data are much harder to work with than transmitted ones since the backscattered signal is significantly weaker than the transmitted one , as well as because some unwanted signals are mixed up with the true one , see figure [ fig2]-a ) for the latter . we refer to our research in @xcite and section 6.9 of @xcite for the case of backscattering real data in 1-d . in the current paper we present results of reconstruction of the 3-d version of the method of @xcite . the approximately globally convergent method of @xcite relies on the structure of the underlying pde operator and does not use optimization techniques . each iterative step consists of solutions of two problems : the dirichlet boundary value problem for an elliptic pde and the cauchy problem for the underlying hyperbolic pde . approximate global convergence " ( global convergence in short ) means that we use a certain reasonable approximate mathematical model . approximation is used because of one inevitably faces with substantial challenges when trying to develop globally convergent numerical methods for multidimensional cips for hyperbolic pdes with single source . it is rigorously established in the framework of this model that the method of @xcite results in obtaining some points in a small neighborhood of the exact coefficient without a priori knowledge of any point in this neighborhood , see theorem 2.9.4 in @xcite and theorem 5.1 in  @xcite . the distance between those points and the exact solution depends on the error in the data , the step size @xmath1 of a certain discretization of the pseudo - frequency interval and the computational domain @xmath2 where the inverse problem is solved ( see section [ sec:4.3 ] for definition of @xmath1 ) . a knowledge of the background medium in @xmath3 is also not required by this method . because of these theorems , convergence analysis is not presented here . a rigorous definition of the approximate global convergence property can be found in section 1.1.2 of @xcite and in @xcite . we use a mild approximation , since it amounts only to the truncation of a certain asymptotic series , and it is used only on the first iterative step ( section [ sec:4.2 ] ) . the validity of this approximate model was verified computationally on both synthetic and transmitted experimental data in @xcite as well as in the current work in the case of experimental backscattering data . different imaging methods are used to compute geometrical information of targets , such as their shapes , sizes and locations , see , e.g.  @xcite . on the other hand , refractive indices , which is our main interest , characterize constituent materials of targets , and they are much more difficult to compute . as to the gradient - like methods , we refer to , e.g. @xcite and references therein . convergence of these methods is guaranteed only if the starting point of iterations is chosen to be sufficiently close to the correct solution . on the other hand , it was shown in section 5.8.4 of @xcite that the gradient method failed to work for transmitted experimental data of @xcite in the case when its starting point was the background medium . an outline of this paper is as follows . in section 2 we state forward and inverse problems . in section 3 we describe the experimental data and briefly outline the data pre - processing procedure . in section 4 we briefly outline the method of @xcite : for reader s convenience . in section 5 we describe a version of the fem which works for our case . in section 6 we describe our algorithm . in section 7 we outline some details of our numerical implementation . results are presented in section 8 and summary is in section 9 . let @xmath4 be a convex bounded domain with the boundary @xmath5 denote by @xmath6 we model the electromagnetic wave propagation in an isotropic and non - magnetic space @xmath7 with the dimensionless coefficient @xmath8 which describes the spatially distributed dielectric constant of the medium . we consider the following cauchy problem for the hyperbolic equation @xmath9@xmath10we assume that the coefficient @xmath11 of equation ( [ 2.1 ] ) is such that @xmath12,~~\varepsilon _ { r}(\mathbf{x})=1\text { for } \mathbf{x}\in \mathbb{r}^{3}\diagdown \omega ,   \label{2.20}\]]where @xmath13 we _ a priori _ assume knowledge of the constant @xmath14 which amounts to the knowledge of the set of admissible coefficients in ( 2.20 ) . however , we do not assume that the number @xmath15 is small , i.e. we do not impose smallness assumptions on the unknown coefficient @xmath16 . below @xmath17 are hlder spaces , where @xmath18 is an integer and @xmath19 let @xmath20 be a part of the boundary @xmath21 later we will designate @xmath22 as the backscattering side of @xmath2 and will explain how we deal with the absence of the data at @xmath23    * coefficient inverse problem ( cip ) . * _ suppose that the coefficient _ @xmath24 _ _  satisfies . determine the function _ _ @xmath24 _ _  for _ _ @xmath25 _ _ , assuming that the following function _ _ @xmath26 _ _  is known for a single source position _ _ @xmath27 @xmath28    the function @xmath26 in ( [ 2.5 ] ) models time dependent measurements of the wave field at the part @xmath22 of the boundary @xmath29 of the domain of interest @xmath2 . we assume below that the source position is fixed and @xmath30 . this assumption allows us to simplify the resulting integral - differential equation because @xmath31 in @xmath32 . the assumption @xmath33 for @xmath34 means that the coefficient @xmath11 has a known constant value outside of the domain of interest @xmath35    this is a cip with single measurement data . uniqueness theorem for such cips in the multidimensional case are currently known only if the function @xmath36 in ( [ 2.2 ] ) is replaced with a function @xmath37 such that @xmath38 @xmath39 a proper example of such function @xmath40 is a narrow gaussian centered around @xmath41 , which approximates the function @xmath36 in the distribution sense .  from the physics standpoint this gaussian is equivalent to @xmath42 that uniqueness theorem can be proved by the method , which was originated in @xcite . this method is based on carleman estimates , also see , e.g. sections 1.10 , 1.11 of the book @xcite about this method . the authors believe that , because of applications , it still makes sense to develop numerical methods for this cip without completely addressing the uniqueness question . the function @xmath43 in ( [ 2.1 ] ) represents the voltage of one component of the electric field @xmath44 in our computer simulations the incident field has only one non - zero component @xmath45 . this component propagates along the @xmath46axis until it reaches the target , where it is scattered . so , we assume that in our experiment @xmath47 we now comment on five main discrepancies between our mathematical model ( [ 2.1])- ( [ 2.20 ] ) and the reality . the first discrepancy which causes the main difficulties , is the aforementioned huge mismatch between experimental data and computational simulations . the second one is that , although we realize that equation ( [ 2.1 ] ) can be derived from maxwell equations only in the 2-d case , we use it to model the full 3-d case . the reason is that our current receiver can measure only one of the polarization components of the scattered electric field @xmath48 . in addition , if using a more complicated mathematical model than the one of ( [ 2.1 ] ) , for example the one that includes vector scattering and thus depolarization effects on scattering , then one would need to develop a globally convergent inverse method for this case . the latter is a quite time consuming task with yet unknown outcome . * *  * * equation ( [ 2.1 ] ) was used in chapter 5 of @xcite for the case of transmitted experimental data , and accurate solutions were obtained . a partial explanation of the latter can be found in @xcite , where the maxwell s system in a non - magnetic and non - conductive medium was solved numerically in time domain . it was shown numerically in section 7.2.2 of @xcite that the component of the vector @xmath49 which was initially incident upon the medium , dominates two other components . this is true for at least a rather simple medium such as ours . therefore , the function @xmath43 in ( [ 2.1 ] ) represents the voltage of the computed component @xmath50 of the electric field , which is emitted and measured by our antennas . the third discrepancy is that the condition @xmath51 is violated on the inclusion / background interface in our experiments . the fourth discrepancy is that formally equation ( [ 2.1 ] ) is invalid for the case when metallic targets are present . on the other hand , it was demonstrated computationally in kbksnf that one can treat metallic targets as dielectrics with large dielectric constants , which we call _ appearing dielectric constant _ , @xmath52modeling metallic targets as integral parts of the unknown coefficient @xmath24 is convenient for the above application to imaging of explosives . indeed , ieds usually consist of mixtures of some dielectrics with a number of metallic parts . such targets are heterogeneous ones , and we consider three heterogeneous cases in section [ sec:8.2 ] . on the other hand , modeling metallic parts of heterogeneous targets as a separate matter than the rest of an a priori unknown background medium would result in significant additional complications of the already difficult problem with yet unknown outcome . the fifth discrepancy is that we use the incident plane wave instead of the point source in our computations . we have discovered that the plane wave case works better in image reconstructions than the point source , while the point source case is more convenient for the convergence analysis in bk , bkjiip12 . in addition , since the distance between our measurement plane and targets is much larger than the wavelength of our signal , it is reasonable to approximate the incident wave as a plane wave . thus , our results of section [ sec:8.2 ] demonstrate the well known fact that computational results are often less pessimistic than the theory , since the theory can not grasp all nuances of the reality . in summary , we believe that accurate solutions of the above cip for experimental data justify our mathematical model . figure [ fig1]-a ) is a photograph explaining the data collection . the data collection is done in a regular room , which contains office furniture , computers , etc . keeping in mind our desired application ( see introduction ) , we intentionally did not arrange a special waveguide , which would protect our data from unwanted signals caused by reflections from various objects in the room . below @xmath53 and @xmath54 are horizontal and vertical axis respectively and the @xmath55 axis is perpendicular to the measurement plane , the positive direction of @xmath55 axis is in the direction from the target to the measurement plane . we dimensionalize our coordinates as @xmath56 where m " stands for meter . however , we do not change notations of coordinates for brevity . hence , below , e.g. 0.05 of length actually means 5 centimeters . `` the transmitter sends the pulse into the medium which contains targets of interest . the electric wave caused by the pulse is scattered by the targets , and the backscattered signal is detected by the detector . the detected signal is recorded by the real time oscilloscope . '' two main pieces of our device are picosecond pulse generator ( figure [ fig1]-b ) ) and textronix oscilloscope ( figure [ fig1]-c ) ) . the picosecond pulse generator generates electric pulses . the duration of each pulse is 300 picoseconds . this pulse goes to the transmitter , which is a horn antenna ( source ) . the transmitter sends the pulse into the medium which contains targets of interest . the electric wave caused by the pulse is scattered by the targets , and the backscattered signal is detected by the detector . the detected signal is recorded by the real time oscilloscope . the oscilloscope produces a digitized time resolved signal with the step size in time of 10 picoseconds . the total time of measurements for one pulse is 10 nanoseconds=10@xmath57 picoseconds=10@xmath58 second .    to decrease the measurement noise , the pulse is generated 800 times for each position of the detector , the backscattering signal is also measured 800 times and resulting signals are averaged . the detector moves in both horizontal and vertical directions covering the square @xmath59 on the measurement plane . we have chosen the step size of this movement to be 0.02 . although we can choose any step size , we found that 0.02 provides a good compromise between the precision of measurements and the total time spent on data collection . the distance between our targets and the measurement plane is approximately 0.8 with about 0.05 deviations , and the wavelength of our signal is about 0.03 . therefore , the distance between the measurement plane and our targets is of about 26 wavelengths . this is in the far field zone . the main difficulty working with experimental data is that there is a huge mismatch between these data and computationally simulated ones . indeed , figure [ fig2]-a ) depicts a sample of experimentally measured data for a wooden block at one position of the detector , see figure [ fig1]-a ) for data collection scheme . on this figure , the direct signal is the signal going directly to the receiver . we use this direct signal as the time reference for data pre - processing . unwanted signals are due to reflections of the electric wave from several objects present in the room . figure [ fig2]-b ) presents the computationally simulated data for the same target , see section [ sec:7.1 ] for data simulations . these figures show a huge mismatch between real and computationally simulated data . therefore , data pre - processing is necessary . we refer to @xcite for details of our data pre - processing procedure . the main steps of this procedure include :    1 . _ time - zero correction_. the time - zero correction is to shift the measured data in time . so that its starting time is the same as when the incident pulse is emitted from the transmitter this is done using the direct signals from the transmitter to the detector as the time reference . 2 .   _ extraction of scattered signals_. apart from the backscattered wave by the targets , our measured data also contain various types of signals , e.g.  direct signals from the horn to the detector , scattered signals from structures inside the room , etc . what we need , however , is the scattered signals by the targets only . to obtain them , we single out the scattered signals caused by the targets only and remove all unwanted signals . data propagation_. after getting the scattered signals , the next step of data pre - processing is to propagate the data closer to the targets , i.e.  to approximate the scattered wave on a plane which is much closer to the targets then the measurement plane . the distance between that propagated plane and the front surface of a target is usually between 0.02 and 0.06 ( compare with the 0.8 distance from the measurement plane ) . there are two reasons for doing this . the first one is that the method of @xcite works with the laplace transform of the function @xmath43 ( section [ sec:4 ] ) . that laplace transform decays exponentially in terms of the time delay , which is proportional to the distance from the target to the measurement plane . hence , the amplitude of the laplace transformed experimental data on the measurement plane is very small and can be dominated by computational round - off error . the second reason is that this propagation procedure helps to substantially reduce the computational cost since the computational domain for the inverse problem is reduced . 4 .   _ data calibration . _ finally , since the amplitudes of the experimental incident and scattered waves are usually significantly different from simulations , we need to bring the former to the same level of the amplitude as the latter . this is done using a known target referred to as _ calibrating object_.    in this paper , the result of data pre - processing is used as the measured data @xmath60 on the backscattering boundary @xmath22 of our computational domain @xmath2 for the inverse problem . in this section we briefly outline the globally convergent method for reader s convenience . we refer to sections 2.3 , 2.5 , 2.6.1 and 2.9.2 of bk as well as @xcite for details . the first step of our inverse algorithm is the laplace transform of the function @xmath61@xmath62where @xmath63 is a certain number . we assume that the number @xmath63 is sufficiently large , and we call the parameter @xmath64 _ pseudo frequency_. it follows from ( [ 2.1 ] ) , ( [ 2.2 ] ) and ( [ 2.6 ] ) that the function @xmath65 is the solution of the following problem @xmath66@xmath67the limit ( [ 2.8 ] ) is proved in theorem 2.7.1 of @xcite . in addition , it was proven in theorem 2.7.2 of @xcite that for the function @xmath24 satisfying ( [ 2.20 ] ) there exists unique solution @xmath68 of the problem ( 2.7 ) , ( [ 2.8 ] ) for every @xmath69 such that @xmath70where @xmath71 is the solution of the problem ( 2.7 ) , ( [ 2.8 ] ) for the case @xmath72@xmath73      it follows from theorem 2.7.2 of @xcite that @xmath74 hence , we can consider the functions @xmath75 , @xmath76 @xmath77substituting @xmath78 in ( [ 2.7 ] ) and keeping in mind that the source @xmath30 , we obtain @xmath79using ( [ transform ] ) we obtain @xmath80where the truncation pseudo frequency @xmath81 is a large number , which is chosen numerically , see section [ sec:8 ] for details . we call @xmath82 the _ tail function , _ and it is unknown . it follows from ( [ transform ] ) and ( 3.7 ) that @xmath83it follows from @xcite ( section 2.3 ) that , under some conditions , there exists a function @xmath84 such that the following asymptotic behavior with respect to @xmath85 holds for functions @xmath86 and @xmath87 @xmath88@xmath89 differentiating both sides of equation ( [ new3.2 ] ) with respect to @xmath64 then using ( [ transform ] ) and ( [ 3.7 ] ) , we obtain the following nonlinear integral differential equation @xmath90 . \end{split } \label{3.8}\]]in addition , ( [ 2.5 ] ) and ( [ transform ] ) lead to the following dirichlet boundary condition for the function @xmath87 @xmath91 ,   \label{3.9}\]]@xmath92here @xmath93 is the laplace transform ( 2.6 ) of the function @xmath60 in ( [ 2.5 ] ) . we now need to complement the boundary data ( [ 3.9 ] ) at the backscattering side @xmath22 with the boundary data at the rest of the boundary @xmath94 using computationally simulated data , it was shown numerically in section 6.8.5 of @xcite as well as in @xcite that it is reasonable to approximate the boundary data on @xmath95 by the solution of the forward problem for the homogeneous medium for the case @xmath96 : recall that this equality holds outside of the domain @xmath97 see ( [ 2.20 ] ) . thus , we use below the following dirichlet boundary condition for the function @xmath98@xmath99 ,   \label{3.91}\]]@xmath100 , \\   \psi ^{0}\left ( \mathbf{x},s\right ) , \forall \left ( \mathbf{x},s\right ) \in \left ( \partial \omega \setminus \gamma \right ) \times \left [ \underline{s},\overline{s}\right ] .\end{array}\right . \label{3.92}\]]where the function @xmath101 is the function @xmath102 in ( [ 3.90 ] ) computed for the case @xmath103 .    even though equation ( [ 3.8 ] ) with the boundary condition ( [ 3.91 ] ) has two unknown functions @xmath87 and @xmath86 , we can approximate both of them because approximation procedures for them are different , see section 7.1 . suppose for a moment that functions @xmath87 and @xmath86 are approximated in @xmath2 together with their derivatives @xmath104 then the corresponding approximation for the coefficient @xmath11 can be found via backwards calculation using ( [ new3.2 ] ) .      to start iterations , we need the first approximation @xmath105 for the tail function . in this section we show how to calculate @xmath106 this is the same choice as the one in section 2.9.2 of the book @xcite as well as in @xcite . let the function @xmath107 satisfying ( 2.20 ) be the exact solution of our cip for the exact data @xmath108 in ( [ 2.5 ] ) . let @xmath109 be the exact tail function  defined as @xmath110let @xmath111 $ ] be the corresponding exact function @xmath98 satisfying equation ( [ 3.8 ] ) . let @xmath112 $ ] be the corresponding exact dirichlet boundary condition for @xmath113 as defined in ( 3.91 ) .  following ( [ 3.92 ] ) , we assume that @xmath114 for @xmath115 .$]hence , ( [ 3.8 ] ) and ( [ 3.91 ] ) hold for functions @xmath116 setting in ( 3.8 ) @xmath117 we obtain @xmath118next , truncating the second term in each of the asymptotics ( [ 2.b ] ) and ( [ 2c ] ) , we obtain that there exists a function @xmath119 such that @xmath120substituting formulae ( [ asymptotic ] ) into ( [ 4.18 * * ] ) , we obtain the following approximate dirichlet boundary value problem for the function @xmath121 @xmath122@xmath123thus , using ( [ 4.17 * ] ) and ( [ asymptotic ] ) , we obtain the following approximate mathematical model .    * approximate mathematical model . *    _ we assume that there exists a function _ @xmath124 _ _  such that the exact tail function _ _ @xmath125 _ _  has the form _ _ @xmath126__and the function _ _ @xmath127 _ _  is__@xmath128    because of ( [ 2.917 ] ) , ( [ 2.918 ] ) and ( [ 2.914 ] ) , we set for the first tail @xmath129where the function @xmath130 is the solution of the following dirichlet boundary value problem @xmath131@xmath132we point out that we calculate @xmath105 without any advanced knowledge of a small neighborhood of the exact coefficient @xmath133 using ( [ asymptotic])-([2.925 ] ) and schauder theorem @xcite , we obtain @xmath134where the number @xmath135 depends only from the domain @xmath35 hence , the error in the calculation of @xmath105 depends only on the error in the boundary data @xmath136 on the other hand , since the boundary function @xmath137 is generated by the function @xmath26 in ( [ 2.5 ] ) , then the error in @xmath138 is generated by the error in measurements . the estimate ( [ 2.926 ] ) is one of elements of the proof of the approximate global convergence theorem for this numerical method , see theorem 2.9.4 in @xcite and theorem 5.1 in @xcite . although a good approximation for the exact solution @xmath107 can be derived from the function @xmath139 we have observed computationally that better approximations are delivered via iterations described below in sections [ sec:6.1 ] , [ sec:6.2 ] .      to approximate both functions @xmath87 and @xmath86 using ( [ 3.8 ] ) and ( [ 3.91 ] ) , we consider the layer stripping procedure with respect to @xmath64 . we divide the interval @xmath140 $ ] into @xmath141 small subintervals with the uniform step size @xmath142 . here , @xmath143 we approximate the function @xmath144 as a piecewise constant function with respect to @xmath145 i.e. we assume that @xmath146 for @xmath147 hence , using ( [ 3.7 ] ) , we approximate the function @xmath148 as @xmath149to obtain a sequence of dirichlet boundary value problems for elliptic pdes for functions @xmath150 we introduce the @xmath151dependent carleman weight function ( cwf ) @xmath152 , $ ] where @xmath153 is a large parameter . in our numerical studies we take @xmath154 . this function mitigates the influence of the nonlinear term in the resulting integral - differential equations on every pseudo - frequency interval @xmath155 . multiply both sides of equation ( [ 3.8 ] ) by @xmath156 and integrate with respect to @xmath157 we obtain @xmath158here @xmath159 is such an approximation of the tail function @xmath160 which corresponds to the function @xmath161 ( section [ sec:6.1 ] ) . numbers @xmath162 are computed explicitly . furthermore , @xmath163 for this reason we ignore the nonlinear term in ( [ 4.21 ] ) , thus setting@xmath164note that ( [ 4.22 ] ) is not a linearization , since ( [ 4.21 ] ) contains products @xmath165 and also because the tail function @xmath166 depends nonlinearly on functions @xmath167 see ( [ 2.a ] ) and step 6 in section [ sec:6.1 ] . in this section we explain how we compute functions @xmath169 on every pseudo - frequency interval @xmath170 using the fem . once the functions @xmath171 along with the function @xmath172 in ( [ 4.21 ] ) are calculated , we compute the function @xmath173 using the direct analog of ( [ 4.20 ] ) , @xmath174using ( [ transform ] ) , we set @xmath175 .   \label{3.104}\ ] ] to find the function @xmath176 we note that the function @xmath177 is the solution of the following analog of the problem ( [ 2.7 ] ) , ( [ 2.8])@xmath178@xmath179where @xmath180 \text { for } \mathbf{x}\in \partial \omega . \label{3.107}\]]to compute the function @xmath181 from ( [ 3.105 ] ) , ( [ 3.106 ] ) and ( [ 3.107 ] ) , we apply a version of the fem as described below in sections [ sec:5.1 ] , [ sec:5.2 ] .      following @xcite we discretize in computations our bounded domain @xmath4 by an unstructured tetrahedral mesh @xmath182 using non - overlapping tetrahedral elements @xmath183 . the elements @xmath184 are such that @xmath185 , where @xmath186 is the total number of elements in @xmath2 , and @xmath187we associate with the mesh @xmath182 the mesh function @xmath188 as a piecewise - constant function such that @xmath189where @xmath190 is the diameter of @xmath184 which we define as the longest side of @xmath184 . we impose the following shape regularity assumption of the mesh @xmath182 for every element @xmath191 @xmath192where @xmath193 is the radius of the maximal sphere contained in the element @xmath184 . define the set of polynomials @xmath194 as @xmath195we introduce now the finite element space @xmath196 as @xmath197where @xmath198 denotes the set of linear functions on @xmath184 defined by ( polynom ) for @xmath199 . hence , the finite element space @xmath196 consists of continuous piecewise linear functions in @xmath2 . to approximate functions @xmath200 we introduce the space of piecewise constant functions @xmath201 , @xmath202where @xmath203 is the piecewise constant function on @xmath184 defined by ( polynom ) for @xmath204 .      to compute the function @xmath205 from ( [ 3.105 ] ) , we formulate the finite element method for the problem ( [ 3.106])-([3.107 ] ) as : find the function @xmath206 for the known function @xmath207 such that @xmath208where @xmath209 is the scalar product in @xmath210 . we expand @xmath211 in terms of the standard continuous piecewise linear functions @xmath212 in the space @xmath196 as @xmath213where @xmath214 denote the nodal values of the function @xmath211 at the nodes @xmath215 of the elements @xmath184 in the mesh @xmath182 . we can determine @xmath214 by knowing already computed functions @xmath216 using the following analog of ( [ 3.104 ] ) @xmath217 , \forall \mathbf{x}\in \omega .\]]substitute ( [ 3.200 ] ) into ( [ 3.107_1 ] ) and choose @xmath218 then we obtain the following linear algebraic system of equations @xmath219 , \label{3.108}\]]where @xmath220 $ ] is the scalar product in @xmath221 the system ( [ 3.108 ] ) can be rewritten in the matrix form for the unknown vector @xmath222 and known vector @xmath223 as @xmath224here @xmath225 is the block mass matrix in space , @xmath226 is the stiffness matrix corresponding to the term containing @xmath227 in ( [ 3.108 ] ) and @xmath228 is the load vector . at the element @xmath184 the matrix entries in ( [ 3.108_1 ] ) are explicitly given by : @xmath229    to obtain an explicit scheme for the computation of coefficients @xmath205 , we approximate the matrix @xmath225 by the lumped mass matrix @xmath230 in space , i.e. , the diagonal approximation is obtained by taking the row sum of @xmath225 @xcite . we obtain @xmath231note that for the case of linear lagrange elements which are used in our computations in section [ sec:8 ] we have @xmath232 . thus , the lumping procedure does not include approximation errors in this case . we present now our algorithm for the numerical solution of equations ( 4.21 ) and computing the functions @xmath205 using the equation([3.109 ] ) . in this algorithm the index @xmath233 denotes the number of inner iterations inside every pseudo - frequency interval @xmath234 when we update tails . step 0 : :    set @xmath235 . compute the initial tail function    @xmath236    as in ( 2.923)-([2.925 ] ) . step 1 : :    here we describe iterations which update tails inside every    pseudo - frequency interval @xmath155 . let @xmath237 suppose that functions @xmath238    are computed . solve the dirichlet boundary value problem for the    function    @xmath239    @xmath240 step 2 : :    compute functions @xmath241 and    @xmath242    @xmath243@xmath244 .\ ] ] step 3 : :    compute the function    @xmath245 via backwards    calculations , using the finite element formulation of equation    ( [ 3.109 ] ) as    @xmath246since    by ( [ 2.20 ] ) we should have    @xmath247    and also since we need to extend the function    @xmath248    outside of the domain @xmath2 by unity , we set    @xmath249 step 4 : :    solve the forward problem ( [ 2.1])-([2.2 ] ) with    @xmath250    and compute the laplace transform ( [ 2.6 ] ) for @xmath251 . we    obtain the function    @xmath252 . step 5 : :    update the tail function    as@xmath253continue    inner iterations with respect to @xmath233 until the stopping    criterion of step 1 of section [ sec:6.2 ] is met at    @xmath254 . step 6 : :    set for the pseudo - frequency interval @xmath255    @xmath256 step 7 : :    if either the stopping criterion with respect to @xmath257 of    step 4 of section [ sec:6.2 ] is met , or @xmath258 then set the    resulting function @xmath259 as    the solution of our cip . otherwise , set @xmath260 and go to    step 1 . when testing the algorithm of section [ sec:6.1 ] on experimental data , we have developed a reliable stopping criterion for iterations @xmath261 in this algorithm . on every pseudo - frequency interval @xmath155 we define first norms  @xmath262 as @xmath263 in ( [ first_norms ] ) the function @xmath264 is the computed tail functions at the inner iteration @xmath265 as in ( [ 7.104 ] ) . functions @xmath266 in ( [ first_norms ] ) are obtained from the known measured function @xmath26 in ( [ 2.5 ] ) as @xmath267where @xmath268 is the laplace transform of the function @xmath26 at @xmath269 .    we have observed that computed first norms  @xmath262 always achieve only one minimum at a certain @xmath270 , where the number @xmath271 depends on the specific set of experimental data . furthermore , in non - blind cases of non - metallic targets , the corresponding values of @xmath272 were in a good agreement with a priori known ones . however , in the cases of non - blind metallic targets we have observed that @xmath273 this contradicts with ( [ 2.51 ] ) . therefore , we have developed the following stopping criterion which consists of four steps .    * the stopping criterion *    the first step in our criterion is for stopping inner iterations with respect to @xmath233 in step 5 of section [ sec:6.1 ] . as to steps 2 - 4 , they are for stopping outer iterations with respect to @xmath257 ( step 7 in section [ sec:6.1 ] ) . first , we define numbers @xmath274 and @xmath275 as@xmath276@xmath277 in ( [ 7.6_1 ] ) functions @xmath278 are computed tail functions corresponding to @xmath279 ( step 6 in section [ sec:6.1 ] ) and functions @xmath280 are calculated using ( [ obstail ] ) . * * step 1 . * iterate with respect to @xmath233 and stop iterations at @xmath281 such that@xmath282or @xmath283 where @xmath284 is a chosen tolerance . * * for every @xmath257 * *  * * compute final norms  @xmath285 as @xmath286 in ( [ final_norms ] ) functions @xmath287 are computed as in ( [ 7.104 ] ) . * * step 3*. compute the number @xmath288 of the pseudo frequency interval such that the first norms @xmath262 in ( [ first_norms ] ) achieve its first minimum with respect to @xmath257 and get corresponding @xmath289 on this interval . compute the number @xmath290 of the pseudo frequency interval such that the final norms @xmath291 in ( [ final_norms ] ) achieve its first minimum or they are stabilized with respect to @xmath257 , and get corresponding @xmath292 on this interval . next , compute the number @xmath293 @xmath294 * * step 4*. if @xmath295 or @xmath296 then take the final reconstructed value of the refractive index @xmath297 as the computed function @xmath298 take @xmath299 and stop iterations . however , if @xmath300 then continue iterations and compute the number @xmath301 $ ] of the pseudo frequency interval such that the global minimum with respect to @xmath257 of final norms @xmath285 in ( [ final_norms ] ) is achieved . then , similarly with ( [ 7.7 ] ) , compute the number @xmath302 @xmath303and take @xmath304 as the final reconstructed value of the refractive index . also , take the function @xmath305 as the computed coefficient @xmath24 and stop iterations . we have observed in all our computations that conditions of our stopping criterion are always achieved . more precisely , one of conditions ( [ 7.5 ] ) , ( [ 7.5_1 ] ) is always achieved for iterations with respect to @xmath233 and the minimal values mentioned in steps 3 and 4 are always achieved . figure [ fig : finalnorms ] displays a typical @xmath306dependence of sequences @xmath307and @xmath308 in this section we present some additional details of our numerical implementation . because of ( [ 2.51 ] ) , we define in all our tests the upper value of the function @xmath309 as @xmath310 see ( [ 2.20 ] ) . thus , we set lower and upper bounds for the reconstructed function @xmath311 in @xmath3 as @xmath312 \}. \label{4.30}\]]as to the lower bound , we ensure it via ( [ 7.101 ] ) . we ensure the upper bound 15 similarly via truncating to 15 those values of @xmath313 which exceed this number . to solve dirichlet boundary value problems ( [ 7.102 ] ) , we use fem . we reconstruct refractive indices rather than dielectric constants of material since they can be directly measured .    to compare our computational results with directly measured refractive indices @xmath314 of dielectric targets and with appearing dielectric constants of metallic targets ( see ( [ 2.51 ] ) ) , we consider maximal values of computed functions @xmath315 , @xmath316 see step 4 of section [ sec:6.2 ] for the definition of @xmath317 using experimental data for non - blind targets and comparing reconstruction results with cases of synthetic data , we have found that our algorithm provided accurate results with the following pseudo frequency interval , which we use in all our computations @xmath318,\underline{s}=8,\overline{s}=10,h=0.05.\ ] ]      as it is clear from step 4 of section [ sec:6.1 ] , we need to solve the forward problem ( [ 2.1 ] ) , ( [ 2.2 ] ) on each iterative step of inner iterations to update the tail via ( [ 7.103 ] ) . since it is impossible to computationally solve equation ( [ 2.1 ] ) in the infinite space @xmath319 we work with a truncated domain . namely , we choose the domain @xmath226 as @xmath320 we use the hybrid fem / fdm method described in @xcite and the software package waves @xcite . we split @xmath226 into two subdomains @xmath321 and @xmath322 so that @xmath323 . we solve the forward problem in @xmath226 and the inverse problem via the algorithm of section [ sec:6.1 ] in @xmath324 the space mesh in @xmath325 and in @xmath322 consists of tetrahedral and cubes , respectively . below @xmath326 since by ( [ 2.20 ] ) @xmath33 in @xmath327 then it is computationally efficient to use fdm in @xmath322 and to use fem in @xmath328 as it is done in the hybrid method of @xcite . the front and back sides of the rectangular prism @xmath226 are @xmath329 and @xmath330 , respectively . the boundary of the domain @xmath226 is @xmath331 here , @xmath332 and @xmath333 are , respectively , front and back sides of the domain @xmath226 , and @xmath334 is the union of left , right , top and bottom sides of this domain . the front side @xmath22 of the rectangular prism @xmath2 where the propagated data @xmath60 in ( [ 2.5 ] ) are given , is @xmath335 now we describe the forward problem which is used in our computations . to compute tail functions @xmath336 via steps 4 , 5 of the algorithm of section [ sec:6.1 ] , we computationally solve the following forward problem in our tests : @xmath337 , \\ { n}u & = -\partial _ { t}u,~\mbox{on}~\partial _ { 1}g\times ( t_{1},t ) , \\ \partial _ { n}u & = -\partial _ { t}u,~\mbox{on}~\partial _ { 2}g\times ( 0,t ) , \\ \partial _ { n}u & = 0,~\mbox{on}~\partial _ { 3}g\times ( 0,t ) , \end{split } \label{8.2}\]]where @xmath338 is the amplitude of the initialized plane wave , @xmath339we use @xmath340 and @xmath341 we solve the problem ( [ 8.2 ] ) using the explicit scheme with the time step size @xmath342 which satisfies the cfl condition . our reconstruction procedure is done in two stages described in this section .      in the first stage we follow the algorithm of section [ sec:6.1 ] . we have observed that this stage provides accurate locations of targets of interest . it also provides accurate values of refractive indices @xmath343 of dielectric targets and large values of appearing dielectric constants @xmath344 for metallic targets , see ( [ 7.7 ] ) and ( [ 7.8 ] ) . however , the algorithm of section [ sec:6.1 ] does not reconstruct well sizes / shapes of targets . thus , we need a postprocessing procedure , which is done in the second stage . let @xmath345 be the function in ( 7.101 ) . then we set @xmath346next , we determine minimal @xmath347 and maximal @xmath348 values in @xmath53 and @xmath54 directions , where the function @xmath349 next , we set@xmath350 , y\in \left [ y_{\min } , y_{\max } \right ] , \\   1 & \text { otherwise}\end{cases}\]]and proceed with step 5 of the algorithm of section [ sec:6.1 ] . in this second stage we perform the same number of iterations with respect to both indices @xmath351 as ones of the first stage . we are concerned in the second stage only with sized and shapes of targets , and we are not concerned with values of @xmath352 rather , we take these values from the first stage . let @xmath353 be the function @xmath24 obtained at the last iteration of the second stage . then we form the image of the target based on the function @xmath354@xmath355    .__object names . _ _ [ cols="^,<,<",options="header " , ] we collected experimental backscattering time resolved data of electrical wave propagation and have applied the approximately globally convergent numerical method of @xcite to these data . results for four non - blind and seven blind cases show a good accuracy of reconstruction of refractive indices of dielectric targets and appearing dielectric constants of metallic targets . in the case of dielectrics , the average reconstruction error is at least three times less than the error of direct measurements . we confidently differentiate between metallic and dielectric targets . in particular , we have accurately computed maximal values of refractive indices / dielectric constants of three blind heterogeneous targets . these targets represent simplified models of improvised explosive devices ieds , which are heterogeneous ones . locations of targets and their sizes in @xmath356 directions are accurately reconstructed . the most difficult cases of sizes in the @xmath46direction ( depth ) are well reconstructed in some cases . in addition , shapes of some targets are well reconstructed in some cases . we believe that a follow up application of the locally convergent adaptivity technique might improve reconstructions of shapes of targets . the adaptivity takes the solution obtained by the approximately globally convergent method as the starting point for the minimization of the tikhonov functional on a sequence of adaptively refined meshes . a significant refinement via the adaptivity was demonstrated in section 5.9 of @xcite for the case of transmitted experimental data , see figures 5.13 and 5.16 in @xcite . this research was supported by us army research laboratory and us army research office grants w911nf-11 - 1 - 0325 and w911nf-11 - 1 - 0399 , the swedish research council , the swedish foundation for strategic research ( ssf ) through the gothenburg mathematical modelling centre ( gmmc ) and by the swedish institute , visby program . the authors are grateful to mr . steven kitchin for his excellent work on data collection . l. beilina , k. samuelsson and k. ahlander , efficiency of a hybrid method for the wave equation . in _ international conference on finite element methods _ , gakuto international series mathematical sciences and applications , gakkotosho co. , ltd , 2001 . l. beilina and m.v . klibanov , a new approximate mathematical model for global convergence for a coefficient inverse problem with backscattering data , _ j. inverse and ill - posed problems _ , 20 , 513 - 565 , 2012 . l.  beilina , energy estimates and numerical verification of the stabilized domain decomposition finite element / finite difference approach for the maxwell s system in time domain , _ central european journal of mathematics _ , 11 , 702 - 733 , 2013 . kuzhuget , l. beilina , m.v . klibanov , a. sullivan ,  l. nguyen and m.a . fiddy , blind experimental data collected in the field and an approximately globally convergent inverse algorithm , _ inverse problems _ , 28 , 095007 , 2012 . n.  t. thnh , l.  beilina , m.  v. klibanov , and m.  a. fiddy . reconstruction of the refractive index from experimental backscattering data using a globally convergent inverse method , preprint , arxiv:1306.3150 [ math - ph ] , 2013 .
Abstract:

Answer: we consider the problem of reconstruction of dielectrics from blind backscattered experimental data . experimental data were collected by a device , which was built at university of north carolina at charlotte . this device sends electrical pulses into the medium and collects the time resolved backscattered data on a part of a plane . the spatially distributed dielectric constant @xmath0 is the unknown coefficient of a wave - like pde . this coefficient is reconstructed from those data in blind cases . to do this , a globally convergent numerical method is used .    * keywords * : coefficient inverse problem ( cip ) , finite element method , globally convergent numerical method for cip , experimental backscattered data . * ams classification codes : * 65n15 , 65n30 , 35j25 .