HomerReid / scuff-em

A comprehensive and full-featured computational physics suite for boundary-element analysis of electromagnetic scattering, fluctuation-induced phenomena (Casimir forces and radiative heat transfer), nanophotonics, RF device engineering, electrostatics, and more. Includes a core library with C++ and python APIs as well as many command-line applications.
http://www.homerreid.com/scuff-em
GNU General Public License v2.0
126 stars 50 forks source link

RF bug #122

Closed akzy closed 7 years ago

akzy commented 8 years ago

Hi Homer,

Seemed to have found a bug in the RF module - Depending on the plot range i decide (keeping the sample frequency the same), the results of the electric field seem to differ. This first plot is when i chose the range 0-200um in z. rf bug z 200um

...and the 0-300um doesnt seem to show these here but instead at 280um. rf bug z 300um rf bug z 300um_numab I noticed this whilst doing contour plots in mathematica - whilst it could be my interpolation or some other maths im doing, i would argue against this as other methods of producing this field data (such a K. Singers Essential Numerical Tools http://www.uni-ulm.de/~ksinger/ent/) do not produce this on the same mathematica sheet.

I'll also provide a copy of my commands im using for scuff and the settings https://github.com/akzy/scuff/tree/patch-3 I would upload the actual data files but they're few hundred MB. Just quick clarification of my settings, I used grid.py to make the .epfile (s) based upon the RAM i have available. compute.sh runs scuff with those settings and Convert.py is used to make it friendly with my mathematica sheet (probably irrelevant to you). The .portcurents file is made in such a way to simulate 1V being applied to the port.

If i can be of anymore help please let me know.