# Import statements
# ============================================================================
import numpy as np
import matplotlib.pyplot as plt
# File written for class
from Utilities.Utilities import Nuclide
from Utilities.Utilities import Background_Cross_Section
from Removal_Matrix import Removal_Matrix
from Scatter_Matrix import Scatter_Matrix
from Fission_Matrix import Fission_Matrix
# ============================================================================
# Load Interpolated Interpreted Plotted Data Files and Doppler-Broadened Data
# ============================================================================
H1 = Nuclide('H1', 1, 1.008, [1, 1, 0])
H1.Load_Doppler_Data([600, 900, 1200])
O16 = Nuclide('O16', 16, 15.995, [1, 1, 0])
O16.Load_Doppler_Data([600, 900, 1200])
U235 = Nuclide('U235', 235, 235.044, [1, 1, 1])
U235.Load_Doppler_Data([600, 900, 1200])
U238 = Nuclide('U238', 238, 238.051, [1, 1, 1])
U238.Load_Doppler_Data([600, 900, 1200])
# ============================================================================
# Testing
# ============================================================================
if True:
print('Saving Broadened Data to .txt files\n')
for j in range(3):
data = np.vstack(np.transpose([NEM[j], NXS[j]]))
if Nuclide.B[j]:
file_name = 'Data/Doppler/' + Nuclide.N + '_' + \
data_type[j] + '_' + str(i) + '.txt'
np.savetxt(file_name, data, delimiter=',')
print('Finished Broadening ' + Nuclide.N + ' Data\n')
# ============================================================================
# Load and Doppler-Broaden Data
# ============================================================================
# Note, the data we are performing Doppler Broadening on are the interpreted
# plot data from the BNL website, which are at 300 K.
if __name__ == "__main__":
# Temperature values to Doppler-Broaden to
Temps = [600, 900, 1200]
# Define the objects for containing data for H1, O16, U_235, U_238
H1 = Nuclide('H1', 1, 1.008, [1, 1, 0])
O16 = Nuclide('O16', 16, 15.995, [1, 1, 0])
U235 = Nuclide('U235', 235, 235.044, [1, 1, 1])
U238 = Nuclide('U238', 238, 238.051, [1, 1, 1])
# Doppler-Broaden Data
Make_Doppler_Data(H1, Temps)
Make_Doppler_Data(O16, Temps)
Make_Doppler_Data(U235, Temps)
Make_Doppler_Data(U238, Temps)
color='b',
label=r'$\sigma_{t}(E), T=1200\degree K$')
line3 = Mlines.Line2D([], [],
color='r',
label=r'$\sigma_{tg},\ T=300\degree K$')
line4 = Mlines.Line2D([], [],
color='g',
label=r'$\sigma_{tg},\ T=1200\degree K$')
ax1.legend(handles=[line1, line2, line3, line4])
return
# ============================================================================
# Loading data, and defining parameters
# ============================================================================
H1 = Nuclide('H1', 1, 1.008, [1, 1, 0])
H1.Load_Doppler_Data([600, 900, 1200])
# O16 = Nuclide('O16', 16, 15.995, [1, 1, 0])
# O16.Load_Doppler_Data([600, 900, 1200])
U235 = Nuclide('U235', 235, 235.044, [1, 1, 1])
U235.Load_Doppler_Data([600, 900, 1200])
# U238 = Nuclide('U238', 238, 238.051, [1, 1, 1])
# U238.Load_Doppler_Data([600, 900, 1200])
Casmo_16 = np.array([
1.00e1, 8.21e-1, 5.53e-3, 4.00e-6, 1.30e-6, 1.15e-6, 1.097e-6, 1.02e-6,
9.71e-7, 8.50e-7, 6.25e-7, 3.50e-7, 2.80e-7, 1.40e-7, 5.80e-8, 3.00e-8,
1.00e-11