Dear meep users,
I am trying to calculate reflection from glass in the
wavelength range 400- 700nm.I used 1.5 as refractive index of glass.but
reflection from glass is increasing as wavelength increases .Here is my
.ctl file
*****************************************************************************
(define-param fcen 0.2222) ; pulse center frequency a=100nm
(define-param df 0.13636)
(define-param resl 20);
(define-param thick 1); thickness of glass in meep units
(define-param dpml 1) ; thickness of PML layers
(define-param sx (* 10 thick) ) ;
(define sx0 (+ sx (* 2 dpml)))
(define sy sx);
(define sy0 (+ sy (* 2 dpml))) ;
(define sz sx)
(define sz0 (+ sz (* 2 dpml))) ;
(define-param sample 10)
(define-param nfreq 100)
(set! k-point (vector3 0 0 0))
(define-param decayby 1e-8)
(define-param structure? false) ;
(define background (make dielectric (epsilon 1.0001 )))
(set! eps-averaging? false)
(set! force-complex-fields? false)
;(set! Courant 0.25)
(set! geometry-lattice (make lattice (size sx0 sy0 no-size)))
(set! geometry (if structure?
(list
;(make block (center 0 0) (size sx sy) (material background))
(make block (center 0 0) (size infinity thick infinity)
(material (make dielectric (index 1.5))))
)
(list
(make block (center 0 0) (size sx sy) (material
background))
)
))
(set! pml-layers (list
(make pml (thickness dpml) (direction X))
(make pml (thickness dpml) (direction Y))
)
)
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ez)
(center 0 (* thick -2)) (size sx 0)
)
)
)
(set! resolution resl)
(define incident; transmitted
flux
(add-flux fcen df nfreq
(make flux-region
(center 0 (* thick 3)) (size sx 0)
)
)
)
;)
(define bottom
(add-flux fcen df nfreq
(make flux-region
(center 0 (* thick -3)) (size sx 0)
)
)
)
(if structure?
(load-minus-flux "bottom_flux" bottom)
)
(run-sources+ (stop-when-fields-decayed 50 Ez (vector3 0 (* thick 3) 0)
decayby))
(if (not structure?)
(save-flux "bottom_flux" bottom)
)
(display-fluxes incident bottom)
Where could the problem be ?
Thank you !
_______________________________________________
meep-discuss mailing list
meep-discuss@ab-initio.mit.edu
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
