conf/elasticity/benchmarks/LandBenchmarkBiVentricle.conf

@@ -8,20 +8,14 @@ h0=1e-4
 h_min=1e-6
 #logfile = CG_DG_sig1000#AS_T4_1_deg1_gamma07
 
-multiDiscs  = 0
-# CAUTION: if useDifferentMeshes = 0, it will use the entry in Mesh!!!
-useDifferentMeshes = 0 # possibility to change meshes isn't implemented right now
-MechDiscretizationList =  Conforming,DG  # chosing different discs for homotopy
-MechPolynomialDegreeList = 1 #if there is only one entry, it usess the same degree for all discs
-MeshList = TestBiventricle #if there is only one entry, it usess the same mesh for all discs
-# If |DGPenaltyList| != |MechDiscretizationList| DGPenalty = 0 für fehlende Einträge
-DGPenaltyList= 0,10000 #10,100,1000,2000, 10000 (geht auch ohne scale nicht)
-
 activeStress = 0 #35
 Model = ActiveStrainElasticity
 MechRuntype = Default
 MechDynamics = Static
-MechDiscretization = DG#Conforming
+MechDiscretization = Conforming
+
+UseExtrapolation=true
+dampingExtrapolation = 0.5
 
 DGSign=-1.0
 DGPenalty= 1000# Penalty for EG=200

src/elasticity/assemble/LagrangeElasticity.hpp

@@ -144,7 +144,7 @@ public:
   void Scale (double s) const {
     //TODO warum konvergiert der LandBenchmark mit mehreren PressureSteps nicht, wenn einkommentiert?
     Material &cellMat = eProblem.GetMaterial();
-// cellMat.ScalePenalty(s*s);
+   cellMat.ScalePenalty(s*s);
   }
       
 };

src/elasticity/solvers/IterativePressureSolver.cpp

@@ -25,7 +25,9 @@ bool IterativePressureSolver::Method(IElasticity &assemble, Vector &u) {
   Vector deltaU(u);
 
   bool useExtrapolation = false;
+  double dampingExtrapolation = 1;
   Config::Get("UseExtrapolation", useExtrapolation);
+  Config::Get("dampingExtrapolation", dampingExtrapolation);
 
   Initialize(assemble, uNew);
   while (!iteration.IsFinished()) {
@@ -55,8 +57,8 @@ bool IterativePressureSolver::Method(IElasticity &assemble, Vector &u) {
     }
     // Extrapolation ab n = 2
     if (useExtrapolation && n >= 2) {
-      deltaU = uNew - uOld;  // Δu = u^n - u^{n-1}
-      uNew += deltaU;        // Extrapolation: u^n = 2 u^n - u^{n-1}
+      deltaU = uNew - uOld;
+      uNew += dampingExtrapolation * deltaU;
     }
   }
   mout.EndBlock();