Elmer Discussion Forum

Dear Gleb,

The form of the friction law here is linear, i.e. \tau_nt = C u_t where C is the "Porous Slip Coef t", for t = 1 and 2. If you want to use more complex friction laws, you should have a look on the wiki in User Functions (Weertman, Coulomb, Budd).

Hope it helps,
Olivier

Hi Clara,

Did you try to decrease your time step?

Also, having a water pressure dependent friction law might help by removing the discontinuity at the GL.

Hope it helps
Olivier

Helo,

Can you be more specific? Are you running a Stokes or a SSA simulation? 2d or 3d? From which setup did you started?

Regards,
Olivier

Hello, When executing a solver on the boundary of the bulk mesh (e.g. the free surface solver), the Dirichlet boundary conditions imposed in the Body Force section associated to this second body are not imposed. The problem is in the SolverUtils routine where the loop is only done over NumberOfBulkE...

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...