hi Hodjat df is different in your code df =0.5 in Fig.3.c flux freq. range is (fcen-df/2) to ( fcen + df/2)
Onder 2011/6/5 <[email protected]> > Send meep-discuss mailing list submissions to > [email protected] > > To subscribe or unsubscribe via the World Wide Web, visit > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss > or, via email, send a message with subject or body 'help' to > [email protected] > > You can reach the person managing the list at > [email protected] > > When replying, please edit your Subject line so it is more specific > than "Re: Contents of meep-discuss digest..." > > > Today's Topics: > > 1. correction of the mistake in calculating normalized > transmission through a cavity in a square lattice PC (hodjat hajian) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 5 Jun 2011 08:45:59 -0700 (PDT) > From: hodjat hajian <[email protected]> > Subject: [Meep-discuss] correction of the mistake in calculating > normalized transmission through a cavity in a square lattice PC > To: Meep Discussion <[email protected]> > Message-ID: <[email protected]> > Content-Type: text/plain; charset="iso-8859-1" > > Hi dear users. > I want to calculate the transmission (Fig.3c) that is illustrated in this > paper (Phys. Rev. B, Vol. 54, Num. 11 , pages:7837-7842), but my result is > different. > The program that I've used in Meep is as below, but the result is different > from Fig. 3c. Help me to correct the mistake please. > Best. > Hodjat > > ?;some parameters to describe the geometry: > (define-param eps 1.843) > (define-param r 0.2) ; radius of holes > (define-param Nx 7) ; number of holes in X direction > (define-param Ny 11) ; number of holes in Y direction > (define-param rd 0.2) ; r of defect > > ?; the cell dimensions: > (define-param pad 4) ; padding from last hole and pml edge > (define-param dpml 1) ; pml thickness > (define sx (+ (* 2 (+ pad dpml r)) Nx -1)) > (define sy (+ (* 2 (+ dpml r)) Ny > -1)) > (define wx (+ Nx -1 (* 2 r))) > (define wy (+ Ny -1 (* 2 r))) > > ?;******************************* > > (set! geometry-lattice ( make lattice (size sx sy no-size))) > > (set! geometry > (list > (make cylinder (center 0 0 0) (radius r) (height infinity) (material (make > dielectric (epsilon eps)))) > ) > ) > > (set! geometry (geometric-objects-lattice-duplicates geometry)) > > (set! geometry > (append geometry > > (list > (make cylinder (center 0 0 0) (radius rd) (height infinity) (material air)) > ) > (list > (make block (center (/ (+ sx wx) 4) 0) (size (* 0.5 (- sx wx)) infinity > infinity) (material air)) > ) > (list > (make block (center (/ (+ (* -1 wx) (* -1 sx)) 4) 0) (size (* 0.5 (- sx > wx)) infinity infinity) (material air)) > ) > > (list > (make block (center 0 (/ (+ sy wy) 4)) (size infinity (* 0.5 (- sy wy)) > infinity) (material air)) > ) > > (list > (make block (center 0 (/ (+ (* -1 wy) > (* -1 sy)) 4)) (size infinity (* 0.5 (- sy wy)) infinity) (material air)) > ) > > > ) > ) > > ?;********************************* > (set! pml-layers (list (make pml (thickness dpml)))) > (set-param! resolution 20) > > ?;some parameters to describe the source: > (define-param fcen 0.35) > (define-param df 0.20) > (define-param nfreq 500) > (set! sources (list > (make > source (src (make gaussian-src (frequency fcen) (fwidth df))) > (component Ez) (center (+ dpml (* -0.5 sx)) 0) (size 0 sy)) ; > )) > > (set! symmetries (list > (make mirror-sym (direction Y) (phase -1)) > )) > > (define trans (add-flux fcen df nfreq > (make flux-region (center (- (* 0.5 sx) dpml 0.5) 0) (size 0 sy)) > )) > > (run-sources+ > (stop-when-fields-decayed 50 Ez (vector3 (- (* 0.5 sx) dpml 0.5) 0) 1e-3) > (at-beginning output-epsilon) > ) > > (display-fluxes trans) > -------------- next part -------------- > An HTML attachment was scrubbed... > URL: < > http://ab-initio.mit.edu/pipermail/meep-discuss/attachments/20110605/0849804f/attachment.html > > > > ------------------------------ > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss > > End of meep-discuss Digest, Vol 64, Issue 4 > ******************************************* > -- Pamukkale University Electrical-Electronics Eng.Dept. 20070 Denizli TURKEY Tel: (+90 258) 296 3213
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
