HI every one,

I am an undergraduate student and doing my final year project on Photonic
Crystal with Heterostructure. I am new user of MEEP. when i calculate the
transmission spectra, i have got some problems. the structure is something
like :

ooooxxxxoooo
ooooxxxxoooo
ooooxxxxoooo

o and x are all rods only with different epsilons. i do a flux-region at
the end of the structure but before PML. run CTL-file but the output just
gave me  0 flux at every time step. i dont understand why. have i done
something wrong?

Could somebody please help me in this problem?

============
CTL-file
============

(reset-meep)

(set! geometry-lattice (make lattice (size 21 21 no-size)))

(set! geometry
    (append
     (geometric-object-duplicates (vector3 1 0) 0 17
        (make cylinder (center -8.5 2.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 12)))))
     (geometric-object-duplicates (vector3 1 0) 0 17
        (make cylinder (center -8.5 1.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 12)))))
     (geometric-object-duplicates (vector3 1 0) 0 17
        (make cylinder (center -8.5 0.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 12)))))
     (geometric-object-duplicates (vector3 1 0) 0 17
        (make cylinder (center -8.5 -0.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 12)))))
     (geometric-object-duplicates (vector3 1 0) 0 17
        (make cylinder (center -8.5 -1.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 12)))))
     (geometric-object-duplicates (vector3 1 0) 0 17
         (make cylinder (center -8.5 -2.5) (radius 0.4) (height infinity)
          (material (make dielectric(epsilon 12)))))

     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 2.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))
     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 1.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))
     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 0.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))
     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 -0.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))
     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 -1.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))
     (geometric-object-duplicates (vector3 1 0) 0 5
        (make cylinder (center -2.5 -2.5) (radius 0.4) (height infinity)
         (material (make dielectric(epsilon 6)))))))

(set! sources (list
              (make source
                (src (make gaussian-src (frequency 0.25) (fwidth 0.4)))
                 (component Ez)
                 (center -9 0)
                 (size 0 6))))

(set! pml-layers (list (make pml (thickness 1.0))))
(set! resolution 10)
(define trans (add-flux 0.25 0.4 100 (make flux-region (center 9 0) (size
0 12))))
(use-output-directory)

(display-fluxes trans)

(run-sources+
        (stop-when-fields-decayed 50 Ez
                (vector3 9 0) 1e-3)
        (at-beginning output-epsilon)
        (at-end output-efield-z))




_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss

Reply via email to