A strength implicit correction scheme for the viscous-plastic sea ice model

Abstract With increasing awareness of the role of the polar regions in global climate, efficient and accurate simulation of sea ice is an important issue. The numerical solution procedure for the viscous-plastic sea ice model is examined. Recent developments are drawn upon and new discretisation practices introduced that improve efficiency and the accuracy of the simulation. The main components of the new solution procedure are bounded and conservative finite volume discretisation and a reinterpretation of the model structure, including the introduction of a method to resolve momentum-ice strength coupling. It is found that traditional segregated momentum solution procedures do not ensure mass conservation, which may introduce substantial thickness errors in Arctic sea ice simulations. The strength implicit correction scheme ensures mass conservation for a similar computational expense as previous, non-conservative methods.

[1]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[2]  M. Leppäranta,et al.  The role of plastic ice interaction in marginal ice zone dynamics , 1985 .

[3]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[4]  A. D. Gosman,et al.  Automatic Resolution Control for the Finite-Volume Method, Part 2: Adaptive Mesh Refinement and Coarsening , 2000 .

[5]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[6]  Hrvoje Jasak,et al.  Application of the finite volume method and unstructured meshes to linear elasticity , 2000 .

[7]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[8]  Marika M. Holland,et al.  Simulating the ice‐thickness distribution in a coupled climate model , 2001 .

[9]  M. Perić,et al.  A collocated finite volume method for predicting flows at all speeds , 1993 .

[10]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[11]  Hrvoje Jasak,et al.  A tensorial approach to computational continuum mechanics using object-oriented techniques , 1998 .

[12]  W. Hibler.  A Dynamic Thermodynamic Sea Ice Model , 1979 .

[13]  William D. Hibler,et al.  On an efficient numerical method for modeling sea ice dynamics , 1997 .

[14]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[15]  Peter Lemke,et al.  Results of the Sea Ice Model Intercomparison Project: Evaluation of sea ice rheology schemes for use in climate simulations , 2000 .

[16]  I. Demirdžić,et al.  Finite volume method for thermo-elasto-plastic stress analysis , 1993 .

[17]  Josef M. Oberhuber,et al.  Simulation of the atlantic circulation with a coupled sea ice-mixed layer-isopycnal general circulat , 1993 .

[18]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[19]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[20]  R. Colony,et al.  The thickness distribution of sea ice , 1975 .

[21]  I. Demirdzic,et al.  Numerical method for coupled fluid flow, heat transfer and stress analysis using unstructured moving meshes with cells of arbitrary topology , 1995 .

[22]  E. Hunke,et al.  An Elastic–Viscous–Plastic Model for Sea Ice Dynamics , 1996 .