I have been facing an issue with building a lattice & plotting, but I don’t know whether the issue here is plot or lattice related, when I create a lattice out of a fuel plate that has half the spacin of water between plates, I get inconsistency in the placement of the fuel meat, for some plates it is deviated forward and for some it is deviated backwards, the image shows that clearly, something tells me that it is just a plotting issue as the fuel meet is only a few millimeters thick and the OpenMC cant capture these fine details when plotting, but I feel like it is a bigger issue and has to do with the way I defined the lattice, I will also publish another post about how that lattice produces an error when running, so please take a look at that if you can. The specific portion of the code that builds the geometry will be included below, in addition toa photo of the plot + the full code.
fuel_top = openmc.ZPlane(z0=30)
fuel_bottom = openmc.ZPlane(z0=-30)
fuel_right = openmc.XPlane(x0=2.98)
fuel_left = openmc.XPlane(x0=-2.98)
fuel_front = openmc.YPlane(y0=0.0255)
fuel_back = openmc.YPlane(y0=-0.0255)
fuel_region = +fuel_bottom & -fuel_top & +fuel_left & -fuel_right & -fuel_front & +fuel_back
fuel_cell = openmc.Cell(region=fuel_region, fill=U3Si2Al)
clad_top = openmc.ZPlane(z0=32.55)
clad_bottom = openmc.ZPlane(z0=-32.55)
clad_right = openmc.XPlane(x0=3.615)
clad_left = openmc.XPlane(x0=-3.615)
clad_front = openmc.YPlane(y0=0.0635)
clad_back = openmc.YPlane(y0=-0.0635)
for surf in [clad_top, clad_bottom, clad_right, clad_left]:
surf.boundary_type = ‘vacuum’
clad_region = +clad_bottom & -clad_top & +clad_left & -clad_right & -clad_front & +clad_back & ~fuel_region
clad_cell = openmc.Cell(region=clad_region, fill=Al6061)
coolant_front = openmc.YPlane(y0=0.2045)
coolant_back = openmc.YPlane(y0=-0.2045)
coolant_region = -clad_top & + clad_bottom & -coolant_front & +coolant_back & +clad_left & -clad_right & ~clad_region & ~ fuel_region
coolant_cell = openmc.Cell(region=coolant_region, fill=water)
infinite_water = openmc.Cell(fill= water)
infinite_water_universe = openmc.Universe(cells=[infinite_water])
fuel_plate = openmc.Universe(cells=[fuel_cell, clad_cell, coolant_cell])
fuel_assembly = openmc.RectLattice(name = ‘fuel assembly’)
fuel_assembly.lower_left = (-3.615, -2.863)
fuel_assembly.pitch = (7.23, 0.409)
fuel_assembly.universes = [[fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate], [fuel_plate]]
fuel_assembly.outer = infinite_water_universe
fuel_assembly_cell = openmc.Cell()
fuel_assembly_cell.fill = fuel_assembly
fuel_assembly_universe = openmc.Universe()
fuel_assembly_universe.add_cell(fuel_assembly_cell)
geometry = openmc.Geometry(root=fuel_assembly_universe)
geometry.export_to_xml()
main code.py (4.6 KB)
Please zoom in to the image and see the inconsistency, the first fuel plate from the top is placed forward not exactly in the middle, as we go to the bottom, plates start getting deviated to the bottom, so that is a pattern that indicates a problem with the lattice definition, but I am still unable to figure out what the issue is exactly.