However, I found that the field was oscillated when it was small. Then, the
computation process spent lots of time.
In my opinion, the computation cell acts as a cavity due to a small
reflection of PML.
The problem is that you didn't use the PML correctly: it needs to
*overlap* your structure in order to absorb fields propagating in the
dielectric material. But you did:
On Tue, 25 Apr 2006, Jun She wrote:
(define-param sx (+ x (* 2 dpml)))
(define-param sy (+ y (* 2 dpml)))
(define-param sz (+ 3 (* 2 dpml)))
(set! geometry-lattice (make lattice (size sx sy sz)))
(set! geometry
(list
(make block
(center 0 0 0)
(size x y z)
(material (make dielectric (epsilon eps))))))
That is, your block ends just as the PML is beginning. This means that
there is a reflection at the boundary of the block and the PML.
Change the block size to, say, (size infinity infinity z).
Cordially,
Steven G. Johnson
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
