A mean-risk mixed integer nonlinear program for transportation network protection

This paper focuses on transportation network protection to hedge against extreme events such as earthquakes. Traditional two-stage stochastic programming has been widely adopted to obtain solutions under a risk-neutral preference through the use of expectations in the recourse function. In reality, decision makers hold different risk preferences. We develop a mean-risk two-stage stochastic programming model that allows for greater flexibility in handling risk preferences when allocating limited resources. In particular, the first stage minimizes the retrofitting cost by making strategic retrofit decisions whereas the second stage minimizes the travel cost. The conditional value-at-risk (CVaR) is included as the risk measure for the total system cost. The two-stage model is equivalent to a nonconvex mixed integer nonlinear program (MINLP). To solve this model using the Generalized Benders Decomposition (GBD) method, we derive a convex reformulation of the second-stage problem to overcome algorithmic challenges embedded in the non-convexity, nonlinearity, and non-separability of first- and second-stage variables. The model is used for developing retrofit strategies for networked highway bridges, which is one of the research areas that can significantly benefit from mean-risk models. We first justify the model using a hypothetical nine-node network. Then we evaluate our decomposition algorithm by applying the model to the Sioux Falls network, which is a large-scale benchmark network in the transportation research community. The effects of the chosen risk measure and critical parameters on optimal solutions are empirically explored

Approximate and exact convexification approaches for solving two-stage mixed-integer recourse models

Many practical decision-making problems are subject to uncertainty. A powerful class of mathematical models designed for these problems is the class of mixed-integer recourse models. Such models have a wide range of applications in, e.g., healthcare, energy, and finance. They permit integer decision variables to accurately model, e.g., on/off restrictions or natural indivisibilities. The additional modelling flexibility of integer decision variables, however, comes at the expense of models that are significantly harder to solve. The reason is that including integer decision variables introduces non-convexity in the model, which poses a significant challenge for state-of-the-art solvers.In this thesis, we contribute to better decision making under uncertainty by designing efficient solution methods for mixed-integer recourse models. Our approach is to address the non-convexity caused by integer decision variables by using convexification. That is, we construct convex approximating models that closely approximate the original model. In addition, we derive performance guarantees for the solution obtained by solving the approximating model. Finally, we extensively test the solution methods that we propose and we find that they consistently outperform traditional solution methods on a wide range of benchmark instances

Algorithms for Stochastic Integer Programs Using Fenchel Cutting Planes

This dissertation develops theory and methodology based on Fenchel cutting planes for solving stochastic integer programs (SIPs) with binary or general integer variables in the second-stage. The methodology is applied to auto-carrier loading problem under uncertainty. The motivation is that many applications can be modeled as SIPs, but this class of problems is hard to solve. In this dissertation, the underlying parameter distributions are assumed to be discrete so that the original problem can be formulated as a deterministic equivalent mixed-integer program. The developed methods are evaluated based on computational experiments using both real and randomly generated instances from the literature. We begin with studying a methodology using Fenchel cutting planes for SIPs with binary variables and implement an algorithm to improve runtime performance. We then introduce the stochastic auto-carrier loading problem where we present a mathematical model for tactical decision making regarding the number and types of auto-carriers needed based on the uncertainty of availability of vehicles. This involves the auto-carrier loading problem for which actual dimensions of the vehicles, regulations on total height of the auto-carriers and maximum weight of the axles, and safety requirements are considered. The problem is modeled as a two-stage SIP, and computational experiments are performed using test instances based on real data. Next, we develop theory and a methodology for Fenchel cutting planes for mixed integer programs with special structure. Integer programs have to be solved to generate a Fenchel cutting plane and this poses a challenge. Therefore, we propose a new methodology for constructing a reduced set of integer points so that the generation of Fenchel cutting planes is computationally favorable. We then present the computational results based on randomly generated instances from the literature and discuss the limitations of the methodology. We finally extend the methodology to SIPs with general integer variables in the second-stage with special structure, and study different normalizations for Fenchel cut generation and report their computational performance

Modelling and solution methods for stochastic optimisation

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University.In this thesis we consider two research problems, namely, (i) language constructs for modelling stochastic programming (SP) problems and (ii) solution methods for processing instances of different classes of SP problems. We first describe a new design of an SP modelling system which provides greater extensibility and reuse. We implement this enhanced system and develop solver connections. We also investigate in detail the following important classes of SP problems: singlestage SP with risk constraints, two-stage linear and stochastic integer programming problems. We report improvements to solution methods for single-stage problems with second-order stochastic dominance constraints and two-stage SP problems. In both cases we use the level method as a regularisation mechanism. We also develop novel heuristic methods for stochastic integer programming based on variable neighbourhood search. We describe an algorithmic framework for implementing decomposition methods such as the L-shaped method within our SP solver system. Based on this framework we implement a number of established solution algorithms as well as a new regularisation method for stochastic linear programming. We compare the performance of these methods and their scale-up properties on an extensive set of benchmark problems. We also implement several solution methods for stochastic integer programming and report a computational study comparing their performance. The three solution methods, (a) processing of a single-stage problem with second-order stochastic dominance constraints, (b) regularisation by the level method for two-stage SP and (c) method for solving integer SP problems, are novel approaches and each of these makes a contribution to knowledge.Financial support was obtained from OptiRisk Systems