diff --git a/src/foam/foam-2D.cfg b/src/foam/foam-2D.cfg index 31ca8ebd25648cef2598819b3e384cabd5a2e613..d5028e7fbe6091dc556ef5585b4c71f7fc906153 100644 --- a/src/foam/foam-2D.cfg +++ b/src/foam/foam-2D.cfg @@ -6,7 +6,7 @@ smallestDiameter = 6.25e-2 # 2e-3 [m] smallestDiameter = 6.25e-2 # 2e-3 [m] [timeSteps] -refinementTolerance = 1e-3 # 1e-5 +refinementTolerance = 1e-5 # 1e-5 [u0.solver] tolerance = 1e-8 @@ -18,7 +18,7 @@ tolerance = 1e-8 tolerance = 1e-8 [v.fpi] -tolerance = 1e-3 +tolerance = 1e-10 [solver.tnnmg.preconditioner.basesolver] tolerance = 1e-10 diff --git a/src/foam/foam.cfg b/src/foam/foam.cfg index dc15cc635841b432f2268e4c9fff60540970fcbf..85411b67803e6e1a51c12f2a3fc37a9f5261eec7 100644 --- a/src/foam/foam.cfg +++ b/src/foam/foam.cfg @@ -1,5 +1,5 @@ # -*- mode:conf -*- -outPath = test # output written to ./output/outPath +outPath = no-friction # output written to ./output/outPath gravity = 9.81 # [m/s^2] @@ -10,12 +10,12 @@ bulkModulus = 4.12e7 # [Pa] #2190 poissonRatio = 0.3 # [1] #0.11 [body0.elastic] density = 5e3 # [kg/m^3] #750 -shearViscosity = 0 # [Pas] -bulkViscosity = 0 # [Pas] +shearViscosity = 0.0 # [Pas] +bulkViscosity = 0.0 # [Pas] [body0.viscoelastic] density = 5e3 # [kg/m^3] #750 -shearViscosity = 0 # [Pas] -bulkViscosity = 0 # [Pas] +shearViscosity = 0.0 # [Pas] +bulkViscosity = 0.0 # [Pas] [body1] length = 5.00 # [m] @@ -24,12 +24,12 @@ bulkModulus = 4.12e7 # [Pa] poissonRatio = 0.3 # [1] [body1.elastic] density = 5e3 # [kg/m^3] -shearViscosity = 0 # [Pas] -bulkViscosity = 0 # [Pas] +shearViscosity = 0.0 # [Pas] +bulkViscosity = 0.0 # [Pas] [body1.viscoelastic] density = 5e3 # [kg/m^3] -shearViscosity = 0 # [Pas] -bulkViscosity = 0 # [Pas] +shearViscosity = 0.0 # [Pas] +bulkViscosity = 0.0 # [Pas] [boundary.friction] @@ -37,9 +37,9 @@ C = 0.0 # [Pa] mu0 = 0.6 # [ ] V0 = 1e-6 # [m/s] L = 1e-5 # [m] -initialAlpha = -10.0 # [ ] +initialAlpha = 0.0 #-10.0 # [ ] stateModel = AgeingLaw -frictionModel = Truncated #Regularised +frictionModel = None #Truncated #Regularised [boundary.friction.weakening] a = 0.010 # [ ] b = 0.015 # [ ] @@ -52,28 +52,28 @@ b = 0.015 # [ ] sigmaN = 0.0 # [Pa] [boundary.dirichlet] -finalVelocity = 2e-4 # [m/s] +finalVelocity = 1e-1 #2e-4 # [m/s] [io] data.write = true printProgress = true restarts.first = 0 -restarts.spacing= 20 +restarts.spacing= 50 restarts.write = true #true vtk.write = true [problem] -finalTime = 15 # [s] #1000 +finalTime = 10 # [s] #1000 bodyCount = 2 [initialTime] timeStep = 0 relativeTime = 0.0 -relativeTau = 1e-4 # 1e-6 +relativeTau = 0.1 #1e-4 # 1e-6 [timeSteps] scheme = newmark -timeSteps = 1e4 # 1e4 +timeSteps = 10 # 1e4 [u0.solver] maximumIterations = 100