Hello Xinyan Lu, I don't know if you have mistakes in the code, but the transmitted field may be the evanescent wave due to total internal reflection. This field is decaying over distance, so it will disappear when you have a lot larger computation cell. You can refer to optics book in the chapter about TIR （全反射，倏逝波）. Hope this helps.
Allan On Jan 19, 2013, at 10:58 AM, Xinyan Lu <[email protected]> wrote: > Hi, everyone. I'm trying to simulate a phenomenon of light dispersion > throught a photonic crystal, Before that, I want to do a test for a simple > reflection from SiO2 to air. > > But I can never get a total reflection no matter how I augment the incident > angle. Why this happens? > > Thanks! > > And here is my ctl file. > > > (define-param sx 50) > (define-param sy 11) ; the size of the computational cell, not including PML > (define-param dpml 1) ; thickness of PML layers > > (define sxx (+ sx (* 2 dpml))) > (define syy (+ sy (* 2 dpml))) ; cell size, including PML > (set! geometry-lattice (make lattice (size sxx syy no-size))) > (set! geometry (list > (make block (center 0 (* -0.25 syy)) (size sxx (/ syy 2) infinity) > (material (make dielectric (epsilon 2.1609)))) > )) > > (set! pml-layers (list (make pml (thickness dpml)))) > (set! resolution 10) > (define ((pw-amp k x0) x) > (exp (* 0+1i (vector3-dot k (vector3+ x x0))))) > > (define-param fcen 0.8) ; pulse center frequency > (define-param df 0.02) ; turn-on bandwidth > (define-param kdir (vector3 1 20)) ; direction of k (length is irrelevant) > (define k (vector3-scale (* 2 pi fcen) > (unit-vector3 kdir))) ; k with correct length > (define kxcos (vector3-x (unit-vector3 k))) ; direction cosine of k in x > (define kycos (vector3-y (unit-vector3 k))) ; direction cosine of k in y > > (set! sources > (list > > > (make source > (src (make continuous-src (frequency fcen) (fwidth df))) > (component Ez) (center 0 (* -0.25 syy)) (size sx 0) > (amp-func (pw-amp k (vector3 0 (* -0.5 sy)))) > ) > > > )) > > > > (define-param T 400) ; run time > (run-until T > (at-beginning output-epsilon) > (at-end output-efield-z)) > > -- > 鲁欣妍 > 北京航空航天大学中法工程师学院 > 电话：+86 1581 054 1921 > 电话：+33 6 18 48 37 20 > > LU Xinyan > Beihang University - ECPkn > Tel: +86 1581 054 1921 > Tel: +33 6 18 48 37 20 > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
