[petsc-users] Staggered solver phase field

[PERRIER-MICHON Augustin]

Dear Mr Bourdin,

thank you for your answer and the remarks.

I will performed time dependent multi-physics analysis including crack 
propagation afterward. To anticipate this time dependency, I chose to 
use TS solver instead of SNES or TAO. Plus, I thought that TS solver can 
be used for quasi-static problems as well.

In my previous simulations with a monolithic TS solver, I controlled the 
time step during all the calculation. In my opinion I could do the same 
in this framework and not let TS solvers adapt the step time. A 
synchronization of the two solvers is necessary.

With these informations, is this framework and especially TSSTEP 
function compatible with my problem ?

Thanks a lot
Augustin

Le 2025-04-23 16:58, Blaise Bourdin a écrit :
> Augustin,
> 
>  Out of curiosity, why TS and not SNES? At the very least the damage
> problem should be a constrained minimization problem so that you can
> model criticality with respect to the phase-field variable.
> Secondly, I would be very wary about letting TS adapt the time step by
> itself. In quasi-static phase-field fracture, the time step affects
> the crack path, not the order of the approximation in time. I doubt
> that any of the mechanisms in TS are appropriate here.
> 
> You are welcome to dig into my implementation for inspiration, or
> reuse it for your problems https://urldefense.us/v3/__https://github.com/bourdin/mef90__;!!G_uCfscf7eWS!YNLdSyA8YkLHqmi8DisU8Hz_g4lAFvUm-N7RZHxO_u1CkWa6ZEMYoUBG2so6IcRb57XYSVCvOLtnD-Fx0ic5BOuyGj5kvwHxvw$ 
> 
> Blaise
> 
>> On Apr 23, 2025, at 10:20 AM, PERRIER-MICHON Augustin
>> <augustin.perrier-michon at ensma.fr> wrote:
>> 
>> [You don't often get email from augustin.perrier-michon at ensma.fr.
>> Learn why this is important at
>> https://urldefense.us/v3/__https://aka.ms/LearnAboutSenderIdentification__;!!G_uCfscf7eWS!YNLdSyA8YkLHqmi8DisU8Hz_g4lAFvUm-N7RZHxO_u1CkWa6ZEMYoUBG2so6IcRb57XYSVCvOLtnD-Fx0ic5BOuyGj4A2oxZ2g$  ]
>> 
>> Caution: External email.
>> 
>> 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
> 
> —
> Canada Research Chair in Mathematical and Computational Aspects of
> Solid Mechanics (Tier 1)
> Professor, Department of Mathematics & Statistics
> Hamilton Hall room 409A, McMaster University
> 1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada
> https://urldefense.us/v3/__https://www.math.mcmaster.ca/bourdin__;!!G_uCfscf7eWS!YNLdSyA8YkLHqmi8DisU8Hz_g4lAFvUm-N7RZHxO_u1CkWa6ZEMYoUBG2so6IcRb57XYSVCvOLtnD-Fx0ic5BOuyGj7SJK5sRQ$  | +1 (905) 525 9140 ext. 27243