This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
problems:parstokes [2017/11/16 13:03] tzwinger [Introduction] |
problems:parstokes [2019/04/13 13:12] (current) tzwinger [Introduction] |
||
---|---|---|---|
Line 1: | Line 1: | ||
====== ParStokes: Iterative Solver for Large-Scale Glaciological Applications ====== | ====== ParStokes: Iterative Solver for Large-Scale Glaciological Applications ====== | ||
+ | |||
+ | ==== Alternatives ==== | ||
+ | Since March 2019, a new [[problems: | ||
==== Introduction ==== | ==== Introduction ==== | ||
Line 46: | Line 49: | ||
A sif-file for the ParStokes module should have (at least) three Solver sections, one | A sif-file for the ParStokes module should have (at least) three Solver sections, one | ||
- | for the coupled problem, one for the velocity system with A, and one for the pressure | + | for the coupled problem, one for the velocity system with '' |
- | system with M. An example could look like this: | + | system with '' |
Solver 1 | Solver 1 | ||
Equation = " | Equation = " | ||
+ | Exec Solver = " | ||
Procedure = " | Procedure = " | ||
Variable = -dofs 3 " | Variable = -dofs 3 " | ||
Line 57: | Line 61: | ||
Solver 2 | Solver 2 | ||
Equation = " | Equation = " | ||
+ | Exec Solver = " | ||
Procedure = " | Procedure = " | ||
Variable = -dofs 1 " | Variable = -dofs 1 " | ||
Line 103: | Line 108: | ||
In the case of sliding conditions, the keyword "Slip Coefficient" | In the case of sliding conditions, the keyword "Slip Coefficient" | ||
| | ||
- | (the same as the keyword | + | (the same as the keyword |
| | ||
- | (same as "External Pressure" | + | (same as '' |
== Remarks == | == Remarks == | ||
Line 117: | Line 122: | ||
Element = "p:1 b: | Element = "p:1 b: | ||
| | ||
- | NB.: The setting | + | The setting |
Bubbles in Global System = False | Bubbles in Global System = False | ||
is not needed, as it is the default in any of the solvers. | is not needed, as it is the default in any of the solvers. | ||
Line 130: | Line 135: | ||
Linear System Row Equilibration = Logical True | Linear System Row Equilibration = Logical True | ||
- | d) For the systems with A and M it is possible to preserve symmetry by setting | + | d) For the systems with '' |
- | Linear System Scaling = Logical True | + | Linear System Scaling = Logical True |
- | Linear System Row Equilibration = Logical False | + | Linear System Row Equilibration = Logical False |
In that case CG can be used instead of BiCGStab. | In that case CG can be used instead of BiCGStab. | ||
Line 166: | Line 171: | ||
outer iteration decreases. | outer iteration decreases. | ||
- | f) ParStokes by default does only solve the Stokes (and not Navier-Stokes) equation. If one wishes to include acceleration terms, setting the keyword " | + | f) ParStokes by default does only solve the Stokes (and not Navier-Stokes) equation. If one wishes to include acceleration terms, setting the keyword " |
+ | |||
+ | Bubbles in Global System = True | ||
+ | |||
+ | to avoid the by default implemented elimination of bubbles from the global system. In that case one also should take care that the '' | ||
g) It should be possible to restart a ParStokes simulation | g) It should be possible to restart a ParStokes simulation |