Search found 80 matches

by gagliar
25 Sep 2014, 08:37
Forum: Elmer/Ice
Topic: Setting basal pressure til fraction of overburden
Replies: 5
Views: 11977

Re: Setting basal pressure til fraction of overburden

Hi Christian, The Weertman friction law doesn't read the 'External Pressure' keyword. This friction law doesn't dependent on the water pressure. If you want to have the friction modulated by the water pressure you will have to use the Coulomb friction law. Have a look on the wiki: http://elmerice.el...
by gagliar
18 Jul 2014, 18:18
Forum: Elmer/Ice
Topic: 2d Ice Shelf without Grounding Line
Replies: 3
Views: 7728

Re: 2d Ice Shelf without Grounding Line

Hi Reinhard,

I had a look in your sif file. Did you check that you really have a constant flux in your inflow boundary?
If you really have constant BC then i would expect that the steady solution is not a function of the time step...

Regards
Olivier
by gagliar
23 May 2014, 09:05
Forum: Elmer/Ice
Topic: Problems with Elmer/ice tests with newest revision
Replies: 2
Views: 6510

Re: Problems with Elmer/ice tests with newest revision

Hello David,

We are aware of this problem and the tests should be modified very soon on the svn to account for this new feature of Elmer. The good news is that we can see the interest of such tests!

Regards
Olivier
by gagliar
10 Jan 2014, 18:37
Forum: Elmer/Ice
Topic: Small basal velocity value patch from Adjoin inverse method
Replies: 7
Views: 15987

Re: Small basal velocity value patch from Adjoin inverse method

Hi Yongmei,

No value or a value = 0 is completely different. If you have no value for the surface velocity, it should not be imposed as a zero value. You should impose velocity only where you have measured value.

Hope it helps
Olivier
by gagliar
04 Dec 2013, 10:20
Forum: Elmer/Ice
Topic: Grid2DInterpolator updated
Replies: 7
Views: 9966

Re: Grid2DInterpolator updated

Hi Rupert, I'm replying to your two last posts. I found these two lines in the InterpolateDEM function: ix = INT((x-xb0)/dbx)+1 iy = INT((y-yb0)/dbx)+1 Surely the second line should be divide by dby not dbx? This would be a problem for grids with non-square cells with variation in the y-direction. Y...
by gagliar
29 Nov 2013, 20:55
Forum: Elmer/Ice
Topic: Troubles with grounding line simulations
Replies: 4
Views: 7077

Re: Troubles with grounding line simulations

Hi Rupert, I think your troubles came from the fact that you are using the Newton linearization of the Stokes equation. When the contact is tested, the boundary conditions are changed during the non-linear iterations of the Stokes solution (some nodes are released) and therefore if you want to use t...
by gagliar
27 Nov 2013, 15:39
Forum: Elmer/Ice
Topic: Problem with solving temperature after inverse method
Replies: 15
Views: 21852

Re: Problem with solving temperature after inverse method

Hi Denis and Martina, Two comments on the discussion above. The fact that the Mesh Velocities have to be set to zero after a relaxation problem when looking to the steady solution for temperature comes from the fact that the TemperateIce solver is using Mesh Velocities in the ALE formulation. For a ...
by gagliar
27 Nov 2013, 09:33
Forum: Elmer/Ice
Topic: Troubles with grounding line simulations
Replies: 4
Views: 7077

Re: Troubles with grounding line simulations

Hi Rupert,

Can you send me your 2D setup and I can have a look. From which test case did you start?

Regards
Olivier
by gagliar
15 Nov 2013, 12:43
Forum: Elmer/Ice
Topic: Grid2DInterpolator updated
Replies: 7
Views: 9966

Re: Grid2DInterpolator updated

Hi Ruppet,
Would you like to have the noData value and a tolerance available as optional keywords in the sif solver syntax? Let me know if so and I'll add that.
Yes, that would be nice.

regards
Olivier
by gagliar
15 Nov 2013, 12:40
Forum: Bug reports
Topic: Bug corrected in AdvectionReaction solver
Replies: 0
Views: 2555

Bug corrected in AdvectionReaction solver

Hello, I have corrected a bug in the AdvectionReaction (version 6431, 15/11/2013). There was a missing term in the ALE formulation of the advectionreaction solver (should modify results if you have non-zero mesh velocity). See explanations below. The AdvectionReaction solver solves the following equ...