From c1e2387695aba90577c920a812422e5458200467 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 28 Oct 2011 10:15:25 -0400 Subject: Added meep code for FEM disk analysis --- meep/wg.ctl | 45 +++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 45 insertions(+) create mode 100644 meep/wg.ctl (limited to 'meep') diff --git a/meep/wg.ctl b/meep/wg.ctl new file mode 100644 index 0000000..6f51f36 --- /dev/null +++ b/meep/wg.ctl @@ -0,0 +1,45 @@ +(set! geometry-lattice (make lattice (size 26 26 no-size))) + +(define-param r 10) ; inner radius of ring +(define-param w 1) ; width of waveguide +(define-param sh (- 0 r) ) ; shift from center + +(define-param fcen 0.15) ; pulse center frequency +(define-param df 0.1) ; pulse width (in frequency) + +(set! symmetries (list (make mirror-sym (direction Y)))) + +(set! geometry (list + (make cylinder (center 0 0) (height infinity) + (radius r) (material (make dielectric (index 3.4)))) +; (make block (center 0 (+ sh 0 (/ w -2))) (size infinity 1 infinity) +; (material (make dielectric (epsilon 12)))) + )) + +(set! sources (list +; (make source +; (src (make gaussian-src (frequency fcen) (fwidth df))) +; (component Ez) +; (center (+ r -1) 0)) + (make source + (src (make continuous-src (frequency fcen))) + (component Ez) + (center (+ r -1) 0)) + )) + +(set! pml-layers (list (make pml (thickness 1.0)))) + +(set! resolution 10) + +(use-output-directory) + +;(run-sources+ 300 +; (at-beginning output-epsilon) +; (after-sources (harminv Ez (vector3 (+ sh 2)) fcen df)) +; ) + +(run-until 2600 + (at-beginning output-epsilon) + (at-every (/ 10 fcen) (output-png Ez "-S3 -Zc bluered -a yarg -A ./wg-out/eps-000000.00.h5")) + ) + -- cgit v1.2.3