Hessian-Based Model Reduction for Large-Scale Data Assimilation Problems

Assimilation of spatially- and temporally-distributed state observations into simulations of dynamical systems stemming from discretized PDEs leads to inverse problems with high-dimensional control spaces in the form of discretized initial conditions. Solution of such inverse problems in “real-time” is often intractable. This motivates the construction of reduced-order models that can be used as surrogates of the high-fidelity simulations during inverse solution. For the surrogates to be useful, they must be able to approximate the observable quantities over a wide range of initial conditions. Construction of the reduced models entails sampling the initial condition space to generate an appropriate training set, which is an intractable proposition for high dimensional initial condition spaces unless the problem structure can be exploited. Here, we present a method that extracts the dominant spectrum of the inputoutput map (i.e. the Hessian of the least squares optimization problem...

• Type: Conference
• Year: 2007
• Where: ICCS
• Authors: Omar Bashir, Omar Ghattas, Judith Hill, Bart G. van Bloemen Waanders, Karen Willcox
• Keywords: ICCS 2007, Initial Condition, Initial Condition Spaces, Inverse Problems