On 06/06/2018 09:06 PM, Ben Melosh wrote:
I am trying to model Stokes flow at very high
viscosities (~10^20), from my results it looks like I cannot achieve this and am better
suited using viscosities of ~10^5 or so. I know there is this older post  which
"A word of warning though: SfePy does not
have a dedicated Navier-Stokes solver (the Newton (+ and direct) solver use
that is in the navier_stokes.py example is rather naive, and works for small
problems with only, with viscosity "high enough")"
Has there been a different solver implemented since this was posted? If not, what do you
recommend to work with such large viscosities?
I have only some general recommendations:
- having a matrix with very different magnitudes of items in different blocks
is problematic - you might need to use some scaling, or choose different
physical units for your input.
- You can try using the FieldSplit preconditioner from PETSc for the Stokes
problem, if it is so large that it cannot be solved using a direct solver
(Umfpack, or, in the git version, also MUMPS). Some of the relevant options of
the solver are mentioned in the docstring of
examples/multi_physics/biot_parallel_interactive.py, see PETSc docs for more, and
examples/multi_physics/biot_short_syntax.py shows how to pass any additional
options to the PETSc solvers from a problem description file.
Does it help?
 "Problem with the mesh while trying to solve the Navier Stokes equations with
SfePy mailing list -- sfepy(a)python.org
To unsubscribe send an email to sfepy-leave(a)python.org