Атомная доля Pu239 в плутонии

Эта страница предназначена для демонстрации того, как вычислить влияние неопределённости на атомную долю $^{239}\text{Pu}$ в $\text{Pu}$ составе.

import sauna

Реализация опирается на чувствительности функционала к полной макросечению. Следовательно, необходимо получить их из соответствующих результатов расчёта нейтронного транспорта. В этом примере это чувствительности MET1000, полученные с помощью Serpent.

sensitivities = sauna.Sensitivities()
sensitivities.from_serpent('../../models/MET1000_FC_28_sens0.m')

После импорта данных нужно указать нуклиды, которые используются при нормировке общей атомной доли к единице. В данном случае это $\text{Pu}$ изотопы, но те, которые нас интересуют, должны быть исключены, то есть $^{239}\text{Pu}$.

pu_nuclides = [942360, 942380, 942400, 942410, 942420]

Затем атомная доля $^{239}\text{Pu}$ должна быть задана для возмущения. В примере значение примерно равно 0,5244:

ao_239 = 0.5244

Имея значения, теперь можно вызвать функцию для оценки влияния неопределённости, предполагая, что её значение равно 1%.

uncertainty_ao = sauna.Analysis.get_concentration_uncertainty(sensitivities    = sensitivities,
                                                              uncertainty      = 0.01,
                                                              targets          = [942390],
                                                              background_zams  = pu_nuclides,
                                                              fraction         = ao_239,
                                                              fraction_type    = 'ao',
                                                              uncertainty_type = 'normal')

Ожидаемый результат составляет примерно 0,31% для собственно­го значения (eigenvalue).