The project investigates important aspects of the Arbitrary Lagrangian-Eulerian (ALE) layer motion framework and high-order weighted essentially non-oscillatory (WENO) advection schemes, in order to fully exploit the potential of these new concepts in realistic ocean climate modelling applications. During the first phase, these concepts have been identified as the most promising techniques to significantly reduce spurious mixing in ocean models. However, now the efforts in basic research have to be extended to address emerged challenges related to general robustness and efficiency as well as other further mandatory model adjustments. The main goals are:
- Development of a robust generalized layer motion algorithm based on Lagrangian layer motion and a combination of different regridding strategies.
- Adaptation and optimization of high-order numerical schemes for remapping, internal pressure gradient and WENO advection to the resulting unstructured mesh layout with sloping layers in FESOM.
- Development of new diagnostics for diapycnal mixing and internal pressure gradient errors to assess the energetic consistency of the newly designed model components.
With these efforts, we aim for enabling a new era of energy-consistent climate simulations, which will not be dominated by spurious numerical mixing anymore, but by the advanced and wellcalibrated physically-motivated mixing parameterizations developed in other subprojects of this CRC.
The proposed project aims to further develop, assess and analyse numerical algorithms leading to reduction in spurious diapycnal mixing in ocean circulation models. This goal will be achieved by (i) the design and implementation of vertical mesh motion algorithms that reduce spurious mixing; (ii) use of advective schemes with isopycnal diffusion and special design of limiters; (iii) development and analysis of high-order advection algorithms relying on high-order flux evaluation.
Reduction of spurious mixing by Lagrangian layer motion
- realistic applications
- with different dynamical regimes
- combination of individual
- layer motion techniques
- triggering of regridding
- efficient mesh regularization
- analysis of diapycnal mixing
- interpretation of mean
- (thickness-weighted) quantities
Reduction of spurious mixing by new advection schemes and by stabilization with isoneutral dissipation
Improved understanding of solvers for generalized Riemann problems
Fast and robust solvers available, but only few rigorous analysis
What do approximate solvers
actually compute from an analytical perspective?
What is the common analytical structure of different solvers?
Two new insights, important steps towards closing the gap
Reducing spurious diapycnal mixing in ocean models
My part of work is studying the newly proposed methods of rotating the diffusive part of the advection schemes into the isoneutral plane.
Hi everyone, my name is Margarita and I’m a PhD student of the subproject M5 “Reducing spurious diapycnal mixing in ocean models.” This project is about development and analysis of algorithms leading to reducing spurious mixing in ocean models.
Particularly my part of work is studying the newly proposed methods of rotating the diffusive part of the advection schemes into the isoneutral plane. I work in AWI under supervision of Sergey Danilov.
By the moment adapted to triangular meshes with vertex-based and cell-based placement of scalar variables algorithm described by Lemarie at al. (2012) was implemented on the sea-ice model FESOM2.0. I started doing reference potential energy (RPE) analysis of the implemented harmonic version of the algorithm. Also biharmonic version is to be implemented soon. RPE analysis will be held for both versions and also on dissorted meshes. This analysis allows to determine spurious mixing depending on advection schemes and meshes type.
The future work includes analysis of spurious mixing with variance decay technique (by Knut Klingbeil et al.), analysis of regular and irregular meshes, carrying of realistic ocean simulations and analysis of spurious mixing under real conditions.
Solving equations faster and more accurately
Designing high-order numerical methods needs a good understanding of the mathematical equations we want to solve.
Hi, I’m Claus and I am a PostDoc in Project M5. The focus of my research in this project is the development and analysis of high-order advection methods and high-oder flux evaluation techniques for ocean models. Our goal is to reduce spurious diapycnal mixing in ocean models and I’m working on the intersection between mathematical analysis and numerical methods to help with that.
High-order numerical methods for flow computations are becoming increasingly more popular in computational engineering, but may be not as widespread in climate and ocean science. So what’s the deal with high and low order?
Classical finite differences and finite volume methods for the discretization of partial differential equations use one piece of information in each cell (say, a function value at the cell-center or an integral average over the cell) and put this information into a discrete version of the PDE. This allows for fast and robust algorithms, but in order to resolve small scale features, we often need very fine grids. On the other hand, the high-order methods we are interested in, use higher degree polynomials or other nonlinear functions in each cell. This extra information allow us resolve more features of the solution, even on coarser girds. In many computational fluid dynamics applications this leads to a smaller overall computational time and we would like to show that this is also true for problems in ocean science.
However, carelessly throwing high-order polynomials at your problem is a sure recipe for failure. If we want to include more features of the analytical solution in our numerical solution, we need a good understanding of the analytical properties and how they can be translated into our numerical scheme. For our particular application, e.g., we want tight control over diffusion properties and need to tune our methods to avoid artificial diffusion without losing stability properties that numerical diffusion brings.
In short, designing high-order numerical methods needs a good understanding of the mathematical equations we want to solve. I’m happy to be involved with learning more about the equations in ocean science, so that we can solve them faster and more accurately.
Schulz, K., Klingbeil, K., Morys, C., & Gerkema, T. (2020). The fate of mud nourishment in response to short-term wind forcing. Estuar. Coast, 10.1007/s12237-020-00767-4.
Smolentseva, M., & Danilov, S. (2020). Comparison of several high-order advection schemes for vertex-based triangular discretization. Ocean Dyn., 70(4), 463-479, https://doi.org/10.1007/s10236-019-01337-4 .
Lorenz, M., Klingbeil, K., & Burchard, H. (2020). Numerical study of the exchange flow of the Persian Gulf using an extended Total Exchange Flow analysis framework. J. Geophys. Res.: Oceans, 125(2), e2019JC015527, https://doi.org/10.1029/2019JC015527 .
Scholz, P., Sidorenko, D., Gurses, O., Danilov, S., Koldunov, N., Wang, Q., Sein, D., Smolentseva, M., Rakowsky, N. & Jung, T. (2019). Assessment of the Finite VolumE Sea Ice Ocean Model (FESOM2.0), Part I: Description of selected key model elements and comparison to its predecessor version, Geosci. Model Dev., https://doi.org/10.5194/gmd-2018-329.
Klingbeil, K., J. Becherer, E. Schulz, H. E. de Swart, H. M. Schuttelaars, A. Valle-Levinson and H. Burchard (2019). Thickness-weighted averaging in tidal estuaries and the vertical distribution of the Eulerian residual transport. J. Phys. Oceanogr., doi: https://doi.org/10.1175/JPO-D-18-0083.1.
Stähler, S. C., Panning, M. P., Hadziioannou, C., Lorenz, R. D., Vance, S., Klingbeil, K., & Kedar, S. (2019). Seismic signal from waves on Titan's seas. Earth Planet Sc. Lett., 520, 250-259, doi: https://doi.org/10.1016/j.epsl.2019.05.043.
Lorenz, M., K. Klingbeil, P. MacCready, and H. Burchard (2019). Numerical issues of the Total Exchange Flow (TEF) analysis framework for quantifying estuarine circulation, Ocean Sci., 15, 601-614.
Gräwe, U., K. Klingbeil, J. Kelln, and S. Dangendorf (2019). Decomposing mean sea level rise in a semi-enclosed basin, the Baltic Sea. J. Climate, doi: https://doi.org/10.1175/JCLI-D-18-0174.1.
Iske, A. (2019). Approximation Theory and Algorithms for Data Analysis. Texts App. Math., 68, Springer, doi: 10.1007/978-3-030-05228-7.
Klingbeil, K., Burchard, H., Danilov, S., Goetz, C. & Iske, A. (2019). Reducing spurious diapycnal mixing in ocean models. In Energy Transfers in Atmosphere and Ocean (pp. 245-286). Springer, Cham., doi: https://doi.org/10.1007/978-3-030-05704-6_8.
Rybicki, M., Moldaenke, C., Rinke, K., Dahlhaus, H., Klingbeil, K., Holtermann, P. L. ... & J. Zhu (2019). WP-C: A Step Towards Secured Drinking Water: Development of an Early Warning System for Lakes. In Chinese Water Systems (pp. 159-205). Springer, Cham, doi: https://doi.org/10.1007/978-3-319-97568-9_5.
Burchard, H., Bolding, K., Feistel, R., Gräwe, U., MacCready, P., Klingbeil, K., Mohrholz, V., Umlauf, L., & van der Lee, E. M. , (2018). The Knudsen theorem and the Total Exchange Flow analysis framework applied to the Baltic Sea, Prog. Oceanogr., 165, 268-286, doi: https://doi.org/10.1016/j.pocean.2018.04.004.
Slavik, K., Lemmen, C., Zhang, W., Kerimoglu, O., Klingbeil, K. & Wirtz, K. W. (2018). The large-scale impact of offshore wind farm structures on pelagic primary productivity in the southern North Sea. Hydrobiologia, 1-19, doi: 10.1175/JAS-D-17-0114.1.
Klingbeil, K., Debreu, L., Lemarié, F. & Burchard, H. (2018). The numerics of hydrostatic structured-grid coastal ocean models: state of the art and future perspectives. Ocean Model., Vol. 125, 80-105, doi: https://doi.org/10.1016/j.ocemod.2018.01.007.
Frassl, M., B. Boehrer, P. Holtermann, W. Hu, K. Klingbeil, Z. Peng, ... & K. Rinke (2018). Opportunities and Limits of Using Meteorological Reanalysis Data for Simulating Seasonal to Sub-Daily Water Temperature Dynamics in a Large Shallow Lake. Water-Sui., 10(5), 594, doi: https://doi.org/10.3390/w10050594.
Lemmen, C., Hofmeister, R., Klingbeil, K., Nasermoaddeli, M. H., Kerimoglu, O., Burchard, H., Kösters, F. & Wirtz, K. W. (2018). Modular System for Shelves and Coasts (MOSSCO v1.0) – a flexible and multi-component framework for coupled coastal ocean ecosystem modelling, Geosci. Model Dev., 10.5194/gmd-2017-138 .
Nasermoaddeli, M. H., Lemmen, C., Stigge, G., Kerimoglu, O., Burchard, H., Klingbeil, K., Hofmeister, R., Kreus, M., Wirtz, K. W. & Kösters, F. A (2018). A model study on the large-scale effect of macrofauna on the suspended sediment concentration in a shallow shelf sea Estuarine, Coastal and Shelf Science, Geosci. Model Dev., https://doi.org/10.1016/j.ecss.2017.11.002.