[petsc-users] FieldSplit and Biot's poroelasticity

[Karin&NiKo]

Dear Petsc-gurus,

I am solving Biot's poroelasticity problem :
 [image: Images intégrées 1]

I am using a mixed P2-P1 finite element discretization.

I am using the fieldsplit framework to solve the linear systems. Here are
the options I am using :
-pc_type fieldsplit
-pc_field_split_type schur
-fieldsplit_0_pc_type gamg
-fieldsplit_0_pc_gamg_threshold -1.0
-fieldsplit_0_ksp_type gmres
-fieldsplit_0_ksp_monitor
-fieldsplit_1_pc_type sor
-fieldsplit_1_ksp_type gmres
-pc_fieldsplit_schur_factorization_type upper


By increasing the mesh size, I get increasing numbers of outer iterations.

According to your own experience, among all the features of fieldsplit, was
is the "best" set of preconditioners for this rather classical problem  in
order to get an extensible solver (I would like to solve this problem on
some tens millions of unknowns of some hundreds of procs)?

Thanks,
Nicolas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161213/608b755f/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 9086 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20161213/608b755f/attachment.png>