From bc7226a97b4f36f29619ef7bc0779d049cf0d055 Mon Sep 17 00:00:00 2001 From: anu1217 Date: Wed, 31 Jul 2024 15:22:17 -0500 Subject: [PATCH 01/11] Comparison of ALARA and OpenMC SS Results --- ALARA_OpenMC_Comparison.py | 164 +++++++++++++++++++++++++++++++++++++ 1 file changed, 164 insertions(+) create mode 100644 ALARA_OpenMC_Comparison.py diff --git a/ALARA_OpenMC_Comparison.py b/ALARA_OpenMC_Comparison.py new file mode 100644 index 0000000..055bec4 --- /dev/null +++ b/ALARA_OpenMC_Comparison.py @@ -0,0 +1,164 @@ +# -*- coding: utf-8 -*- +""" +Created on Wed Jul 10 12:54:10 2024 + +@author: Anupama Rajendra +""" + +import matplotlib.pyplot as plt +import csv + +#Initializing arrays to store data from OpenMC results .csv file: +OpenMC_Nuclides = [] +OpenMC_Day = [] +OpenMC_Month = [] +OpenMC_Shutdown = [] + +with open('Densities_CSV.csv', 'r') as OpenMC_Data: + OpenMC_csv_reader = csv.reader(OpenMC_Data) + for OpenMC_line in OpenMC_csv_reader: + OpenMC_Nuclides.append(OpenMC_line[0]) + OpenMC_Day.append(OpenMC_line[2]) + OpenMC_Month.append(OpenMC_line[3]) + OpenMC_Shutdown.append(OpenMC_line[4]) + +# Reading content from ALARA Output file: +with open('ALARA_Data.txt', 'r') as Data_File: + ALARA_Data = Data_File.readlines() + +#Identifying the part of the ALARA output file that contains the relevant density data: +ALARA_FileBounds = [] +for Density_Index, Density_Content in enumerate(ALARA_Data): + if '==' in Density_Content: + ALARA_FileBounds.append(Density_Index) + #If '==' is found twice, end the loop + if len(ALARA_FileBounds) == 2: + break + +#Process the data in between the two lines that begin with '==' +ALARA_Nuclide_Data = ALARA_Data[ALARA_FileBounds[0] + 1:ALARA_FileBounds[1]] + +# Stable W nuclides present at beginning of operation: +Stable_Nuc = ['w-180', 'w-182', 'w-183', 'w-184', 'w-186'] + +#Initializing arrays to store data from ALARA and data common to both ALARA & OpenMC +ALARA_Nuc_no_W = [] +ALARA_no_W_Day = [] +ALARA_no_W_Month = [] +ALARA_no_W_Shutdown = [] +ALARA_List = [ALARA_no_W_Day, ALARA_no_W_Month, ALARA_no_W_Shutdown] + +Common_Nuclides = [] +Common_Day_Diff = [] +Common_Month_Diff = [] +Common_Shutdown_Diff = [] + +#Nuclide data found in ALARA only +ALARA_Nuclides_Only = [] +ALARA_no_W_Day_Only = [] +ALARA_no_W_Month_Only = [] +ALARA_no_W_Shutdown_Only = [] + +for ALARA_Filtered_Lines in ALARA_Nuclide_Data: + #For any nuclide that is not one of the stable W nuclides: + if not any(ALARA_Filtered_Lines.startswith(Nuclide) for Nuclide in Stable_Nuc): + #Formatting change to make OpenMC and ALARA nuclide formats match up + ALARA_Filtered_Lines = ALARA_Filtered_Lines.replace('-','').capitalize() + #Storing all density information from ALARA output + + #The first part of each line is the nuclide + Nuc_Each_Line = (ALARA_Filtered_Lines.strip().split()[0]) + ALARA_Nuc_no_W.append(Nuc_Each_Line) + ALARA_no_W_Day.append(ALARA_Filtered_Lines.strip().split()[1]) + ALARA_no_W_Month.append(ALARA_Filtered_Lines.strip().split()[2]) + ALARA_no_W_Shutdown.append(ALARA_Filtered_Lines.strip().split()[3]) + #Identifying nuclides from ALARA that are also found from OpenMC data: + if Nuc_Each_Line in OpenMC_Nuclides: + Common_Nuclides.append(Nuc_Each_Line) + else: + #Add to lists that contain nuclides/densities only found in ALARA output + ALARA_Nuclides_Only.append(Nuc_Each_Line) + ALARA_no_W_Day_Only.append(ALARA_Filtered_Lines.split()[1]) + ALARA_no_W_Month_Only.append(ALARA_Filtered_Lines.split()[2]) + ALARA_no_W_Shutdown_Only.append(ALARA_Filtered_Lines.split()[3]) + +# Times after shutdown [s] +time_steps = [0, 86400, 2.6864e+06] + +# Plot data for each isotope found in ALARA +for i, isotope in enumerate(ALARA_Nuc_no_W): + densities = [ALARA_no_W_Day[i], ALARA_no_W_Month[i], ALARA_no_W_Shutdown[i]] + plt.plot(time_steps, densities, marker='o', label=isotope) + +plt.xlabel('Time (seconds)') +plt.ylabel('Number Density (atoms/cm^3)') +plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') +plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) +plt.subplots_adjust(right=0.7) +plt.xscale('log') +plt.yscale('log') +plt.grid(True) +plt.savefig('Nuclide_density_ALARA') +plt.show() + +#Initializing arrays to store ratios of nuclide densities from both codes +Ratio_Day = [] +Ratio_Month = [] +Ratio_Shutdown = [] + +#print("All nuclides found in OpenMC", OpenMC_Nuclides) +print("Nuclides common to ALARA and OpenMC", Common_Nuclides) +#print("Nuclides found only in ALARA:", ALARA_Nuclides_Only) + +#Iterating over all nuclides common to both sets of results +for Common_Name in Common_Nuclides: + #Identifying the location of each common nuclide in the OpenMC nuclide list + Index_OpenMC = OpenMC_Nuclides.index(Common_Name) + Density_OpenMC_Day = float(OpenMC_Day[Index_OpenMC]) + Density_OpenMC_Month = float(OpenMC_Month[Index_OpenMC]) + Density_OpenMC_Shutdown = float(OpenMC_Shutdown[Index_OpenMC]) + + #Identifying the location of each common nuclide in the ALARA nuclide list + Index_ALARA = ALARA_Nuc_no_W.index(Common_Name) + Density_ALARA_Day = float(ALARA_no_W_Day[Index_ALARA]) + Density_ALARA_Month = float(ALARA_no_W_Month[Index_ALARA]) + Density_ALARA_Shutdown = float(ALARA_no_W_Shutdown[Index_ALARA]) + + #Ratios between the values from ALARA and OpenMC: + Ratio_Day.append(Density_OpenMC_Day / Density_ALARA_Day) + Ratio_Month.append(Density_OpenMC_Month / Density_ALARA_Month) + Ratio_Shutdown.append(Density_OpenMC_Shutdown / Density_ALARA_Shutdown) + + #Absolute values of the differences between ALARA and OpenMC + Common_Day_Diff.append(abs(Density_OpenMC_Day - Density_ALARA_Day)) + Common_Month_Diff.append(abs(Density_OpenMC_Month - Density_ALARA_Month)) + Common_Shutdown_Diff.append(abs(Density_OpenMC_Shutdown - Density_ALARA_Shutdown)) + +for j, com_iso in enumerate(Common_Nuclides): + common_diff = [Common_Day_Diff[j], Common_Month_Diff[j], Common_Shutdown_Diff[j]] + common_ratios = [Ratio_Day[j], Ratio_Month[j], Ratio_Shutdown[j]] + #plt.plot(time_steps, common_diff, marker='o', label=com_iso) + plt.plot(time_steps, common_ratios, marker = 'o', label = com_iso) + +# plt.xlabel('Time [seconds]') +# plt.xlim(1, 1E+7) +# #plt.ylim(1E+12, 1E+15) +# plt.ylabel('Number Density Difference [atoms/cm^3]') +# plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') +# plt.subplots_adjust(right=0.7) +# plt.xscale('log') +# plt.yscale('log') +# plt.grid(True) +#plt.savefig('Nuclide_density_diff') +# plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) +# plt.show() + +plt.xlabel('Time [seconds]') +plt.ylabel('Ratio Between OpenMC and ALARA') +plt.title('Number Density Ratio vs. Time After Shutdown [s]') +plt.subplots_adjust(right=0.7) +plt.legend(loc='upper left', bbox_to_anchor=(1.15, 1.1)) +plt.grid(True) +plt.savefig('Nuclide_density_ratio') +plt.show() +plt.close() \ No newline at end of file From 06b4ad1642b3366507a2a89c01bf2710954f51c3 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 15:28:30 -0500 Subject: [PATCH 02/11] Rename ALARA_OpenMC_Comparison.py to Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py --- .../Comparisons/ALARA_OpenMC_Comparison.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) rename ALARA_OpenMC_Comparison.py => Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py (99%) diff --git a/ALARA_OpenMC_Comparison.py b/Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py similarity index 99% rename from ALARA_OpenMC_Comparison.py rename to Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py index 055bec4..5dfdba0 100644 --- a/ALARA_OpenMC_Comparison.py +++ b/Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py @@ -161,4 +161,4 @@ plt.grid(True) plt.savefig('Nuclide_density_ratio') plt.show() -plt.close() \ No newline at end of file +plt.close() From 7e905c52a7e37bf717023699e1382af26c6c8c36 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 15:31:29 -0500 Subject: [PATCH 03/11] Delete Spherical_Shell/Comparisons directory --- .../Comparisons/ALARA_OpenMC_Comparison.py | 164 ------------------ 1 file changed, 164 deletions(-) delete mode 100644 Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py diff --git a/Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py b/Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py deleted file mode 100644 index 5dfdba0..0000000 --- a/Spherical_Shell/Comparisons/ALARA_OpenMC_Comparison.py +++ /dev/null @@ -1,164 +0,0 @@ -# -*- coding: utf-8 -*- -""" -Created on Wed Jul 10 12:54:10 2024 - -@author: Anupama Rajendra -""" - -import matplotlib.pyplot as plt -import csv - -#Initializing arrays to store data from OpenMC results .csv file: -OpenMC_Nuclides = [] -OpenMC_Day = [] -OpenMC_Month = [] -OpenMC_Shutdown = [] - -with open('Densities_CSV.csv', 'r') as OpenMC_Data: - OpenMC_csv_reader = csv.reader(OpenMC_Data) - for OpenMC_line in OpenMC_csv_reader: - OpenMC_Nuclides.append(OpenMC_line[0]) - OpenMC_Day.append(OpenMC_line[2]) - OpenMC_Month.append(OpenMC_line[3]) - OpenMC_Shutdown.append(OpenMC_line[4]) - -# Reading content from ALARA Output file: -with open('ALARA_Data.txt', 'r') as Data_File: - ALARA_Data = Data_File.readlines() - -#Identifying the part of the ALARA output file that contains the relevant density data: -ALARA_FileBounds = [] -for Density_Index, Density_Content in enumerate(ALARA_Data): - if '==' in Density_Content: - ALARA_FileBounds.append(Density_Index) - #If '==' is found twice, end the loop - if len(ALARA_FileBounds) == 2: - break - -#Process the data in between the two lines that begin with '==' -ALARA_Nuclide_Data = ALARA_Data[ALARA_FileBounds[0] + 1:ALARA_FileBounds[1]] - -# Stable W nuclides present at beginning of operation: -Stable_Nuc = ['w-180', 'w-182', 'w-183', 'w-184', 'w-186'] - -#Initializing arrays to store data from ALARA and data common to both ALARA & OpenMC -ALARA_Nuc_no_W = [] -ALARA_no_W_Day = [] -ALARA_no_W_Month = [] -ALARA_no_W_Shutdown = [] -ALARA_List = [ALARA_no_W_Day, ALARA_no_W_Month, ALARA_no_W_Shutdown] - -Common_Nuclides = [] -Common_Day_Diff = [] -Common_Month_Diff = [] -Common_Shutdown_Diff = [] - -#Nuclide data found in ALARA only -ALARA_Nuclides_Only = [] -ALARA_no_W_Day_Only = [] -ALARA_no_W_Month_Only = [] -ALARA_no_W_Shutdown_Only = [] - -for ALARA_Filtered_Lines in ALARA_Nuclide_Data: - #For any nuclide that is not one of the stable W nuclides: - if not any(ALARA_Filtered_Lines.startswith(Nuclide) for Nuclide in Stable_Nuc): - #Formatting change to make OpenMC and ALARA nuclide formats match up - ALARA_Filtered_Lines = ALARA_Filtered_Lines.replace('-','').capitalize() - #Storing all density information from ALARA output - - #The first part of each line is the nuclide - Nuc_Each_Line = (ALARA_Filtered_Lines.strip().split()[0]) - ALARA_Nuc_no_W.append(Nuc_Each_Line) - ALARA_no_W_Day.append(ALARA_Filtered_Lines.strip().split()[1]) - ALARA_no_W_Month.append(ALARA_Filtered_Lines.strip().split()[2]) - ALARA_no_W_Shutdown.append(ALARA_Filtered_Lines.strip().split()[3]) - #Identifying nuclides from ALARA that are also found from OpenMC data: - if Nuc_Each_Line in OpenMC_Nuclides: - Common_Nuclides.append(Nuc_Each_Line) - else: - #Add to lists that contain nuclides/densities only found in ALARA output - ALARA_Nuclides_Only.append(Nuc_Each_Line) - ALARA_no_W_Day_Only.append(ALARA_Filtered_Lines.split()[1]) - ALARA_no_W_Month_Only.append(ALARA_Filtered_Lines.split()[2]) - ALARA_no_W_Shutdown_Only.append(ALARA_Filtered_Lines.split()[3]) - -# Times after shutdown [s] -time_steps = [0, 86400, 2.6864e+06] - -# Plot data for each isotope found in ALARA -for i, isotope in enumerate(ALARA_Nuc_no_W): - densities = [ALARA_no_W_Day[i], ALARA_no_W_Month[i], ALARA_no_W_Shutdown[i]] - plt.plot(time_steps, densities, marker='o', label=isotope) - -plt.xlabel('Time (seconds)') -plt.ylabel('Number Density (atoms/cm^3)') -plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') -plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) -plt.subplots_adjust(right=0.7) -plt.xscale('log') -plt.yscale('log') -plt.grid(True) -plt.savefig('Nuclide_density_ALARA') -plt.show() - -#Initializing arrays to store ratios of nuclide densities from both codes -Ratio_Day = [] -Ratio_Month = [] -Ratio_Shutdown = [] - -#print("All nuclides found in OpenMC", OpenMC_Nuclides) -print("Nuclides common to ALARA and OpenMC", Common_Nuclides) -#print("Nuclides found only in ALARA:", ALARA_Nuclides_Only) - -#Iterating over all nuclides common to both sets of results -for Common_Name in Common_Nuclides: - #Identifying the location of each common nuclide in the OpenMC nuclide list - Index_OpenMC = OpenMC_Nuclides.index(Common_Name) - Density_OpenMC_Day = float(OpenMC_Day[Index_OpenMC]) - Density_OpenMC_Month = float(OpenMC_Month[Index_OpenMC]) - Density_OpenMC_Shutdown = float(OpenMC_Shutdown[Index_OpenMC]) - - #Identifying the location of each common nuclide in the ALARA nuclide list - Index_ALARA = ALARA_Nuc_no_W.index(Common_Name) - Density_ALARA_Day = float(ALARA_no_W_Day[Index_ALARA]) - Density_ALARA_Month = float(ALARA_no_W_Month[Index_ALARA]) - Density_ALARA_Shutdown = float(ALARA_no_W_Shutdown[Index_ALARA]) - - #Ratios between the values from ALARA and OpenMC: - Ratio_Day.append(Density_OpenMC_Day / Density_ALARA_Day) - Ratio_Month.append(Density_OpenMC_Month / Density_ALARA_Month) - Ratio_Shutdown.append(Density_OpenMC_Shutdown / Density_ALARA_Shutdown) - - #Absolute values of the differences between ALARA and OpenMC - Common_Day_Diff.append(abs(Density_OpenMC_Day - Density_ALARA_Day)) - Common_Month_Diff.append(abs(Density_OpenMC_Month - Density_ALARA_Month)) - Common_Shutdown_Diff.append(abs(Density_OpenMC_Shutdown - Density_ALARA_Shutdown)) - -for j, com_iso in enumerate(Common_Nuclides): - common_diff = [Common_Day_Diff[j], Common_Month_Diff[j], Common_Shutdown_Diff[j]] - common_ratios = [Ratio_Day[j], Ratio_Month[j], Ratio_Shutdown[j]] - #plt.plot(time_steps, common_diff, marker='o', label=com_iso) - plt.plot(time_steps, common_ratios, marker = 'o', label = com_iso) - -# plt.xlabel('Time [seconds]') -# plt.xlim(1, 1E+7) -# #plt.ylim(1E+12, 1E+15) -# plt.ylabel('Number Density Difference [atoms/cm^3]') -# plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') -# plt.subplots_adjust(right=0.7) -# plt.xscale('log') -# plt.yscale('log') -# plt.grid(True) -#plt.savefig('Nuclide_density_diff') -# plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) -# plt.show() - -plt.xlabel('Time [seconds]') -plt.ylabel('Ratio Between OpenMC and ALARA') -plt.title('Number Density Ratio vs. Time After Shutdown [s]') -plt.subplots_adjust(right=0.7) -plt.legend(loc='upper left', bbox_to_anchor=(1.15, 1.1)) -plt.grid(True) -plt.savefig('Nuclide_density_ratio') -plt.show() -plt.close() From 69bc34429ed8e2475ca158dd3798d483e91af341 Mon Sep 17 00:00:00 2001 From: anu1217 Date: Wed, 31 Jul 2024 15:32:53 -0500 Subject: [PATCH 04/11] Comparison between ALARA and OpenMC results --- ALARA_OpenMC_Comparison.py | 164 +++++++++++++++++++++++++++++++++++++ 1 file changed, 164 insertions(+) create mode 100644 ALARA_OpenMC_Comparison.py diff --git a/ALARA_OpenMC_Comparison.py b/ALARA_OpenMC_Comparison.py new file mode 100644 index 0000000..055bec4 --- /dev/null +++ b/ALARA_OpenMC_Comparison.py @@ -0,0 +1,164 @@ +# -*- coding: utf-8 -*- +""" +Created on Wed Jul 10 12:54:10 2024 + +@author: Anupama Rajendra +""" + +import matplotlib.pyplot as plt +import csv + +#Initializing arrays to store data from OpenMC results .csv file: +OpenMC_Nuclides = [] +OpenMC_Day = [] +OpenMC_Month = [] +OpenMC_Shutdown = [] + +with open('Densities_CSV.csv', 'r') as OpenMC_Data: + OpenMC_csv_reader = csv.reader(OpenMC_Data) + for OpenMC_line in OpenMC_csv_reader: + OpenMC_Nuclides.append(OpenMC_line[0]) + OpenMC_Day.append(OpenMC_line[2]) + OpenMC_Month.append(OpenMC_line[3]) + OpenMC_Shutdown.append(OpenMC_line[4]) + +# Reading content from ALARA Output file: +with open('ALARA_Data.txt', 'r') as Data_File: + ALARA_Data = Data_File.readlines() + +#Identifying the part of the ALARA output file that contains the relevant density data: +ALARA_FileBounds = [] +for Density_Index, Density_Content in enumerate(ALARA_Data): + if '==' in Density_Content: + ALARA_FileBounds.append(Density_Index) + #If '==' is found twice, end the loop + if len(ALARA_FileBounds) == 2: + break + +#Process the data in between the two lines that begin with '==' +ALARA_Nuclide_Data = ALARA_Data[ALARA_FileBounds[0] + 1:ALARA_FileBounds[1]] + +# Stable W nuclides present at beginning of operation: +Stable_Nuc = ['w-180', 'w-182', 'w-183', 'w-184', 'w-186'] + +#Initializing arrays to store data from ALARA and data common to both ALARA & OpenMC +ALARA_Nuc_no_W = [] +ALARA_no_W_Day = [] +ALARA_no_W_Month = [] +ALARA_no_W_Shutdown = [] +ALARA_List = [ALARA_no_W_Day, ALARA_no_W_Month, ALARA_no_W_Shutdown] + +Common_Nuclides = [] +Common_Day_Diff = [] +Common_Month_Diff = [] +Common_Shutdown_Diff = [] + +#Nuclide data found in ALARA only +ALARA_Nuclides_Only = [] +ALARA_no_W_Day_Only = [] +ALARA_no_W_Month_Only = [] +ALARA_no_W_Shutdown_Only = [] + +for ALARA_Filtered_Lines in ALARA_Nuclide_Data: + #For any nuclide that is not one of the stable W nuclides: + if not any(ALARA_Filtered_Lines.startswith(Nuclide) for Nuclide in Stable_Nuc): + #Formatting change to make OpenMC and ALARA nuclide formats match up + ALARA_Filtered_Lines = ALARA_Filtered_Lines.replace('-','').capitalize() + #Storing all density information from ALARA output + + #The first part of each line is the nuclide + Nuc_Each_Line = (ALARA_Filtered_Lines.strip().split()[0]) + ALARA_Nuc_no_W.append(Nuc_Each_Line) + ALARA_no_W_Day.append(ALARA_Filtered_Lines.strip().split()[1]) + ALARA_no_W_Month.append(ALARA_Filtered_Lines.strip().split()[2]) + ALARA_no_W_Shutdown.append(ALARA_Filtered_Lines.strip().split()[3]) + #Identifying nuclides from ALARA that are also found from OpenMC data: + if Nuc_Each_Line in OpenMC_Nuclides: + Common_Nuclides.append(Nuc_Each_Line) + else: + #Add to lists that contain nuclides/densities only found in ALARA output + ALARA_Nuclides_Only.append(Nuc_Each_Line) + ALARA_no_W_Day_Only.append(ALARA_Filtered_Lines.split()[1]) + ALARA_no_W_Month_Only.append(ALARA_Filtered_Lines.split()[2]) + ALARA_no_W_Shutdown_Only.append(ALARA_Filtered_Lines.split()[3]) + +# Times after shutdown [s] +time_steps = [0, 86400, 2.6864e+06] + +# Plot data for each isotope found in ALARA +for i, isotope in enumerate(ALARA_Nuc_no_W): + densities = [ALARA_no_W_Day[i], ALARA_no_W_Month[i], ALARA_no_W_Shutdown[i]] + plt.plot(time_steps, densities, marker='o', label=isotope) + +plt.xlabel('Time (seconds)') +plt.ylabel('Number Density (atoms/cm^3)') +plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') +plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) +plt.subplots_adjust(right=0.7) +plt.xscale('log') +plt.yscale('log') +plt.grid(True) +plt.savefig('Nuclide_density_ALARA') +plt.show() + +#Initializing arrays to store ratios of nuclide densities from both codes +Ratio_Day = [] +Ratio_Month = [] +Ratio_Shutdown = [] + +#print("All nuclides found in OpenMC", OpenMC_Nuclides) +print("Nuclides common to ALARA and OpenMC", Common_Nuclides) +#print("Nuclides found only in ALARA:", ALARA_Nuclides_Only) + +#Iterating over all nuclides common to both sets of results +for Common_Name in Common_Nuclides: + #Identifying the location of each common nuclide in the OpenMC nuclide list + Index_OpenMC = OpenMC_Nuclides.index(Common_Name) + Density_OpenMC_Day = float(OpenMC_Day[Index_OpenMC]) + Density_OpenMC_Month = float(OpenMC_Month[Index_OpenMC]) + Density_OpenMC_Shutdown = float(OpenMC_Shutdown[Index_OpenMC]) + + #Identifying the location of each common nuclide in the ALARA nuclide list + Index_ALARA = ALARA_Nuc_no_W.index(Common_Name) + Density_ALARA_Day = float(ALARA_no_W_Day[Index_ALARA]) + Density_ALARA_Month = float(ALARA_no_W_Month[Index_ALARA]) + Density_ALARA_Shutdown = float(ALARA_no_W_Shutdown[Index_ALARA]) + + #Ratios between the values from ALARA and OpenMC: + Ratio_Day.append(Density_OpenMC_Day / Density_ALARA_Day) + Ratio_Month.append(Density_OpenMC_Month / Density_ALARA_Month) + Ratio_Shutdown.append(Density_OpenMC_Shutdown / Density_ALARA_Shutdown) + + #Absolute values of the differences between ALARA and OpenMC + Common_Day_Diff.append(abs(Density_OpenMC_Day - Density_ALARA_Day)) + Common_Month_Diff.append(abs(Density_OpenMC_Month - Density_ALARA_Month)) + Common_Shutdown_Diff.append(abs(Density_OpenMC_Shutdown - Density_ALARA_Shutdown)) + +for j, com_iso in enumerate(Common_Nuclides): + common_diff = [Common_Day_Diff[j], Common_Month_Diff[j], Common_Shutdown_Diff[j]] + common_ratios = [Ratio_Day[j], Ratio_Month[j], Ratio_Shutdown[j]] + #plt.plot(time_steps, common_diff, marker='o', label=com_iso) + plt.plot(time_steps, common_ratios, marker = 'o', label = com_iso) + +# plt.xlabel('Time [seconds]') +# plt.xlim(1, 1E+7) +# #plt.ylim(1E+12, 1E+15) +# plt.ylabel('Number Density Difference [atoms/cm^3]') +# plt.title('Number Density [atoms/cm^3] vs. Time After Shutdown [s]') +# plt.subplots_adjust(right=0.7) +# plt.xscale('log') +# plt.yscale('log') +# plt.grid(True) +#plt.savefig('Nuclide_density_diff') +# plt.legend(loc='upper left', bbox_to_anchor=(1.05, 1)) +# plt.show() + +plt.xlabel('Time [seconds]') +plt.ylabel('Ratio Between OpenMC and ALARA') +plt.title('Number Density Ratio vs. Time After Shutdown [s]') +plt.subplots_adjust(right=0.7) +plt.legend(loc='upper left', bbox_to_anchor=(1.15, 1.1)) +plt.grid(True) +plt.savefig('Nuclide_density_ratio') +plt.show() +plt.close() \ No newline at end of file From e9115e3ac82bbe2dae293ad63c016ee05f3c8cce Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 15:33:53 -0500 Subject: [PATCH 05/11] Rename ALARA_OpenMC_Comparison.py to SphericalShell/ALARA_OpenMC_Comparison.py --- .../ALARA_OpenMC_Comparison.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) rename ALARA_OpenMC_Comparison.py => SphericalShell/ALARA_OpenMC_Comparison.py (99%) diff --git a/ALARA_OpenMC_Comparison.py b/SphericalShell/ALARA_OpenMC_Comparison.py similarity index 99% rename from ALARA_OpenMC_Comparison.py rename to SphericalShell/ALARA_OpenMC_Comparison.py index 055bec4..5dfdba0 100644 --- a/ALARA_OpenMC_Comparison.py +++ b/SphericalShell/ALARA_OpenMC_Comparison.py @@ -161,4 +161,4 @@ plt.grid(True) plt.savefig('Nuclide_density_ratio') plt.show() -plt.close() \ No newline at end of file +plt.close() From 476663c6914d82b58c27e42345bed57835eef00f Mon Sep 17 00:00:00 2001 From: anu1217 Date: Wed, 31 Jul 2024 15:39:37 -0500 Subject: [PATCH 06/11] Switching to CoupledOperator for greater precision --- SS_Activation_OpenMC.py | 12 +++++++----- 1 file changed, 7 insertions(+), 5 deletions(-) diff --git a/SS_Activation_OpenMC.py b/SS_Activation_OpenMC.py index 7e2e185..79d2aa6 100644 --- a/SS_Activation_OpenMC.py +++ b/SS_Activation_OpenMC.py @@ -85,8 +85,8 @@ W.depletable = True W.volume = 4.0/3.0 * np.pi * (R_2**3 - R_1**3) #volume of W wall material fluxes, micros = openmc.deplete.get_microxs_and_flux(model, Cells) -operator = openmc.deplete.IndependentOperator(materials, fluxes[0:1], micros[0:1],normalization_mode='source-rate') -# operator = openmc.deplete.CoupledOperator(model, normalization_mode='source-rate') +#operator = openmc.deplete.IndependentOperator(materials, fluxes[0:1], micros[0:1],normalization_mode='source-rate') +operator = openmc.deplete.CoupledOperator(model, normalization_mode='source-rate') time_steps = [3e8, 86400, 2.6e6] source_rates = [1E+18, 0, 0] integrator = openmc.deplete.PredictorIntegrator(operator=operator, timesteps=time_steps, source_rates=source_rates, timestep_units='s') @@ -154,13 +154,14 @@ plt.plot(time, num_dens[nuclide], marker='.', linestyle='solid', color=plot_color, label=nuclide) # Adding labels and title -plt.xlabel('Time after beginning of operation [s]') -plt.xlim(1, sum(time_steps) +plt.xlabel('Time after start of operation [s]') +plt.xlim(1, sum(time_steps)) +#plt.ylim(1e-09, 1e+20) #plt.gca().set_ylim(bottom=0) plt.ylabel('Nuclide density [atoms/cm^3]') plt.xscale("log") plt.yscale("log") -plt.title('Plot of number density vs time') +plt.title('Plot of number density vs time after operation') # Adding a legend plt.legend() @@ -168,3 +169,4 @@ plt.savefig('Nuclide_density_OpenMC') # Display the plot plt.show() +plt.close() \ No newline at end of file From 47d31df8248a597451609ffda6e98ce35413b4f6 Mon Sep 17 00:00:00 2001 From: anu1217 Date: Wed, 31 Jul 2024 15:59:58 -0500 Subject: [PATCH 07/11] Storing nuclides from all timesteps --- SS_Activation_OpenMC.py | 26 +++++++++++++++----------- 1 file changed, 15 insertions(+), 11 deletions(-) diff --git a/SS_Activation_OpenMC.py b/SS_Activation_OpenMC.py index 79d2aa6..f348eff 100644 --- a/SS_Activation_OpenMC.py +++ b/SS_Activation_OpenMC.py @@ -84,7 +84,7 @@ #Depletion calculation W.depletable = True W.volume = 4.0/3.0 * np.pi * (R_2**3 - R_1**3) #volume of W wall material -fluxes, micros = openmc.deplete.get_microxs_and_flux(model, Cells) +#fluxes, micros = openmc.deplete.get_microxs_and_flux(model, Cells) #operator = openmc.deplete.IndependentOperator(materials, fluxes[0:1], micros[0:1],normalization_mode='source-rate') operator = openmc.deplete.CoupledOperator(model, normalization_mode='source-rate') time_steps = [3e8, 86400, 2.6e6] @@ -127,25 +127,29 @@ # Stable W nuclides present at beginning of operation (will not be plotted) stable_nuc = ['W180', 'W182', 'W183', 'W184', 'W186'] -#Store list of nuclides from last timestep as a Materials object -materials_last = results.export_to_materials(-1) -# Storing depletion data from 1st material -mat_dep = materials_last[0] -# Obtaining the list of nuclides from the results file -nuc_last = mat_dep.get_nuclides() +materials_list=[] +#Store list of nuclides from each timestep as a Materials object +for step in range(len(time_steps)): + materials_object = results.export_to_materials(step) + # Storing depletion data from 1st material + mat_dep = materials_object[0] + # Obtaining the list of nuclides in the results file + rad_nuc = mat_dep.get_nuclides() + materials_list.append(rad_nuc) + # Removing stable W nuclides from list so that they do not appear in the plot for j in stable_nuc : - nuc_last.remove(j) -print(nuc_last) + rad_nuc.remove(j) +print(rad_nuc) colors = list(mcolors.CSS4_COLORS.keys()) num_dens= {} pair_list = {} -with open(r'Densities_CSV.csv', 'a') as density_file: +with open(r'Densities_CSV.csv', 'w') as density_file: - for nuclide in nuc_last: + for nuclide in rad_nuc: plot_color = random.choice(colors) time, num_dens[nuclide] = results.get_atoms('1', nuclide, nuc_units = 'atom/cm3') print(time, num_dens[nuclide]) From f753c3df85810efbc7ea57c6b42b1a6d99f389fc Mon Sep 17 00:00:00 2001 From: anu1217 Date: Wed, 31 Jul 2024 16:01:24 -0500 Subject: [PATCH 08/11] Changing ALARA flux to go from high-low energy --- flux_ss.txt | 206 ++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 206 insertions(+) create mode 100644 flux_ss.txt diff --git a/flux_ss.txt b/flux_ss.txt new file mode 100644 index 0000000..5b63432 --- /dev/null +++ b/flux_ss.txt @@ -0,0 +1,206 @@ + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 + 0.00000E+00 + +0 +0 +0 +0 +0 +0 +0 +5.62E-08 +1.78E-09 +2.44E-10 +1.33E-10 +1.43E-10 +2.57E-10 +2.84E-10 +2.57E-10 +2.73E-10 +3.06E-10 +2.97E-10 +3.05E-10 +2.66E-10 +2.76E-10 +2.43E-10 +2.25E-10 +2.56E-10 +5.81E-11 +1.48E-10 +2.36E-10 +2.80E-10 +2.94E-10 +3.07E-10 +3.08E-10 +3.38E-10 +3.69E-10 +8.40E-10 +1.05E-09 +1.18E-09 +7.44E-10 +7.84E-10 +8.35E-10 +8.89E-10 +1.05E-09 +1.04E-09 +8.19E-10 +1.89E-10 +1.92E-10 +3.67E-10 +7.91E-10 +1.24E-09 +1.41E-09 +1.55E-09 +1.58E-09 +1.67E-09 +1.74E-09 +1.97E-09 +2.07E-09 +2.24E-09 +2.30E-09 +2.56E-09 +2.77E-09 +2.91E-09 +3.17E-09 +7.13E-09 +3.18E-09 +4.72E-09 +4.39E-09 +4.65E-09 +4.76E-09 +4.93E-09 +5.26E-09 +5.58E-09 +5.55E-09 +5.80E-09 +6.06E-09 +6.03E-09 +6.17E-09 +6.30E-09 +1.26E-08 +1.26E-08 +6.15E-09 +6.11E-09 +1.20E-08 +1.14E-08 +1.33E-09 +5.03E-10 +1.00E-09 +2.75E-09 +5.60E-09 +1.04E-08 +5.05E-09 +4.98E-09 +4.92E-09 +4.65E-09 +4.39E-09 +4.24E-09 +4.16E-09 +4.04E-09 +3.87E-09 +3.83E-09 +3.74E-09 +3.53E-09 +3.36E-09 +3.25E-09 +3.15E-09 +3.00E-09 +7.20E-09 +6.55E-09 +2.12E-09 +1.52E-09 +3.72E-09 +2.21E-09 +5.04E-09 +1.80E-09 +2.61E-09 +2.30E-09 +2.66E-09 +1.00E-09 +1.23E-09 +4.86E-10 +3.21E-10 +3.99E-10 +1.91E-10 +1.60E-10 +4.95E-10 +7.19E-10 +1.04E-09 +7.33E-10 +1.79E-10 +2.44E-10 +3.10E-10 +1.68E-10 +9.88E-11 +4.64E-11 +2.15E-11 +1.61E-11 +2.11E-11 +6.80E-12 +4.26E-12 +1.45E-11 +1.24E-11 +1.83E-11 +1.47E-11 +6.21E-12 +4.00E-13 +0 +0 +0 +4.64E-13 +2.68E-13 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 +0 From 123b824843a34f05102665dc4fe4aabd74830a3f Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 16:10:41 -0500 Subject: [PATCH 09/11] Changing flux file format for readability --- flux_ss.txt | 207 ++++++++-------------------------------------------- 1 file changed, 31 insertions(+), 176 deletions(-) diff --git a/flux_ss.txt b/flux_ss.txt index 5b63432..3c45ca9 100644 --- a/flux_ss.txt +++ b/flux_ss.txt @@ -28,179 +28,34 @@ 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - -0 -0 -0 -0 -0 -0 -0 -5.62E-08 -1.78E-09 -2.44E-10 -1.33E-10 -1.43E-10 -2.57E-10 -2.84E-10 -2.57E-10 -2.73E-10 -3.06E-10 -2.97E-10 -3.05E-10 -2.66E-10 -2.76E-10 -2.43E-10 -2.25E-10 -2.56E-10 -5.81E-11 -1.48E-10 -2.36E-10 -2.80E-10 -2.94E-10 -3.07E-10 -3.08E-10 -3.38E-10 -3.69E-10 -8.40E-10 -1.05E-09 -1.18E-09 -7.44E-10 -7.84E-10 -8.35E-10 -8.89E-10 -1.05E-09 -1.04E-09 -8.19E-10 -1.89E-10 -1.92E-10 -3.67E-10 -7.91E-10 -1.24E-09 -1.41E-09 -1.55E-09 -1.58E-09 -1.67E-09 -1.74E-09 -1.97E-09 -2.07E-09 -2.24E-09 -2.30E-09 -2.56E-09 -2.77E-09 -2.91E-09 -3.17E-09 -7.13E-09 -3.18E-09 -4.72E-09 -4.39E-09 -4.65E-09 -4.76E-09 -4.93E-09 -5.26E-09 -5.58E-09 -5.55E-09 -5.80E-09 -6.06E-09 -6.03E-09 -6.17E-09 -6.30E-09 -1.26E-08 -1.26E-08 -6.15E-09 -6.11E-09 -1.20E-08 -1.14E-08 -1.33E-09 -5.03E-10 -1.00E-09 -2.75E-09 -5.60E-09 -1.04E-08 -5.05E-09 -4.98E-09 -4.92E-09 -4.65E-09 -4.39E-09 -4.24E-09 -4.16E-09 -4.04E-09 -3.87E-09 -3.83E-09 -3.74E-09 -3.53E-09 -3.36E-09 -3.25E-09 -3.15E-09 -3.00E-09 -7.20E-09 -6.55E-09 -2.12E-09 -1.52E-09 -3.72E-09 -2.21E-09 -5.04E-09 -1.80E-09 -2.61E-09 -2.30E-09 -2.66E-09 -1.00E-09 -1.23E-09 -4.86E-10 -3.21E-10 -3.99E-10 -1.91E-10 -1.60E-10 -4.95E-10 -7.19E-10 -1.04E-09 -7.33E-10 -1.79E-10 -2.44E-10 -3.10E-10 -1.68E-10 -9.88E-11 -4.64E-11 -2.15E-11 -1.61E-11 -2.11E-11 -6.80E-12 -4.26E-12 -1.45E-11 -1.24E-11 -1.83E-11 -1.47E-11 -6.21E-12 -4.00E-13 -0 -0 -0 -4.64E-13 -2.68E-13 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 + +0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 +0.00000E+00 5.61859E-08 1.77985E-09 2.44430E-10 1.33249E-10 1.43070E-10 +2.57039E-10 2.84265E-10 2.57053E-10 2.73441E-10 3.05985E-10 2.97163E-10 +3.05423E-10 2.66059E-10 2.75828E-10 2.42649E-10 2.25309E-10 2.55894E-10 +5.80846E-11 1.47979E-10 2.36268E-10 2.80374E-10 2.94313E-10 3.06798E-10 +3.07944E-10 3.38472E-10 3.68752E-10 8.39627E-10 1.05107E-09 1.18463E-09 +7.43989E-10 7.84362E-10 8.34827E-10 8.88988E-10 1.05383E-09 1.03639E-09 +8.19308E-10 1.89473E-10 1.91778E-10 3.67294E-10 7.90709E-10 1.23886E-09 +1.40678E-09 1.54582E-09 1.58122E-09 1.66747E-09 1.74264E-09 1.97076E-09 +2.07143E-09 2.23929E-09 2.29999E-09 2.56112E-09 2.77486E-09 2.90754E-09 +3.17344E-09 7.13059E-09 3.18010E-09 4.71682E-09 4.39266E-09 4.64557E-09 +4.76033E-09 4.93324E-09 5.26116E-09 5.58171E-09 5.54849E-09 5.79688E-09 +6.06466E-09 6.03362E-09 6.17182E-09 6.30089E-09 1.26341E-08 1.26458E-08 +6.14676E-09 6.11037E-09 1.19974E-08 1.14091E-08 1.32644E-09 5.03253E-10 +1.00192E-09 2.75104E-09 5.60392E-09 1.03880E-08 5.04515E-09 4.97501E-09 +4.91826E-09 4.65095E-09 4.38533E-09 4.24375E-09 4.15808E-09 4.04411E-09 +3.87115E-09 3.82862E-09 3.74139E-09 3.52514E-09 3.35997E-09 3.25314E-09 +3.14582E-09 3.00398E-09 7.19754E-09 6.54630E-09 2.11976E-09 1.52330E-09 +3.72074E-09 2.21419E-09 5.04412E-09 1.80260E-09 2.61338E-09 2.30091E-09 +2.66260E-09 1.00339E-09 1.23306E-09 4.85876E-10 3.21014E-10 3.98677E-10 +1.90782E-10 1.60292E-10 4.95211E-10 7.19330E-10 1.03528E-09 7.32532E-10 +1.79365E-10 2.43513E-10 3.10382E-10 1.68471E-10 9.87651E-11 4.64443E-11 +2.14687E-11 1.61128E-11 2.10588E-11 6.80051E-12 4.26141E-12 1.44660E-11 +1.24312E-11 1.82837E-11 1.47184E-11 6.20704E-12 4.00378E-13 0.00000E+00 +0.00000E+00 0.00000E+00 4.64253E-13 2.67558E-13 0.00000E+00 0.00000E+00 +0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 +0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 +0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 +0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 +0.00000E+00 From d66cd1b9909bb27a5bfc5883b499f01be5ffe8b9 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 16:12:54 -0500 Subject: [PATCH 10/11] Rename flux_ss.txt to ALARA_flux_ss.txt --- flux_ss.txt => ALARA_flux_ss.txt | 0 1 file changed, 0 insertions(+), 0 deletions(-) rename flux_ss.txt => ALARA_flux_ss.txt (100%) diff --git a/flux_ss.txt b/ALARA_flux_ss.txt similarity index 100% rename from flux_ss.txt rename to ALARA_flux_ss.txt From c1a46ec461fe6775119d281514ba7ac88a7fa5b0 Mon Sep 17 00:00:00 2001 From: Anupama Rajendra <113371601+anu1217@users.noreply.github.com> Date: Wed, 31 Jul 2024 16:19:39 -0500 Subject: [PATCH 11/11] Delete ALARA_flux_ss.txt --- ALARA_flux_ss.txt | 61 ----------------------------------------------- 1 file changed, 61 deletions(-) delete mode 100644 ALARA_flux_ss.txt diff --git a/ALARA_flux_ss.txt b/ALARA_flux_ss.txt deleted file mode 100644 index 3c45ca9..0000000 --- a/ALARA_flux_ss.txt +++ /dev/null @@ -1,61 +0,0 @@ - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 - 0.00000E+00 - -0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -0.00000E+00 5.61859E-08 1.77985E-09 2.44430E-10 1.33249E-10 1.43070E-10 -2.57039E-10 2.84265E-10 2.57053E-10 2.73441E-10 3.05985E-10 2.97163E-10 -3.05423E-10 2.66059E-10 2.75828E-10 2.42649E-10 2.25309E-10 2.55894E-10 -5.80846E-11 1.47979E-10 2.36268E-10 2.80374E-10 2.94313E-10 3.06798E-10 -3.07944E-10 3.38472E-10 3.68752E-10 8.39627E-10 1.05107E-09 1.18463E-09 -7.43989E-10 7.84362E-10 8.34827E-10 8.88988E-10 1.05383E-09 1.03639E-09 -8.19308E-10 1.89473E-10 1.91778E-10 3.67294E-10 7.90709E-10 1.23886E-09 -1.40678E-09 1.54582E-09 1.58122E-09 1.66747E-09 1.74264E-09 1.97076E-09 -2.07143E-09 2.23929E-09 2.29999E-09 2.56112E-09 2.77486E-09 2.90754E-09 -3.17344E-09 7.13059E-09 3.18010E-09 4.71682E-09 4.39266E-09 4.64557E-09 -4.76033E-09 4.93324E-09 5.26116E-09 5.58171E-09 5.54849E-09 5.79688E-09 -6.06466E-09 6.03362E-09 6.17182E-09 6.30089E-09 1.26341E-08 1.26458E-08 -6.14676E-09 6.11037E-09 1.19974E-08 1.14091E-08 1.32644E-09 5.03253E-10 -1.00192E-09 2.75104E-09 5.60392E-09 1.03880E-08 5.04515E-09 4.97501E-09 -4.91826E-09 4.65095E-09 4.38533E-09 4.24375E-09 4.15808E-09 4.04411E-09 -3.87115E-09 3.82862E-09 3.74139E-09 3.52514E-09 3.35997E-09 3.25314E-09 -3.14582E-09 3.00398E-09 7.19754E-09 6.54630E-09 2.11976E-09 1.52330E-09 -3.72074E-09 2.21419E-09 5.04412E-09 1.80260E-09 2.61338E-09 2.30091E-09 -2.66260E-09 1.00339E-09 1.23306E-09 4.85876E-10 3.21014E-10 3.98677E-10 -1.90782E-10 1.60292E-10 4.95211E-10 7.19330E-10 1.03528E-09 7.32532E-10 -1.79365E-10 2.43513E-10 3.10382E-10 1.68471E-10 9.87651E-11 4.64443E-11 -2.14687E-11 1.61128E-11 2.10588E-11 6.80051E-12 4.26141E-12 1.44660E-11 -1.24312E-11 1.82837E-11 1.47184E-11 6.20704E-12 4.00378E-13 0.00000E+00 -0.00000E+00 0.00000E+00 4.64253E-13 2.67558E-13 0.00000E+00 0.00000E+00 -0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -0.00000E+00