Lost-sales inventory systems with a service level criterion

Competitive retail environments are characterized by service levels and lost sales in case of excess demand. We contribute to research on lost-sales models with a service level criterion in multiple ways. First, we study the optimal replenishment policy for this type of inventory system as well as base-stock policies and (R, s, S) policies. Furthermore, we derive lower and upper bounds on the order-up-to level, and we propose efficient approximation procedures to determine the order-up-to level. The procedures find values of the inventory control variables that are close to the best (R, s, S) policy and comply to the service level restriction for most of the instances, with an average cost increase of 2.3% and 1.2% for the case without and with fixed order costs, respectively.     

Economic application in a finite capacity multi-channel queue with second optional channel

We consider a finite capacity M/M/R queue with second optional channel. The interarrival times of arriving customers follow an exponential distribution. The service times of the first essential channel and the second optional channel are assumed to follow an exponential distribution. As soon as the first essential service of a customer is completed, a customer may leave the system with probability (1 − θ) or may opt for the second optional service with probability θ (0 ? θ ? 1). Using the matrix-geometric method, we obtain the steady-state probability distributions and various system performance measures. A cost model is established to determine the optimal solutions at the minimum cost. Finally, numerical results are provided to illustrate how the direct search method and the tabu search can be applied to obtain the optimal solutions. Sensitivity analysis is also investigated.     

Vacation policies for machine repair problem with two type spares

This paper studies the machine repair problem consisting of M operating machines with two types of spare machines (S = S₁ + S₂), and R servers (repairmen) who leave for a vacation of random length when there are no failed machines queuing up for repair in the repair facility. At the end of the vacation the servers return and operate two vacation policies. First, the servers take vacations repeatedly until they find the repair facility has at least one waiting failed machine in the queue. Second, the servers do not take a vacation again and remain idle until the first arriving failed machine arrives, which starts a busy period in the repair facility. For both policies, the servers have two service rates for repair-slow and fast. The matrix geometric theory is used to find the steady-state probabilities of the number of failed machines in the system as well as the performance measures. Some special cases are given. A direct search algorithm is used to simultaneously determine the optimal values of the number of two types of spares and the number of servers while maintaining a minimum specified level of system availability.     

Three formulas for eigenfunctions of integrable Schrödinger operators

We give three formulas for meromorphic eigenfunctions (scatteringstates) of Sutherlandsintegrable N-body Schrödinger operators and their generalizations.The first is an explicit computation of the Etingof–Kirillov tracesof intertwining operators, the second an integral representationof hypergeometric type, and the third is a formula of Bethe ansatz type.The last two formulas are degenerations of elliptic formulasobtained previously in connection with theKnizhnik–Zamolodchikov–Bernardequation. The Bethe ansatz formulas in the elliptic case are reviewed and discussed in more detail here: Eigenfunctionsare parametrized by a Hermite–Bethe variety, a generalizationof the spectral variety of the Lamé operator.We also give the q-deformed version of ourfirst formula. In the scalar slN case, this gives common eigenfunctionsof the commuting Macdonald–Rujsenaars difference operators.     

Two-dimensional cargo overbooking models

This paper introduces two-dimensional (weight and volume) overbooking problems arising mainly in the cargo revenue management, and compares them with one-dimensional problems. It considers capacity spoilage and cargo offloading costs, and minimizes their sum. For one-dimensional problems, it shows that the optimal overbooking limit does not change with the magnitude of the booking requests. In two-dimensional problems, the overbooking limit is replaced by a curve. The curve, along with the volume and weight axes, encircles the acceptance region. The booking requests are accepted if they fall within this region. We present Curve (Cab) and Rectangle (Rab) models. The boundary of the acceptance region in the Cab (resp. Rab) model is a curve (resp. rectangle). The optimal curve for the Cab model is shown to be unique and continuous. Moreover, it can be obtained by solving a series of simple equations. Finding the optimal rectangle for the Rab model is more challenging, so we propose an approximate rectangle. The approximate rectangle is a limiting solution in the sense that it converges to the optimal rectangle as the booking requests increase. The approximate rectangle is numerically shown to yield costs that are very close to the optimal costs.     

节能情形下制造企业动态批量生产策略与节能方式选择

随着诸多严厉的环保政策的实行,节能管理已成为我国高耗能制造企业运营管理中的一个重要内容。通过构建自我节能、节能效益分享和节能量保证三种节能方式下制造企业的动态批量生产决策模型,首先分析了给定节能方式下制造企业的最优动态批量生产策略,其次探讨了制造企业最优节能方式选择。结果表明,最优动态批量生产策略满足“零库存生产”规则。进一步的数值实验分析结果表明:面临前两种节能方式时,当节能服务公司与制造企业的节能投资成本系数比率和采用合同能源管理的预计节能率较小时,制造企业应选择后者,反之应选择前者;面临后两种节能方式时,当采用合同能源管理的预计节能率较小时,制造企业选择两者无差异,反之应选择前者。     

Hölder and Minkowski type inequalities for pseudo-integral

There are proven generalizations of the Hölder’s and Minkowski’s inequalities for the pseudo-integral. There are considered two cases of the real semiring with pseudo-operations: one, when pseudo-operations are defined by monotone and continuous function g, the second semiring ([a, b], sup, ⊙), where ⊙ is generated and the third semiring where both pseudo-operations are idempotent, i.e., ⊕ = sup and ⊙ = inf.     

Quasi-steady flight to quasi-steady flight transition in a windshear: Trajectory optimization and guidance

This paper is concerned with the near-optimum guidance of an aircraft from quasi-steady flight to quasi-steady flight in a windshear. The take-off problem is considered with reference to flight in a vertical plane. In addition to the horizontal shear, the presence of a downdraft is considered. It is assumed that the power setting is held at the maximum value and that the aircraft is controlled through the angle of attack. Inequality constraints are imposed on both the angle of attack and its time derivative.First, trajectory optimization is considered. The optimal transition problem is formulated as a Chebyshev problem of optimal control: the performance index being minimized is the peak value of the modulus of the difference between the absolute path inclination and a reference value, assumed constant. Two types of optimal trajectories are studied: type 1 is concerned with gamma recovery (recovery of the initial value of the relative path inclination); and type 2 is concerned with quasisteady flight recovery (recovery of the initial values of the relative velocity, the relative path inclination, and the relative angle of attack). The numerical results show that the type 1 trajectory and the type 2 trajectory are nearly the same in the shear portion, while they diverge to a considerable degree in the aftershear portion of the optimal trajectory.Next, trajectory guidance is considered. A guidance scheme is developed so as to achieve near-optimum quasi-steady flight recovery in a windshear. The guidance scheme for quasi-steady flight recovery includes three parts in sequence. The first part refers to the shear portion of the trajectory and is based on the result that this portion of the trajectory depends only mildly on the boundary conditions; therefore, any of the guidance schemes already developed for type 1 trajectories can be employed (for instance, variable gamma guidance). The second part (constant gamma guidance) refers to the initial aftershear portion of the trajectory and is designed to achieve almost velocity recovery. The third part (constant rate of climb guidance) refers to the final aftershear portion of the trajectory and is designed to achieve almost complete restoration of the initial quasi-steady state.While the shear guidance and the initial aftershear guidance employ constant gain coefficients, the final aftershear guidance employs a variable gain coefficient. This is done in order to obtain accuracy and prompt response, while avoiding oscillations and overshoots. The numerical results show that the guidance scheme for quasi-steady flight recovery yields a transition from quasi-steady flight to quasi-steady flight which is close to that of the optimal trajectory, ensures the restoration of the initial quasi-steady state, and has good stability properties.This paper is based on Refs. 1 and 2.This research was supported by NASA-Langley Research Center, Grant No. NAG-1-516, and by Boeing Commercial Aircraft Company. The authors are indebted to Dr. R. L. Bowles, NASA-Langley Research Center, for helpful discussions.     

Difference schemes of exponential type for singularly perturbed Volterra integro-differential problems

Extended one-step schemes of exponential type are introduced for solving singularly perturbed Volterra integro-differential problems. These schemes are of order (m + 1), m = 0, 1, 2, …, when the perturbation parameter, ε, is fixed. These schemes have the property that if ε is of order h they reduced to first order of accuracy and optimal when ε → 0. Stability analysis of these schemes are presented. Numerical results and comparisons with other schemes are presented.     

Preemptive multiprocessor order scheduling to minimize total weighted flowtime

Consider m identical machines in parallel, each of which can produce k different product types. There is no setup cost when the machines switch from producing one product type to another. There are n orders each of which requests various quantities of the different product types. All orders are available for processing at time t = 0, and preemption is allowed. Order i has a weight w_i and its completion time is the time when its last requested product type finishes. Our goal is to find a preemptive schedule such that the total weighted completion time ∑_wiCi$\sum w_i{C}_i$ is minimized. We show that this problem is NP-hard even when all jobs have identical weights and there are only two machines. Motivated by the computational complexity of the problem, we propose a simple heuristic and show that it obeys a worst-case bound of 2 − 1/m. Finally, empirical studies show that our heuristic performs very well when compared with a lower bound of the optimal cost.     

机票预定问题的多阶段决策模型

本文讨论了机票预定的一种方法 .通过建立多阶段决策模型 ,将订票时期分成若干个阶段 ,在每一个阶段航空公司对乘客要求订票作出不同的反应 ,保证了检票时准备登机的人数与飞机上的座位数目相当接近     

Batch delivery scheduling with batch delivery cost on a single machine

We consider a scheduling problem in which n independent and simultaneously available jobs are to be processed on a single machine. The jobs are delivered in batches and the delivery date of a batch equals the completion time of the last job in the batch. The delivery cost depends on the number of deliveries. The objective is to minimize the sum of the total weighted flow time and delivery cost. We first show that the problem is strongly NP-hard. Then we show that, if the number of batches is B, the problem remains strongly NP-hard when B ? U for a variable U ? 2 or B ? U for any constant U ? 2. For the case of B ? U, we present a dynamic programming algorithm that runs in pseudo-polynomial time for any constant U ? 2. Furthermore, optimal algorithms are provided for two special cases: (i) jobs have a linear precedence constraint, and (ii) jobs satisfy the agreeable ratio assumption, which is valid, for example, when all the weights or all the processing times are equal.     

Dynamic cruise ship revenue management

Recently, it has been recognized that revenue management of cruise ships is different from that of airlines or hotels. Among the main differences is the presence of multiple capacity constraints in cruise ships, i.e., the number of cabins in different categories and the number of lifeboat seats, versus a single constraint in airlines and hotels (i.e., number of seats or rooms). We develop a discrete-time dynamic capacity control model for a cruise ship characterized by multiple constraints on cabin and lifeboat capacities. Customers (families) arrive sequentially according to a stochastic process and request one cabin of a certain category and one or more lifeboat seats. The cruise ship revenue manager decides which requests to accept based on the remaining cabin and lifeboat capacities at the time of an arrival as well as the type of the arrival. We show that the opportunity cost of accepting a customer is not always monotone in the reservation levels or time. This non-monotone behavior implies that “conventional” booking limits or critical time periods capacity control policies are not optimal. We provide analysis and insights justifying the non-monotone behavior in our cruise ship context. In the absence of monotonicity, and with the optimal solution requiring heavy storage for “large” (industry-size) problems, we develop several heuristics and thoroughly test their performance, via simulation, against the optimal solution, well-crafted upper bounds, and a first-come first-served lower bound. Our heuristics are based on rolling-up the multi-dimensional state space into one or two dimensions and solving the resulting dynamic program (DP). This is a strength of our approach since our DP-based heuristics are easy to understand, solve and analyze. We find that single-dimensional heuristics based on decoupling the cabins and lifeboat problems perform quite well in most cases.     

Inventory control in multi-channel retail

In recent years multi-channel retail systems have received increasing interest. Partly due to growing online business that serves as a second sales channel for many firms, offering channel specific prices has become a common form of revenue management. We analyze conditions for known inventory control policies to be optimal in presence of two different sales channels. We propose a single item lost sales model with a lead time of zero, periodic review and nonlinear non-stationary cost components without rationing to realistically represent a typical web-based retail scenario. We analyze three variants of the model with different arrival processes: demand not following any particular distribution, Poisson distributed demand and a batch arrival process where demand follows a Pòlya frequency type distribution. We show that without further assumptions on the arrival process, relatively strict conditions must be imposed on the penalty cost in order to achieve optimality of the base stock policy. We also show that for a Poisson arrival process with fixed ordering costs the model with two sales channels can be transformed into the well known model with a single channel where mild conditions yield optimality of an (s, S) policy. Conditions for optimality of the base stock and (s, S) policy for the batch arrival process with and without fixed ordering costs, respectively, are presented together with a proof that the batch arrival process provides valid upper and lower bounds for the optimal value function.     

Optimal access control for broadband services: Stochastic knapsack with advance information

Suppose that call reservation requests of K different types arrive randomly at a single capacitated link where each in turn must be accepted or rejected. Each request requires an amount of bandwidth and a random duration, both depending on its type. An accepted call generates an immediate fixed revenue followed by a variable revenue per unit time. We assume that each call's requested duration is known when it arrives and that the state of the link is constantly monitored. The problem is to find a call acceptance policy that maximizes the long-run average revenue per unit time generated by the accepted calls. We propose and investigate threshold-type control policies that use the known duration of arriving calls as well as the current link status. Two main contributions result. First are interpretable analytical results for the case of K = 1 that can also be applied to the case of K > 1 using complete partitioning scheme. Second, we illustrate how to apply stochastic optimization techniques to compute optimal thresholds in the general case. We also include the results of initial simulation experiments.     

Large time behavior of differential equations with drifted periodic coefficients modeling carbon storage in soil

This paper is concerned with the linear ODE in the form y′(t) = λρ(t)y(t) + b(t), λ < 0 which represents a simplified storage model of the carbon in the soil. In the first part, we show that, for a periodic function ρ(t), a linear drift in the coefficient b(t) involves a linear drift for the solution of this ODE. In the second part, we extend the previous results to a classical heat non-homogeneous equation. The connection with an analytic semi-group associated to the ODE equation is considered in the third part. Numerical examples are given.     

Optimal Interplanetary Orbital Transfers via Electrical Engines

The Hohmann transfer theory, developed in the 19th century, is the kernel of orbital transfer with minimum propellant mass by means of chemical engines. The success of the Deep Space 1 spacecraft has paved the way toward using advanced electrical engines in space. While chemical engines are characterized by high thrust and low specific impulse, electrical engines are characterized by low thrust and hight specific impulse. In this paper, we focus on four issues of optimal interplanetary transfer for a spacecraft powered by an electrical engine controlled via the thrust direction and thrust setting: (a) trajectories of compromise between transfer time and propellant mass, (b) trajectories of minimum time, (c) trajectories of minimum propellant mass, and (d) relations with the Hohmann transfer trajectory. The resulting fundamental properties are as follows:

New types of 3-D systems of quadratic differential equations with chaotic dynamics based on Ricker discrete population model

The wide class of 3-D autonomous systems of quadratic differential equations, in each of which either there is a couple of coexisting limit cycles or there is a couple of coexisting chaotic attractors, is found. In the second case the couple consists of either Lorentz-type attractor and another attractor of a new type or two Lorentz-type attractors. It is shown that the chaotic behavior of any system of the indicated class can be described by the Ricker discrete population model: z_i+1 = z_i exp(r − z_i), r > 0, z_i > 0, i = 0, 1, … . The values of parameters, at which in the 3-D system appears either the couple of limit cycles or the couple of chaotic attractors, or only one limit cycle, or only one sphere-shaped chaotic attractor, are indicated. Examples are given.     

Revenue Models and Policies for the Car Rental Industry

In this paper, we consider the application of revenue management techniques in the context of the car rental industry. In particular, the paper presents a dynamic programming formulation for the problem of assigning cars of several categories to different segments of customers, with rental requests arising dynamically and randomly with time. Customers make a rental request for a given type of car, for a given number of days at a given pickup time. The rental firm can satisfy the demand for a given product with either the product requested or with a car of at most one category superior to that initially required, in this case an “upgrade” can take place. The one-way rental scenario, which allows the possibility of the rental starting and ending at different locations, is also addressed. In the framework considered, the logistic operator has to decide whether to accept or reject a rental request. Since the proposed dynamic programming formulations are impractical due to the curse of dimensionality, linear programming approximations are used to derive revenue management decision policies for the operator. Indeed, primal and dual acceptance policies are developed (i.e. booking limits, bid prices) and their effectiveness is assessed on the basis of an extensive computational phase.     

Online-optimization of multi-elevator transport systems with reoptimization algorithms based on set-partitioning models

We develop and experimentally compare policies for the control of a system of k elevators with capacity one in a transport environment with ? floors, an idealized version of a pallet elevator system in a large distribution center of the Herlitz PBS AG in Falkensee. Each elevator in the idealized system has an individual waiting queue of infinite capacity. On each floor, requests arrive over time in global waiting queues of infinite capacity. The goal is to find a policy that, without any knowledge about future requests, assigns an elevator to each request and a schedule to each elevator so that certain expected cost functions (e.g., the average or the maximal flow times) are minimized. We show that a reoptimization policy for minimizing average squared waiting times can be implemented to run in real-time (1 s) using dynamic column generation. Moreover, in discrete event simulations with Poisson input it outperforms other commonly used policies like multi-server variants of greedy and nearest neighbor.