Dear All There seems to be some problem when defining sources and structures in cylindrical dimensions: When defining a source or a structure that is centered somewhere on the r-z plane and stretches until r=0, then either the dielectric function (for a structure) or the field (for a source) are lower around r=0 than they are defined
This problem can be solved for structures http://www.mail-archive.com/[email protected]/msg01354.html But, this solution doesn't work for sources, because sources cannot be defined to be larger than the computational cell. Perhaps someone has a hint how to overcome this problem? I suspect that this discontinuity leads to inaccurate results, in certain conditions. (I give here a small code that shows the problem, and attached an image created from it. Zooming lets you see the discontinuity on the top part of the image (that is, r=0) (use-output-directory) (define-param wlength 4) (define-param field_comp Er) (set! dimensions CYLINDRICAL) (set! geometry-lattice (make lattice (size 20 infinity 20))) (set-param! m 0) (set! eps-averaging? false) (set-param! resolution 5) (set! sources (list (make source (src (make continuous-src (wavelength wlength))) (component field_comp) (center 0 0 0) (size 20 infinity no-size) ))) (define sample_T (* 0.1 wlength)) (run-until (* 20 wlength) (at-beginning output-epsilon) (to-appended "er" (at-every sample_T output-efield-r)) (to-appended "ep" (at-every sample_T output-efield-z)) )
<<attachment: er.t000.png>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
