------
Steven, et al.,or more likely i'm doing something silly.simple simulation trying to calculate reflectance, transmission and absorption through a cubic layer of spheres (i.e. a single plane of spheres arranged in a square lattice). send a gaussian pulse down the z axis, place the sphere at the center of the simulation and use periodic boundary conditions in x & y and PML at the start and end of the z ordinate. instead of collecting flux in one plane, it was collected at 4 planes at varying distance from the sphere (err, layer of spheres) for both the reflected and transmitted.
what i'm seeing is that R and T are NOT the same at each plane. in particular, at large frequencies, there is some spurious (?) oscillations in R & T which depend on the position of the flux planes. since there is nothing but air between the various flux planes, it appears energy is not conserved.
attached is the ctl file used and the results of the simulation using pc_metal?=true. it's true that if you calculate n & k for this dielectric function, you find 0.036 < n < 0.54 and 2.39 < k < 15 for the frequency range of interest, so the interaction w/ the metal sphere may be violating the Courant condition. i did reduce the resolution from 64 to 32, but these 'wiggles' still start appearing at the same place (i.e. w~1). trying to run w/ higher resolution but computer resources are a problem. btw, if a constant eps is used, say 4, all is acceptable.
i do realize one could use symmetry to reduce the computational burden, but i wanted to make sure i didn't introduce (another?) silly error.
pulling hair,
gp
------
test.ctl
Description: Binary data
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
