| Download | - View final version: Smoothed particle hydrodynamics implementation of the standard viscous–plastic sea-ice model and validation in simple idealized experiments (PDF, 6.8 MiB)
|
|---|
| DOI | Resolve DOI: https://doi.org/10.5194/tc-18-1013-2024 |
|---|
| Author | Search for: Marquis, OresteORCID identifier: https://orcid.org/0000-0003-4771-3801; Search for: Tremblay, Bruno; Search for: Lemieux, Jean-FrançoisORCID identifier: https://orcid.org/0000-0003-2084-5759; Search for: Islam, Mohammed1ORCID identifier: https://orcid.org/0000-0002-2129-5333 |
|---|
| Affiliation | - National Research Council of Canada. Ocean, Coastal and River Engineering
|
|---|
| Funder | Search for: Natural Sciences and Engineering Research Council of Canada; Search for: National Research Council Canada |
|---|
| Format | Text, Article |
|---|
| Abstract | The viscous–plastic (VP) rheology with an elliptical yield curve and normal flow rule is implemented in a Lagrangian modelling framework using the smoothed particle hydrodynamics (SPH) meshfree method. Results show, from a perturbation analysis of SPH sea-ice dynamic equations, that the classical SPH particle density formulation expressed as a function of sea-ice concentration and mean ice thickness leads to incorrect plastic wave speed. We propose a new formulation for particle density that gives a plastic wave speed in line with theory. In all cases, the plastic wave in the SPH framework is dispersive and depends on the smoothing length (i.e., the spatial resolution) and on the SPH kernel employed in contrast to its finite-difference method (FDM) implementation counterpart. The steady-state solution for the simple 1D ridging experiment is in agreement with the analytical solution within an error of 1 %. SPH is also able to simulate a stable upstream ice arch in an idealized domain representing the Nares Strait in a low-wind regime (5.3 m s⁻¹) with an ellipse aspect ratio of 2, an average thickness of 1 m and free-slip boundary conditions in opposition to the FDM implementation that requires higher shear strength to simulate it. In higher-wind regimes (7.5 m s⁻¹) no stable ice arches are simulated – unless the thickness is increased – and the ice arch formation showed no dependence on the size of particles, in contrast to what is observed in the discrete-element framework. Finally, the SPH framework is explicit, can take full advantage of parallel processing capabilities and shows potential for pan-Arctic climate simulations. |
|---|
| Publication date | 2024-03-04 |
|---|
| Publisher | Copernicus Publications
European Geosciences Union |
|---|
| Licence | |
|---|
| In | |
|---|
| Related data | |
|---|
| Language | English |
|---|
| Peer reviewed | Yes |
|---|
| Export citation | Export as RIS |
|---|
| Report a correction | Report a correction (opens in a new tab) |
|---|
| Record identifier | 8823d658-e6aa-4a33-a5cb-0363194d769f |
|---|
| Record created | 2025-04-16 |
|---|
| Record modified | 2025-11-03 |
|---|
