### Abstract

Original language | English |
---|---|

Title of host publication | PHYSOR 2016: Unifying Theory and Experiments in the 21st Century |

Publisher | American Nuclear Society ANS |

Pages | 255-269 |

ISBN (Print) | 978-0-89448-762-2 |

Publication status | Published - 2016 |

MoE publication type | A4 Article in a conference publication |

Event | International Conference on the Physics of Reactors, PHYSOR 2016: Unifying Theory and Experiments in the 21st Century - Sun Valley, United States Duration: 1 May 2016 → 5 May 2016 |

### Conference

Conference | International Conference on the Physics of Reactors, PHYSOR 2016 |
---|---|

Abbreviated title | PHYSOR 2016 |

Country | United States |

City | Sun Valley |

Period | 1/05/16 → 5/05/16 |

### Fingerprint

### Keywords

- monte carlo neutronics
- CFD
- multiphysics

### Cite this

*PHYSOR 2016: Unifying Theory and Experiments in the 21st Century*(pp. 255-269). American Nuclear Society ANS.

}

*PHYSOR 2016: Unifying Theory and Experiments in the 21st Century.*American Nuclear Society ANS, pp. 255-269, International Conference on the Physics of Reactors, PHYSOR 2016, Sun Valley, United States, 1/05/16.

**Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations.** / Tuominen, Riku; Valtavirta, Ville; Peltola, Juho; Leppänen, Jaakko.

Research output: Chapter in Book/Report/Conference proceeding › Conference article in proceedings › Scientific › peer-review

TY - GEN

T1 - Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations

AU - Tuominen, Riku

AU - Valtavirta, Ville

AU - Peltola, Juho

AU - Leppänen, Jaakko

N1 - SDA: SHP: SASUNE Nuclear Project : 100502

PY - 2016

Y1 - 2016

N2 - The main goal of this work was to couple the Monte Carlo neutronics code Serpent~2 with a CFD solver from the OpenFOAM toolbox. The coupling was implemented with the already available multi-physics interface of Serpent. The interface allows the passing of high fidelity density and temperature distributions from an external solver to Serpent and also the passing of fission power distribution from Serpent to the external solver. The coupled CFD-neutronics problem was solved by iteration. The coupling was tested by modelling a mock-up 5x5 fuel assembly cooled with water in a steady state condition at full power. The effects of boiling were neglected and the flow was modelled as single phase. The main results of the coupled calculation were high fidelity temperature and density distributions. The effect of distribution fidelity on neutronics was studied by running Serpent simulations with varying level of axial refinement. Most obvious differences were observed in the axial power density and in the collision density in the coolant, with minor differences in e.g. the effective multiplication factor. In addition the convergence of the coupled calculation was briefly studied. The default Serpent relaxation scheme based on the stochastic approximation performed well in the coupled calculations. This work can be continued in the future for example by including the effects of boiling in the CFD solution.

AB - The main goal of this work was to couple the Monte Carlo neutronics code Serpent~2 with a CFD solver from the OpenFOAM toolbox. The coupling was implemented with the already available multi-physics interface of Serpent. The interface allows the passing of high fidelity density and temperature distributions from an external solver to Serpent and also the passing of fission power distribution from Serpent to the external solver. The coupled CFD-neutronics problem was solved by iteration. The coupling was tested by modelling a mock-up 5x5 fuel assembly cooled with water in a steady state condition at full power. The effects of boiling were neglected and the flow was modelled as single phase. The main results of the coupled calculation were high fidelity temperature and density distributions. The effect of distribution fidelity on neutronics was studied by running Serpent simulations with varying level of axial refinement. Most obvious differences were observed in the axial power density and in the collision density in the coolant, with minor differences in e.g. the effective multiplication factor. In addition the convergence of the coupled calculation was briefly studied. The default Serpent relaxation scheme based on the stochastic approximation performed well in the coupled calculations. This work can be continued in the future for example by including the effects of boiling in the CFD solution.

KW - monte carlo neutronics

KW - CFD

KW - multiphysics

M3 - Conference article in proceedings

SN - 978-0-89448-762-2

SP - 255

EP - 269

BT - PHYSOR 2016: Unifying Theory and Experiments in the 21st Century

PB - American Nuclear Society ANS

ER -