Chapter9-Radiation Doses and Hazard Assessment

Ex1-pg241

In [1]:
##Chapter 9, Example 9.1, Page 241
import math
## Iron kerma and absorbed dose rates
Sp = 10**14
r = 100.
mew = 0.03031
mtr = 0.02112 ## mew/pro
men = 0.01983 ## mew/pro
p0 = 10**-6*Sp*math.exp(-mew*r)/(4.*math.pi*r**2)
K0 = 1.602*10**-10*mtr*p0
D0 = 1.602*10**-10*men*p0
print'%s %.2e %s'%("p0 = ",p0," cm^-2s^-1\n")
print'%s %.2e %s'%(" K0 = ",K0," Gy/s\n")
print'%s %.2e %s'%(" D0 = ",D0," Gy/s\n")
## Answers may vary due to round off error

p0 =  3.84e+01  cm^-2s^-1

 K0 =  1.30e-10  Gy/s

 D0 =  1.22e-10  Gy/s

Ex2-pg242

In [2]:
##Chapter 9, Example 9.2, Page 242
import math
## kerma rate 
fsMs = (0.6022/18.)*((2.*12.8*0.5)+(3.5*0.1107))
K = 1.602*10**-10*fsMs*10**10*0.1
print'%s %.2f %s'%("fsUs/p = ",fsMs," cm^2/g\n")
print'%s %.2f %s'%(" K = ",K," Gy/s\n")
## Answers may vary due to round off error

fsUs/p =  0.44  cm^2/g

 K =  0.07  Gy/s

Ex3-pg245

In [3]:
##Chapter 9, Example 9.3, Page 245
import math
##Find fluence and H
Sp = 10**9
dt = 600.
r = 1500.
E = 0.03103
phi = Sp*dt/(4.*math.pi*r**2)
H = 1.602*10**-10*E*phi
print'%s %.2f %s'%("fluence = ",phi," cm^2\n")
print'%s %.2f %s'%(" H = ",H*10**8," microSv\n")
## Answer may vary due to round off error

fluence =  21220.66  cm^2

 H =  10.55  microSv