[EMAIL PROTECTED]
MIME-Version: 1.0
Content-Type: text/plain; charset=iso-8859-1
Content-Transfer-Encoding: 8bit
Message-ID: <[EMAIL PROTECTED]>
Dear Meep Users and Steven
I am simulating a periodic dielectric structure with
periodicity 0.24um and indices as 3.248+_ Delta.
However, In my results, I donot find a change in
behaviour with altering delta or the time width of
gaussian source. Is there something, i need to
correct? I desperately need a solution. I thank you
for your kind review of the problem.
I am also attaching the ctl file.
;Simulate propagation of Gausian Pulse in dielctric
slab
; AND Calculating Electric Field Component at the end
of a dielectric slab
;Consider propagation in Z direction and width of slab
in X direction
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
(reset-meep)
(define-param a 1e-6) ; base unit
in meters (1e-6 = 1 um)
(define-param wavelength-center 1.535) ; central
wavelength in micrometer
(set! eps-averaging? false)
;define the largest wavelength in the gaussian
spectrum, in base units
(define fc (/ 1 wavelength-center)) ; source center
frequency
(define-param delta_index 0.003) ; the difference in
refractive index
(define base_index 3.248) ; the base refractive index
(define-param eps-lo 1) ; the surrounding
low-dielectric material
(define-param w 1.1) ; the thickness of the waveguide
(set! Courant 0.5)
(define mat1 (make dielectric (index (+ base_index
delta_index))))
(define mat2 (make dielectric (index (- base_index
delta_index))))
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Define the structure and the computational cell
; Here we define the size of the computational cell.
Since it is 2d,
; it has no-size in the z direction.
(set! geometry-lattice (make lattice (size 450 5
no-size)))
; the default-material is what fills space where we
haven't placed objects
(set! default-material (make dielectric (epsilon
eps-lo)))
; a list of geometric objects to create structures in
our computational cell:
; (in this case,we have a periodic medium with high
and low refractive indices)
(set! geometry(append
(list
(make block
(center -200.2 0 0)
(size 50 w infinity)
(material (make dielectric (index
base_index))))
(make block
(center 200.1 0 0)
(size 49.8 w infinity)
(material (make dielectric (index
base_index)))))
(geometric-object-duplicates (vector3
0.24 0 0) 1 1460
(make block ; a dielectric block (a
rectangle)
(center -175.14 0 0)
(size 0.12 w infinity) ; block is finite only
in x & y direction
(material mat1)))
(geometric-object-duplicates(vector3 0.24
0 0) 1 1460
(make block
(center -175.02 0 0)
(size 0.12 w infinity)
(material mat2)))
))
(set-param! resolution 20)
(set! pml-layers (list (make pml (thickness 1))))
(set! sources (list (make source
(src (make gaussian-src (frequency fc) (width (
* 1.5 303.03 ))))
(component Ez) (center -224 0)
(amplitude 10) (size 0 w))))
(use-output-directory
"/home/shiv/Desktop/examples/FDTD/Output")
(define (Power)
(print "Power:," (meep-time)", " (flux-in-box X
(volume (center -220 0) (size 0 2)))", " (flux-in-box
X (volume (center 224 0) (size 0 2))) "\n" ))
(define trans ;transmitted flux
(add-flux fc 0.0033 2000
(make flux-region
(center 220 0) (size 0 (* w 2)))))
(run-until 5500 (at-every 1 Power))
(display-fluxes trans)
Regards
Shiv
Bring your gang together - do your thing. Go to
http://in.promos.yahoo.com/groups
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
