Teams to exploit spatial locality among agents

Size: px
Start display at page:

Download "Teams to exploit spatial locality among agents"

Transcription

1 Teams to exploit spatial locality amon aents James Parker and Maria Gini Department of Computer Science and Enineerin, University of Minnesota Abstract In many situations, task allocation is hihly dependent on the spatial locality of the aents and the tasks. In this paper we explore task allocation in complex spatial environments, where the cost of completin a task can vary over time. This work focuses on a subset of cost functions for tasks that row over time and presents the real-time Latest Finishin First (RT- LFF) alorithm, which strikes a balance between the optimal zero travel time solution and minimizin the amount of time aents spend transferrin to tasks in different locations. We also present how spatial knowlede can be used to infer missin knowlede in partially observable spaces. I. INTRODUCTION Multi-aent systems offer robustness, flexibility, and efficiency over sinle-aent systems. The fields of wide-area surveillance, exploration and mappin, transportation, and search & rescue are all bein enriched by incorporatin multiple aents. In all these examples, not only must aents travel throuh the spatial topoloy of the environment to accomplish their oals, but the topoloy itself can effect how the aent s oals can chane over time. However, the benefits of multiaent systems increases the burden on coordination. In dynamic environments where the set of aents miht chane over time, as aents fail or new aents are added to the system, and where new tasks can appear and old tasks can expire, task allocation has to be done in real-time to handle the uncertainty as to location, completion requirements and size of tasks. Also aents travel throuh unknown and possibly unsafe areas to reach the tasks. Unlike most task allocation problems [1], we focus on domains where the amount of resources that a task requires chanes over time. We present the Latest Finishin First (LFF) alorithm for task allocation in domains where the cost of tasks increases with time (Section IV). These types of problems can occur in nature, such as invasive species or forest fires, where if a task is not completed quickly, then it can become difficult or impossible later. We then extend our LFF alorithm to a real-time heuristic version for partially observable dynamic environments where new tasks can appear as time proresses (Section V). Finally, we propose a novel way of incorporatin partial information in RoboCup Rescue to accurately predict a model for clustered tasks (Section VI), and we evaluate experimentally our work aainst other methods (Section VII). II. RELATED WORK Multi aent task allocation is well known to be an NP hard problem, ivin rise to many different techniques for findin an approximate solution. One example is swarm robotics, where aents often use threshold based methods to individually assess the constraints and their ability to complete each task. If an aent s abilities surpass a threshold on the constraints, then the aent allocates itself to the task. If not, the aent passes the information to other aents. An example is [2], which uses distributed constraint optimization (DCOP) as a basis for task allocation. A comparison between DCOP and other swarm techniques is provided in [3]. In comparison, market inspired auction methods typically require more communication and are more centralized. Zhan et al. [4] present an auction based approach to form executable coalitions, allowin multiple aents from different locations to reach a task and compete it efficiently. Recently, decentralized applications have been desined that add flexibility to the system (e.., [5]). Our work strikes a balance between swarm and centralization, each aent is directed to an area by a central authority, but upon reachin the destination, aents act on their own individual abilities. Other approaches have been developed, such as modelin task allocation as a potential ame [6]. Sandholm et al. [7] present a eneralized coalition formation alorithm which produces solutions within a bound from the optimal via prunin a subset of the search tree. Dan et al. [8] improve on it by further prunin, but the search time is still exponential. The work in [9] focuses on tasks that require more than one aent to do them, while simultaneously tryin to use efficiently the aent s resources and time. Our approach also assumes multiple aents are required, but we also allow the requirements of tasks to chane over time. Our work specifically applies to urban search & rescue usin the RoboCup Search and Rescue Simulator [10], providin simulations based off street and buildin maps of real cities. Emerency situations are very time critical and often lackin in information, as demonstrated by [11]. Most notably, when the emerency occurs aents are spatially spread out in a disoranized fashion and must quickly coordinate and form teams to accomplish tasks. III. PROBLEM DESCRIPTION Our problem is task allocation for multi-aent systems, where aents must travel on desinated pathways. The aents are scattered at random initially and need to convere on tasks. These tasks require more work to be completed as time proresses, which necessitates multiple aents bein allocated to them. Thus, simply oin to the nearest task can cause overallocation and other tasks miht row out of control. We denote the set of homoeneous aents by A = {a 1,..., a A } and the set of tasks by B = {b 1,..., b B }. The set of assinments of aents to a task is denoted by N(t) = {n 1 (t),..., n B (t)}, where n i (t) is the set of aents from A that are assined to task b i at time unit t. An aent a i A can only work on one task, b j, at a time and at most one assinment n j (t) N can contain a i. All aents and tasks

2 have a spatial location in the environment and the travel time between two locations is assumed to be computable. Each aent provides the amount of work w per time unit towards the task it is assined. Every task b i B has a cost function fi t t 1 : (fi, n i (t 1) ) R with the relationship: where f t i f t+1 i (fi t, n i (t) ) = fi t (f t 1 i, n i (t 1) ) + fi t (1) has the form: f t i = i (f t i ) w n i (t) (2) and i : R >0 R >0 is a monotonically increasin function with initial value fi 0. For simplicity, we will treat f i t as a value rather than a function if the values of f t 1 i and n i (t 1) are clear in the context. This means fi t is strictly monotonically increasin when i (fi t) > w n i(t) and strictly monotonically decreasin when i (fi t) < w n i(t). When the cost function fi t reaches or passes zero at some time t i (0), the task is considered complete and is removed from the problem. In our problem the cost of completin a task increases over time, therefore the oal of task allocation is to minimize the time the last task is completed. IV. LATEST FINISHING FIRST With perfect knowlede about how a task rows, we can predict what allocation of aents would be ideal. This ideal solution is not reachable due to travel time bein required to reach tasks, but even in this ideal settin task allocation is not as simple as it seems. Even nelectin the spatial limitations, when tasks row over time some allocations of aents are better than others. Takin into account the spatial limitations, we derive the Latest Finishin First (LFF) alorithm, which we present after first describin properties of the dynamic cost functions we will use. Throuhout this section, we assume that time t is sufficiently small to enable us to treat t as continuous and we assume a solution exists. T T (x, y) is defined to be the travel time between location x and location y. A solution is found at some time t s if and only if all b i B have f ts i = 0. Usin (2), we can write the previous equation as fi 0 + t<t s fi t = 0. Since this is true for every b i B a solution is reached if and only if: ( fi 0 + ( i (fi t ) w n i (t) )) = 0 (3) t<t s Note that b f 0 i B i is a constant and t<t s w n i (t) = p t s w A, where p is the overall percent of time the aents are workin. We can then rewrite (3) as: fi 0 + i (fi t ) p t s w A = 0 (4) t<t s This means that t<t s i (fi t ) is the amount of work added to the system from time t = 0, which we call R ts or the lobal reret. We define R t = R t R t 1 = b i B i(fi t). Theorem 1: Minimizin R ts is an optimal solution when T T (x, y) = 0 x, y. Proof: When T T (x, y) = 0 x, y we can assume p = 1, since aents can instantly move between tasks. Suppose a better solution ˆf i t exists, then ˆR tŝ = t<tŝ i ( ˆf i t) where Ŝ < S with ˆR tŝ > R ts. Rewritin (4) for R ts yields: = w A t s R ts Solvin (4) for ˆR tŝ in terms of t s : ˆR tŝ + w A (t s tŝ) = w A t s Usin our two assumptions, we can then write: f 0 i R ts + w A (t s tŝ) < ˆR tŝ + w A (t s tŝ) which implies that R ts satisfied (3) at time tŝ. This is a contradiction. Theorem 2: Minimizin R ts for every time unit t also minimizes the lobal reret, R ts, when 1. T T (x, y) = 0 x, y, 2. i = j i, j, 3. lim x 0 + (x) = 0, and 4. 2 t 2 i (fi t) 0. Proof: Since i = j i, j, we will denote i with to simplify notation. Assume there exists a better solution ˆF i = { ˆf i 1, ˆf i 2,...} which differs from the reedy minimization F i = {fi 1, f i 2,...}. This means that at some time t d the better solution must assin aents differently than the reedy solution. By definition the reedy minimization is a minimum R td at time t d, thus R td ˆR td, where ˆR t = b ( ˆf t i B i ). After time t d, the reedy minimization will allocate aents exactly in the same way as the better solution and because T T (x, y) = 0 x, y this is possible from any confiuration. Next we prove by induction that b f t i B i ˆf b t i B i. At time t d, b f t d i B i = ˆf t d i combined with the fact that fi x = fi 0 + t<x f i x alon with F i is a reedy choice implies b f t d+1 i B i ˆf t d+1 i, which is the base case in the induction. If we write ˆf i t = f i t+c i, then we can conclude b c i B i 0. We can then compute for f t+1 i as: and ˆf t+1 i f t+1 i = f t i + (f i ) w n i (t) as (n i (t) is the same since assinments are copied): ˆf t+1 i = (f t i + c i ) + (f i + c i ) w n i (t) If we use the monotonicity of, then we can see (f i + c i ) = t+1 (f i ) + δ i c i for some δ i > 0. This means ˆf i f t+1 i = c i + δ i c i or rearranin: f t+1 t+1 i + c i + δ i c i = ˆf i. Since 2 t 2 i (fi t) 0 we know δ i δ j when c i > c j i, j. Thus, b δ i B i c i 0 since b c i B i 0. We can conclude that b f t+1 i B i t+1 ˆf i, the inductive step. Usin a similar loic to extractin the definition of R t from (3), we drop b f 0 i B i and t<t s b w n i B i(t) from both sides and et: R ts ˆR ts where ˆR ts = t<t s ( ˆf i t ). This is a contradiction to the fact that ˆF is a better solution, therefore minimizin R t at every time t also minimizes R ts. Before concludin the proof, we must consider the cases when F completes a task that ˆF did not and vice versa. If F completes a task that ˆF did not, then the reedy minimization f 0 i

3 will assin aents from an already completed task to a random unfinished tasks. This does not invalidate any of the inequalities above. When ˆF completes a task that F has not, this task will never be completed by direct mimicry from the reedy minimization solution, instead this will be finished by the random assinment described above. There is no discontinuity in the sums since we require lim x 0 + (x) = 0, so when a task is completed it simply disappears from the equations. Unfortunately, the assumption that T T (x, y) = 0 x, y is rather stron. When aents require time to move between different task clusters, every time an aent is allocated to a different task, p decreases in (4). To address this issue, we introduce Latest Finishin First (LFF) in Alorithm 1. The inspiration behind LFF is to create a stable assinment that will try to maintain p close to 1 to maximize the overall output of aents in the system. To do this, each task is prioritized and the closest aent to the hihest priority task is assined to that task. After an aent has been assined, the priority of the hihest priority task is recomputed and this process is repeated. The priority of a task is the expected time to complete that task, so tasks that take loner have a hiher priority. Since i increases work needed to complete a task, all tasks initially start out never completin (or t i (0) =, i). In case of ties, the task which has the larest current cost is allocated the closest aent. This causes the completion times of all the tasks to be as balanced as possible, which in turn reduces the amount of travel needed for aents. Alorithm 1 Latest Finishin First 1: while a i A and j such that a i n j (0) do 2: Find b i for armax t i (0) {Task b i ends last} 3: if b i, b j where t i (0) = and t j (0) = then 4: Find b k for armax fk 0 with t k(0) = b k B 5: Assin armin T T (a j, b k ) a j A and a j unassined 6: else 7: Assin armin T T (a j, b i ) a j A and a j unassined 8: end if 9: end while V. REAL-TIME LATEST FINISHING FIRST This section extends LFF outlined in Alorithm 1 to a decentralized real-time solution for dynamic partially observable environments in Alorithm 2. The LFF alorithm assumed all the tasks were known initially, but without this assumption it becomes impractical to try and minimize aent movement. When new tasks are discovered, it is crucial that aents are reallocated to the new tasks. The Real-Time Latest Finishin First (RT-LFF) alorithm assumes that new tasks are identified from an oracle, and when a new task is identified aents are reallocated to try finish all tasks at the same time. When the task allocation first starts, the normal LFF is used to compute the assinments of aents. While time is proressin, if a new task, b j, is identified then RT-LFF uses a reedy heuristic to reallocate aents to b j. Tasks are ordered based on spatial proximity to b j, and tasks that are closest to the new task attempt to reallocate aents first. Let t i (0) be the time for task b i to complete with one fewer aent and t + j (0) be the time for task b j to complete with one additional aent includin the delay from the aent travelin. An aent is only reallocated from b i to b j if t i (0) < t+ j (0), namely a reedy heuristic which only transfers aents if the oriinal task will still finish first even after reallocatin an aent. If the oracle informs all aents of the new task, each aent can compute this alorithm and aree on which aents should transfer without any communication. In the case where aents have limited computational power, a sinle more powerful aent at or nearby each task could direct this alorithm. If the oracle only informs the aents closest to the newly identified task, these aents would run the RT-LFF alorithm, since they are the hihest priority, and after any assinments pass on updated information to the next closest aents. Theorem 3: RT-LFF has at most B 1 misallocated aents from LFF at any time assumin zero travel time. Proof: We show that usin RT-LFF a task b 1 will never be assined more than one aent when usin LFF by breakin it down into two cases. First, assume that we have two tasks b 1 and b 2 and let t ++ j (0) be the time task b j is completed with two additional aents assined to it and similarly t i (0) be the time task b i is completed with two less aents assined to it. Without loss of enerality assume t 2 (0) < t 1 (0). If RT-LFF did not reallocate an aent from b 2 to b 1, then t 2 (0) > t+ 1 (0). We notice that reallocatin two aents from b 2 to b 1 means t 2 (0) > t 2 (0) and t++ 1 (0) < t + 1 (0), which is a contradiction to the RT-LFF alorithm since b 1 would allocate an aent back to b 2 since t 2 (0) > t+ 1 (0). Next we show that a task b 1 would never receive a sinle aent from more than one cluster. The arument is similar to the first part, only it now involves task b 3. If b 2 and b 3 did not reallocate an aent to b 1 under RT-LFF, then t + 1 (0) < t 2 (0) and t + 1 (0) < t 3 (0). Suppose LFF did have both b 2 and b 3 reallocate an aent to b 1, then without a loss of enerality assume t 2 (0) < t 3 (0). Then t++ 1 (0) < t + 1 (0) < t 2 (0) < t 3 (0). This is a contradiction to how LFF works. Currently t 3 (0) is the last finishin and t++ 1 (0) is the fastest finishin, and from the equation above if b 1 ave back an aent to b 3 then t + 1 (0) < t 2 (0). This means the transfer back has reduced the time of latest finishin task, thus under RT-LFF b 1 will be not allocated an aent from more than one cluster compared to LFF. This implies a task in RT-LFF cannot have more than one aent under the allocation of LFF for each time step. The worst case is when B 1 tasks have one aent fewer than LFF s allocation with the last task over allocated. T T (a, b ) Alorithm 2 Real-Time Latest Finishin First 1: Use LF F for initial allocation 2: while t < t s do 3: for b i B do 4: if b i b j where t i (0) < t+ j (0) with armin b i T T (b i, b j ) then 5: Reassin a k N i (t) to b j with armin 6: end if a k k j 7: end for 8: end while

4 VI. MODELING COST OF TASK COMPLETION FOR EXTINGUISHING FIRES In this section we focus only on the subproblem of dealin with fires in RoboCup Rescue and present a way of modelin how fire spreads and assinin aents to fires. The RoboCup Rescue simulator is desined for urban search & rescue after an earthquake, which fires to start in buildins. The environment is lare (see Fiure 1) with upward of 100 aents. The full simulator uses heteroeneous aents, but for this work we focus only on the aents that can extinuish fires, i.e. firetrucks. The other aents simulate the oracle and constantly search for new fires, and relay the location and estimated size when a new fire is found. Fires are the most danerous hazard in RoboCup Rescue. While a sinle buildin on fire can be dealt with quickly, if too many buildins inite the fire becomes very difficult to tackle both due to its size and re-initions of buildins. We present two novel contributions. First, fire clusters are modeled as sinle tasks that have a cost which increases as time passes. Second, we present a method to estimate the number of buildins on fire in a cluster when only a few buildins from that cluster have been observed. The RT-LFF alorithm is then shown to out-perform more naïve heuristics. A. Growth of Fire Clusters Each buildin on fire individually has a chance to inite nearby buildins based primarily on distance. This means the rate of rowth,, is proportional to the number of current buildins on fire, x: δx (a) δt = x (b) x = C e t (5) Ten simulations with a sinle fire startin in various locations were used to empirically evaluate to be rouhly This is a first order linear differential equation that can be explicitly solved by separation of variables to yield the well known exponential rowth function iven by (5b). A constant C is introduced by interation to satisfy the initial conditions. Eq. (5a) can be modified to incorporate fire aents extinuishin effect on a fire. If there are n fire aents assined to a fire cluster that each extinuish at a rate w, then the rate of rowth of a fire of a fire cluster will be reduced by n w: δx δt = x n w (6) The constant w was also empirically calculated to be about by runnin ten simulations with a sinle fire aent and trackin the total number of fires extinuished after 100 cycles. Matters are further complicated since the intensity of the fire has an effect on the amount of time it takes to extinuish, thus, w is biased hiher than the real value. Aents can extinuish small fires much more quickly than larer fires, which often causes the aent to repeatedly put out small fires as they are reinited from nearby larer fires. Nevertheless, w ives a reasonable estimate for the aent s capabilities as shown in Section VII. Since both n and w are constants, this still is a linear differential equation which simplifies to a slihtly modified exponential rowth function shown in (7). This satisfies all the conditions in Theorem 2 except the zero travel-time assumption. x = n w + C e t (7) Fi. 1. Pink and blue buildins averae position determine the two points on the red fire cluster circle. B. Estimatin the Size of a Cluster Clusterin is commonly done to roup similar objects into a sinle abstract object to reduce the complexity and to enable a macro-level analysis. Runnin a probabilistic model for every sinle buildin to predict the fire spread would require a lare amount of computation and have a low probability for every possible state. Furthermore, the model will chane based on which fires aents are assined to extinuish, so ideally the model should be recomputed after every assinment. This method will not scale and is too complex for a scenario in which a lare amount of information is already missin. For that reason fires are abstracted into clusters and the macro-level behaviors of these clusters are analyzed instead. The lack of full information in RoboCup Rescue makes clusterin difficult and requires some assumptions to be made. If a lare number of aents were available to circle around the fire cluster and monitor its rowth, then direct clusterin usin (6) with current known information would be a ood estimate. However, normally aents are not able to dedicate this much time to information atherin, so it is necessary to come up with a way of estimatin the size of a fire cluster while only bein able to see a few buildins. To do this we need to make one assumption: on averae fire propaates equally in all directions. From this assumption, a circle is a ood estimate of the area on fire. This circle is identified by two points on its ede in the followin manner: First, the most recently seen buildin and the 9 closest burnin buildins (or all known nearby burnin buildins if fewer than 9 are on fire) are included in a set. Then k-means clusterin with k = 2 is run to find two subsets. The averae position of each subset is then used as the two points to fit the circle upon, shown in Fiure 1. The method for findin the radius can be extended from the exponential model iven in Section VI-A. If x is the number of buildins on fire, then we can estimate π r 2 = A x, where r is the circle radius and A is the averae buildin area (estimated over all buildins on the map). This can be rewritten

5 as: x = π r2 A which can then be substituted into (5a) yieldin: δr δt = 2 r, and r(t) = D e 2 t (8) To overestimate the radius of the circle, we initially assume the fire started at the beinnin of the simulation. For example if a fire is found at time t = 50, we would assume this fire started as a sinle burnin buildin, D = 1, at time t = 0. Thus the radius would simply be r(50). With a radius and two points on a circle, there are two possible circles to choose from. The circle with the hihest ratio of burnin buildins to total buildins is the one chosen as the fire cluster. Since this radius is the worst case estimate, we should incorporate current information to refine the estimate. For all buildins in the circle, we check their status at the last time viewed. If a buildin has never been seen, we nelect it. If there is a conflict between past information and the assumption that the fire started at time t = 0, we assume the center of the circle is still the same and recompute D for a radius that satisfies the information. This new initial condition ives a smaller radius to be fit on the two circle points and this time we choose the circle whose center was inside the old overestimated circle. This process is repeated until there is no conflictin information. C. Assinment to Fire Clusters The oal of assinin aents to clusters is to extinuish fires as quickly as possible. Eq. (7) can be solved for t when x = 0 as shown in (9). If the constant C is nonneative, then the fire is rowin faster than the n aents can extinuish it. Thus the fire will row indefinitely and t(0) is set to. When C is neative, (9) will be computable and will ive the time when the fire will be fully extinuished. t(0) = (ln n w )/ (9) C Alorithm 3 shows the implementation of the RT-LFF alorithm for the case when (x) = x in RoboCup Rescue. Normally in RoboCup Rescue no fire clusters are known initially, so when the first fire cluster is found all aents are assined to that cluster. When a new fire cluster is found, multiple aents may need to be transfered from the old clusters. This is why any time an aent was transfered, we should run the alorithm aain until no more useful transfers exist. The order in which the fire clusters are labeled is important, since every cluster first attempts to transfer an aent to b 1. When a new fire is discovered, it will initially have no aents assined and will need to receive aents from other clusters. Since this fire cluster is in the reatest need of aents, it is assined b 1. All other clusters are then reordered based on their distance from cluster b 1, the newly discovered fire. Thus b 2 will be the closest fire cluster to b 1, b 3 the second closest and so forth. This helps reduce the travel time T T (i, j) in order to more efficiently use fire trucks. VII. RESULTS In this section, we show the validity of our exponential model and RT-LFF by empirical evaluation in RoboCup Rescue. First, we show how the exponential model accurately fits the real fire rowth data. Then the model is used to estimate Alorithm 3 RT-LFF for (x) = x Require: C i, C j, n i, n j, w,, t {C i is the interation constant and n i (t) is the number of aents on task b i, C j and n j correspond to similar thins on task b j, w is the task completion rate of aents, is the rowth rate of fires and t is the current time.} 1: chane true 2: while chane do 3: chane false 4: for i = 1 to k do 5: t i (ni(t) 1) w (0) (ln C i )/ {Compute the time to extinuish fire i with one less aent than it currently is assined.} 6: for j = 1 to k do 7: ˆx j nj(t) w + C j e (t+t T (i,j)) {Before the aent from fire cluster i arrives to cluster j, compute the effect of the aents.} 8: Ĉ j ( ˆx j (nj(t)+1) w )/(e (t+t T (i,j)) ) {Once the aent from cluster i arrives, we need to recompute C j.} 9: t + (nj(t)+1) w j (0) (ln )/ {Finally compute when fire j is fully Ĉj extinuished.} 10: if t i (0) < t+ j (0) then 11: Transfer an aent from i to j 12: n i n i 1 13: n j n j : chane true 15: end if 16: end for 17: end for 18: end while the time a fire is extinuished and compared aainst the real time it took for aents to extinuish the fire. Finally RT-LFF is compared aainst a random assinment method and a closest distance reedy heuristic. A. Model Fittin A sinle fire was tracked over 5 simulations on two maps (Virtual City and Berlin) and the actual number of fires is compared aainst the exponential function best fit in Fiure 2. The exponential function slihtly under estimates the fire cluster between t = 40 and t = 60 but is otherwise fairly accurate. Fi. 2. Averae fire spread estimated by best exponential fit to data. The estimate of fire cluster size described in Section VI-B

6 was tested by assinin a static amount of fire trucks to extinuish the fire cluster and comparin the estimated and real extinuish time. The estimated size of the cluster is recomputed every time step and increases in accuracy as a larer number of buildins in the cluster is observed. For worst case analysis the real extinuish time is compared aainst the estimated extinuish time when the fire cluster is first discovered. A historam of 54 fires extinuished with the normalized error is shown in Fiure 3. The error in estimation had a mean of and standard deviation of The distribution is close to Gaussian and fairly unbiased at overestimatin or underestimatin. Some error may be due to imprecise empirically derived constants in the modelin equation or a lack of incorporatin the intensity in the model. Fi. 4. Averae finish time (no fires left) per map for each confiuration. First alorithm, which attempts to minimize the time when the last task finishes. We documented properties exhibited by dynamic tasks and went on to present a real-time solution, which approximates LFF for partially observable spaces. The alorithms currently work only for homoeneous aents with identical abilities. We hope to expand the current model to include other aent types. The choice of metric for evaluatin an alorithm for dynamic cost functions is also an open question, especially when no solution exists. Completin fewer tasks at the expense of inorin others miht be beneficial more than keepin the current total cost low. REFERENCES Fi. 3. Relative error of extinuish time when a cluster is first discovered. B. Performance Comparison of Different Allocation Strateies Each map was run with 5 simulations and the time the last fire was extinuished is reported in Fiure 4. In each case, the oracle is simulated by a fixed number of aents randomly searchin the environment for new tasks and broadcastin the location of the tasks when they are discovered. The fixed number of aents searchin the environment do not interact with any of the tasks directly. RT-LFF method does much better in Virtual City (8.3% better than the closest heuristic approach) than in Berlin (4.0% better). This is probably due to Virtual City bein an artificial map with similar buildin sizes and a compact buildin arranement. Berlin has open spaces where there is a river or park that throws off the accuracy of the circle estimatin method. There is also a larer discrepancy in buildin size in Berlin which can aain effect the circle estimate. If a lare buildin is chosen for one of the two clusters to estimate the two points on the circle, the buildins in this cluster will be more spread out and would ive an inaccurate point estimate since the center of the lare buildin is far away from the others. The random method performs poorly on both maps, due to lack of cooperation. VIII. CONCLUSIONS AND FUTURE WORK This paper addresses task allocation with multiple aents, where each task has a cost that chanes over time. This adds a substantial amount of complexity and requires more coordination between aents. We limited the investiation of dynamic cost tasks to a eneral family of functions in order to reduce the complexity. We presented the Latest Finishin [1] S. D. Ramchurn, M. Polukarov, A. Farinelli, C. Truon, and N. R. Jennins, Coalition formation with spatial and temporal constraints, in Proc. Int l Conf. on Autonomous Aents and Multi-Aent Systems, 2010, pp [2] P. Scerri, A. Farinelli, S. Okamoto, and M. Tambe, Allocatin tasks in extreme teams, in Proc. Int l Conf. on Autonomous Aents and Multi- Aent Systems, 2005, pp [3] P. R. Ferreira, Jr., F. dos Santos, A. L. C. Bazzan, D. Epstein, and S. J. Waskow, RoboCup Rescue as multiaent task allocation amon teams: experiments with task interdependencies, Journal of Autonomous Aents and Multi-Aent Systems, vol. 20, no. 3, pp , [4] Y. Zhan and L. E. Parker, Task allocation with executable coalitions in multirobot tasks, in Proc. IEEE Int l Conf. on Robotics and Automation, [5] M. Nanjanath, A. Erlandson, S. Andrist, A. Raipindi, A. Mohammed, A. Sharma, and M. Gini, Decision and coordination strateies for robocup rescue aents, in Proc. SIMPAR, 2010, pp [6] A. Chapman, R. A. Micillo, R. Kota, and N. Jennins, Decentralised dynamic task allocation usin overlappin potential ames, The Computer Journal, [7] T. Sandholm, K. Larson, M. Andersson, O. Shehory, and F. Tohmé, Generation with worst case uarantees, Artificial Intellience, vol. 111, no. 1 2, pp , [8] V. D. Dan and N. R. Jennins, Generatin coalition structures with finite bound from the optimal uarantees, in Proc. Int l Conf. on Autonomous Aents and Multi-Aent Systems, 2004, pp [9] X. Zhen and S. Koeni, Reaction functions for task allocation to cooperative aents, in Proc. Int l Conf. on Autonomous Aents and Multi-Aent Systems, 2008, pp [10] H. Kitano and S. Tadokoro, RoboCup rescue: A rand challene for multiaent and intellient systems, AI Maazine, vol. 22, no. 1, pp , [11] Álvaro Monares, S. F. Ochoa, J. A. Pino, V. Herskovic, J. Rodriuez- Covili, and A. Neyem, Mobile computin in urban emerency situations: Improvin the support to firefihters in the field, Expert Systems with Applications, vol. 38, no. 2, pp , 2011.

PHY 133 Lab 1 - The Pendulum

PHY 133 Lab 1 - The Pendulum 3/20/2017 PHY 133 Lab 1 The Pendulum [Stony Brook Physics Laboratory Manuals] Stony Brook Physics Laboratory Manuals PHY 133 Lab 1 - The Pendulum The purpose of this lab is to measure the period of a simple

More information

Simple Optimization (SOPT) for Nonlinear Constrained Optimization Problem

Simple Optimization (SOPT) for Nonlinear Constrained Optimization Problem (ISSN 4-6) Journal of Science & Enineerin Education (ISSN 4-6) Vol.,, Pae-3-39, Year-7 Simple Optimization (SOPT) for Nonlinear Constrained Optimization Vivek Kumar Chouhan *, Joji Thomas **, S. S. Mahapatra

More information

5. Network Analysis. 5.1 Introduction

5. Network Analysis. 5.1 Introduction 5. Network Analysis 5.1 Introduction With the continued rowth o this country as it enters the next century comes the inevitable increase in the number o vehicles tryin to use the already overtaxed transportation

More information

LP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time

LP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time LP Roundin and Combinatorial Alorithms for Minimizin Active and Busy Time Jessica Chan, Samir Khuller, and Koyel Mukherjee University of Maryland, Collee Park {jschan,samir,koyelm}@cs.umd.edu Abstract.

More information

A Constant Complexity Fair Scheduler with O(log N) Delay Guarantee

A Constant Complexity Fair Scheduler with O(log N) Delay Guarantee A Constant Complexity Fair Scheduler with O(lo N) Delay Guarantee Den Pan and Yuanyuan Yan 2 Deptment of Computer Science State University of New York at Stony Brook, Stony Brook, NY 79 denpan@cs.sunysb.edu

More information

2.2 Differentiation and Integration of Vector-Valued Functions

2.2 Differentiation and Integration of Vector-Valued Functions .. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin

More information

Strong Interference and Spectrum Warfare

Strong Interference and Spectrum Warfare Stron Interference and Spectrum Warfare Otilia opescu and Christopher Rose WILAB Ruters University 73 Brett Rd., iscataway, J 8854-86 Email: {otilia,crose}@winlab.ruters.edu Dimitrie C. opescu Department

More information

Wire antenna model of the vertical grounding electrode

Wire antenna model of the vertical grounding electrode Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna

More information

Lecture 5 Processing microarray data

Lecture 5 Processing microarray data Lecture 5 Processin microarray data (1)Transform the data into a scale suitable for analysis ()Remove the effects of systematic and obfuscatin sources of variation (3)Identify discrepant observations Preprocessin

More information

Efficient method for obtaining parameters of stable pulse in grating compensated dispersion-managed communication systems

Efficient method for obtaining parameters of stable pulse in grating compensated dispersion-managed communication systems 3 Conference on Information Sciences and Systems, The Johns Hopkins University, March 12 14, 3 Efficient method for obtainin parameters of stable pulse in ratin compensated dispersion-manaed communication

More information

Causal Bayesian Networks

Causal Bayesian Networks Causal Bayesian Networks () Ste7 (2) (3) Kss Fus3 Ste2 () Fiure : Simple Example While Bayesian networks should typically be viewed as acausal, it is possible to impose a causal interpretation on these

More information

A Probabilistic Analysis of Propositional STRIPS. Planning. Tom Bylander. Division of Mathematics, Computer Science, and Statistics

A Probabilistic Analysis of Propositional STRIPS. Planning. Tom Bylander. Division of Mathematics, Computer Science, and Statistics A Probabilistic Analysis of Propositional STRIPS Plannin Tom Bylander Division of Mathematics, Computer Science, and Statistics The University of Texas at San Antonio San Antonio, Texas 78249 USA bylander@riner.cs.utsa.edu

More information

Phase Diagrams: construction and comparative statics

Phase Diagrams: construction and comparative statics 1 / 11 Phase Diarams: construction and comparative statics November 13, 215 Alecos Papadopoulos PhD Candidate Department of Economics, Athens University of Economics and Business papadopalex@aueb.r, https://alecospapadopoulos.wordpress.com

More information

Learning for Decentralized Control of Multiagent Systems in Large, Partially-Observable Stochastic Environments

Learning for Decentralized Control of Multiagent Systems in Large, Partially-Observable Stochastic Environments Learnin for Decentralized Control of Multiaent Systems in Lare, Partially-Observable Stochastic Environments Miao Liu LIDS, MIT Cambride, MA miaoliu@mit.edu Christopher Amato Dept.of CS, U. of New Hampshire

More information

Physics 20 Homework 1 SIMS 2016

Physics 20 Homework 1 SIMS 2016 Physics 20 Homework 1 SIMS 2016 Due: Wednesday, Auust 17 th Problem 1 The idea of this problem is to et some practice in approachin a situation where you miht not initially know how to proceed, and need

More information

Optimization of Mechanical Design Problems Using Improved Differential Evolution Algorithm

Optimization of Mechanical Design Problems Using Improved Differential Evolution Algorithm International Journal of Recent Trends in Enineerin Vol. No. 5 May 009 Optimization of Mechanical Desin Problems Usin Improved Differential Evolution Alorithm Millie Pant Radha Thanaraj and V. P. Sinh

More information

Habits and Multiple Equilibria

Habits and Multiple Equilibria Habits and Multiple Equilibria Lorenz Kuen and Eveny Yakovlev February 5, 2018 Abstract Abstract TBA. A Introduction Brief introduction TBA B A Structural Model of Taste Chanes Several structural models

More information

Generalized Distance Metric as a Robust Similarity Measure for Mobile Object Trajectories

Generalized Distance Metric as a Robust Similarity Measure for Mobile Object Trajectories Generalized Distance Metric as a Robust Similarity Measure for Mobile Object rajectories Garima Pathak, Sanjay Madria Department of Computer Science University Of Missouri-Rolla Missouri-6541, USA {madrias}@umr.edu

More information

Linearized optimal power flow

Linearized optimal power flow Linearized optimal power flow. Some introductory comments The advantae of the economic dispatch formulation to obtain minimum cost allocation of demand to the eneration units is that it is computationally

More information

Ning Wu Institute for Traffic Engineering Ruhr University Bochum, Germany Tel: ; Fax: ;

Ning Wu Institute for Traffic Engineering Ruhr University Bochum, Germany Tel: ; Fax: ; MODELLING THE IMPACT OF SIDE-STREET TRAFFIC VOLUME ON MAJOR- STREET GREEN TIME AT ISOLATED SEMI-ACTUATED INTERSECTIONS FOR SIGNAL COORDINATION DECISIONS Donmei Lin, Correspondin Author Center for Advanced

More information

1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance

1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance CHAPTER 7 PROJECTILES 7 No Air Resistance We suppose that a particle is projected from a point O at the oriin of a coordinate system, the y-axis bein vertical and the x-axis directed alon the round The

More information

IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1

IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1 Intervention in Gene Reulatory Networks via a Stationary Mean-First-Passae-Time Control Policy Golnaz Vahedi, Student Member, IEEE, Babak Faryabi, Student

More information

Problem Set 2 Solutions

Problem Set 2 Solutions UNIVERSITY OF ALABAMA Department of Physics and Astronomy PH 125 / LeClair Sprin 2009 Problem Set 2 Solutions The followin three problems are due 20 January 2009 at the beinnin of class. 1. (H,R,&W 4.39)

More information

A Rate-Splitting Strategy for Max-Min Fair Multigroup Multicasting

A Rate-Splitting Strategy for Max-Min Fair Multigroup Multicasting A Rate-Splittin Stratey for Max-Min Fair Multiroup Multicastin Hamdi Joudeh and Bruno Clercx Department of Electrical and Electronic Enineerin, Imperial Collee London, United Kindom School of Electrical

More information

Altitude measurement for model rocketry

Altitude measurement for model rocketry Altitude measurement for model rocketry David A. Cauhey Sibley School of Mechanical Aerospace Enineerin, Cornell University, Ithaca, New York 14853 I. INTRODUCTION In his book, Rocket Boys, 1 Homer Hickam

More information

Experiment 3 The Simple Pendulum

Experiment 3 The Simple Pendulum PHY191 Fall003 Experiment 3: The Simple Pendulum 10/7/004 Pae 1 Suested Readin for this lab Experiment 3 The Simple Pendulum Read Taylor chapter 5. (You can skip section 5.6.IV if you aren't comfortable

More information

V DD. M 1 M 2 V i2. V o2 R 1 R 2 C C

V DD. M 1 M 2 V i2. V o2 R 1 R 2 C C UNVERSTY OF CALFORNA Collee of Enineerin Department of Electrical Enineerin and Computer Sciences E. Alon Homework #3 Solutions EECS 40 P. Nuzzo Use the EECS40 90nm CMOS process in all home works and projects

More information

Matrix multiplication: a group-theoretic approach

Matrix multiplication: a group-theoretic approach CSG399: Gems of Theoretical Computer Science. Lec. 21-23. Mar. 27-Apr. 3, 2009. Instructor: Emanuele Viola Scribe: Ravi Sundaram Matrix multiplication: a roup-theoretic approach Given two n n matrices

More information

THE SIGNAL ESTIMATOR LIMIT SETTING METHOD

THE SIGNAL ESTIMATOR LIMIT SETTING METHOD ' THE SIGNAL ESTIMATOR LIMIT SETTING METHOD Shan Jin, Peter McNamara Department of Physics, University of Wisconsin Madison, Madison, WI 53706 Abstract A new method of backround subtraction is presented

More information

(a) Find the function that describes the fraction of light bulbs failing by time t e (0.1)x dx = [ e (0.1)x ] t 0 = 1 e (0.1)t.

(a) Find the function that describes the fraction of light bulbs failing by time t e (0.1)x dx = [ e (0.1)x ] t 0 = 1 e (0.1)t. 1 M 13-Lecture March 8, 216 Contents: 1) Differential Equations 2) Unlimited Population Growth 3) Terminal velocity and stea states Voluntary Quiz: The probability density function of a liht bulb failin

More information

MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM. Apurba Sakti EGEE 520, Mathematical Modeling of EGEE systems Spring 2007

MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM. Apurba Sakti EGEE 520, Mathematical Modeling of EGEE systems Spring 2007 MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM Apurba Sakti EGEE 520, Mathematical Modelin of EGEE systems Sprin 2007 Table of Contents Abstract Introduction Governin Equations Enery balance Momentum

More information

Generalized Least-Squares Regressions V: Multiple Variables

Generalized Least-Squares Regressions V: Multiple Variables City University of New York (CUNY) CUNY Academic Works Publications Research Kinsborouh Community Collee -05 Generalized Least-Squares Reressions V: Multiple Variables Nataniel Greene CUNY Kinsborouh Community

More information

Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013

Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013 Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR Andrew Hall 11/7/2013 Outline RBWR Motivation and Desin Why use Serpent Cross Sections? Modelin the RBWR Axial Discontinuity

More information

LP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time

LP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time LP Roundin and Combinatorial Alorithms for Minimizin Active and Busy Time Jessica Chan University of Maryland Collee Park, MD, USA jschan@umiacs.umd.edu Samir Khuller University of Maryland Collee Park,

More information

arxiv: v1 [math.oc] 26 Sep 2017

arxiv: v1 [math.oc] 26 Sep 2017 Data-Driven Analysis of Mass-Action Kinetics Camilo Garcia-Tenorio, Duvan Tellez-Castro, Eduardo Mojica-Nava, and Jore Sofrony arxiv:1709.09273v1 [math.oc] 26 Sep 2017 Abstract The physical interconnection

More information

Using Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems

Using Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems Usin Quasi-Newton Methos to Fin Optimal Solutions to Problematic Kriin Systems Steven Lyster Centre for Computational Geostatistics Department of Civil & Environmental Enineerin University of Alberta Solvin

More information

Extended Target Poisson Multi-Bernoulli Filter

Extended Target Poisson Multi-Bernoulli Filter 1 Extended Taret Poisson Multi-Bernoulli Filter Yuxuan Xia, Karl Granström, Member, IEEE, Lennart Svensson, Senior Member, IEEE, and Maryam Fatemi arxiv:1801.01353v1 [eess.sp] 4 Jan 2018 Abstract In this

More information

First- and Second Order Phase Transitions in the Holstein- Hubbard Model

First- and Second Order Phase Transitions in the Holstein- Hubbard Model Europhysics Letters PREPRINT First- and Second Order Phase Transitions in the Holstein- Hubbard Model W. Koller 1, D. Meyer 1,Y.Ōno 2 and A. C. Hewson 1 1 Department of Mathematics, Imperial Collee, London

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.01x Fall Term 2001 EXAM 1 SOLUTIONS

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.01x Fall Term 2001 EXAM 1 SOLUTIONS MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01x Fall Term 2001 EXAM 1 SOLUTIONS Problem 1: We define a vertical coordinate system with positive upwards. The only forces actin

More information

Robust Semiparametric Optimal Testing Procedure for Multiple Normal Means

Robust Semiparametric Optimal Testing Procedure for Multiple Normal Means Veterinary Dianostic and Production Animal Medicine Publications Veterinary Dianostic and Production Animal Medicine 01 Robust Semiparametric Optimal Testin Procedure for Multiple Normal Means Pen Liu

More information

Disclaimer: This lab write-up is not

Disclaimer: This lab write-up is not Disclaimer: This lab write-up is not to be copied, in whole or in part, unless a proper reference is made as to the source. (It is stronly recommended that you use this document only to enerate ideas,

More information

Round-off Error Free Fixed-Point Design of Polynomial FIR Predictors

Round-off Error Free Fixed-Point Design of Polynomial FIR Predictors Round-off Error Free Fixed-Point Desin of Polynomial FIR Predictors Jarno. A. Tanskanen and Vassil S. Dimitrov Institute of Intellient Power Electronics Department of Electrical and Communications Enineerin

More information

4.3. Solving Friction Problems. Static Friction Problems. Tutorial 1 Static Friction Acting on Several Objects. Sample Problem 1.

4.3. Solving Friction Problems. Static Friction Problems. Tutorial 1 Static Friction Acting on Several Objects. Sample Problem 1. Solvin Friction Problems Sometimes friction is desirable and we want to increase the coefficient of friction to help keep objects at rest. For example, a runnin shoe is typically desined to have a lare

More information

PhysicsAndMathsTutor.com

PhysicsAndMathsTutor.com . A particle of mass m is projected vertically upwards, at time t =, with speed. The particle is mv subject to air resistance of manitude, where v is the speed of the particle at time t and is a positive

More information

Motion in Two Dimensions Sections Covered in the Text: Chapters 6 & 7, except 7.5 & 7.6

Motion in Two Dimensions Sections Covered in the Text: Chapters 6 & 7, except 7.5 & 7.6 Motion in Two Dimensions Sections Covered in the Tet: Chapters 6 & 7, ecept 7.5 & 7.6 It is time to etend the definitions we developed in Note 03 to describe motion in 2D space. In doin so we shall find

More information

Energizing Math with Engineering Applications

Energizing Math with Engineering Applications Enerizin Math with Enineerin Applications Understandin the Math behind Launchin a Straw-Rocket throuh the use of Simulations. Activity created by Ira Rosenthal (rosenthi@palmbeachstate.edu) as part of

More information

Emotional Optimized Design of Electro-hydraulic Actuators

Emotional Optimized Design of Electro-hydraulic Actuators Sensors & Transducers, Vol. 77, Issue 8, Auust, pp. 9-9 Sensors & Transducers by IFSA Publishin, S. L. http://www.sensorsportal.com Emotional Optimized Desin of Electro-hydraulic Actuators Shi Boqian,

More information

11 Free vibrations: one degree of freedom

11 Free vibrations: one degree of freedom 11 Free vibrations: one deree of freedom 11.1 A uniform riid disk of radius r and mass m rolls without slippin inside a circular track of radius R, as shown in the fiure. The centroidal moment of inertia

More information

Prince Sultan University Physics Department First Semester 2012 /2013. PHY 105 First Major Exam Allowed Time: 60 min

Prince Sultan University Physics Department First Semester 2012 /2013. PHY 105 First Major Exam Allowed Time: 60 min Prince Sultan University Physics Department First Semester 01 /01 PHY 105 First Major Exam Allowed Time: 60 min Student Name: 1. Write your name in the specified space NOW.. Any paper without name will

More information

Bidirectional Clustering of Weights for Finding Succinct Multivariate Polynomials

Bidirectional Clustering of Weights for Finding Succinct Multivariate Polynomials IJCSNS International Journal of Computer Science and Network Security, VOL.8 No.5, May 28 85 Bidirectional Clusterin of Weihts for Findin Succinct Multivariate Polynomials Yusuke Tanahashi and Ryohei Nakano

More information

2.3. PBL Equations for Mean Flow and Their Applications

2.3. PBL Equations for Mean Flow and Their Applications .3. PBL Equations for Mean Flow and Their Applications Read Holton Section 5.3!.3.1. The PBL Momentum Equations We have derived the Reynolds averaed equations in the previous section, and they describe

More information

Statistical Mechanics: A Possible Model for Market-based Electric Power Control

Statistical Mechanics: A Possible Model for Market-based Electric Power Control Proceedins of the 37th Hawaii International Conference on System Sciences - 4 Statistical Mechanics: A Possible Model for Market-based Electric Power Control David P Chassin Pacific orthwest ational Laboratory

More information

A Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere

A Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere A Mathematical Model for the Fire-extinuishin Rocket Fliht in a Turbulent Atmosphere CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tiroviste, Km 8 ROMANIA crismihailescu@yahoo.com http://www.elmec.ro

More information

RATE OPTIMIZATION FOR MASSIVE MIMO RELAY NETWORKS: A MINORIZATION-MAXIMIZATION APPROACH

RATE OPTIMIZATION FOR MASSIVE MIMO RELAY NETWORKS: A MINORIZATION-MAXIMIZATION APPROACH RATE OPTIMIZATION FOR MASSIVE MIMO RELAY NETWORKS: A MINORIZATION-MAXIMIZATION APPROACH Mohammad Mahdi Nahsh, Mojtaba Soltanalian, Petre Stoica +, Maryam Masjedi, and Björn Ottersten Department of Electrical

More information

Slopes of Lefschetz fibrations and separating vanishing cycles

Slopes of Lefschetz fibrations and separating vanishing cycles Slopes of Lefschetz fibrations and separatin vanishin cycles Yusuf Z. Gürtaş Abstract. In this article we find upper and lower bounds for the slope of enus hyperelliptic Lefschetz fibrations. We demonstrate

More information

5 Shallow water Q-G theory.

5 Shallow water Q-G theory. 5 Shallow water Q-G theory. So far we have discussed the fact that lare scale motions in the extra-tropical atmosphere are close to eostrophic balance i.e. the Rossby number is small. We have examined

More information

Calculating Well Deliverability in Gas Condensate Reservoirs

Calculating Well Deliverability in Gas Condensate Reservoirs 1 14 Calculatin Well Deliverability in Gas Condensate Reservoirs ROBERT MOTT AEA Technoloy, Winfrith, Dorchester, Dorset DT2 8DH, United Kindom Abstract Well deliverability in most as condensate reservoirs

More information

Assessment of the MCNP-ACAB code system for burnup credit analyses. N. García-Herranz 1, O. Cabellos 1, J. Sanz 2

Assessment of the MCNP-ACAB code system for burnup credit analyses. N. García-Herranz 1, O. Cabellos 1, J. Sanz 2 Assessment of the MCNP-ACAB code system for burnup credit analyses N. García-Herranz, O. Cabellos, J. Sanz 2 Departamento de Ineniería Nuclear, Universidad Politécnica de Madrid 2 Departamento de Ineniería

More information

Physics 18 Spring 2011 Homework 2 - Solutions Wednesday January 26, 2011

Physics 18 Spring 2011 Homework 2 - Solutions Wednesday January 26, 2011 Physics 18 Sprin 011 Homework - s Wednesday January 6, 011 Make sure your name is on your homework, and please box your final answer. Because we will be ivin partial credit, be sure to attempt all the

More information

Assessment of the MCNP-ACAB code system for burnup credit analyses

Assessment of the MCNP-ACAB code system for burnup credit analyses Assessment of the MCNP-ACAB code system for burnup credit analyses N. García-Herranz, O. Cabellos, J. Sanz UPM - UNED International Workshop on Advances in Applications of Burnup Credit for Spent Fuel

More information

ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor

ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor Ade afar Abdullah Electrical Enineerin Department, Faculty of Technoloy and Vocational Education Indonesia University of Education

More information

From Social Trust Assisted Reciprocity (STAR) to Utility-Optimal Mobile Crowdsensing

From Social Trust Assisted Reciprocity (STAR) to Utility-Optimal Mobile Crowdsensing From ocial Trust Assisted eciprocity (TA) to Utility-Optimal Mobile Crowdsensin Xiaowen Gon, Xu Chen, Junshan Zhan, H. Vincent Poor chool of Electrical, Computer and Enery Enineerin Arizona tate University,

More information

Solution to the take home exam for ECON 3150/4150

Solution to the take home exam for ECON 3150/4150 Solution to the tae home exam for ECO 350/450 Jia Zhiyan and Jo Thori Lind April 2004 General comments Most of the copies we ot were quite ood, and it seems most of you have done a real effort on the problem

More information

A Performance Comparison Study with Information Criteria for MaxEnt Distributions

A Performance Comparison Study with Information Criteria for MaxEnt Distributions A Performance Comparison Study with nformation Criteria for MaxEnt Distributions Ozer OZDEMR and Aslı KAYA Abstract n statistical modelin, the beinnin problem that has to be solved is the parameter estimation

More information

Asymptotic Behavior of a t Test Robust to Cluster Heterogeneity

Asymptotic Behavior of a t Test Robust to Cluster Heterogeneity Asymptotic Behavior of a t est Robust to Cluster Heteroeneity Andrew V. Carter Department of Statistics University of California, Santa Barbara Kevin. Schnepel and Doulas G. Steierwald Department of Economics

More information

Termination criteria in the Moore-Skelboe Algorithm for Global Optimization by Interval Arithmetic

Termination criteria in the Moore-Skelboe Algorithm for Global Optimization by Interval Arithmetic To appear in series Nonconvex Global Optimization and Its Applications Kluwer Academic Publishers Termination criteria in the Moore-Skelboe Alorithm for Global Optimization by Interval Arithmetic M.H.

More information

The 10 th international Energy Conference (IEC 2014)

The 10 th international Energy Conference (IEC 2014) The 10 th international Enery Conference (IEC 2014) Wind Turbine Interated Control durin Full Load Operation 1. Hamed Habibi and 2. Ahil Yousefi-Koma 1. MSc, Centre of Advanced Systems and Technoloies

More information

Stochastic simulations of genetic switch systems

Stochastic simulations of genetic switch systems Stochastic simulations of enetic switch systems Adiel Loiner, 1 Azi Lipshtat, 2 Nathalie Q. Balaban, 1 and Ofer Biham 1 1 Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 2 Department

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Journal of Inequalities in Pure and Applied Mathematics L HOSPITAL TYPE RULES FOR OSCILLATION, WITH APPLICATIONS IOSIF PINELIS Department of Mathematical Sciences, Michian Technoloical University, Houhton,

More information

Security Constrained Optimal Power Flow

Security Constrained Optimal Power Flow Security Constrained Optimal Power Flow 1. Introduction and notation Fiure 1 below compares te optimal power flow (OPF wit te security-constrained optimal power flow (SCOPF. Fi. 1 Some comments about tese

More information

This is a repository copy of The effects of gas-phase and in-depth radiation absorption on ignition and steady burning rate of PMMA.

This is a repository copy of The effects of gas-phase and in-depth radiation absorption on ignition and steady burning rate of PMMA. This is a repository copy of The effects of as-phase and in-depth radiation absorption on inition and steady burnin rate of PMMA. White ose esearch Online UL for this paper: http://eprints.whiterose.ac.uk/9515/

More information

Path Loss Prediction in Urban Environment Using Learning Machines and Dimensionality Reduction Techniques

Path Loss Prediction in Urban Environment Using Learning Machines and Dimensionality Reduction Techniques Path Loss Prediction in Urban Environment Usin Learnin Machines and Dimensionality Reduction Techniques Mauro Piacentini Francesco Rinaldi Technical Report n. 11, 2009 Path Loss Prediction in Urban Environment

More information

Finite Lifetimes, Patentsʹ Length and Breadth, and Growth

Finite Lifetimes, Patentsʹ Length and Breadth, and Growth Auburn University Department of Economics Workin Paper Series Finite Lifetimes, Patentsʹ Lenth and Breadth, and Growth Bharat Diwakar and Gilad Sorek Auburn University AUWP 06 08 This paper can be downloaded

More information

11. Excitons in Nanowires and Nanotubes

11. Excitons in Nanowires and Nanotubes Excitons in Nanowires and Nanotubes We have seen that confinement in quantum wells leads to enhanced excitonic effects in the optical response of semiconductors The bindin enery of the stronest bound excitons

More information

Convergence of DFT eigenvalues with cell volume and vacuum level

Convergence of DFT eigenvalues with cell volume and vacuum level Converence of DFT eienvalues with cell volume and vacuum level Sohrab Ismail-Beii October 4, 2013 Computin work functions or absolute DFT eienvalues (e.. ionization potentials) requires some care. Obviously,

More information

Graph Entropy, Network Coding and Guessing games

Graph Entropy, Network Coding and Guessing games Graph Entropy, Network Codin and Guessin ames Søren Riis November 25, 2007 Abstract We introduce the (private) entropy of a directed raph (in a new network codin sense) as well as a number of related concepts.

More information

Solution sheet 9. λ < g u = 0, λ = g α 0 : u = αλ. C(u, λ) = gλ gλ = 0

Solution sheet 9. λ < g u = 0, λ = g α 0 : u = αλ. C(u, λ) = gλ gλ = 0 Advanced Finite Elements MA5337 - WS17/18 Solution sheet 9 On this final exercise sheet about variational inequalities we deal with nonlinear complementarity functions as a way to rewrite a variational

More information

Modeling for control of a three degrees-of-freedom Magnetic. Levitation System

Modeling for control of a three degrees-of-freedom Magnetic. Levitation System Modelin for control of a three derees-of-freedom Manetic evitation System Rafael Becerril-Arreola Dept. of Electrical and Computer En. University of Toronto Manfredi Maiore Dept. of Electrical and Computer

More information

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions

Ground Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions PC11 Fundamentals of Physics I Lectures 7 and 8 Motion in Two Dimensions Dr Tay Sen Chuan 1 Ground Rules Switch off your handphone and paer Switch off your laptop computer and keep it No talkin while lecture

More information

INFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW

INFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW HEFAT7 5 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics Sun City, South Africa Paper number: GJ1 INFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW Gylys

More information

On K-Means Cluster Preservation using Quantization Schemes

On K-Means Cluster Preservation using Quantization Schemes On K-Means Cluster Preservation usin Quantization Schemes Deepak S. Turaa Michail Vlachos Olivier Verscheure IBM T.J. Watson Research Center, Hawthorne, Y, USA IBM Zürich Research Laboratory, Switzerland

More information

An Optimization Problem With a Surprisingly Simple Solution

An Optimization Problem With a Surprisingly Simple Solution An Optimization Problem With a Surprisinly Simple Solution D. Drinen, K. G. Kennedy, and W. M. Priestley 1 INTRODUCTION. Suppose you and n of your friends play the followin ame. A random number from the

More information

MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION

MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION = = = MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION Mauro Cettolo and Marcello Federico ITC-irst - Centro per la Ricerca Scientifica e Tecnoloica I-385 Povo, Trento, Italy ABSTRACT Robust acoustic

More information

THE IMPACT OF ELEVATED MASTS ON LIGHTNING INCIDENCE AT THE RADIO-BASE-STATION VICINITIES

THE IMPACT OF ELEVATED MASTS ON LIGHTNING INCIDENCE AT THE RADIO-BASE-STATION VICINITIES X International Symosium on Lihtnin Protection 9 th -13 th November, 9 Curitiba, Brazil THE IMPACT OF ELEVATED MASTS ON LIGHTNING INCIDENCE AT THE RADIO-BASE-STATION VICINITIES Rosilene Nietzsch Dias 1,

More information

Design of Chevron Gusset Plates

Design of Chevron Gusset Plates 017 SEAOC CONENTION PROCEEDINGS Desin of Chevron Gusset Plates Rafael Sali, Director of Seismic Desin Walter P Moore San Francisco, California Leih Arber, Senior Enineer American Institute of Steel Construction

More information

Indeterminacy in discrete-time infinite-horizon models with non linear utility and endogenous labor

Indeterminacy in discrete-time infinite-horizon models with non linear utility and endogenous labor Indeterminacy in discrete-time infinite-horizon models with non linear utility and endoenous labor Kazuo NISHIMURA Institute of Economic Research, Kyoto University, Japan and Alain VENDII CNRS - GREQAM,

More information

REAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION. Scott Rickard, Radu Balan, Justinian Rosca

REAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION. Scott Rickard, Radu Balan, Justinian Rosca REAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEARATION Scott Rickard, Radu Balan, Justinian Rosca Siemens Corporate Research rinceton, NJ scott.rickard,radu.balan,justinian.rosca @scr.siemens.com ABSTRACT

More information

Stoichiometry of the reaction of sodium carbonate with hydrochloric acid

Stoichiometry of the reaction of sodium carbonate with hydrochloric acid Stoichiometry of the reaction of sodium carbonate with hydrochloric acid Purpose: To calculate the theoretical (expected) yield of product in a reaction. To weih the actual (experimental) mass of product

More information

arxiv:cond-mat/ v1 [cond-mat.str-el] 15 Dec 2003

arxiv:cond-mat/ v1 [cond-mat.str-el] 15 Dec 2003 Imperial Collee Preprints PREPRINT arxiv:cond-mat/0312367v1 [cond-mat.str-el] 15 Dec 2003 First- and Second Order Phase Transitions in the Holstein- Hubbard Model W. Koller 1, D. Meyer 1, Y. Ōno2 and A.

More information

OPTIMAL CONTROL OF AN HIV MODEL

OPTIMAL CONTROL OF AN HIV MODEL H. R. ERFANIAN and M. H. NOORI SKANDARI/ TMCS Vol.2 No. (211) 65-658 The ournal of Mathematics and Computer Science Available online at http://www.tmcs.com The ournal of Mathematics and Computer Science

More information

Stochastic learning feedback hybrid automata for dynamic power management in embedded systems

Stochastic learning feedback hybrid automata for dynamic power management in embedded systems Electrical and Computer Enineerin Faculty Publications Electrical & Computer Enineerin 2005 Stochastic learnin feedback hybrid automata for dynamic power manaement in embedded systems T. Erbes Eurecom

More information

Renormalization Group Theory

Renormalization Group Theory Chapter 16 Renormalization Group Theory In the previous chapter a procedure was developed where hiher order 2 n cycles were related to lower order cycles throuh a functional composition and rescalin procedure.

More information

Fluid Modeling of Pollution Proliferation in P2P Networks

Fluid Modeling of Pollution Proliferation in P2P Networks Fluid Modelin of Pollution Proliferation in P2P Networks Rakesh Kumar David D. Yao Amitabha Bachi Keith W. Ross Dan Rubenstein ABSTRACT P2P systems are hihly vulnerable to pollution attacks in which attackers

More information

A matching of. matrix elements and parton showers. J. Andre and T. Sjostrand 1. Department of Theoretical Physics, Lund University, Lund, Sweden

A matching of. matrix elements and parton showers. J. Andre and T. Sjostrand 1. Department of Theoretical Physics, Lund University, Lund, Sweden LU TP 97{18 Auust 1997 A matchin of matrix elements and parton showers J. Andre and T. Sjostrand 1 Department of Theoretical Physics, Lund University, Lund, Sweden Abstract We propose a simple scheme to

More information

arxiv: v2 [math.oc] 11 Jun 2018

arxiv: v2 [math.oc] 11 Jun 2018 : Tiht Automated Converence Guarantees Adrien Taylor * Bryan Van Scoy * 2 Laurent Lessard * 2 3 arxiv:83.673v2 [math.oc] Jun 28 Abstract We present a novel way of eneratin Lyapunov functions for provin

More information

Whither probabilistic security management for real-time operation of power systems?

Whither probabilistic security management for real-time operation of power systems? 2013 IREP Symposium-Bulk Power Systems Dynamics and Control-IX (IREP), Auust 25-30,2013, Rethymnon, Crete, Greece 1 Whither probabilistic security manaement for real-time operation of power systems? Efthymios

More information

Thermo-Mechanical Damage Modeling of Polymer Matrix Composite Structures in Fire

Thermo-Mechanical Damage Modeling of Polymer Matrix Composite Structures in Fire Thermo-Mechanical Damae Modelin of Polymer Matrix Composite Structures in Fire CHANGSONG LUO, and JIM LUA Global Enineerin and Materials, Inc. One Airport Place, Suite One Princeton, NJ 08540 USA ABSTRACT

More information

8. Compartmental Modeling of Two Phase Systems: Reactions in Gas and Liquid

8. Compartmental Modeling of Two Phase Systems: Reactions in Gas and Liquid ChE 505 Chapter 8N Updated 0/0/07 8. Compartmental Modelin of Two Phase Systems: Reactions in Gas and Liquid In many situations, includin the atmosphere, reactions occur in more than one phase e.. as/air)

More information

On the Chromatic Art Gallery Problem

On the Chromatic Art Gallery Problem CCCG 2014, Halifax, Nova Scotia, Auust 11 13, 2014 On the Chromatic Art Gallery Problem Sándor P. Fekete Stephan Friedrichs Michael Hemmer Joseph B. M. Mitchell Christiane Schmidt Abstract For a polyonal

More information

Health Monitoring of a Truss Bridge using Adaptive Identification

Health Monitoring of a Truss Bridge using Adaptive Identification Proceedins of the 27 IEEE Intellient Transportation Systems Conference Seattle, WA, USA, Sept. 3 - Oct. 3, 27 WeC.4 Health Monitorin of a Truss Bride usin Adaptive Identification James W. Fonda, Student

More information