Beginner question: FreeSurfaceSolver in two dimensions

Numerical methods and mathematical models of Elmer
Post Reply
MTrevers
Posts: 10
Joined: 16 Nov 2015, 20:19
Antispam: Yes

Beginner question: FreeSurfaceSolver in two dimensions

Post by MTrevers »

Hi All,

I'm currently trying to get to grips with Elmer/Ice. To help me learn how to put together .sif files, I've been trying to set up the very simple case of a narrow 2D column of ice deforming under the pressure of it's own weight, with no external forcings. I've decided to adapt Thomas' simple Test_Glacier example from the beginner Elmer/Ice training course.

So far I've been unable to get my case to run. I receive multiple instances of the following:

Code: Select all

ERROR:: ElementMetric: Degenerate 1D element: 8041
... followed by:

Code: Select all

ERROR:: LUDecomp: Matrix is singluar.
*** Error in `ElmerSolver': free(): invalid pointer: 0x0000000002831e20 ***
Program received signal SIGABRT: Process abort signal.
As far as I can see, the problem seems to be that I'm trying to use two instances of the FreeSurfaceSolver to change the mesh in two dimensions. This is why I've posted the question here instead of under the Elmer/Ice forum.

I've attached my full example directory. Is anyone able to take a look at my .sif and suggest what I'm doing wrong/how I might be able to put it right?

Thanks
Matt
Attachments
ice_column.tar.gz
(451.9 KiB) Downloaded 256 times
Post Reply