No abstract

No abstract

No abstract

Abstract. Spectral deferred correction (SDC) methods for solving ordinary differential equations (ODEs) were introduced by Dutt, Greengard and Rokhlin (2000). It was shown in that paper that SDC methods can achieve arbitrary high order accuracy and possess nice stability properties. Their SDC methods are constructed with low order integrators, such as forward Euler or backward Euler, and are able to handle stiff and non-stiff terms in the ODEs. In this paper, we use high order Runge-Kutta (RK) integrators to construct a family of related methods, which we refer to as integral deferred correction (IDC) methods. The distribution of quadrature nodes is assumed to be uniform, and the corresponding local error analysis is given. The smoothness of the error vector associated with an IDC method, measured by the discrete Sobolev norm, is a crucial tool in our analysis. The expected order of accuracy is demonstrated through several numerical examples. Superior numerical stability and accuracy regions are observed when high order RK integrators are used to construct IDC methods.

Abstract. In this paper, we develop a class of high order conservative semi-Lagrangian (SL) discontinuous Galerkin (DG) methods for solving multi-dimensional linear transport equations. The methods rely on a characteristic Galerkin weak formulation, leading to L 2 stable discretizations for linear problems. Unlike many existing SL methods, the high order accuracy and mass conservation of the proposed methods are realized in a non-splitting manner. Thus, the detrimental splitting error, which is known to significantly contaminate long term transport simulations, will be not incurred. One key ingredient in the scheme formulation, borrowed from CSLAM [Lauritzen, Nair & Ullrich, 2010], is the use of Green's theorem which allows us to convert volume integrals into a set of line integrals. The resulting line integrals are much easier to approximate with high order accuracy, hence facilitating the implementation. Another novel ingredient is the construction of quadratic curves in approximating sides of upstream cell, leading to quadratic-curved quadrilateral upstream cells. Formal third order accuracy is obtained by such a construction. The desired positivity-preserving property is further attained by incorporating a high order bound-preserving filter. To assess the performance of the proposed methods, we test and compare the numerical schemes with a variety of configurations for solving several benchmark transport problems with both smooth and nonsmooth solutions. The efficiency and efficacy are numerically verified.

The discontinuous Galerkin (DG) methods designed for hyperbolic problems arising from a wide range of applications are known to enjoy many computational advantages. DG methods coupled with strong-stability-preserving explicit Runge–Kutta discontinuous Galerkin (RKDG) time discretizations provide a robust numerical approach suitable for geoscience applications including atmospheric modeling. However, a major drawback of the RKDG method is its stringent Courant–Friedrichs–Lewy (CFL) stability restriction associated with explicit time stepping. To address this issue, the authors adopt a dimension-splitting approach where a semi-Lagrangian (SL) time-stepping strategy is combined with the DG method. The resulting SLDG scheme employs a sequence of 1D operations for solving multidimensional transport equations. The SLDG scheme is inherently conservative and has the option to incorporate a local positivity-preserving filter for tracers. A novel feature of the SLDG algorithm is that it can be used for multitracer transport for global models employing spectral-element grids, without using an additional finite-volume grid system. The quality of the proposed method is demonstrated via benchmark tests on Cartesian and cubed-sphere geometry, which employs nonorthogonal, curvilinear coordinates.

No abstract

No abstract