Hey all, I am working on generalising a 2D model into 3D. To start, I am just trying a single dielectric slab acting as a waveguide. The source is put in the middle of the slab. However, the guide does not confine the waves.
Any help is greatly appreciated as I have been fumbling around with this for about a month now. Here is the code: (define-param fcen .6452) ; frequency of source (define-param fspread .01) ; spread in gaussian freq range (define-param guide_z .2) (define-param guide_r .2) ;this is the radial thickness of ring (define-param size_z 1) (define-param size_x 8) (define-param size_y 1) (define-param pml_thickness 2) (set! geometry-lattice (make lattice (size (+ pml_thickness size_x) (+ pml_thickness size_y) (+ pml_thickness size_z)))) (set! geometry (list (make block (center 0 0 0) (size infinity guide_r guide_z) ; size is given as (size x y z) (material (make dielectric (epsilon 11.7)))) )) (set! sources (list (make source (src (make gaussian-src (frequency fcen) (fwidth fspread))) (component Ez) (center 0 0 0 ) (size .05 .05 .05) (cutoff 15)))) (set! pml-layers (list (make pml (thickness 1)))) (set! resolution 20) (use-output-directory) (run-until 200 (at-beginning output-epsilon) (to-appended "ez" (at-every 5 output-efield-z))) ------------------------------------------------------------------------------------------------------------------ Also attached is a gif showing what I mean. Thanks so much, John Bordelon, Louisiana State University
<<attachment: ezslice.gif>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
