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://github.com/bourdin/mef90

Blaise


On Apr 23, 2025, at 10:20 AM, PERRIER-MICHON Augustin <augustin.perrier-mic...@ensma.fr> wrote:

[You don't often get email from augustin.perrier-mic...@ensma.fr. Learn why this is important at https://aka.ms/LearnAboutSenderIdentification ]

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://www.math.mcmaster.ca/bourdin | +1 (905) 525 9140 ext. 27243

Reply via email to