Merge pull request #28 from wdhawkins/spherical_benchmark_boundary_update

Updating spherical benchmarks to use new boundary interface

Author: wdhawkins

Six_1g_spherical_benchmarks/Problem_1.py

@@ -32,6 +32,7 @@
         partitioner=PETScGraphPartitioner(type='parmetis'),
     )
     grid = meshgen.Execute()
+    grid.SetUniformBoundaryID("sphere_surface")
 
     # Get measured volumes
     volumes_per_block = grid.ComputeVolumePerBlockID()
@@ -103,12 +104,7 @@
         ],
         scattering_order=0,
         boundary_conditions=[
-            {"name": "xmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "xmax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmax", "type": "isotropic", "group_strength": bsrc},
+            {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
     )
 

Six_1g_spherical_benchmarks/Problem_2.py

@@ -68,14 +68,6 @@
         ],
         scattering_order=0,
         volumetric_sources=[mg_src],
-        boundary_conditions=[
-            {"name": "xmin", "type": "vacuum"},
-            {"name": "xmax", "type": "vacuum"},
-            {"name": "ymin", "type": "vacuum"},
-            {"name": "ymax", "type": "vacuum"},
-            {"name": "zmin", "type": "vacuum"},
-            {"name": "zmax", "type": "vacuum"}
-        ],
     )
 
     ss_solver = SteadyStateSourceSolver(problem=phys)

Six_1g_spherical_benchmarks/Problem_3.py

@@ -65,14 +65,6 @@
         ],
         scattering_order=0,
         volumetric_sources=[mg_src],
-        boundary_conditions=[
-            {"name": "xmin", "type": "vacuum"},
-            {"name": "xmax", "type": "vacuum"},
-            {"name": "ymin", "type": "vacuum"},
-            {"name": "ymax", "type": "vacuum"},
-            {"name": "zmin", "type": "vacuum"},
-            {"name": "zmax", "type": "vacuum"},
-        ],
     )
 
     ss_solver = SteadyStateSourceSolver(problem=phys)

Six_1g_spherical_benchmarks/Problem_4.py

@@ -31,6 +31,7 @@
         partitioner=PETScGraphPartitioner(type='parmetis'),
     )
     grid = meshgen.Execute()
+    grid.SetUniformBoundaryID("sphere_surface")
 
     # Create XS
     num_groups = 1
@@ -67,12 +68,7 @@
         ],
         scattering_order=0,
         boundary_conditions=[
-            {"name": "xmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "xmax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmax", "type": "isotropic", "group_strength": bsrc},
+            {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
     )
 

Six_1g_spherical_benchmarks/Problem_5.py

@@ -31,6 +31,7 @@
         partitioner=PETScGraphPartitioner(type='parmetis'),
     )
     grid = meshgen.Execute()
+    grid.SetUniformBoundaryID("sphere_surface")
 
     # Create XS
     num_groups = 1
@@ -67,12 +68,7 @@
         ],
         scattering_order=0,
         boundary_conditions=[
-            {"name": "xmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "xmax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "ymax", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmin", "type": "isotropic", "group_strength": bsrc},
-            {"name": "zmax", "type": "isotropic", "group_strength": bsrc},
+            {"name": "sphere_surface", "type": "isotropic", "group_strength": bsrc},
         ],
     )
 

Six_1g_spherical_benchmarks/Problem_6.py

@@ -68,14 +68,6 @@
         ],
         scattering_order=0,
         volumetric_sources=[mg_src],
-        boundary_conditions=[
-            {"name": "xmin", "type": "vacuum"},
-            {"name": "xmax", "type": "vacuum"},
-            {"name": "ymin", "type": "vacuum"},
-            {"name": "ymax", "type": "vacuum"},
-            {"name": "zmin", "type": "vacuum"},
-            {"name": "zmax", "type": "vacuum"},
-        ],
     )
 
     ss_solver = SteadyStateSourceSolver(problem=phys)