@article{LeKoch2023, author = {Le, Thi Thai and Koch, Thorsten}, title = {Effect of inertia force on the interface stability of a tangential-velocity discontinuity in porous media}, volume = {169}, journal = {International Journal of Multiphase Flow}, doi = {10.1016/j.ijmultiphaseflow.2023.104612}, year = {2023}, abstract = {The present study investigates the stability of a tangential-velocity discontinuity in porous media during the withdrawing and injecting processes of natural gases from and into an underground gas storage. The focus is placed on analyzing the impact of inertia forces on the interface stability using the Forchheimer equations. Other publications have relied primarily on Darcy's law to describe flow stability in porous media. However, Darcy's law only adequately describes flows in which viscous forces dominate over inertia forces. As the flow rate increases, the significance of inertia forces becomes more pronounced, and Darcy's law becomes insufficient for considering such flows. Our findings indicate that even a slight consideration of the inertia effect leads to permanent destabilization of the discontinuity interface, regardless of the fluid viscosity or the Mach number. In contrast, when the inertia effect is neglected, the interface is stabilized across the entire Mach number range if the fluid viscosity is strong enough.}, language = {en} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method}, journal = {The 35th International Conference on Efficiency, Cost, Optimization, Simulation and Environmental Impact of Energy Systems (ECOS 2022)}, pages = {1891 -- 1902}, year = {2022}, abstract = {Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. In addition, some strategies have been proposed to enhance the computation efficiency furthermore. As one of the strategies, a method of reducing model by time aggregation has been proposed to search design candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method with the strategies is extendedly applied to the optimal design of energy supply systems with storage units. Especially, the method of reducing model is extended by aggregating both representative days and sampling times separately in consideration of the characteristics of energy storage units. A case study is conducted on the optimal design of a cogeneration system with a thermal storage tank. Through the study, it turns out the hierarchical MILP method is effective to derive the optimal solutions in short computation times. It also turns out that the model reduction with day and time aggregations is effec- tive to shorten the computation times furthermore when the number of candidates for equipment capacities is relatively small.}, language = {en} } @article{YokoyamaShinanoWakui2022, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {Day and time aggregations for optimal design of energy supply systems with storage units by a hierarchical MILP method}, journal = {PROCEEDINGS OF ECOS 2022 - THE 35 TH INTERNATIONAL CONFERENCE ON EFFICIENCY, COST, OPTIMIZATION, SIMULATION AND ENVIRONMENTAL IMPACT OF ENERGY SYSTEMS}, pages = {1891 -- 1902}, year = {2022}, abstract = {Mixed-integer linear programming (MILP) methods have been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems efficiently. In addition, some strategies have been proposed to enhance the computation efficiency furthermore. As one of the strategies, a method of reducing model by time aggregation has been proposed to search design candidates efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method with the strategies is extendedly applied to the optimal design of energy supply systems with storage units. Especially, the method of reducing model is extended by aggregating both representative days and sampling times separately in consideration of the characteristics of energy storage units. A case study is conducted on the optimal design of a cogeneration system with a thermal storage tank. Through the study, it turns out the hierarchical MILP method is effective to derive the optimal solutions in short computation times. It also turns out that the model reduction with day and time aggregations is effective to shorten the computation times furthermore when the number of candidates for equipment capacities is relatively small.}, language = {en} } @article{Rehfeldt2022, author = {Rehfeldt, Daniel}, title = {Faster Algorithms for Steiner Tree and related Problems: From Theory to Practice}, journal = {OR News: Das Magazin der GOR}, year = {2022}, language = {en} } @inproceedings{ZakiyevaPetkovic2022, author = {Zakiyeva, Nazgul and Petkovic, Milena}, title = {Modeling and Forecasting Gas Network Flows with Multivariate Time Series and Mathematical Programming Approach}, booktitle = {Operations Research Proceedings 2021. OR 2021.}, publisher = {Springer, Cham}, doi = {10.1007/978-3-031-08623-6_31}, pages = {200 -- 205}, year = {2022}, abstract = {With annual consumption of approx. 95 billion cubic meters and similar amounts of gas just transshipped through Germany to other EU states, Germany's gas transport system plays a vital role in European energy supply. The complex, more than 40,000 km long high-pressure transmission network is controlled by several transmission system operators (TSOs) whose main task is to provide security of supply in a cost-efficient way. Given the slow speed of gas flows through the gas transmission network pipelines, it has been an essential task for the gas network operators to enhance the forecast tools to build an accurate and effective gas flow prediction model for the whole network. By incorporating the recent progress in mathematical programming and time series modeling, we aim to model natural gas network and predict gas in- and out-flows at multiple supply and demand nodes for different forecasting horizons. Our model is able to describe the dynamics in the network by detecting the key nodes, which may help to build an optimal management strategy for transmission system operators.}, language = {en} } @misc{ShinanoVigerske2025, author = {Shinano, Yuji and Vigerske, Stefan}, title = {Smoothie: Mixing the strongest MIP solvers to solve hard MIP instances on supercomputers - Phase I development}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-101906}, year = {2025}, abstract = {Mixed-Integer Linear Programming (MIP) is applicable to such a wide range of real-world decision problems that the competition for the best code to solve such problems has lead to tremendous progress over the last decades. While current solvers can solve some of the problems that seemed completely out-of-reach just 10 years ago, there are always relevant MIP problems that currently cannot be solved. With the Smoothie solver we intend to solve extremely hard MIP problems by building on the many years that went into the development of several state-of-the-art MIP solvers and by utilizing some of the largest computing resources available. The high-level task parallelization framework UG (Ubiquity Generator) is used and extended by Smoothie to build a solver that uses large-scale parallelization to distribute the solution of a single MIP on a shared- or distributed-memory computing infrastructure, thereby employing several established MIP solvers simultaneously. For the first development phase, which is the topic of this report, both FICO Xpress and Gurobi are used in concurrent mode on a single machine, while information on incumbent solutions and explored branch-and-bound subtrees is exchanged. A dynamic restarting mechanism ensures that solver configurations are selected that promise most suitable for the MIP to be solved. We report on initial findings using this early version of Smoothie on unsolved problems from MIPLIB 2017.}, language = {en} } @article{LePetkovic2025, author = {Le, Thi Thai and Petkovic, Milena}, title = {Predicting Fluid Interface Instability in Energy Systems for Sustainable Energy Transition}, journal = {Operations Research Proceedings 2023. OR 2023.}, publisher = {Springer}, doi = {10.1007/978-3-031-58405-3_45}, pages = {349 -- 356}, year = {2025}, abstract = {Due to the coexistence of different gases in underground storage, this work explores the interface stability's impact on energy storage, specifically during the injection and withdrawal of gases such as hydrogen and natural gas. A new approach of combing simulation and time series analysis is used to accurately predict instability modes in energy systems. Our simulation is based on the 2D Euler equations, solved using a second-order finite volume method with a staggered grid. The solution is validated by comparing them to experimental data and analytical solutions, accurately predicting the instability's behavior. We use time series analysis and state-of-the-art regime-switching methods to identify critical features of the interface dynamics, providing crucial insights into system optimization and design.}, language = {en} } @misc{ZakiyevaPetkovic2023, author = {Zakiyeva, Nazgul and Petkovic, Milena}, title = {High-dimensional high-frequency time series prediction with a mixed integer optimisation method}, journal = {Operations Research Proceedings}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-93114}, year = {2023}, language = {en} } @article{LeTranPhung2016, author = {Le, Thi Thai and Tran, Van Cuc and Phung, Dang Hieu}, title = {Numerical simulation of tidal flow in Danang Bay Based on non-hydrostatic shallow water equations}, volume = {8}, journal = {Pacific Journal of Mathematics for Industry}, number = {1}, doi = {10.1186/s40736-015-0020-6}, year = {2016}, abstract = {This paper presents a numerical simulation of the tidal flow in Danang Bay (Vietnam) based on the non-hydrostatic shallow water equations. First, to test the simulation capability of the non-hydrostatic model, we have made a test simulation comparing it with the experiment by Beji and Battjes 1993. Simulation results for this case are compared with both the experimental data and calculations obtained from the traditional hydrostatic model. It is shown that the non-hydrostatic model is better than the hydrostatic model when the seabed topography variation is complex. The usefulness of the non-hydrostatic model is father shown by successfully simulating the tidal flow of Danang Bay.}, language = {en} } @article{LeFukumotoKoch2024, author = {Le, Thi Thai and Fukumoto, Yasuhide and Koch, Thorsten}, title = {Linear stability of a simple shear layer between two parallel streams in a shallow water flow}, volume = {493}, journal = {Physics Letters A}, publisher = {Elservier}, doi = {10.1016/j.physleta.2023.129264}, year = {2024}, abstract = {The stability of shear layers in fluid flows is a crucial factor in forming vortices and jets and plays a fundamental role in the development of turbulence. Such shear layer instabilities are ubiquitous in natural phenomena, such as atmospheric and oceanic flows, contributing to the formation of weather systems and predicting tsunamis. This study specifically focuses on the stability of a shear layer sandwiched between two semi-infinite layers within a two-dimensional flow. The velocity profile of the shear layer is assumed to be linearly dependent on the vertical coordinate, while the velocity of the other layers remains uniform with differing strengths. The effect of viscosity and surface tension is ignored to simplify the analysis. The shallow water equations are used to analyze the interface stability of the shear layer, and the resulting dispersion relation between wave frequency and other wave characteristics is obtained. This relation incorporates Whittaker functions and their first derivatives and is used to derive appropriate limits corresponding to various physical conditions. Our study thus contributes to a deeper understanding of the stability of shear layers and their implications for natural phenomena.}, language = {en} } @article{YokoyamaShinanoWakui2023, author = {Yokoyama, Ryohei and Shinano, Yuji and Wakui, Tetsuya}, title = {An effective approach for deriving and evaluating approximate optimal design solutions of energy supply systems by time series aggregation}, volume = {11}, journal = {Frontiers in Energy Research}, doi = {10.3389/fenrg.2023.1128681}, pages = {1 -- 14}, year = {2023}, abstract = {It is important to design multi-energy supply systems optimally in consideration of their operations for variations in energy demands. An approach for efficiently solving such an optimal design problem with a large number of periods for variations in energy demands is to derive an approximate optimal design solution by time series aggregation. However, such an approach does not provide any information on the accuracy for the optimal value of the objective function. In this paper, an effective approach for time series aggregation is proposed to derive an approximate optimal design solution and evaluate a proper gap between the upper and lower bounds for the optimal value of the objective function based on a mixed-integer linear model. In accordance with aggregation, energy demands are relaxed to uncertain parameters and the problem for deriving an approximate optimal design solution and evaluating it is transformed to a three-level optimization problem, and it is solved by applying both the robust and hierarchical optimization methods. A case study is conducted on a cogeneration system with a practical configuration, and it turns out that the proposed approach enables one to derive much smaller gaps as compared with those obtained by a conventional approach.}, language = {en} } @article{品野2020, author = {品野, 勇治}, title = {MIPLIB と Hans Mittelmann's benchmarks}, volume = {65}, journal = {機関誌「オペレーションズ・リサーチ」}, number = {1}, pages = {49 -- 56}, year = {2020}, abstract = {混合整数計画法 (Mixed Integer Programming: MIP) は,MIP を解くソフトウェアである MIP ソルバが大規模な現実問題を解けるようになったこともあり,現実問題を解く有用な OR の手法として広く知られるようになった.しかしながら,MIP ソルバの開発に欠かせないベンチマーク・データセットおよび性能測定方法についてはそれほど広く知られているとは言い難い.ベンチマーク・データセットは注意を払って作成しないと,多くのバイアスがかかってしまう.それらのバイアスを可能な限りのぞき,真に有用なベンチマーク・テストの結果を得るためには複数の人数で多大な労力を割く必要がある.本稿では,そのような MIP ソルバ開発の背景として重要な役割を果たしてきた MIPLIB と Hans Mittelmann's benchmarks について解説する.また,本稿において Hans Mittelmann's benchmarks は,BENCHMARKS FOR OPTIMIZATION SOFTWAREのページ (http://plato.asu.edu/bench.html) に示されているベンチマークである.}, language = {ja} } @article{品野藤井2019, author = {品野, 勇治 and 藤井, 浩一}, title = {使ってみよう線形計画ソルバ}, volume = {64}, journal = {機関誌「オペレーションズ・リサーチ」}, number = {4}, pages = {238 -- 245}, year = {2019}, abstract = {本稿では,線形計画問題を解くソフトウェアである線形計画ソルバの標準的な利用方法と,原稿執筆時点で利用可能なソルバを可能な限り紹介する.筆者らはソルバ開発コミュニティの中で仕事をしているので,この機会に現在のソルバ開発現場の様子も紹介するとともに,各ソルバの特徴として何を紹介するべきかは,できる限り開発者,または,開発者に近い研究者に問い合わせて記述した.本稿は,多くの線形計画ソルバの選択肢がある中で,読者がもつ問題を解くのに最も適当なソルバを選択するための指針を与えることを主たる目的としている.}, language = {ja} } @inproceedings{PetkovicZittel2024, author = {Petkovic, Milena and Zittel, Janina}, title = {Forecasting and modeling the dynamics of large-scale energy networks under the supply and demand balance constraint}, volume = {12}, booktitle = {Optimization in Green Sustainability and Ecological Transition. ODS 2023}, doi = {10.1007/978-3-031-47686-0_29}, year = {2024}, abstract = {With the emergence of "Big Data" the analysis of large data sets of high-dimensional energy time series in network structures have become feasible. However, building large-scale data-driven and computationally efficient models to accurately capture the underlying spatial and temporal dynamics and forecast the multivariate time series data remains a great challenge. Additional constraints make the problem more challenging to solve with conventional methods. For example, to ensure the security of supply, energy networks require the demand and supply to be balanced. This paper introduces a novel large-scale Hierarchical Network Regression model with Relaxed Balance constraint (HNR-RB) to investigate the network dynamics and predict multistep-ahead flows in the natural gas transmission network, where the total in- and out-flows of the network have to be balanced over a period of time. We concurrently address three main challenges: high dimensionality of networks with more than 100 nodes, unknown network dynamics, and constraint of balanced supply and demand in the network. The effectiveness of the proposed model is demonstrated through a real-world case study of forecasting demand and supply in a large-scale natural gas transmission network. The results demonstrate that HNR-RB outperforms alternative models for short- and mid-term horizons.}, language = {en} } @article{TateiwaShinanoYasudaetal.2023, author = {Tateiwa, Nariaki and Shinano, Yuji and Yasuda, Masaya and Kaji, Shizuo and Yamamura, Keiichiro and Fujisawa, Katsuki}, title = {Development and analysis of massive parallelization of a lattice basis reduction algorithm}, volume = {40}, journal = {Japan Journal of Industrial and Applied Mathematics}, number = {1}, doi = {10.1007/s13160-023-00580-z}, pages = {13 -- 56}, year = {2023}, abstract = {The security of lattice-based cryptography relies on the hardness of solving lattice problems. Lattice basis reduction is a strong tool for solving lattice problems, and the block Korkine-Zolotarev (BKZ) reduction algorithm is the de facto standard in cryptanalysis. We propose a parallel algorithm of BKZ-type reduction based on randomization. Randomized copies of an input lattice basis are independently reduced in parallel, while several basis vectors are shared asynchronously among all processes. There is a trade-off between randomization and information sharing; if a substantial amount of information is shared, all processes might work on the same problem, which diminishes the benefit of parallelization. To monitor the balance between randomness and sharing, we propose a new metric to quantify the variety of lattice bases, and we empirically find an optimal parameter of sharing for high-dimensional lattices. We also demonstrate the effectiveness of our parallel algorithm and metric through experiments from multiple perspectives.}, language = {en} } @article{LeKoch2021, author = {Le, Thi Thai and Koch, Thorsten}, title = {Instability of a tangential discontinuity surface in a three-dimensional compressible medium}, volume = {33}, journal = {Physics of Fluids}, number = {1}, doi = {10.1063/5.0033753}, pages = {016106}, year = {2021}, abstract = {Compressible flows appear in many natural and technological processes, for instance, the flow of natural gases in a pipe system. Thus, a detailed study of the stability of tangential velocity discontinuity in compressible media is relevant and necessary. The first early investigation in two-dimensional (2D) media was given more than 70 years ago. In this article, we continue investigating the stability in three-dimensional (3D) media. The idealized statement of this problem in an infinite spatial space was studied by Syrovatskii in 1954. However, the omission of the absolute sign of cos θ with θ being the angle between vectors of velocity and wave number in a certain inequality produced the inaccurate conclusion that the flow is always unstable for entire values of the Mach number M. First, we revisit this case to arrive at the correct conclusion, namely that the discontinuity surface is stabilized for a large Mach number with a given value of the angle θ. Next, we introduce a real finite spatial system such that it is bounded by solid walls along the flow direction. We show that the discontinuity surface is stable if and only if the dispersion relation equation has only real roots, with a large value of the Mach number; otherwise, the surface is always unstable. In particular, we show that a smaller critical value of the Mach number is required to make the flow in a narrow channel stable.}, language = {en} } @inproceedings{YokoyamaKamadaShinanoetal.2020, author = {Yokoyama, Ryohei and Kamada, Hiroki and Shinano, Yuji and Wakui, Tetsuya}, title = {A Hierarchical Optimization Approach to Robust Design of Energy Supply Systems Based on a Mixed-Integer Linear Model}, booktitle = {Conference}, pages = {601 -- 613}, year = {2020}, abstract = {In designing energy supply systems, designers should heighten the robustness in performance criteria against the uncertainty in energy demands. In this paper, a robust optimal design method using a hierarchi- cal mixed-integer linear programming (MILP) method is proposed to maximize the robustness of energy sup- ply systems under uncertain energy demands based on a mixed-integer linear model. A robust optimal design problem is formulated as a three-level min-max-min MILP one by expressing uncertain energy demands by intervals, evaluating the robustness in a performance criterion based on the minimax regret cri- terion, and considering relationships among integer design variables, uncertain energy demands, and inte- ger and continuous operation variables. This problem is solved by evaluating upper and lower bounds for the minimum of the maximum regret of the performance criterion repeatedly outside, and evaluating lower and upper bounds for the maximum regret repeatedly inside. Since these different types of optimization problems are difficult to solve even using commercial MILP solvers, they are solved by applying a hierarchi- cal MILP method developed for ordinary optimal design problems with its modifications. In a case study, the proposed approach is applied to the robust optimal design of a cogeneration system. Through the study, its validity and effectiveness are ascertained, and some features of the obtained robust designs are clarified.}, language = {en} } @inproceedings{YokoyamaTakeuchiShinanoetal.2020, author = {Yokoyama, Ryohei and Takeuchi, Kotaro and Shinano, Yuji and Wakui, Tetsuya}, title = {Effect of Model Reduction by Time Aggregation in Multiobjective Optimal Design of Energy Supply Systems by a Hierarchical MILP Method}, booktitle = {Conference}, pages = {627 -- 639}, year = {2020}, abstract = {The mixed-integer linear programming (MILP) method has been applied widely to optimal design of energy supply systems. A hierarchical MILP method has been proposed to solve such optimal design problems effi- ciently. As one of the strategies to enhance the computation efficiency furthermore, a method of reducing model by time aggregation has been proposed to search design candidates accurately and efficiently in the relaxed optimal design problem at the upper level. In this paper, the hierarchical MILP method and model reduction by time aggregation are applied to the multiobjective optimal design. In applying the model reduc- tion, the methods of clustering periods by the order of time series, based on an operational strategy, and by the k-medoids method are applied. As a case study, the multiobjective optimal design of a gas turbine cogeneration system with a practical configuration is investigated by adopting the annual total cost and pri- mary energy consumption as the objective functions to be minimized simultaneously, and the clustering methods are compared with one another in terms of the computation efficiency. It turns out that the model reduction by any clustering method is effective to enhance the computation efficiency when importance is given to minimizing the first objective function. It also turns out that the model reduction only by the k- medoids method is effective very limitedly when importance is given to minimizing the second objective function.}, language = {en} } @misc{PetkovicZakiyeva2022, author = {Petkovic, Milena and Zakiyeva, Nazgul}, title = {Mathematical Optimization for Analyzing and Forecasting Nonlinear Network Time Series}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-88037}, year = {2022}, abstract = {This work presents an innovative short to mid-term forecasting model that analyzes nonlinear complex spatial and temporal dynamics in energy networks under demand and supply balance constraints using Network Nonlinear Time Series (TS) and Mathematical Programming (MP) approach. We address three challenges simultaneously, namely, the adjacency matrix is unknown; the total amount in the network has to be balanced; dependence is unnecessarily linear. We use a nonparametric approach to handle the nonlinearity and estimate the adjacency matrix under the sparsity assumption. The estimation is conducted with the Mathematical Optimisation method. We illustrate the accuracy and effectiveness of the model on the example of the natural gas transmission network of one of the largest transmission system operators (TSOs) in Germany, Open Grid Europe. The obtained results show that, especially for shorter forecasting horizons, the proposed method outperforms all considered benchmark models, improving the average nMAPE for 5.1\% and average RMSE for 79.6\% compared to the second-best model. The model is capable of capturing the nonlinear dependencies in the complex spatial-temporal network dynamics and benefits from both sparsity assumption and the demand and supply balance constraint.}, language = {en} } @misc{HosodaMaherShinano2024, author = {Hosoda, Junko and Maher, Stephen J. and Shinano, Yuji}, title = {Combining Benders' decomposition and parallelisation to solve large-scale adaptive network supply chain management problems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-98744}, year = {2024}, abstract = {Integrated supply chain network design—combining warehouse clustering with inter- and intra-cluster vehicle routing—is a fundamental feature of supply chain management problems. In this paper, a Benders' decomposition-based parallel heuristic approach for the integrated long-haul and local vehicle routing problem is proposed. The proposed algorithm comprises lower- and upper-bound search algorithms running in parallel to find high-quality solutions for challenging supply chain management problems. A novel Benders' decomposition approach is proposed to aid the improvement in the lower bound for large-scale problems. The results demonstrate that the Benders' decomposition approach is effective in finding strong lower bounds and is beneficial for obtaining high-quality solutions.}, language = {en} }