IMIS

Publications | Institutes | Persons | Datasets | Projects | Maps
[ report an error in this record ]basket (0): add | show Print this page

A split-explicit second order Runge-Kutta method for solving 3D hydrodynamic equations
Ishimwe, A.P.; Deleersnijder, E.; Legat, V.; Lambrechts, J. (2023). A split-explicit second order Runge-Kutta method for solving 3D hydrodynamic equations. Ocean Modelling 186: 102273. https://dx.doi.org/10.1016/j.ocemod.2023.102273
In: Ocean Modelling. Elsevier: Oxford. ISSN 1463-5003; e-ISSN 1463-5011, more
Peer reviewed article  

Available in  Authors 

Keyword
    Marine/Coastal
Author keywords
    Unstructure mesh; Split-explicit Runge-Kutta; Hydrodynamic equations; Timestep methods

Authors  Top 
  • Ishimwe, A.P., more
  • Deleersnijder, E., more
  • Legat, V., more
  • Lambrechts, J., more

Abstract
    Numerical models of marine hydrodynamics have to deal with processes exhibiting a wide range of timescales. These processes include fast external gravity waves and slower internal fully three-dimensional motions. In order to be both time-efficient and numerically stable, the time stepping scheme has to be chosen carefully to cope with the characteristic time scale of each phenomenon. An usual approach is to split the fast and slow dynamics into separate modes. The fast waves are modeled with a two-dimensional system through depth averaging while the other motions, where characteristic times are much longer, are dealt with in three dimensions. However, if the splitting is inexact, for instance in projecting the fields in a new 3D mesh, this procedure can lead to improper results regarding to the physical properties such as mass conservation and tracer consistency. In this work, a new split-explicit Runge-Kutta scheme is adapted and developed for the Discontinuous-Galerkin Finite Element method in order to obtain a new second-order time stepping, yielding to more accurate results. This method combines a three-stage low-storage Runge-Kutta for the slow processes and a two-stage low-storage for the fast ones. The 3D iterations are not affecting the surface elevation, hence an Arbitrary Lagrangian Eulerian implementation is straightforward. Water volume and tracers are conserved. A set of test cases for baroclinic flows as well as a laboratory application demonstrate the performance of the scheme. They suggest that the new scheme has little numerical diffusion.

All data in the Integrated Marine Information System (IMIS) is subject to the VLIZ privacy policy Top | Authors