Dear Steven Johnson and MEEP users,
The first email did not distribute so I am sending this again...
I encounter a convergence problem at the interface of PEC and anisotropic
media. I defined a 2D simulation with PEC triangle on PEC substrate in an
anisotropic medium. When I apply the following code, it ends up with exploded
values at the interface of PEC and anisotropic medium. I have tried to set
'eps-averaging' to false and decreased the courant factor, but it is still not
working.
Please let me know if you have any idea on this problem. Thanks!!
====Code
(set! eps-averaging? false)
(set! Courant 0.01)
(set! geometry-lattice (make lattice (size 10 10 no-size)))
(set! geometry (list
(make block (center -2.5 0) (size 5 10 infinity)
(material (make dielectric (epsilon-diag 3 4 1)
(epsilon-offdiag 0.5 0 0))))
(make block (center 2.5 0) (size 5 10 infinity)
(material perfect-electric-conductor))
(make cone (center -1 0) (radius 3) (height 2) (axis -1 0 0)
(material perfect-electric-conductor))
))
(set! sources (list
(make source
(src (make continuous-src (frequency 1)))
(component Hz)
(center -3 0)
(size 0 8)
)
))
(set! resolution 48)
(set! pml-layers (list (make pml (thickness 1.0))))
(run-until 10
(at-beginning output-epsilon)
(to-appended "hz" (after-time 0 (at-every 1 output-hfield-z)))
)
Best Wishes,
Hanhong Gao
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
