Dear Prof. Johnson and all meep users,
I have tried EM field wavepacket propagation one sublayer of FCC stacked
spheres with a ctl input file given below as input, to ensure validation of
FDTD propagation and to continue to perform FDTD calculation with a lot of
FCC layers photonic crystal.
(set-param! resolution 20)
(set! geometry-lattice (make lattice (size 2 2 22)))
(define nshell (make dielectric (index 3.5)))
(define ncore (make dielectric (index 1.46)))
(set! geometry (append
(list
(make sphere (center 0 1 0) (radius (sqrt 0.5)) (material nshell))
(make sphere (center 0 -1 0) (radius (sqrt 0.5)) (material nshell))
(make sphere (center 1 0 0) (radius (sqrt 0.5)) (material nshell))
(make sphere (center -1 0 0) (radius (sqrt 0.5)) (material nshell))
)))
(set! k-point (vector3 0 0 0))
(set! pml-layers (list (make pml (thickness 1.0) (direction Z))))
(define-param nfreq 501)
(define-param fcen (/ (* 2 3.5)))
(define-param df (* 2 fcen))
(set! sources
(list (make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey)
(center 0 0 10) (size 2 2 0))))
(define trans
(add-flux fcen df nfreq
(make flux-region (center 0 0 -9.4) (size 2 2 0))))
(define refl
(add-flux fcen df nfreq
(make flux-region (center 0 0 9.4) (size 2 2 0))))
(load-minus-flux "refl" refl)
(run-sources+ 500 (at-beginning output-epsilon) (at-every 10
output-efield-y))
;(save-flux "refl" refl)
(display-fluxes trans refl)
I also used reference ctl file in which (set! geometry ) parentheses was
eliminated, and comment of (load- ) and (save- ) was changed.
With this input, I could obtained somewhat valid trans. and refl. spectrum
in the wavelength coordinate about from 3.5 to 30,
but z coordinate of objects was changed to any nonzero value, values of
propagated EM field and following spectrum were diverged,
thinking theoretically, propagation and follwing results should not be
affected by parallel transfomation along any z direction.
What cause this abnormal behavior?
I sincerely hope any help or comment from dear everyone.
Dr. ChangGyun Shin
Professional Researcher
Reformable Display Group, Display Lab.
Samsung Advanced Institute of Technology, Korea
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
