Elmer Discussion Forum

Hello, In section "Equation 1", there is no solver indicated! You should have a line of that type in this section Active Solvers (1)= 1 But also you should have at least one solver section which tells which equation you are solving! See examples in /elmerfem/fem/tests of the source code. R...

Hi Cyrille and Martin, What would be great is that you test the fixes proposed by Brandon Berg and Jeremy Bassis on your setups to see if it changes the dependence to time step of your simulations? My understanding of this fix is that you simply need to reintroduce the acceleration term in the force...

Hi all, You may have seen this paper in TCD ( https://www.the-cryosphere-discuss.net/tc-2019-315/ ) by Brandon Berg and Jeremy Bassis regarding the "time step dependence (and fixes) in Stokes simulations of calving ice shelves". As they are solving the contact problem in the same way as us...

Hi Norbert, By convention, the boundary conditions in your sif file should write: Boundary Condition 1 Target Boundaries = 1 ... End Boundary Condition 2 Target Boundaries = 2 ... End Boundary Condition 3 Target Boundaries = 3 ... End Boundary Condition 4 Target Boundaries = 4 ... End ! cavity roof ...

Hello,

I can already see one problem:

Variable 1 lx = Real 6512.347200
Variable 1 ly = Real 2022.023200

are not the correct lengths (lx = 6.5013280e+03 and ly = 2.0110040e+03 in the DEM files).

Does it improve results to have the correct values?

Regards,
Olivier

Dear Nordert,

Can you post your sif and 2d footprint mesh so that I have a look? Did you try to visualize the resulting 3d mesh in paraview instead of ElmerGUI?

Regards,
Olivier

Hi Arran,

I only see your message now. Did you solve your problem or is it still not working?

Regards,
Olivier

It is not possible to switch off the cavity (just the channels). I guess you will not be able to do what you want directly with the GlaDS solver but you will have to write your own solver with the appropriate h=h(N) function implemented.
Regards,
Olivier

Yes, this is m_b from Eq. (12). So in fact, contrary to what I was saying in my previous post, m_b can be directly evaluated as a function of the geothermal heat flux and the frictional heating.

Dear Arran, GlaDS doesn't see the geothermal heat flux, but only melt, trough the source term: Hydraulic Potential Volume Source. If you want to include melt from geothermal heat flux (and basal friction and internal deformation) you will have to compute the temperature field and evaluate how much w...