Search found 94 matches

by tzwinger
12 Nov 2014, 08:55
Forum: Elmer/Ice
Topic: Help to simulate Sheet/Shelf problems
Replies: 4
Views: 6445

Re: Help to simulate Sheet/Shelf problems

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
by tzwinger
26 Sep 2014, 11:22
Forum: Elmer/Ice
Topic: Setting basal pressure til fraction of overburden
Replies: 5
Views: 8819

Re: Setting basal pressure til fraction of overburden

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....
by tzwinger
23 Sep 2014, 12:15
Forum: Elmer/Ice
Topic: Setting basal pressure til fraction of overburden
Replies: 5
Views: 8819

Re: Setting basal pressure til fraction of overburden

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...
by tzwinger
05 Sep 2014, 12:30
Forum: Elmer/Ice
Topic: Test failed for test case 3 (Contact) in Elmer/Ice
Replies: 5
Views: 7766

Re: Test failed for test case 3 (Contact) in Elmer/Ice

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...
by tzwinger
04 Sep 2014, 22:50
Forum: Elmer/Ice
Topic: Test failed for test case 3 (Contact) in Elmer/Ice
Replies: 5
Views: 7766

Re: Test failed for test case 3 (Contact) in Elmer/Ice

Hi, unfortunately, I am not able to reproduce your error. Could you try to go into the subdirectory of this test and run it manually and send the output. To that end, you have to declare the following variables (provided you are running bash) export ELMER_GRID="ElmerGrid" export FORT="...
by tzwinger
07 Aug 2014, 09:54
Forum: Elmer/Ice
Topic: non-reproducible machine-specific convergence failure
Replies: 3
Views: 6048

Re: non-reproducible machine-specific convergence failure

Hi Rupert,
I send you the setup scripts of MUMPS on taito - although they might not be of particular help, as system differs. If I interpret them right, we use parmetis there.

Cheers,

Thomas
by tzwinger
07 Aug 2014, 08:17
Forum: Elmer/Ice
Topic: non-reproducible machine-specific convergence failure
Replies: 3
Views: 6048

Re: non-reproducible machine-specific convergence failure

Hi Rupert, Checking the commit history of the SVN repository, the most likely changes after revision 6798 that could influence your simulation where actually committed by you (Grid2DInterpolator.f90 and FlowDepth,f90). I don't know if those influence something set for your Stokes solver. Else, my la...
by tzwinger
03 May 2014, 11:50
Forum: Elmer/Ice
Topic: Mask to restrict sliding in ice free areas
Replies: 10
Views: 16936

Re: Mask to restrict sliding in ice free areas

Hi Christian, I apparently missed a few potential issues in your setup: !! bedrock: Boundary Condition 3 Name = "bedrock" Height = Real 0.0 Bottom Surface = Equals bedrockDEM Normal-Tangential Velocity = Logical True Flow Force BC = Logical True Flow Solution 1 = Real 0.0 Flow Solution 1 C...
by tzwinger
03 May 2014, 10:51
Forum: Elmer/Ice
Topic: Mask to restrict sliding in ice free areas
Replies: 10
Views: 16936

Re: Mask to restrict sliding in ice free areas

Christian, I didn't realise that you use Block Pre-Conditioner (BPC). As you apparently were able to run with no-slip condition (Dirichlet) the issue is not directly connected to your Mesh or Solver rather than the boundary condition. Solver settings for ParStokes are well documented here: http://el...
by tzwinger
29 Apr 2014, 17:58
Forum: Elmer/Ice
Topic: Mask to restrict sliding in ice free areas
Replies: 10
Views: 16936

Re: Mask to restrict sliding in ice free areas

Hi Christian, I made an approach with a mountain glacier using Weertman law. I didn't switch between a Dirichlet condition and the sliding law, but simply, depending on the depth, switched from one to another sliding coefficient. Something like that: $ function depthwise(D) {\ if (D > 10.0) {\ _dept...