D3DEnergetic / FIDASIM

A Neutral Beam and Fast-ion Diagnostic Modeling Suite
http://d3denergetic.github.io/FIDASIM/
Other
29 stars 19 forks source link

Updated 3D functionality #155

Closed alvin-garcia closed 6 years ago

alvin-garcia commented 6 years ago

Luke, I fixed a few bugs for 3D FIDASIM, please have a look. The spectra and NPA outputs are fine with these changes, but the neutron rate is a bit off. I am using nphi = 10 and phi bounds = [4pi/3, 5pi/3]. Below are my results.

test neutron rate = 5.9759e+14 neutrons/s

test3D_np10 neutron rate = 8.9640e+13 neutrons/s

ratio = 6.6665549

Please let me know if you have any suggestions. Thanks.

lstagner commented 6 years ago

was the test case using the same phi bounds? if not then the neutron rate would be for the range phi=[0,2pi) and you would expect the neutron rate to be about 1/6 the test rate.

alvin-garcia commented 6 years ago

Yes, my concern is that there is a small discrepancy and it is off by about 0.67. I varied nphi and n_pfida and the rate did not change. The flux comparisons improved, however, as one would expect.