arrow
Volume 27, Issue 1
Bayesian Model Calibration with Interpolating Polynomials Based on Adaptively Weighted Leja Nodes

Laurent van den Bos, Benjamin Sanderse, Wim Bierbooms & Gerard van Bussel

Commun. Comput. Phys., 27 (2020), pp. 33-69.

Published online: 2019-10

Export citation
  • Abstract

An efficient algorithm is proposed for Bayesian model calibration, which is commonly used to estimate the model parameters of non-linear, computationally expensive models using measurement data. The approach is based on Bayesian statistics: using a prior distribution and a likelihood, the posterior distribution is obtained through application of Bayes' law. Our novel algorithm to accurately determine this posterior requires significantly fewer discrete model evaluations than traditional Monte Carlo methods. The key idea is to replace the expensive model by an interpolating surrogate model and to construct the interpolating nodal set maximizing the accuracy of the posterior. To determine such a nodal set an extension to weighted Leja nodes is introduced, based on a new weighting function. We prove that the convergence of the posterior has the same rate as the convergence of the model. If the convergence of the posterior is measured in the Kullback–Leibler divergence, the rate doubles. The algorithm and its theoretical properties are verified in three different test cases: analytical cases that confirm the correctness of the theoretical findings, Burgers' equation to show its applicability in implicit problems, and finally the calibration of the closure parameters of a turbulence model to show the effectiveness for computationally expensive problems.

  • AMS Subject Headings

62F15, 65D0

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address

l.m.m.van.den.bos@cwi.nl (Laurent van den Bos)

b.sanderse@cwi.nl (Benjamin Sanderse)

w.a.a.m.bierbooms@tudelft.nl (Wim Bierbooms)

g.j.w.vanbussel@tudelft.nl (Gerard van Bussel)

  • BibTex
  • RIS
  • TXT
@Article{CiCP-27-33, author = {van den Bos , LaurentSanderse , BenjaminBierbooms , Wim and van Bussel , Gerard}, title = {Bayesian Model Calibration with Interpolating Polynomials Based on Adaptively Weighted Leja Nodes}, journal = {Communications in Computational Physics}, year = {2019}, volume = {27}, number = {1}, pages = {33--69}, abstract = {

An efficient algorithm is proposed for Bayesian model calibration, which is commonly used to estimate the model parameters of non-linear, computationally expensive models using measurement data. The approach is based on Bayesian statistics: using a prior distribution and a likelihood, the posterior distribution is obtained through application of Bayes' law. Our novel algorithm to accurately determine this posterior requires significantly fewer discrete model evaluations than traditional Monte Carlo methods. The key idea is to replace the expensive model by an interpolating surrogate model and to construct the interpolating nodal set maximizing the accuracy of the posterior. To determine such a nodal set an extension to weighted Leja nodes is introduced, based on a new weighting function. We prove that the convergence of the posterior has the same rate as the convergence of the model. If the convergence of the posterior is measured in the Kullback–Leibler divergence, the rate doubles. The algorithm and its theoretical properties are verified in three different test cases: analytical cases that confirm the correctness of the theoretical findings, Burgers' equation to show its applicability in implicit problems, and finally the calibration of the closure parameters of a turbulence model to show the effectiveness for computationally expensive problems.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.OA-2018-0218}, url = {http://global-sci.org/intro/article_detail/cicp/13313.html} }
TY - JOUR T1 - Bayesian Model Calibration with Interpolating Polynomials Based on Adaptively Weighted Leja Nodes AU - van den Bos , Laurent AU - Sanderse , Benjamin AU - Bierbooms , Wim AU - van Bussel , Gerard JO - Communications in Computational Physics VL - 1 SP - 33 EP - 69 PY - 2019 DA - 2019/10 SN - 27 DO - http://doi.org/10.4208/cicp.OA-2018-0218 UR - https://global-sci.org/intro/article_detail/cicp/13313.html KW - Bayesian model calibration, interpolation, Leja nodes, surrogate modeling. AB -

An efficient algorithm is proposed for Bayesian model calibration, which is commonly used to estimate the model parameters of non-linear, computationally expensive models using measurement data. The approach is based on Bayesian statistics: using a prior distribution and a likelihood, the posterior distribution is obtained through application of Bayes' law. Our novel algorithm to accurately determine this posterior requires significantly fewer discrete model evaluations than traditional Monte Carlo methods. The key idea is to replace the expensive model by an interpolating surrogate model and to construct the interpolating nodal set maximizing the accuracy of the posterior. To determine such a nodal set an extension to weighted Leja nodes is introduced, based on a new weighting function. We prove that the convergence of the posterior has the same rate as the convergence of the model. If the convergence of the posterior is measured in the Kullback–Leibler divergence, the rate doubles. The algorithm and its theoretical properties are verified in three different test cases: analytical cases that confirm the correctness of the theoretical findings, Burgers' equation to show its applicability in implicit problems, and finally the calibration of the closure parameters of a turbulence model to show the effectiveness for computationally expensive problems.

Laurent van den Bos, Benjamin Sanderse, Wim Bierbooms & Gerard van Bussel. (2019). Bayesian Model Calibration with Interpolating Polynomials Based on Adaptively Weighted Leja Nodes. Communications in Computational Physics. 27 (1). 33-69. doi:10.4208/cicp.OA-2018-0218
Copy to clipboard
The citation has been copied to your clipboard