foam.cfg

# -*- mode:conf -*-
[general]
outPath = no-friction # output written to ./output/outPath
gravity         = 9.81     # [m/s^2]

[body0]
length          = 6.0      # [m]
height          = 1.0     # [m]
bulkModulus     = 0.0 #4.12e9   # [Pa] #2190
poissonRatio    = 0.0 #0.3     # [1]  #0.11
[body0.elastic]
density         = 5e3      # [kg/m^3] #750
shearViscosity  = 0.0     # [Pas]
bulkViscosity   = 0.0     # [Pas]
[body0.viscoelastic]
density         = 5e3      # [kg/m^3] #750
shearViscosity  = 0.0     # [Pas]
bulkViscosity   = 0.0     # [Pas]

[body1]
length          = 5.00     # [m]
height          = 1.00     # [m]
bulkModulus     = 4.12e7    # [Pa]
poissonRatio    = 0.3     # [1]
[body1.elastic]
density         = 5e3      # [kg/m^3]
shearViscosity  = 0.0     # [Pas]
bulkViscosity   = 0.0     # [Pas]
[body1.viscoelastic]
density         = 5e3      # [kg/m^3]
shearViscosity  = 0.0     # [Pas]
bulkViscosity   = 0.0     # [Pas]


[boundary.friction]
C               = 0.0       # [Pa]
mu0             = 0.6      # [ ]
V0              = 1e-6     # [m/s]
L               = 1e-5  # [m]
initialAlpha    = -10.0  # [ ]
stateModel      = AgeingLaw
frictionModel   = None #Tresca #None #Truncated #Regularised
[boundary.friction.weakening]
a               = 0.010    # [ ]
b               = 0.015    # [ ]
[boundary.friction.strengthening]
a               = 0.010    # [ ]
b               = 0.015    # [ ]


[boundary.neumann]
sigmaN          = 0.0      # [Pa]

[boundary.dirichlet]
finalVelocity   = 2e-3 #2e-4     # [m/s]

[io]
data.write      = true
printProgress   = true
restarts.first  = 0
restarts.spacing= 50
restarts.write  = true #true
vtk.write       = true

[problem]
finalTime       = 15     # [s] #1000
bodyCount       = 2

[initialTime]
timeStep = 0
relativeTime = 0.0
relativeTau = 5e-4 # 1e-6

[timeSteps]
scheme = newmark
timeSteps = 2 #2e3

[u0.solver]
maximumIterations = 100
verbosity         = full

[a0.solver]
maximumIterations = 100
verbosity         = full

[v.solver]
maximumIterations = 100
verbosity         = quiet

[v.fpi]
maximumIterations = 10000
lambda            = 0.5

[solver.tnnmg.preconditioner]
mode         = additive
patchDepth   = 1
maximumIterations = 2
verbosity         = quiet
[solver.tnnmg.preconditioner.patchsolver]
maximumIterations = 100
verbosity         = quiet
[solver.tnnmg.preconditioner.basesolver]
maximumIterations = 10000
verbosity         = quiet

[solver.tnnmg.main]
pre   = 1
multi = 5 # number of multigrid steps
post  = 0