Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread PERRIER-MICHON Augustin
Dear all, I agree that TAO or SNES should be better solutions for fracture analysis using phase field models alone. In my case, the use of TS is not real a choice. It is motivating by later adding new time-dependent physics (like thermal or species diffusion). To be fair, I chose to use grad

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread Pierre Jolivet
On 23 Apr 2025, at 7:28 PM, neil liu wrote:MMG only supports serial execution, whereas ParMMG supports parallel mode (although ParMMG is not as robust or mature as MMG). Given this, could you please provide some guidance on how to handle this in the code? Here are my current thoughts; please let

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread Matthew Knepley
On Wed, Apr 23, 2025 at 2:20 PM Blaise Bourdin wrote: > Hi, > > Typically, phase-field models are formulated as rate independent > unilateral minimization problems of the form > > u_i,\alpha_i = \argmin_{u,\alpha \le \alpha_{i-1}} F(u,\alpha) > > Where i denotes the time step. These are technical

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread Blaise Bourdin
Hi, Typically, phase-field models are formulated as rate independent unilateral minimization problems of the form u_i,\alpha_i = \argmin_{u,\alpha \le \alpha_{i-1}} F(u,\alpha) Where i denotes the time step. These are technically neither DAE nor ODE since there is the only time derivative

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread neil liu
*MMG only supports serial execution, whereas ParMMG supports parallel mode (although ParMMG is not as robust or mature as MMG).* Given this, could you please provide some guidance on how to handle this in the code? Here are my current thoughts; please let know whether it could work as a temporary

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread Stefano Zampini
If mmg does not support parallel communicators, we should handle it internally in the code, always use commself, and raise an error if there are two or more processes in the comm that have cEnd - cStart > 0 Il giorno mer 23 apr 2025 alle ore 20:05 neil liu ha scritto: > Thanks a lot. Pierre. > D

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread neil liu
Thanks a lot. Pierre. Do you have any suggestions to build a real serial DM from this gatherDM? I tried several ways, which don't work. DMClone? Thanks, On Wed, Apr 23, 2025 at 11:39 AM Pierre Jolivet wrote: > > > On 23 Apr 2025, at 5:31 PM, neil liu wrote: > > Thanks a lot, Stefano. > I tried

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread Barry Smith
> On Apr 23, 2025, at 11:47 AM, PERRIER-MICHON Augustin > wrote: > > Dear Mr Smith, > > Thank you for your answer. > > Maybe I do not understand your interrogation because I am not familiar with > fluid flow simulations. Each physics are solved on an entire time step at > every time step.

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread neil liu
Thanks a lot, Stefano. I tried DMPlexGetGatherDM and DMPlexDistributeField. It can give what we expected. The final gatherDM is listed as follows, rank 0 has all information (which is right) while rank 1 has nothing. Then I tried to feed this gatherDM into adaptMMG on rank 0 only (it seems MMG wor

Re: [petsc-users] Mixing PETSc Parallelism with Serial MMG3D Workflow

2025-04-23 Thread Pierre Jolivet
> On 23 Apr 2025, at 5:31 PM, neil liu wrote: > > Thanks a lot, Stefano. > I tried DMPlexGetGatherDM and DMPlexDistributeField. It can give what we > expected. > The final gatherDM is listed as follows, rank 0 has all information (which is > right) while rank 1 has nothing. > Then I tried

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread PERRIER-MICHON Augustin
Dear Mr Smith, Thank you for your answer. Maybe I do not understand your interrogation because I am not familiar with fluid flow simulations. Each physics are solved on an entire time step at every time step. This is a quasi static problem with an incremental loading. Later I will add time de

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread Barry Smith
I am unfamiliar with this simulation type, so I have some elementary questions. Is each "physics" taking a fractional time-step so that the total time integrated (in a time-step) is 1/2dt + 1/2dt or is one physics stepping the entire time-step and the other physics "fixing up some of the v

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread 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-s

Re: [petsc-users] Staggered solver phase field

2025-04-23 Thread Blaise Bourdin
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 quas

[petsc-users] Staggered solver phase field

2025-04-23 Thread 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 so