It means the back-reflectance is of the inverse direction compared to you incident beam.
roy On Nov 6, 2007 11:04 AM, fanguofang <[EMAIL PROTECTED]> wrote: > Hi: > Why the reflected flux has the negative and positive? > > creating output file "./bend-flux-denergy-001497.45.h5"... > run 0 finished at t = 1497.45 (29949 timesteps) > flux1:, 0.645, 0.127093598718564, -0.102789149726399 > flux1:, 0.64510101010101, 0.113700231828847 , -0.054823836101576 > flux1:, 0.64520202020202, 0.0987521784729818, 0.0216289088302877 > flux1:, 0.64530303030303, 0.0973125826225362, 0.101522556030181 > flux1:, 0.64540404040404, 0.129662980886524, 0.144391028816318 > flux1:, 0.64550505050505, 0.216171239988684, 0.108368038923313 > flux1:, 0.64560606060606, 0.3700750104319, -0.0275032375188931 > flux1:, 0.64570707070707, 0.590580218163083, -0.240648529190736 > flux1:, 0.64580808080808 , 0.859649223524078, -0.454151126045829 > flux1:, 0.64590909090909, 1.14572538087916, -0.549298050635704 > flux1:, 0.6460101010101, 1.4160014627803, -0.405801661728238 > flux1:, 0.646111111111111, 1.65580772502274, 0.0389828877641945 > flux1:, 0.646212121212121, 1.8900668630813, 0.734790154879933 > flux1:, 0.646313131313131, 2.19889827048105, 1.49840073936785 > > control file: > > (define-param n 3.4) ; Dielectric index of waveguide > > (define-param sx 20) ; size of cell in X direction > (define-param sy 20) ; size of cell in Y direction > > (set! geometry-lattice (make lattice (size sx sy no-size))) > (set! pml-layers (list (make pml (thickness 1.0)))) > (set-param! resolution 10) > > (define-param nfreq 100) ; number of frequencies at which to compute flux > (define-param w 0.5) ; width of waveguide > > (define-param no-bend? false) ; if true, have straight waveguide, not bend > > (set! geometry > (if no-bend? > (list > (make block > (center 0 0) > (size infinity w infinity) > (material (make dielectric (index n))))) > (list > (make block > (center 0 0) > (size 12 w infinity) > (material (make dielectric (index n)))) > ))) > > (define-param fcen 0.65) ; pulse center frequency > (define-param df 0.01) ; pulse width (in frequency) > (set! sources (list > (make source > (src (make gaussian-src (frequency fcen) (fwidth df))) > (component Ez) > (center -5 0) > ))) > > > (define trans ; transmitted flux > (add-flux fcen df nfreq > (if no-bend? > (make flux-region > (center 6.5 0) (size 0 (* w 2))) > (make flux-region > (center 5.5 0) (size 0 (* w 2) ))))) > (define refl ; reflected flux > (add-flux fcen df nfreq > (make flux-region > (center -4.5 0) (size 0 (* w 2))))) > > > (if (not no-bend?) (load-minus-flux "refl-flux" refl)) > > (run-sources+ > (stop-when-fields-decayed 10 Ez > (if no-bend? > (vector3 6.5 0) > (vector3 5.5 0)) > 1e-3) > (at-beginning output-epsilon) > (at-end output-efield-x) > (at-end output-dpwr)) > > > (if no-bend? (save-flux "refl-flux" refl)) > > (display-fluxes trans refl) > > _______________________________________________ > meep-discuss mailing list > [email protected] > http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss > -- Photonics Lab Rochester Inst of Tech Rochester, NY 14623
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
