Dear Steven and meep useres.
I want to simulate the GDV (group velocity dispersion) of PCF (photonic
crystal fiber) with meep. Could you please tell me whether meep can do?
If it can do, I want to try
it with the following geometric objects. How should I add this ctl-file?
; Parameters =========================================================
(define-param eps 12) ; dielectric constant of core
(define-param a 2.4) ; lattice constant
(define-param r (* 0.315 a)) ; radius of holes
(define-param N 4) ; layer of clad
(define-param pad 1) ; padding between clad and edge of PML
(define-param dpml 0.5) ; PML thickness
;=========================================================== Prameters
; Primitive cell =====================================================
(define-param sx (+ (* (* 2 N) a) (* 2 r) pad (* 2 dpml)))
; size of cell in
; x direction
(define-param sy (+ (* (* 2 N) (* (sqrt 3) (/ a 2)))
(* 2 r) pad (* 2 dpml))) ; size of
cell in
; y direction
(set! geometry-lattice (make lattice (size sx sy no-size)))
(set! default-material (make dielectric (epsilon eps)))
;====================================================== Primitive cell
; Strucrure ==========================================================
(define (clad v w x y)
(geometric-objects-duplicates (vector3 x y) 0 N
(geometric-object-duplicates (vector3 v w) 1 N
(make cylinder (center 0)
(radius r)
(height infinity)
(material air)))))
(set! geometry
(append
(append
(clad a
0
(/ -1 (/ 2 a))
(/ (sqrt 3) (/ 2 a)))
(clad (/ -1 (/ 2 a))
(/ (sqrt 3) (/ -2 a))
a
0)
(clad (/ -1 (/ 2 a))
(/ (sqrt 3) (/ 2 a))
(/ -1 (/ 2 a))
(/ (sqrt 3) (/ -2 a))))
))
(set! pml-layers (list (make pml (thickness dpml))))
(set! resolution 30)
;=========================================================== Strucrure
; Source =============================================================
(set! sources
(list (make source
(src (make continuous-src
(wavelength 0.88)))
(component Ez)
(center 0 0))))
(set! symmetries
(list (make mirror-sym (direction X) (phase 1))
(make mirror-sym (direction Y) (phase 1))
))
;============================================================== Source
; Run=================================================================
(run-until 100
(at-beginning output-epsilon)
(at-every 1.0 (output-png Ez "-Zc bluered")
))
;================================================================= Run
Thank you very much!
Yuta
Yuta Wakayama
E-mail: [EMAIL PROTECTED]
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
