Mathematical Sciences: Numerical Solution of Differential Equations in Mechanics

数学科学：力学微分方程的数值解

• 批准号: 9500672
• 负责人: Douglas Arnold
• 金额: $ 12.4万
• 依托单位: Pennsylvania State Univ University Park
• 依托单位国家: 美国
• 项目类别: Standard Grant
• 财政年份: 1995
• 资助国家: 美国
• 起止时间: 1995-07-01 至 1998-06-30
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=9500672&HistoricalAwards=false
• 关键词: Mathematical; Sciences; Numerical; Solution; Differential

9500672 Arnold This project will involve devising, improving, and analyzing methods for the numerical simulation of complex physical phenomena modeled by partial differential equations. A major emphasis will be on the modeling of thin structures (plates and shells). Specifically, an innovative new finite element method recently proposed by a team of computational engineers for plate bending will be analyzed and a new approach the investigator introduced in collaboration with Brezzi to develop and analyze new finite element methods for modeling the deformation of elastic shells will be pursued. This latter approach has produced the best theoretical results to date for avoiding the phenomenon known as membrane locking which has plagued shell computations. In addition to studying numerical methods and the errors they introduce, the investigator will consider the errors introduced by mathematical models of thin structures. Specifically, a new technique of analysis to study the convergence of two-dimensional models of thin plates will be used. This technique gives the possibility of obtaining sharp convergence estimates for plate modeling errors both globally and away from the lateral boundary (where convergence may be higher). It is hoped that this research will clarify in what situations the Reissner-Mindlin plate model affords better approximation than the classical Kirchhoff model, an important question that eluded resolution for decades. A second area of emphasis will be on the development, implementation, validation, and dissemination of a finite element code to simulate the collision of black holes. This problem involves the solution of the Einstein field equations on an unbounded domain in four dimensional space-time and is extremely demanding computationally. Work will focus, at least initially, on the solution of the initial data problem for black hole collisions, i.e., the determination of physical ly consistent and meaningful initial data. The initial data problem is a nonlinear elliptic system posed on an unbounded three-dimensional domain. For this problem the computer code to be produced will involve a variety of sophisticated modern techniques that have not yet been applied in computational relativity, including adaptive tetrahedral mesh generation and multigrid preconditioning. The investigator's work is directly related to several critical technological issues. It is generally agreed by practitioners that the current state of the art of numerical modeling of shells is inadequate, and this has been implicated in the failure of engineering structures such as stadium roofs and airplane fuselages. The investigator's work in both the area of the development and analysis of finite element methods for plates and shells, and on plate modeling errors is directed at overcoming these problems. The work on the simulation of black hole collisions is directly motivated by the ongoing construction of the LIGO gravitational wave observatory. LIGO is expected to be able to detect radiation from black hole collisions early in the next decade. Massive numerical simulations are needed to identify the cosmological events which are the sources for the detected radiation. It is clear that these simulations push the limits of available hardware and software technology. The investigator's work is directed at improving the speed and accuracy of these simulations by employing and adapting recent advances in numerical algorithms. The resulting computer codes will be adapted to parallel computer architectures and will provide an excellent opportunity for the graduate students and postdoctoral fellows working under the supervision of the investigator to learn high performance computing techniques.

这个项目将涉及设计、改进和分析用偏微分方程模拟的复杂物理现象的数值模拟方法。主要的重点将放在薄结构（板和壳）的建模上。具体来说，将分析一组计算工程师最近提出的一种用于板弯曲的创新有限元方法，以及研究人员与Brezzi合作开发和分析用于模拟弹性壳变形的新有限元方法的新方法。后一种方法产生了迄今为止最好的理论结果，避免了困扰壳计算的称为膜锁定的现象。除了研究数值方法及其引入的误差外，研究者还将考虑薄结构数学模型引入的误差。具体地说，将采用一种新的分析技术来研究薄板二维模型的收敛性。这种技术提供了对全球和远离横向边界（收敛可能更高）的板块建模误差获得尖锐收敛估计的可能性。希望这项研究能够澄清在什么情况下Reissner-Mindlin板模型比经典Kirchhoff模型提供更好的近似，这是一个几十年来未能解决的重要问题。第二个重点领域将是开发、实现、验证和传播模拟黑洞碰撞的有限元代码。这个问题涉及四维时空无界域上爱因斯坦场方程的解，计算要求极高。至少在最初阶段，工作将集中于解决黑洞碰撞的初始数据问题，即确定物理上一致和有意义的初始数据。初始数据问题是一个无界三维区域上的非线性椭圆系统。对于这个问题，生成的计算机代码将涉及各种尚未在计算相对论中应用的复杂现代技术，包括自适应四面体网格生成和多网格预处理。研究者的工作直接关系到几个关键的技术问题。从业人员普遍认为，目前的炮弹数值模拟技术是不充分的，这已经涉及到工程结构的失败，如体育场屋顶和飞机机身。研究者在板壳有限元方法的开发和分析领域的工作，以及在板建模误差方面的工作，都是为了克服这些问题。模拟黑洞碰撞的工作直接受到正在进行的LIGO引力波天文台建设的推动。LIGO有望在未来十年内探测到黑洞碰撞产生的辐射。需要大量的数值模拟来确定作为探测到的辐射源的宇宙事件。很明显，这些模拟突破了现有硬件和软件技术的极限。研究者的工作是通过采用和适应数值算法的最新进展来提高这些模拟的速度和准确性。由此产生的计算机代码将适用于并行计算机体系结构，并将为在研究者监督下工作的研究生和博士后提供学习高性能计算技术的绝佳机会。