Total atomic density change in fuels with different id

User Support
1

Dear experts, I have determined two different volumes of fuel in the openmc input file, let’s say both of them contain U235 nuclide. How can I observe the total atomic density change for these?

uo2 = openmc.Material(material_id=1,name="uo2")
uo2.add_element("U", 1, percent_type="ao", enrichment=4.25)
uo2.add_element("O", 2)
uo2.set_density("g/cc", 10.4)
uo2.volume = 252 * pi * fuel_radius**2

uo3 = openmc.Material(material_id=6, name='UO2/60+ThO2/40 %5 enrichment')
uo3.add_nuclide('U234', 0.00000298656706683134)
uo3.add_nuclide('U235', 0.000734315535331843)
uo3.add_nuclide('U238', 0.0139490086042382)
uo3.add_nuclide('Th232', 0.00912327894861938)
uo3.add_nuclide('O16', 0.0476191793105125)
uo3.volume = 78 * pi * fuel_radius**2

dif
dif640×960 26.8 KB

1 Like
2

@pshriwise @paulromano

1 Like
3

As far as I know, openmc does not do this process. You can handle the total calculation with the help of post processing. Or through a formula. These experts can help you with this issue. @gridley @Daedalus @Shimwell

4

Not really sure from the question but perhaps a depletion simulation

Here is an example

Code block 24 in the example notebook is similar to your plots

1 Like
5

Here, the change in total atomic density in depletion of any nuclide for materials with two different IDs is desired.He wants to find the total thorium depletion in material id 1 and 6 to show them both in a single graph.For example, they stated this as TOT_ADENS in Serpent.

6

Yes you can track multiple materials in a depletion simulation.

Be sure to set the material.depletable =True for the materials you want to track. Then the results will be in the h5 file for both materials

1 Like
7

Yes, it saves the material with two different IDs to the h5 file. For example, 78 pin 1 id, 234 pin 3 id. But here we want to find the total change in these. That is, the total Th232 depletion in the number of pins 78 + 234. Do we need to produce a different formula to calculate this?

8

I guess if you have the change for the two individual pins you can combine those two results to get the total. Perhaps I’ve misunderstood something here. But would you not just add them up?

2 Likes
9

Yes, I think it is a logical solution to add them together to find the change in Th232 atomic density on pins with two different IDs.I will combine the list of atomic density changes at id 1 and id 6 with the addition process.