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
RightsThis publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.
Showing items related by title, author, creator and subject.
Multigrid preconditioners for the hybridised discontinuous Galerkin discretisation of the shallow water equations Gibson, Thomas H.; Müller, Eike H.; Betteridge, Jack; Graham, Ivan G. (Elsevier, 2020-10-26);Numerical climate- and weather-prediction models require the fast solution of the equations of fluid dynamics. Discontinuous Galerkin (DG) discretisations have several advantageous properties. They can be used for arbitrary ...
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 ...
Mugg, Patrick (Monterey, CA; Naval Postgraduate School, 2021-09);Within atmospheric earth system models, accurately resolving the various physical multi-scale processes is critical. For example, the single-layered shallow water equations are widely used to model ocean waves in shallow ...