Determining white noise forcing from Eulerian observations in the Navier-Stokes equation

The Bayesian approach to inverse problems is of paramount importance in quantifying uncertainty about the input to, and the state of, a system of interest given noisy observations. Herein we consider the forward problem of the forced 2D Navier-Stokes equation. The inverse problem is to make inference concerning the forcing, and possibly the initial condition, given noisy observations of the velocity field. We place a prior on the forcing which is in the form of a spatially-correlated and temporally-white Gaussian process, and formulate the inverse problem for the posterior distribution. Given appropriate spatial regularity conditions, we show that the solution is a continuous function of the forcing. Hence, for appropriately chosen spatial regularity in the prior, the posterior distribution on the forcing is absolutely continuous with respect to the prior and is hence well-defined. Furthermore, it may then be shown that the posterior distribution is a continuous function of the data. We complement these theoretical results with numerical simulations showing the feasibility of computing the posterior distribution, and illustrating its properties.

[1]  Andrew C. Lorenc,et al.  The potential of the ensemble Kalman filter for NWP—a comparison with 4D‐Var , 2003 .

[2]  Sari Lasanen,et al.  Posterior convergence for approximated unknowns in non-Gaussian statistical inverse problems , 2011, 1112.0906.

[3]  Sebastian J. Vollmer,et al.  Dimension-Independent MCMC Sampling for Inverse Problems with Non-Gaussian Priors , 2013, SIAM/ASA J. Uncertain. Quantification.

[4]  Kody J. H. Law Proposals which speed up function-space MCMC , 2014, J. Comput. Appl. Math..

[5]  G. Roberts,et al.  MCMC Methods for Functions: ModifyingOld Algorithms to Make Them Faster , 2012, 1202.0709.

[6]  Hantaek Bae Navier-Stokes equations , 1992 .

[7]  Andrew M. Stuart,et al.  Inverse problems: A Bayesian perspective , 2010, Acta Numerica.

[8]  A. Stuart,et al.  MAP estimators and their consistency in Bayesian nonparametric inverse problems , 2013, 1303.4795.

[9]  D. Menemenlis Inverse Modeling of the Ocean and Atmosphere , 2002 .

[10]  Faming Liang,et al.  Statistical and Computational Inverse Problems , 2006, Technometrics.

[11]  D. Zupanski A General Weak Constraint Applicable to Operational 4DVAR Data Assimilation Systems , 1997 .

[12]  Franco Flandoli,et al.  Dissipativity and invariant measures for stochastic Navier-Stokes equations , 1994 .

[13]  A. Stuart,et al.  Signal processing problems on function space: Bayesian formulation, stochastic PDEs and effective MCMC methods , 2011 .

[14]  L. Tierney A note on Metropolis-Hastings kernels for general state spaces , 1998 .

[15]  Jonathan C. Mattingly Ergodicity of 2D Navier–Stokes Equations with¶Random Forcing and Large Viscosity , 1999 .

[16]  James C. Robinson,et al.  Bayesian inverse problems for functions and applications to fluid mechanics , 2009 .

[17]  P. Kloeden,et al.  TAYLOR EXPANSIONS OF SOLUTIONS OF STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS WITH ADDITIVE NOISE , 2009, 1010.0161.

[18]  S. Lasanen,et al.  Measurements and infinite‐dimensional statistical inverse theory , 2007 .

[19]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[20]  Sari Lasanen,et al.  Non-Gaussian statistical inverse problems. Part I: Posterior distributions , 2012 .

[21]  T. Kurtz,et al.  Stochastic equations in infinite dimensions , 2006 .