Dear Pr. Steven and all Meep users
First I want to thank you all for your postes because it really helps.
I am beginner user of Meep and I am trying to find the diagram band of a simple
rectangular waveguide. The problem is that i got the same frequencies for all K
points! I checked my units, the resolution, the source and I haven't found the
error yet.
I worked on 2D and i defined K in z direction (because the waveguide is uniform
in Z direction). The unit of distance that I chose is a=1mm. the dimension of
the waveguide are width=22.86mm length=11.43mm. The frequency used is in order
of GHz so in Meep it has to be in order of 0.01 and K is in order of 100, so in
Meep it has to be order of 0.01 also.
Please I need help so if any one can help me I would appriciate that so much!
Thank you in advance:)
The code :
(define-param dpml 4); width of pml layer
(define-param w 6) ; width of metallic boundaries
(define-param a 22.86);
(define-param b 11.43);
(define-param sx (* 2 (+ (/ a 2) w dpml))); size of cell in x direction
(define-param sy (* 2 (+ (/ b 2) w dpml))); size of cell in y direction
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! geometry
(list
(make block (center 0 0) (size (+ a w) (+ b w) infinity)
(material perfect-electric-conductor))
(make block (center 0 0) (size a b infinity)
(material (make dielectric (epsilon 1))))))
(set-param! resolution 5)
(set! pml-layers (list (make pml (direction X) (thickness dpml))))
(set! pml-layers (list (make pml (direction Y) (thickness dpml))))
(define-param fcen 0.06) ; pulse center frequency
(define-param df 0.1) ;
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey) (center 0 0) (size 0 b)))))
(define-param kmin 0)
(define-param kmax 0.1)
(define-param k-interp 10)
(run-k-points 300 (interpolate k-interp (list (vector3 0 0 kmin) (vector3 0 0
kmax))))
RESULTS
For K= 0 0 0
harminv0 0.0215337528554796
harminv0 0.06466384783274905
harminv0 0.0876000615675477
For K= 0 0 0.01 to K= 0 0 0.1 i got the same results
harminv 0.02253547268004267
harminv 0.06479032574739023
harminv 0.08764276787837490 _______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
