Many supply chain problems involve discrete decisions in a dynamic environment. The inventory routing problem is an example that combines the dynamic control of inventory at various facilities in a supply chain with the discrete routing decisions of a fleet of vehicles that moves product between the facilities.
We study these problems modeled as mixed-integer programs and propose a time decomposition based on approximate inventory valuation. We generate the approximate value function with an algorithm that combines data fitting, discrete optimization and dynamic programming methodology. Our framework allows the user to specify a class of piecewise linear, concave functions from which the algorithm chooses the value function. The use of piecewise linear concave functions is motivated by intuition, theory and practice. Intuitively, concavity reflects the notion that inventory is marginally more valuable the closer one is to a stock-out. Theoretically, piecewise linear concave functions have certain structural properties that also hold for finite mixed-integer program value functions. (Whether the same properties hold in the infinite case is an open question, to our knowledge.) Practically, piecewise linear concave functions are easily embedded in the objective function of a maximization mixed-integer or linear program, with only a few additional auxiliary continuous variables. We evaluate the solutions generated by our value functions in a case study using maritime inventory routing instances inspired by the petrochemical industry.
The thesis also includes two other contributions. First, we review various data fitting optimization models related to piecewise linear concave functions, and introduce new mixed-integer programming formulations for some cases. The formulations may be of independent interest, with applications in engineering, mixed-integer non-linear programming, and other areas. Second, we study a discounted, infinite-horizon version of the canonical single-item lot-sizing problem and characterize its value function, proving that it inherits all properties of interest from its finite counterpart. We then compare its optimal policies to our algorithm's solutions as a proof of concept.PhDCommittee Chair: George Nemhauser; Committee Member: Ahmet Keha; Committee Member: Martin Savelsbergh; Committee Member: Santanu Dey; Committee Member: Shabbir Ahme
