Elmer Discussion Forum

Dear Kevin, the "side" boundaries in the test example for the Reykjavik course simply set the horizontal velocity component (what ever direction this is in your model, in ours it is X or direction #1) to zero while leaving the vertical (in our model Y or #2) without a Dirichlet condition, ...

Checking the code, it seems that the dimension of the problem is being taken care of - i.e., I found no reason why it shouldn't work in 2D. Question is - do you really need it in a 2D problem or are you equally happy with the standard Navier-Stokes Solver and a direct solver?

Regards,

Thomas

Hi Rupert,
your issue might be connected to the fact that the default branch is "devel" and not "elmerice".
A might be at place.

Regards,

Thomas

Hello, I didn't realize your MeshUpdate in the first instance - sorry. Then you can do it such that you just initialize, imposing the data directly to Top Surface and then couple the FS to the MeshUpdate variable, sacrificing the advantage of having a vertically structured mesh. Or, you do as I sugg...

Initial Condition 2 Zs = Variable Coordinate 2 Real cubic include "icescarp70degground.dat" End RefZs = Equals Coordinate 2 End should do it. So no "Equals" for the Zs-initialization. And then, of course, you have to take care that the solver that updates the mesh (StructuredMes...

The main issue is conected to the fact that your "Top Surface" keyword points to the initial data all the time, but rather should point to the result of your free-surface evolution equation. So, the good strategy would be to let the initial shape define the initial condition for the free-s...

Hello,
seems you want to solve the Grounding Line Problem, which for 2D is explained (including link to example) under http://elmerice.elmerfem.org/wiki/doku. ... undingline

Good luck!

Best Regards,

Thomas

Hi, perhaps the name "External Pressure" is misleading. It is basically the normal component of the Cauchy stress vector, which is the Cauchy stress tensor dot product with surface normal. It follows from the weak formulation of the Stokes equation. It is a Neumann condition to the system....

Hi Christian, I guess you are using this 0.6*overburden as the effective normal pressure in some sliding law, right? Did you check that your values are dimensional consistent with the rest of your material parameters? Simply, if you do your computations in MPa-a-m system (which usually is the case i...

So your case went through - something then went wrong outside the run. Reading your error message from before, my suspicion is, that perhaps your permissions for the test subdirectory are such that you are not allowed to write anything into it. If you do a chmod -R u+rwX contact does it help? Best R...