Help regarding understanding calculation of One-group cross section

151 views
Skip to first unread message

Khurrum

unread,
Dec 28, 2015, 1:12:14 PM12/28/15
to OpenMC Users Group
Hi every one,

sorry if the question I am going to ask seems too straight forward but I am asking for your help.

for any geometry, let's say pincell, I want to calculate the MICROSCOPIC cross section for fission for U-235 isotope.

I have set up a tally for fission score on nuclide U-235 (let us assume that it gives a value of 0.52 in tallies.out file) and another tally to calculate the flux for whole fuel (assuming it gives a value of 12.5 in tallies.out) material (as OpenMC gives an error if I try to apply filter of nuclide U-235 on fuel material to tally the flux). where to go from here?

looking forward to your reply.

Khurrum

unread,
Dec 30, 2015, 6:42:15 AM12/30/15
to OpenMC Users Group
Ping!!!!!!

Will Boyd

unread,
Dec 30, 2015, 11:20:53 AM12/30/15
to Khurrum, OpenMC Users Group
Khurrum,

If you are interested in computing multi-group cross sections using OpenMC, you may consider using the "openmc.mgxs" Python module that was developed specifically for this purpose and included in the v0.7.1 release. The "openmc.mgxs" module can compute macroscopic and/or microscopic cross sections for one, some or all nuclides in a material, cell, distribcell, or universe tally filter. To get started using "openmc.mgxs," consider looking at the example IPython Notebooks and the API documentation on the website here. The Notebooks in particular provide a lot of examples of how to compute macro/micro cross sections for real pin cell and assembly problems.

If instead you wish to compute your microscopic cross sections manually using your own inputs and data processing scripts, then I think you are following the correct direction with the process you described earlier. You will need a reaction rate tally and flux tally in the *same* "domain" (material, cell, or universe filter). The flux tally cannot be defined on a per-nuclide basis - only the reaction rate tally should be defined for one or more nuclides - the flux tally cannot be defined for a single nuclide and is instead tallied for all nuclides across the tallied "domain."

After running a simulation with the reaction rate and flux tallies, you can then divide their values to compute the *macroscopic* multi-group cross section for each nuclide in the reaction rate tally. If you wish to obtain microscopic cross sections, you must divide out the nuclide density (in units of at/b-cm) which you can find in the "summary.h5" output file.

Again, all of this can be done for you with just a few lines of Python code if you choose to the "openmc.mgxs" module included in the Python API.

Best,
Will

On Wed, Dec 30, 2015 at 6:42 AM, Khurrum <khurrum...@gmail.com> wrote:
Ping!!!!!!

--
You received this message because you are subscribed to the Google Groups "OpenMC Users Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to openmc-users...@googlegroups.com.
To post to this group, send email to openmc...@googlegroups.com.
Visit this group at https://groups.google.com/group/openmc-users.
For more options, visit https://groups.google.com/d/optout.



--
Will Boyd
Nuclear Science & Engineering
Massachusetts Institute of Technology
Reply all
Reply to author
Forward
0 new messages