High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model
Giraldo. Francis X.
MetadataShow full item record
We extend the explicit in time high-order triangular discontinuous Galerkin (DG) method to semi-implicit (SI) and then apply the algorithm to the two-dimensional oceanic shallow water equations; we implement high-order SI time-integrators using the backward difference formulas from orders one to six. The reason for changing the time-integration method from explicit to SI is that explicit methods require a very small time step in order to maintain stability, especially for high-order DG methods. Changing the time- integration method to SI allows one to circumvent the stability criterion due to the gravity waves, which for most shallow water applications are the fastest waves in the system (the exception being supercritical flow where the Froude number is greater than one). The challenge of constructing a SI method for a DG model is that the DG machinery requires not only the standard finite element-type area integrals, but also the finite volume-type boundary integrals as well. These boundary integrals pose the biggest challenge in a SI discretization because they require the construction of a Riemann solver that is the true linear representation of the nonlinear Riemann problem; if this condition is not satisfied then the resulting numerical method will not be consistent with the continuous equations. In this paper we couple the SI time-integrators with the DG method while maintaining most of the usual attributes associated with DG methods such as: high-order accuracy (in both space and time), parallel efficiency, excellent stability, and conservation. The only property lost is that of a compact communication stencil typical of time-explicit DG methods; implicit methods will always require a much larger communication stencil. We apply the new high-order SI DG method to the shallow water equations and show results for many standard test cases of oceanic interest such as: standing, Kelvin and Rossby soliton waves, and the Stommel problem. The results show that the new high-order SI DG model, that has already been shown to yield exponentially convergent solutions in space for smooth problems, results in a more efficient model than its explicit counterpart. Furthermore, for those problems where the spatial resolution is sufficiently high compared with the length scales of the flow, the capacity to use high-order (HO) time-integrators is a necessary complement to the employment of HO space discretizations, since the total numerical error would be otherwise dominated by the time discretization error. In fact, in the limit of increasing spatial resolution, it makes little sense to use HO spatial discretizations coupled with low-order time discretizations.
The Discontinuous Galerkin Coastal Ocean Model (DGCOM) is a two-dimensional shallow water model for simulating coastal ocean processes.The article of record may be found at http://www.doi.org/10.1002/fld.2118
Showing items related by title, author, creator and subject.
Hood, Karoline (Monterey, California: Naval Postgraduate School, 2016-06);Storm surges have a devastating impact on coastlines throughout the United States. In order to accurately understand the impacts of storm surges there needs to be an effective model. One of the governing systems of equations ...
Acceleration of the implicit-explicit non-hydrostatic unified model of the atmosphere (NUMA) on Manycore processors Abdi, Daniel S.; Giraldo, Francis X.; Constantinescu, Emil M.; Carr, Lester E., III; Wilcox, Lucas C.; Warburton, Timothy C. (Sage, 2017);We present the acceleration of an IMplicit-EXplicit (IMEX) non-hydrostatic atmospheric model on manycore processors such as GPUs and Intel’s MIC architecture. IMEX time integration methods sidestep the constraint imposed ...
Giraldo, F.X. (2000);The weak Lagrange-Galerkin finite element method fur the 2D shallow water equations on the sphere is presented. This method offers stable and accurate solutions because the equations are integrated along thr characteristics. ...