Hi John. I think the gif shows right. In your code, you defines a waveguide with epsilon 11.7. But the background material is air. So, the guide will not confine the wave, it is ok. Some energy will flow away the guide. What is the problem?
2011/12/29 John Bordelon <[email protected]> > 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 > > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss > -- Cícero Julião da Silva Junior ____________________________________ Institute of Physics - UFAL www.ufal.br
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
