Building a Framework for Predictive Science

Key questions that scientists and engineers typically want to address can be formulated in terms of predictive science. Questions such as: "How well does my computational model represent reality?", "What are the most important parameters in the problem?", and "What is the best next experiment to perform?" are fundamental in solving scientific problems. Mystic is a framework for massively-parallel optimization and rigorous sensitivity analysis that enables these motivating questions to be addressed quantitatively as global optimization problems. Often realistic physics, engineering, and materials models may have hundreds of input parameters, hundreds of constraints, and may require execution times of seconds or longer. In more extreme cases, realistic models may be multi-scale, and require the use of high-performance computing clusters for their evaluation. Predictive calculations, formulated as a global optimization over a potential surface in design parameter space, may require an already prohibitively large simulation to be performed hundreds, if not thousands, of times. The need to prepare, schedule, and monitor thousands of model evaluations, and dynamically explore and analyze results, is a challenging problem that requires a software infrastructure capable of distributing and managing computations on large-scale heterogeneous resources. In this paper, we present the design behind an optimization framework, and also a framework for heterogeneous computing, that when utilized together, can make computationally intractable sensitivity and optimization problems much more tractable.

[1]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[2]  P. Venkataraman,et al.  Applied Optimization with MATLAB Programming , 2001 .

[3]  R. Storn,et al.  Differential Evolution - A simple and efficient adaptive scheme for global optimization over continuous spaces , 2004 .

[4]  Colin McDiarmid,et al.  Surveys in Combinatorics, 1989: On the method of bounded differences , 1989 .

[5]  Gábor Lugosi,et al.  Concentration Inequalities , 2008, COLT.

[6]  Michael S. Eldred,et al.  DAKOTA, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis Version 3.0 Reference Manual , 2001 .

[7]  I. N. Egorov,et al.  IOSO Optimization Toolkit - Novel Software to Create Better Design , 2002 .

[8]  Guido van Rossum,et al.  Python Programming Language , 2007, USENIX Annual Technical Conference.

[9]  S. N. Kramer,et al.  An Augmented Lagrange Multiplier Based Method for Mixed Integer Discrete Continuous Optimization and Its Applications to Mechanical Design , 1994 .

[10]  Houman Owhadi,et al.  Optimal Uncertainty Quantification , 2010, SIAM Rev..

[11]  J. F. A. de O. Pantoja Algorithms for constrained optimization , 1984 .

[12]  Ufuk Topcu,et al.  Uncertainty quantification via codimension‐one partitioning , 2011 .

[13]  G. Vetrovec DES , 2021, Encyclopedia of Systems and Control.

[14]  M. Ledoux The concentration of measure phenomenon , 2001 .

[15]  Houman Owhadi,et al.  Rigorous model-based uncertainty quantification with application to terminal ballistics—Part II. Systems with uncontrollable inputs and large scatter , 2012 .

[16]  Michael S. Eldred,et al.  DAKOTA, A Multilevel Parallel Object-Oriented Framework for Design Optimization, Parameter Estimation, Uncertainty Quantification, and Sensitivity Analysis Version 3.0 Developers Manual (title change from electronic posting) , 2002 .

[17]  Michael S. Eldred,et al.  DAKOTA : a multilevel parallel object-oriented framework for design optimization, parameter estimation, uncertainty quantification, and sensitivity analysis. Version 5.0, user's reference manual. , 2010 .

[18]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[19]  Houman Owhadi,et al.  Rigorous model-based uncertainty quantification with application to terminal ballistics—Part II. Systems with uncontrollable inputs and large scatter , 2012 .

[20]  Houman Owhadi,et al.  Rigorous verification, validation, uncertainty quantification and certification through concentration-of-measure inequalities , 2008 .

[21]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[22]  Edwin K. P. Chong,et al.  Algorithms for Constrained Optimization , 2011 .