[petsc-users] Solvers for phase field fracture

[Matthew Knepley]

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

> 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:

  https://bitbucket.org/litianyi/dynamic-gradient-damage

    Matt


>
> 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>