Multifidelity Monte Carlo methods often rely on a preprocessing phase consisting of standard Monte Carlo sampling to estimate correlation coefficients between models of different fidelity to determine the weights and number of samples for each level. For computationally intensive models, as are often encountered in simulations of chaotic systems, this up-front cost can be prohibitive. In this work, a correlation estimation procedure is developed for the case in which the highest and next highest fidelity models are generated via discretizing the same mathematical model using different resolution. The procedure uses discretization error estimates to estimate the required correlation coefficient without the need to sample the highest fidelity model, which can dramatically decrease the cost of the preprocessing phase. The method is extended to chaotic problems by using discretization error estimates that account for the statistical nature of common quantities of interest and the accompanying finite sampling errors that pollute estimates of such quantities of interest. The methodology is then demonstrated on a model problem based on the Kuramoto-Sivashinsky equation.
翻译:Monte Carlo 方法往往依赖由标准 Monte Carlo 抽样抽样组成的预处理阶段来估计不同忠诚模型之间的相关系数,以确定每一等级的重量和样品数量。对于计算密集模型,正如在模拟混乱系统时经常遇到的那样,这种前期成本可能令人望而却步。在这项工作中,为以下情况制定了一个相关估计程序:最高和下一级忠诚模型是通过使用不同分辨率分离同一数学模型生成的。该程序使用离散错误估计来估计所需的相关系数,而不必采样最高忠诚模型,这可以大大降低预处理阶段的成本。这种方法通过使用离散错误估计方法扩大到混乱问题,即考虑到共同利益的统计性质以及伴随的有限抽样错误,从而对此类利益的数量进行估算。然后根据Kuramoto-Sivashinsky 方程式的模型问题来证明该方法。