Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 4.252 IF 4.252
  • IF 5-year value: 4.890 IF 5-year 4.890
  • CiteScore value: 4.49 CiteScore 4.49
  • SNIP value: 1.539 SNIP 1.539
  • SJR value: 2.404 SJR 2.404
  • IPP value: 4.28 IPP 4.28
  • h5-index value: 40 h5-index 40
  • Scimago H index value: 51 Scimago H index 51
Volume 11, issue 3
Geosci. Model Dev., 11, 1161-1179, 2018
https://doi.org/10.5194/gmd-11-1161-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Special issue: Modelling lakes in the climate system (GMD/HESS inter-journal...

Geosci. Model Dev., 11, 1161-1179, 2018
https://doi.org/10.5194/gmd-11-1161-2018
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

Development and technical paper 29 Mar 2018

Development and technical paper | 29 Mar 2018

A fully consistent and conservative vertically adaptive coordinate system for SLIM 3D v0.4 with an application to the thermocline oscillations of Lake Tanganyika

Philippe Delandmeter1,2, Jonathan Lambrechts1, Vincent Legat1, Valentin Vallaeys1, Jaya Naithani1, Wim Thiery3,4, Jean-François Remacle1, and Eric Deleersnijder5,6 Philippe Delandmeter et al.
  • 1Université catholique de Louvain, Institute of Mechanics, Materials and Civil Engineering (IMMC), Avenue Georges Lemaître 4, 1348 Louvain-la-Neuve, Belgium
  • 2Utrecht University, Institute for Marine and Atmospheric Research, Princetonplein 5, 3584 CC Utrecht, the Netherlands
  • 3ETH Zürich, Institute for Atmospheric and Climate Sciences, Universitätstrasse 16, 8092 Zürich, Switzerland
  • 4Vrije Universiteit Brussel, Department of Hydrology and Hydraulic Engineering, Pleinlaan 2, 1050 Brussels, Belgium
  • 5Université catholique de Louvain, Institute of Mechanics, Materials and Civil Engineering (IMMC) & Earth and Life Institute (ELI), Avenue Georges Lemaître 4, 1348 Louvain-la-Neuve, Belgium
  • 6Delft University of Technology, Delft Institute of Applied Mathematics (DIAM), Mekelweg 4, 2628 CD Delft, the Netherlands

Abstract. The discontinuous Galerkin (DG) finite element method is well suited for the modelling, with a relatively small number of elements, of three-dimensional flows exhibiting strong velocity or density gradients. Its performance can be highly enhanced by having recourse to r-adaptivity. Here, a vertical adaptive mesh method is developed for DG finite elements. This method, originally designed for finite difference schemes, is based on the vertical diffusion of the mesh nodes, with the diffusivity controlled by the density jumps at the mesh element interfaces.

The mesh vertical movement is determined by means of a conservative arbitrary Lagrangian–Eulerian (ALE) formulation. Though conservativity is naturally achieved, tracer consistency is obtained by a suitable construction of the mesh vertical velocity field, which is defined in such a way that it is fully compatible with the tracer and continuity equations at a discrete level.

The vertically adaptive mesh approach is implemented in the three-dimensional version of the geophysical and environmental flow Second-generation Louvain-la-Neuve Ice-ocean Model (SLIM 3D; www.climate.be/slim). Idealised benchmarks, aimed at simulating the oscillations of a sharp thermocline, are dealt with. Then, the relevance of the vertical adaptivity technique is assessed by simulating thermocline oscillations of Lake Tanganyika. The results are compared to measured vertical profiles of temperature, showing similar stratification and outcropping events.

Publications Copernicus
Special issue
Download
Short summary
The discontinuous Galerkin (DG) finite element method is well suited for the modelling of three-dimensional flows exhibiting strong density gradients. Here, a vertical adaptive mesh method is developed for DG finite element methods and implemented into SLIM 3D. This technique increases drastically the accuracy of simulations including strong stratification, without affecting the simulation cost. SLIM 3D is then used to simulate the thermocline oscillations of Lake Tanganyika.
The discontinuous Galerkin (DG) finite element method is well suited for the modelling of...
Citation
Share