GitLab maintenance scheduled form Friday, 2021-06-18 5:00pm to Satursday, 2021-06-19 10:00pm CT - Services will be unavailable during this time.

Commit 520914e6 authored by Patrick Shriwise's avatar Patrick Shriwise

Updates to openmc model. Adjusting additional reflector thickness to maintain k-eff.

parent 559c61cd
......@@ -51,7 +51,7 @@ radius_pebble_flibe = radius_pebble_outer+tolerance
# Pebble centers coordinates (x,y,z)
print("Reading pseudo-random pebble centers from file pebble_centers.txt")
pebble_centers = np.loadtxt('list_pebbles_1568.csv', delimiter=',', skiprows=1)[:, 1:]
pebble_centers = np.loadtxt('pebble_centers.txt', delimiter=' ')
print("File pebble_centers.txt reading completed")
# Vessel Dimensions
......@@ -62,7 +62,7 @@ vessel_radius = np.max(np.linalg.norm(pebble_centers[:, :-1], axis=1)) + radius_
# Add additional moderation surrounding the pebble bed to reduce the
# artificially high k-eff due to reflecting boundary conditions
extra_thickness = 12.125 # cm
extra_thickness = 16.125 # cm
vessel_z_min -= extra_thickness # bottom plane vessel ; scaled by TAMU experiment factor
vessel_z_max += extra_thickness # top plane vessel ; scaled by TAMU experiment factor
vessel_radius += extra_thickness # scaled by TAMU experiment factor
......@@ -256,11 +256,12 @@ if verbose:
# -------------- Settings --------------------
settings = openmc.Settings()
settings.source = openmc.Source(space=openmc.stats.Box(*vessel_cell.bounding_box, only_fissionable=True))
settings.source = openmc.Source(space=openmc.stats.Box(*vessel_cell.bounding_box))
settings.particles = 50000
settings.inactive = 50
settings.batches = 250
settings.temperature = dict(default=573, method='interpolation', multipole=True, range=(300.0, 1500.0), tolerance=1000.0)
settings.material_cell_offsets = False
# Fuel volume calculation
volume_fuel = openmc.VolumeCalculation([m_fuel], 10000000, *c_pebble_central.bounding_box)
settings.volume_calculations = [volume_fuel]
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment