Error Analysis and Uncertainty Quantification for the Heterogeneous Transport Equation in Slab Geometry
IMA journal of numerical analysis(2020)
摘要
We present an analysis of multilevel Monte Carlo (MLMC) techniques for the forward problem of uncertainty quantification for the radiative transport equation, when the coefficients (cross-sections) are heterogenous random fields. To do this we first give a new error analysis for the combined spatial and angular discretisation in the deterministic case, with error estimates that are explicit in the coefficients (and allow for very low regularity and jumps). This detailed error analysis is done for the one-dimensional space-one-dimensional angle slab-geometry case with classical diamond differencing. Under reasonable assumptions on the statistics of the coefficients, we then prove an error estimate for the random problem in a suitable Bochner space. Because the problem is not self-adjoint, stability can only be proved under a path-dependent mesh resolution condition. This means that, while the Bochner space error estimate is of order O(h(eta)) for some eta where h is a (deterministically chosen) mesh diameter, smaller mesh sizes might be needed for some realisations. We also show that the expected cost for computing a typical quantity of interest remains of the same order as for a single sample. This leads to rigorous complexity estimates for Monte Carlo (MC) and MLMC: for particular linear solvers, the multilevel version gives up to two orders of magnitude improvement over MC. We provide numerical results supporting the theory.
更多查看译文
关键词
radiative transport,neutron transport,spatial heterogeneity,random coefficients,error estimate,multilevel Monte Carlo
AI 理解论文
溯源树
样例
生成溯源树,研究论文发展脉络
Chat Paper
正在生成论文摘要