Abstract:Abstract. Time dependent simulations of ice sheets require two equations to be solved: the mass transport equation, derived from the conservation of mass, and the stress balance equation, derived from the conservation of momentum. The mass transport equation controls the advection of ice from the interior of the ice sheet towards its periphery, thereby changing its geometry. Because it is based on a hyperbolic partial differential equation, a stabilization scheme needs to be employed when solved using the finite element method. Several stabilization schemes exist in the finite element method framework, but their respective accuracy and robustness have not yet been systematically assessed for glaciological applications. Here, we compare classical schemes used in the context of the finite element method: (i) Artificial Diffusion, (ii) Streamline Upwinding, (iii) Streamline Upwind Petrov-Galerkin, (iv) Discontinuous Galerkin, and (v) Flux Corrected Transport. We also look at the stress balance equation, which is responsible for computing the ice velocity that `advects' the ice dowstream. To improve the velocity computation accuracy, the ice sheet modeling community employs several sub-element parameterizations of physical processes at the grounding line, the point where the grounded ice starts to float onto the ocean. Here, we introduce a new sub-element parameterization for the driving stress, the force that drives the ice sheet flow. We analyze the response of each stabilization scheme by running transient simulations forced by ice shelf basal melt. The simulations are based on an idealized ice sheet geometry for which there is no influence of bedrock topography. We also perform transient simulations of the Amundsen Sea Sector, West Antarctica, where real bedrock and surface elevations are employed. In both idealized and real ice sheet experiments, stabilization schemes based on artificial diffusion lead systematically to a bias towards more mass loss in comparison to the other schemes, and therefore, should be avoided or employed with a sufficiently high mesh resolution in the vicinity of the grounding line. We also run diagnostic simulations to assess the accuracy of the driving stress parameterization, which in combination with an adequate parameterization for basal stress, provides improved numerical convergence in ice speed computations and more accurate results.

Laboratory Tests to Investigate the Initial Phase of Ice Ridge Consolidation

The Determination of Heat Transfer Coefficient on Water-ice Surface in a Free Convection

Method for Testing Shear and Tensile Strengths of Freshwater/Seawater Ice

Investigation of Heat and Mass Transfer Characteristics in the Ice Rink: Ice Making, Maintaining and Resurfacing Processes

How the Growth of Ice Depends on the Fluid Dynamics Underneath

Investigation of the temperature field of the ice sheet in a prefabricated curling ice rink: Experiment and finite element simulation

Turbulent heat exchange between water and ice at an evolving ice-water interface

Effect of compact fin structure on heat transfer performance of ice storage unit during freezing

Laboratory Studies on the Parametrization Scheme of the Melting Rate of Ice–Air and Ice–Water Interfaces

An Experimental Investigation on the Unsteady Heat Transfer Process over an Ice Accreting NACA 0012 Airfoil

Equilibrium states of the ice-water front in a differentially heated rectangular cell

Analysis of Solidification Heat Transfer Process on Double-Convection Condition

Identification and Optimization for Effective Thermal Diffusivity of Fresh Ice

Experimental investigation and scale analysis on melting of salty ice in a 3D-printed cavity filled with porous media

Changes of Water/Ice Morphological, Thermodynamic, and Mechanical Parameters During the Freezing Process

STUDY ON THERMAL ICE STRESSES OF A POND IN COLD REGION

Numerical research on the solidification heat transfer characteristics of ice thermal storage device based on a compact multichannel flat tube-closed rectangular fin heat exchanger

HEAT TRANSFER ANALYSIS OF PLANE FREEZE PHASE TRANSITION ON DOUBLE-CONVECTION CONDITION

Experimental Study on Current Velocity Profile under Ice Floe

Experimental Study on the Melting Characteristics Modulation of Cubic Ice Cubes with Different Trace Air Contents under Natural Convection Condition

Assessment of numerical schemes for transient, finite-element ice flow models using ISSM v4.18