Accurate high performance computing for nonlinear collisional kinetic theory

非线性碰撞动力学理论的精确高性能计算

• 批准号: 1217154
• 负责人: Irene Gamba
• 金额: $ 16.77万
• 依托单位: University of Texas at Austin
• 依托单位国家: 美国
• 项目类别: Standard Grant
• 财政年份: 2012
• 资助国家: 美国
• 起止时间: 2012-09-15 至 2015-08-31
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=1217154&HistoricalAwards=false
• 关键词: Accurate; performance; computing; nonlinear; collisional

This work will implement, analyze, and benchmark accurate numerical schemes for nonlinear collisional kinetic equations and their extension to high performance computing. Kinetic equations describe physical systems at the mesoscopic level, where there are large numbers of interactions between particles but they are not frequent enough to be described as a continuum. At each point in phase space a nonlinear integral term models the effect of interactions with all of the other particles in the system. This integral is difficult and expensive to evaluate due to the delicate conservation structure of the interactions and its dimensionality, however it has many properties that makes it well-suited for parallel computation. This work seeks to overcome the computational bottleneck of the collisional integral term in kinetic models by adapting a conservative spectral numerical method to massively parallel computer architectures, and will investigate how this scales with increasing computational nodes. The proposed work will further enhance the method by applying high order space and time discretization to the transport terms in the system, which was previously infeasible due to the expense of evaluating the collision term, as well as investigating methods for reducing the complexity of the integration for further efficiency. It will also carefully investigate the order of accuracy of computation with regards to the velocity domain cutoff and quadrature. This work will extend the conservative spectral method to the case of anisotropic in angle collisional models, in particular the grazing collisions (Landau) limit and applications to collisional plasma models. Finally, this work will attempt to extend the ideas of asymptotic-preserving schemes to efficiently compute the collision term in stiff regimes.Many important problems in science can be described at the molecular level as individual particles bouncing around, each with its speed and direction. These particles interact with each other as well as other objects through collisions that exchange energy and momentum, and the average behavior of these interactions can be felt as wind, for example. Kinetic models describe problems where particle density is low enough that collisional effects matter in the dynamics of the problem, but is not so strong that one can simply use the average behavior to describe the evolution. Kinetic models can characterize a wide variety of applications such as design of nano- and micro-scale devices, energy applications in plasma physics such as semiconductors and nuclear fusion, and atmospheric re-entry for spacecraft and satellites. In addition to physical applications, these models can be used for the modeling of biological and social dynamics, such as the dynamics of crowds in confined spaces, or modeling the flocking and herding of animals, which could be of security and environmental interest. This work will develop new simulation tools that are able to leverage high performance computing resources to perform high accuracy simulation of problems that were previously computationally infeasible, which can provide a broader view of kinetic applications.

这项工作将实现，分析和基准精确的数值方案的非线性碰撞动力学方程及其扩展到高性能计算。动力学方程描述了介观层次的物理系统，其中粒子之间存在大量的相互作用，但它们不够频繁，不能被描述为连续体。在相空间中的每个点处，非线性积分项模拟与系统中所有其他粒子的相互作用的效果。由于相互作用及其维数的微妙守恒结构，这种积分很难评估，但它有许多属性，使其非常适合并行计算。这项工作旨在克服动力学模型中的碰撞积分项的计算瓶颈，通过将保守的谱数值方法应用于大规模并行计算机架构，并将研究如何随着计算节点的增加而扩展。所提出的工作将进一步加强该方法，通过应用高阶空间和时间离散化的系统中的运输条款，这是以前不可行的，由于费用的碰撞项，以及调查的方法，以降低复杂性的集成，以进一步提高效率。它还将仔细研究关于速度域截止和求积的计算精度的顺序。这项工作将保守谱方法的情况下，各向异性的角度碰撞模型，特别是掠碰撞（朗道）限制和应用碰撞等离子体模型。最后，本文将尝试扩展渐近保持格式的思想，以有效地计算刚性区域中的碰撞项。科学中的许多重要问题可以在分子水平上描述为单个粒子的反弹，每个粒子都有自己的速度和方向。这些粒子通过碰撞相互作用，交换能量和动量，这些相互作用的平均行为可以被感觉为风。动力学模型描述的问题是，粒子密度足够低，碰撞效应在问题的动力学中很重要，但又没有那么强，以至于人们可以简单地使用平均行为来描述演化。动力学模型可以描述各种各样的应用，如纳米和微米级器件的设计，等离子体物理学中的能源应用，如半导体和核聚变，以及航天器和卫星重返大气层。除了物理应用外，这些模型还可用于生物和社会动态的建模，例如密闭空间中人群的动态，或对动物的群集和放牧进行建模，这可能具有安全性和环境利益。这项工作将开发新的模拟工具，能够利用高性能计算资源对以前计算上不可行的问题进行高精度模拟，这可以提供更广泛的动力学应用。