Abstract
Realizations of the numerical solution of the scalar transport equation on the sphere, written in divergent form, are presented. Various temporal discretizations are considered: the one-step Taylor–Galerkin method (TG2), the two-step Taylor–Galerkin method of the second (TTG2), third (TTG3), and fourth (TTG4) orders. The standard Finite-Element Galerkin method with linear basis functions on a triangle is applied as spatial discretization. The flux correction technique (FCT) is implemented. Test runs are carried out with different initial profiles: a function from C∞ (Gaussian profile) and a discontinuous function (slotted cylinder). The profiles are advected by reversible, nondivergent velocity fields, therefore the initial distribution coincides with the final one. The case of a divergent velocity field is also considered to test the conservation and positivity properties of the schemes. It is demonstrated that TG2, TTG3, and TTG4 schemes with FCT applied give the best result for small Courant numbers, and TTG2, TTG4 are preferable in case of large Courant number. However, TTG2+FCT scheme has the worst stability. The use of FCT increases the integral errors, but ensures that the solution is positive with high accuracy. The implemented schemes are included in the dynamic core of a new sea ice model developed using the INMOST package. The acceleration of the parallel program and solution convergence with spatial resolution are demonstrated.
Acknowledgment
The authors are grateful to the staff of the INM RAS: Gordey Goyman, Pavel Perezhogin, Vladimir Shashkin, Vasily Kramarenko, and Alexander Danilov for valuable advice and comments, help with the construction of grids and with mastering the INMOST and Ani3D software packages.
-
Funding: The study was carried out at the INM RAS with the support of the Russian Science Foundation (project 21-71-30023).
References
[1] R. Comblen, S. Legrand, E. Deleersnijder, and V. Legat, A finite element method for solving the shallow water equations on the sphere. Ocean Modelling28 (2009), No. 1-3,12-23.10.1016/j.ocemod.2008.05.004Search in Google Scholar
[2] A. Danilov, Unstructured tetrahedral mesh generation technology. Comput. Math. Math. Phys. 50 (2010), 139-156.10.1134/S0965542510010124Search in Google Scholar
[3] A. Danilov, K. Terekhov, I. Konshin, and Y. Vassilevski, Parallel software platform INMOST: a framework for numerical modeling. Supercomp. Frontiers Innov. 2 (2015), No. 4, 55-66.Search in Google Scholar
[4] S. Danilov, Q. Wang, R. Timmermann, N. Iakovlev, D. Sidorenko, M. Kimmritz, T. Jung, and J. Schroter, Finite-element sea ice model (FESIM), version 2. Geoscientific Model Development 8 (2015), No. 6,1747-1761.10.5194/gmd-8-1747-2015Search in Google Scholar
[5] Error values for Finite Element advection tests on sphere. https://github.com/chuck97/FEtransport_resuitsSearch in Google Scholar
[6] INMOST: a toolkit for distributed mathematical modelling. http://www.inmost.orgSearch in Google Scholar
[7] R. Kwok and D. Sulsky, Arctic Ocean sea ice thickness and kinematics: satellite retrievals and modeling. Oceanography 23 (2010), No. 4.134-143.10.5670/oceanog.2010.11Search in Google Scholar
[8] K. Lipnikov, Yu. Vassilevski, A. Danilovet al., AdvancedNumericallnstruments3D. http://sourceforge.net/projects/ani3dSearch in Google Scholar
[9] R. Löhner, K. Morgan, J. Peraire, and M. Vahdati, Finite element flux-corrected transport (FEM-FCT) for the Euler and Navier–Stokes equations. Int.J. Numer. Meth. Fluids 7 (1987), 1093-1109.10.1002/fld.1650071007Search in Google Scholar
[10] S. Marcq and J. Weiss, Influence of sea ice lead-width distribution on turbulent heat transfer between the ocean and the atmosphere. The Cryosphere 8 (2012), 143-156.10.5194/tc-6-143-2012Search in Google Scholar
[11] C. Mehlmann, S. Danilov, M. Losch, J. F. Lemieux, N. Hutter, T. Richter, P. Blain, E. C. Hunke, and P. Korn, Simulating linear kinematic features in viscous-plastic sea ice models on quadrilateral and triangular grids. arXiv:2103.04431v. 7 Mar 2021.10.1029/2021MS002523Search in Google Scholar
[12] R. D. Nair and P. H. Lauritzen, A class of deformational flow test cases for linear transport problems on the sphere. J. Comput. Phys. 229 (2010), 8868-8887.10.1016/j.jcp.2010.08.014Search in Google Scholar
[13] S. Petrovand N. Iakovlev, The optimized Finite Element dynamical core of the Arctic Ocean sea ice model. Commun. Computer Inform. Sci., 2021 (in press).10.1007/978-3-030-78759-2_32Search in Google Scholar
[14] L. Quartapelle, Numerical Solution of the Incompressible Navier–Stokes Equations, Vol. 113. Birkhäuser, 1993.10.1007/978-3-0348-8579-9Search in Google Scholar
[15] Y. Vassilevski, K. Terekhov, K. Nikitin, and I. Kapyrin, Parallel Finite Volume Computation on General Meshes. Springer, Cham, 2020.10.1007/978-3-030-47232-0Search in Google Scholar
[16] O. C. Zienkiewicz and R. L. Taylor, The Finite Element Method, Fifth Edition. Vol. 3: Fluid Dynamics. Oxford, Butterworth-Heinemann, 2000.Search in Google Scholar
© 2021 Walter de Gruyter GmbH, Berlin/Boston
