Skip to content



Adaptive mesh refinement for a 2D unified continuous/discontinuous Galerkin Non-hydrostatic Atmospheric Model

Michal A. Kopera, ; Frank X. Giraldo, (Naval Postgraduate School)
Thursday 25 October 2012, 10:30-10:55



The adaptive mesh refinement techniques for element-based Galerkin methods are becoming a strong candidate for future numerical weather prediction models. Particular attention has been paid to the discontinuous Galerkin method [1], [2], [3] as it avoids global assembly of data and makes the implementation of the algorithm easier. In this presentation we will focus on the extension of the 2D discontinuous Galerkin, quad-based non-conforming adaptive mesh refinement algorithm to a continuous Galerkin formulation. The novelty of this approach is that we propose to do this within a unified CG/DG nonhydrostatic atmospheric model that we call NUMA (Nonhydrostatic Unified Model of the Atmosphere). NUMA is equipped to handle AMR at various levels: IMEX time-integrators are used to be able to use large time-steps and a new class of preconditioners [4] have been specifically designed to handle the IMEX methods with AMR.

[1] A. Muller, J. Behrens, F.X. Giraldo, V. Wirth (2011). An Adaptive Discontinuous Galerkin Method for Modelling Atmospheric Convection. Defense Technical Information Center Report,

[2] S. Blaise, and A. St-Cyr (2011). A Dynamic hp-Adaptive Discontinuous Galerkin Method for Shallow-Water Flows on the Sphere with Application to a Global Tsunami Simulation, Monthly

[3] M. A. Kopera and F.X. Giraldo (2012). AMR for a 2d DG Nonhydrostatic atmospheric model, in preparation.

[4] L.E. Carr, C.F. Borges, and F.X. Giraldo (2012). An element-based spectrally-optimized approximate inverse preconditioner for the Euler equations, SIAM J. Sci. Comp. (in press).


[zip ] [pdf ]

Back to top ∧