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: a carousel is an automated storage and retrieval system , widely used in modern warehouses . it consists of a number of shelves or drawers , which are linked together and are rotating in a closed loop . it is operated by a picker ( human or robotic ) that has a fixed position in front of the carousel . a typical vertical carousel is given in figure  [ fig : carousel ] . carousels are widely used for storage and retrieval of small and medium - sized items , such as health and beauty products , repair parts of boilers for space heating , parts of vacuum cleaners and sewing machines , books , shoes and many other goods . in e - commerce companies use carousel to store small items and manage small individual orders . an order is defined as a set of items that must be picked together ( for instance , for a single customer ) . carousels are highly versatile , and come in a huge variety of configurations , sizes , and types . they can be horizontal or vertical and rotate in either one or both directions . although both unidirectional ( one - way rotating ) or bidirectional ( two - way rotating ) carousels are encountered in practice , the bidirectional types are the most common ( as well as being the most efficient )  @xcite . one of the main advantages of carousels is that , rather than having the picker travel to an item ( as is the case in a warehouse where items are stored on shelves ) , the carousel rotates the items to the picker . while the carousel is travelling , the picker has the time to perform other tasks , such as pack or label the retrieved items , or serve another carousel . this practice enhances the operational efficiency of the warehouse . carousel models have received much attention in the literature and continue to pose interesting problems . there is a rich literature on carousels that dates back to 1980  @xcite . in section  [ s : literature ] we shall review some of the main research topics that have been of interest to the research community so far . to name a few , one may wish to study various ways of storing the items on a carousel ( storage arrangements ) so as to minimise the total time needed until an order is completed ( response time ) or the strategy that should be followed in rotating the carousel so as the total time the carousel travels between items of one order is minimised ( travel time for a single order ) . one may also consider design issues , for instance , the problem of pre - positioning the carousel in anticipation of storage or retrieval requests ( choosing a _ dwell point _ ) in order to improve the average response time of the system . the list of references presented here is by no means exhaustive ; it rather serves the purpose of indicating the continuing interest in carousels .    in this review paper we focus on the modelling and the performance of carousel systems . usually a carousel is modelled as a circle , either as a discrete model @xcite , where the circle consists of a fixed number of locations , or as a continuous one @xcite , where the circle has unit length and the locations of the required items are represented as arbitrary points on the circle . throughout this paper we shall view the carousel as a continuous loop of unit length . beyond this initial assumption , we shall examine modelling issues such as how to model travel times or picking times of items in a system of several carousels so as to be able to derive approximations of various performance characteristics . under `` performance '' one may understand a variety of notions . for example , in single - carousel single - order problems ( cf.section  [ s : nelly ] ) , the performance measure under consideration is the travel time of the carousel until all items in an order are picked . on the other hand , in section  [ s : maria ] , performance may be measured by the time the picker is idle between picking items from various carousels , i.e.by the picker s utilisation .    in this paper we consider two research topics in detail . in section  [ s : nelly ] , we discuss the problem of choosing a reasonable picking strategy for one order and a single carousel , where the order is represented as a list of items , and by order pick strategy we mean an algorithm that prescribes in which sequence the items are to be retrieved . we present a general probabilistic approach developed by litvak _ et al_.  @xcite to analytically derive the probability distribution of the travel time in case when items locations are independent and uniformly distributed . this line of research seems to be the only example in the literature where exact statistical characteristics of the travel time have been obtained by means of a systematic mathematical approach . the presented technique is based on properties of uniform spacings and their relations to exponential distributions . we demonstrate the effectiveness of this method by considering several relevant order - picking strategies , such as the greedy nearest - item strategies and so - called @xmath2-step strategies that provide a good approximation for the optimal ( shortest ) route . in section  [ s : maria ] we consider the second topic that relates to multiple - carousel settings and the modelling challenges that appear in such problems . having optimised the travel time of a single carousel for a single order , one wonders if optimising locally every time each order on each carousel leads to the best solution ( fastest , cheapest , or with the largest picker utilisation ) for a complicated system . as is mentioned later on , multiple - carousel problems become too complicated too quickly , and often exact analysis is not possible . therefore , we discuss which concessions have to be made in order to be able to obtain estimates of the performance measures we are interested in , and we give in detail the impact that these concessions have on our estimations . there exist a few exact results for two - carousel models and related models in healthcare logistics ; see boxma and vlasiou @xcite and vlasiou _ et al . _  @xcite@xcite . however , to the best of our knowledge , no exact results exist for systems involving more than two carousels . preferably , these two research topics that we consider in this paper should be studied in parallel . however , establishing any exact results , say on determining the optimal retrieval and travelling strategy for a multiple - carousel model , without any restrictions to the sequence the items in an order are picked or the sequence the carousels are served , seems to be intractable . nonetheless , quite a few research opportunities related to the optimal design and control of carousel systems are still available . we elaborate on further research topics in section  [ s : research ] . we conclude with section  [ s : literature ] , which outlines the problems examined so far on carousels and related storage and retrieval systems . performance analysis of single units is a necessary step in structural design of order pick systems  @xcite . in a setting of a single order on a single carousel , the major performance characteristic is the response time , that is , the total time it takes to retrieve an order . the response time consists of pick times needed to collect the items from their locations by an operator , and the travel ( rotation ) time of the carousel . while pick times can hardly be improved , the travel time depends on the location of each item and the order picking sequence , and thus , it is subject to analysis and optimisation . therefore , in this section , we discuss properties of the travel time needed to collect an order of @xmath3 items . in this section , our focus is on the case when the item locations are randomly distributed on a carousel circumference . this model allows one to compute statistical characteristics of the travel time such as the average travel time or the travel time distribution . later on , in section  [ ss : picking a single order ] we discuss some results from the literature on evaluating the travel times under different assumptions on the items locations , in particular , the case when the pick positions are fixed . we note that in case of a single carousel , it is natural to assume that the pick times and the travel time are independent . the situation , however , is quite different in the systems of two or more carousels , where pick times on one carousel affect the travel times on other carousels . this issue will be discussed in detail in section  [ s : maria ] . the model addressed in this section is as follows . we model a carousel as a circle of length  1 . the order is represented by the list of @xmath3 items whose positions are independent and uniformly distributed on @xmath4 . for ease of presentation , we act as if the picker travels to the pick positions instead of the other way around . also , we assume that the acceleration / deceleration time of the carousel is negligible or that it is assigned to the pick time , and that the carousel rotates at unit speed . therefore the travel distance can be identified with the travel time ( see also section  [ ss : design issues ] ) . obviously , the travel time depends heavily on the pick strategy . here by order pick strategy _ we mean an algorithm that prescribes the sequence in which the items are collected . for example , assume that the picker always proceeds in the clockwise ( cw ) direction and denote by @xmath5 the time needed to collect @xmath3 items under this simple strategy . then , clearly , the distribution function @xmath6 of @xmath5 simply equals @xmath7 , @xmath8 . however , we would like to study strategies that provide smaller travel times . in this sense , a better algorithm that one can think of is the ` greedy ' strategy , also called the _ nearest - item heuristic _ : always travel to the nearest item to be picked ( as in figure  [ fig_ni_route ] ) . , width=144 ]    the nearest - item strategy indeed performs very well and is often used in practice , but the question is : `` what is the distribution of the travel time under the nearest - item heuristic ? '' . this problem is not at all trivial . for example , straightforward methods , such as conditioning on possible item locations , do not lead to feasible calculations . the same applies to the optimal strategy . bartholdi and platzman  @xcite showed that the shortest route admits at most one turn . intuitively , this follows merely by observing figure  [ fig_ni_route ] , where the displayed route can be shortened by collecting the first item in the counterclockwise direction and then collecting the rest of the items rotating clockwise . thus , the shortest route is merely the minimum among the @xmath9 candidate routes than have at most one turn . however , in spite of this simple structure of the shortest route , its distribution function is hard to derive .    below we discuss in detail a general methodology developed by litvak _ et al . _ @xcite to obtain the distribution of the travel time under various order pick strategies . the proposed technique is based on properties of uniform spacings and their connection with exponential random variables . we show how this approach allows us to derive exact and often counterintuitive results on several relevant order pick strategies . some other methods from the literature are described in section  [ ss : picking a single order ] . we start with introducing the notation and presenting some background results . let the random variable @xmath10 be the picker s starting point and the random variable @xmath11 , where @xmath12 , be the position of the @xmath13th item . we suppose that the @xmath11 s , @xmath14 , are independent and uniformly distributed on @xmath4 . let @xmath15 denote the order statistics of @xmath16 and set @xmath17 , @xmath18 . then the uniform spacings are defined as @xmath19 if we consider @xmath3 items randomly located on a circle , then the spacings @xmath20 are the distances between two neighbouring items , and the spacings @xmath21 and @xmath22 are the distances between the starting point and the two items adjacent to it . whatever strategy the picker takes , he always has to cover one or more uniform spacings on his way from one location to another . hence , in general , the travel time can be expressed as a function of the uniform spacings . uniform spacings have been analysed extensively in two classical review papers by pyke  @xcite . the author gives four useful constructions that establish a connection between uniform spacings and exponential random variables . we will use such a connection in the following form . let @xmath23 , @xmath24 be independent exponential random variables with mean  1 . moreover , define the random variables @xmath25 then , according to pyke  @xcite , uniform spacings can be represented as follows : @xmath26 here and throughout this paper @xmath27 means that @xmath28 and @xmath29 have the same probability distribution . linear combinations of uniform spacing have nice properties . in particular , the moments of linear combinations with non - negative coefficients can be easily computed , and their distribution function has been derived by ali  @xcite , ali and obaidullah  @xcite . now , let @xmath30 and @xmath31 be independent exponential random variables with parameters @xmath32 and @xmath33 , respectively . we write @xmath34 , @xmath35 , where @xmath23 and @xmath36 are independent exponential random variables with parameter  1 . then , given the event @xmath37 $ ] , we obtain the following useful statements : ( i )  the distribution of @xmath38 is exponential with parameter @xmath39 ( property of the minimum of two exponentials ) , which is distributed as @xmath40 ; ( ii )  since @xmath41 $ ] , then , according to the memoryless property , @xmath31 can be written as a sum of two terms : @xmath42 and another independent exponential with parameter @xmath33 , so @xmath31 is distributed as @xmath43 . ( iii )  it is easy to check that the distribution of @xmath44 is independent of the event @xmath37 $ ] because according to ( i ) and ( ii ) , given @xmath37 $ ] , @xmath45 is again distributed as @xmath46 ( see also chapter  2 of @xcite ) . based on the above - mentioned properties of exponential random variables , and their connections to uniform spacings and travel times , one may adopt the following methodology for analysing the travel times under various strategies  @xcite :    1 . express the travel time under a given strategy as a function of uniform spacings . 2 .   by conditioning on linear inequalities between the spacings and employing the above mentioned properties of exponential random variables , rewrite the travel time as a linear combination of uniform spacings or as a probabilistic mixture of such linear combinations . 3 .   use the results from  @xcite to obtain the moments and the distribution of the travel time . below we show how this approach works in case of the nearest - item heuristic  @xcite and so - called @xmath2-step strategies  @xcite .      under the nearest - item heuristic , the picker always moves towards the nearest item to be retrieved . the positions of the items partition the circle in @xmath47 uniform spacings @xmath48 defined by . under the nearest - item heuristic , the picker first considers the two spacings adjacent to his starting position and then travels to the nearest item . next he also looks at the other spacing adjacent to that item and compares the distance to the item located at the endpoint of that spacing and the distance to the first item in the other direction , which is the sum of the spacings previously considered . then he travels again to the nearest item , and so on . furthermore , by employing  , we may act as if the picker faces non - normalised exponential spacings , and afterwards divide the travel time ( which is equal to the travel distance ) by the sum of all spacings . then it is clear that each new spacing faced by the picker is independent of the ones already observed . now let @xmath49 , @xmath50 , denote the @xmath13-th non - normalised exponential spacing faced by the picker . that is , the spacings are numbered as observed by the picker operating under the nearest - item heuristic ( see figure  [ fig_travel ] ) . , width=340 ]    then @xmath51 can be expressed as @xmath52    we first provide an informal explanation of how the proposed methodology can be applied to . to start with , note that first term in the right - hand side of is @xmath53 , which is distributed simply as @xmath54 . moreover , under the event @xmath55 $ ] the rest of the sum remains unaltered . further , consider the term @xmath56 let @xmath57 , @xmath58 , @xmath59 be auxiliary independent exponential random variables with mean  1 . given @xmath60 $ ] , the random variable @xmath61 is distributed as @xmath62 , @xmath23 is distributed as @xmath63 and @xmath64 is distributed as @xmath59 . then the term in is distributed as @xmath65 . furthermore , given the event @xmath66 $ ] , we obtain that @xmath23 is distributed as @xmath62 , @xmath61 is distributed as @xmath67 and @xmath64 is distributed as @xmath68 . substituting the above in , we obtain again @xmath65 ! remarkably , under the event @xmath69 $ ] , again transforms into @xmath65 . furthermore , the sum @xmath70 becomes simply @xmath71 . we may now rename @xmath72 back to @xmath73 since the two 3-dimensional vectors are identically distributed . then the term becomes @xmath74 , and the rest of the terms in the right - hand side of remain unaltered in all three cases . proceeding further , we obtain the next statement which is proved rigorously in  @xcite . [ litvak and adan  @xcite ] [ th_nearest_item ] for all @xmath75 , @xmath76 and @xmath77 where @xmath78 if @xmath79 and @xmath80 otherwise .    here follows directly from and the result by ali  @xcite , which we applied in the form given by theorem  2 in @xcite . the above theorem is surprising because it provides an elegant solution for a problem that looks intractable at first . an interesting by - product is the distribution of the number of turns under the nearest - item heuristics and the counterintuitive result that _ the travel time and the number of turns are independent _  @xcite ! the latter can be seen directly from . indeed , a turn after step @xmath13 is equivalent to the event @xmath81 $ ] . however , as we saw earlier , the form of the distribution of the travel time is given by and it is independent of this sort of events .      under the @xmath2-step strategies , the picker chooses the shortest route among the @xmath82 routes that change direction at most once , and only do so after collecting no more than @xmath2 items . note that the optimal strategy is in fact an @xmath83-step strategy since it is never optimal to turn more than once , and the maximal possible number of items collected before a turn is @xmath84 . the @xmath2-step strategies give a good approximation for the shortest travel time . in fact , they often provide the optimal route even for moderate values of @xmath2 , as in figure  [ fig : m - step ] . et al . _ @xcite were the first to propose these strategies as an upper bound for the optimal route . in case of independent uniformly distributed pick positions , they obtained the distribution of the travel time under the @xmath2-step strategy for @xmath85 using analytical methods . later on , litvak and adan  @xcite applied the described methodology based on the properties of uniform spacings to completely analyse the travel time under the @xmath2-step strategies , provided @xmath86 . the travel under the @xmath2-step strategy can be expressed as follows @xmath87 indeed , the term @xmath88 is the gain in travel time ( compared to one full rotation ) obtained by skipping the spacing @xmath89 and going back instead , ending in a clockwise direction . on the other hand , @xmath90 is the gain obtained by skipping the spacing @xmath91 and going back ending counterclockwise . under the @xmath2-step strategy the picker skips the spacing that provides the largest possible gain ( see figure  [ fig : m - step ] ) . -step strategy.[fig : m - step],width=340 ]    using property , and after appropriate manipulations of exponential random variables , one can prove the following result . [ th_m - step ] for any @xmath92 , with @xmath86 , @xmath93    the maximum in the right - hand side of implies that @xmath94 is distributed as a complicated probabilistic mixture of linear combinations of uniform spacings  @xcite . the number of terms in this mixture is the well - known catalan number @xmath95 which grows extremely fast with @xmath2 . computing the expectations , we conclude that on average , the @xmath2-step strategy performs better than the nearest - item heuristic already for @xmath96 provided @xmath97 .    again , as a by - product , we can obtain the distribution of the number of steps before the turn . moreover , the latter random variable turns out to be independent of the travel time . this surprising statement follows from a similar reasoning as the independence of the travel time and the number of turns under the nearest - item heuristic . furthermore , when @xmath3 goes to infinity , the number of steps before the turn converges to a shifted geometric distribution with parameter @xmath98 . that is , in the limit , with probability @xmath98 there will be no turn , with probability @xmath99 there will be one step before a turn , etc . also , in the limit , the @xmath2-step strategy with @xmath86 coincides with the optimal strategy since the probability of achieving the minimal travel time by making more than @xmath100 steps before a turn will converge to zero . thus , for large enough @xmath3 , the probability that a 2-step strategy provides an optimal route is about @xmath101 . this explains the remarkably good performance of the @xmath2-step strategies .    as a side remark , we would like to note that @xcite provides slightly more general results than those presented in and .      since the optimal strategy simply coincides with the @xmath83-step strategy ( at most one turn after collecting at most @xmath84 items ) it can be analysed by methods from section  [ ss : m - step ] . however , the condition @xmath86 is violated for @xmath102 , and hence , does not hold . in fact , the proposed methodology applied to the optimal travel time @xmath103 very soon results in analytically infeasible calculations . litvak and van zwet  @xcite analysed the optimal route . they employed the results on the @xmath2-step strategy to derive a recursive expression for the distribution of the minimal travel time . we would like to also note that the process of comparing the lengths of the spacings and deriving corresponding linear combinations of normalised exponentials can be easily translated into a computer program . then , for moderate values of @xmath3 the exact distribution of the optimal travel time can be obtained numerically . the result will be a complicated mixture of linear combinations of uniform spacings . for large values of @xmath3 such exact calculations will require too much computer capacity . however , in this case , the knowledge of the exact distribution is not very important since one can apply approximations based on asymptotic results discussed in the next section . when the order is large , we can model this situation by letting @xmath104 . then the expressions in and for the travel time allow us to obtain asymptotic results that are of independent mathematical interest . obviously , if @xmath104 then the travel time under any strategy goes to one with probability  1 . however , with linear scaling , we obtain non - trivial distributions that we present below for the nearest - item heuristic and for the optimal travel time . [ theorem2 ] let @xmath23 , @xmath64 , @xmath105,@xmath57 , @xmath58 , @xmath105 , be independent exponentials with mean 1 . then @xmath106 as @xmath104 . result is also generalised to the case when items positions are independent and have some positive density @xmath107  @xcite . the expression in the right - hand side of is a well - known functional of the poisson process , which has been extensively studied in the literature . we will briefly discuss this topic in section  [ ss : exponential_functionals ] . the problems examined so far relate to one - carousel models . in industry though , one rarely meets a facility where only one carousel is used . multiple - carousel systems tend to have a higher level of throughput ; however , they increase the investment cost due to the extra driving and control mechanisms @xcite . a natural question is how much the throughput of a standard carousel can be improved by the corresponding multiple - carousel system that has the same number of shelves as the standard carousel . thus , the question we would like to examine in this section is the following : given a setup , i.e.a specific storage scheme of the items stored on the carousel _ and _ a specific travelling strategy , such as those described in the previous section , how much can we increase the utilisation of the picker ( by assigning to him more carousels to handle ) without increasing the response time of an order above some chosen level ? in other words , how do we reach a _ quality and efficiency _ regime in a real situation ?    to illustrate things better , consider the following simple example . a facility assigns @xmath3 carousels to a single picker . each carousel is assigned to an order of a single customer , and each order consists of exactly one item . moreover , each carousel rotates independently until the desired item reaches the picker , who is standing at a fixed point , the origin . once this position is reached , the carousel stops until the item is picked . only then will the next order be given to the carousel , which will start rotating the new order to the origin . the picker serves the carousels in a fixed order , visiting each carousel only once in every cycle . clearly , as @xmath3 goes to infinity , the utilisation of the picker in steady state tends to one , since almost surely he will never have to wait . the carousels will have brought each of their respective items to the origin by the time the picker is ready to serve them . on the other hand , the time until the picker returns to the first carousel tends to infinity ; i.e.  each individual customer suffers long waiting times . multiple carousel problems differ intrinsically from single - carousel problems in a number of ways . such systems tend to be more complicated . the system can not be viewed as a number of independently operating carousels ( cf . @xcite and section  [ ss : design issues ] ) , since there may be some interaction between two separate carousels by means of the picker that is assigned to them . namely , if the number of pickers is less than the number of carousels , then the picking strategy that is chosen for an isolated carousel may affect significantly the waiting time of another carousel . thus , one can not guarantee that minimising the travel time of a single carousel maximises the total throughput of the system ; the outcome may be quite the contrary because of the system s interdependency . another point is that in multiple - carousel problems , the i.i.d . assumption of the time needed to pick each of two consecutive orders with random item storage is in principle invalid . characteristics such as the time needed to reach the optimal point or the travel time for each carousel depend on one another through the picker s movements . for all these reasons , multiple - carousel systems merit a special reference . ideally , the problems of minimising the travel time of all carousels and maximising the picker s utilisation without surpassing certain levels of each order s response time should be studied together . however , the interdependence that appears in multiple carousel problems usually leads to complicated mathematical structures that can hardly be analysed exactly . one will have to resort to simplifications . one technique that can help overcome some of these difficulties is the setting proposed in vlasiou _ et al . _ the system we consider below consists of two carousels operated by a single picker . given a setting , i.e.  a storage scheme and a travel strategy , one first needs to obtain an estimate of the travel time needed in order to collect all items under this setting . for example , if the items are stored in random positions on the carousel , then the distribution of the travel time under the nearest - item heuristic is given by . in most settings though , this distribution can not be computed analytically , in which cases the empirical distribution or simulation may provide a partial answer . subsequently , one may need to approximate this distribution by a phase - type distribution ; see e.g.  @xcite . then , the following modelling assumption is made . we aggregate all items in one . that is , we consider an order that consists of exactly one item . it is assumed that the travel time of the carousel until that single item is reached is uniformly distributed ( i.e.  it is assumed that the item is located randomly on the carousel ) , while the distribution of the pick time for that item is taken to be equal to the phase - type distribution computed previously . under these assumptions , one can compute the utilisation of the picker by applying the results developed in vlasiou _ et al . _ this procedure can be repeated until the desired quality and efficiency regime is reached .    to describe things concretely , we consider a system consisting of two identical carousels and one picker . at each carousel there is an infinite supply of pick orders that need to be processed . the picker alternates between the two carousels , picking one order at a time . there are two ways one can view this . either , as mentioned above , one aggregates all items in an order in one super - item ( i.e.  we consider an order that consists of exactly one item ) or under the term `` picking time '' we understand the total time needed for the actual picking and travelling from the moment the picker is about to pick the first item in an order until the time the last item is picked . for ease of presentation , we will opt for the first solution , considering orders consisting of exactly one item .    as in section [ s : nelly ] , we model a carousel as a circle of length @xmath108 and we assume that it rotates in one direction at a constant speed . the picking process may be visualised as follows . when the picker is about to pick an item at one of the carousels , he may have to wait until the item is rotated in front of him . in the meantime , the other carousel rotates towards the position of the next item . after completion of the first pick the carousel is instantaneously replenished and the picker turns to the other carousel , where he may have to wait again , and so on . let the random variables @xmath109 , @xmath110 and @xmath111 , @xmath112 , denote the pick time , rotation time and waiting time for the @xmath3-th item . clearly , the waiting times @xmath111 satisfy the recursion @xmath113 where @xmath114 . we assume that both @xmath115 and @xmath116 , @xmath112 , are sequences of independent identically distributed random variables , also independent of each other . the pick times @xmath109 follow a phase - type distribution and the rotation times @xmath110 are uniformly distributed on @xmath4 ( which means that the items are randomly located on the carousels ) . then @xmath117 is a markov chain , with state space @xmath4 . moreover , it can be shown that @xmath117 is an aperiodic , recurrent harris chain , which possesses a unique equilibrium distribution . in equilibrium , equation becomes @xmath118 once the distribution of @xmath119 is computed from , we can compute @xmath120 $ ] and thus also the throughput of the system @xmath121 from @xmath122+{\mathbb{e}}[a]}.\ ] ] equation   with a plus sign instead of minus sign in front of @xmath119 at the right - hand side , is precisely lindley s equation for the stationary waiting time in a ph / u/1 queue . the equation for the standard ph / u/1 queue has no simple solution , but in vlasiou _ et al_.  @xcite we show that the waiting time of the picker in our problem can be solved for explicitly . for example , assume that the service times follow an erlang distribution with scale parameter @xmath32 and @xmath3 stages ; that is , @xmath123 and define @xmath124 $ ] . then , for the laplace transform @xmath125 of @xmath119 , i.e.@xmath126 where @xmath127 is the density of @xmath119 , the following theorem holds ( recall that the domain of integration is bounded by the length of the carousel ) .    [ 3th : erlang / u ] for all @xmath128 , the transform @xmath125 satisfies @xmath129 where @xmath130    in we still need to determine the @xmath47 unknowns @xmath131 and @xmath132 for @xmath133 . note that for any zero of the polynomial @xmath134 , the left - hand side of vanishes ( since @xmath135 is analytic everywhere ) . this implies that the right - hand side should also vanish . hence , the zeros of @xmath134 provide the equations necessary to determine the unknowns . in @xcite it is explained how to determine these unknown parameters ( which incidentally form the unique solution to a linear system of equations ) and how to invert the transform . qualitatively , the result is as follows . [ 3th : densityold ] the density of @xmath119 on @xmath136 $ ] is given by @xmath137 and @xmath138 = 1-\sum_{i=1}^{2n+2 } \frac{c_i}{r_i}(\mathrm{e}^{r_i}-1),\ ] ] where @xmath139 is a zero of the polynomial @xmath134 appearing in theorem  [ 3eq : erl_polywn ] , and where the coefficients @xmath140 are known explicitly . as a by - product , we have that    [ 3cor : tau ] the throughput @xmath121 satisfies @xmath141 + { \mathbb{e}}[w ] = \frac{n}{\lambda } + { \sum\limits}_{i=1}^{2n+2 } \frac{c_i}{r^2_i } [ 1 + ( r_i -1)\mathrm{e}^{r_i } ] .\ ] ]    the same qualitative result holds in case the pick times follow a mixed - erlang distribution . in this case , the waiting time density is again a mixture of exponentials , where all parameters can be computed explicitly ; cf . @xcite .    in a series of papers , et al_.  @xcite have relaxed several of the assumptions made above for the two - carousel setting . for example , the travel time needed to pick all items in an order can have any general distribution ( e.g.  depending on the pick strategy that is followed ) . in such cases , one can compute the distribution of the waiting time of the picker by approximating the distribution of the travel time by an appropriate phase - type distribution . phase - type distributions may be used to approximate any given distribution on @xmath4 for the travel times arbitrarily close ; see for example asmussen  @xcite . as an illustrative example , we give below the steady - state distribution of the waiting time of the picker in case the pick times follow some general distribution with laplace - stieltjes transform ( lst ) @xmath142 , and the travel times follow an erlang distribution with parameter @xmath33 and @xmath3 stages . here , @xmath135 denotes the ( unknown ) lst of the waiting time of the picker . in this case we have the following :    [ 4th : density ] the waiting - time distribution has a mass @xmath131 at the origin , which is given by @xmath143=1-\sum_{i=0}^{n-1}\frac{(-\mu)^i}{i!}\,\phi^{(i)}(\mu)\ ] ] and has a density @xmath144 on @xmath145 that is given by @xmath146 in the above expression , we have that @xmath147 and that the parameters @xmath148 for @xmath149 are the unique solution to the system of equations @xmath150    as a final curiosity , we present figure  [ fig : throughput ] . for single - server queuing models it is well - known that the mean waiting time depends ( approximately linearly ) on the squared coefficients of variation of the interarrival ( and service ) times ; see also section  [ sss : lindley ] for connections of this model to the classical single - server queue . the results in figure  [ fig : throughput ] , however , show that for this two - carousel model , the throughput @xmath121 , and thus the mean waiting time , is not very sensitive to the squared coefficient of variation of the pick time ; it indeed decreases as @xmath151 increases , but very slowly . this phenomenon may be explained by the fact that the waiting time of the server is bounded by one , that is , the time needed for a full rotation of the carousel . we refrain from giving all results derived for the waiting time distribution in this setting , as they can be found in the papers mentioned so far . one point needs to be stressed though . this technique makes usage of several simplifications ( e.g.  aggregating orders in one item ) and approximations ( e.g.  modelling various distributions as a phase - type distribution ) . some of them are almost unavoidable . for example , a carousel storing items in separate drawers should be evidently modelled with a discrete travel - time distribution ; for the application of these results though , one should approximate this distribution by a ( continuous ) phase - type distribution . however , the effect that some of these assumptions have to the final result is marginal , or at least fully controlled . as was shown in vlasiou and adan  @xcite , the error made in computing the distribution of the time the picker has to wait ( is not utilised ) is bounded . error bounds have been studied widely . the main question is to define an upper bound of the distance between the distribution in question and its approximation , that depends on the distance between the governing distributions . for our model , recall that @xmath152 , @xmath153 , and @xmath119 denote respectively the pick time needed for an item , the travel time of the carousel until this item is reached , and the waiting time of the picker until the carousel stops for the pick . moreover , @xmath154 represents the distribution of @xmath153 ( and similarly also for @xmath119 ) and @xmath155 is its approximation ( such as the phase - type approximation mentioned above ) . using this approximation , @xmath155 , one can derive analytically an exact solution that is obtained for this case for the distribution of @xmath119 . denote this solution by @xmath156 . then the following error bound holds . [ 6th : error bound ] let @xmath157 . then @xmath158)$ ] .    in the theorem above , the norm under consideration is the uniform norm . the main ingredient of the proof relies on the fact that the density for the stationary waiting time of can be described in terms of an integral equation that is a contraction mapping . as a result , approximation errors can be bounded . an almost identical result can be derived in case one approximates the pick time , rather than the travel time . thus , as this theorem indicates , resorting to approximations yields results of validity that can be controlled , provided that one has an estimation of the error that is being made by the original approximation . other results derived for the two - carousel setting include the study of the conditions under which there exists a steady - state distribution  @xcite , the study of the tail behaviour of this distribution under general assumptions for the pick and travel times  @xcite , the derivation of the steady - state distribution for various cases for the distributions of the pick and travel times  @xcite , as well as the time - dependent distribution of the waiting times of the picker for a specific setting for the distributions of the pick and travel times  @xcite . moreover , certain types of dependencies between the pick and travel times have also been studied , and the steady - state distribution has been derived for these cases as well  @xcite . it is worth a mention that such multiple - carousel systems , their mathematical peculiarities , their results and the way those are derived are not limited only to carousel , warehousing , or manufacturing problems . the same equation describing the dynamics of a two - carousel setting describes also the dynamics of a queuing model with two nodes that is applied to situations varying from a university canteen to a surgeon s operating room . for a description of such systems and detailed analysis see vlasiou _ et al . _ @xcite . what we have discussed so far on multiple - carousel problems is summarised as follows . multiple - carousel problems are intrinsically different from their single - carousel counterparts . what is of interest in such problems is to strike a balance between the utilisation of the picker and the response time of an order . to date , not much is known about such systems ; see section  [ ss : problems involving multiple carousels ] for an exhaustive literature review . a few of these results are simulation studies . however , it is almost inevitable to make use of some simulation or approximations in these problems . the results developed in vlasiou _ et al . _ @xcite help predict the performance of two - carousel systems and ultimately , combined with the results on e.g nearest - item heuristic or @xmath2-step strategies discussed in section  [ s : nelly ] , they help design a facility having a specific quality and efficiency target . however , such results are still far from accurate . more research is needed on the subject ; specific directions are provided in the next section . the mathematics and models involved in the research regarding carousel systems have surprisingly many connections to broader areas in queuing theory and applied probability . other than the relation to polling systems which will be explained in detail in section  [ polling systems ] , the subjects we have presented so far are connected to the classical single - server queue , to rendezvous networks and layered queues and even to graph theory . in the following , we highlight few of these connections . the uniform spacings defined in constitute a classical mathematical construction which is very well studied . uniform spacings have been analysed extensively in two classical review papers by pyke  @xcite . in particular , @xcite discusses the connections between uniform spacings and exponential random variables that are a main concept in the methodology presented in section  [ s : nelly ] . the markovian property ( which is also called the memoryless property ) of the exponential distribution is systematically exploited in operations research and in particular in queuing theory  @xcite . uniform spacings play an important role in mathematical statistics . mainly , they are used for goodness - of - fit tests which examine how well a sample of data agrees with a given distribution @xmath159 as its population . the idea of using uniform spacings is based on the integral transformation @xmath160 which reduces the problem to the testing of uniformity of the transformed sample . there is a vast literature on the distributions , limiting behaviour , approximations and bounds for various goodness - of - fit test statistics and empirical processes based on uniform spacings . these investigations are of great mathematical and practical interest . considerable progress in the area has been achieved in the eighties , but there are still many open problems motivating new studies .    in his detailed review , pyke  @xcite distinguishes two main types of goodness - of - fit statistics based on a function of uniform spacings , namely a sum of the form @xmath161 or a function of the ordered spacings and their ranks . the analysis of the first kind of tests goes back to le  cam  @xcite and gives rise to an extensive literature , see e.g.  @xcite and references therein . recent progress on multivariate spacings has been reported in @xcite . the second type of tests requires the knowledge of the properties of ordered spacings . this subject has been extensively studied ; we refer the interested reader to the work by deheuvels  @xcite and devroye  @xcite . an original discrete version of the problem is analysed by henze  @xcite who derives the distribution of the maximal and minimal spacings in lottery tickets .    apart from the tests mentioned above , there are also tests based on @xmath2-spacings which are the gaps between the order statistics @xmath162 and @xmath163 . for the analysis of such test statistics and their asymptotic properties as the number of observations goes to infinity , see , e.g. , del pino  @xcite , hall  @xcite , and references therein . the tests based on ordered @xmath2-spacings have been also analysed , see , e.g. ,  @xcite . more references on this subject and results on the approximations for @xmath2-spacings can be found in @xcite . for further analysis and applications of various empirical processes based on spacings see pyke  @xcite , beirlant et al . @xcite , einmahl and van zuijlen  @xcite and references therein . let @xmath164 be i.i.d . exponential random variables with mean  1 . for any @xmath165 , define @xmath166 @xmath167 note that the right - hand side of is exactly @xmath168 with @xmath169 . likewise , the right - hand side of is the minimum of two independent random variables distributed as @xmath170 . we see that the sums of independent exponentials with exponentially decreasing coefficients play an important role in the limiting results for the travel time in carousel systems as the number of items goes to infinity . specifically , these random variables appear if we consider the difference between the travel time and one complete carousel rotation , and then scale this quantity linearly with the number of items . now let @xmath171 be a standard poisson process . then we can write @xmath168 as an exponential functional associated with @xmath171 : @xmath172 the functional @xmath168 has been intensively studied in the literature . its density was obtained independently in @xcite , and in @xcite for @xmath169 . carmona  _ et al . _ @xcite derived a density of @xmath173 for a large class of functions @xmath174 , in particular , for @xmath175 . bertoin and yor  @xcite found the fractional moments of @xmath168 . if @xmath176 is a density of @xmath168 , then @xmath176 and all its derivatives equal 0 at point @xmath177 . this implies that all moments of @xmath178 are finite . however , for @xmath179 , it was proved in @xcite that @xmath180 is not determined by its moments . et al . _ @xcite found the distribution and the fractional moments of the exponential functional @xmath181 where @xmath182 is a compound poisson process . the distribution function of @xmath168 and @xmath183 has an interesting asymptotic behaviour in the neighbourhood of zero . bertoin and yor  @xcite obtained the following logarithmic asymptotics : @xmath184 where @xmath185 is a density of @xmath186 the exact asymptotic behaviour has been derived by litvak and van zwet  @xcite . compared to the logarithmic asymptotics , their formula contains several additional terms and reveals an unexpected oscillating behaviour involving theta - functions . the explanation of why the oscillations appear seems to lie in the sort of a ` binary tree structure ' of the functional @xmath187 , whose coefficients are negative powers of @xmath188 . later on , robert  @xcite and mohamed and robert  @xcite found that such oscillating asymptotic behaviour is a typical feature of algorithms with a tree structure . this phenomenon is compelling and deserves further studies . exponential functionals of poisson process and , more generally , of lvy processes , appear in a number important applications . for instance , they are relevant to the analysis of randomised algorithms  @xcite and in mathematical finance  @xcite . in @xcite and @xcite the exponential functionals of poisson processes , and , respectively , of compound poisson processes , play a key role in the analysis of the limiting behaviour of a transmission control protocol connection for the internet . we refer to the survey @xcite for further applications , results and references . the study of exponential functionals of lvy processes are a current subject of research , see e.g. @xcite , @xcite . one of the most intriguing mathematical observations that arise when studying the two - carousel model presented in section  [ s : maria ] is that recursion   differs from the original lindley s recursion  @xcite , which is @xmath189 , only in the change of a plus sign into a minus sign . at the right - hand side of these two recursions , the sign in front of @xmath111 is reversed . lindley s recursion describes the waiting time @xmath190 of a customer in a single - server queue in terms of the waiting time of the previous customer , his or her service time @xmath191 , and the interarrival time @xmath109 between them . it is one of the fundamental and most well - studied equations in queuing theory . for a detailed study of lindley s equation we refer to asmussen  @xcite , cohen  @xcite , and the references therein .    in the applied probability literature there has been a considerable amount of interest in generalisations of lindley s recursion , namely the class of markov chains , which are described by the recursion @xmath192 . the model in section  [ s : maria ] is a special case of this general recursion and it is obtained by taking @xmath193 . many structural properties of the recursion @xmath192 have been derived . for example asmussen and sigman  @xcite develop a duality theory , relating the steady - state distribution to a ruin probability associated with a risk process . for more references in this domain , see for example borovkov  @xcite and kalashnikov  @xcite . an important assumption which is often made in these studies is that the function @xmath194 is non - decreasing in its main argument @xmath195 . for example , in @xcite this assumption is crucial for their duality theory to hold . clearly , in the special case of @xmath193 which is discussed in section  [ s : maria ] , this assumption does not hold . this fact produces some surprising results when analysing the equation . the implications of this ` minor ' difference in sign are rather far reaching . for example , in section  [ s : maria ] we have presented two results in theorems  [ 3th : erlang / u ] and [ 3th : densityold ] , where we have seen that the waiting time of the picker can be solved for explicitly . for lindley s recursion , i.e.  with a plus sign instead of minus sign for @xmath119 in stationarity , this case correspond to the stationary waiting time in a classical single - server ph / u/1 queue . however , this equation has no simple solution for lindley s recursion , while we have derived a closed - form expression for the carousel recursion . furthermore , numerical results ( see also figure  [ fig : throughput ] ) show that for this carousel model the mean waiting time is not very sensitive to the coefficient of variation of the pick time , which is in complete contrast to lindley s recursion . for these reasons , we believe that it is interesting to investigate in detail the impact on the analysis of such a ` slight ' modification to the original equation . in this section , we highlight some of the differences of these two models .      for the single - server queue , i.e.  lindley s recursion , it is well - known ( * ? ? ? iii.6 ) that the random variables representing waiting times of customers converge in distribution ( and in total variation ) when the mean of the associated random walk is less than zero , or equivalently when the traffic intensity @xmath196 is less than 1 ; i.e. , when @xmath197<{\mathbb{e}}[a]$ ] , where we recall that @xmath153 is the generic service - time random variable , and @xmath152 is the generic interarrival - time random variable .    for the two - carousel model , though , which is given by recursion  , the situation is slightly different . in case @xmath198>0 $ ] , the stochastic process @xmath199 is an aperiodic , ( possibly delayed ) regenerative process with the time points where @xmath200 being the regeneration points . moreover the process has a finite mean cycle length . to see this , let @xmath201 , define the stopping time @xmath202 , and observe that a generic cycle length is stochastically bounded by @xmath121 and that @xmath203\leqslant{\mathbb{p}}[x_k > 0 \mbox { for all } k=2,\dotsc , n+1]={\mathbb{p}}[x_2 > 0]^{n}.\ ] ] moreover , we have that @xmath204<1 $ ] because of the condition @xmath198>0 \leftrightarrow { \mathbb{p}}[x<0]>0 $ ] . therefore , from the standard theory on regenerative processes it follows that the limiting distribution exists and the process converges to it in total variation . through coupling , stability can be shown also for the case where @xmath205=0 $ ] ; see @xcite for details . we see thus that while for lindley s recursion the stability condition is given by @xmath206<0 $ ] , for recursion   _ stability always holds _ ; moreover , excluding the deterministic case , we have convergence in total variation .      for lindley s recursion , there has been a substantial amount of investigations on the behaviour of @xmath207 $ ] as @xmath208 , the state of the art can be found in @xcite . results of this type for recursion   have been derived in @xcite . if the right tail of @xmath209 is regularly varying of index @xmath210 ( see @xcite for background ) , then @xmath211 { \mathbb{p}}[x > x].\ ] ] if the right tail of @xmath209 is of rapid variation ( see again @xcite ) , then @xmath212 { \mathbb{p}}[x > x].\ ] ] in both equations , we use the notational convention @xmath213 to denote @xmath214 as @xmath215 . note that the class of distributions covering these results include all phase - type distributions , as well as the weibull , gamma , lognormal and pareto distributions . moreover , these results indicate that large values of @xmath119 are caused by a single large value of @xmath30 . this is contrasting with the qualitative picture for lindley s recursion , where a large value of @xmath119 is most likely caused only by a single big jump only in the case where @xmath30 is heavy - tailed . if @xmath30 is light - tailed ( for example phase type ) , then a large value of @xmath119 is the cause of a more intricate event involving a change of measure ; see @xcite for background . a natural question is whether it is possible to unify the results for lindley s recursion and recursion  . this is possible by considering a recursion that has a minus before @xmath111 ( cf . recursion   too ) only with probability @xmath216 , @xmath217 $ ] , and has a plus before @xmath111 ( i.e.  equal to lindley s recursion ) with probability @xmath218 . for this recursion , the tail behaviour has been studied in @xcite under assumptions similar to the ones made in @xcite . to summarise the qualitative picture emerging from that paper , the tail behaviour for the unified recursion with @xmath217 $ ] converges continuously to the results for recursion   ( i.e.  if @xmath219 ) for the heavy - tailed case , while it has a discontinuity for @xmath220 ; for the so - called cramr case the result is reversed : the unified recursion is continuous for @xmath220 and discontinuous for @xmath219 , while for the intermediate case ( where @xmath30 is light tailed but does not satisfy the cramr condition ) the results for the unified recursion are continuous at both end - points . it is well known that for lindley s recursion , the time - dependent waiting - time distribution is determined by the solution of a wiener - hopf problem , see for example @xcite and @xcite . recursion   though , regularly gives rise to _ generalised _ wiener - hopf equation . for example , in @xcite we have derived a generalised wiener - hopf equation for the density of the stationary waiting time , while @xcite contains an integral equation for the generating function of the distribution of @xmath111 that is equivalent to a generalised wiener - hopf equation , which can not be solved in general . in noble  @xcite it is shown that such equations can sometimes be solved , but a general solution , as is possible for the classical wiener - hopf problem ( arising in lindley s recursion ) , seems to be absent . this makes it appear that may have a more complicated time - dependent behaviour than lindley s recursion . however , a point we make in @xcite is that this is not necessarily the case . thus , equation is a rare example of a stochastic model which allows for an _ explicit _ time - dependent analysis . the reason is that , if @xmath221 has a phase - type distribution , we can completely describe in terms of the evolution of a finite - state markov chain . we shall refrain from giving all results on the time - dependent behaviour of or their differences from the classical lindley recursion for the single - server queue , as these results have been well documented elsewhere @xcite . here , we simply list the major findings . other than deriving the time - dependent waiting time distribution for under the assumption that the random variables @xmath222 are phase - type distributed , one can derive explicit expressions for the correlation between two waiting times . it results that the covariance function @xmath223 between two waiting times with lag @xmath224 converges to zero geometrically fast in @xmath224 . this is consistent with the fact that the distribution of @xmath111 converges geometrically fast to that of @xmath119 , cf . vlasiou  @xcite . one of the properties of @xmath223 is that it is non - negative if @xmath224 is even and non - positive if @xmath224 is odd . if in addition , the random variable @xmath225 has a strictly positive density on an arbitrary interval , then the inequalities are strict . in contrast , the literature on the covariance function of the waiting times for the single - server queue seems to be sporadic . for the g / g/1 queue , daley  @xcite and blomqvist  @xcite give some general properties . in particular , in @xcite it is shown that the serial correlation coefficients of a stationary sequence of waiting times are non - negative and decrease monotonically to zero .    as we have mentioned before , @xmath117 , as given by , is a regenerative process ; regeneration occurs at times when @xmath200 . other transient results relate to the length of a generic regeneration cycle @xmath226 . for recursion , we do not need to resort to the usage of generating functions , as is necessary when analysing the corresponding quantity in lindley s recursion . note that the interpretation of @xmath226 for the carousel model is completely different from the corresponding quantity for lindley s recursion . there , @xmath226 represents the number of customers that arrived during a busy period . in the carousel setting , @xmath226 represents the number of pauses a picker has until he needs to pick two consecutive orders without any pause . in this sense @xmath226 can be seen as a `` non - busy period '' .      when deriving equation  , one of the main assumptions we have made , which led to this particular form for the equation is that the picker is not allowed to pick two consecutive orders at the same carousel and must alternate between the two carousels ( thus picking all odd - numbered orders from one carousel and all even - numbered orders from the other ) . this condition is crucial . if we remove this condition , then under certain distributional assumptions , the problem turns out to be the classical machine repair problem , and certain analogies between these two models arise .    in the machine repair problem , there is a number of machines working in parallel ( two in our situation , corresponding to the two carousels ) and one repairman ( corresponding to the picker ) , who serves the machines when they fail . the machines are working independently and as soon as a machine fails , it joins a queue formed in front of the repairman where it is served in order of arrival . a machine that is repaired is assumed to be as good as new . the machine repair problem , also known as the computer terminal model ( see for example bertsekas and gallager  @xcite ) or as the time sharing system ( see , e.g. , asmussen  @xcite or kleinrock  ( * ? ? ? * section 4.11 ) ) is a well studied problem in the literature . it is one of the key models to describe problems with a finite input population . a fairly extensive analysis of the machine repair problem can be found in takcs  ( * ? ? ? * chapter 5 ) . in @xcite we compare the two models and discuss their performance . the issue that is usually investigated in the machine repair problem is the waiting time of a machine until it becomes again operational . in the situation described in section [ s : maria ] though , we are concerned with the waiting time of the repairman . it is quite surprising that although the machine repair problem under general assumptions is thoroughly treated in the literature , this question remains unanswered . in the machine repair problem the operating time of the machine is usually more valuable than the utilisation of the repairman , which might explain why the classical literature has been mainly focused on performance measures related to the machines .    in @xcite the waiting time of the repairman is derived under the assumption that ` rotation ' times follow a phase - type distribution while ` pick ' times are generally distributed . moreover , it is shown that the random variables for the waiting time for the picker / repairman in the two models _ are not stochastically ordered_. however , on average , the alternating strategy connected to the two - carousel model leads to longer waiting times for the picker , which readily implies that the _ throughput of the machine repair model is bigger_. furthermore it is shown that the probability that the picker does _ not _ have to wait is _ larger _ in the two - carousel alternating system than in the machine repair ( i.e.  non - alternating ) model one . this result is perhaps counterintuitive , since the inequality for the mean waiting times of the picker in the two situations is reversed . regarding the relationship between the @xmath13-th waiting time of the picker in the two - carousel alternating model ( denote this by @xmath227 ) , and that of the repairman in the machine repair problem ( let this be given by @xmath228 ) , an immediate corollary of the results stated above is as follows .    for all @xmath13 , @xmath229 so , although the stationary random variables @xmath230 and @xmath231 are not stochastically ordered , the partial sums of the sequences @xmath232 and @xmath233 are . moreover , a conjecture stated in @xcite suggests that a direct application of the karlin - novikoff cut - criterion ( cf . szekli  @xcite ) leads to an _ increasing convex ordering _ , namely :    for all increasing convex functions @xmath234 , for which the mean exists , we have that @xmath235\leqslant{\mathbb{e}}[\phi(w^{\mathrm{a}})].\ ] ]      the essence of layered queueing ( a special case of which is rendezvous networks ) is a form of simultaneous resource possession @xcite .    in its most simple form in computer science applications , in a rendezvous network , a task may serve requests in two rounds ( phases ) of service . in computer applications , tasks or applications may act both as customers that needs service from other tasks and as servers to other tasks too . as a naive example , think of an application that adds up numbers . it acts both as a server , accepting requests from other applications that need numbers added , and as a customer , requiring service from the central processing unit . one can imagine that tasks are ordered in several levels or layers . tasks have directed arcs to other tasks at lower layers to represent service requests . a request from an task ( client ) to a lower - layered task ( server ) may return a reply to the requester ( a synchronous request , or rendezvous ) . while in the first phase ( i.e.  in the rendezvous ) the client is blocked and the server merely continues the thread of control of its client . however , in the second phase the client has an independent thread of control of its own . for example , task a makes a request to task b which then makes a request to task c. while task c is servicing the request from task b , tasks a and b are both blocked @xcite . among the advantages of the rendezvous is efficiency , since it provides communication without the effort of buffer management and the message copying associated with asynchronous communication . however , some potential for concurrency is lost , and there may be performance - impairing bottlenecks when a key task spends long periods send - blocked @xcite . special approximations are needed to solve queueing models which contain a two - phase server , because the second phase effectively creates a new customer in the queueing network , violating the conditions of product form queueing @xcite .    distributed or parallel software with synchronous communication via rendezvous is found in client - server systems and in proposed open distributed systems , in implementation environments such as ada , v , remote procedure call systems , in transputer systems , and in specification techniques such as csp , ccs and lotos . the delays induced by rendezvous can cause serious performance problems , which are not easy to estimate using conventional models which focus on hardware contention , or on a restricted view of the parallelism which ignores implementation constraints . stochastic rendezvous networks are queueing networks of a new type which have been proposed as a modelling framework for these systems . they incorporate the two key phenomena of _ included service _ and the _ second phase of service _ mentioned above . the main work on rendezvous networks focuses on mean value analysis and gives approximate performance estimates . this method has been applied to moderately large industrial software systems @xcite .    a layered queuing network ( lqn ) model is a canonical form for extended queueing networks that represent layered service . in a layered queue a server , while executing a service , may request a lower layer service and wait for it to complete . thus , in lqns there exist entities that have a dual role ; they act as servers to other entities of a lower layer and as customers to higher layered entities . the service time of the upper server includes the queueing delay and service time of the lower server , and this may extend through multiple layers . lqn was developed for modelling software servers , with for example blocking remote procedure calls to lower layer software servers , however it applies to any extended queueing network in which resource usages are nested , lower layer usages within higher layer usages @xcite . the two - carousel model we have presented in section  [ s : maria ] is a layered queue , and in particular a rendezvous network . to see this , organise the system as follows . the items that are stored on the carousel and have to be picked comprise the lowest layer . carousels are in the middle layer , while the picker is put in the highest layer . one may view the rotation time of a carousel as a first phase of service for the item that will be picked . the carousel ( middle layer ) acts in this case as a server . however , the second phase of service ( the actual picking ) does not necessarily happen immediately ( rendezvous ) . the item might have to wait for the picker to return from the previous carousel  cf . recursion  . at this stage , the carousels act as customers waiting to be served by the higher layer , the picker . we see thus that each carousel acts both as a server ( rotating items to the picking location ) and as a customer ( waiting until the picker completes his task before the carousel can resume its role as a server , bringing the next item to the picking location ) . layered systems are quite unknown outside the computer - science community . e.g. , in @xcite it is mentioned that `` this paper presents a model , never studied before in the queueing literature , of a system of two connected queues where customers of one queue act as the servers of the other queue ''  a comment that may very well be valid outside the computer - science literature . the analysis of recursion  , as it developed in @xcite as well as @xcite are the only papers we are aware of that deal with lqns using analytic and probabilistic tools , and admittedly all the aforementioned work on the two - carousel model had not made the connection between this model and layered queues . however unusual it might be in queuing theory to encounter a non - increasing lindley - type recursion , recursion   appears in problems involving the computation of the distribution of the maximum weight of an independent set in a sparse random graph . consider an @xmath3-node sparse random @xmath236-regular graph ( i.e.  a graph selected uniformly at random from the set of all graphs on @xmath3 nodes in which every node has degree @xmath236 ) . an independent set is a set of nodes of the graph where no two nodes in the set are connected by an edge . suppose that the nodes of the graph are equipped with some nonnegative weights @xmath237 which are generated independently according to some common distribution @xmath238 . one may be interested for example in the limits of maximum weight independent sets and matchings in sparse random graphs for some types of i.i.d . weight distributions . then recursion   corresponds exactly to the one related to the weight distribution in an 1-regular graph ; see @xcite . moreover , if one considers @xmath236-regular graphs , then the corresponding recursion giving the weight distribution in this case is similar to the one corresponding to the waiting time of a picker serving @xmath236 carousels ; see . the crucial difference in this case is in the random variables @xmath190 and @xmath111 appearing at the right - hand side of the recursion are not independent , while the corresponding variables in the recursion related to @xmath236-regular graphs are independent ; see ( * ? ? ? * eq .  ( 3 ) ) . it would be interesting to investigate the connections between the research areas of warehouse logistics and graph theory . as mentioned in section  [ s : nelly ] , as of yet the case of independent uniformly distributed items locations is the only known scenario where the travel time can be evaluated analytically by applying a systematic mathematical approach . it is important to develop methods to obtain statistical characteristics of the travel time under more realistic assumptions on the items locations . as we discuss below in section  [ ss : picking a single order ] , there are not many results in this direction in the literature . the non - uniform distributions of pick positions and especially the correlations between the items in an order lead to challenging mathematical problems . we believe that no feasible analytical solutions can be obtained in most of the realistic models . thus , the problem calls for well justified heuristics and efficient numerical methods .      the model we have considered in section  [ s : maria ] applies to a two - carousel system that is operated by a single picker . two - carousel systems have received some attention in the literature ( cf . section  [ ss : problems involving multiple carousels ] ) but many questions remain open . a line of research is directed towards studying the performance of two - carousel systems under various storage - assignment policies ( randomised or not ) , for various pick / travel time strategies and heuristics ( sequential picking , nearest - item heuristic , @xmath2-step strategies , etc . ) , for single- or dual - command cycles , and for open- and closed - loop strategies . here a single command cycle assumes a single operation , such as only storage or only retrieval . in a dual - command cycle , a storage and retrieval are combined to efficiently use the time of the operator . furthermore , an open - loop strategy implies that the carousel remains stationary at the point where the last item was retrieved ( awaiting the next order to be fed ) , while under the closed - loop strategy the carousel returns to a predefined point after the retrieval of an order is completed .    as explained in section  [ s : maria ] , two - carousel systems differ in nature and in analysis from the corresponding one - carousel problems even when studied under the same assumptions on the various storage , pick , cycle , and starting - point strategies that are followed . since two - carousel systems perform in broad terms better than single - carousel systems @xcite , studying the expected increase of the throughput of the system can help answer questions of financial nature , such as whether the benefits from the increased throughput justify the increased cost of building and operating a two - carousel system . the model discussed in section  [ s : maria ] can be extended to the case of multiple carousels as follows . for instance , consider the situation where a single picker operates three carousels . apart from the number of carousels , all other characteristics of the model remain the same as in section  [ s : maria ] . that is , we consider again an infinite queue of orders that need to be picked , we have again a rotation stage and a picking stage for each item . moreover , as before , the picker serves all carousels cyclically . for three carousels , this leads to the recursion @xmath239 where now the variables appearing at the right - hand side are not independent of one another , as was the case for all variables appearing at the right - hand side of recursion  . we may assume for convenience that the sequences @xmath115 and @xmath116 are independent among them and between them . furthermore , we note that the waiting times @xmath240 and @xmath190 are not independent . the state of the system can be modelled e.g.  as a two - dimensional markov chain , where apart from the waiting time of the picker for the @xmath3-th item that will be picked we also need to incorporate the remaining rotation time of the next carousel to be served . evidently , if the rotation times are assumed to be exponentially distributed , the system ( for three or more carousels ) can be analysed explicitly by similar techniques as the ones applied in chapter  4 of  @xcite , although it is doubtful how realistic such an assumption is . naturally , if one considers a system with multiple carousels or stations , one can think about optimisation questions . namely , as the number of carousels increases , the waiting time of the picker is expected to decrease . after serving a long series of carousels cyclically , when you return to the beginning of the cycle , with high probability the item to be picked will have reached the origin . this implies that an item will have to wait for the picker at the origin more frequently than in the two - carousel system , which means that the throughput of a single carousel decreases . intuitively , as the number of carousels increases to infinity , the utilisation of the picker increases to one , while the throughput of each individual carousel decreases to zero . given a setting , one might wonder how many carousels a single picker can operate so that we maximise both the throughput of the carousels and the utilisation of the picker simultaneously . the ultimate goal of the analysis of carousel systems is to provide a mathematical model that adequately describes the reality and , at the same time , can be efficiently evaluated either analytically or numerically . at the moment , the literature on a single carousel has advanced enough to characterise the travel time with great precision , at least for independent uniform items locations . however , as mentioned above , single carousel systems are rarely used in modern warehouses . clearly , multiple carousel models are more relevant from a practical point of view . the drawback is that such models tend to become extremely complex . until now the studies of multiple carousel systems were either solely based on simulations or employed analytical models that involved simplifying assumption on the order picking time . for instance , in section  [ s : maria ] we assumed that each order is collected within a random time that has the same distribution for each order . this is definitely a simplifying assumption , because , for instance , the orders may differ in size , and as we saw in section  [ s : nelly ] , the distribution of the travel time depends on the number of items to be collected . further literature on multiple carousels discussed in detail in section  [ ss : problems involving multiple carousels ] also involves significant simplifications of the real - life situation .    in this respect , a major challenge for future studies is to develop a unified approach for rigorous studies of real - life automated storage and retrieval systems . such an approach is expected to involve the methods proposed so far for single and multiple carousels . in sections  [ s : nelly ]  and  [ s : maria ] we presented well - developed methodologies for analytical studies of order picking in one and two carousel units . thus , an important topic for further research is to combine these two problems in one integrated study of multiple carousel systems . one may hope to obtain interesting analytical results in this direction because of the analytical nature of both methodologies . however , the problems of combining these two settings are challenging . in section  [ s : nelly ] we have seen that the travel time distribution can be of a complicated form , while the results in section  [ s : maria ] often rely on assumption such as exponential or phase - type pick times ( recall that the travel time needed to pick all items corresponds to the pick time for orders aggregated in one item ) . also , as mentioned above , the travel time depends on the size of the order , while the technique of aggregating orders in one item has made use of the assumed independence between pick times and rotation times ( while one might expect that in orders with multiple items , long travel times might be correlated to orders with multiple items and thus to shorter rotation times to the first item in the order ) . eventually , one will have to resort to the development of reasonable algorithms rather than the derivation of exact distributions . in this respect , we emphasise again that algorithmic studies of realistic carousel models constitute an important part of further research . it is also interesting to study if single or multiple carousel systems can be analysed in case there is an arrival process according to which the orders arrive . if orders arrive according to a poisson process in front of the carousel , what can be said for the waiting time of the picker ? this question can also be combined with a non - alternating system , where the picker serves the first carousel that has completed the rotation to the next item on that carousel that needs to be picked , or with bernoulli - type requests , where the picker has to serve with a certain probability the `` first '' carousel and with the complementary probability the `` other '' carousel ( potentially waiting for an item if none is present at the designated carousel ) . for each case , one should also consider the stability of the system in case the arrival rate of the orders is less than the throughput of the system with an infinite queue of orders . a polling system is comprised of a number of customer queues that are served in an order by a single server . in the literature on polling systems , the polling system with two queues where at each queue the server serves exactly one customer before switching to the other queue is often referred to as the _ 1-limited alternating - service _ model . the model described in section  [ s : maria ] is closely related to such polling systems . the two main differences are the existence of an extra stage , the rotation time of the carousel , and the absence of an arrival process for the orders . in polling systems one deals only with one stage , which in the terminology of section  [ s : maria ] is represented by the picking stage . extending the model of section  [ s : maria ] by introducing an arrival process of the orders as suggested above , is equivalent to studying an 1-limited alternating - service model with switch - over times between the stations ( which can be seen as being equivalent to the rotation time towards the single item ) .    the polling model with two queues , poisson arrivals , and no switch - over times has first been studied by eisenberg  @xcite , where the main question studied is the queue - length distribution , as is often the case in the literature on polling systems . eisenberg  @xcite gives the generating function for the stationary joint distribution of the two queue sizes . cohen and boxma  @xcite study the single server queue with two poissonian arrival streams and no switch - over times . the server handles alternatingly a customer of each queue if the queues are not empty and it is assumed that customers of the same arrival stream have the same service time distribution . it is shown that the determination of the joint queue - length distribution at the departure epoch can be formulated as a riemann - hilbert boundary problem that can be completely solved for general service time distributions . introducing switch - over times increases the complexity of the problem . in boxma @xcite the analysis is extended to include switch - over times of the server between queues , under the restriction that both queues have identical characteristics . this work is further extended in boxma and groenendijk  @xcite , where the authors no longer request that both queues have identical characteristics . it is assumed that service times and switch - over times are generally distributed .    the literature on polling systems with alternating service is not limited to the references above but is rather extensive ; see @xcite for some references . it seems though , that the question regarding the waiting time of the picker for the 1-limited polling system with two carousels has not been considered outside the scope of  @xcite . thus , introducing an arrival process for the orders in the model of section  [ s : maria ] complements the existing literature on polling systems and forms a challenging problem . the interesting feature then is that the switch - over time between two queues depends on the current picking time . again , the results from section  [ s : nelly ] can be incorporated into the model for adequate description of order picking times . an extension considered in polling systems is the _ @xmath224-limited _ service policy , where the server switches queues after having served at most @xmath224 customers in one queue . for an extensive list of references on @xmath224-limited polling systems see van vuuren and winands  @xcite . the main focus of the existing literature is again on the queue - length distribution of all stations . as the authors note in @xcite , `` to this very day , not only hardly any exact results for polling systems with the @xmath224-limited service policy have been obtained , but also their derivations give little hope for extensions to more realistic systems '' . it is worth considering the @xmath224-limited service discipline under the exact setting we have established in section  [ s : maria ] , where now the focus is on the distribution of the waiting time of the server . in the following , we classify the literature on carousels according to the main theme handled . this taxonomy allows for a better overview of the variety of the subjects examined . a crucial distinction is made between systems that involve a single carousel and systems with multiple carousels . the first four categories presented relate to single - carousel systems , while systems with multiple carousels are examined later on . the performance of a carousel system depends greatly upon the way it is loaded and the demand frequency of the items placed on it . an effective storage scheme may reduce significantly the travel time of the carousel . several strategies have been followed in practice to store items on a carousel . the simplest strategy is to place the items randomly on the carousel . _ randomised policies _ have been examined extensively @xcite , and various performance characteristics have been derived under the assumption that the items are uniformly distributed on the carousel .    one way to improve the throughput of a carousel system is to adopt a storage policy other than the randomised assignment policy . ha and hwang  @xcite have studied what they call the `` two - class - based storage '' , which is a storage scheme that divides the items in two classes based on their demand frequency . the items with a higher turnover are randomly assigned to one continuous region of the carousel , while the less frequently asked items occupy the complementary region . the authors show by simulation that the two - class - based storage can reduce significantly the expected cycle time , both in the case where a cycle is a single pick or storage of an item ( single - command cycle ) , and in the case where a cycle consists of the paired operations of storing and retrieving ( dual - command cycle ) . the same authors in @xcite examine the effects of the two - class - based storage policy on the throughput of the system , and present a case where there is a 16.29% improvement of this policy over the randomised policy . another storage scheme is suggested by stern  @xcite . assignments are made using a _ maximal adjacency principle _ , that is , two items are placed closely if their probability of appearing in the same order is high . the author evaluates this storage assignment analytically by using a markov chain model he develops . the _ organ pipe arrangement _ for a carousel system is introduced in lim _ et al . _ @xcite and is proven to be optimal in beng  @xcite and in vickson and fujimoto  @xcite under a wide variety of settings . the organ pipe arrangement has been widely used in storage units , such as magnetic tapes @xcite and warehouses @xcite . this arrangement is based on the classical mathematical work of hardy , littlewood and polya  @xcite . their concept is used in @xcite to minimise the expected distance travelled by an access head as it travels from one record to another . various optimality properties of this arrangement have been proven ; see for example keane _ et al . _ @xcite and references therein .    in carousel systems , the organ pipe arrangement places the item with the highest demand in an arbitrary bin , the items with the second and third highest demands in the bin next to the first one but from opposite sides , and sequentially all other items next to the previous ones , where the odd - numbered items according to their frequency are grouped together and placed next to one another in a decreasing order from the one side of the most frequent item ( and similarly the even - numbered items are grouped together and placed to the other side ) . figure  [ test ] illustrates the organ pipe arrangement . the numbers at the top indicate the ranking of an item in a decreasing order of frequency . park and rhee  @xcite study the system throughput and the job sojourn times under the organ pipe arrangement , where independent one - item orders arrive according to a poisson process . they explicitly quantified the gain of the organ pipe arrangement compared to random assignment and showed that this gain grows with the ` skewness ' in the demand distribution . abdel - malek and tang  @xcite study the travel times in carousels with @xmath241 bins and the organ pipe arrangement under the assumption that each order consists of one item and a sequence of orders forms a markov chain : if the current order requires bin @xmath218 then the next order requires bin @xmath242 with probability @xmath243 . the objective is to find the optimal assignment , which minimises the average travel time . their extensive numerical experiments show that although the organ - pipe arrangement is not optimal in this setting , it performs very close to optimality in a wide range of system parameters . the optimal solution in  @xcite is determined by solving a _ quadratic assignment problem_. the quadratic assignment problem is a well - known optimisation problem on choosing an optimal permutation of @xmath3 coordinates of a vector @xmath244 in order to minimise @xmath245 , where @xmath226 is a cost matrix . such problems have a long history started with the work of koopmans and beckmann  @xcite . litvak  @xcite shows by experimental studies and by providing asymptotic results for large orders that in general the optimal storage depends on the order size . moreover , the organ - pipe storage is disadvantageous when an order is large . another question related to storage is about the number of items of each type that should be stored on the carousel in order to maximise the number of orders that can be retrieved without having to reload . this question is examined in jacobs _ et al . _ @xcite , where the authors propose a heuristic that yields a reasonable solution , the error of which can be bounded . this method has been improved by yeh  @xcite , where a more accurate solution is obtained , and further on by kim  @xcite , where it is observed that the heuristic described in @xcite does not always lead to the optimal solution . the author constructs an algorithm that yields the optimal solution . this algorithm is further improved in li and wan  @xcite . this line of research has been continued in the recent paper by hassini  @xcite . in the formulation used in jacobs _ et al . _ @xcite , the author determines the optimal allocation . along with exact optimal solutions for deterministic and stochastic demand , @xcite also provides heuristics that perform close to optimal . one of the most important performance characteristics of a carousel system is the total time to pick an order . the total time to retrieve all items of an order may be expressed as a sum of the total time that the carousel is travelling plus the total time that the carousel is stopped for picking . the latter is effectively the total pick time , and it is not affected by the sequence in which we choose to retrieve the objects . however , the total travelling time greatly depends upon the retrieval sequence . the analysis of the travel time under various strategies is , in general , a non - trivial problem . this problem , however , has been resolved for independent and uniformly distributed item locations  @xcite , as we discussed in detail in section  [ s : nelly ] . various picking strategies have been proposed . bartholdi and platzman  @xcite assume a discrete model and study the performance of an algorithm and three heuristics that determine an efficient , but not necessarily optimal , sequence of retrieving all items . a heuristic is a simpler , non - optimal procedure that is based on a specific strategy . the heuristic methods proposed are the _ nearest - item heuristic _ , where the next item to be picked is always the one that is closer to the picker at any given moment , the _ shorter - direction heuristic _ , where the carousel chooses the shortest direction between the route that simply rotates clockwise and the route that rotates counterclockwise , and the _ monomaniacal heuristic _ , that always chooses to rotate to the right and pick items sequentially . the _ optimal retrieval _ algorithm that is presented enumerates all possible paths ; therefore , it is guaranteed to find the quickest sequence in which to retrieve a single order . in @xcite the authors prove among other things that the travel time under the nearest - item heuristic is never greater than one rotation of the carousel . et al . _ @xcite provide the upper bound of @xmath246 full rotations , where @xmath3 is the number of items in the order , and show that the new upper bound is tight . litvak and adan  @xcite obtained the distribution and the asymptotic properties of the travel time under the nearest - item heuristic for uniformly distributed independent items locations . these results , based on properties of uniform spacings , have been discussed in detail in sections  [ ss : narest - item ] ,  [ ss : asymptotic results ] . in @xcite , the first two moments of the travel time and the distribution of the number of turns are computed recursively by conditioning on the event that there is an empty space of size @xmath247 on one side of the picker s current position . we presume that such methods may lead to the travel time distribution in some special cases with non - uniform items locations . another interesting picking strategy that has been already discussed in section  [ ss : m - step ] is the so - called _ @xmath2-step strategy _ , where the carousel chooses the shortest route among the ones that change direction at most once , and only do so after collecting at most @xmath2 items . in case of independent uniformly distributed items locations the average travel time under the @xmath2-step strategy is smaller than the one under the nearest - item heuristic already for @xmath96 ; see @xcite . the results by litvak and adan  @xcite on the @xmath2-step strategies have been presented in section  [ ss : m - step ] . in an earlier paper , rouwenhorst _ et al . _ @xcite apply analytical methods to study the case when @xmath85 . this means that the carousel changes direction after collecting at most two items . they interpret @xmath2-step strategies as stochastic upper bounds for the minimal travel time and present convincing numerical results on the excellent performance of such strategies . wen and chang  @xcite model the carousel as a discrete bidirectional loop and assume that the time to move between the bins of a shelf is _ not _ negligible . they propose three heuristic solution procedures and compare their performance . an earlier version of this work can be found in wen  @xcite . ghosh and wells  @xcite model the carousel as a continuum of _ clusters _ and _ gaps _ , where a cluster is a segment on the circle that corresponds to a series of locations that have to be visited for the retrieval of an order , while a gap is the segment of the circle between two clusters . the authors develop two algorithms to find optimal retrieval strategies . in particular , to find an optimal path , they avoid a complete enumeration by noticing that a turn can never be made after covering more than 1/3 length of the carousel . stern  @xcite studies properties of the optimal , i.e.minimal , picking sequence both for the _ open - loop _ strategy , where the carousel remains stationary at the point where the last item was retrieved ( awaiting the next order to be fed ) , and for the _ closed - loop _ strategy , where the carousel returns to a predefined point after the retrieval of an order is completed . he formally shows that under the open - loop strategy the carousel will change its direction at most once when following the optimal picking sequence , while under the closed - loop strategy the carousel will turn at most twice . a recursive expression for the distribution of the minimal travel time needed to collect one order of @xmath3 randomly distributed items in the open - loop scenario is given explicitly by litvak and van zwet  @xcite . the case when positions of the items in an order are dependent has not received much studies . one way to model the dependencies is described by abdel - malek and tang  @xcite who assume that the positions of successive items form a markov chain . in this setting , they study the performance of the organ - pipe storage rule . stern  @xcite introduces correlations between items in an order by considering several order types , where each type corresponds to a fixed list of items . the work of wan and wolff  @xcite focuses on minimising the travel time for `` clumpy '' orders , that is , orders concentrated on a relatively small segment of the carousel , and introduces the nearest - endpoint heuristic for which they obtain conditions for it to be optimal . in this setting , one can no longer assume that the items locations are uniformly distributed . moreover , there is clearly a strong dependence between items positions .    the model with non - uniform items locations reflects a relevant situation when some of the items are required more frequently than others . most of the papers that assume distinct frequencies assume the orders of one item ( see e.g.  @xcite ) . an interesting work on non - uniformly distributed items is given by litvak  @xcite , where the focus is on the length of the shortest rotation time needed to collect a single order when the order size is large and the items locations have a non - uniform continuous distribution with a positive density @xmath107 on @xmath136 $ ] . a popular strategy for reducing the mean travel time per order in carousel storage and retrieval systems is batching together a number of orders and then picking them sequentially . batch _ is a set of orders that is picked in a single tour . two consecutively picked items do not necessarily belong to the same order . an excellent literature survey by van den berg  @xcite on planning and control of warehousing systems addresses this issue and the problems that arise if large batches are formed . apart from the questions mentioned before , stern  @xcite also considers the performance of a carousel for a fixed set of order types ( for example , big orders with many items , and small ones ) . bartholdi and platzman  @xcite are mainly concerned with sequencing batches of requests in a bidirectional carousel . they specify the number of orders to be retrieved ( ignoring any new arrivals ) and propose three heuristic methods to solve this static problem . orders may be picked in any sequence ( and not necessarily at the order they arrive ) , but picks within the same order are performed consecutively . they define the _ minimum spanning interval _ , which is the shortest interval containing all the items of an order and , by assuming that the picker always begins and finishes retrieving an order at one of the endpoints of this interval , they construct the shortest matching chain by ordering the orders accordingly . this procedure may fail to give an uninterrupted sequence in which to pick the orders ; therefore , they propose the following heuristics . the first one , called the _ hierarchical _ heuristic , picks any order that happens to have a common endpoint with another order , and then travels clockwise until an unpicked endpoint is encountered , and repeats the procedure . the _ nearest - order _ heuristic is practically an extension of the nearest - item heuristic described earlier in the paper , as is the case with the _ second monomaniacal _ heuristic they propose . under these heuristics , they obtain upper bounds for the travel time . ghosh and wells  @xcite assume that the orders have to be picked under a fifo sequencing restriction , which means that the first order to arrive at the warehouse is the first order that will be picked , and so on . since the orders are retrieved in a fifo fashion , the problem is reduced to finding how to retrieve each individual order so that the best overall retrieval is achieved . they develop an algorithm for the optimal retrieval path of @xmath3 orders via dynamic programming , and show how to update dynamically the solution when new orders arrive . rouwenhorst _ et al . _  @xcite model the carousel as an m / g/1 queuing system , where the orders are the `` customers '' that require service , and the service they get depends on the pick strategy that is followed . this approach permits the derivation of various queuing characteristics such as the mean response time and the waiting time when orders arrive randomly . the authors mention that the tight upper bounds for the mean response time can be further exploited to obtain also good approximations for excess probabilities of the response time . van den berg  @xcite assumes either a fixed or an arbitrary sequence of orders . when the sequence of the orders is given , he presents an efficient dynamic programming algorithm that finds an optimum path that visits all orders in the specified sequence . furthermore , when there is no given order sequence , he simplifies the problem to a _ rural postman problem _ on a circle and solves this problem to optimality . the rural postman problem is the problem of finding the shortest route in an undirected graph which includes all edges at least one time . van den berg  @xcite concludes that the obtained solution requires at most 1.5 revolutions more than a lower bound of an optimal solution to the original problem . simulation results suggest that the average rotation time may be reduced up to 25% when allowing a free order sequence . lee and kuo  @xcite formulate the problem of optimal sequencing of items and orders as a _ multi - travelling salesman problem_. in the multi - travelling salesman problem , there are several salesmen in a home city , and each of the other cities has to be visited only by one salesman . using this formulation , lee and kuo  @xcite provide efficient heuristics for optimal picking of several orders consisting of multiple items . all research papers mentioned so far that deal with travel time models of carousel systems assume average uniform velocity of the carousel . in other words , the main assumption is that the carousel travels with constant speed and the acceleration from the stationary position ( when a pick is performed ) to its full speed , as well as the deceleration from the maximum speed to zero speed , are negligible factors when computing the travel time of the carousel . guenov and raeside  @xcite give some empirical evidence that the error induced when neglecting acceleration and deceleration of an order picking vehicle is indeed negligible . thus the problem of minimising retrieval times can be considered to be equivalent to the problem of minimising the average distance travelled by the carousel per retrieval . hwang _ et al . _ @xcite , however , develop strategies for picking that take into consideration the variation in speed of the carousel . for unit - load automated storage and retrieval systems there are several travel - time models that consider the speed profiles of the storage and retrieval robot . in @xcite some relevant references are given . unlike the unit - load automated storage and retrieval systems , almost all the existing travel - time models for carousel systems assume that the effects of the variation in speed are negligible . in @xcite the authors try to bridge this gap in the literature . they assume that the items are randomly distributed on the carousel and derive the expected travel time both in the case of a single command cycle and in the case of a dual command cycle . they verify the accuracy of the proposed models by comparing the results to results directly obtained from discrete racks . egbelu and wu  @xcite study the problem of pre - positioning the carousel in anticipation of storage or retrieval requests in order to improve the average response time of the system . choosing the right starting point of a carousel in anticipation of an order is also referred as the _ dwell point selection problem . this strategy becomes relevant when the items are stored under the organ pipe arrangement . in this situation the dwell point should be chosen to be the location of the most popular item ; see , e.g. , @xcite . spee  @xcite is concerned with developing design criteria for carousels . he states the basic conditions for designing an automatic order picking system with carousels and comments on the optimal storage design . namely , he is interested in finding the right number of picking robots and the right number and dimensions of a carousel so that the investment is minimised , provided that the size of the orders that need to be retrieved is given . mcginnis  @xcite studies some of the design and control issues relevant to _ rotary racks_. a rotary rack is an automated storage and retrieval system that strongly resembles carousels . in fact , conceptually , a rotary rack is simply a carousel , where the only difference is that each level or shelf of this carousel can rotate independently of the others . the author concludes that , while rotary racks appear to be a simple generalisation of conventional carousels , the control strategies that have been shown effective for carousels do not appear to be as effective for these systems . rotary racks can be viewed as a multiple - carousel system ( where each level is considered as a sub - carousel ) with a single picker . while almost all work mentioned in this section concerns one - carousel models , real applications have triggered the study of models involving multiple carousels . the study of such models is not as developed yet as the study of models involving a single carousel . the list of references that follows seems to be complete . perhaps the first academic study that investigates the performance of a system involving several carousels is that of emerson and schmatz  @xcite . the authors simulated the operation of the warehouse of rockwell s collins telecommunications products . the system consists of twenty - two carousels , where each pair of carousels had a single - operator station ( so there are in total eleven operator stations ) . the questions they are concerned with are how big the batch size of orders should be so as to complete the week s work ( which is used as a performance measure ) and keep all operators busy , what happens when a carousel or a station is down , and how is an overload or an imbalance ( for example , unequal operator performance , unequal carousel loading , or large orders ) handled . in order to investigate potential solutions to these three imbalance conditions , the authors investigate two operating rules . the first operating rule studies six different storage schemes with seven carousel pairs ( and thus seven operators ) . it uses simulation models to study simple storage schemes such as random storage , sequential alternating storage , and storage in the carousel with the largest number of openings . the aim in @xcite is to study the degree of carousel usage . the authors find that there is no significant difference between the carousel loads among the storage schemes . however , they do not treat the problem of optimally assigning items to carousel bins , and do not present any analytical models to help investigate the problem . the second operating rule they investigate is a floating operator . this is an operator who is trained to work at any station , and who is moving to different stations according to specific needs ( for example , depending on the size of the queue at a particular station ) . they conclude that this solution seems advantageous for the purposes of the warehouse they investigate . koenigsberg  @xcite presents analytic solutions for evaluating the performance of a single carousel , and discusses the ways in which his approach can be extended to a system involving two unidirectional carousels both served by a single robotic operator . the carousels are related only through the state of the robot , which means that each carousel is independent of the other except for the time it waits for an operation to commence ( such as pick , storage , or repair ) because the robot is busy at the other carousel . the author concludes that under some conditions , it is often more advantageous to have two carousels of identical length instead of one carousel of double the length . furthermore , going to three carousels of equal length ( i.e.  one third of the length of the single carousel ) will offer little further improvement . hwang and ha  @xcite study the throughput performance both of a single and of a double carousel system . based on a randomised storage assignment policy , cycle time models are developed for single and dual commands . furthermore , they examine the value of the information on the succeeding orders in terms of system efficiency , which may lead to better scheduling of the orders to be processed .    in a later work , et al . _ @xcite attempt to measure analytically the effects of double shuttles of the storage and retrieval machine ( i.e.  the robotic picker ) on the throughput both of the standard and of the double carousel system . storage and retrieval machines with double shuttles are machines that have space for two items . thus , for example , an item can be retrieved from the carousel and stored on one shuttle , while the other shuttle has an item that needs to be stored to the carousel . after this item is stored , a second item can be retrieved from the carousel and placed on the empty shuttle . all these operations occur during a single cycle of the carousel operation . for the double carousel system , a retrieval sequence rule is proposed which utilises the characteristics of the two independently rotating carousels . from the test results , double shuttles are shown to have a substantial improvement over single shuttles . this improvement tends to be more prominent in the double carousel system . due to cost concerns , the authors note that an economic evaluation will be needed to justify the extra cost of double carousel systems and double shuttles before implementing them in real world situations . wen _ et al . _ @xcite consider a system comprised of two carousels and a single retrieval machine . their main assumption is that every order must be picked in a single tour , i.e. , an order can not be divided into two or more sub - tours . batching orders together is also not allowed . they analyse the retrieval time and propose four heuristic algorithms for the scheduling sequence of retrieving items from the system to satisfy an order . their method is an extension of the algorithm presented in @xcite and @xcite . meller and klote study the throughput of a group of several carousels , a so - called carousel pod  @xcite . they use approximations to evaluate the order pick time in one carousel and then evaluate the throughput of a pod by plugging in the average response times of each unit and modelling the pod of @xmath248 carousels as a queuing system where @xmath249 picker operates one carousel . further , they derive an approximation for the system s throughput using a diffusion approximation by gelenbe  @xcite which was earlier applied by bozer and white  @xcite in the analysis of end - of - aisle order - picking systems . recently , hassini and vickson  @xcite studied storage locations for items , aiming to minimise the long - run expected travel time in a two - carousel setting with a single picker . they assume that the products are available at all times ( so as to be able to ignore possible delays due to lack of stock ) , and that orders are not batched ; that is , the carousel system processes only single - item orders . this is applicable in situations where individual product orders are processed in a first - come - first - served policy , or when the next item to be retrieved is known only after the present one has been picked . the authors compare the performance of three heuristic storage schemes and a genetic algorithm  @xcite that for small - sized problems completely enumerates the solution space . they conclude that none of the heuristic approaches leads to a solution that outperforms the algorithmic solution they provide . the same model is also studied by park _ et al . _ as is the case in @xcite , in @xcite the basic assumptions are that there is an infinite number of items to be picked and that an order consists of a single item . the authors , however , are not interested in storage issues . they further assume that the single operator , the picker , is alternately serving the two carousels . this may cause the picker to have to wait for an amount of time until the item at the carousel he is currently serving is rotated in front of him . they derive the distribution of the waiting time of the picker under specific assumptions for the pick times . this allows them to derive expressions for the system throughput and the picker utilisation . the model presented in @xcite has been extended further in vlasiou _ _  @xcite by removing all assumptions related to the pick times or rotation times . in related work , vlasiou _ et al . _ @xcite have shown that the two - carousel model studied in @xcite is equivalent to an alternating service queue , if one allows for rotation times with an infinite support . some of these results have been presented in section  [ s : maria ] . finally , we would like to mention that there is a broad literature on automated storage and retrieval systems ( see e.g. the survey by le - duc  @xcite ) . an extensive list of references has been also made available on - line by roodbergen  @xcite . the authors would like to thank ivo adan for the suggestion to write this survey paper and for his active involvement in the research described in sections  [ s : nelly ]  and  [ s : maria ] . in the netherlands , the three universities of technology have formed the 3tu federation . this article is the result of joint research in the 3tu centre of competence nirict ( netherlands institute for research on ict ) . poissonian exponential functionals , @xmath242-series , @xmath242-integrals , and the moment problem for log - normal distributions . in _ seminar on stochastic analysis , random fields and applications iv_. vol . 58 of _ progress in probability_. birkhuser , basel pp . 4556 . ( 1988 ) . two queues with alternating service and switching times . in _ queueing theory and its applications  liber amicorum for j.w . cohen_. ed . o.  boxma and r.  syski . vol .  7 of _ cwi monographs_. north - holland , amsterdam pp .  261282 . ( 1997 ) . on the distribution and asymptotic results for exponential functionals of lvy processes . in _ exponential functionals and principal values related to brownian motion rev . mat . iberoamericana , rev . mat . iberoamericana , madrid pp . 73130 . real time optimisation of man on board order picking . in _ proceedings of the 10th international conference on automation in warehousing_. ed . j.  white . ifs(publications ) , dallas , texas . .  8793 . stability bounds for queueing models in terms of weighted metrics . in _ analytic methods in applied probability_. ed . y.  suhov . vol . 207 of _ american mathematical society translations ser .  2_. american mathematical society , providence , ri pp . 7790 . analysis of the efficiency of carousel and tote - stacker performance . in _ proceedings of the 7th international conference on automation in warehousing_. ed . j.  white . springer , san francisco , california . .  173183 . analysis of multi - server systems via dimensionality reduction of markov chains . . school of computer science , carnegie mellon university pittsburgh , pa 15213 . available at http://www.cs.cmu.edu/~osogami/[http://www.cs . cmu.edu/@xmath252osogami/ ] . spacings revisited . in _ proceedings of the sixth berkeley symposium on mathematical statistics and probability ( univ . california , berkeley , calif . , 1970/1971 ) , berkeley , calif . , 1972 . california press , vol . i : theory of statistics _ , pp . 417427 . performance analysis of a carousel system . in _ progress in material handling research : 1996_. ed . r.  j. graven , l.  f. mcginnis , d.  j. medeiros , r.  e. ward , and m.  r. wilhelm . the material handling institute , charlotte , nc pp .  495511 . automatic order picking system with horizontal racks . in _ progress in material handling research : 1996_. ed . r.  j. graven , l.  f. mcginnis , d.  j. medeiros , r.  e. ward , and m.  r. wilhelm . the material handling institute , charlotte , nc pp .  545550 . . parts location and optimal picking rules for a carousel conveyor automatic storage and retrieval system . in _ proceedings of the 7th international conference on automation in warehousing_. ed . j.  white . springer , san francisco , california . .  185193 . the order picking problem in the carousel systems . in _ proceedings of the 7th international conference on automation in warehousing_. ed . j.  white . springer , san francisco , california . .  195199 . order picking for a two - carousel - single - server system in an automated warehouse . in proceedings of the 10th international conference on automation in warehousing_. ed . j.  white . ifs(publications ) , dallas , texas . .  8793 .
Abstract:

Answer: this paper gives an overview of recent research on the performance evaluation and design of carousel systems . we discuss picking strategies for problems involving one carousel , consider the throughput of the system for problems involving two carousels , give an overview of related problems in this area , and present an extensive literature review . emphasis has been given on future research directions in this area . @xmath0 faculty of electrical engineering , mathematics and computer science , + department of applied mathematics , university of twente , + 7500 ae enschede , the netherlands . @xmath1 eurandom and department of mathematics & computer science , + eindhoven university of technology , + p.o . box 513 , 5600 mb eindhoven , the netherlands . + n.litvak@ewi.utwente.nl , m.vlasiou@tue.nl    * keywords : * order picking , carousels systems , travel time , throughput    * ams subject classification : * 90b05 , 90b15