Best way to define distributed source

I’m trying to run a fixed source problem for a PWR where the source is defined for every pin at various axial levels. My original attempt was to use a MeshSpatial distribution using a RectilinearMesh but it seems that MeshSpatial is only implemented for unstructured mesh.

What would be the best way to define this? I could make an unstructured mesh that is rectilinear but it seems fairly time consuming. I could also define a Box source for every pin and axial mesh but I’m afraid this might not be very performant.

Thanks,
Ben

Hi Ben! I think @pshriwise may have plans for extending MeshSpatial to cover the other mesh types. Not sure what the ETA is on that but hopefully he can chime in. In the meantime, the other options would be a bit painful – like you said you could define a box source for each region. Right now I think that’s going to be O(n) in the number of sources. Depending on the problem though, maybe it’s not a big deal.

I can confirm that I’m working on a mesh source class that will enable additional mesh types for MeshSpatial and add a MeshSource class that will allow for unique angle/energy distributions for each element if desired. I’m currently on leave for the next couple of weeks, but once I return I’d like to submit that feature asap.

Hi, I’m a beginner with OpenMC. I want to create a core source with pins having different axial and radial power profiles, as shown in the diagram. Is it correct to set the yellow boxes in the figure as mesh, where the source strength of the two smaller areas inside mesh is defined as 0?Could you please provide me with a simple demonstration?