Code: Select all

```
div <mat> [sep <1/0>] [subz <nz> <zmin> <zmax>] [subr <nr> <rmin> <rmax>] [subs <ns> <s0>]
```

The remaining parameters divide the materials into sub-regions. In the current update the sub-division is based on cylindrical objects. Parameter "subz" divides the material into <nz> axial regions between z = <zmin> and z = <zmax>. Parameter "subr" divides the region radially into <nr> rings between r = <rmin> and r = <rmax>. The last paramter "subs" divides the region into <ns> angular sectors, with the zero position given by angle <s0>.

Here's an example...

Imagine a BWR assembly geometry consisting of a 10 x 10 lattice of fuel pins with 13 different types:

Defining a div-card for each fuel composition:

Code: Select all

```
div fuel1 sep 1
div fuel2 sep 1
...
div fuel13 sep 1
```

(some of the colors in the plot may be same or similar, but each regions now contains a different material)

Fuel types 9-13 contain gadolinium, so let's divide them into 10 rings:

Code: Select all

```
div fuel9 sep 1 subr 10 0.0 4.38E-01
div fuel10 sep 1 subr 10 0.0 4.38E-01
div fuel11 sep 1 subr 10 0.0 4.38E-01
div fuel12 sep 1 subr 10 0.0 4.38E-01
div fuel13 sep 1 subr 10 0.0 4.38E-01
```

Finally let's divide every pin into 4 angular sectors:

Code: Select all

```
div fuel1 sep 1 subs 4 45.0
div fuel2 sep 1 subs 4 45.0
...
div fuel13 sep 1 subr 10 0.0 4.38E-01 subs 4 45.0
```

Here's a closer view of 4 divided pins:

The axial division works in a similar manner. The code should be able to calculate material volumes correctly for this type of geometries. The division can be used with any configuration, but if the volumes cannot be calculated for the original cells, they cannot be calculated for divisions either. I'll post more about new volume division features later...