[Gfs-devel] Boundary Conditions for particles

[Gaurav T. <gau...@gm...>]

HI Stephane,

I and Daniel have been working on particulate object a bit.
Specifically we were thinking of implementing boundary conditions for
particles.

We are starting with periodic boundary conditions.
Here is what I have been planning to do:

1. Create an object: GfsBcParticle with the name of the eventlist as
an input so that this object may have access to the GfsParticleList
2. Construct a list of particles to which BC has to be applied and
remove these particles from the particlelist and create a
GtsSListContainer containing these particle objects. (This may be
avoided by using the third point directly)
3. Constructing a GSList of Particle_send structures containing
box->id through which they are migrating to the neighboring boundary
box, direction and the particle object.
4. Write methods such as box_send_bc and box_send_rcv computing
information such as how many particles and what particles are passing
through box-1 to the neighboring box-2. If box->neighbor[d] is a
boundary then method send_particles is called.
5. In send_particles if box_neighbor[d] exists and
GFS_BOUNDARY_PERIODIC(box->neighbor[d])->matching is the box to which
it must go. Particle is again added to the ParticleList with
coordinates modified to account for periodicity.

Please suggest if this is alright or a better algo could be devised.

Sincerely,
Gaurav.