What is the uncertainty in the activity of Cesium 137 over time?

AI Thread Summary
The discussion focuses on calculating the uncertainty in the activity of Cesium 137 over time, specifically from an initial measurement of 60 mCi in 1981 to a present measurement of 26.84 mCi in 2016. Participants suggest using the error propagation rule to estimate uncertainty, with a reasonable assumption of 1% error in the original activity. There is some debate about the appropriate uncertainty in time, with estimates ranging from 1 month to 6 months based on the lack of specific dates. Ultimately, a conservative estimate of uncertainty leads to reporting the present activity as 27 ± 0.5 mCi. The conversation emphasizes the importance of understanding the context of measurements and applying statistical principles correctly.
Smigglet
Messages
2
Reaction score
0
Hello, one of my friends came over to get some help on a homework problem, and it went past over my head I am pretty good with stats, just don't know how to approach this to help him out.

So a radionuclides' activity (Cesium 137) is measured in the year of 1981, (no date just year), and its 60mCi, and its present activity is measured on 02/05/2016.

so i calculated the present act to be 26.84mCi.

however, how do i find totaI uncertainty in the activity.

I was thinking may be it was (InitialA^2 + FinalA^2)^2? However, there does seem to be error in the decay time due to not having a date opposed to just a year.

any help?
 
Engineering news on Phys.org
Hello Smigglet, :welcome:

If the original activity is given as 60 mCi it will be hard to give an answer in 4 digits. 1% is a reasonable estimate of the error in original activity.
And the uncertainty in time is 6 months.
Apply the error propagation rule $$\left (\Delta f(x,y) \right )^2 = \left ( \partial f\over \partial x\right )^2 \left (\Delta x\right )^2 + \left ( \partial f\over \partial y\right )^2 \left (\Delta y\right )^2$$

And PF greetings to your friend
 
Thank you for the reply :)

How would i go about applying the error propogation rule this scenario, i normally have detector data, or number of counts and time that use i the propagation rule. Not too sure how i would apply it here, and excuse me for dumb my question, i cannot see how the uncertainty in time is 6 months if anything id think it would be the difference from date to date, assuming the measurement was taken 02/05/1981, so i would think the uncertainty would be 1 month since we don't know if the measurement was on january.

I think I am confusing myself here lol
 
Like most years, 1981 runs from January 1 to December 31; so if you use ##t_0## June 1 you're at most 6 months off.
The text as you render it (no date just year) gives me no reason to expect that ##t_0## is May 2, 1981.
It's a conservative estimate for the error in the ##t_0## (a better estimate might be the sigma for a uniform distribution -- 3.5 months)

If someone gives you 60 for the activity, you may assume it's not 61 and not 59, so ##60.0 \pm 0.5## is a reasonable interpretation.

You are confusing yourself if you interpret the 60 as a number of counts (in which case ##\sqrt {60}## would be the estimated standard deviation)

Anyway, you have two contributions of about 1 % so reporting ##\ 27\pm 0.5## mCi or just plain ##27## mCi would be reasonable IMHO.

Like the pirates say: it's more what you'd call "guidelines" than actual rules.

Smigglet said:
Not too sure how i would apply it here
If ##\ \displaystyle {\ f = A_0\, 2^{t-t_0\over t_{1/2}}} \ ## then you can differentiate and calculate. I just did the differentiation numerically by calculating the power of 2 for ##t_0 = ## 1-1-81, 1-6-81, 1-12-81 and adding 1% (in ##A_0## ) and 1% (in the power of 2) in quadrature: about 1.4% of 27, so 0.5 :smile:

--
 
Hello everyone, I am currently working on a burnup calculation for a fuel assembly with repeated geometric structures using MCNP6. I have defined two materials (Material 1 and Material 2) which are actually the same material but located in different positions. However, after running the calculation with the BURN card, I am encountering an issue where all burnup information(power fraction(Initial input is 1,but output file is 0), burnup, mass, etc.) for Material 2 is zero, while Material 1...
Hi everyone, I'm a complete beginner with MCNP and trying to learn how to perform burnup calculations. Right now, I'm feeling a bit lost and not sure where to start. I found the OECD-NEA Burnup Credit Calculational Criticality Benchmark (Phase I-B) and was wondering if anyone has worked through this specific benchmark using MCNP6? If so, would you be willing to share your MCNP input file for it? Seeing an actual working example would be incredibly helpful for my learning. I'd be really...

Similar threads

Back
Top