Dear Steven and other meep user, I simulated a photonic crystal waveguide. The field is good at the beginning, but blows up in the end. The structure is just made of dielectric. I don't understand why the field blow up in the end. Could you help me? Thank you in advance.
Li
The snapshots of Hz field at 100 time unit and 5000 time unit are attached. The
control file is as follows,
________________________________________________________________________________________________________________________________
(define-param sx 30) ; size of cell in X direction
(define-param sy 20) ; size of cell in Y direction
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! default-material (make dielectric (index 3.45)))
(define-param r 0.29) ; the rod radius in the bulk crystal
(define obj-list1 (geometric-object-duplicates (vector3 0.5 (* 0.5 (sqrt 3))) 0
40
(make cylinder (material air) (center 0 0) (radius r) (height
infinity))))
(define obj-list2 (geometric-object-duplicates (vector3 0.5 (* -0.5 (sqrt 3)))
0 40
(make cylinder (material air) (center 0 0) (radius r) (height
infinity))))
(use-output-directory)
(set! geometry
(append
(geometric-objects-duplicates (vector3 1 0) 0 20 obj-list1)
(geometric-objects-duplicates (vector3 -1 0) 0 20 obj-list1)
(geometric-objects-duplicates (vector3 1 0) 0 20 obj-list2)
(geometric-objects-duplicates (vector3 -1 0) 0 20 obj-list2)
(list (make block (material (make dielectric (index 3.45))) (center 0 0)
(size sx 0.7)))
))
(define-param fcen 0.236) ; pulse center frequency
(define-param df 0.2) ; turn-on bandwidth
(set! sources (list
(make source
(src (make continuous-src (frequency fcen) (fwidth df)))
(component Hz) (amplitude 1e10)
(center -13 0)
(size 0 1)
)))
(set! pml-layers (list (make pml (thickness 2.0))))
(set-param! resolution 20)
(run-until 30000
(at-beginning output-epsilon)
(at-every 100 output-hfield-z)
)
<<attachment: hz-000200.00.jpg>>
<<attachment: hz-005000.00.jpg>>
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
