I am using the mutliphysics interface with a completely enclosed, rectangular primse OpenFOAM mesh created with blockMesh. The simulation typically runs without a hitch, but every once in awhile it crashed with lost particles:

Code: Select all

```
Tracking error: infinite loop (possibly lost particle)
Particle type : neutron
(x,y,z) = 3.36331E+01 9.04791E+02 3.46109E+00
(u,v,w) = -3.91149E-01 5.34645E-06 9.20327E-01
E = 5.24396E-01
material = UO2
totxs = 1.94493E-02
majorant = 1.94493E-02
Checking geometry errors (this may fail)...
No geometry errors.
Geometry structure at last position is printed below.
Some levels may not be printed correctly. Look for the
first level that fails: undefined region, suspicious
local coordinates, etc., and check the corresponding
geometry definition in the input.
Level 0 / universe 0 :
Cell universe / cell 1
(x,y,z) = 3.36331E+01 9.04791E+02 3.46109E+00
(u,v,w) = -3.91149E-01 5.34645E-06 9.20327E-01
Level 1 / universe 1 :
Unstructured mesh universe
(x,y,z) = 3.36331E+01 9.04791E+02 3.46109E+00
(u,v,w) = -3.91149E-01 5.34645E-06 9.20327E-01
Simulation aborted.
```

This error is similar in nature to http://ttuki.vtt.fi/serpent/viewtopic.php?f=25&t=2143

I've defined the geometry as such, very simply:

Code: Select all

```
cell 1 0 fill 1 -1
cell 2 0 outside 1
surf 1 cuboid 0 150 0 1000 -10 10
surf 3 inf
cell 4 2 void -3
solid 1 1 2
1000 4 50 20 10 2
"constant/polyMesh/points"
"constant/polyMesh/faces"
"constant/polyMesh/owner"
"constant/polyMesh/neighbour"
"constant/polyMesh/materials"
```