Commit 193d250e authored by Gurecky, William's avatar Gurecky, William
Browse files

change rodMap to include nonheated and heated rods

parent 03b8380d

SubKit/build/SquareLatticeLWR_Nodal.py

+5 −4
Original line number Diff line number Diff line
@@ -87,6 +87,7 @@ class NodalChannelGeom(object): 
      """ Returns the size of nodal channel in y-direction"""

      return me.nodalChanSizes.getYsize(row, col)





class NodalChannelBaseGeom(object):

   """ Calculates geometry of a nodal channel
@@ -859,14 +860,14 @@ class SquareLatticeLWR_Nodal(CoreBuilder.Core): 
               if me.modelGuideTubes:

                   # top nodes in asm

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 0] = int(assemID)

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 1] = int(0)

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 1] = int(assemID)

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 2] = int(assemID)

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 3] = int(0)

                   rodMap[asm_i * 2 + 0][asm_j * 4 + 3] = int(assemID)

                   # bottom nodes in asm

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 0] = int(assemID)

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 1] = int(0)

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 1] = int(assemID)

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 2] = int(assemID)

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 3] = int(0)

                   rodMap[asm_i * 2 + 1][asm_j * 4 + 3] = int(assemID)

               else:

                   # top nodes in asm

                   rodMap[asm_i * 2 + 0][asm_j * 2 + 0] = int(assemID)

tests/nodal_p7_pow3d/buildDeck.py

+2 −2
Original line number Diff line number Diff line
@@ -247,8 +247,8 @@ def main(power_3d_flag=1): 
               0.78, 0.77, 0.76, 0.75, 0.00, 0.00, 0.00, 0.00]]



    builder = SquareLatticeLWR_Nodal(model, coreMap, pwr17x17, pitch, dz, solidGeoms, assemPitch=assemPitch,

                                     baffleGap=baffleGap, assemLosses=formLoss, assemLossLevels=gridLevels, rodDomains=rodDomains,

                                     modelGuideTubes=True, setCoreMaps=True)

                                     baffleGap=baffleGap, assemLosses=formLoss, assemLossLevels=gridLevels,

                                     rodDomains=rodDomains, modelGuideTubes=True, setCoreMaps=True)



    # Run to 1 second to let the solution reach steady state

    model.addTransientGroup(tEnd=1.0)