Dealing with various sources of uncertainty in the operational control of water systems using ensemble based MPC with convex optimization

Author(s):  
Klaudia Horvath ◽  
Maarten Smoorenburg ◽  
Diederik Vreeken ◽  
Ruben Sinnige ◽  
Rodolfo Alvarado Montero ◽  
...  

<p>Model Predictive Control (MPC) can be an effective tool for the operational control of water systems, but there are still many open questions about how this technique can effectively take into uncertainties of forecasts, initial states or the model setup. Moreover, computational cost and robustness often prohibit the use of existing methods in practice. We here report recent developments in the open source RTC-Tools software framework that allow representing these uncertainties through ensembles and computing the optimal control strategy with convex optimization techniques in combination with lexicographical goal programming. Convex optimization is required to have robust mathematical solutions within the short computation times that are feasible in operational practice. Goal programming is here used to facilitate straightforward optimization of competing objectives with results understandable for end-users. Adaptations of Raso’s Tree-Based MPC (e.g. Raso et al., 2014) are used to represent the possibilities offered in future control steps, permitting a realistic moving horizon control strategy while not being excessively conservative.</p><p>The developments are illustrated with applications in different water systems using methods for convex optimization of linear Mixed Integer problems as well as quadratically constrained problems with both open source and commercial solvers. We also demonstrate how RTC-Tools build-in methods can be used for linearization of system equations and objectives. The applications were evaluated in controlled experiments to learn about strengths and weaknesses in comparison with other ensemble and deterministic MPC methods.</p><p>Exploration of the added value of selected uncertainty representation techniques within MPC solutions is presented in a separate contribution (Smoorenburg et al. 2020, session HS4.3 “Ensemble hydrological forecasting: Decision making under uncertainty”).</p><p>Raso, L., D. Schwanenberg, N. C. van de Giesen, and P. J. van Overloop. 2014. “Short-Term Optimal Operation of Water Systems Using Ensemble Forecasts.” Advances in Water Resources 71 (September): 200–208.</p>

2019 ◽  
Vol 9 (24) ◽  
pp. 5524
Author(s):  
Dongju Cao ◽  
Wendong Yang ◽  
Gangyi Xu

With the rapid development and evolvement of unmanned aerial vehicle (UAV) technology, UAV aided wireless communication technology has been widely studied recently. In this paper, a buffer aided multi-UAV relaying network is investigated to assist blocked ground communication. According to the mobility and implementation flexibility of UAV relays, it is assumed that the communication link between air-to-ground is the Rician fading channel. On the basis of information causality, we derive the state change of the information in the buffer of UAV relays and maximize the end-to-end average throughput by join the relay selection, UAV transmit power, and UAV trajectory optimization. However, the considered problem is a mixed integer non-convex optimization problem, and therefore, it is difficult to solve directly with general optimization methods. In order to make the problem tractable, an efficient iterative algorithm based on the block coordinate descent and the successive convex optimization techniques is proposed. The convergence of the proposed algorithm will be verified analytically at the end of this paper. The simulation results show that by alternately optimizing the relay selection, UAV transmit power, and UAV trajectory, the proposed algorithm is able to achieve convergence quickly and significantly improve the average throughput, as compared to other benchmark schemes.


2018 ◽  
Vol 62 ◽  
pp. 579-664 ◽  
Author(s):  
Enrique Fernandez-Gonzalez ◽  
Brian Williams ◽  
Erez Karpas

The state of the art practice in robotics planning is to script behaviors manually, where each behavior is typically generated using trajectory optimization. However, in order for robots to be able to act robustly and adapt to novel situations, they need to plan these activity sequences autonomously. Since the conditions and effects of these behaviors are tightly coupled through time, state and control variables, many problems require that the tasks of activity planning and trajectory optimization are considered together. There are two key issues underlying effective hybrid activity and trajectory planning: the sufficiently accurate modeling of robot dynamics and the capability of planning over long horizons. Hybrid activity and trajectory planners that employ mixed integer programming within a discrete time formulation are able to accurately model complex dynamics for robot vehicles, but are often restricted to relatively short horizons. On the other hand, current hybrid activity planners that employ continuous time formulations can handle longer horizons but they only allow actions to have continuous effects with constant rate of change, and restrict the allowed state constraints to linear inequalities. This is insufficient for many robotic applications and it greatly limits the expressivity of the problems that these approaches can solve. In this work we present the ScottyActivity planner, that is able to generate practical hybrid activity and motion plans over long horizons by employing recent methods in convex optimization combined with methods for planning with relaxed plan graphs and heuristic forward search. Unlike other continuous time planners, ScottyActivity can solve a broad class of robotic planning problems by supporting convex quadratic constraints on state variables and control variables that are jointly constrained and that affect multiple state variables simultaneously. In order to support planning over long horizons, ScottyActivity does not resort to time, state or control variable discretization. While straightforward formulations of consistency checks are not convex and do not scale, we present an efficient convex formulation, in the form of a Second Order Cone Program (SOCP), that is very fast to solve. We also introduce several new realistic domains that demonstrate the capabilities and scalability of our approach, and their simplified linear versions, that we use to compare with other state of the art planners. This work demonstrates the power of integrating advanced convex optimization techniques with discrete search methods and paves the way for extensions dealing with non-convex disjoint constraints, such as obstacle avoidance.


Author(s):  
Sarat Chandra Nayak ◽  
Subhranginee Das ◽  
Mohammad Dilsad Ansari

Background and Objective: Stock closing price prediction is enormously complicated. Artificial Neural Networks (ANN) are excellent approximation algorithms applied to this area. Several nature-inspired evolutionary optimization techniques are proposed and used in the literature to search the optimum parameters of ANN based forecasting models. However, most of them need fine-tuning of several control parameters as well as algorithm specific parameters to achieve optimal performance. Improper tuning of such parameters either leads toward additional computational cost or local optima. Methods: Teaching Learning Based Optimization (TLBO) is a newly proposed algorithm which does not necessitate any parameters specific to it. The intrinsic capability of Functional Link Artificial Neural Network (FLANN) to recognize the multifaceted nonlinear relationship present in the historical stock data made it popular and got wide applications in the stock market prediction. This article presents a hybrid model termed as Teaching Learning Based Optimization of Functional Neural Networks (TLBO-FLN) by combining the advantages of both TLBO and FLANN. Results and Conclusion: The model is evaluated by predicting the short, medium, and long-term closing prices of four emerging stock markets. The performance of the TLBO-FLN model is measured through Mean Absolute Percentage of Error (MAPE), Average Relative Variance (ARV), and coefficient of determination (R2); compared with that of few other state-of-the-art models similarly trained and found superior.


Water ◽  
2019 ◽  
Vol 11 (7) ◽  
pp. 1392 ◽  
Author(s):  
Iram Parvez ◽  
JianJian Shen ◽  
Mehran Khan ◽  
Chuntian Cheng

The hydro generation scheduling problem has a unit commitment sub-problem which deals with start-up/shut-down costs related hydropower units. Hydro power is the only renewable energy source for many countries, so there is a need to find better methods which give optimal hydro scheduling. In this paper, the different optimization techniques like lagrange relaxation, augmented lagrange relaxation, mixed integer programming methods, heuristic methods like genetic algorithm, fuzzy logics, nonlinear approach, stochastic programming and dynamic programming techniques are discussed. The lagrange relaxation approach deals with constraints of pumped storage hydro plants and gives efficient results. Dynamic programming handles simple constraints and it is easily adaptable but its major drawback is curse of dimensionality. However, the mixed integer nonlinear programming, mixed integer linear programming, sequential lagrange and non-linear approach deals with network constraints and head sensitive cascaded hydropower plants. The stochastic programming, fuzzy logics and simulated annealing is helpful in satisfying the ramping rate, spinning reserve and power balance constraints. Genetic algorithm has the ability to obtain the results in a short interval. Fuzzy logic never needs a mathematical formulation but it is very complex. Future work is also suggested.


Author(s):  
Álinson S. Xavier ◽  
Ricardo Fukasawa ◽  
Laurent Poirrier

When generating multirow intersection cuts for mixed-integer linear optimization problems, an important practical question is deciding which intersection cuts to use. Even when restricted to cuts that are facet defining for the corner relaxation, the number of potential candidates is still very large, especially for instances of large size. In this paper, we introduce a subset of intersection cuts based on the infinity norm that is very small, works for relaxations having arbitrary number of rows and, unlike many subclasses studied in the literature, takes into account the entire data from the simplex tableau. We describe an algorithm for generating these inequalities and run extensive computational experiments in order to evaluate their practical effectiveness in real-world instances. We conclude that this subset of inequalities yields, in terms of gap closure, around 50% of the benefits of using all valid inequalities for the corner relaxation simultaneously, but at a small fraction of the computational cost, and with a very small number of cuts. Summary of Contribution: Cutting planes are one of the most important techniques used by modern mixed-integer linear programming solvers when solving a variety of challenging operations research problems. The paper advances the state of the art on general-purpose multirow intersection cuts by proposing a practical and computationally friendly method to generate them.


2018 ◽  
Author(s):  
Fabien Maussion ◽  
Anton Butenko ◽  
Julia Eis ◽  
Kévin Fourteau ◽  
Alexander H. Jarosch ◽  
...  

Abstract. Despite of their importance for sea-level rise, seasonal water availability, and as source of geohazards, mountain glaciers are one of the few remaining sub-systems of the global climate system for which no globally applicable, open source, community-driven model exists. Here we present the Open Global Glacier Model (OGGM, http://www.oggm.org), developed to provide a modular and open source numerical model framework for simulating past and future change of any glacier in the world. The modelling chain comprises data downloading tools (glacier outlines, topography, climate, validation data), a preprocessing module, a mass-balance model, a distributed ice thickness estimation model, and an ice flow model. The monthly mass-balance is obtained from gridded climate data and a temperature index melt model. To our knowledge, OGGM is the first global model explicitly simulating glacier dynamics: the model relies on the shallow ice approximation to compute the depth-integrated flux of ice along multiple connected flowlines. In this paper, we describe and illustrate each processing step by applying the model to a selection of glaciers before running global simulations under idealized climate forcings. Even without an in-depth calibration, the model shows a very realistic behaviour. We are able to reproduce earlier estimates of global glacier volume by varying the ice dynamical parameters within a range of plausible values. At the same time, the increased complexity of OGGM compared to other prevalent global glacier models comes at a reasonable computational cost: several dozens of glaciers can be simulated on a personal computer, while global simulations realized in a supercomputing environment take up to a few hours per century. Thanks to the modular framework, modules of various complexity can be added to the codebase, allowing to run new kinds of model intercomparisons in a controlled environment. Future developments will add new physical processes to the model as well as tools to calibrate the model in a more comprehensive way. OGGM spans a wide range of applications, from ice-climate interaction studies at millenial time scales to estimates of the contribution of glaciers to past and future sea-level change. It has the potential to become a self-sustained, community driven model for global and regional glacier evolution.


Sign in / Sign up

Export Citation Format

Share Document