[petsc-users] Periodic boundary conditions using swarm

[MIGUEL MOLINOS PEREZ]


Dear all,

I am working on the implementation of periodic bcc using a discretisation (PIC-style). I am working with a test case which consists on solving the advection of a set of particles inside of a box (DMDA mesh) with periodic bcc on the x axis.

My implementation updates the position of each particle with a velocity field, afterwards I check if the particle is inside, or not, the supercell (periodic box). If not, I correct the position using bcc conditions. Once this step is done, I call Dmswarmmigrate.

It works in serial, but crashes in parallel with MPI (see attached nohup file). I have checked some of the Dmswarmmigrate examples, and they looks similar to my implementation. However they do not use periodic bcc.

I am missing any step in addition to Dmswarmmigrate?

Best regards
Miguel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240418/03299c9f/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nohup.out
Type: application/octet-stream
Size: 7966 bytes
Desc: nohup.out
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240418/03299c9f/attachment.obj>