Dear Meep users:
I wonder if anyone has this experience.
If I obtain one resonant frequency from one broadband pulse source(use
harminv and source +) and then use it as the center wavelength of the next
narrowband pulse source, sometimes the resonant frequencies from this run don't
have this frequency. Does it mean it is related to the convergence issue and
that frequency is not the real resonant frequency?
In addition, I found if I change my resolution to 0.1, then there is no
output fields(the png files show black) even the program works properly.
Regards
Adrian
Attach code as follows:
(define-param Rd 200) ; radius of hexagonal
(define-param L (* 1 Rd)) ; radius of hexagonal
(define-param pad 800) ; padding between waveguide and edge of PML
(define-param dpml 150) ; thickness of PML
(define-param Lz 400) ; z direction length
(define-param Lsub 1200) ; z direction length
(define sxy (* 2 (+ Rd pad dpml))) ; radial size (cell is from 0 to sr)
(define-param ang (/ (* 120 pi) 180)) ; angle of triangle
(define sz (+ Lz Lsub pad)) ; radial size (cell is from 0 to sr)
(set! geometry-lattice (make lattice (size sxy sxy sz)))
(define lam (/ 1 0.0026810081089669))
(define ZnO (make dielectric (index 2.4)))
(define AlO (make dielectric (index 1.8)))
(define SiO2 (make dielectric (index 1.5)))
(set! pml-layers (list (make pml (thickness dpml) )))
(define-param H_tri (* (* 2 Rd) (sin (* 0.5 ang))))
(define-param ay (sin (* 0.5 ang)))
(define-param ax (cos (* 0.5 ang)))
(set! geometry (list
(make block (center (* 0.5 L) 0 (* 0.5 Lz)) (size L L Lz)
(material ZnO) (e1 ax ay 0) (e2 ax (* -1 ay) 0) (e3 0 0 1) )
(make block (center (* -0.5 L) 0 (* 0.5 Lz)) (size L L Lz)
(material ZnO) (e1 ax ay 0) (e2 ax (* -1 ay) 0) (e3 0 0 1) )
(make block (center 0 0 (* 0.5 Lz)) (size L H_tri Lz) (material
ZnO) )
(make block (center 0 0 (* -0.5 Lsub)) (size sxy sxy Lsub)
(material SiO2) )
))
(set-param! resolution 0.05)
(define-param fcen (/ 1 lam)) ; pulse center
frequency
(define-param df 10e-4) ; pulse width (in frequency)
(define ((pw-amp k x0) x)
(exp (* 0+1i (vector3-dot k (vector3+ x x0)))))
(define-param kdir (vector3 0 1 1)) ; direction of k (length is irrelevant)
(define k (vector3-scale (* 2 pi fcen)
(unit-vector3 kdir))) ; k with correct length
(set! sources
(list
; left
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez) (center 0 (* -0.5 sxy) 0) (size sxy 0 sz)
(amp-func (pw-amp k (vector3 0 (* -0.5 sxy) 0))))
; bottom
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey) (center 0 0 (* -0.5 sz)) (size sxy sxy 0 )
(amp-func (pw-amp k (vector3 0 0 (* 0.5 sz)))))
))
(run-sources+ 70000
(at-beginning output-epsilon)
(at-end output-efield-x output-efield-z output-efield-y)
(after-sources (harminv Ez (vector3 (* 0.6 L) 0 (* 0.5 Lz)) fcen
df))
)
______________________________________________________________________________________________________
付費才容量無上限？Yahoo!奇摩電子信箱2.0免費給你，信件永遠不必刪！ http://tw.mg0.mail.yahoo.com/dc/landing_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
