AF: Small: Fully Implicit Methods for Partial Differential Equations and Software for Hybrid Architecture

该项目的第一个目标是开发新的区域分解算法和软件，用于多物理应用中产生的一些高度非线性、耦合的偏微分方程系统的数值解。 对于简单问题，隐式方法相对容易从给定的显式或半隐式方法发展，但对于一些多物理场问题，很难发展一个允许高性能实现的全隐式方法。 大多数现有的技术是非光滑的，因此很难用牛顿型求解器求解。在这个项目中，将开发一些离散化技术，包括用于捕获域信息的高阶非光滑离散化，以及用于构建牛顿迭代的雅可比系统的低阶局部光滑离散化。低阶离散作为一个非线性预条件，加快收敛速度，但不改变解决方案的精度。该项目的第二个目标是开发一个有效的实现建议的线性和非线性预处理方法的高性能计算机上有大量的处理器。对于许多类型的单物理场问题，已有相对成熟的算法和软件技术，但对于多物理场耦合问题，特别是对于大规模并行计算机，需要健壮的、可扩展的技术。所提出的算法和软件将对几个重要的应用领域产生重大影响，如全球大气流动和生物流体的模拟，也将对需要求解大型线性和非线性方程的计算科学的其他领域产生重大影响。为了扩大研究的影响，该软件将与广泛使用的PETSc软件包完全兼容。该研究为对高性能计算和一般计算科学与工程感兴趣的研究生和本科生提供了丰富的机会。

The first goal of the project is to develop new domain decomposition algorithms and software for the numerical solution of some highly nonlinear, coupled systems of partial differential equations arising from multi-physics applications. For simple problems, implicit methods are relatively easy to develop from a given explicit or semi-implicit method, but for some multi-physics problems it is quite difficult to develop a fully implicit method that allows a high performance implementation. Most of the existing techniques are non-smooth and therefore difficult to solve with Newton type solver. In the project, some discretization techniques will be developed that involve high order non-smooth discretization for capturing the domain information, and low order locally smooth discretization for building the Jacobian system of the Newton iterations. The low order discretization serves as a nonlinear preconditioner that speeds up the convergence, but doesn't change the accuracy of the solution. The second goal of the project is to develop an efficient implementation of the proposed linear and nonlinear preconditioning approaches on high performance computers with a large number of processors. Relatively mature technologies including algorithms and software are available for solving many types of single physics problems, but for coupled multi-physics problems, robust and scalable techniques are badly needed, especially for large scale parallel computers with accelerators. The proposed algorithms and software will have a great impact on several important application areas, such as the simulation of global atmospheric flows and the bio-fluids, and will also have substantial influence on other areas of computational sciences where large linear and nonlinear equations need to be solved. To broaden the impact of the research, the software will be made fully compatible with the widely used PETSc package. The research is a rich area in opportunities for both graduate and undergraduate students interested in high performance computing and general Computational Science and Engineering.