How to handle ElectrostaticParticleContainer after mesh refinement (AMR)?

[Afraid31]

Hi,

I'm currently studying the ElectrostaticPIC example to understand how the interaction between fluids and particles is handled.

I noticed that the particle container is constructed like this:
ElectrostaticParticleContainer myPC(geom, dm, grids, rr);
This makes me wonder: in an adaptive mesh refinement (AMR) simulation, when the grid hierarchy is updated (e.g., after regridding), the myPC object will not automatically reflect those updated grids or dm, since they are passed by value during construction.
Does this mean that, in AMR particle simulations, I need to reconstruct the ElectrostaticParticleContainer each time the grid hierarchy changes?
Is this the recommended way to handle particles in AMReX when using adaptive mesh refinement? Or is there a way to update the internal mesh information of an existing particle container?

Thanks in advance for your help!

[atmyers]

If you use that constructor, then you will need to change the grids every time the hierarchy changes. This is mostly for single level simulations or simulations where the particles use a different set of grids than the fluid. However, there is another constructor that takes an AmrCore*, and that one will automatically track changes to the fluid grids. You can see it used here. You'll still need to call Redistribute() on the particle container after the grids change. The AMR Advection tutorial with tracer particles linked above is a good example to look at for using particles in an AMR context.