(119a) A Branch-and-Bound Algorithm with Growing Datasets for Large-Scale Parameter Estimation | AIChE

(119a) A Branch-and-Bound Algorithm with Growing Datasets for Large-Scale Parameter Estimation

Authors 

Mitsos, A. - Presenter, RWTH Aachen University
Sass, S., RWTH Aachen University
Bongartz, D., KU Leuven
Bell, I. H., National Institute of Standards and Technology
Nikolov, N. I., Institute of Statistics, RWTH Aachen University
Accurate results of parameter estimation problems are the basis for model-based optimization of process design and control. Only deterministic global optimization can ensure that suboptimal solutions are excluded as a reason for model-data mismatch [1,2]. However, general nonconvex parameter estimation problems considering large measurement datasets are challenging for state-of-the-art solvers. Thus, we propose to extend a common optimization method, the spatial branch-and-bound (B&B) algorithm [3,4], by using growing datasets. In detail, we propose to start with a reduced dataset in the root node of the B&B tree and augment it successively until converging to the original full dataset. We present different rules determining for each B&B node whether to augment the dataset or branch the parameter domain. By extending a proof of Locatelli and Schoen [5], we can show for nonlinear programs (NLPs) that our extension retains the convergence of the B&B algorithm. We implement the proposed extension in our open-source solver MAiNGO [6] and test the algorithm for a real-world case study stemming from chemical engineering, namely fitting an equation of state of propane [7]. The numerical results for this case study show significant CPU time savings.

Acknowledgement

This work was partially funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under grant MI 1851/10-1 “Parameter estimation with (almost) deterministic global optimization”.

Susanne Sass is grateful for her association to the International Research Training Group (DFG) IRTG-2379 Modern Inverse Problems under grant 333849990/GRK2379.

Disclaimer: Commercial equipment, instruments, or materials are identified only in order to adequately specify certain procedures. In no case does such identification imply recommendation or endorsement by the National Institute of Standards and Technology, nor does it imply that the products identified are necessarily the best available for the purpose.

References

[1] Singer, A.B., Barton, P.I., 2006. Global Optimization with Nonlinear Ordinary Differential Equations. Journal of Global Optimization 34, 159–190. https://doi.org/10.1007/s10898-005-7074-4

[2] Mitsos, A., Chachuat, B., Barton, P.I., 2009. McCormick-Based Relaxations of Algorithms. SIAM Journal on Optimization 20, 573–601. https://doi.org/10.1137/080717341

[3] Horst, R., Tuy, H., 1996. Global Optimization: Deterministic Approaches. Springer Berlin. https://doi.org/10.1007/978-3-662-03199-5

[4] Tawarmalani, M., Sahinidis, N.V., 2002. Convexification and global optimization in continuous and mixed-integer nonlinear programming: Theory, Algorithms, Software, and Applications. volume 65. Springer US. https://doi.org/10.1007/978-1-4757-3532-1

[5] Locatelli, M., Schoen, F., 2013. Global optimization: Theory, Algorithms, and Applications. MOS-SIAM series on optimization. https://doi.org/10.1137/1.9781611972672

[6] Bongartz, D., Najman, J., Sass, S., Mitsos, A., 2018. MAiNGO - McCormick-based Algorithm for mixed-integer Nonlinear Global Optimization. http://permalink.avt.rwth-aachen.de/?id=729717 (Accessed March 20, 2023)

[7] Lemmon, E.W., McLinden, M.O.,Wagner, W., 2009. Thermodynamic Properties of Propane. III. A Reference Equation of State for Temperatures from the Melting Line to 650 K and Pressures up to 1000 MPa. J. Chem. Eng. Data 54, 3141–3180. https://doi.org/10.1021/je900217v