Dear Steven,
I want to calculate the transmission spectrum in PC slabs: (Refer to the Meep
tutorial)
- First, I define the structure: My structure is PC slab with triangular
lattice of air holes (refractive index of slab - n; radius of hole - r; slab
thickness - h). The input/output waveguide is formed by removing a row of air
holes. Between input and output waveguide is a resonator.
- The center frequency of Gaussian source, that is a resonant frequency of the
resonator, is located at one edge of the cell just outside the PML.
- Finally, I compute the flux spectrum at the other end of the computational
cell, but before the PML.
My computational cell: 30x20x0.5; After making the structure, I have some codes
are below. Could you help to check some errors?
(set! pml-layers (list (make pml (thickness dpml))))
(set-param! resolution 20)
(define-param fcen 0.27915) ; pulse center frequency
(define-param df 0.01) ; pulse width
(define-param nfreq 500) ; number of frequencies at which to compute flux
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth
df)))
(component Hz)
(center -14.5 0 0)
(size 0 w h)) ; h: slab thickness, w: width
of waveguide
))
(set! symmetries (list (make mirror-sym (direction Y) (phase -1))))
(define trans
(add-flux fcen df nfreq
(make flux-region
(center 14 0 0) (size 0 (* 2 w) h))))
(run-sources+ (stop-when-fields-decayed
50 Hz
(vector3 14 0 0)
1e-3)
(at-beginning output-epsilon)
(during-sources
(in-volume (volume (center 0 0 0) (size 30 0 h))
(to-appended "hz-slice" (at-every 0.4 output-hfield-z)))))
(display-fluxes)
I thank you in advance,
Chad Huskow
____________________________________________________________________________________
The fish are biting.
Get more visitors on your site using Yahoo! Search Marketing.
http://searchmarketing.yahoo.com/arp/sponsoredsearch_v2.php_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
