Updating side boundaries with two free surfaces
Posted: 26 Apr 2018, 21:49
Hello All,
I am attempting to model a 2D section (x and z, across flow) of a floating ice shelf.
In my model set up, I define both the ice shelf surface (hs) and ice shelf base (hb) as free surfaces, whose position in z changes with accumulation and basal melt, respectively. I keep the x-location of the right and left hand sides fixed, as I want my 2D slice of the ice shelf to not be able to grow outwards to the sides. However, I do want the right and left boundaries to be able to change their position in z to match the ice shelf surface and ice shelf base.
I am not able to achieve this at present by updating the mesh along the RHS and LHS (but I could be doing the update incorrectly). The extent of the LHS and RHS stays fixed in z- while the free surface ice shelf base, for example, moves down in z, resulting in a discontinuity at the corner which is compensated for across a few corner elements.
I have attached the .sif file (BCs at the bottom) and .grd files I am using below, as well as a figure showing the issues at the corners:
How can I update the left and right boundaries in z such that they span the distance from the ice shelf surface to the ice shelf base?
Thanks for your help and time,
Laura
I am attempting to model a 2D section (x and z, across flow) of a floating ice shelf.
In my model set up, I define both the ice shelf surface (hs) and ice shelf base (hb) as free surfaces, whose position in z changes with accumulation and basal melt, respectively. I keep the x-location of the right and left hand sides fixed, as I want my 2D slice of the ice shelf to not be able to grow outwards to the sides. However, I do want the right and left boundaries to be able to change their position in z to match the ice shelf surface and ice shelf base.
I am not able to achieve this at present by updating the mesh along the RHS and LHS (but I could be doing the update incorrectly). The extent of the LHS and RHS stays fixed in z- while the free surface ice shelf base, for example, moves down in z, resulting in a discontinuity at the corner which is compensated for across a few corner elements.
I have attached the .sif file (BCs at the bottom) and .grd files I am using below, as well as a figure showing the issues at the corners:
How can I update the left and right boundaries in z such that they span the distance from the ice shelf surface to the ice shelf base?
Thanks for your help and time,
Laura