Dear Meep users,
I have a problem with divergence of the field at source point
inside material with low index of refraction. The field at this point
increases at every period. I tried to decrease Courant number but it does
not help, as a determination of dispersive material instead of simple case.
(set! geometry-lattice (make lattice (size 15 15 no-size)))
(define-param f 0.5)
(define-param eps_inf 0.2)
(define-param mu_inf 0.2)
(set! Courant 0.05)
(set! geometry (list
(make sphere (radius 3) (center 0 0 0)
(material (make dielectric
(epsilon eps_inf) (mu mu_inf)
)))))
(set! sources (list
(make source
(src (make continuous-src
(frequency f)
(width 20)))
(component Ey)
(center 0 0 0)
)))
(set! pml-layers (list (make pml (thickness 2.0))))
(set! resolution 25)
(define T (/ 1 (* 4 f)))
(use-output-directory "meta-out")
(define (Field)
(print "Field:, " (meep-time) ", " (get-field-point Ey (vector3 0
0 0))", "(get-field-point Ex (vector3 0 0 0)) "\n"))
(run-until 1000
(at-every 0.02 Field))
Could someone explain me how to avoid this divergence?
I would be much appreciate your answers and suggestions
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
