lwr1g.py

from numpy import logspace
from xsgen.nuc_track import transmute

#############################
### General specifcations ###
#############################
reactor = "lwr1g"
plugins = ['xsgen.pre', 'xsgen.buk']
solver = 'openmc+origen'
formats = ('brightlite',)
burn_regions = 1     # Number of burnup annular regions.
burn_time = 365*10  # Number of days to burn the material [days]
time_step = 100      # Time step by which to increment the burn [days]
burn_times = [0, 3]
burn_times.extend(range(100, 7001, 100))  # we now have [0, 3, 100, 200 .. 4000]

batches = 3


################################
### Unit Cell Sepcifications ###
################################
fuel_cell_radius = 0.410
void_cell_radius = 0.4185
clad_cell_radius = 0.475
unit_cell_pitch  = 0.65635 * 2.0
unit_cell_height = 10.0

fuel_density = 10.7
clad_density = 5.87                         # Cladding Density
cool_density = 0.73                         # Coolant Density

fuel_specific_power = 40.0 / 1000.0   # Power garnered from fuel [W / g]

# LEU
initial_heavy_metal = {     # Initial heavy metal mass fraction distribution
    922350: 0.04,
    922380: 0.96,
    }

enrichment = 0.04

pnl = 0.96

# UOX
fuel_chemical_form = {                 # Dictionary of initial fuel loading.
    80160: 2.0,
    "IHM": 1.0,
    }

k_particles   = 500       # Number of particles to run per kcode cycle
k_cycles      = 130       # Number of kcode cycles to run
k_cycles_skip = 30        # Number of kcode cycles to run but not tally at the begining.

# group_structure = [1.0e-9, 10]
group_structure = logspace(-9, 1, 10)

# Temperature
# Should be a positive multiple of 300 K (ie 300, 600, 900, etc)
temperature = 600

track_nucs = ["Ac227",
              "Am241",
              "AM242",
              "BA140",
              "C14",
              "CM251",
              "CS141",
              "CS142",
              "CS147",
              "H1",
              "H3",
              "PU236",
              "PU237",
              "Pu238",
              "Pu239",
              "Pu240",
              "Pu241",
              "Th228",
              "Th229",
              "Th230",
              "Th232",
              "U230",
              "U231",
              "U232",
              "U233",
              "U234",
              "U235",
              "U236",
              "U237",
              "U238",
              "U239",
              "Zr93",
              "ZR95",
]