Solving a Huff-like competitive location and design model for profit maximization in the plane

A chain wants to set up a single new facility in a planar market where similar facilities of competitors, and possibly of its own chain, are already present. Fixed demand points split their demand probabilistically over all facilities in the market proportionally with their attraction to each facility, determined by the different perceived qualities of the facilities and the distances to them, through a gravitational or logit type model. Both the location and the quality (design) of the new facility are to be found so as to maximise the profit obtained for the chain. Several types of constraints and costs are considered.     

An automated planning engine for biopharmaceutical production

We introduce an optimization-based production planning tool for the biotechnology industry. The industry’s planning problem is unusually challenging because the entire production process is regulated by multiple external agencies – such as the US Food and Drug Administration – representing countries where the biopharmaceutical is to be sold. The model is structured to precisely capture the constraints imposed by current and projected regulatory approvals of processes and facilities, as well as capturing the outcomes of quality testing and processing options, facility capacities and initial status of work-in-process. The result is a supply chain “Planning Engine” that generates capacity-feasible batch processing schedules for each production facility within the biomanufacturing supply chain and an availability schedule for finished product against a known set of demands and regulations. Developing the formulation based on distinct time grids tailored for each facility, planning problems with more than 27,000 boolean variables, more than 130,000 linear variables and more than 80,000 constraints are automatically formulated and solved within a few hours. The Planning Engine’s development and implementation at Bayer Healthcare’s Berkeley, CA manufacturing site is described.     

Strategic facility location: A review

Facility location decisions are a critical element in strategic planning for a wide range of private and public firms. The ramifications of siting facilities are broadly based and long-lasting, impacting numerous operational and logistical decisions. High costs associated with property acquisition and facility construction make facility location or relocation projects long-term investments. To make such undertakings profitable, firms plan for new facilities to remain in place and in operation for an extended time period. Thus, decision makers must select sites that will not simply perform well according to the current system state, but that will continue to be profitable for the facility's lifetime, even as environmental factors change, populations shift, and market trends evolve. Finding robust facility locations is thus a difficult task, demanding that decision makers account for uncertain future events. The complexity of this problem has limited much of the facility location literature to simplified static and deterministic models. Although a few researchers initiated the study of stochastic and dynamic aspects of facility location many years ago, most of the research dedicated to these issues has been published in recent years. In this review, we report on literature which explicitly addresses the strategic nature of facility location problems by considering either stochastic or dynamic problem characteristics. Dynamic formulations focus on the difficult timing issues involved in locating a facility (or facilities) over an extended horizon. Stochastic formulations attempt to capture the uncertainty in problem input parameters such as forecast demand or distance values. The stochastic literature is divided into two classes: that which explicitly considers the probability distribution of uncertain parameters, and that which captures uncertainty through scenario planning. A wide range of model formulations and solution approaches are discussed, with applications ranging across numerous industries.     

An analysis of p-median location problem: Effects of backup service level and demand assignment policy

Any solution to facility location problems will consider determining the best suitable locations with respect to certain criteria. Among different types of location problems, involving emergency service system (ESSs) are one of the most widely studied in the literature, and solutions to these problems will mostly aim to minimize the mean response time to demands. In practice, however, a demand may not be served from its nearest facility if that facility is engaged in serving other demands. This makes it a requirement to assign backup services so as to improve response time and service quality. The level of backup service is a key, strategic-level planning factor, and must be taken into consideration carefully. Moreover, in emergency service operations conducted in congested demand regions, demand assignment policy is another important factor that affects the system performance. Models failing to adopt sufficient levels of backup service and realistic demand assignment policies may significantly deteriorate the system performance.Considering the classic p-median problem (pMP) location model, this paper investigates the effects of backup service level, demand assignment policy, demand density, and number of facilities and their locations on the solution performance in terms of multiple metrics. For this purpose, we adopt a combined optimization and simulation approach. We will first modify the classic pMP to account for distances to backup services. Next, we employ a discrete event simulation to evaluate the performance of location schemes obtained from the deterministic mathematical model. Our results provide insights for decision-makers while planning ESS operations.     

A discrete competitive facility location model with variable attractiveness

We consider the discrete version of the competitive facility location problem in which new facilities have to be located by a new market entrant firm to compete against already existing facilities that may belong to one or more competitors. The demand is assumed to be aggregated at certain points in the plane and the new facilities can be located at predetermined candidate sites. We employ Huff's gravity-based rule in modelling the behaviour of the customers where the probability that customers at a demand point patronize a certain facility is proportional to the facility attractiveness and inversely proportional to the distance between the facility site and demand point. The objective of the firm is to determine the locations of the new facilities and their attractiveness levels so as to maximize the profit, which is calculated as the revenue from the customers less the fixed cost of opening the facilities and variable cost of setting their attractiveness levels. We formulate a mixed-integer nonlinear programming model for this problem and propose three methods for its solution: a Lagrangean heuristic, a branch-and-bound method with Lagrangean relaxation, and another branch-and-bound method with nonlinear programming relaxation. Computational results obtained on a set of randomly generated instances show that the last method outperforms the others in terms of accuracy and efficiency and can provide an optimal solution in a reasonable amount of time.     

Competitive facility location problem with attractiveness adjustment of the follower: A bilevel programming model and its solution

We are concerned with a problem in which a firm or franchise enters a market by locating new facilities where there are existing facilities belonging to a competitor. The firm aims at finding the location and attractiveness of each facility to be opened so as to maximize its profit. The competitor, on the other hand, can react by adjusting the attractiveness of its existing facilities with the objective of maximizing its own profit. The demand is assumed to be aggregated at certain points in the plane and the facilities of the firm can be located at predetermined candidate sites. We employ Huff’s gravity-based rule in modeling the behavior of the customers where the fraction of customers at a demand point that visit a certain facility is proportional to the facility attractiveness and inversely proportional to the distance between the facility site and demand point. We formulate a bilevel mixed-integer nonlinear programming model where the firm entering the market is the leader and the competitor is the follower. In order to find the optimal solution of this model, we convert it into an equivalent one-level mixed-integer nonlinear program so that it can be solved by global optimization methods. Apart from reporting computational results obtained on a set of randomly generated instances, we also compute the benefit the leader firm derives from anticipating the competitor’s reaction of adjusting the attractiveness levels of its facilities. The results on the test instances indicate that the benefit is 58.33% on the average.     

The capacitated maximal covering location problem with backup service

The maximal covering location problem has been shown to be a useful tool in siting emergency services. In this paper we expand the model along two dimensions — workload capacities on facilities and the allocation of multiple levels of backup or prioritized service for all demand points. In emergency service facility location decisions such as ambulance sitting, when all of a facility's resources are needed to meet each call for service and the demand cannot be queued, the need for a backup unit may be required. This need is especially significant in areas of high demand. These areas also will often result in excessive workload for some facilities. Effective siting decisions, therefore, must address both the need for a backup response facility for each demand point and a reasonable limit on each facility's workload. In this paper, we develop a model which captures these concerns as well as present an efficient solution procedure using Lagrangian relaxation. Results of extensive computational experiments are presented to demonstrate the viability of the approach.     

Planning for meals-on-wheels: algorithms and application

Home-delivered meals provision, also known as meals-on-wheels, is a volunteer-staffed activity for which little strategic planning is performed. We develop a Memetic Algorithm to solve the Home Delivered Meals Location-Routing Problem. This planning model addresses facility location, allocation of demand to facilities, and design of delivery routes, while balancing efficiency and effectiveness considerations. The case study presented on a large data set shows how trade-off curves, which are very useful for decision making, can be obtained by the method developed.     

The budget constrained r-interdiction median problem with capacity expansion

In this article, we elaborate on a budget constrained extension of the r-interdiction median problem with fortification (RIMF). The objective in the RIMF is to find the optimal allocation of protection resources to a given service system consisting of p facilities so that the disruptive effects of r possible attacks to the system are minimized. The defender of the system needs to fortify q facilities of the present system to offset the worst-case loss of r non-fortified facilities due to an interdiction in which the attacker’s objective is to cause the maximum possible disruption in the service level of the system. The defender-attacker relationship fits a bilevel integer programming (BIP) formulation where the defender and attacker take on the respective roles of the leader and the follower. We adopt this BIP formulation and augment it with a budget constraint instead of a predetermined number of facilities to be fortified. In addition, we also assume that each facility has a flexible service capacity, which can be expanded at a unit cost to accommodate the demand of customers who were serviced by some other interdicted facility before the attack. First, we provide a discrete optimization model for this new facility protection planning scenario with a novel set of closest assignment constraints. Then, to tackle this BIP problem we use an implicit enumeration algorithm performed on a binary tree. For each node representing a different fortification scheme, the attacker’s problem is solved to optimality using Cplex 11. We report computational results obtained on a test bed of 96 randomly generated instances. The article concludes with suggestions for future research.     

Hybrid heuristics for the capacitated lot sizing and loading problem with setup times and overtime decisions

The capacitated lot sizing and loading problem (CLSLP) deals with the issue of determining the lot sizes of product families/end items and loading them on parallel facilities to satisfy dynamic demand over a given planning horizon. The capacity restrictions in the CLSLP are imposed by constraints specific to the production environment considered. When a lot size is positive in a specific period, it is loaded on a facility without exceeding the sum of the regular and overtime capacity limits. Each family may have a different process time on each facility and furthermore, it may be technologically feasible to load a family only on a subset of existing facilities. So, in the most general case, the loading problem may involve unrelated parallel facilities of different classes. Once loaded on a facility, a family may consume capacity during setup time. Inventory holding and overtime costs are minimized in the objective function. Setup costs can be included if setups incur costs other than lost production capacity. The CLSLP is relevant in many industrial applications and may be generalized to multi-stage production planning and loading models. The CLSLP is a synthesis of three different planning and loading problems, i.e., the capacitated lot sizing problem (CLSP) with overtime decisions and setup times, minimizing total tardiness on unrelated parallel processors, and, the class scheduling problem, each of which is NP in the feasibility and optimality problems. Consequently, we develop hybrid heuristics involving powerful search techniques such as simulated annealing (SA), tabu search (TS) and genetic algorithms (GA) to deal with the CLSLP. Results are compared with optimal solutions for 108 randomly generated small test problems. The procedures developed here are also compared against each other in 36 larger size problems.     

Multi-objective competitive location problem with distance-based attractiveness and its best non-dominated solution

The multi-objective competitive location problem (MOCLP) with distance-based attractiveness is introduced. There are m potential competitive facilities and n demand points on the same plane. All potential facilities can provide attractiveness to the demand point which the facility attractiveness is represented as distance-based coverage of a facility, which is “full coverage” within the maximum full coverage radius, “no coverage” outside the maximum partial coverage radius, and “partial coverage” between those two radii. Each demand point covered by one of m potential facilities is determined by the greatest accumulated attractiveness provided the selected facilities and least accumulated distances between each demand point and selected facility, simultaneously. The tradeoff of maximum accumulated attractiveness and minimum accumulated distances is represented as a multi-objective optimization model. A proposed solution procedure to find the best non-dominated solution set for MOCLP is introduced. Several numerical examples and instances comparing with introduced and exhaustive method demonstrates the good performance and efficiency for the proposed solution procedure.     

Location of Multiple-Server Congestible Facilities for Maximizing Expected Demand, when Services are Non-Essential

We formulate a model for locating multiple-server, congestible facilities. Locations of these facilities maximize total expected demand attended over the region. The effective demand at each node is elastic to the travel time to the facility, and to the congestion at that facility. The facilities to be located are fixed, so customers travel to them in order to receive service or goods, and the demand curves at each demand node (which depend on the travel time and the queue length at the facility), are known. We propose a heuristic for the resulting integer, nonlinear formulation, and provide computational experience.     

The gravity p-median model

In this paper we propose a new model for the p-median problem. In the standard p-median problem it is assumed that each demand point is served by the closest facility. In many situations (for example, when demand points are communities of customers and each customer makes his own selection of the facility) demand is divided among the facilities. Each customer selects a facility which is not necessarily the closest one. In the gravity p-median problem it is assumed that customers divide their patronage among the facilities with the probability that a customer patronizes a facility being proportional to the attractiveness of that facility and to a decreasing utility function of the distance to the facility.     

Competitive facility location on decentralized supply chains

This paper addresses a novel competitive facility location problem about a firm that intends to enter an existing decentralized supply chain comprised of three tiers of players with competition: manufacturers, retailers and consumers. It first proposes a variational inequality for the supply chain network equilibrium model with production capacity constraints, and then employs the logarithmic-quadratic proximal prediction–correction method as a solution algorithm. Based on this model, this paper develops a generic mathematical program with equilibrium constraints for the competitive facility location problem, which can simultaneously determine facility locations of the entering firm and the production levels of these facilities so as to optimize an objective. Subsequently, a hybrid genetic algorithm that incorporates with the logarithmic-quadratic proximal prediction–correction method is developed for solving the proposed mathematical program with an equilibrium constraint. Finally, this paper carries out some numerical examples to evaluate proposed models and solution algorithms.     

竞争环境下截流设施选址问题

研究企业新建设施时,市场上已有设施存在的情况下,使本企业总体利润最大的截流设施选址问题。在一般截留设施选址模型的基础上引入引力模型,消费者到某个设施接受服务的概率与偏离距离及设施的吸引力相关,同时设施的建设费用与设施吸引力正相关,建立非线性整数规划模型并使用贪婪算法进行求解。数值分析表明,该算法求解速度快,模型计算精度较高。     

Facility location for market capture when users rank facilities by shorter travel and waiting times

A firm wants to locate several multi-server facilities in a region where there is already a competitor operating. We propose a model for locating these facilities in such a way as to maximize market capture by the entering firm, when customers choose the facilities they patronize, by the travel time to the facility and the waiting time at the facility. Each customer can obtain the service or goods from several (rather than only one) facilities, according to a probabilistic distribution. We show that in these conditions, there is demand equilibrium, and we design an ad hoc heuristic to solve the problem, since finding the solution to the model involves finding the demand equilibrium given by a nonlinear equation. We show that by using our heuristic, the locations are better than those obtained by utilizing several other methods, including MAXCAP, p-median and location on the nodes with the largest demand.     

IP modeling of the survivable hop constrained connected facility location problem

We consider a generalized version of the rooted connected facility location problem which occurs in planning of telecommunication networks with both survivability and hop-length constraints. Given a set of client nodes, a set of potential facility nodes including one predetermined root facility, a set of optional Steiner nodes, and the set of the potential connections among these nodes, that task is to decide which facilities to open, how to assign the clients to the open facilities, and how to interconnect the open facilities in such a way, that the resulting network contains at least λ edge-disjoint paths, each containing at most H edges, between the root and each open facility and that the total cost for opening facilities and installing connections is minimal. We study two IP models for this problem and present a branch-and-cut algorithm based on Benders decomposition for finding its solution. Finally, we report computational results.     

Optimizing pricing and location decisions for competitive service facilities charging uniform price

In this paper, we present the problem of optimizing the location and pricing for a set of new service facilities entering a competitive marketplace. We assume that the new facilities must charge the same (uniform) price and the objective is to optimize the overall profit for the new facilities. Demand for service is assumed to be concentrated at discrete demand points (customer markets); customers in each market patronize the facility providing the highest utility. Customer demand function is assumed to be elastic; the demand is affected by the price, facility attractiveness, and the travel cost for the highest-utility facility. We provide both structural and algorithmic results, as well as some managerial insights for this problem. We show that the optimal price can be selected from a certain finite set of values that can be computed in advance; this fact is used to develop an efficient mathematical programming formulation for our model.     

The equitable location problem on the plane

This paper considers the problem of locating M facilities on the unit square so as to minimize the maximal demand faced by each facility subject to closest assignments and coverage constraints. Focusing on uniform demand over the unit square, we develop upper and lower bounds on feasibility of the problem for a given number of facilities and coverage radius. Based on these bounds and numerical experiments we suggest a heuristic to solve the problem. Our computational results show that the heuristic is very efficient, as the average gap between its solutions and the lower bound is 4.34%.     

Environmental site evaluation of waste management facilities embedded into EUGÈNE model: A multicriteria approach

EUGÈNE is a sophisticated mixed integer linear programming model developed to help regional decision makers on long-term planning for solid waste management activities. The model removes almost every limitations encountered in other waste management models and contains a large quantity of variables and constraints. The method used to embed waste management environmental parameters in the EUGÈNE model consists in building global impact index (GII) for all site/facility combinations. First, an environmental and spatial evaluation of waste management facilities over sites is based on qualitative and quantitative criteria measuring biophysical and social impacts. Spatial analysis is carried out by geographical information system routines. Then, a multicriteria analysis ranks all site/facility combinations, according to their global performance based on all criteria. The net flow, computed by the PROMETHEE multicriteria outranking method, is considered as a GII to be embedded into EUGÈNE. The model objective function is thus modified to minimize total system cost and GII. Some practical results obtained for the City of Montreal are discussed.