Problem with DAGMC Simulation

Hello,

I wanted to perform a DAGMC simulation of a geometry as follows:

It is a very simple geometry with two wall thicknesses and a floor, all made of concrete.

Here is how I proceed to convert and simulate the geometry (.step file) created in a CAD software:

First, I use GMSH to create the mesh, and then I use Python:

First function:

convert_file.py

Conversion from .msh to .h5m, applying the different materials to the 3 volumes (WALL1, WALL2, floor).

Second function:

gmsh_sim.py

The goal of my import is to have a box surrounding my CAD geometry by about 1 meter on all sides. I also want to fill the box with air, knowing that the concrete materials have already been assigned to the volumes.

I perform my simulation and I do not get anything that resembles what might be expected. There is no attenuation in the flux at the walls.

In my opinion, the function that collects data from the output file tallies.out is correct as I have used it in other simulations and it displayed expected results. I think that the way I coded the geometry import in gmsh_sim.py and the setup of the cells and universes might be incorrect, and this could be the source of the wrong results.

I have attached to this email my code, an image of my geometry, and the plot obtained which is not the expected one.

Could you please help me?

Have a good day.

Best regards.
wrong_plot

convert_files.py (1.8 KB)
gmsh_sim.py (8.9 KB)