Collaborative Research: A posteriori error analysis and adaptivity for discontinuous interface problems

协作研究：后验误差分析和不连续界面问题的自适应性

• 批准号: 1016232
• 负责人: Haiying Wang
• 金额: $ 8万
• 依托单位: Michigan Technological University
• 依托单位国家: 美国
• 项目类别: Standard Grant
• 财政年份: 2010
• 资助国家: 美国
• 起止时间: 2010-10-01 至 2013-09-30
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=1016232&HistoricalAwards=false
• 关键词: Collaborative; Research; posteriori; error; analysis

There are many practical physical situations that can be modeled by partial differential equations in which some of the coefficients in the problem, e.g. those describing material properties, are discontinuous across an interface. Such problems arise in a very broad range of scientific and engineering disciplines, including computational biology, ground water flow and reservoir simulation, environmental remediation studies, crystal growth, wave propagation, sedimentation phenomena, and the preparation of nuclear fuel rods. Such problems cause notorious difficulties for classic numerical methods as a direct result of the lack of smoothness. At the same time, numerical methods for efficiently solving interface problems using a fixed Cartesian grid have attracted considerable attention because they offer a number of computational advantages and a number of approaches have been pursued. Unfortunately, regular-shape discretizations are generally problematic because the interface ?cuts? through the cells without respecting the regular geometry of the discretization, which has a strongly negative impact on the accuracy of the resulting approximations. Consequently, it is critically important to provide computational error estimates that quantify the accuracy of a computed quantity of interest in terms of various sources of discretization and modeling error. This project will develop variational finite element frameworks for several discrete interface methods that identify both a discretization and a modeling component to the model and use the variational framework to derive accurate a posteriori error estimates for specified quantities of interest. Both stationary and evolutionary problems in two and three space dimensions will be considered. The development of efficient adaptive discretization methods for both stationary and evolution problems will be addressed. The methodology and analytic tools to be developed in this proposal will provide a powerful tool for the systematic treatment of problems in which interfaces have complex geometry and the material properties vary considerably on a scale smaller than the overall scale of the discretization. The project will yield a systematic approach to deriving computable and accurate error estimates for quantities of interest and further, provide detailed information about the relative contributions to the error arising from discretization and modeling.

有许多实际的物理情况可以通过偏微分方程进行建模，其中问题中的一些系数，例如那些描述材料属性的数据在界面上是不连续的。这些问题出现在非常广泛的科学和工程学科中，包括计算生物学、地下水流和水库模拟、环境修复研究、晶体生长、波传播、沉降现象以及核燃料棒的制备。由于缺乏平滑性，此类问题给经典数值方法带来了众所周知的困难。与此同时，使用固定笛卡尔网格有效解决界面问题的数值方法引起了相当大的关注，因为它们提供了许多计算优势，并且已经研究了许多方法。不幸的是，规则形状的离散化通常会出现问题，因为界面会“切割”。穿过细胞而不考虑离散化的规则几何形状，这对所得近似值的准确性产生强烈的负面影响。因此，提供计算误差估计至关重要，该估计可以根据离散化和建模误差的各种来源来量化感兴趣的计算量的准确性。该项目将为几种离散接口方法开发变分有限元框架，这些方法识别模型的离散化和建模组件，并使用变分框架为指定的感兴趣量导出准确的后验误差估计。将考虑两个和三个空间维度的稳态和演化问题。将解决针对平稳和演化问题的有效自适应离散化方法的开发。本提案中要开发的方法和分析工具将为系统处理问题提供强大的工具，在这些问题中，界面具有复杂的几何形状，并且材料属性在小于离散化总体规模的规模上变化很大。该项目将产生一种系统方法来导出感兴趣数量的可计算且准确的误差估计，并进一步提供有关离散化和建模所产生的误差的相对贡献的详细信息。