CSG and DAGMC geometry

User Support
1

Hi,
As previously announced as one of last release new features, it should be now possible to combine DAGMC and CSG in one single model. I’ve tested it in a simple model, and even though the geometry seems correct from plotting, when I run openmc the CSG geometry gets ignored. I attach below the part of the script that defines the two geometry cells and few images. Am I missing something here? Thank you for your help.

graveyard = openmc.Sphere(surface_id=991, r=201, boundary_type='vacuum')
fuel_cyl = openmc.Cylinder(surface_id=992, r=1.5)
fuel_top = openmc.ZPlane(surface_id=993, z0=0)
fuel_bot = openmc.ZPlane(surface_id=994, z0=-100)
csg_cell = openmc.Cell(cell_id=995 , region= -fuel_cyl & -fuel_top & +fuel_bot, fill= fuel)

cad_univ = openmc.DAGMCUniverse(filename=h5m,universe_id=996 )
cad_cell = openmc.Cell(cell_id=997 , region= -graveyard, fill= cad_univ)

root = openmc.Universe(universe_id=998)
root.add_cells([csg_cell,cad_cell])
geometry = openmc.Geometry(root)
geometry.export_to_xml()
2

Hi again,
The problem was indeed I was not excluding the csg region from the cad region.
Corrected script and output below.
Thanks.

graveyard = openmc.Sphere(surface_id=991, r=201, boundary_type='vacuum')
fuel_cyl = openmc.Cylinder(surface_id=992, r=1.5)
fuel_top = openmc.ZPlane(surface_id=993, z0=0)
fuel_bot = openmc.ZPlane(surface_id=994, z0=-100)
csg_cell = openmc.Cell(cell_id=995 , region= -fuel_cyl & -fuel_top & +fuel_bot, fill= fuel)

cad_univ = openmc.DAGMCUniverse(filename=h5m,universe_id=996 )
cad_cell = openmc.Cell(cell_id=997 , region= -graveyard | + fuel_cyl | -fuel_bot | +fuel_top, fill= cad_univ)

root = openmc.Universe(universe_id=998)
root.add_cells([csg_cell,cad_cell])
geometry = openmc.Geometry(root)
geometry.export_to_xml()

flux
flux1920×1440 64.5 KB

1 Like