Glacier width in FreeSurfaceSolver

Extension of Elmer in computational glaciology

Glacier width in FreeSurfaceSolver

Postby vanGeffen » 04 Jun 2015, 12:08

Hi,

I have a question: is there an easy way to account for the glacier width W(x) in the FreeSurfaceSolver when running a flow-line model? I want to determine the free surface following the equation for the ice thickness as depicted in the attachment.

Thanks,

Sharon

ChangeIceThickness.png
ChangeIceThickness.png (24.96 KiB) Viewed 2133 times

with H the ice thickness, u the depth-averaged horizontal velocity, W the glacier width and b the surface mass balance
vanGeffen
 
Posts: 1
Joined: 03 Apr 2015, 13:16

Re: Glacier width in FreeSurfaceSolver

Postby joeatodd » 09 Jun 2015, 15:58

Hi Sharon,

It's possible to do this as an accumulation term in the freesurface equation, but I think that it's more accurate to directly account for the addition of mass in the stokes equations.

I wrote some code to do this in a flowline case by adding a "Pressure Load" following Peter's advice in this thread:

viewtopic.php?f=3&t=2273

The plan is for this code to make its way into the elmer/ice repository before too long. If you like I can pass it on to you for the time being, but it comes without any warranty ;)

Joe
joeatodd
 
Posts: 34
Joined: 02 Feb 2012, 18:49


Return to Elmer/Ice

Who is online

Users browsing this forum: No registered users and 1 guest