Hello!
I am currently looking at what happens with various shapes on control rods when they are inserted in a fuel mixture. The control rods are deformed cylinders, and are imported as dagmc.h5m files whereas the fuel mixture is just a cylinder.
I want the control rod to be submerged in the fuel mixture, so I have tried creating just that in the code down below:
import openmc
# Define materials
fuelm = openmc.Material(name="fuelm")
fuelm.set_density('g/cm3', 10.0)
fuelm.add_nuclide('U235', 0.00229)
fuelm.add_nuclide('U238', 0.13847)
fuelm.add_nuclide('O16', 0.29)
fuelm.add_nuclide('H1', 0.57)
fuelm.add_s_alpha_beta('c_H_in_H2O')
crod = openmc.Material(name="crod")
crod.add_element('Fe', 1.0) # 100% natural iron
crod.set_density('g/cc', 7.87) # Density of iron
# Create Materials collection
mats = openmc.Materials([fuelm, crod])
mats.export_to_xml()
# Read the DAGMC geometry for the crod
# Create geometry using CSG for the fuel region
# Define surfaces
fuel_surface = openmc.openmc.model.RightCircularCylinder(center_base = (0,0,0), radius = 2, height = 5, axis = 'x')
fuel_cell = openmc.Cell(region = -fuel_surface, name='fuelm', fill = fuelm)
dagmc_univ = openmc.DAGMCUniverse(r"dagmc.h5m")
dagmc_cell = openmc.Cell(fill = crod)
# Create a root universe
root = openmc.Universe()
root.add_cells([fuel_cell,dagmc_cell])
# Create geometry and export
geometry = openmc.Geometry(root)
geometry.export_to_xml()
# Set up OpenMC settings
settings = openmc.Settings()
settings.dagmc = True
settings.batches = 10
settings.inactive = 2
settings.particles = 5000
settings.export_to_xml()
However, as I do so and plot the results I see that the fuel region is contained within the dagmc region (the total opposite). I know how to define regions, but I don’t know how to define the fuel region so that it it’s contained between the fuel surface and the control rod (dagmc model). I don’t see any way to approximate the control rod region with a simple cylinder, since the geometry is slightly more complicated than that.
Do you have any ideas. And by the way, this is what the deformed control rod looks like:
Thanks in advance.