Closed-form discovery of structural errors in models of chaotic systems by integrating Bayesian sparse regression and data assimilation

10/01/2021
by   Rambod Mojgani, et al.
13

Models used for many important engineering and natural systems are imperfect. The discrepancy between the mathematical representations of a true physical system and its imperfect model is called the model error. These model errors can lead to substantial difference between the numerical solutions of the model and the observations of the system, particularly in those involving nonlinear, multi-scale phenomena. Thus, there is substantial interest in reducing model errors, particularly through understanding their physics and sources and leveraging the rapid growth of observational data. Here we introduce a framework named MEDIDA: Model Error Discovery with Interpretability and Data Assimilation. MEDIDA only requires a working numerical solver of the model and a small number of noise-free or noisy sporadic observations of the system. In MEDIDA, first the model error is estimated from differences between the observed states and model-predicted states (the latter are obtained from a number of one-time-step numerical integrations from the previous observed states). If observations are noisy, a data assimilation (DA) technique such as ensemble Kalman filter (EnKF) is first used to provide a noise-free analysis state of the system, which is then used in estimating the model error. Finally, an equation-discovery technique, such as the relevance vector machine (RVM), a sparsity-promoting Bayesian method, is used to identify an interpretable, parsimonious, closed-form representation of the model error. Using the chaotic Kuramoto-Sivashinsky (KS) system as the test case, we demonstrate the excellent performance of MEDIDA in discovering different types of structural/parametric model errors, representing different types of missing physics, using noise-free and noisy observations.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/04/2022

A multi-model ensemble Kalman filter for data assimilation and forecasting

Data assimilation (DA) aims to optimally combine model forecasts and noi...
research
06/25/2020

Combining Ensemble Kalman Filter and Reservoir Computing to predict spatio-temporal chaotic systems from imperfect observations and models

Prediction of spatio-temporal chaotic systems is important in various fi...
research
06/08/2023

Learning Closed-form Equations for Subgrid-scale Closures from High-fidelity Data: Promises and Challenges

There is growing interest in discovering interpretable, closed-form equa...
research
07/24/2018

Model uncertainty estimation in data assimilation for multi-scale systems with partially observed resolved variables

Model uncertainty quantification is an essential component of effective ...
research
09/07/2021

Combining data assimilation and machine learning to estimate parameters of a convective-scale model

Errors in the representation of clouds in convection-permitting numerica...
research
11/08/2021

A toolkit for data-driven discovery of governing equations in high-noise regimes

We consider the data-driven discovery of governing equations from time-s...

Please sign up or login with your details

Forgot password? Click here to reset