PML boundary conditions for embedded boundary

[Omstavan]

Hi,

I am trying to use warpx.2d to make a rectangular simulation domain with a circular embedded boundary in the middle which I can make by using the warpx.eb_implicit_function in the input parameters. I need all the boundaries to be PML and charge absorbing. I am aware that the default boundary condition is PML but I am not sure if it applies to embedded boundaries as well. I also need it to be particle absorbing which I am aware that the embedded boundaries are so by default. I am also open to getting this done in the rz version of warpx if possible, but again I am not sure we can have the rectangular outer boundaries PML in the rz version. Further, I am hoping for slightly more complex structure in the middle in the future, so having a PML condition for it would be useful.

Additionally, the warpx.eb_implicit_function description says "the physics simulation area is where the function value is negative" but the following input gives an error "amrex::Abort::0::AMReX_EB2_Level.H: Domain is completely covered !!!" which is odd given the circle is much smaller. Using the negative of that does seem to run.

#################################
####### GENERAL PARAMETERS ######
#################################
max_step = 1000
amr.n_cell =  736 736
amr.max_grid_size = 256
amr.blocking_factor = 512
amr.max_level = 0
geometry.dims = 2
geometry.prob_lo = -1.e1   -1.e1    # physical domain
geometry.prob_hi =  1.e1    1.e1
my_constants.ev_to_K = 11604.51812
my_constants.T = 5
#################################
####### Boundary condition ######
#################################
boundary.field_lo = pml pml
boundary.field_hi = pml pml
warpx.eb_implicit_function = 0.25 - x*x - y*y

Can someone please clarify these for me?

[ax3l]

@roelof-groenewald just explained this to me: you need to quote the spaces:

warpx.eb_implicit_function = "0.25 - x*x - y*y"

otherwise it evaluates to 0.25
Does that help?