Uncertainty Quantification and Data Assimilation

In many applications, initial data or model parameters of partial differential equations

are not known precisely. This issue can be addressed in several ways:

If the probability distribution of the unknown data is available one can investigate how this probability distribution is propagated via the solution operator of the PDE. This is possible by intrusive methods (e.g. stochastic Galerkin) as well as by non-intrusive methods (e.g. Monte-Carlo or stochastic collocation). We construct reliable and efficient numerical schemes for this task and derive error estimates for them.

Another scenario is that measurement data of the system state (e.g. values of the solution of the PDE at certain points) are available. Data assimilation aims at combining measurement data and models in such a way that reliable estimates of the system state can be obtained. We approach this task based on so-called observers, i.e. modified versions of the PDE that include measurement data. Our goal is to show that for long times the observer state is close to the state of the original system.

  • Gugat, Giesselmann, Kunkel: Exponential synchronization of a nodal observer for a semilinear model for the flow in gas networks, IMA J. Math. Control Inform. (2021)
  • Lang, Scheichl, Silvester: A Fully Adaptive Multilevel Stochastic Collocation Strategy for Solving Elliptic PDEs with Random Data, J. Comput. Phys. 419 (2020)
  • Meyer, Rohde, Giesselmann: A posteriori error analysis for random scalar conservation laws using the stochastic Galerkin method, IMA J. Numer. Anal. (2019)
  • Müller, Ullmann, Lang: A Bramble-Pasciak conjugate gradient method for discrete Stokes problems with random viscosity, SIAM/ASA J. Uncertainty Quantification 7 (2019)
  • Uncertainty quantification for phase-field models (Lang)