Adding lead material

Changing TRIGA fuel radius (was using clad radius)

src/indicators.py

@@ -205,6 +205,7 @@ def do_all(inputDict):
             materials.append(mat.get_triga_moderator_material())
             materials.append(mat.get_triga_air_material())
             materials.append(mat.get_triga_grid_plate_material())
+            materials.append(mat.get_triga_lead_material())
     elif materialOpt == 'trigamore':
         materials.append(mat.get_triga_fuel_material())
         materials.append(mat.get_triga_clad_material())
@@ -216,6 +217,7 @@ def do_all(inputDict):
             materials.append(mat.get_triga_moderator_material())
             materials.append(mat.get_triga_air_material())
             materials.append(mat.get_triga_grid_plate_material())
+            materials.append(mat.get_triga_lead_material())
     elif materialOpt == 'deb':
         materials.append(mat.get_bruss_enriched_rod_fuel_material())
     else:

src/indicators_clustering.py

@@ -403,6 +403,8 @@ def auto_apportion(observations, numElementsRRR, groupBdrs, coarseBdrs, apportio
     metric = np.zeros(numCoarseGroups)
     numFineGroupsPerCoarseGroup = np.zeros(numCoarseGroups, dtype=np.int)
     timeBirch = 0.0
+    if numCoarseGroups == 1:
+        return np.array([numElementsRRR])
     for coarseGroup in range(len(coarseBdrs)-1):
         # Only look at the fine groups within the current coarse group
         strt = np.argmin(np.abs(groupBdrs - coarseBdrs[coarseGroup+1]))

src/materials_materials.py

@@ -65,6 +65,7 @@ def get_materials_name2function_dict():
         'tGRAPHITE': get_triga_graphite_material,
         'tBORATEDGRAPHITE': get_triga_borated_graphite_material,
         'tB4C': get_triga_b4c_material,
+        'tLEAD': get_triga_lead_material,
         #'tAIRTUBE': get_triga_air_tube_material,
         #'tIRRADIATIONTUBE': get_triga_irradiation_tube_material,
         # Iron (for time-dependent dissertation problem)
@@ -664,7 +665,7 @@ def get_triga_fuel_material():
     shortName = 'tFUEL'
     longName = 'U-ZrH fuel'
     atomDensity = 8.71115E-2
-    fuelRadius = 1.7920 #cm
+    fuelRadius = 1.7411 #cm
     temperature = 296. #K
     temperatureIndex = 0 # X in .9Xc
     thermalOpt = 'zrh'
@@ -936,6 +937,29 @@ def get_triga_air_material():
         temperatureIndex=temperatureIndex)
     return material
 
+def get_triga_lead_material():
+    shortName = 'tLEAD'
+    longName = 'lead'
+    massDensity = 11.34 # g/cc
+    fuelRadius = 10. #cm (complete guess)
+    temperature = 293.6 #K
+    temperatureIndex = 0 # X in .9Xc
+    thermalOpt = 'free'
+    elemAtomFracDict = {'Pb': 1.0}
+    #
+    chordLength = calc_chord_length(fuelRadius)
+    symDict, ZList, ZAList = get_all_isotopes(elemAtomFracDict)
+    abundanceDict = lookup_natl_abundances(ZAList)
+    #
+    material = Material(
+        shortName=shortName, longName=longName,
+        temperature=temperature, thermalOpt=thermalOpt,
+        symDict=symDict, ZList=ZList, ZAList=ZAList,
+        abundanceDict=abundanceDict, chordLength=chordLength,
+        elemAtomFracDict=elemAtomFracDict, massDensity=massDensity,
+        temperatureIndex=temperatureIndex)
+    return material
+
 ###############################################################################
 def get_bruss_enriched_rod_fuel_material():
     shortName = 'debFUEL'