Elmer Discussion Forum

Thanks Peter!

One more question regarding elemental variable. Are they supported for Restart procedure (writing and reading in .result)?

Thanks again
Olivier

Hi, When creating a variable from a solver, one can use the VariableAddVector subroutine: SUBROUTINE VariableAddVector( Variables,Mesh,Solver,Name,DOFs,Values,& Perm,Output,Secondary,VarType,Global,InitValue,IpPoints) By default it will create a nodal value. Which value should be given to VarTyp...

Hello, The problem is that you have zero thickness at the front of the glacier and that you should take care that the mesh has a minimum height by setting Minimum Mesh Height = Real 1.0 Correct Surface = Logical True in the Solver section of StructuredMeshMapper. Moreover, for a 2d flow line problem...

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