particle tracking on a DA

[Berend van Wachem]

Hi,

Yes I use PETSc to do particle tracking with DA's. But it is not very
simple and I get bad scaling. I have done geometric decomposition of the
problem and are keeping "ghost cells". The cells are the size of the
lookup for collisions.

Questions you should address are,

1. What kind of particle-particle interactions?
2. How many particles?
3. Is there a fluid phase as well?

I think the communication with the scheme you propose would be major if
the particle numbers go over 500.000

Berend.

> Hi All,
>
> Is anyone aware of PETSc particle tracking code?  I'm looking for
> something that will do Lagrangian particle tracking in parallel for
> 2D fluid flow on a regular cartesian mesh.
>
> If this doesn't exist, I may try to create it---quick and dirty
> style.  My thinking was to a master-slaves parallelism: keep a
> complete directory of the particle locations on node 0 then send the
> entire directory to each node after each timestep to do an update.
> Each node updates any points that fall into it's subdomain, flags
> those particle entries, then returns the entire directory to node 0.
> Node 0 merges the updates into the master directory.
>
> This isn't scalable but it is simpler to program.  Does anyone have
> any suggestions about strategy, implementation, PETSc tools that
> would facilitate this, etc?
>
> Cheers
> Rich
>
> Richard Foa Katz     NSF IRFP Postdoc
> http://www.damtp.cam.ac.uk/user/rfk22
>
>
>