Dear colleague:
I just use meep to calculate the loss of the 300micrometer
straight waveguide.The loss is very big.about just 2e-5 of the source
power.Is it normal?
What is wrong with it? I think it the loss should be very very
small just for the straight waveguide.
The following is my mode.
(define-param n 2.7) ;
(define-param pad 2) ; padding between waveguide and edge of PML
(define-param dpml 1) ; thickness of PML
(define-param fcen 0.6452) ; pulse center frequency
(define-param df 0.01) ; pulse width (in frequency)
(define-param w 0.5) ;
(define-param l 300) ;
(set! resolution 10)
(define sx (+ w (* 2 (+ pad dpml)))) ; cell size
(define sy (+ l (* 2 (+ pad dpml)))) ; cell size
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! pml-layers (list (make pml (thickness dpml))))
(set! geometry
(list
(make block (center 0 0 ) (size w infinity )
(material (make dielectric (index n)) ))
))
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ex)
(center 0 (- (/ sy 2) 2) ))))
(define-param sum1 0)
(define ((do-sum1 box))
(set! sum1 (+ sum1
(electric-energy-in-box box))) )
(define (output1) (print "electric-energy= " sum1 "\n"))
(define-param sum2 0)
(define ((do-sum2 box))
(set! sum2 (+ sum2
(electric-energy-in-box box))) )
(define (output2) (print "electric-energy= " sum2 "\n"))
(run-sources+
(stop-when-fields-decayed 50 Ez
(vector3 0 (- 2 (/ sy 2) ) )
1e-3)
(at-every 10
(do-sum1
(volume (center 0 (- 2 (/ sy 2) ) )
(size w 0) ) )
(do-sum2
(volume (center 0 (- (/ sy 2) 5 ))
(size w 0) ) )
)
(at-end output1 output2 )
)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
