i..

Z

Closes #216 (closed)

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

-activeDeformation = 1
+#VolumetricPenalty = 10000
+#PrestressSteps = 20
+#CoupledProblem = BiventricleCoarse
+
+CoupledProblem = DeformedBiventricleCoarse
+ActiveMaterial = Bonet
+PassiveMaterial = Bonet
+GuccioneMat_C = 0.29
+GuccioneMat_bf = 42.0
+GuccioneMat_bs = 16.8 
+GuccioneMat_bfs = 29.4			
+
+# = Removes all files in the data folder upon start ======
+ClearData = true
+# ========================================================
+# = Location of .geo files ===============================
+GeoPath = ../geo/
+# ========================================================
+h0=1e-4
+h_min=1e-6
+
+
+activeStress = 0
+Mesh = TestBiventricle
+
+ProblemDimension = 3
+ProblemGeometry = Tet
+
 WithPrestress = true
-PrestressSteps = 1
-withPressure = 1
+PrestressSteps = 10
+PressureIterations = 20
+PressureDepth = 0
+
+LVPressure = 1.06658
+RVPressure = 0.53329
+RAPressure = 0
+LAPressure = 0
 
-pressureED_LV = 0.001
-pressureED_RV = 0.001
-pressureES_RV = 0.001
-pressureES_RV = 0.001
 
-WithPrestress = false
+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
 ElphyDeltaTime = 0.001 # 0.001 #0.000001
@@ -19,24 +52,22 @@ ElphyDeltaTime = 0.001 # 0.001 #0.000001
 linearImplicitSolverCell = 1
 cell_values = 1
 deformConcentrationStretch = 1 #1
-withPressure = 0
 timeStepExtrapolation = 0
 DampExtrapolation = 0.5;
-Mesh = KovachevaBiventricle
-#Mesh = KovachevaBiventricleInstantActivation
-ClearData = true
 
 Model = Coupled
 CoupledModel = Segregated
 
-CoupledProblem = DeformedBiventricleCoarse
-#CoupledProblem = BiventricleCoarse
+#CoupledProblem = DeformedBiventricleCoarse
+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
@@ -53,9 +84,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
@@ -73,11 +101,18 @@ ClearDistribution = 0
 MechDiscretization = Conforming
 MechRuntype = Default
 #MechDynamics = Newmark
+
+## Rayleigh
+
 MechDynamics = RayleighNewmark
 RayleighAlpha = 100
+
+##  KelvinVoigt
+
 #MechDynamics = ViscoNewmark
 #ViscoelasticDamping = KelvinVoigt
 #KelvinVoigtEta = 2
+
 #------------------------
 # Newmark Parameters
 #------------------------
@@ -88,27 +123,16 @@ RayleighAlpha = 100
 NewmarkBeta = 0.25
 NewmarkGamma = 0.5
 
-
-
-# Amount of pressure steps per iteration
-PressureSteps = 0
 # If Iteration with Pressuresteps fails, recursively restart PressureSolver with Steps*2 until iteration limit is reached
 PressureIterations = 10
 PressureDepth = 0
 
-WithPrestress=false
-#WithPrestress=true
-PrestressSteps=10
-
 ######################
 #       Meshes       #
 ######################
-MechMeshPart = with_vessels
+#MechMeshPart = with_vessels
 MechMeshPart = vessel_base
 
-#TractionK = 0
-#TractionC = 0
-
 TractionK = 0#0.1
 TractionC = 0#5e-3
 
@@ -132,7 +156,7 @@ VolumetricSplit = 0
 Incompressible = false
 QuasiCompressiblePenalty = Ciarlet #None, Ciarlet, Mixed
 VolumetricPenalty = 1000 # ?
-PermeabilityPenalty = 1000 # ?
+PermeabilityPenalty = 0 #1000 # ?
 
 
 ElphyModel = SemiImplicit
@@ -204,9 +228,9 @@ MechReduction = 1e-10 #1e-12
 MechSteps = 20000
 
 # === Newton Method === #
-NewtonEpsilon = 1e-8 #1e-8
+NewtonEpsilon = 1e-5 #1e-8
 NewtonReduction = 1e-12 #1e-5
-NewtonSteps = 300
+NewtonSteps = 100
 NewtonLineSearchSteps = 10
 NewtonLinearizationReduction = 1e-4 #1e-12
 
@@ -216,7 +240,7 @@ NewtonLinearizationReduction = 1e-4 #1e-12
 presmoothing = 3
 postsmoothing = 3
 
-ConfigVerbose = 10
+ConfigVerbose = 1
 AssembleVerbose = 1
 
 # = Solver Verbose =
@@ -246,5 +270,5 @@ PressureSolverVTK = 0
 DynamicSolverVTK = 10
 ElphySolverVTK = 10
 CoupledSolverVTK = 10
-PlotCalcium=1
+PlotCalcium=0
 PlotVTK = 1
\ No newline at end of file