The Level Set Ensemble Kalman Filter: Sequential Data Assimilation for Flows With Shocks

Abstract

The ensemble Kalman filter is commonly used to perform sequential data assimilation in science and engineering because of its robust performance and scalability when the state space dimension is high. However, the standard EnKF generates spurious oscillations when applied to compressible flows with shocks; the cause is the uncertainty, across the ensemble, in the shock location. This paper develops the \emph{level set EnKF}, which overcomes this challenge by employing a nonlinear mapping from the original state space into a latent space in which the analysis step is undertaken. The latent representation comprises a level set function that encodes the discontinuity location and smooth state extensions, defined over the full physical domain, that represent the solution on either side of the discontinuity. We first formulate the level set EnKF as a specific instance of a more general approach of performing data assimilation in a latent space. We then demonstrate the feasibility of the level set EnKF for several one-dimensional compressible flows and a two-dimensional blast wave.