A Tabu Search Approach for Delivering Pet Food and Flour in Switzerland

In this paper, we consider a real-life vehicle routeing problem that occurs in a major Swiss company producing pet food and flour. In contrast with usual hypothetical problems, a large variety of restrictions has to be considered. The main constraints are relative to the accessibility and the time windows at customers, the carrying capacities of vehicles, the total duration of routes and the drivers' breaks. To find good solutions to this problem, we propose two heuristic methods: a fast straightforward insertion procedure and a method based on tabu search techniques. Next, the produced solutions are compared with the routes actually covered by the company. Our outcomes indicate that the total distance travelled can be reduced significantly when such methods are used.     

Exact and heuristic methodologies for scheduling in hospitals: problems, formulations and algorithms

This text summarizes the PhD thesis defended by the author in January 2006 under the supervision of Professor Erik Demeulemeester at the Katholieke Universiteit Leuven. The thesis is written in English and is available from the author’s website (http://www.econ.kuleuven.be/jeroen.belien). In this research we propose a number of exact and heuristic algorithms for various scheduling problems encountered in hospitals. The emphasis lies on the design of new methodologies as well as on the applicability of the algorithms in real-life environments. The main contributions include a new decomposition approach for a particular class of staff scheduling problems, an extensive study of master surgery scheduling algorithms that aim at leveling the resultant bed occupancy and an innovative method for integrating nurse and surgery scheduling.      

Optimal portfolio selection for general provisioning and terminal wealth problems

In Dhaene et al. (2005), multiperiod portfolio selection problems are discussed, using an analytical approach to find optimal constant mix investment strategies in a provisioning or a savings context. In this paper we extend some of these results, investigating some specific, real-life situations. The problems that we consider in the first section of this paper are general in the sense that they allow for liabilities that can be both positive or negative, as opposed to Dhaene et al. (2005), where all liabilities have to be of the same sign. Secondly, we generalize portfolio selection problems to the case where a minimal return requirement is imposed. We derive an intuitive formula that can be used in provisioning and terminal wealth problems as a constraint on the admissible investment portfolios, in order to guarantee a minimal annualized return. We apply our results to optimal portfolio selection.     

Optimal fleet composition and periodic routing of offshore supply vessels

The supply vessel planning problem is a maritime transportation problem faced by amongst others the energy company Statoil. A set of offshore installations requires supplies from an onshore supply depot on a regular basis, a service performed by a fleet of offshore supply vessels. The problem consists of determining the optimal fleet composition of offshore supply vessels and their corresponding weekly voyages and schedules. We present a voyage-based solution method for the supply vessel planning problem. A computational study shows how the solution method can be used to solve real-life problems. Statoil has implemented a planning tool based on the voyage-based solution method and reports significant savings.     

Supply chain design considering economies of scale and transport frequencies

In this paper we consider a 3-echelon, multi-product supply chain design model with economies of scale in transport and warehousing that explicitly takes transport frequencies into consideration. Our model simultaneously optimizes locations and sizes of tank farms, material flows, and transport frequencies within the network. We consider all relevant costs: product cost, transport cost, tank rental cost, tank throughput cost, and inventory cost. The problem is based on a real-life example from a chemical company. We show that considering economies of scale and transport frequencies in the design stage is crucial and failing to do so can lead to substantially higher costs than optimal. We solve a wide variety of problems with branch-and-bound and with the efficient solution heuristics based on iterative linearization techniques we develop. We show that the heuristics are superior to the standard branch-and-bound technique for large problems like the one of the chemical company that motivated our research.     

A Posteriori Error Estimates for Distributed Convex Optimal Control Problems

In this paper, we present an a posteriori error analysis for finite element approximation of distributed convex optimal control problems. We derive a posteriori error estimates for the coupled state and control approximations under some assumptions which hold in many applications. Such estimates, which are apparently not available in the literature, can be used to construct reliable adaptive finite element approximation schemes for control problems. Explicit estimates are obtained for some model problems which frequently appear in real-life applications.     

Weak solutions to a nonlinear variational wave equation and some related problems

In this paper we present some results on the global existence of weak solutions to a nonlinear variational wave equation and some related problems. We first introduce the main tools, the L ^p Young measure theory and related compactness results, in the first section. Then we use the L ^p Young measure theory to prove the global existence of dissipative weak solutions to the asymptotic equation of the nonlinear wave equation, and comment on its relation to Camassa-Holm equations in the second section. In the third section, we prove the global existence of weak solutions to the original nonlinear wave equation under some restrictions on the wave speed. In the last section, we present global existence of renormalized solutions to two-dimensional model equations of the asymptotic equation, which is also the so-called vortex density equation arising from sup-conductivity.     

A Trust-Region Method for Unconstrained Multiobjective Problems with Applications in Satisficing Processes

Multiobjective optimization has a significant number of real-life applications. For this reason, in this paper we consider the problem of finding Pareto critical points for unconstrained multiobjective problems and present a trust-region method to solve it. Under certain assumptions, which are derived in a very natural way from assumptions used to establish convergence results of the scalar trust-region method, we prove that our trust-region method generates a sequence which converges in the Pareto critical way. This means that our generalized marginal function, which generalizes the norm of the gradient for the multiobjective case, converges to zero. In the last section of this paper, we give an application to satisficing processes in Behavioral Sciences. Multiobjective trust-region methods appear to be remarkable specimens of much more abstract satisficing processes, based on “variational rationality” concepts. One of their important merits is to allow for efficient computations. This is a striking result in Behavioral Sciences.     

A flexible adaptive memory-based algorithm for real-life transportation operations: Two case studies from dairy and construction sector

Effective routing of vehicles remains a focal goal of all modern enterprises, thriving for excellence in project management with minimal investment and operational costs. This paper proposes a metaheuristic methodology for solving a practical variant of the well-known Vehicle Routing Problem, called Heterogeneous Fixed Fleet VRP (HFFVRP). Using a two-phase construction heuristic, called GEneralized ROute Construction Algorithm (GEROCA), the proposed metaheuristic approach enhances its flexibility to easily adopt various operational constraints. Via this approach, two real-life distribution problems faced by a dairy and a construction company were tackled and formulated as HFFVRP. Computational results on the aforementioned case studies show that the proposed metaheuristic approach (a) consistently outperforms previous published metaheuristic approaches we have developed to solve the HFFVRP, and (b) substantially improves upon the current practice of the company. The key result that impressed both companies’ management was the improvement over the bi-objective character of their problems: the minimization of the total distribution cost as well as the minimization of the number of the given heterogeneous number of vehicles used.     

A survey for the quadratic assignment problem

The quadratic assignment problem (QAP), one of the most difficult problems in the NP-hard class, models many real-life problems in several areas such as facilities location, parallel and distributed computing, and combinatorial data analysis. Combinatorial optimization problems, such as the traveling salesman problem, maximal clique and graph partitioning can be formulated as a QAP. In this paper, we present some of the most important QAP formulations and classify them according to their mathematical sources. We also present a discussion on the theoretical resources used to define lower bounds for exact and heuristic algorithms. We then give a detailed discussion of the progress made in both exact and heuristic solution methods, including those formulated according to metaheuristic strategies. Finally, we analyze the contributions brought about by the study of different approaches.     

Improved Steplength by More Practical Information in the Extragradient Method for Monotone Variational Inequalities

The extragradient method uses only the values of the function in the variational inequality to solve it. In many real-life problems, the functions do not have explicit expressions and their evaluations are expensive. Therefore, it is important and practical to reduce the number of function evaluations used in these problems. In order to do this, we present some modifications to one existing extragradient method. First, we analyze the limit of the projection along a direction onto a closed convex set. Then, using the obtained result, we give a horizontal asymptotic property of an estimated function to the one measuring the progress in each iteration. From this property, a Newton search as well a self-adjusted relaxation procedure are introduced with an improved steplength into the extragradient method. Besides the theoretical background, numerical results are given to show the progress of the modified method. This work was supported by National Natural Science Foundation of China 10571083 and Doctoral Fund of Ministry of Education of China 20060284001.     

One Approach to the Comparison of Error Bounds at a Point and on a Set in the Solution of Ill-Posed Problems

The approximate solution of ill-posed problems by the regularization method always involves the issue of estimating the error. It is a common practice to use uniform bounds on the whole class of well-posedness in terms of the modulus of continuity of the inverse operator on this class. Local error bounds, which are also called error bounds at a point, have been studied much less. Since the solution of a real-life ill-posed problem is unique, an error bound obtained on the whole class of well-posedness roughens to a great extent the true error bound. In the present paper, we study the difference between error bounds on the class of well-posedness and error bounds at a point for a special class of ill-posed problems. Assuming that the exact solution is a piecewise smooth function, we prove that an error bound at a point is infinitely smaller than the exact bound on the class of well-posedness.     

Location-based techniques for the synergy approximation in combinatorial transportation auctions

The use of combinatorial auctions for the procurement of transportation services is investigated in this paper. We focus on the carrier viewpoint who is interested in submitting to the auction a selected bundle of loads that avoids some of the empty movement of the vehicles in his transportation network and that increases his profits. For this purpose we develop an optimization approach based on the use of the location techniques, usually used in the context of facility planning. Mathematically, this means maximizing the synergy among the bundle’s auctioned loads from one side and between the auctioned and the pre-existing loads from the other side. We show the validity of our approach by using first an illustrative example and then by applying it to solve a real-life problem related to a logistics company installed in the Arabic Gulf region.     

Optimal Control of Continuous Casting by Nondifferentiable Multiobjective Optimization

A new version of an interactive NIMBUS method for nondifferentiable multiobjective optimization is described. It is based on the reference point idea and the classification of the objective functions. The original problem is transformed into a single objective form according to the classification information. NIMBUS has been designed especially to be able to handle complicated real-life problems in a user-friendly way.The NIMBUS method is used for solving an optimal control problem related to the continuous casting of steel. The main goal is to minimize the defects in the final product. Conflicting objective functions are constructed according to certain metallurgical criteria and some technological constraints. Due to the phase changes during the cooling process there exist discontinuities in the derivative of the temperature distribution. Thus, the problem is nondifferentiable.Like many real-life problems, the casting model is large and complicated and numerically demanding. NIMBUS provides an efficient way of handling the difficulties and, at the same time, aids the user in finding a satisficing solution. In the end, some numerical experiments are reported and compared with earlier results.     

Wavelet methods for continuous-time prediction using Hilbert-valued autoregressive processes

We consider the prediction problem of a continuous-time stochastic process on an entire time-interval in terms of its recent past. The approach we adopt is based on the notion of autoregressive Hilbert processes that represent a generalization of the classical autoregressive processes to random variables with values in a Hilbert space. A careful analysis reveals, in particular, that this approach is related to the theory of function estimation in linear ill-posed inverse problems. In the deterministic literature, such problems are usually solved by suitable regularization techniques. We describe some recent approaches from the deterministic literature that can be adapted to obtain fast and feasible predictions. For large sample sizes, however, these approaches are not computationally efficient.With this in mind, we propose three linear wavelet methods to efficiently address the aforementioned prediction problem. We present regularization techniques for the sample paths of the stochastic process and obtain consistency results of the resulting prediction estimators. We illustrate the performance of the proposed methods in finite sample situations by means of a real-life data example which concerns with the prediction of the entire annual cycle of climatological El Niño-Southern Oscillation time series 1 year ahead. We also compare the resulting predictions with those obtained by other methods available in the literature, in particular with a smoothing spline interpolation method and with a SARIMA model.     

Algorithm engineering for optimal alignment of protein structure distance matrices

Protein structural alignment is an important problem in computational biology. In this paper, we present first successes on provably optimal pairwise alignment of protein inter-residue distance matrices, using the popular dali scoring function. We introduce the structural alignment problem formally, which enables us to express a variety of scoring functions used in previous work as special cases in a unified framework. Further, we propose the first mathematical model for computing optimal structural alignments based on dense inter-residue distance matrices. We therefore reformulate the problem as a special graph problem and give a tight integer linear programming model. We then present algorithm engineering techniques to handle the huge integer linear programs of real-life distance matrix alignment problems. Applying these techniques, we can compute provably optimal dali alignments for the very first time.     

Optimal Value and Post Optimal Solution in a Transportation Problem

In this work, we analyse the transportation problem of a real-life situation by obtaining the optimal feasible solutions, thus carrying out the sensitivity analysis of the problem. The work utilises the data obtained from the Asejire and Ikeja plants of Coca-Cola company, aiming to aid decision-making regarding the best possible options to satisfy customers at the barest minimum cost of transportation. Rerunning the optimization of a problem is an expensive scheme for gathering and obtaining enough data required for a problem. Thus, to minimize the transportation cost, the sensitivity analysis of parameters is a good tool to determine the behaviour of some input parameters where the values of these parameters are varied arbitrarily such that optimal results are verified. Maple 18 Software is used to solve the problem and the result obtained is compared with the values evaluated from northwest corner method, least cost method and Vogel''s approximation method. The study critically shows how a little change in a unit or more of any model parameter affects the expected results.     

A new two-point scheme for multiple roots of nonlinear equations

In most of the earlier research for multiple zeros, in order to obtain a new iteration function from the existing scheme, the usual practice is to make no change at the first substep. In this paper, we explore the idea that what are the advantages if the flexibility of choice is also given at the first substep. Therefore, we present a new two-point sixth-order scheme for multiple roots (m>1). The main advantages of our scheme over the existing schemes are flexibility at both substeps, simple body structure, smaller residual error, smaller error difference between two consecutive iterations, and smaller asymptotic error constant. The development of the scheme is based on midpoint formula and weight functions of two variables. We compare our methods with the existing methods of the same order with real-life applications as well as standard test problems. From the numerical results, we find that our methods can be considered as better alternates for the existing methods of the same order. Finally, dynamical study of the proposed schemes is presented that confirms the theoretical results.     

An algorithm for addressing the real interval eigenvalue problem

In this paper we present an algorithm for approximating the range of the real eigenvalues of interval matrices. Such matrices could be used to model real-life problems, where data sets suffer from bounded variations such as uncertainties (e.g. tolerances on parameters, measurement errors), or to study problems for given states.The algorithm that we propose is a subdivision algorithm that exploits sophisticated techniques from interval analysis. The quality of the computed approximation and the running time of the algorithm depend on a given input accuracy. We also present an efficient C++ implementation and illustrate its efficiency on various data sets. In most of the cases we manage to compute efficiently the exact boundary points (limited by floating point representation).     

Wall effects on the slow steady motion of a particle in a viscous incompressible fluid

In this paper, asymptotic expansions with respect to small Reynolds numbers are proved for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid past a single plane wall. The flow problem is modelled by a certain boundary value problem for the stationary, nonlinear Navier-Stokes equations. The coefficients of these expansions are the solutions of various, linear Stokes problems which can be constructed by single layer potentials and corresponding boundary integral equations on the boundary surface of the particle. Furthermore, some asymptotic estimates at small Reynolds numbers are presented for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid between two parallel, plane walls and in an infinitely long, rectilinear cylinder of arbitrary cross section. In the case of the flow problem with a single plane wall, the paper is based on a thesis, which the author has written under the guidance of Professor Dr. Wolfgang L. Wendland.