Hello, My simple meep example is coming up with an incorrect result.I found some posts on the list regarding similar problems (incorrect reflection coefficient calculation), where the recommendation was to move the flux planes away from the source, by at least a wavelength. This did not solve the problem for me.
In my example, I have a 2D simulation, with periodicity in X, and PML in Y. The scatterer is a dielectric slab with epsilon 9. The source is a plane wave source at the upper PML interface (the angle of incidence is 30 degrees).
The field output looks great (see attached png; resized smaller for posting to list).
To calculate the fluxes, I simple setup flux planes in front of and beyond the slab, and run the computation twice; with and without the slab.
The coefficents are clearly wrong: transmission: 0.647382725227634, reflection: 3.21374994463727, loss: -2.8611326698649 code: http://www.pastebin.us/19072 The command to run that code, which also computes the coefficients is: meep no-scatterer\?=true compute-flux\?=true test.ctl && meep no-scatterer\?=false output-finalfield\?=true compute-flux\?=true output-finalfield\?=true test.ctl | tee lastrun.out Best Regards, Matt
test2.png
Description: test2.png
_______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
