Hi ,all,
I use meep to simulate grating structure with oblique plane wave
illumination. In order to normalize, the electric field distribution without
structure is calculated. However, I found that the amplitude of electric field
is uniform along x direction but not in y direction. Distribution of amplitude
in y direction is sinusoidal and peak-valley value increase with incident
angle. As I interpreted, the amplitude of electric field of plane wave should
be uniform spatially.
Why is not amplitude in y direction uniform? Any advice will be
appreciated.
(set! force-complex-fields? true)
(define-param withstructure? false)
(define a 100 )
(define wave_length (/ 632.8 a))
(define sx 15.8)
(define sy 50)
(define pml_depth 5)
(define my_material (make dielectric (epsilon 2.25)))
(if withstructure? (define-param grating_material my_material)
(define-param grating_material default-material))
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set-param! resolution 10)
(set! geometry
(list
(make block
(center 0 -12.5)
(size 15.8 25)
(material grating_material))
(make block
(center 0 5)
(size 7.9 10)
(material grating_material))))
(set! pml-layers (list (make pml (thickness pml_depth) (direction Y))))
(define theta (* (/ 75 180) pi))
(define freq (/ 1 wave_length))
(define kx (* freq (sin theta)))
(define (my-amp-func p)(exp (* 0+2i pi kx (vector3-x p))))
(set! k-point (vector3 kx 0 0))
(set! sources (list (make source
(src (make continuous-src (frequency (/ 1 wave_length)) (width (/ 1
wave_length))))
(component Ez)
(center 0 (- (/ sy 2) pml_depth))
(size sx 0 )
(amp-func my-amp-func))))
(run-until (* 200 wave_length)
(at-beginning output-epsilon)
(at-end output-efield-z))
feng
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
