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

Research output: Chapter in Book/Report/Conference proceedingConference article in proceedingsScientificpeer-review

5 Citations (Scopus)

Abstract

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.
Original languageEnglish
Title of host publicationPHYSOR 2016: Unifying Theory and Experiments in the 21st Century
PublisherAmerican Nuclear Society ANS
Pages255-269
ISBN (Print)978-0-89448-762-2
Publication statusPublished - 2016
MoE publication typeA4 Article in a conference publication
EventInternational Conference on the Physics of Reactors, PHYSOR 2016: Unifying Theory and Experiments in the 21st Century - Sun Valley, United States
Duration: 1 May 20165 May 2016

Conference

ConferenceInternational Conference on the Physics of Reactors, PHYSOR 2016
Abbreviated titlePHYSOR 2016
CountryUnited States
CitySun Valley
Period1/05/165/05/16

Fingerprint

charge flow devices
boiling
physics
density distribution
temperature distribution
coolants
multiplication
iteration
fission
radiant flux density
assembly
collisions
approximation
water
simulation

Keywords

  • monte carlo neutronics
  • CFD
  • multiphysics

Cite this

Tuominen, R., Valtavirta, V., Peltola, J., & Leppänen, J. (2016). Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations. In PHYSOR 2016: Unifying Theory and Experiments in the 21st Century (pp. 255-269). American Nuclear Society ANS.
Tuominen, Riku ; Valtavirta, Ville ; Peltola, Juho ; Leppänen, Jaakko. / Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations. PHYSOR 2016: Unifying Theory and Experiments in the 21st Century. American Nuclear Society ANS, 2016. pp. 255-269
@inproceedings{d609741931eb4b5383bb058f8e3ba914,
title = "Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations",
abstract = "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.",
keywords = "monte carlo neutronics, CFD, multiphysics",
author = "Riku Tuominen and Ville Valtavirta and Juho Peltola and Jaakko Lepp{\"a}nen",
note = "SDA: SHP: SASUNE Nuclear Project : 100502",
year = "2016",
language = "English",
isbn = "978-0-89448-762-2",
pages = "255--269",
booktitle = "PHYSOR 2016: Unifying Theory and Experiments in the 21st Century",
publisher = "American Nuclear Society ANS",
address = "United States",

}

Tuominen, R, Valtavirta, V, Peltola, J & Leppänen, J 2016, Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations. in 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.

PHYSOR 2016: Unifying Theory and Experiments in the 21st Century. American Nuclear Society ANS, 2016. p. 255-269.

Research output: Chapter in Book/Report/Conference proceedingConference article in proceedingsScientificpeer-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 -

Tuominen R, Valtavirta V, Peltola J, Leppänen J. Coupling Serpent and OpenFOAM for neutronics - CFD multi-physics calculations. In PHYSOR 2016: Unifying Theory and Experiments in the 21st Century. American Nuclear Society ANS. 2016. p. 255-269