I'm working on a comparison study using SERPENT 1.1.14 and SCALE6.0. I'm basing my comparison on the report "Burnup Calculation Capability in the PSG2/SERPENT Monte Carlo Reactor Physics Code" by Jaakko Leppanen and Maria Pusa.

As part of the study, I am trying to compare the mass of various isotopes such as PU-239, Xe-135, Zp-237, ect. The isotopic mass is being compared as burnup increases from 0-51000MWd. My difficulty is determining the form of output in the SERPENT file. The output "<input>_dep.m" lists a MAT_fuelp1r1_ADENS and MAT_fuelp1r1_MDENS matrix. It also shows a formula for total mass as being TOT_MASS(:,j) = TOT_MASS(:,j) + MAT_fuelp1r1_VOLUME(j).*MAT_fuelp1r1_MDENS(:,j);

So for each burnup step, do I determine the total mass of the nuclide by multiplying the volume (cm^3) with the mass density (gm/cm^3)? This does not seem to give the result I am looking for.

Thanks for your help.

Respectfully,

Tim Wyant