Search

The ANSWERS® WIMS reactor physics code is being developed for whole core multiphysics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding and pin-by-pin differential depletion. An integrated thermal hydraulic solver permits differential temperature and density variations to feedback to the neutronics calculation. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver. Two coupling routes are presented and compared using a challenging PWR SMR benchmark. The first route, called GEOM, dynamically calculates the resonance shielding and homogenisation with the whole core flux solution. The second coupling route, called CAMELOT, separates the resonance shielding and pincell homogenisation from the whole core solution via generating tabulated cross sections. Both routes can use the MERLIN homogenised pin-by-pin whole core flux solver and couple to the same integrated thermal hydraulic solver, called ARTHUR. Heterogeneous differences between the neutronics and thermal hydraulics are mapped via thermal identifiers for neutronics materials and thermal regions. The ability for the integrated thermal hydraulic solver to call an external code via a Fortran-C-Python (FCP) interface is also summarised. This flexible external coupling permits one way coupling to an external fuel performance code or two way coupling to an external thermal hydraulic code.

For liquid metal-cooled fast reactors (LMFRs), improved predictive modelling is desirable to facilitate reactor licensing and operation and move towards a best estimate plus uncertainty (BEPU) approach. A key source of uncertainty in fast reactor calculations arises from the underlying nuclear data. Addressing the propagation of such uncertainties through multiphysics calculations schemes is therefore of importance, and is being addressed through international projects such as the Sodium-cooled Fast Reactor Uncertainty Analysis in Modelling (SFR-UAM) benchmark. In this paper, a methodology for propagation of nuclear data uncertainties within WIMS is presented. Uncertainties on key reactor physics parameters are calculated for selected SFR-UAM benchmark exercises, with good agreement with previous results. A methodology for coupled neutronic-thermal-hydraulic calculations within WIMS is developed, where thermal feedback is introduced to the neutronic solution through coupling with the ARTHUR subchannel code within WIMS and applied to steady-state analysis of the Horizon 2020 ESFR-SMART project reference core. Finally, integration of reactor physics and fuel performance calculations is demonstrated through linking of the WIMS reactor physics code to the TRAFIC fast reactor fuel performance code, through a Fortran-C-Python (FCP) interface. Given the 3D multiphysics calculation methodology, thermal-hydraulic and fuel performance uncertainties can ultimately be sampled alongside the nuclear data uncertainties. Together, these developments are therefore an important step towards enabling propagation of uncertainties through high fidelity, multiphysics SFR calculations and hence facilitate BEPU methodologies.