Boundary conditions

[Tissot11]

I used following for particle boundary conditions (#5131)

boundary.particle_lo = Absorbing periodic
boundary.particle_hi = thermal periodic

using following for fields

boundary.field_lo = PML periodic
boundary.field_hi = PML periodic

results in error saying that I should choose neumann pec . Choosing following

boundary.field_lo = neumann periodic
boundary.field_hi = neumann periodic

works. But I do not know why PML or absorbing_silver_mueller won't work with particle boundary conditions above. Moreover, I know neumann from mathematics but I wonder what it means in WarpX e.g. in the context of PML, open or silver_mueller. I could not reproduce my results using open boundary conditions from another code.

[RemiLehe]

Thanks for your question.

In order to answer your question, we would need a few more details about the type of field solver that you are using. In particular:

• is this a purely electrostatic simulation? (i.e. does your simulation script use do_electrostatic)
• is this an electromagnetic simulation in which you are asking WarpX to compute the initial fields using a Poisson solver? (i.e. does your simulation script use initialize_self_fields)

If you could post your input script (or a modified version thereof, if it contains sensitive information), that would be even easier.

[Tissot11]

As you can in the input deck in #5131 that it is an electromagnetic simulation and I do not use do_electrostatic also initialize_self_fields. But noting the initialize_self_fields could be on by default, I specifically put in another simulations species.initialize_self_fields=0 but it still gives the same error (see the attached file)

errWarpX-2770081.txt outWarpX-2770081.txt