Dear Meep users,
I want to simulate total internal reflection （TIR） effects using Meep. The
attached code was modified from the
plane wave sample codes provided in previous posts.
The code is a 633nm laser plane wave illuminated from glass to air in a
larger incident angle than TIR critical angle.
What I should see is that all the incoming wave will be reflected along the
air/glass interface and a evanencent wave will propagate along it.
But the plane wave still exists in the simulation results.
Could you please give me some idea what is wrong, and I would be appreciate
if you can provide a right code.
Thank you very much.
Ryan
======= code =====
(define-param sy 10) ; the size of the computational cell NOT including PML
inthe y direction
(define-param sx 2) ; the size of the computational cell including PML in
the x direction
(define-param dpml 2) ; the thickness of PML layer
(define-param svert (+ sy (* 2 dpml)))
(define-param frec (/ 1 0.633))
(define-param res 25)
(define-param theta (/ (* 1 pi) 8))
(define kx (* frec (sin theta)))
(define (my-amp-func p)(exp (* 0+2i pi kx (vector3-x p))))
(set! k-point (vector3 kx 0 0))
(set! geometry-lattice (make lattice (size sx svert no-size)))
(set! geometry (list
(make block (center 0 -3) (size 4 4 infinity)
(material (make dielectric (epsilon 2.31)))))); a
glass inserted
(set! sources (list
(make source
(src (make continuous-src (frequency frec)))
(component Ez)
(center 0 0)
(size sx 0)
(amp-func my-amp-func))))
(set! pml-layers (list (make pml (thickness dpml) (direction Y))))
(set! resolution res)
(run-until 200
(at-end (output-png Ez "-Zc
/usr/share/h5utils/colormaps/dkbluered")))
====end code====
--
―――――――――――
[EMAIL PROTECTED]
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
