Hi Steven and Meep's user,
When i calculate a surface plamson waveguide with metal-air-metal
structure, fields always deverge at the metal-air interface nearby (and
insider) the PML.
The problem also exists in some other FDTD programs. Could you give some
comment and advice?
The following ctl file runs in Meep 0.9.
Best regards
Zhichao
(define-param sx 6) ; size of cell in X direction
(define-param sy 4) ; size of cell in Y direction
(set! geometry-lattice (make lattice (size sx sy no-size)))
(define-param wy 0.5) ; the size of x window
(define-param shifty 0)
(define-param ir 0.5)
(define-param a0 100e-9)
(define-param c0 3e8)
(define-param wavelgth 1550e-9)
(define-param fcen (/ a0 wavelgth)) ; pulse center frequency
(define-param df 0.1) ; pulse width (in frequency)
(define normalfreq (* (* 2 pi) (/ c0 a0)))
(define-param EpsA1 1.000000e+000)
(define-param EpsB1 1.214794556259851e+014)
(define-param EpsC1 1.0e5)
(define-param EpsDelta1 (* 1.143359534478220e+016 1.143359534478220e+016))
(define (epsraw->delta-epsilon a b c delta f)
(/ delta c))
(define (epsraw->gamma a b c delta f)
(/ (/ b a) f))
(define (epsraw->omega a b c delta f)
(/ (sqrt (/ c a)) f))
(define x-material
(make dielectric (epsilon 1)
(polarizations
(make polarizability
(omega (epsraw->omega EpsA1 EpsB1 EpsC1 EpsDelta1
normalfreq))
(gamma (epsraw->gamma EpsA1 EpsB1 EpsC1 EpsDelta1 normalfreq))
(delta-epsilon (epsraw->delta-epsilon EpsA1 EpsB1 EpsC1 EpsDelta1
normalfreq)))
)))
(set! geometry (list
(make block
(center 0 0 0)
(size sx sy infinity)
(material x-material))
(make block
(center 0 shifty 0)
(size sx wy infinity)
(material air))
; (make block
; (center 0 0 0)
; (size sx sy infinity)
; (material air))
; (make block
; (center 0 shifty 0)
; (size sx wy infinity)
; (material (make dielectric (epsilon 12))))
))
(set! sources (list
(make source
(src (make continuous-src (frequency fcen)))
(component Hz)
(center 0 shifty))))
(set! pml-layers (list (make pml (thickness 1.0))))
(set-param! resolution 32)
(run-until 200
(at-beginning output-epsilon)
(to-appended "hz" (at-every 0.6 output-hfield-z))
; (at-end output-hfield-z)
)
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
