hi every one
I have written a ctl file to calculate Q factor in a pc cavity . I want to
calculate TE mode in the cavity and I used Ey for src component . but
harminv give a minus q factor. what is wrong?
my cod is :(reset-meep)
(define-param diel_hole 1);dielectric constant of the holes
(define-param diel_slab 7.43);effective dielectric constant of fundamental
mode in the bulk
(define-param res 10);resolution pixels per unit of distance
(define-param rad 0.29);radius of the holes
(define-param numlayers 10);number of layers (10)
(define-param fcen_h 0.3);for harminv - bandgap for h=0.5 r=0.29
triangularhole struct. is between ~0.27 and 0.35
(define-param df_h 0.1);for harminv
(define-param sx 20); x-size of the cell
(define-param sy 20); y-size of the cell
(define-param runtime 5000)
(set! default-material (make dielectric (epsilon
diel_slab)));defaultmaterial: high dielectric
(set! geometry-lattice (make lattice (size sx sy no-size)));lattice sx*sy*sz
(set! geometry
(append;combining list of objects
(geometric-objects-duplicates (vector3 1 0 0 ) (* -1 numlayers)
numlayers
(geometric-objects-duplicates (vector3
0.5 (/ (sqrt 3) 2) 0 ) (* -1 numlayers) numlayers (list
(make
cylinder (center 0 0 0)(radius rad)(height infinity)(material (make
dielectric(epsilon diel_hole)))))))));duplicate whithin duplicate to
makeholes in the slab
(set! geometry (append geometry (list
(make cylinder (center 0)(radius
rad)(height infinity)(material (make dielectric (epsilon
diel_slab)))))));defect in the middle
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen_h )(fwidth 0.1)))
(component Ey)
(center 0.1 0.3 0))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution res)
(run-sources+ runtime;run until sources off and "runtime" more
(at-beginning output-epsilon)
(after-sources
(harminv Ey (vector3 0.1 0.3 0) 0.3 0.1)))
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
