def test_sets(self):
xs = CrossSection()
xs.set_d(0, 1.0)
xs.set_siga(0, 2.0)
xs.set_nusigf(0, 3.0)
xs.set_xi(0, 1.0)
xs.set_d(1, 11.0)
xs.set_siga(1, 12.0)
xs.set_nusigf(1, 13.0)
xs.set_xi(1, 0.0)
xs.set_smat([[1.0, 2.0], [3.0, 4.0]])
# sm(kg, kkg)
#
# kkg
# 0 1
# kg 0 1.0 2.0
# 1 3.0 4.0
self.assertEqual(xs.dif(0), 1.0)
self.assertEqual(xs.siga(0), 2.0)
self.assertEqual(xs.nusigf(0), 3.0)
self.assertEqual(xs.xi(0), 1.0)
self.assertEqual(xs.dif(1), 11.0)
self.assertEqual(xs.siga(1), 12.0)
self.assertEqual(xs.nusigf(1), 13.0)
self.assertEqual(xs.xi(1), 0.0)
self.assertEqual(xs.sigs(0, 0), 1.0)
self.assertEqual(xs.sigs(0, 1), 2.0)
self.assertEqual(xs.sigs(1, 0), 3.0)
self.assertEqual(xs.sigs(1, 1), 4.0)
def test_sets(self):
xs = CrossSection()
xs.set_d(1.0)
xs.set_siga(2.0)
xs.set_nusigf(3.0)
self.assertEqual(xs.dif(), 1.0)
self.assertEqual(xs.siga(), 2.0)
self.assertEqual(xs.nusigf(), 3.0)
def test_container(self):
xs_fuel = CrossSection([1.36, 0.0181, 0.0279])
delta = 1.0
albedo = -1.0
geom = [{'xs':xs_fuel, 'width':100}]
cont = Container(geom, delta, albedo)
#cont.debug()
keff = 1.0
keff_old = 1.0
total_fis_src_old = 1.0
conv = 1.0e-7
for idx_outer in range(100):
total_fis_src = cont.get_total_fis_src()
norm_factor = 1.0 / (total_fis_src/keff)
cont.normalize_fis_src(norm_factor)
for idx_inner in range(4):
for color in range(2):
cont.calc(color)
cont.calc_fis_src()
total_fis_src = cont.get_total_fis_src()
keff = total_fis_src / (total_fis_src_old/keff_old)
diff = abs((keff - keff_old)/keff)
#print( keff, diff)
if(diff < conv):
break
keff_old = keff
total_fis_src_old = total_fis_src
cont.set_keff(keff)
kana = xs_fuel.nusigf() / (xs_fuel.dif() * math.pi ** 2 / 100**2 + xs_fuel.siga())
#print( 'kana = ', kana)
self.assertAlmostEqual(keff, kana, places=4)
flux = cont.get_flux_dist()
self.assertEqual(len(flux), 2) # x, y
self.assertEqual(len(flux[0]), int(geom[0]['width']/delta))
self.assertEqual(flux[0][0], delta/2.0)
self.assertEqual(flux[0][-1], geom[0]['width']-delta/2.0)
self.assertEqual(len(flux[1]), int(geom[0]['width']/delta))
def test_calculation_manager(self):
xs_fuel = CrossSection([1.36, 0.0181, 0.0279])
delta = 1.0
albedo = -1.0
geom = [{'xs': xs_fuel, 'width': 100}]
config = {'geometry': geom, 'mesh_width': delta, "albedo": albedo}
calc_man = CalculationManager(config)
calc_man.run()
keff = calc_man.get_keff()
kana = xs_fuel.nusigf() / (xs_fuel.dif() * math.pi**2 / 100**2 +
xs_fuel.siga())
#print( 'kana = ', kana)
self.assertAlmostEqual(keff, kana, places=4)
def test_onenode(self):
node = Node()
xs = CrossSection([1.36, 0.0181, 0.0279])
#xs.debug()
node.set_xs(xs)
node.set_keff( xs.nusigf() / xs.siga() )
node.calc()
for k in range(1000):
jout_xm = node.get_jout(XM)
jout_xp = node.get_jout(XP)
node.set_jin(XM, jout_xm)
node.set_jin(XP, jout_xp)
node.calc()
self.assertEqual(node.get_jout(XM), node.get_jout(XP))
self.assertEqual(node.get_jin(XM), node.get_jout(XM))
self.assertEqual(node.get_jin(XP), node.get_jout(XP))
self.assertEqual(node.get_jout(XM)+node.get_jin(XM), node.get_flux() / 2.0)
def test_container_controller(self):
xs_fuel = CrossSection([1.36, 0.0181, 0.0279])
delta = 1.0
albedo = -1.0
geom = [{'xs': xs_fuel, 'width': 100}]
container = Container(geom, delta, albedo)
#cont.debug()
controller = ContainerController(container)
controller.calc()
keff = controller.get_keff()
kana = xs_fuel.nusigf() / (xs_fuel.dif() * math.pi**2 / 100**2 +
xs_fuel.siga())
#print( 'kana = ', kana)
self.assertAlmostEqual(keff, kana, places=4)
def test_uniform_zeroflux_bc(self):
xs_fuel = CrossSection([1.36, 0.0181, 0.0279])
delta = 1.0
geom = [{'xs': xs_fuel, 'width': 100}]
nodes = []
for r in geom:
for k in range(int(r['width'] / delta)):
the_node = Node(r['xs'])
the_node.set_width(delta)
nodes.append(the_node)
keff = 1.0
keff_old = 1.0
total_fis_src_old = 1.0
conv = 1.0e-7
for idx_outer in range(2000): # outer iteration
# calculation of total fission source
total_fis_src = 0.0
for the_node in nodes:
total_fis_src += the_node.get_fis_src()
# normalize fis src to make total fis src unity
norm_factor = 1.0 / (total_fis_src / keff)
for the_node in nodes:
the_node.normalize_fis_src(norm_factor)
# inner iteration with fixed fis src
for istart in range(2): # start color (0: red, 1:black)
for ix in range(istart, len(nodes), 2):
if (ix == 0): # left boundary
jin_xm = -nodes[ix].get_jout(XM)
else:
jin_xm = nodes[ix - 1].get_jout(XP)
if (ix == len(nodes) - 1): # right boundary
jin_xp = -nodes[ix].get_jout(XP)
else:
jin_xp = nodes[ix + 1].get_jout(XM)
nodes[ix].set_jin(XM, jin_xm)
nodes[ix].set_jin(XP, jin_xp)
nodes[ix].calc()
# calculation of new fission source
for the_node in nodes:
the_node.calc_fis_src()
# calculation of total fission source
total_fis_src = 0.0
for the_node in nodes:
total_fis_src += the_node.get_fis_src()
# estimation of eigen value as a ratio of generations
keff = total_fis_src / (total_fis_src_old / keff_old)
diff = abs((keff - keff_old) / keff)
#print( idx_outer, keff, diff)
if (diff < conv):
break
keff_old = keff
total_fis_src_old = total_fis_src
# set new eigen value to all the nodes
for the_node in nodes:
the_node.set_keff(keff)
# reference as analytical solution: Eq.(6) on p.107 of the textbook
kana = xs_fuel.nusigf() / (xs_fuel.dif() * math.pi**2 / 100**2 +
xs_fuel.siga())
#print( 'kana = ', kana)
self.assertAlmostEqual(keff, kana, places=5)