Optimal Bounds on Nonlinear Partial Differential Equations in Model Certification, Validation, and Experimental Design
Authors:
M. McKerns,
F. J. Alexander,
K. S. Hickmann,
T. J. Sullivan,
D. E. Vaughan
Abstract:
We demonstrate that the recently developed Optimal Uncertainty Quantification (OUQ) theory, combined with recent software enabling fast global solutions of constrained non-convex optimization problems, provides a methodology for rigorous model certification, validation, and optimal design under uncertainty. In particular, we show the utility of the OUQ approach to understanding the behavior of a s…
▽ More
We demonstrate that the recently developed Optimal Uncertainty Quantification (OUQ) theory, combined with recent software enabling fast global solutions of constrained non-convex optimization problems, provides a methodology for rigorous model certification, validation, and optimal design under uncertainty. In particular, we show the utility of the OUQ approach to understanding the behavior of a system that is governed by a partial differential equation -- Burgers' equation. We solve the problem of predicting shock location when we only know bounds on viscosity and on the initial conditions. Through this example, we demonstrate the potential to apply OUQ to complex physical systems, such as systems governed by coupled partial differential equations. We compare our results to those obtained using a standard Monte Carlo approach, and show that OUQ provides more accurate bounds at a lower computational cost. We discuss briefly about how to extend this approach to more complex systems, and how to integrate our approach into a more ambitious program of optimal experimental design.
△ Less
Submitted 13 September, 2020;
originally announced September 2020.
A Multiresolution Ensemble Kalman Filter using Wavelet Decomposition
Authors:
Kyle S. Hickmann,
Humberto C. Godinez
Abstract:
We present a method of using classical wavelet based multiresolution analysis to separate scales in model and observations during data assimilation with the ensemble Kalman filter. In many applications, the underlying physics of a phenomena involve the interaction of features at multiple scales. Blending of observational and model error across scales can result in large forecast inaccuracies since…
▽ More
We present a method of using classical wavelet based multiresolution analysis to separate scales in model and observations during data assimilation with the ensemble Kalman filter. In many applications, the underlying physics of a phenomena involve the interaction of features at multiple scales. Blending of observational and model error across scales can result in large forecast inaccuracies since large errors at one scale are interpreted as inexact data at all scales. Our method uses a transformation of the observation operator in order to separate the information from different scales of the observations. This naturally induces a transformation of the observation covariance and we put forward several algorithms to efficiently compute the transformed covariance. Another advantage of our multiresolution ensemble Kalman filter is that scales can be weighted independently to adjust each scale's effect on the forecast. To demonstrate feasibility we present applications to a one dimensional Kuramoto-Sivashinsky (K-S) model with scale dependent observation noise and an application involving the forecasting of solar photospheric flux. The latter example demonstrates the multiresolution ensemble Kalman filter's ability to account for scale dependent model error. Modeling of photospheric magnetic flux transport is accomplished by the Air Force Data Assimilative Photospheric Transport (ADAPT) model.
△ Less
Submitted 5 November, 2015;
originally announced November 2015.