Search found 41 matches

by fgillet
02 Jul 2020, 16:23
Forum: Elmer/Ice
Topic: getFrictionLoads User Function resulting in segmentation fault error
Replies: 8
Views: 1458

Re: getFrictionLoads User Function resulting in segmentation fault error

Hi Thomas, Yes I agree that core issue is with the GetEquation in the USF and this should be changed. However, this does not crash for me even if I don't have a body for the bed; and I think this is why it has been un-noticed so far, espacially when the Aruco test case was designed. So to play safe ...
by fgillet
02 Jul 2020, 15:09
Forum: Elmer/Ice
Topic: getFrictionLoads User Function resulting in segmentation fault error
Replies: 8
Views: 1458

Re: getFrictionLoads User Function resulting in segmentation fault error

Hello,
code is ok in getFrictionLoads if Equation is not assigned, which should be the case if element do not belong to a body with a given equation?

I think this could be fixed in GetEquationId (DefUtils.F90) which should return with Found=.False. if body_id < 1 ?

Cheers
Fab
by fgillet
04 Jun 2020, 10:53
Forum: Elmer/Ice
Topic: Spiky Ice Sheet Surface in transient flow line model
Replies: 6
Views: 1384

Re: Spiky Ice Sheet Surface in transient flow line model

Hello; I think it is also sensitive to the stabilisation method. By default the Free surafce solver uses the Bubbles; SUPG stabilisation is also possible; However in this case the stabilisation parameter is optimal in steady-state. There is a lot of litterature on this subject and possible stabilisa...
by fgillet
10 Apr 2018, 11:04
Forum: Elmer/Ice
Topic: GL dynamics and Stokes in Elmer v8.3
Replies: 2
Views: 2642

GL dynamics and Stokes in Elmer v8.3

Dear ElmerIce users, Following important changes on how Dirichlet Conditions are handle in Elmer v8.3, we found that the conditional dirichlet condition for the normal velocity in the grounded part of the ice sheet was not properly released during grounding line retreat. This has now be fixed in the...
by fgillet
10 Apr 2018, 10:48
Forum: Elmer/Ice
Topic: Elmer/Ice tests
Replies: 2
Views: 3024

Re: Elmer/Ice tests

Hello Rupert, thanks for the report and sorry for the late reply. At Grenoble we are running all the Elmer/Ice tests, after each code revision. I try to keep an eye on the results as much as I can. I have updated most of the tests with failure. On the 52 elmerice tests (ctest -L elmerice); There is ...
by fgillet
04 Apr 2018, 16:12
Forum: Bug reports
Topic: Pb with conditional Dirichlet
Replies: 4
Views: 2464

Re: Pb with conditional Dirichlet

ok.
Thanks!

Fabien
by fgillet
04 Apr 2018, 13:00
Forum: Bug reports
Topic: Pb with conditional Dirichlet
Replies: 4
Views: 2464

Re: Pb with conditional Dirichlet

Hi Peter; Yes a call to DefaultInitialize() will also solve our problem. Currently we have this problem with the Navier-Stokes solver, I see that it is calling InitializeToZero(). Replacing the InitializeToZero() line by CALL DefaultInitialize() will sove our problem. I see that there is few other s...
by fgillet
03 Apr 2018, 19:10
Forum: Bug reports
Topic: Pb with conditional Dirichlet
Replies: 4
Views: 2464

Pb with conditional Dirichlet

Hello; It seems that Solver % Matrix % ConstrainedDOF is never reinitailised to FALSE so that when a node get a Dirichlet condition, it will always apply a Dirichlet condition. This could be a problem when using a "Condition" for the dirichlet condition and when the condition change from TRUE to FAL...
by fgillet
02 Feb 2018, 14:00
Forum: Bug reports
Topic: Passive condition
Replies: 3
Views: 2290

Re: Passive condition

Hello Peter; Here is two tests for the passive condition keyword. I have been a bit confused as I get the expected behaviour for solvers with DOFs > 1 (e.g. Stokes); this is because the values for each component are set in a second round where it does not visit the section that reset the values to t...
by fgillet
25 Jan 2018, 14:12
Forum: Bug reports
Topic: Passive condition
Replies: 3
Views: 2290

Passive condition

Hello, Dirichlet conditions for passive elements (set with the keywords VarName Condition Passive = Logical in Body Forces) are not taken into account anymore due to changes in how Dirichlet conditions are set. Lines 4700 to 4728 of SolverUtils.F90 take care of the matrix entries for the passive ele...