Towards Petascale Simulation of Urban Earthquake Impacts

该奖项是NSF 07 - 559计划征集的结果，“通过Petascale模拟和分析加速科学和工程发现”竞赛。 领导机构是卡内基梅隆大学，分奖项给加州大学伯克利分校、戴维斯分校和圣地亚哥分校。该奖项将开发用于整个高地震活动城市区域（大洛杉矶盆地（GLAB））的高保真、基于物理的千万亿次模拟的方法、能力和软件，以评估大震级地震对建筑物、交通系统组件和地下民用基础设施的工程影响。 地震是美国面临的最严重的自然灾害之一。 模拟它们对构建环境的影响是千万亿次高性能计算（HPC）系统的独特科学驱动力，因为具有不同物理特性的多个长度尺度，大数据量，以及对高度可扩展的并行可视化和数据查询的需求。目前的三维地震模拟仅限于线性土壤特性和孤立的单个结构。超越目前的方法，该项目将开发模拟能力，包括土壤，基础和大型库存的结构之间的相互作用;密集建筑区的结构之间的相互作用;和结构对自由场运动的影响;以及非线性土壤行为和饱和土壤中导致液化的孔隙水压力的影响。 目前的HPC能力不足以进行这些复杂的模拟。该项目将使用新兴的千万亿次系统来执行和集成（1）大型沉积盆地的地面运动模拟，（2）土壤的非线性行为模拟，（3）建筑物，桥梁和其他基础设施系统的大型库存模拟，（4）计算数据库，以及（5）可扩展的可视化技术。 该项目将在千万亿次模拟的分层，多尺度方法学方面取得新的进展。从一个大的区域，宽带地震模拟到多个子区域，将允许以高度可扩展的方式进行详细建模，以准确地捕捉整个库存的场地响应和结构响应。域缩减方法将首次使用区域模拟的地面运动作为多个人口密集子区域的输入。 这些子区域包括高度非线性土壤模型和数百个建筑物和桥梁的详细模型，沿着嵌入式基础。将开发一种新的可扩展隐式-显式时间积分方法，为一个分区的复杂地震模拟提供最佳计算性能。 数据分析和可视化能力将被开发为在与模拟相同的并行处理器上运行，从而大大减少移动数据的需求。使用这种方法进行数据密集型超级计算，现场可视化和新的计算数据库系统将使人们能够前所未有地了解地震影响。新功能将提供有助于理解和决策的信息，通过深入到感兴趣的细节，同时保持全球背景。GLAB之所以被选中，是因为它具有很高的地震风险，重要的公共政策需求，以及它为验证和确认提供的丰富信息。 GLAB测试台的年度里程碑经过校准，以尽早利用HPC资源，并随着千万亿次系统的上线而扩展模拟。 模拟的方法，应用和结果将有助于灾害规划和管理，因为它是一个城市系统如何在大地震，需要改善备灾和减灾的详细知识。 从事该项目的研究生将有机会通过土木工程，计算工程和计算机科学的多学科研究创造新知识。

This award is an outcome of the NSF 07-559 program solicitation, "accelerating Discovery in Science and Engineering through Petascale Simulations and Analysis" competition. The lead institution is Carnegie Mellon University, with subawards to the Universities of California at Berkeley, Davis, and San Diego. This award will develop methodologies, capability, and software for high-fidelity, physics-based petascale simulations of an entire high seismicity urban region, the Greater Los Angeles Basin (GLAB), to assess the engineering impacts of large magnitude earthquakes on buildings, transportation system components, and the underground civil infrastructure. Earthquakes are one of the most severe natural hazards facing the United States. Simulating their effects on the built environment is a unique science driver for petascale high-performance computing (HPC) systems because of the multiple length scales with different physics, large data volumes, and need for highly scalable parallel visualization and data querying. Current three-dimensional earthquake simulations have been limited to linear soil behavior and isolated, individual structures. Advancing beyond current approaches, this project will develop simulation capability that includes the interaction between the soil, foundations, and large inventories of structures; the interaction between structures in densely built areas; and the effect of the structures on the free-field motion; as well as the nonlinear soil behavior and the effect of the pore water pressure in saturated soils leading to liquefaction. Current HPC capability is inadequate for these complex simulations. This project will use emerging petascale systems to perform and integrate (1) ground motion simulation of large sedimentary basins, (2) simulation of the nonlinear behavior of soil, (3) simulation of large inventories of buildings, bridges, and other infrastructure systems, (4) computational databases, and (5) scalable visualization techniques. The project will make new advances in a hierarchical, multi-scale methodology for petascale simulation. Going from a large regional, broadband earthquake simulation to multiple subregions will allow detailed modeling in a highly scalable manner to capture site response and structural response accurately for entire inventories. The Domain Reduction Method will be applied for the first time using ground motion from a regional simulation as input for multiple highly populated subregions. These subregions include models of highly nonlinear soil and detailed models of hundreds of buildings and bridges along with embedded foundations. A new scalable implicit-explicit time integration method will be developed to provide optimal computational performance for the complex earthquake simulations for a subregion. Data analysis and visualization capability will be developed to run on the same parallel processors as the simulation, drastically reducing the need to move data. Using this approach for data-intensive supercomputing, in-situ visualization and a new computational database system will allow unprecedented ability to understand earthquake impacts. The new capability will present information that facilitates understanding and decision-making by drilling down to the detail of interest while maintaining the global context. The GLAB was selected because of the high seismic risk, important public policy need, and the wealth of information that it provides for verification and validation. The annual milestones for the GLAB test bed are calibrated to take early advantage of HPC resources, with scaling of the simulations to take place as petascale systems are brought online. The methodology, applications, and the results of the simulations will be useful for disaster planning and management, since it is the detailed knowledge of how an urban system performs in a large earthquake that is needed for improving disaster preparedness and mitigation. Graduate students working on the project will have the opportunity to create new knowledge through multidisciplinary research in civil engineering, computational engineering, and computer science.