Numerical instabilities (spikes in the ice surface) at grounding line
Numerical instabilities (spikes in the ice surface) at grounding line
I'm running a simulation with a transient grounding line and am observing numerical instabilities in the form of spikes on the upper ice surface and inverted spikes on the lower ice surface. Has anyone else experienced this and has a solution to fix these numerical instabilities?
Re: Numerical instabilities (spikes in the ice surface) at grounding line
Helo,
Can you be more specific? Are you running a Stokes or a SSA simulation? 2d or 3d? From which setup did you started?
Regards,
Olivier
Can you be more specific? Are you running a Stokes or a SSA simulation? 2d or 3d? From which setup did you started?
Regards,
Olivier
Re: Numerical instabilities (spikes in the ice surface) at grounding line
Hi Olivier,
Thanks for your reply. I am running a full Stokes 3D simulation and the grounding line in question is in a collision zone. The basal friction coefficient is set to 7.624e7 and the basal melt rate is dependent on the ice thickness and distance from the GL (Favier et al., 2016). Our current thinking is that the numerical instabilities could be caused by the resulting high melt rates in the collision zone. I have attached a couple of photos for reference.
Cheers,
Clara
Thanks for your reply. I am running a full Stokes 3D simulation and the grounding line in question is in a collision zone. The basal friction coefficient is set to 7.624e7 and the basal melt rate is dependent on the ice thickness and distance from the GL (Favier et al., 2016). Our current thinking is that the numerical instabilities could be caused by the resulting high melt rates in the collision zone. I have attached a couple of photos for reference.
Cheers,
Clara
Re: Numerical instabilities (spikes in the ice surface) at grounding line
Hi Clara,
Did you try to decrease your time step?
Also, having a water pressure dependent friction law might help by removing the discontinuity at the GL.
Hope it helps
Olivier
Did you try to decrease your time step?
Also, having a water pressure dependent friction law might help by removing the discontinuity at the GL.
Hope it helps
Olivier
Re: Numerical instabilities (spikes in the ice surface) at grounding line
Thank you for your help, Olivier. I was able to fix the problem in the end by avoiding a BMB discontinuity at the grounding line.
All the best,
Clara
All the best,
Clara