Collaborative Research: Efficient Coupling of Multilevel Partial Differential Equation Solvers and Advanced Sampling Methods

今天，我们知道如何在计算机上模拟许多过程，例如机翼周围的空气流动或物体在受力时如何变形。然而，从实践的角度来看，通常需要确定我们预测的“量化不确定性”，即，当输入不完全已知时，能够说出计算机模拟的准确性或不准确性。一个例子是当前面的空气是湍流时翼型的行为。在这些情况下，需要模拟许多场景，这是一个非常昂贵的提议。该项目旨在开发数学工具，使这一过程更加可行和计算负担得起的利用多级层次的采样和解决方案的approach.Many过程在科学和工程描述的参数依赖偏微分方程（PDE）的数值解是昂贵的。虽然它们的正解可能是可行的，但计算成本使得难以或不可能统计地评估问题，例如（i）在给定参数值的已知或假设分布的情况下量化输出变量中的不确定性，或（ii）解决从噪声测量推断参数之间的概率分布的统计逆问题。这两个问题通常通过高维的，因此非常昂贵的采样方法来回答。该项目将使用PDE求解器的分层和多级分解，加上误差估计，以获得更快的方法来获得足够数量的样本，以实现规定精度的统计抽样。该方法意味着，只要有可能，样本和PDE解决方案可以发生在其各自的层次结构的粗糙，廉价的水平，昂贵的解决方案只在必要时执行。因此，这些层次结构的结合使用有望以更低的成本提供更多、更丰富的样本，从而实现以前不可能实现的科学和工程应用。该奖项反映了NSF的法定使命，并通过使用基金会的知识价值和更广泛的影响审查标准进行评估，被认为值得支持。

We know today how to simulate many processes on computers, such as air flow around airfoils or how objects deform when a force is applied. However, from a practical perspective, it is often desirable to determine to "quantify the uncertainty" of our predictions, i.e., to be able to say how accurate or inaccurate a computer simulation is when the inputs are not exactly known. An example is how an airfoil behaves when the air ahead is turbulent. In these cases, one needs to simulate many scenarios, an exceedingly expensive proposition. This project aims to develop mathematical tools to make this process more feasible and computationally affordable by taking advantage of multilevel hierarchies in the sampling and solution approach.Many processes in the sciences and engineering are described by parameter-dependent partial differential equations (PDEs) whose numerical solution is expensive. While their forward solution may be feasible, the computational cost makes it difficult or impossible to statistically evaluate questions such as (i) to quantify the uncertainty in output variables given a known or assumed distribution of parameter values, or (ii) to solve the statistical inverse problem of inferring a probability distribution among parameters from noisy measurements. Both of these questions are typically answered through high-dimensional and consequently very expensive sampling methods. This project will use hierarchical and multilevel decompositions of PDE solvers, coupled with error estimates, to derive vastly faster methods to obtain a sufficient number of samples to achieve statistical sampling with prescribed accuracy. The approach implies that whenever possible, samples and PDE solutions can happen on coarse, cheap levels of their respective hierarchies, and expensive solves are performed only when necessary. The combined use of these hierarchies therefore promises to deliver more, and more informative, samples at a cheaper cost - thereby enabling applications in the sciences and engineering that were not possible before.This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.