# The Library

### Approximation of Bayesian inverse problems for PDEs

Tools

Cotter, Simon L., Dashti, M. and Stuart, A. M..
(2010)
*Approximation of Bayesian inverse problems for PDEs.*
SIAM Journal on Numerical Analysis, Vol.48
(No.1).
pp. 322-345.
ISSN 0036-1429

PDF
WRAP_Cotter_approximation_Bayesian_inverse.pdf - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader Download (577Kb) |

Official URL: http://dx.doi.org/10.1137/090770734

## Abstract

Inverse problems are often ill posed, with solutions that depend sensitively on data. In any numerical approach to the solution of such problems, regularization of some form is needed to counteract the resulting instability. This paper is based on an approach to regularization, employing a Bayesian formulation of the problem, which leads to a notion of well posedness for inverse problems, at the level of probability measures. The stability which results from this well posedness may be used as the basis for quantifying the approximation, in finite dimensional spaces, of inverse problems for functions. This paper contains a theory which utilizes this stability property to estimate the distance between the true and approximate posterior distributions, in the Hellinger metric, in terms of error estimates for approximation of the underlying forward problem. This is potentially useful as it allows for the transfer of estimates from the numerical analysis of forward problems into estimates for the solution of the related inverse problem. It is noteworthy that, when the prior is a Gaussian random field model, controlling differences in the Hellinger metric leads to control on the differences between expected values of polynomially bounded functions and operators, including the mean and covariance operator. The ideas are applied to some non-Gaussian inverse problems where the goal is determination of the initial condition for the Stokes or Navier–Stokes equation from Lagrangian and Eulerian observations, respectively.

Item Type: | Journal Article |
---|---|

Subjects: | Q Science > QA Mathematics |

Divisions: | Faculty of Science > Mathematics |

Library of Congress Subject Headings (LCSH): | Differential equations, Partial -- Improperly posed problems, Inverse problems (Differential equations), Bayesian statistical decision theory, Stokes flow |

Journal or Publication Title: | SIAM Journal on Numerical Analysis |

Publisher: | Society for Industrial and Applied Mathematics |

ISSN: | 0036-1429 |

Date: | 2010 |

Volume: | Vol.48 |

Number: | No.1 |

Number of Pages: | 24 |

Page Range: | pp. 322-345 |

Identification Number: | 10.1137/090770734 |

Status: | Peer Reviewed |

Publication Status: | Published |

Access rights to Published version: | Open Access |

Funder: | Engineering and Physical Sciences Research Council (EPSRC), European Research Council (ERC) |

References: | [1] A. Apte, C. K. R. T. Jones, A. M. Stuart, and J. Voss, Data assimilation: Mathematical and statistical perspectives, Internat. J. Numer. Methods Fluids, 56 (2008), pp. 1033–1046. [2] S. R. Arridge, J. P. Kaipio, V. Kolehmainen, M. Schweiger, E. Somersalo, and T. Tarvainen, Approximation errors and model reduction with an application to optical diffusion tomography, Inverse Problems, 22 (2006), pp. 175–195. [3] V. I. Bogachev, Gaussian Measures, American Mathematical Society, Providence, RI, 1998. [4] D. Calvetti and E. Somersalo, Introduction to bayesian scientific computing, Surveys and Tutorials in the Applied Mathematical Sciences 2, Springer, New York, 2007. [5] S. L. Cotter, M. Dashti, J. C. Robinson, and A. M. Stuart, Bayesian inverse problems for functions and applications to fluid mechanics, Inverse Problems, 25 (2010), p. 115008. [6] S. L. Cotter, M. Dashti, J. C. Robinson, and A. M. Stuart, MCMC Methods on Function Space and Applications to Fluid Mechanics, in preparation, 2010. [7] M. Dashti and J. C. Robinson, A simple proof of uniqueness of the particle trajectories for solutions of the Navier-Stokes equations, Nonlinearity, 22 (2009), pp. 735–746. [8] H. K. Engl, M. Hanke, and A. Neubauer, Regularization of inverse problems, Kluwer, Dordrecht, the Netherlands, 1996. [9] J. N. Franklin, Well-posed stochastic extensions of ill posed linear problems, J. Math. Anal. Appl., 31 (1970), pp. 682–716. [10] A. Hofinger and H. K. Pikkarainen, Convergence rates for the Bayesian approach to linear inverse problems, Inverse Problems, 23 (2007), pp. 2469–2484. [11] K. Ide, M. Courier, M. Ghil, and A. Lorenc, Unified notation for assimilation: Operational, sequential and variational, J. Met. Soc. Japan, 75 (1997), pp. 181–189. [12] J. Kaipio and E. Somersalo, Statistical and computational inverse problems, Applied Mathematical Sciences 160, Springer, New York, 2005. [13] J. P. Kaipio and E. Somersalo, Approximation errors in nonstationary inverse problems, Inverse Prob. Imag., 1 (2007), pp. 77–93. [14] J. P. Kaipio and E. Somersalo, Statistical inverse problems: Discretization, model reduction and inverse crimes, J. Comput. Appl. Math., 198 (2007), pp. 493–504. [15] A. Kirsch, An Introduction to the Mathematical Theory of Inverse Problems, Springer, New York, 1996. [16] M. S. Lehtinen, L. Paivarinta, and E. Somersalo, Linear inverse problems for generalized random variables, Inverse Problems, 5 (1989), pp. 599–612. [17] J. Liu, Monte Carlo strategies in scientific computing, Springer Texts in Statistics, Springer, New York, 2001. [18] A. Neubauer and H. K. Pikkarainen, Convergence results for the Bayesian inversion theory, J. Inverse Ill-Posed Probl., 16 (2008), pp. 601–613. [19] N. K. Nichols, Data assimilation: Aims and basic concepts, Data Assimilation for the Earth System, R. Swinbank, V. Shutyaev, and W. A. Lahoz, eds., Kluwer, Dordrecht, the Netherlands, 2003, pp. 9–20. [20] A. Pazy, Semigroups of Linear Operators and Applications to Partial Differential Equations, Springer, New York, 1983. [21] H. K. Pikkarainen, State estimation approach to nonstationary inverse problems: Discretization error and filtering problem, Inverse Problems, 22 (2006), pp. 365–379. [22] G. Da Prato and J. Zabczyk, Stochastic equations in infinite dimensions, Encyclopedia of Mathematics and its Applications 44, Cambridge University Press, Cambridge, 1992. [23] C. Schwab and R. A. Todor, Karhunen-Loeve approximation of random fields in domains by generalized fast multipole methods, J. Comput. Phys., 217 (2006), pp. 100–122. [24] I. Sloan and H. Wozniakowski, When are quasi-Monte Carlo algorithms efficient for highdimensional integrals? J. Complexity, 14 (1998), pp. 1–33. [25] P. D. Spanos and R. Ghanem, Stochastic finite element expansion for random media, J. Eng. Mech., 115 (1989), pp. 1035–1053. [26] P. D. Spanos and R. Ghanem, Stochastic Finite Elements: A Spectral Approach, Dover, New York, 2003. [27] A. M. Stuart, Inverse problems: A Bayesian approach, Acta Numer., 19 (2010), to appear. [28] A. Tarantola, Inverse Problem Theory. SIAM, Philadelphia, PA, 2005. [29] R. Temam, Navier-Stokes equations and nonlinear functional analysis, CBMS-NSF Regional Conference Series in Applied Mathematics 66, SIAM, Philadelphia, PA, 1995. [30] R. Temam, Navier-Stokes equations. AMS Chelsea, Providence, RI, 2001. [31] R. A. Todor and C. Schwab, Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients, IMA J. Numer. Anal., 27 (2007), pp. 232–261. [32] C. R. Vogel, Computational Methods for Inverse Problems, SIAM, Philadelphia, PA, 2002. [33] L. W. White, A study of uniqueness for the initialiazation problem for burgers’ equation, J. Math. Anal. Appl., 172 (1993), pp. 412–431. |

URI: | http://wrap.warwick.ac.uk/id/eprint/3313 |

Data sourced from Thomson Reuters' Web of Knowledge

### Actions (login required)

View Item |