[petsc-users] Staggered solver phase field

[PERRIER-MICHON Augustin]

Dear Petsc users,

I am currently dealing with finite element fracture analysis using phase 
field model. To perform such simulations, I have to develop a staggered 
solver : mechanical problem is solved at constant damage and damage 
problem is solved at constant displacement.

I created 2 TS solver and 2 DMPLEX for each "physics".
Each physics's system is built using TSSetIFunction and TSSetIJacobian 
with associated functions.

The TS calls are performed with TSSTEP in order to respect staggered 
solver scheme in iterative loops.

My question : Is the using of TSSTEP function adapted to a staggered 
solver ? How to use this function in my framework ? Have you got any 
other suggestions or advices ?

Thanks a lot
Best regards

-- 
Augustin PERRIER-MICHON
PhD student institut PPRIME
Physics and Mechanics of materials department
ISAE-ENSMA
Téléport 2
1 Avenue Clément ADER
86361 Chasseneuil du Poitou- Futuroscope
Tel : +33-(0)-5-49-49-80-97