Development of Discontinuous Galerkin Methods for Kinetic Transport Models and Control Problems with State Constraints

该项目的目标是开发和分析新的间断Galerkin(DG)方法来求解不同应用领域中的偏微分方程组。DG方法是一类使用完全不连续的分段多项式空间作为数值解和检验函数的有限元方法。这些健壮、紧凑、局部保守的方法可以处理任意非结构网格，是hp自适应策略的理想选择。该格式的良好性质要求在传统的DG方法无法解决的领域进行进一步的研究。在这项拨款提案中，PI计划在以下方向进行研究：(1)求解动力学方程的正性保持DG方法，包括Boltzmann方程和Vlasov方程，(2)所提出的方法在太阳能电池/半导体器件模拟和等离子体物理中的应用，(3)用于Hamilton-Jacobi方程的新型DG求解器及其在状态约束控制问题中的应用。拟议的活动介于算法开发、分析和应用之间。为动力学模型和控制问题开发健壮、高精度、成本效益高的数值算法是非常具有挑战性的，这不仅是因为这些模型的高维，而且还因为需要对基本物理有深刻的理解。最终目标是生成计算效率高且适合应用程序需要的求解器。PI的工作源于现实世界应用程序的计算需求。该提案中提出的许多想法将在半导体器件模拟、高效燃料电池建模、控制问题和等离子体物理中产生直接的应用和影响。PI积极与数学、物理、电气工程和化学系的学生和教职员工互动。此外，国际和平研究所将把该项目与研究生培训结合起来，以便在更广泛的背景下进行交流。

The objective of this project is to develop and analyze novel discontinuous Galerkin (DG) methods for solving partial differential equations arising from various application areas. The DG method is a class of finite element methods using completely discontinuous piecewise polynomial space for the numerical solution and the test functions. Those robust, compact, locally conservative methods can treat arbitrarily unstructured meshes and are ideal for hp-adaptive strategies. The good properties of the scheme call for further research in areas that are traditionally not solved by DG methods. In this grant proposal, the PI plans to conduct research in the following directions: (1) a positivity-preserving DG method for solving the kinetic equations, including the Boltzmann equations and Vlasov equations, (2) application of the proposed method to solar cell/semiconductor device simulations and plasma physics, (3) a novel DG solver for the Hamilton-Jacobi equations and its applications in control problems with state constraints. The proposed activity lies between algorithm development, analysis and applications. Developing robust, high-order accurate, cost-efficient numerical algorithms for kinetic models and control problem is very challenging, not only because of the high dimensionality of such models, but also because of the fact that a deep understanding of the underlying physics is required. The eventual goal is to produce solvers that are computationally efficient and suit the need for applications. The PI's work arises from the computational demand of real world applications. Many ideas developed in this proposal will have straightforward applications and impacts in semiconductor device simulations, high-efficiency fuel cell modeling, control problems and plasma physics. The PI actively interacts with students and faculty members in mathematics, physics, electrical engineering and chemistry departments. In addition, the PI will integrate the project with the training of graduate students in order to communicate in a broader context.