Hello,
I am working with OLEDs and PhCs and running some simulations with the
following ctl file (attached).
I am interested in simulating an impulse like dipole which turns on then
shuts off almost immediately but you are able to see the effects for a
certain length of time. Please see the following animation for further
clarification: http://www.lumerical.com/movs/led_textured_fdtd.mpg
Right now, the dipole source remains on throughout the meep simulation.
I was hoping you could help suggest some ideas that would enable the result
to look more like what is seen in the animation.
Thanks very much.
; WLctl - Coupled cavity system proposed by Bose and Yang w/ Wong
;
; ---------------------------------------------------
;Ranojoy Bose May 12, 2008
;Calculate transmission and cavity modes
; ---------------------------------------------------
(define-param eps1 1) ; dielectric constant of Air (n_air = 1)
(define-param eps2 2.25) ; dielectric constant of Glass
(define-param a 1) ; lattice parameter 420 nm
(define-param eps3 3.61); SiNx
(define-param eps4 3.24); ITO
(define-param eps5 3.0625); Organic
(define-param r 0.36) ; radius of air holes
(define-param t1 0.33) ; glass rod height
(define-param h 3) ; Glass layer
(define-param h1 1) ; SiNx layer
(define-param h2 0.33) ; ITO Layer a = 600 nm
(define-param h3 0.3) ; Olayer
(define-param dpml 1) ; PML thickness
(define-param sx 20) ; size of cell in x direction
(define-param sy 20) ; size of cell in y direction
(define-param sz 10) ; very large to avoid PML layer
(set! progress-interval 100)
(set! eps-averaging? false)
(set! geometry-lattice (make lattice (size sx no-size sz)))
(set! pml-layers (list (make pml (thickness dpml))))
;(meep-field-set-boundary fields Low Z Metallic);
(set! geometry
(append
(list (make block (center 0 0 2.23) (size sx sy 1.5) (material (make
dielectric (epsilon eps2)))))
(list (make block (center 0 0 .95) (size sx sy 1) (material (make
dielectric (epsilon eps3)))))
(list (make block (center 0 0 .315) (size sx sy .33) (material (make
dielectric (epsilon eps4)))))
(list (make block (center 0 0 0) (size sx sy .3) (material (make
dielectric (epsilon eps5)))))
(list (make block (center 0 0 -2.15) (size sx sy 1) (material
metal)))
(geometric-objects-lattice-duplicates
(list
(make cylinder (center 0 0 1.615) (radius r) (height 0.33)
(material (make dielectric (epsilon eps2))))
) 1 1 80)
))
(set-param! resolution 8)
(use-output-directory "OLEDedit6")
(define-param fcen 0.55) ; pulse center frequency
(define-param df 0.02) ; pulse width (in frequency) corresponds to 1 fs
(define-param nfreq 300) ; number of frequencies at which to compute flux
; false = transmission spectrum, true = resonant modes:
(define-param compute-mode? true)
(if compute-mode?
(begin
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)
(start-time 0.5) (end-time 0.500001)))
(component Ex)
(center 0 0 0)
(size 0 0 0)
)))
;(set! sources (list
; (make source
; (src (make gaussian-src (frequency fcen) (fwidth df)
(end-time 0.05)))
; (component Ez)
; (center 0 0 4)
; (size 0 0 0)
; )))
; (set! symmetries
; (list (make mirror-sym (direction Z) (phase +1)) ; z-even
; (list (make mirror-sym (direction Y) (phase +1)) )) ;
y-even
;(make mirror-sym (direction X) (phase -1)); x-odd
;))
(run-sources+ 1
(at-beginning output-epsilon)
(after-sources (harminv Ex (vector3 0) fcen df)))
(run-until (/ 100 fcen) (at-every (/ 1 fcen 10000) output-efield-x
;output-hfield-x
;output-hfield-y
))
))
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
