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
