Exemple #1
0
overland_flux.PhaseSources.water.Geom.domain.Value = 0.0

#-----------------------------------------------------------------------------
# Exact solution specification for error calculations
#-----------------------------------------------------------------------------

overland_flux.KnownSolution = 'NoKnownSolution'


#-----------------------------------------------------------------------------
# Set solver parameters
#-----------------------------------------------------------------------------

overland_flux.Solver = 'Richards'
overland_flux.Solver.MaxIter = 2500
overland_flux.OverlandFlowDiffusive = 0


overland_flux.Solver.Nonlinear.MaxIter = 20
overland_flux.Solver.Nonlinear.ResidualTol = 1e-9
overland_flux.Solver.Nonlinear.EtaChoice = 'EtaConstant'
overland_flux.Solver.Nonlinear.EtaValue = 0.01
overland_flux.Solver.Nonlinear.UseJacobian = False
overland_flux.Solver.Nonlinear.UseJacobian = True
overland_flux.Solver.Nonlinear.DerivativeEpsilon = 1e-8
overland_flux.Solver.Nonlinear.StepTol = 1e-20
overland_flux.Solver.Nonlinear.Globalization = 'LineSearch'
overland_flux.Solver.Linear.KrylovDimension = 20
overland_flux.Solver.Linear.MaxRestart = 2

overland_flux.Solver.Linear.Preconditioner = 'PFMGOctree'
sandtank.PhaseSources.water.GeomNames = 'domain'
sandtank.PhaseSources.water.Geom.domain.Value = 0.0

# -----------------------------------------------------------------------------
# Exact solution specification for error calculations
# -----------------------------------------------------------------------------

sandtank.KnownSolution = 'NoKnownSolution'

# -----------------------------------------------------------------------------
# Set solver parameters
# -----------------------------------------------------------------------------

sandtank.Solver = 'Richards'
sandtank.Solver.MaxIter = 2500000
sandtank.OverlandFlowDiffusive = 0

sandtank.Solver.Nonlinear.MaxIter = 100
sandtank.Solver.Nonlinear.ResidualTol = 1e-5
sandtank.Solver.Nonlinear.EtaValue = 0.01
sandtank.Solver.Nonlinear.UseJacobian = False
sandtank.Solver.Nonlinear.UseJacobian = True
sandtank.Solver.Nonlinear.DerivativeEpsilon = 1e-8
sandtank.Solver.Nonlinear.StepTol = 1e-20
sandtank.Solver.Nonlinear.Globalization = 'LineSearch'
sandtank.Solver.Linear.KrylovDimension = 100
sandtank.Solver.Linear.MaxRestarts = 5

sandtank.Solver.Linear.Preconditioner = 'PFMG'
sandtank.Solver.PrintSubsurf = False
sandtank.Solver.Drop = 1E-20