watts/examples/1App_PyARC_UnitCell/pyarc_template at 4ae5ad72cd7b47619c70cd9eb7c2fdad79941869 · watts-dev/watts · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
=arc
geometry{

    materials{
        material ( LEU ) {
            temp    = {{ temp }}
            adens{
                 aden ( u235 ) = 2.0e-2  % atom density in at/barn-cm
                 aden ( u238 ) = 2.0e-2  % atom density in at/barn-cm
                 aden ( o16 ) = 8.0e-2  % atom density in at/barn-cm
            }
            lumped_element_aden( lu35 ) = 1e-3
        }
        material ( LEU2 ) {
            temp    = {{ temp }}
            adens{
                 aden ( u236 ) = 2.0e-2  % atom density in at/barn-cm
                 aden ( u237 ) = 2.0e-2  % atom density in at/barn-cm
                 aden ( o16 ) = 8.0e-2  % atom density in at/barn-cm
            }
        }
    }
    surfaces{
        hexagon ( hex ){ orientation=y   normal = z  pitch = {{ assembly_pitch }} }
        plane ( z0 ) { z = 0.0  }
        plane ( z10 ) { z = {{ assembly_length }} }
    }

    reactor{
        lower_axial_surf=z0
        upper_axial_surf=z10

        lower_boundary_condition=reflective
        upper_boundary_condition=reflective
        reactor_hex_lattice {
            pitch         = hex
            num_ring      = 1
            fill          = [ assembly ]
        }
        assembly ( assembly ) {
            sub_assembly ( sub_assembly_name ) {
                lower_axial_surf=z0
                upper_axial_surf=z10
                material=LEU
            }
        }
    }
}
calculations{
	lumped_element_text_file( lu35 ) = "lumped.son"
    mcc3{
       xslib            = "endf7.0"
       egroupname       = ANL33
       scattering_order = 1

        cell( a ){
            associated_sub_assembly     = sub_assembly_name
        }
    }
    dif3d{
        power           = 1.0
        geometry_type   = hexagonal_infinite_lattice
        isotxs          = previous
        run_dif3d       = true
        max_axial_mesh_size = 0.05
        dif_fd_options{}
    }
}


end