2018/08/03 by Borggaard, Jeff, Glatt-Holtz, Nathan E., Krometis, Justin A.
#65N21 (Primary) 62F15 #76R05 (Secondary) #FOS: Mathematics #Numerical Analysis (math.NA)
paper · doi:10.48550/arxiv.1808.01084
We consider the statistical inverse problem of estimating a background flow field (e.g., of air or water) from the partial and noisy observation of a passive scalar (e.g., the concentration of a solute), a common experimental approach to visualizing complex fluid flows. Here the unknown is a vector field that is specified by a large or infinite number of degrees of freedom. Since the inverse problem is ill-posed, i.e., there may be many or no background flows that match a given set of observations, we adopt a Bayesian approach to regularize it. In doing so, we leverage frameworks developed in recent years for infinite-dimensional Bayesian inference. The contributions in this work are threefold. First, we lay out a functional analytic and Bayesian framework for approaching this problem. Second, we define an adjoint method for efficient computation of the gradient of the log likelihood, a key ingredient in many numerical methods. Finally, we identify interesting example problems that exhibit posterior measures with simple and complex structure. We use these examples to conduct a large-scale benchmark of Markov Chain Monte Carlo methods developed in recent years for infinite-dimensional settings. Our results indicate that these methods are capable of resolving complex multimodal posteriors in high dimensions.