Dear Dr. Steven and Meep Users,
I am a new Meep user and trying to calculate band diagram of an EBG unit cell
in 3D for a wave incident normal to the surface. It has a dielectric layer (5
mm x 5 mm x 1.2 mm) with a ground plane ( 5mm x 5 mm x 0.2 mm) and a patch (4.8
mm x 4.8 mm x 0.2 mm) and a center via of r = 0.2 mm. I have periodic boundary
on X and Y ; PML on the Z direction. I use a gaussian source (Ey) 5 mm above
the patch. Band diagram looks wrong when compared with the band diagram
calculated using HFSS (shown in the tutorial document contained in the link). I
am trying to calculate band only in X direction.
Can someone point out what is wrong with my code (given below) or point me in
the right direction ?
Thanks a lot !
Arun.
(set! geometry-lattice (make lattice (size 4 4 30)))
(set! pml-layers (list
(make pml (direction Z) (thickness 10))))
(set-param! resolution 10);
(set! k-point (vector3 0 0 0))
(set! ensure-periodicity true)
(set! eps-averaging? false)
(set! geometry (list (make block (center 0 0 0) (size 5 5 1.2)
(material (make dielectric (epsilon 10.2)))) ;
substrate
(make cylinder (center 0 0 0) (radius 0.2) (height 1.2) (material
metal)) ; via
(make block (center 0 0 0.7)
(size 4.8 4.8 0.2) (material metal)) ; patch
(make block (center 0 0 -0.7)
(size 5 5 0.2) (material metal)))) ; ground
(set! sources (list
(make source
(src (make gaussian-src (frequency 0.25) (fwidth 1.5)))
(component Ey) (center 0 1 5))))
(define-param k-interp 19)
(run-k-points 200 (interpolate k-interp (list (vector3 0 )
(vector3 0.5 ) )))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
