Antti
Hi,
I tried to use ORIGEN as a stand-alone code to estimate tritium production in a aluminiun-lithium alloy.
Fluxes for irradiation are based on a separate MCNP model. I input them with THERM, RES and FAST factors as described in the manual.
The input is basically as follows (including only lithium):
#origens
0$$ a4 33 a11 71 e t
tritium production
3$$ 28 a3 0 0 a11 2 a16 2 a33 0 e
4** 0.876 57.9 3.50 0 t
35$$ 0 t
56$$ 10 10 1 a5 1 a10 0 a13 2 a14 3 a15 3 a18 3 e
57** 0 a3 1e-05 0.1118285 e t
Case total production
0.0001 Al-Li
59** 2.07e11 2.07e11 2.07e11 2.07e11 2.07e11 2.07e11
2.07e11 2.07e11 2.07e11 2.07e11
60** 679 1358 2037 2716 3395 4074 4753 5432 6111 6790
66$$ a1 2 a5 2 a9 2 e
73$$ 30070 30060
74** 484 33
75$$ 4 4
t
56$$ 0 0 a10 1 e t
56$$ 0 0 a10 2 e t
56$$ 0 0 a10 3 e t
56$$ 0 0 a10 4 e t
56$$ 0 0 a10 5 e t
56$$ 0 0 a10 6 e t
56$$ 0 0 a10 7 e t
56$$ 0 0 a10 8 e t
56$$ 0 0 a10 9 e t
56$$ 0 0 a10 10 e t
56$$ f0 t
end
=opus
LIBUNIT=33
TYPARAMS=NUCLIDES
UNITS=BEQUERELS
LIBTYPE=ALL
TIME=HOURS
NPOSITION=1 2 3 4 5 6 7 8 9 10 end
end
=opus
LIBUNIT=33
TYPARAMS=NUCLIDES
UNITS=GRAMS
LIBTYPE=ALL
TIME=HOURS
NPOSITION=1 2 3 4 5 6 7 8 9 10 end
end
The code runs fine and results seem logical, but when I compare final results with MCNP production rates for Li-6 and Li-7 ORIGEN calculates over 50 times higher values.
What could be the
explanation? I know that the resonance flux in this case seems very high, but
that is not a mistake here. However, when I calculate the same material in another type of spectrum, MCNP and ORIGEN results are much closer.
Is there any other way of inputting the spectrum?