[petsc-users] Solvers for phase field fracture

Matthew Knepley knepley at gmail.com
Thu Jan 4 04:50:37 CST 2018

On Thu, Jan 4, 2018 at 1:12 AM, Avtar Singh <er.bhullersingh at gmail.com>

> Hello Petsc Users,
> I am solving a fortran framework to simulate phase-field fracture problem.
> Upto the failure point the code is working fine. But as the crack-phase
> starts to propagate, there is abrupt changes in the Stiffness Matrix, so
> determinant becomes nearly equal to zero.  Hence, the code crashes.

I am not sure what formulation you are using, but that does not happen with
the Francfort-Marrigo-Bourdin formulation I am familiar with:



> I tried with mumps package with pc lu. Also tried, superlu, gmres and
> fgmres with jacobi, bjacobi, lu. But the problem still persist.
> Can anyone suggest, Which solver and preconditioner should i use?
> Thank you
> --
> Avtar Singh
> Research Scholar
> Multiscale Mechanics and Multiphysics Lab
> IIT Roorkee, Roorkee

What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.caam.rice.edu/~mk51/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180104/5af766fb/attachment.html>

More information about the petsc-users mailing list