Valid Inequalities for Optimal Transmission Switching
|
|
- Brian Marshall
- 6 years ago
- Views:
Transcription
1 Valid Inequalities for Optimal Transmission Switching Hyemin Jeon Jeff Linderoth Jim Luedtke Dept. of ISyE UW-Madison Burak Kocuk Santanu Dey Andy Sun Dept. of ISyE Georgia Tech 19th Combinatorial Optimization Workshop Aussois, France, January, 2015 Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 1 / 28
2 Power Background Notation This is a Power Systems Talk But I don t know much about power systems Hopefully you will find the mathematics interesting Economic Dispatch Focus today is on a simple problem of meeting demand for power at minimum cost Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 2 / 28
3 Power Background Notation Power Grid Networks Look Weird Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 3 / 28
4 Power Background Notation It s Just a Network Power Network: (N, A) with... G N: generation nodes D N: demand nodes Load forecasts (MW) b i for i D Generation cost ($/MW) c i and capability p i (MW) for i G Peak load rating (MW) u ij for (i, j) A Economic Dispatch Problem Determine power generation levels for i G and power transmission levels for (i, j) A to meet demands b i, i D, at minimum cost Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 4 / 28
5 Power Background Optimal Power Flow Power Flow Electric power grids follow the laws of physics, characterized by nonlinear, nonconvex equations Direct control is difficult We cannot dictate how power will flow. Making many assumptions, we can model the real power flow using the Direct Current (DC) model. DC Power Flow Assumption The (real) power x ij transmit over line (i, j) A is proportional to angle differences (θ i θ j ) at the endpoint nodes i N and j N. x ij = α ij (θ i θ j ) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 5 / 28
6 Power Background Economic Dispatch Linear Program for (DC) Economic Dispatch s.t. j:(i,j) E min c i p i x,p,θ i G x ij j:(j,i) E x ji = u ij x ij u ij p i p i p i p i d i i G i D 0 i N \ G \ D (i, j) E i G x ij = α ij (θ i θ j ) x ij R p i R + θ i R (i, j) E i G i N (i, j) E x, θ need not be 0 Bounds on x, but no a priori bounds on θ (Usually derived from bounds on x) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 6 / 28
7 Power Background Tradeoffs MCNF++ This economic dispatch problem is just a min cost network flow problem with some additional potential constraints The potential drop (θ A θ D ) must be the same aloing the paths: A B D and A C D B A D C Braess Paradox If line (C, D) didn t exist, I wouldn t have to enforce this potential balance constraint. Thus, removing lines of the transmission network may actually increase the efficiency of delivery. Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 7 / 28
8 Power Background Tradeoffs Transmission Switching Tradeoff Having Lines Allows You to Send Flow: U ij x ij U ij (i, j) E Having Lines Induces Constraints in the Network: x ij = α ij (θ i θ j ) (i, j) E Fisher, O Neill & Ferris ( 08) show that efficiency improved by switching off transmission lines Lines Off % Improvement 1 6.3% % % % 24.9% Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 8 / 28
9 Transmission Switching Complexity Problem Complexity DC Transmission Switching Given: A network G = (N, A) with arc capacities and susceptances (u ij, α ij ) (i, j) A, generation levels p i i G, demand levels b i i D. Does there exist a subset of arcs S A such that deactivating arcs in S leads to a feasible DC power flow? Theorem DC Transmission Switching is NP-Complete Reduction from subset-sum The problem remains hard... Even if there are a polynomial number of cycles in the network Even on a series-parallel graph with only one supply/demand pair Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 9 / 28
10 Transmission Switching Formulations Switching Off Lines Regular Flow Constraints x ij = α ij (θ i θ j ) U ij x ij U ij (i, j) E (i, j) E Let z ij {0, 1} (i, j) A, Switched Flow Constraints x ij = α ij z ij (θ i θ j ) (i, j) E If (and only if) there are finite bounds on flow, then one can write a MILP formulation (Fisher, O Neil, and Ferris 08). z ij = 1 line (i, j) A is used Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 10 / 28
11 Transmission Switching Formulations MILP Formulation s.t. j:(i,j) E min x,p,θ,z x ij c i p i i G j:(j,i) E x ij = p i d i U ij z ij x ij U ij z ij α ij (θ i θ j ) x ij + M(1 z ij ) 0 α ij (θ i θ j ) x ij M(1 z ij ) 0 L i θ i L i p i p i p i z ij {0, 1} i G i D 0 i N \ G \ D (i, j) E (i, j) E (i, j) E i N i G (i, j) E Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 11 / 28
12 Transmission Switching The Challenge Throwing Down the Gauntlet Hedman, Ferris, O Neill, Fisher, Oren, (2010) state When solving the transmission switching problem,... the techniques for closing the optimality gap, specifically improving the lower bound, are largely ineffective. So they resort to a variety of heuristic, ad-hoc techniques to get good solutions to the MILP they propose My good colleague and continuous optimizer Michael Ferris impugns the good name of integer programmers by stating that we are not smart enough to solve DC transmission switching You will later see that CPLEX v12 is already orders of magnitude better than CPLEX v9 on DC transmission switching instances But still it s not good enough for large-scale networks... Thus we have... Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 12 / 28
13 Transmission Switching The Challenge Ferris s Challenge to Integer Programmers Solve realistically-sized DC transmission switching instances to provable optimality As integer programmers, we would like to rise to the challenge, and improve these ineffective lower bound techniques. The Padberg Way Let s study the mathematical structure of the problem, create a useful relaxation of the problem, and improve our description of the relaxation through cutting planes (facets) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 13 / 28
14 Shagadelic Cycle Inequalities Intuition Key (Simple) Insight?! A Assume (WLOG) that α ij = 1 We can just set x ij = α ij x ij and scale u ij by α ij Then we have... B C x AB = θ A θ B x BC = θ B θ C x CA = θ C θ A x AB + x BC + x CA = 0 The potential constraints essentially (only) enforce that flow around a cycle is zero. If you didn t forget everything from your introductory electrical engineering class (like I did), then you will recognize this as Kirchoff s Voltage Law. Insight We should focus on what goes on around a cycle and try to model this in a better way Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 14 / 28
15 Shagadelic Cycle Inequalities Intuition The Padberg Way Simple IP People (like me) Like Simple Sets Directed cycle G = (V, C), with V = [n], C = {(i, i + 1) i [n 1]} {(n, 1)}: { C = (x, θ, z) R 2n {0, 1} n u ij x ij u ij (i, j) C } z ij (θ i θ j ) = x ij (i, j) C The inequalities in this set model the potential drop across each arc in a cycle Even though C has the nonlinear equations z ij (θ i θ j ) = x ij, it is the union of 2 n polyhedra, so cl conv(c) is a polyhedron. The Padberg Way Structure, Structure, Structure! Even though C is just a relaxation of the true problem, we hope that by generating valid inequalities for C, we can improve performance of IP approaches Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 15 / 28
16 Shagadelic Cycle Inequalities Theorems Now We Do Math { C = (x, θ, z) R 2n {0, 1} n u ij x ij u ij (i, j) C } z ij (θ i θ j ) = x ij (i, j) C Theorem For S C such that u(s) > u(c\s), the shagadeliccycle inequalities (sci) x(s) + a C β S az a b S (1) x(s) + a C β S az a b S (2) are valid for C, where β S a = u(s \ a) u(c \ S) a C b S = (n 1)(2u(S) u(c)) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 16 / 28
17 Shagadelic Cycle Inequalities Theorems Shagadelic-Cycle Inequalities, Example x 1 2z 1 x 2 4z 2 x 3 3z 3x1 + x2 + z1 z2 + 3z3 6 S = {1, 2} x 1 + x 3 z 1 + z 2 2z 3 2 S = {1, 3} x 2 + x 3 + 5z 1 + z 2 + 2z 3 10 S = {2, 3} x 1 + x 2 + x 3 + 7z 1 + 5z 2 + 6z 3 18 S = {1, 2, 3} Logic Enforced For S = {1, 2}, if z 1 = z 2 = 1, then { 6 z3 = 0 x 1 + x 2 3 z 3 = 1 For S = {1, 3}, if z 1 = z 3 = 1, then { 5 z2 = 0 x 1 + x 3 4 z 2 = 1 Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 17 / 28
18 Shagadelic Cycle Inequalities Theorems Proofs! Yeah, Baby! Theorem If S C, and u(c \ S) < u(s), then the shagadelic-cycle inequalities (sci) are facet-defining for cl conv(c). Thus, all 2 n inequalities are necessary in the description of cl conv(c) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 18 / 28
19 Shagadelic Cycle Inequalities Theorems Even More Proofs Along with some other trivial inequalities, the shagadelic cycle inequalities are sufficient to describe the convex hull of C cl conv(c) = { (x, θ, z) R 3n u ij z ij x ij u ij z ij (i, j) C z ij 1 (i, j) C x(s) + a C β S az a b S S C : u(s) > u(c \ S) x(s) + } β S az a b S S C : u(s) > u(c \ S) a C Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 19 / 28
20 Separation Heuristic Can We Use the sci? Given solution ^x R n +, ^z [0, 1] n, the separation problem for (sci) is where max C A:C is a cycle max {^x(s) + (β S ) ^z b S }, S C:2u(S) u(c) β S a = u(s \ a) u(c \ S) b S = (n 1)(2u(S) u(c)) a C Theorem The separation problem for (sci) is NP-Hard Theorem Something that seems like it must be true, but I can t prove it. Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 20 / 28
21 Separation Heuristic Simple Observations Observation: If a C ^za C 1, then (^x, ^z) cannot be violated by any (sci) This suggests a two-phase separation heuristic. Separation Heuristic 1 Find a necessary cycle C such that a C ^za > C 1 2 Find S C in the given cycle Do (1) by (truncated) enumeration Given C, algebra shows that (2) is a knapsack problem: K C = 1 a C (1 ^za) 0 ^v a = ^x a + u a^z a 2u a( e C\a (1 ^ze)) ν = max ^v ay a u ay a 1 2 u(c) y {0,1} n a C a C If ν u(c)k C > 0, then (sci) is violated by (^x, ^z) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 21 / 28
22 Separation Heuristic P=NP I Can Solve the Knapsack Problem in Polynomial Time! Since I have proved that P = NP, the Clay Mathematics Institute should pay me... Not really, it is just that this specific knapsack problem is easy Take the items: S C = {a C ^x a u a^z a + 2u ak C > 0} Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 22 / 28
23 Computational Results Transmission Switching Standard IEEE Benchmark Instances Optimal Switching can make some difference in generation cost Generation Cost Instance No Switching With Switching case3les case6ww case case case ieee case case case case118b case case Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 23 / 28
24 Computational Results Transmission Switching But These Are (Now) Too Easy no cuts with cuts instance time nodes time nodes case3les case6ww case case case ieee case case case case118b case case Solving the MIP model using CPLEX v12.5 Case118 was the instance that Ferris et al. report not being able to solve with CPLEX (version 9) Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 24 / 28
25 Computational Results Transmission Switching Creating More Instances Modify the 118B instance by modifying the demands randomly. Create 15 new instances Comparing on Many 118B Instances Avg. Time Avg. Nodes No Cuts With Cuts Cuts show some promise We continue to work on pure transmission switching on larger instances (> 2000) nodes. These problems are still way too hard for CPLEX with and without cuts There are many alternative optimal solutions to the linear programming relaxation Which one(s)? should we cut off? Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 25 / 28
26 The End Accomplishments Prove that transmission switching problem is NP-Complete Understand a cycle relaxation derived from the structure of the problem Give a complete description of the convex hull of the set with 2 n inequalities Also have an (integer) extended formulation in dimension 6n + 1 Even with initial implementation, we can significantly improve default CPLEX behavior Still To Do Working on effective mechanisms for using these inequalities for larger instances A special challenge for (pure) transmission switching is the extreme dual generacy of LP solutions so engineering effective cutting plane mechanisms is important Up Next Study more complicated structures besides cycles Try to include demands at nodes, for Flow-(sci) Extend to potential preserved, but nonlinear relationship between potential and flow Gas and Water Network design Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 26 / 28
27 The End The Real Conclusion Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 27 / 28
28 The End Thank you! For Giving Us a Great Foundation Many Authors (UW-Madison & Ga. Tech) IP for Transmission Switching Aussois 28 / 28
Valid Inequalities for Potential-Constrained Network Design
Valid Inequalities for Potential-Constrained Network Design Hyemin Jeon Jeff Linderoth Jim Luedtke Dept. of ISyE UW-Madison jrluedt1@wisc.edu MINLP 2014, CMU, Pittsburgh, PA June 3, 2014 Jeon, Linderoth,
More informationA Cycle-Based Formulation and Valid Inequalities for DC Power Transmission Problems with Switching
A Cycle-Based Formulation and Valid Inequalities for DC Power Transmission Problems with Switching arxiv:1412.6245v1 [math.oc] 19 Dec 2014 Burak Kocuk Hyemin Jeon Santanu S. Dey Jeff Linderoth James Luedtke
More informationComputational Integer Programming. Lecture 2: Modeling and Formulation. Dr. Ted Ralphs
Computational Integer Programming Lecture 2: Modeling and Formulation Dr. Ted Ralphs Computational MILP Lecture 2 1 Reading for This Lecture N&W Sections I.1.1-I.1.6 Wolsey Chapter 1 CCZ Chapter 2 Computational
More informationStochastic Unit Commitment with Topology Control Recourse for Renewables Integration
1 Stochastic Unit Commitment with Topology Control Recourse for Renewables Integration Jiaying Shi and Shmuel Oren University of California, Berkeley IPAM, January 2016 33% RPS - Cumulative expected VERs
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective Vehicle routing problem 1 / 33 Cutting stock
More informationGLOBAL OPTIMIZATION METHODS FOR OPTIMAL POWER FLOW AND TRANSMISSION SWITCHING PROBLEMS IN ELECTRIC POWER SYSTEMS
GLOBAL OPTIMIZATION METHODS FOR OPTIMAL POWER FLOW AND TRANSMISSION SWITCHING PROBLEMS IN ELECTRIC POWER SYSTEMS A Thesis Presented to The Academic Faculty by Burak Kocuk In Partial Fulfillment of the
More informationIntroduction to Integer Programming
Lecture 3/3/2006 p. /27 Introduction to Integer Programming Leo Liberti LIX, École Polytechnique liberti@lix.polytechnique.fr Lecture 3/3/2006 p. 2/27 Contents IP formulations and examples Total unimodularity
More informationOn the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times
CORE DISCUSSION PAPER 2000/52 On the Polyhedral Structure of a Multi Item Production Planning Model with Setup Times Andrew J. Miller 1, George L. Nemhauser 2, and Martin W.P. Savelsbergh 2 November 2000
More informationIntroduction to Bin Packing Problems
Introduction to Bin Packing Problems Fabio Furini March 13, 2015 Outline Origins and applications Applications: Definition: Bin Packing Problem (BPP) Solution techniques for the BPP Heuristic Algorithms
More informationIntroduction to Mathematical Programming IE406. Lecture 21. Dr. Ted Ralphs
Introduction to Mathematical Programming IE406 Lecture 21 Dr. Ted Ralphs IE406 Lecture 21 1 Reading for This Lecture Bertsimas Sections 10.2, 10.3, 11.1, 11.2 IE406 Lecture 21 2 Branch and Bound Branch
More informationIntroduction to Integer Linear Programming
Lecture 7/12/2006 p. 1/30 Introduction to Integer Linear Programming Leo Liberti, Ruslan Sadykov LIX, École Polytechnique liberti@lix.polytechnique.fr sadykov@lix.polytechnique.fr Lecture 7/12/2006 p.
More informationInteger Programming ISE 418. Lecture 8. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 8 Dr. Ted Ralphs ISE 418 Lecture 8 1 Reading for This Lecture Wolsey Chapter 2 Nemhauser and Wolsey Sections II.3.1, II.3.6, II.4.1, II.4.2, II.5.4 Duality for Mixed-Integer
More informationa 1 a 2 a 3 a 4 v i c i c(a 1, a 3 ) = 3
AM 221: Advanced Optimization Spring 2016 Prof. Yaron Singer Lecture 17 March 30th 1 Overview In the previous lecture, we saw examples of combinatorial problems: the Maximal Matching problem and the Minimum
More informationMVE165/MMG630, Applied Optimization Lecture 6 Integer linear programming: models and applications; complexity. Ann-Brith Strömberg
MVE165/MMG630, Integer linear programming: models and applications; complexity Ann-Brith Strömberg 2011 04 01 Modelling with integer variables (Ch. 13.1) Variables Linear programming (LP) uses continuous
More informationOn the knapsack closure of 0-1 Integer Linear Programs. Matteo Fischetti University of Padova, Italy
On the knapsack closure of 0-1 Integer Linear Programs Matteo Fischetti University of Padova, Italy matteo.fischetti@unipd.it Andrea Lodi University of Bologna, Italy alodi@deis.unibo.it Aussois, January
More informationSolving the MWT. Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP:
Solving the MWT Recall the ILP for the MWT. We can obtain a solution to the MWT problem by solving the following ILP: max subject to e i E ω i x i e i C E x i {0, 1} x i C E 1 for all critical mixed cycles
More informationSeparation Techniques for Constrained Nonlinear 0 1 Programming
Separation Techniques for Constrained Nonlinear 0 1 Programming Christoph Buchheim Computer Science Department, University of Cologne and DEIS, University of Bologna MIP 2008, Columbia University, New
More information- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs
LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs
More informationLectures 6, 7 and part of 8
Lectures 6, 7 and part of 8 Uriel Feige April 26, May 3, May 10, 2015 1 Linear programming duality 1.1 The diet problem revisited Recall the diet problem from Lecture 1. There are n foods, m nutrients,
More informationWeek 4. (1) 0 f ij u ij.
Week 4 1 Network Flow Chapter 7 of the book is about optimisation problems on networks. Section 7.1 gives a quick introduction to the definitions of graph theory. In fact I hope these are already known
More information16.1 Min-Cut as an LP
600.469 / 600.669 Approximation Algorithms Lecturer: Michael Dinitz Topic: LPs as Metrics: Min Cut and Multiway Cut Date: 4//5 Scribe: Gabriel Kaptchuk 6. Min-Cut as an LP We recall the basic definition
More informationBBM402-Lecture 20: LP Duality
BBM402-Lecture 20: LP Duality Lecturer: Lale Özkahya Resources for the presentation: https://courses.engr.illinois.edu/cs473/fa2016/lectures.html An easy LP? which is compact form for max cx subject to
More informationON MIXING SETS ARISING IN CHANCE-CONSTRAINED PROGRAMMING
ON MIXING SETS ARISING IN CHANCE-CONSTRAINED PROGRAMMING Abstract. The mixing set with a knapsack constraint arises in deterministic equivalent of chance-constrained programming problems with finite discrete
More information3.4 Relaxations and bounds
3.4 Relaxations and bounds Consider a generic Discrete Optimization problem z = min{c(x) : x X} with an optimal solution x X. In general, the algorithms generate not only a decreasing sequence of upper
More information15.081J/6.251J Introduction to Mathematical Programming. Lecture 24: Discrete Optimization
15.081J/6.251J Introduction to Mathematical Programming Lecture 24: Discrete Optimization 1 Outline Modeling with integer variables Slide 1 What is a good formulation? Theme: The Power of Formulations
More informationSolving Box-Constrained Nonconvex Quadratic Programs
Noname manuscript No. (will be inserted by the editor) Solving Box-Constrained Nonconvex Quadratic Programs Pierre Bonami Oktay Günlük Jeff Linderoth June 13, 2016 Abstract We present effective computational
More informationFixed-charge transportation problems on trees
Fixed-charge transportation problems on trees Gustavo Angulo * Mathieu Van Vyve * gustavo.angulo@uclouvain.be mathieu.vanvyve@uclouvain.be November 23, 2015 Abstract We consider a class of fixed-charge
More informationLagrangian Relaxation in MIP
Lagrangian Relaxation in MIP Bernard Gendron May 28, 2016 Master Class on Decomposition, CPAIOR2016, Banff, Canada CIRRELT and Département d informatique et de recherche opérationnelle, Université de Montréal,
More informationCutting Planes in SCIP
Cutting Planes in SCIP Kati Wolter Zuse-Institute Berlin Department Optimization Berlin, 6th June 2007 Outline 1 Cutting Planes in SCIP 2 Cutting Planes for the 0-1 Knapsack Problem 2.1 Cover Cuts 2.2
More informationExtended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications
Extended Formulations, Lagrangian Relaxation, & Column Generation: tackling large scale applications François Vanderbeck University of Bordeaux INRIA Bordeaux-Sud-Ouest part : Defining Extended Formulations
More informationAn Integer Programming Approach for Linear Programs with Probabilistic Constraints
An Integer Programming Approach for Linear Programs with Probabilistic Constraints James Luedtke Shabbir Ahmed George Nemhauser Georgia Institute of Technology 765 Ferst Drive, Atlanta, GA, USA luedtke@gatech.edu
More informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
More informationBranch-and-Price-and-Cut for the Split Delivery Vehicle Routing Problem with Time Windows
Branch-and-Price-and-Cut for the Split Delivery Vehicle Routing Problem with Time Windows Guy Desaulniers École Polytechnique de Montréal and GERAD Column Generation 2008 Aussois, France Outline Introduction
More information3.7 Strong valid inequalities for structured ILP problems
3.7 Strong valid inequalities for structured ILP problems By studying the problem structure, we can derive strong valid inequalities yielding better approximations of conv(x ) and hence tighter bounds.
More informationLecture 8: Column Generation
Lecture 8: Column Generation (3 units) Outline Cutting stock problem Classical IP formulation Set covering formulation Column generation A dual perspective 1 / 24 Cutting stock problem 2 / 24 Problem description
More informationOn a Cardinality-Constrained Transportation Problem With Market Choice
On a Cardinality-Constrained Transportation Problem With Market Choice Pelin Damcı-Kurt Santanu S. Dey Simge Küçükyavuz Abstract It is well-known that the intersection of the matching polytope with a cardinality
More informationLecture 2. MATH3220 Operations Research and Logistics Jan. 8, Pan Li The Chinese University of Hong Kong. Integer Programming Formulations
Lecture 2 MATH3220 Operations Research and Logistics Jan. 8, 2015 Pan Li The Chinese University of Hong Kong 2.1 Agenda 1 2 3 2.2 : a linear program plus the additional constraints that some or all of
More informationWeek Cuts, Branch & Bound, and Lagrangean Relaxation
Week 11 1 Integer Linear Programming This week we will discuss solution methods for solving integer linear programming problems. I will skip the part on complexity theory, Section 11.8, although this is
More information1 Column Generation and the Cutting Stock Problem
1 Column Generation and the Cutting Stock Problem In the linear programming approach to the traveling salesman problem we used the cutting plane approach. The cutting plane approach is appropriate when
More informationMixed Integer Linear Programming Formulations for Probabilistic Constraints
Mixed Integer Linear Programming Formulations for Probabilistic Constraints J. P. Vielma a,, S. Ahmed b, G. Nemhauser b a Department of Industrial Engineering, University of Pittsburgh 1048 Benedum Hall,
More informationStructure of Valid Inequalities for Mixed Integer Conic Programs
Structure of Valid Inequalities for Mixed Integer Conic Programs Fatma Kılınç-Karzan Tepper School of Business Carnegie Mellon University 18 th Combinatorial Optimization Workshop Aussois, France January
More informationInteger Equal Flows. 1 Introduction. Carol A. Meyers Andreas S. Schulz
Integer Equal Flows Carol A Meyers Andreas S Schulz Abstract We examine an NP-hard generalization of the network flow problem known as the integer equal flow problem The setup is the same as a standard
More informationA new primal-dual framework for European day-ahead electricity auctions
A new primal-dual framework for European day-ahead electricity auctions Mehdi Madani*, Mathieu Van Vyve** *Louvain School of management, ** CORE Catholic University of Louvain Mathematical Models and Methods
More informationMixed-Integer Nonlinear Programming
Mixed-Integer Nonlinear Programming Claudia D Ambrosio CNRS researcher LIX, École Polytechnique, France pictures taken from slides by Leo Liberti MPRO PMA 2016-2017 Motivating Applications Nonlinear Knapsack
More informationNetwork Flows. 6. Lagrangian Relaxation. Programming. Fall 2010 Instructor: Dr. Masoud Yaghini
In the name of God Network Flows 6. Lagrangian Relaxation 6.3 Lagrangian Relaxation and Integer Programming Fall 2010 Instructor: Dr. Masoud Yaghini Integer Programming Outline Branch-and-Bound Technique
More informationAggregation and Mixed Integer Rounding to Solve MILPs (Marchand and Wolsey)
Aggregation and Mixed Integer Rounding to Solve MILPs (Marchand and Wolsey) SAS Institute - Analytical Solutions Lehigh University - Department of Industrial and Systems Engineering July 7, 2005 Classical
More informationAn Integrated Approach to Truss Structure Design
Slide 1 An Integrated Approach to Truss Structure Design J. N. Hooker Tallys Yunes CPAIOR Workshop on Hybrid Methods for Nonlinear Combinatorial Problems Bologna, June 2010 How to Solve Nonlinear Combinatorial
More informationInteger Programming ISE 418. Lecture 12. Dr. Ted Ralphs
Integer Programming ISE 418 Lecture 12 Dr. Ted Ralphs ISE 418 Lecture 12 1 Reading for This Lecture Nemhauser and Wolsey Sections II.2.1 Wolsey Chapter 9 ISE 418 Lecture 12 2 Generating Stronger Valid
More informationPreliminaries Overview OPF and Extensions. Convex Optimization. Lecture 8 - Applications in Smart Grids. Instructor: Yuanzhang Xiao
Convex Optimization Lecture 8 - Applications in Smart Grids Instructor: Yuanzhang Xiao University of Hawaii at Manoa Fall 2017 1 / 32 Today s Lecture 1 Generalized Inequalities and Semidefinite Programming
More informationBranch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems
Branch-and-cut Approaches for Chance-constrained Formulations of Reliable Network Design Problems Yongjia Song James R. Luedtke August 9, 2012 Abstract We study solution approaches for the design of reliably
More informationLecture 5 January 16, 2013
UBC CPSC 536N: Sparse Approximations Winter 2013 Prof. Nick Harvey Lecture 5 January 16, 2013 Scribe: Samira Samadi 1 Combinatorial IPs 1.1 Mathematical programs { min c Linear Program (LP): T x s.t. a
More informationTechnische Universität München, Zentrum Mathematik Lehrstuhl für Angewandte Geometrie und Diskrete Mathematik. Combinatorial Optimization (MA 4502)
Technische Universität München, Zentrum Mathematik Lehrstuhl für Angewandte Geometrie und Diskrete Mathematik Combinatorial Optimization (MA 4502) Dr. Michael Ritter Problem Sheet 1 Homework Problems Exercise
More informationCutting Plane Separators in SCIP
Cutting Plane Separators in SCIP Kati Wolter Zuse Institute Berlin DFG Research Center MATHEON Mathematics for key technologies 1 / 36 General Cutting Plane Method MIP min{c T x : x X MIP }, X MIP := {x
More informationUsing Sparsity to Design Primal Heuristics for MILPs: Two Stories
for MILPs: Two Stories Santanu S. Dey Joint work with: Andres Iroume, Marco Molinaro, Domenico Salvagnin, Qianyi Wang MIP Workshop, 2017 Sparsity in real" Integer Programs (IPs) Real" IPs are sparse: The
More informationOptimization Methods in Management Science
Optimization Methods in Management Science MIT 15.05 Recitation 8 TAs: Giacomo Nannicini, Ebrahim Nasrabadi At the end of this recitation, students should be able to: 1. Derive Gomory cut from fractional
More informationLecture 11 October 7, 2013
CS 4: Advanced Algorithms Fall 03 Prof. Jelani Nelson Lecture October 7, 03 Scribe: David Ding Overview In the last lecture we talked about set cover: Sets S,..., S m {,..., n}. S has cost c S. Goal: Cover
More informationReconnect 04 Introduction to Integer Programming
Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, Reconnect 04 Introduction to Integer Programming Cynthia Phillips, Sandia National Laboratories Integer programming
More informationOptimization Exercise Set n.5 :
Optimization Exercise Set n.5 : Prepared by S. Coniglio translated by O. Jabali 2016/2017 1 5.1 Airport location In air transportation, usually there is not a direct connection between every pair of airports.
More informationInteger Programming Methods LNMB
Integer Programming Methods LNMB 2017 2018 Dion Gijswijt homepage.tudelft.nl/64a8q/intpm/ Dion Gijswijt Intro IntPM 2017-2018 1 / 24 Organisation Webpage: homepage.tudelft.nl/64a8q/intpm/ Book: Integer
More informationOn mathematical programming with indicator constraints
On mathematical programming with indicator constraints Andrea Lodi joint work with P. Bonami & A. Tramontani (IBM), S. Wiese (Unibo) University of Bologna, Italy École Polytechnique de Montréal, Québec,
More informationOrbital Conflict. Jeff Linderoth. Jim Ostrowski. Fabrizio Rossi Stefano Smriglio. When Worlds Collide. Univ. of Wisconsin-Madison
Orbital Conflict When Worlds Collide Jeff Linderoth Univ. of Wisconsin-Madison Jim Ostrowski University of Tennessee Fabrizio Rossi Stefano Smriglio Univ. of L Aquila MIP 2014 Columbus, OH July 23, 2014
More informationCombinatorial Optimization
Combinatorial Optimization Lecture notes, WS 2010/11, TU Munich Prof. Dr. Raymond Hemmecke Version of February 9, 2011 Contents 1 The knapsack problem 1 1.1 Complete enumeration..................................
More informationLinear and Integer Programming - ideas
Linear and Integer Programming - ideas Paweł Zieliński Institute of Mathematics and Computer Science, Wrocław University of Technology, Poland http://www.im.pwr.wroc.pl/ pziel/ Toulouse, France 2012 Literature
More informationOptimization Exercise Set n. 4 :
Optimization Exercise Set n. 4 : Prepared by S. Coniglio and E. Amaldi translated by O. Jabali 2018/2019 1 4.1 Airport location In air transportation, usually there is not a direct connection between every
More informationSoftware for Integer and Nonlinear Optimization
Software for Integer and Nonlinear Optimization Sven Leyffer, leyffer@mcs.anl.gov Mathematics & Computer Science Division Argonne National Laboratory Roger Fletcher & Jeff Linderoth Advanced Methods and
More informationComputer Sciences Department
Computer Sciences Department Valid Inequalities for the Pooling Problem with Binary Variables Claudia D Ambrosio Jeff Linderoth James Luedtke Technical Report #682 November 200 Valid Inequalities for the
More informationSection Notes 9. Midterm 2 Review. Applied Math / Engineering Sciences 121. Week of December 3, 2018
Section Notes 9 Midterm 2 Review Applied Math / Engineering Sciences 121 Week of December 3, 2018 The following list of topics is an overview of the material that was covered in the lectures and sections
More informationStructure in Mixed Integer Conic Optimization: From Minimal Inequalities to Conic Disjunctive Cuts
Structure in Mixed Integer Conic Optimization: From Minimal Inequalities to Conic Disjunctive Cuts Fatma Kılınç-Karzan Tepper School of Business Carnegie Mellon University Joint work with Sercan Yıldız
More informationInteger Linear Programs
Lecture 2: Review, Linear Programming Relaxations Today we will talk about expressing combinatorial problems as mathematical programs, specifically Integer Linear Programs (ILPs). We then see what happens
More informationStrongly Polynomial Algorithm for a Class of Minimum-Cost Flow Problems with Separable Convex Objectives
Strongly Polynomial Algorithm for a Class of Minimum-Cost Flow Problems with Separable Convex Objectives László A. Végh April 12, 2013 Abstract A well-studied nonlinear extension of the minimum-cost flow
More informationTutorial 2: Modelling Transmission
Tutorial 2: Modelling Transmission In our previous example the load and generation were at the same bus. In this tutorial we will see how to model the transmission of power from one bus to another. The
More informationCO 250 Final Exam Guide
Spring 2017 CO 250 Final Exam Guide TABLE OF CONTENTS richardwu.ca CO 250 Final Exam Guide Introduction to Optimization Kanstantsin Pashkovich Spring 2017 University of Waterloo Last Revision: March 4,
More informationReformulation and Sampling to Solve a Stochastic Network Interdiction Problem
Network Interdiction Stochastic Network Interdiction and to Solve a Stochastic Network Interdiction Problem Udom Janjarassuk Jeff Linderoth ISE Department COR@L Lab Lehigh University jtl3@lehigh.edu informs
More informationStochastic Decision Diagrams
Stochastic Decision Diagrams John Hooker CORS/INFORMS Montréal June 2015 Objective Relaxed decision diagrams provide an generalpurpose method for discrete optimization. When the problem has a dynamic programming
More informationA MIXED INTEGER DISJUNCTIVE MODEL FOR TRANSMISSION NETWORK EXPANSION
A MIXED INTEGER DISJUNCTIVE MODEL FOR TRANSMISSION NETWORK EXPANSION Laura Bahiense*, Gerson C. Oliveira (PUC/Rio), Mario Pereira*, Member, Sergio Granville*, Member Abstract: The classical non-linear
More informationIndicator Constraints in Mixed-Integer Programming
Indicator Constraints in Mixed-Integer Programming Andrea Lodi University of Bologna, Italy - andrea.lodi@unibo.it Amaya Nogales-Gómez, Universidad de Sevilla, Spain Pietro Belotti, FICO, UK Matteo Fischetti,
More informationPower Grid Partitioning: Static and Dynamic Approaches
Power Grid Partitioning: Static and Dynamic Approaches Miao Zhang, Zhixin Miao, Lingling Fan Department of Electrical Engineering University of South Florida Tampa FL 3320 miaozhang@mail.usf.edu zmiao,
More informationON THE INTEGRALITY OF THE UNCAPACITATED FACILITY LOCATION POLYTOPE. 1. Introduction
ON THE INTEGRALITY OF THE UNCAPACITATED FACILITY LOCATION POLYTOPE MOURAD BAÏOU AND FRANCISCO BARAHONA Abstract We study a system of linear inequalities associated with the uncapacitated facility location
More informationLecture 9: Dantzig-Wolfe Decomposition
Lecture 9: Dantzig-Wolfe Decomposition (3 units) Outline Dantzig-Wolfe decomposition Column generation algorithm Relation to Lagrangian dual Branch-and-price method Generated assignment problem and multi-commodity
More informationInteger Programming Duality
Integer Programming Duality M. Guzelsoy T. K. Ralphs July, 2010 1 Introduction This article describes what is known about duality for integer programs. It is perhaps surprising that many of the results
More informationThe Fixed Charge Transportation Problem: A Strong Formulation Based On Lagrangian Decomposition and Column Generation
The Fixed Charge Transportation Problem: A Strong Formulation Based On Lagrangian Decomposition and Column Generation Yixin Zhao, Torbjörn Larsson and Department of Mathematics, Linköping University, Sweden
More informationNP-problems continued
NP-problems continued Page 1 Since SAT and INDEPENDENT SET can be reduced to each other we might think that there would be some similarities between the two problems. In fact, there is one such similarity.
More informationOutline. 1 Introduction. 2 Modeling and Applications. 3 Algorithms Convex. 4 Algorithms Nonconvex. 5 The Future of MINLP? 2 / 126
Outline 1 Introduction 2 3 Algorithms Convex 4 Algorithms Nonconvex 5 The Future of MINLP? 2 / 126 Larry Leading Off 3 / 126 Introduction Outline 1 Introduction Software Tools and Online Resources Basic
More informationA Lower Bound on the Split Rank of Intersection Cuts
A Lower Bound on the Split Rank of Intersection Cuts Santanu S. Dey H. Milton Stewart School of Industrial and Systems Engineering, Georgia Institute of Technology. 200 Outline Introduction: split rank,
More informationwith Binary Decision Diagrams Integer Programming J. N. Hooker Tarik Hadzic IT University of Copenhagen Carnegie Mellon University ICS 2007, January
Integer Programming with Binary Decision Diagrams Tarik Hadzic IT University of Copenhagen J. N. Hooker Carnegie Mellon University ICS 2007, January Integer Programming with BDDs Goal: Use binary decision
More informationComputational Integer Programming Universidad de los Andes. Lecture 1. Dr. Ted Ralphs
Computational Integer Programming Universidad de los Andes Lecture 1 Dr. Ted Ralphs MIP Lecture 1 1 Quick Introduction Bio Course web site Course structure http://coral.ie.lehigh.edu/ ted/teaching/mip
More informationDiscrete Optimization 2010 Lecture 7 Introduction to Integer Programming
Discrete Optimization 2010 Lecture 7 Introduction to Integer Programming Marc Uetz University of Twente m.uetz@utwente.nl Lecture 8: sheet 1 / 32 Marc Uetz Discrete Optimization Outline 1 Intro: The Matching
More informationTravelling Salesman Problem
Travelling Salesman Problem Fabio Furini November 10th, 2014 Travelling Salesman Problem 1 Outline 1 Traveling Salesman Problem Separation Travelling Salesman Problem 2 (Asymmetric) Traveling Salesman
More informationBranching Rules for Minimum Congestion Multi- Commodity Flow Problems
Clemson University TigerPrints All Theses Theses 8-2012 Branching Rules for Minimum Congestion Multi- Commodity Flow Problems Cameron Megaw Clemson University, cmegaw@clemson.edu Follow this and additional
More informationTopics in Theoretical Computer Science April 08, Lecture 8
Topics in Theoretical Computer Science April 08, 204 Lecture 8 Lecturer: Ola Svensson Scribes: David Leydier and Samuel Grütter Introduction In this lecture we will introduce Linear Programming. It was
More informationDuality of LPs and Applications
Lecture 6 Duality of LPs and Applications Last lecture we introduced duality of linear programs. We saw how to form duals, and proved both the weak and strong duality theorems. In this lecture we will
More informationSubset selection with sparse matrices
Subset selection with sparse matrices Alberto Del Pia, University of Wisconsin-Madison Santanu S. Dey, Georgia Tech Robert Weismantel, ETH Zürich February 1, 018 Schloss Dagstuhl Subset selection for regression
More informationStochastic Integer Programming
IE 495 Lecture 20 Stochastic Integer Programming Prof. Jeff Linderoth April 14, 2003 April 14, 2002 Stochastic Programming Lecture 20 Slide 1 Outline Stochastic Integer Programming Integer LShaped Method
More informationCS 6783 (Applied Algorithms) Lecture 3
CS 6783 (Applied Algorithms) Lecture 3 Antonina Kolokolova January 14, 2013 1 Representative problems: brief overview of the course In this lecture we will look at several problems which, although look
More informationCS 6820 Fall 2014 Lectures, October 3-20, 2014
Analysis of Algorithms Linear Programming Notes CS 6820 Fall 2014 Lectures, October 3-20, 2014 1 Linear programming The linear programming (LP) problem is the following optimization problem. We are given
More informationMIP reformulations of some chance-constrained mathematical programs
MIP reformulations of some chance-constrained mathematical programs Ricardo Fukasawa Department of Combinatorics & Optimization University of Waterloo December 4th, 2012 FIELDS Industrial Optimization
More informationThe Pooling Problem - Applications, Complexity and Solution Methods
The Pooling Problem - Applications, Complexity and Solution Methods Dag Haugland Dept. of Informatics University of Bergen Norway EECS Seminar on Optimization and Computing NTNU, September 30, 2014 Outline
More informationSection Notes 8. Integer Programming II. Applied Math 121. Week of April 5, expand your knowledge of big M s and logical constraints.
Section Notes 8 Integer Programming II Applied Math 121 Week of April 5, 2010 Goals for the week understand IP relaxations be able to determine the relative strength of formulations understand the branch
More informationLecture 13 March 7, 2017
CS 224: Advanced Algorithms Spring 2017 Prof. Jelani Nelson Lecture 13 March 7, 2017 Scribe: Hongyao Ma Today PTAS/FPTAS/FPRAS examples PTAS: knapsack FPTAS: knapsack FPRAS: DNF counting Approximation
More information