Hi,
I am trying to couple neutronics, heat conduction and fluid dynamics using Cardinal for my HCPB fusion blanket; in the attempt to model conjugate heat transfer. The issue is I cannot define any square pipes using MOOSEs thermal hydrauic module so I had to convert them into cylinders. This made the element count far beyond the element limit, so I came up with solution of stitching blocks together using MOOSEs stitchmeshgenerator. I had thought this worked at first until I saw that there were discontinuities in the temperature field. I believe that this is due to the nodes on many of the surfaces that I had stitched together not being in the same place. Here is an image of the temperature field.
I used openmc2cad (offered by Cardinal) to convert the csg to a journal file, which I then loaded in Cubit. For simplicity, I modelled 2 sets of 12 pipes going through a steel plate, which is supposed to cool another solid that is sandwiched between them. I cut the plates into 3 equal sections to get around the element limit; this can be seen in the discontinuity image. I will add a link to my GitHub that has the resulting exodus file; and some other files that I used to make the mesh. Github link:https://github.com/mias1234/Neutronics-and-tritium-breeding/tree/main/CHT
I do not have access to the cubit files as they are on my personal computer, but I will upload them if need be.