conf for damping methods

conf/coupled/gamm-paper/coarse-biventricle.conf

@@ -9,7 +9,7 @@ h_min=1e-6
 
 
 activeStress = 0
-Mesh = TestBiventricleInstantActivation
+Mesh = TestBiventricle
 
 ProblemDimension = 3
 ProblemGeometry = Tet
@@ -25,20 +25,20 @@ RAPressure = 0
 LAPressure = 0
 
 
-#activeDeformation = Quarteroni2
+activeDeformation = Quarteroni2
 #g = 0.8 #gamma for updateStretch
 #scale_gamma = 0.8 #0.9, 1.0
 
 withPressure = 1
 
-MechDeltaTime  = 0.001 #0.004# 0.001
+MechDeltaTime  = 0.004 #0.004# 0.001
 ElphyDeltaTime = 0.001 # 0.001 #0.000001
 
 #StretchFactorSheet = 2.0
 #StretchFactorNormal = 5.0
 
 linearImplicitSolverCell = 1
-cell_values = 0
+cell_values = 1
 deformConcentrationStretch = 1 #1
 timeStepExtrapolation = 0
 DampExtrapolation = 0.5;
@@ -52,8 +52,10 @@ CoupledProblem = BiventricleCoarse
 ElphyModelClass = MElphyModel
 ElphyModelName = TenTusscher
 TensionModelName = Rossi
-RossiViscosity = 7000
-RossiForce = 7
+RossiViscosity = 500
+#RossiViscosity = 7000
+RossiForce = 10
+#RossiForce = 7
 
 ElphyLevel = 0 #1
 MechLevel = 0 #1
@@ -70,9 +72,6 @@ EndTime = 0.8
 #EndTime = 0.8
 #EndTime = 0.02
 
-MechDeltaTime  = 0.001 #0.004# 0.001
-ElphyDeltaTime = 0.0001 # 0.001 #0.000001
-
 DeltaTime = 0.001
 DeltaTimeMin = 0.00001
 DeltaTimeMax = 0.001
@@ -115,9 +114,6 @@ PressureDepth = 0
 #MechMeshPart = with_vessels
 MechMeshPart = vessel_base
 
-#TractionK = 0
-#TractionC = 0
-
 TractionK = 0#0.1
 TractionC = 0#5e-3
 
@@ -131,7 +127,7 @@ StretchFactorNormal = 4.0
 ######################
 # Material Parameters#
 ######################
-Density = 0#1082e-9
+Density = 1082e-9
 #------------------------
 # Volumetric Penalty
 # Incompressible should be true or false
@@ -141,7 +137,7 @@ VolumetricSplit = 0
 Incompressible = false
 QuasiCompressiblePenalty = Ciarlet #None, Ciarlet, Mixed
 VolumetricPenalty = 1000 # ?
-PermeabilityPenalty = 1000 # ?
+PermeabilityPenalty = 0 #1000 # ?
 
 
 ElphyModel = SemiImplicit
@@ -213,11 +209,11 @@ MechReduction = 1e-10 #1e-12
 MechSteps = 2000
 
 # === Newton Method === #
-NewtonEpsilon = 1e-8 #1e-8
+NewtonEpsilon = 1e-4 #1e-8
 NewtonReduction = 1e-12 #1e-5
 NewtonSteps = 300
 NewtonLineSearchSteps = 10
-NewtonLinearizationReduction = 1e-4 #1e-12
+NewtonLinearizationReduction = 1e-3 #1e-12
 
 #Smoother = GaussSeidel
 #Smoother = PointBlockGaussSeidel
@@ -225,14 +221,14 @@ NewtonLinearizationReduction = 1e-4 #1e-12
 presmoothing = 3
 postsmoothing = 3
 
-ConfigVerbose = -1
+ConfigVerbose = 1
 AssembleVerbose = 1
 
 # = Solver Verbose =
 MultigridVerbose = 0
 BaseSolverVerbose = -1
 LinearVerbose = -1
-NewtonVerbose = -1
+NewtonVerbose = 1
 EulerVerbose = 0
 LinearImplicitVerbose = 1
 ElphyLinearSolverVerbose =1
@@ -255,5 +251,5 @@ PressureSolverVTK = 0
 DynamicSolverVTK = 10
 ElphySolverVTK = 10
 CoupledSolverVTK = 10
-PlotCalcium=1
+PlotCalcium=0
 PlotVTK = 1
\ No newline at end of file