Dear Mr.Steven Johnson and MEEP's users,
I tray to simulate the transmission trough a 2D triangle structure with
dispersive material (Ag).
The problem I have is that the field is blowing up at higher resolutions
i.e. over 28.
I have tried to change resolution 30-900 or Courant factor : 0.4-0.001.
Decreasing Courant factor reduce field value but it's still over e20.
Does anyone have an idea about what may be causing this problem or
what parameters I should adjust? Thanks in advance!
Best Regards,
Piotr Wrobel
[EMAIL PROTECTED]
University of Warsaw
This is my .ctl file
:
;==============================================================================================
(set! geometry-lattice (make lattice (size 16 1 no-size)))
(set-param! resolution 28)
(set! eps-averaging? false)
;(set! Courant 0.001)
(define-param a_value 1e-6)
(define con (/ (* 2 pi 2.99792458e8) a_value))
(define-param e0 3.7)
(define-param g (/ 9.617120e13 con)) ;gamma in meep units
(define-param wp (/ 1.346390e16 con)) ;wp in meep units
(define-param omega1 1e-20)
(define materialdysp
(make dielectric (epsilon e0)
(polarizations
(make polarizability
(omega omega1) (gamma g) (delta-epsilon (/ (* wp wp) (* omega1
omega1))) ;delta eps =wp^2/omega^2
))))
(set! geometry (list
(make block (center 0 -0.39) (size 0.2 1 infinity)
(material materialdysp)
(e1 -1 0.57735 )
(e2 -1 -0.57735 ) )
(make block (center 0 0.39) (size 0.2 1 infinity)
(material materialdysp)
(e1 1 0.57735 )
(e2 1 -0.57735 ) )))
(define-param fcen 2) ; pulse center frequency
(define-param df 0.1) ; pulse width (in frequency)
(define-param nfreq 10) ; number of frequencies at which to compute flux
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey) (center -7 0))))
(set! pml-layers (list (make pml (thickness 1.0))))
(define trans ; transmitted flux
(add-flux fcen df nfreq
(make flux-region
(center 7 0 ) (size 0 0.8 0) (direction X))))
(run-until 300
(at-beginning output-epsilon)
(at-end output-efield-y))
(display-fluxes
trans)
;=================================================================================================
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
