Model UO2(5%)+Th+U233 Fuel in MCNPX for SCWR

  • Thread starter Thread starter Aly_19f
  • Start date Start date
  • Tags Tags
    Fuel Mcnpx Model
Click For Summary

Discussion Overview

The discussion revolves around modeling a mixed-oxide fuel consisting of UO2, Th, and U233 using MCNPX for a Supercritical Water Reactor (SCWR). Participants explore the challenges of calculating number densities for the mixed fuel components and seek guidance on methodologies and resources for handling such calculations.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant describes their successful calculation of number density for UO2(5%) but expresses uncertainty about incorporating Th and U233 into the model.
  • Another participant explains that modeling mixed-oxide fuels like (Th,U)O2 involves understanding the weight percentages of each component and adjusting for oxide density variations.
  • Concerns are raised about the method used for calculating atom fractions in mixed oxides, with one participant noting that their calculations yield results greater than one, which is problematic.
  • Participants discuss the importance of accounting for porosity in ceramic materials when determining densities and suggest that textbooks on nuclear engineering may provide relevant examples for burnup calculations.
  • Specific formulas for calculating atomic densities are mentioned, including the relationship between mass, density, and volume, as well as the rule of mixtures for determining densities of mixed oxides.
  • One participant provides theoretical density values for various compounds and discusses how to approximate the atomic ratios for UO2 based on isotopic compositions.
  • Another participant shares a formula for calculating number densities for U233 and Th232 based on the overall fuel density.

Areas of Agreement / Disagreement

Participants express differing views on the methodology for calculating number densities in mixed-oxide fuels, with no consensus reached on the best approach. Some participants agree on the need for adjustments due to porosity and the complexities introduced by multiple isotopes, while others remain uncertain about specific calculations.

Contextual Notes

Limitations include potential missing assumptions regarding isotopic distributions and the specific methodologies for density calculations. The discussion also highlights the complexity of handling multiple isotopes in mixed-oxide fuels, which may not be fully resolved.

Who May Find This Useful

This discussion may be useful for researchers and practitioners in nuclear engineering, particularly those interested in fuel modeling, mixed-oxide fuels, and computational methods in reactor physics.

Aly_19f
Messages
9
Reaction score
1
TL;DR
How do I calculate the number density of each element in an enriched fuel?
Hey there, I'm working on an MCNPX modelling for SCWR using different clads and fuels, the first fuel was UO2(5%) and I have calculated the number density correctly since there was only one vector U.
But now I don't know how top deal with the Th+U233 due to the existence of Thorium.
Can anyone help me or guide me to a similar solved problem (Uranium + another heavy fuel)?
 

Attachments

  • fuel.png
    fuel.png
    11.2 KB · Views: 203
Engineering news on Phys.org
Aly_19f said:
Summary:: How do I calculate the number density of each element in an enriched fuel?

But now I don't know how top deal with the Th+U233 due to the existence of Thorium.
Well, one is asking how to model a 'mixed-oxide' fuel, which could be (U,Pu)O2, (Th,U)O2, (Th,Pu)O2, or even (U,Th,Pu)O2; they all form MO2, where M = Th, U, and/or Pu. The enrichments are given in percent by weight (wt%). It's relatively easy in UO2, since the lattice is the same regardless of U-isotope, but it is fairly straightforward in MOX fuel. One has to do slight adjustments for oxide density of mixtures.

The 3.66% U235 or U233 would imply that the remainder of the metal (Th) is 96.34%. Note the densities to the mixed-oxide and are less than the theoretical density (due to porosity in the ceramic material) and vary slightly due to the differences in mass of the fissile nuclides.
 
Last edited:
  • Informative
  • Like
Likes   Reactions: Aly_19f and Alex A
Astronuc said:
Well, one is asking how to model a 'mixed-oxide' fuel, which could be (U,Pu)O2, (Th,U)O2, (Th,Pu)O2, or even (U,Th,Pu)O2; they all form MO2, where M = Th, U, and/or Pu. The enrichments are given in percent by weight (wt%). It's relatively easy in UO2, since the lattice is the same regardless of U-isotope, but it is fairly straightforward in MOX fuel. One has to do slight adjustments for oxide density of mixtures.

The 3.66% U235 or U233 would imply that the remainder of the metal (Th) is 96.34%. Note the densities to the mixed-oxide and are less than the theoretical density (due to porosity in the ceramic material) and vary slightly due to the differences in mass of the fissile nuclides.
I still don't know how to deal with the problem, in UO2 the sum of the Uranium atom fractions is equal one whatever the enrichment is.
But by applying the same method in the MOX it is always >1.
Is there any textbook or tutorial about burnup calculations?
 
Astronuc said:
The 3.66% U235 or U233 would imply that the remainder of the metal (Th) is 96.34%. Note the densities to the mixed-oxide and are less than the theoretical density (due to porosity in the ceramic material) and vary slightly due to the differences in mass of the fissile nuclides.
Of course, I'm referring to (U,Th)O2.

Aly_19f said:
But by applying the same method in the MOX it is always >1.
Probably the method assumes a common density, which is a reasonable assumption for UO2, but not MOX, but of course, the atom fractions should add to one.

Aly_19f said:
Is there any textbook or tutorial about burnup calculations?
I would think most textbooks for Introductory Nuclear Engineering or Nuclear Reactor Theory would have some background or examples. This is essentially a basic chemistry problem.

The atom density, N, is given by N = ρ Av / M, where ρ is the density, Av is Avogadro's number and M is the molecular mass.

When calculating atomic densities of mixed oxides, MO2, one may apply the rule of mixtures, but must understand the relationship of mass, density and volume, and it can get complicated when each heavy metal element has multiple isotopes, e.g., 233U, 235U, 238U, and there may be 234U and 236U, and with Pu, 238Pu, 239Pu, 240Pu, 241Pu and 242Pu.

Rule of mixtures - https://en.wikipedia.org/wiki/Rule_of_mixtures
https://www.sciencedirect.com/topics/engineering/rule-of-mixture-equation

With mixed oxides, one has to remember that there is one M (Th,U,Pu) for each O2.
So if we want to calculate atomic densities of M, then we need to keep in mind that we would first calculate the molecular density of MO2 from the mass density, and then multiply by the mass ratio M/MO2.

Also, one must remember the physical density is less than the theoretical density, because ceramics do not sinter to 100% of theoretical density. There is some inherent porosity, which is process specific.

The theoretical densities of the compounds interest are:

ThO2, 10.00 g/cm3
UO2, 10.96 g/cm3, although some use 10.97 or 10.98, which might indicate an influence of the proportions of 235U and 238U in the UO2; getting a pure amount of 235U or 238U is exceedingly difficult, or there could be some porosity or other impurity.
PuO2, 11.46 g/cm3

One should note the ratio of M to MO2. For example, for UO2, one could approximate it by 238/(238+2*16) = 238/270 = 0.8815. If one wishes to be more precise, then one would have to account for the proportions of 235U and 238U. For example, for 5% enrichment, the ratio is (0.05*235+0.95*238)/(0.05*235+0.95*238+2*16) = (237.85)/(237.85+32) = 0.8814. So, it's a small difference, since the atomic masses of the isotopes are very close.

For densities of mixtures of oxides, one can use the atom% of the metals, which are close to the wt%. For example, the theoretical density of (U0.8, Pu0.2)O2 is approximately 11.08 g/cm3. Using the densities above for UO2 = 10.96 and PuO2 = 11.46, then multiplying by the atomic fractions 0.8 and 0.2, respectively, one obtains 10.96 g/cm3 * 0.8 + 11.46 g/cm3 * 0.2 = 11.06 g/cm3.

In the above example, given the 3.66% 233U enrichment in (U,Th)O2, and densities of 10.96 g/cm3 and 10 g/cm3 for UO2 and ThO2, respectively, then the theoretical density of the mixture/compound is 0.0366*10.96 g/cm3 + 0.9634 * 10 g/cm3 = 10.04 g/cm3. So, the theoretical density (TD) is slighly greater than that of ThO2. The table one provided has the density at 9.5 g/cm3, from which one would infer that the mixed oxide ceramic has a density of ~0.95 of theoretical, which is a typical target of oxide ceramic nuclear fuel. In the distant past, ceramic densities were lower, ~0.87-0.94 TD.
 
  • Like
Likes   Reactions: Aly_19f
Aly_19f said:
Summary:: How do I calculate the number density of each element in an enriched fuel?

Hey there, I'm working on an MCNPX modelling for SCWR using different clads and fuels, the first fuel was UO2(5%) and I have calculated the number density correctly since there was only one vector U.
But now I don't know how top deal with the Th+U233 due to the existence of Thorium.
Can anyone help me or guide me to a similar solved problem (Uranium + another heavy fuel)?
N/v for fuel is (rho_fuel*N_a/M(th232*(1-0.0366)+u233*(0.0366)))
N/v for u233 is n/v_fuel*0.0366
N/v for th232 is n/v_fuel*(1-0.0366)

These quantities will be in cm^-3
 
  • Like
Likes   Reactions: Aly_19f

Similar threads

  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
3
Views
8K
Replies
10
Views
13K
  • · Replies 7 ·
Replies
7
Views
3K
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
10K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 4 ·
Replies
4
Views
9K
  • · Replies 1 ·
Replies
1
Views
5K