| Summary: | While there has been a colossal effort in the ongoing decades, the ability to simulate ocean ice has fallen behind various parts of the climate system and most Earth System Models are unable to capture the observed adversities of Arctic sea ice, which is, as it were, attributed to our frailty to determine sea ice dynamics. Viscous Plastic rheology is the most by and large recognized model for sea ice dynamics and it is expressed as a set of partial differential equations that are hard to tackle numerically. Using the 1D sea ice momentum equation as a prototype, we use the method of lines based on Euler's backward method. This results in a nonlinear PDE in space only. At that point, we apply the Damped Newton’s method which has been introduced in Looper and Rapetti et al. and used and generalized to 2D in Saumier et al. to solve the Monge-Ampere equation. However, in our case, we need to solve 2nd order linear equation with discontinuous coefficients during Newton iteration. To overcome this difficulty, we use the Finite element method to solve the linear PDE at each Newton iteration. In this paper, we show that with the adequate smoothing and re-scaling of the linear equation, convergence can be guaranteed and the numerical solution indeed converges efficiently to the continuum solution unlike other numerical approaches that typically solve an alternate set of equations and avoid the difficulty of the Newton method for a large nonlinear algebraic system. The finite element solver failed to converge when the original setting of the smoothed SIME with a smoothing constant $K=2.8 \times 10^8$ was used. A much smaller constant of K=100 was necessary. The large smoothing constant K leads to an ill conditioned mass matrix. Graduate
|
|---|