The procedure can run before,but recently,I run the procedure in a new machine, it doesn't work,errors!
Who knows？please！Error !too many iterations were needed ! in the navier-Stokes equations!

Solver 4
Equation = "Navier-Stokes"
Stabilization Method = String Stabilized
Flow Model = Stokes
Exported Variable 1 = -dofs 1 "dSdt"
Exported Variable 2 = -dofs 1 "dS"
Exported Variable 3 = -dofs 1 "ZbDEM"
Exported Variable 4 = -dofs 1 "ZsDEM"
Linear System Solver = Iterative
Linear System Iterative Method = BICGStab
Linear System Preconditioning = ILU0
Linear System Convergence Tolerance = Real 1.0e-4 !8
Linear System Max Iterations = 500
Linear System Residual Output = 500
Linear System Abort Not Converged = True
Nonlinear System Max Iterations = 50
Nonlinear System Convergence Tolerance = 1.0e-3 !5
Nonlinear System Newton After Iterations = 5
Nonlinear System Newton After Tolerance = 1.0e-02
Nonlinear System Relaxation Factor = 1.00
Steady State Convergence Tolerance = Real 1.0e-3
End

Try with better preconditioners, in fact you could use the IncompressibleNSVec module that comes with block preconditioner and it custom made for the glaciology community. Black-box cheap precondtioners do not work for all problems.

-Peter

PS. Please don't sent the same issue to multiple forums.

Try with better preconditioners, in fact you could use the IncompressibleNSVec module that comes with block preconditioner and it custom made for the glaciology community. Black-box cheap precondtioners do not work for all problems.

-Peter

PS. Please don't sent the same issue to multiple forums.

Thank you!but i can't find the use of the IncompressibleNSVec module,how to use it,where have a case,please!I am a newer!

hello, i have used the block preconditioning strategy in the IncompressibleNSVec,but i met a new error in the TemperateIceSolver,error in the picture,so should i put the block preconditioning strategy in the TemperateIceSolver and even the freesurfacesolver?

Solver 8
!Exec Solver = "Never"
Equation = String "Homologous Temperature Equation"
Procedure = File "ElmerIceSolvers" "TemperateIceSolver"
! Comment next line in parallel, as EliminateDirichlet does
! not work in parallel!
!------------------------------------------------------------
! Before Linsolve = "EliminateDirichlet" "EliminateDirichlet"
Variable = String "Temperature"
Stabilize = True
Optimize Bandwidth = Logical True
Linear System Solver = "Iterative"
Linear System Iterative Method = "GCR"
Linear System Max Iterations = 500
Linear System Convergence Tolerance = 1.0E-07
Linear System Abort Not Converged = False
Linear System Preconditioning = "ILU1"
Linear System Residual Output = 250
Nonlinear System Convergence Tolerance = 1.0E-06
Nonlinear System Max Iterations = 50
!Nonlinear System Relaxation Factor = Real 9.999E-01
Steady State Convergence Tolerance = 1.0E-03
! the contact algorithm (aka Dirichlet algorithm)
!-----------------------------------------------------
!Apply Dirichlet = Logical True
! those two variables are needed in order to store
! the relative or homologous temperature as well
! as the residual
!-------------------------------------------------
Exported Variable 1 = String "Temperature Homologous"
Exported Variable 1 DOFs = 1
Exported Variable 2 = String "Temperature Residual"
Exported Variable 2 DOFs = 1
Apply Dirichlet = True
End