Sylvain Bouillon

A new modeling framework for sea-ice mechanics based on elasto-brittle rheology

Abstract We present a new modeling framework for sea-ice mechanics based on elasto-brittle (EB) behavior. the EB framework considers sea ice as a continuous elastic plate encountering progressive damage, simulating the opening of cracks and leads. As a result of long-range elastic interactions, the stress relaxation following a damage event can induce an avalanche of damage. Damage propagates in narrow linear features, resulting in a very heterogeneous strain field. Idealized simulations of the Arctic sea-ice cover are analyzed in terms of ice strain rates and contrasted to observations and simulations performed with the classical viscous–plastic (VP) rheology. the statistical and scaling properties of ice strain rates are used as the evaluation metric. We show that EB simulations give a good representation of the shear faulting mechanism that accommodates most sea-ice deformation. the distributions of strain rates and the scaling laws of ice deformation are well captured by the EB framework, which is not the case for VP simulations. These results suggest that the properties of ice deformation emerge from elasto-brittle ice-mechanical behavior and motivate the implementation of the EB framework in a global sea-ice model.

Parallel implementation of a Lagrangian-based model on an adaptive mesh in C++: Application to sea-ice

Abstract We present a parallel implementation framework for a new dynamic/thermodynamic sea-ice model, called neXtSIM , based on the Elasto–Brittle rheology and using an adaptive mesh. The spatial discretisation of the model is done using the finite-element method. The temporal discretisation is semi-implicit and the advection is achieved using either a pure Lagrangian scheme or an Arbitrary Lagrangian Eulerian scheme (ALE). The parallel implementation presented here focuses on the distributed-memory approach using the message-passing library MPI. The efficiency and the scalability of the parallel algorithms are illustrated by the numerical experiments performed using up to 500 processor cores of a cluster computing system. The performance obtained by the proposed parallel implementation of the neXtSIM code is shown being sufficient to perform simulations for state-of-the-art sea ice forecasting and geophysical process studies over geographical domain of several millions squared kilometers like the Arctic region.