Publications Details

Modeling impact-induced reactivity changes using DAG-MCNP5

There is a long literature studying the criticality of space reactors immersed in water/sand after a launch accident; however most of these studies evaluate nominal or uniformly compacted system configurations. There is less research on the reactivity consequences of impact, which can cause large structural deformation of reactor components that can result in changes in the reactivity of the system. Predicting these changes is an important component of launch safety analysis. This paper describes new features added to the DAG-MCNP5 neutronics code that allow the criticality analysis of deformed geometries. A CAD-based solid model of the reactor geometry is used to generate an initial mesh for a structural mechanics impact calculation using the PRONTO3D/PRESTO continuum mechanics codes. Boundary conditions and material specifications for the reactivity analysis are attached to the solid model that is then associated with the initial mesh representation. This geometry is then updated with the deformed finite element mesh to perturb node coordinates. DAG-MCNP5 was extended to accommodate two consequences of the large structural deformations: dead elements representing fracture, and small overlaps between adjacent volumes. The dead elements are removed during geometry initialization and adjustments are made to conseve mass. More challenging, small overlaps where adjacent mesh elements contact cause the geometric queries to become unreliable. A new point membership test was developed that is tolerant of self-intersecting volumes, and the particle tracking algorithm was adjusted to enable transport through small overlaps. These new features enable DAG-MCNP5 to perform particle transport and criticality eigenvalue calculations on both deformed mesh geometry and CAD geometry with small geometric defects. Detailed impact simulations were performed on an 85-pin space reactor model. Iin the most realistic model that included NaK coolant and water in the impact simulation, the eigenvalue was determined to increase 2.7% due to impact.