Dear all meep user.
I want to ask about how to define area for flux computation. In examples, I see
tat the length o is twice the wide of wave guide and the wide is 0. Why is it
like that.. Now I am doing this kind of work for photonic crystal . This is my
code
(define notomi (make dielectric (index 12.96)))
(define awal (/ (sqrt 3) -2))
(set! geometry-lattice (make lattice (size 70 24 no-size)))
(set! default-material air) ; everything is airl
(define-param kdir (vector3 1 1))
(set! geometry
(append
(geometric-objects-duplicates (vector3 0 (sqrt 3) 0) 0 4
(geometric-objects-duplicates (vector3 1 0 0) 0 30
(list
(make cylinder (center -14.5 (* -2
(sqrt 3)) 0) (radius 0.35) (height infinity)
(material notomi)))))
(geometric-objects-duplicates (vector3 0 (sqrt 3) 0) 0 3
(geometric-objects-duplicates (vector3 1 0 0) 0 29
(list
(make cylinder (center -14 (+ (* -1
(sqrt 3)) awal) 0) (radius 0.35) (height infinity)
(material notomi)))))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution 20)
(define-param fcen 0.5) ; pulse center frequency
(define-param df 1) ; pulse width (in frequency)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez)
(center 0 -5))))
(define-param nfreq 500) ; number of frequencies at which to compute flux
(define trans ;
(add-flux fcen df nfreq
(make flux-region
(center 0 -2.6) (size 64 0))))
(define-param compute-mode? false)
(run-sources+ (stop-when-fields-decayed
50 Ey
(vector3 0 3.5 0)
1e-3)
(at-beginning output-epsilon)
(during-sources
(in-volume (volume (center 0 0) (size 0 24))
(to-appended "ez" (at-every 0.5 output-efield-z)))))
(display-fluxes trans)
Thank you for the answer
Best regard
andika
____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now.
http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
