Dear Steven and other Friends,
I met a problem.
I calculate the transmittance of W1 PhC waveguide. I randomly change
the radius of each air hole.
To varify the correction, I have added a monitor at a single point to
observe the feild at the output.
I surprisely found that when time is lager than 1000s, the field is
suddenly change to larger than 100000000. but before time=1000s, the field is
normal(<1).
I check my code:
3D structure
the purlse freqs (0.25, 0.3)
the resolution is 16, a is around 400nm,
delta X(Y,Z)=400/16=25nm<32nm=400/0.3/3.45/12
the pml is 1, which take 16 grids
the input light source is putted 16 grids away from the pml layer,
the output monitor is putted 16 grids away from the pml layer
I can not find where is the problem.
Prof. Johnson, could you please give me some help?
thank you very much for your time~
PS:
another question: how to set a guaussian shape source in 3D situation?
here are my code:
================================================
(define ((my-amp1 sigma k) y)
(exp (- (* 0+2i pi (vector3-dot k y))
(/ (vector3-dot y y) (* 2 sigma sigma)))))
(define ((my-amp2 sigma k) z)
(exp (- (* 0+2i pi (vector3-dot k z))
(/ (vector3-dot z z) (* 2 sigma sigma)))))
;;;;define the sources
(set! sources (list
(make source
(src (make gaussian-src (frequency fcen) (fwidth df)))
(component Ey) (amp-func (my-amp1 ny (vector3 (+ 1.5 (* -0.5 size-x)) 0 0)))
(amp-func (my-amp2 h (vector3 (+ 1.5 (* -0.5 size-x)) 0 0))) (center (+ 1.5 (*
-0.5 size-x)) 0 0)
(size 0 ny h)))) ;add source E in z direction
================================================
Sincerely
Ran
_______________________________________________
meep-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
