Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: Limited-area mode

Continuous and discontinuous Galerkin methods for a scalable three-dimensional nonhydrostatic atmospheric model: Limited-area mode
复制标题

用于可扩展三维非静水压大气模型的连续和不连续伽辽金方法:有限区域模式

DOI:
10.1016/j.jcp.2012.04.042
复制
发表时间:
2012
期刊:
J. Comput. Phys.
影响因子:
--
通讯作者:
F. Giraldo
F. Giraldo
中科院分区:
--
文献类型:
--
作者:
J. Kelly;F. Giraldo

文献摘要

被引文献

相似文献

本文介绍了一个统一的,基于元素的伽辽金(EBG)框架的三维,非静力大气模式。EBG方法具有高阶精度、几何灵活性、优良的色散特性和良好的可扩展性。我们的非静力模式,可压缩的欧拉方程的基础上,是适当的有限区域和全球大气模拟。连续伽辽金(CG),或谱元素,和不连续伽辽金(DG)模型被认为是使用六面体元素。该公式适用于全球和有限区域的大气建模,虽然我们限制我们的注意力在这项研究中的3D有限区域的现象,全球大气模拟将在后续文件。我们的CG和DG模型使用的区域分解和通信算法。对CG和DG的通信量和交换算法进行了比较。该模型的数值验证进行了两个测试用例:流过去的3D山和浮力对流的气泡在中性大气中,这些测试表明,CG和DG可以模拟必要的物理干燥的大气动力学。这两种方法的可扩展性显示高达8192个CPU核心,DG接近理想的可扩展性高达32,768个核心。
This paper describes a unified, element based Galerkin (EBG) framework for a three-dimensional, nonhydrostatic model for the atmosphere. In general, EBG methods possess high-order accuracy, geometric flexibility, excellent dispersion properties and good scalability. Our nonhydrostatic model, based on the compressible Euler equations, is appropriate for both limited-area and global atmospheric simulations. Both a continuous Galerkin (CG), or spectral element, and discontinuous Galerkin (DG) model are considered using hexahedral elements. The formulation is suitable for both global and limited-area atmospheric modeling, although we restrict our attention to 3D limited-area phenomena in this study; global atmospheric simulations will be presented in a follow-up paper. Domain decomposition and communication algorithms used by both our CG and DG models are presented. The communication volume and exchange algorithms for CG and DG are compared and contrasted. Numerical verification of the model was performed using two test cases: flow past a 3D mountain and buoyant convection of a bubble in a neutral atmosphere; these tests indicate that both CG and DG can simulate the necessary physics of dry atmospheric dynamics. Scalability of both methods is shown up to 8192 CPU cores, with near ideal scaling for DG up to 32,768 cores.