scholarly journals Simple algorithms to compute meridional overturning and barotropic streamfunction on unstructured meshes

2020 ◽  
Author(s):  
Dmitry Sidorenko ◽  
Sergey Danilov ◽  
Nikolay Koldunov ◽  
Patrick Scholz ◽  
Qiang Wang

Abstract. Computation of barotropic and meridional overturning streamfunctions for models formulated on unstructured meshes is commonly preceded by interpolation to a regular mesh. This operation destroys the original conservation which can be then artificially imposed to make the computation possible. An elementary method is proposed that avoids interpolation and preserves conservation in a strict model sense. The method is described as applied to the discretization of the Finite volumE Sea ice -- Ocean Model (FESOM2) on triangular meshes. It however is generalizable to collocated vertex based discretization on triangular meshes and to both triangular and hexagonal C-grid discretizations.

2020 ◽  
Vol 13 (7) ◽  
pp. 3337-3345
Author(s):  
Dmitry Sidorenko ◽  
Sergey Danilov ◽  
Nikolay Koldunov ◽  
Patrick Scholz ◽  
Qiang Wang

Abstract. Computation of barotropic and meridional overturning streamfunctions for models formulated on unstructured meshes is commonly preceded by interpolation to a regular mesh. This operation destroys the original conservation, which can be then artificially imposed to make the computation possible. An elementary method is proposed that avoids interpolation and preserves conservation in a strict model sense. The method is described as applied to the discretization of the Finite volumE Sea ice – Ocean Model (FESOM2) on triangular meshes. It, however, is generalizable to colocated vertex-based discretization on triangular meshes and to both triangular and hexagonal C-grid discretizations.


2015 ◽  
Vol 8 (2) ◽  
pp. 855-896 ◽  
Author(s):  
S. Danilov ◽  
Q. Wang ◽  
R. Timmermann ◽  
N. Iakovlev ◽  
D. Sidorenko ◽  
...  

Abstract. The Finite-Element Sea-Ice Model, used as a component of the Finite-Element Sea ice Ocean Model, is presented. Version 2 includes the elastic-viscous-plastic (EVP) and viscous-plastic (VP) solvers and employs a flux corrected transport algorithm to advect the ice and snow mean thicknesses and concentration. The EVP part also includes a modified approach proposed recently by Bouillon et al., which is characterized by an improved stability compared to the standard EVP approach. The model is formulated on unstructured triangular meshes. It assumes a collocated placement of ice velocities, mean thicknesses and concentration at mesh vertices, and relies on a piecewise-linear (P1) continuous elements. Simple tests for the modified EVP and VP solvers are presented to show that they may produce very close results provided the number of iterations is sufficiently high.


2021 ◽  
Author(s):  
Ozgur Gurses ◽  
Judith Hauck ◽  
Moritz Zeising ◽  
Laurent Oziel

<p>Marine biogeochemistry models are generally coupled to a physical ocean model. The biases in these coupled models can be attributed to simplified and empirical representation of biogeochemical processes, insufficient spatial mesh resolution which has an impact on the transport and mixing of biogeochemical substances in the ocean, and a deficit of physical parameterizations that intent to mimic unresolved processes such as eddies. Ocean Biogeochemical models based on variable mesh resolution proved to be convenient tools due to their computational efficiency and flexibility. Unlike standard structured-mesh ocean models, the mesh flexibility allows for a realistic representation of eddy dynamics in certain regions. Here, we present preliminary results of the coupling between the Finite-volumE Sea ice-Ocean Model (FESOM2.0) and the biogeochemical model REcoM2 (Regulated Ecosystem Model 2) in a coarse spatial resolution global configuration.<br>Surface maps of the simulated nutrients, chlorophyll a and net primary production (NPP) are comparable to available observational data sets. The control simulation forced with the JRA55-do data set reveals a realistic spatial distribution of nutrients, nanophytoplankton and diatom NPP, carbon stocks and fluxes. <br>FESOM2 utilizes a new dynamical core based on a finite-volume approach. The computational efficiency is about 2-3 times higher than the previous version FESOM1.4, whereas the quality of the simulated ocean and sea ice conditions and representation of biogeochemical variables are comparable in the two models. Thus, the new coupled model FESOM2- REcoM2 is very promising for ocean biogeochemical modelling applications.</p>


2020 ◽  
Author(s):  
Dmitry Sidorenko ◽  
Sergey Danilov ◽  
Nikolay Koldunov ◽  
Patrick Scholz

<p>The Atlantic meridional overturning circulation (AMOC) is the most common diagnostics of numerical simulations. Generally it is computed as a streamfunction of zonally averaged flow along the constant depth. More rarely it is computed as zonally averaged along constant isopycnals. The latter computation, however, allows one to better distinguish between water masses and physical processes contributing to the meridional transport. We analyze the AMOC in global simulations based on the Finite-volumE Sea ice–Ocean Model (FESOM 2.0) using eddy permitting to eddy resolving configurations in the North Atlantic. We (1) split the AMOC computed in density space into the constitutes induced by surface buoyancy fluxes and cross isopycnal transformations, (2) identify the water masses which contribute to the formation of the North Atlantic Deep Water and (3) study the AMOC response to the permitting or resolving eddies in the North Atlantic ocean.</p>


2019 ◽  
Author(s):  
Nikolay V. Koldunov ◽  
Vadym Aizinger ◽  
Natalja Rakowsky ◽  
Patrick Scholz ◽  
Dmitry Sidorenko ◽  
...  

Abstract. A study of the scalability of the Finite-volumE Sea ice-Ocean circulation Model, Version 2.0 (FESOM2), the first mature global model of its kind formulated on unstructured meshes, is presented. This study includes an analysis of main computational kernels with a special focus on bottlenecks in parallel scalability. Several model enhancements, improving this scalability for large numbers of processes, are described and tested. Model grids at different resolutions are used on four HPC systems with differing computation and communication hardware to demonstrate model's scalability and throughput. Furthermore, strategies for improvements in parallel performance are presented and assessed. We show that in terms of throughput FESOM2.0 is on par with the state-of-the-art structured ocean models and in realistic eddy resolving configuration (1/10° resolution) can produce about 16 years per day on 14 000 cores. This suggests that unstructured-mesh models are becoming extremely competitive tools in high-resolution climate modelling. It is shown that main bottlenecks of FESOM parallel scalability are the two-dimensional components of the model, namely the computations of external (barotropic) mode and the sea-ice model. It is argued that these bottlenecks are shared with other general ocean circulation models.


2020 ◽  
Vol 12 (12) ◽  
Author(s):  
Dmitry Sidorenko ◽  
Sergey Danilov ◽  
Vera Fofonova ◽  
William Cabos ◽  
Nikolay Koldunov ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document