A combined immersed body and adaptive mesh method for simulating neutron transport within complex structures

Abstract This article describes a new adaptive immersed body method for the efficient and accurate modelling of neutron transport within geometrically complex domains. It combines two techniques of immersed body projections and self-adaptive mesh refinement to form a unique method that can efficiently resolve problems that contain complex internal structures. The approach allows complete freedom of where mesh resolution is placed and the meshes, which do not need to conform to the problem structure, are optimised for resolving the physics of the problem. Importantly the benefits can be found for problems with many internal structures such as fuel pins, control rods or cooling pipes and channels. Standard finite element meshes typically become large as they capture the geometrical detail and as a result solving times typically increase. This method overcomes this issue and allows for smaller meshes that are optimised to increase solution accuracy, as demonstrated in the numerical examples.

[1]  C C Pain,et al.  Anisotropic mesh adaptivity for multi-scale ocean modelling , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[2]  Matthew D. Eaton,et al.  The application of isogeometric analysis to the neutron diffusion equation for a pincell problem with an analytic benchmark , 2012 .

[3]  J. Duderstadt,et al.  Nuclear reactor analysis , 1976 .

[4]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[5]  Christopher C. Pain,et al.  Minimising the error in eigenvalue calculations involving the Boltzmann transport equation using goal-based adaptivity on unstructured meshes , 2013, J. Comput. Phys..

[6]  W. H. Reed,et al.  Triangular mesh methods for the neutron transport equation , 1973 .

[7]  C.R.E. de Oliveira,et al.  Tetrahedral mesh optimisation and adaptivity for steady-state and transient finite element calculations , 2001 .

[8]  R. P. Smedley-Stevenson,et al.  Goal-based angular adaptivity applied to the spherical harmonics discretisation of the neutral particle transport equation , 2014 .

[9]  T. Hughes,et al.  Isogeometric analysis : CAD, finite elements, NURBS, exact geometry and mesh refinement , 2005 .

[10]  Dimitrios Pavlidis,et al.  Modelling of fluid–solid interactions using an adaptive mesh fluid model coupled with a combined finite–discrete element model , 2012, Ocean Dynamics.

[11]  Axelle Viré,et al.  How to float a wind turbine , 2012, Reviews in Environmental Science and Bio/Technology.

[12]  Liangzhi Cao,et al.  A spherical harmonics—Finite element discretization of the self-adjoint angular flux neutron transport equation , 2007 .

[13]  Jim E. Morel,et al.  A Self-Adjoint Angular Flux Equation , 1999 .

[14]  C. C. Pain,et al.  An immersed body method for coupled neutron transport and thermal hydraulic simulations of PWR assemblies , 2014 .

[15]  Matthew D. Piggott,et al.  Conservative interpolation between unstructured meshes via supermesh construction , 2009 .

[16]  Wolfgang Bangerth,et al.  Three-dimensional h-adaptivity for the multigroup neutron diffusion equations , 2009 .

[17]  Antony J. H. Goddard,et al.  Simulated transient dynamics and heat transfer characteristics of the water boiler nuclear reactor – SUPO – with cooling coil heat extraction , 2012 .

[18]  R. P. Smedley-Stevenson,et al.  The Inner-Element Subgrid Scale Finite Element Method for the Boltzmann Transport Equation , 2010 .

[19]  S. Sherwin,et al.  An isoparametric approach to high-order curvilinear boundary-layer meshing , 2015 .

[20]  Miguel O. Bernabeu,et al.  Anisotropic mesh adaptivity for cardiac electrophysiology , 2010, ICCS.

[21]  C.R.E. de Oliveira,et al.  An arbitrary geometry finite element method for multigroup neutron transport with anisotropic scattering , 1986 .

[22]  Jean C. Ragusa,et al.  A simple Hessian-based 3D mesh adaptation technique with applications to the multigroup diffusion equations , 2008 .

[23]  Patrick E. Farrell.  The addition of fields on different meshes , 2011, J. Comput. Phys..

[24]  Danny Lathouwers,et al.  Goal-oriented spatial adaptivity for the SN equations on unstructured triangular meshes , 2011 .

[25]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[26]  R. T. Ackroyd,et al.  Finite Element Methods for Particle Transport: Applications to Reactor and Radiation Physics , 1997 .

[27]  R. P. Smedley-Stevenson,et al.  The immersed body supermeshing method for modelling reactor physics problems with complex internal structures , 2014 .

[28]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[29]  A. R. Owens,et al.  Discontinuous isogeometric analysis methods for the first-order form of the neutron transport equation with discrete ordinate (SN) angular discretisation , 2016, J. Comput. Phys..

[30]  Christopher C. Pain,et al.  Goal based mesh adaptivity for fixed source radiation transport calculations , 2013 .

[31]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[32]  Danny Lathouwers,et al.  Spatially adaptive eigenvalue estimation for the SN equations on unstructured triangular meshes , 2011 .

[33]  Christopher C. Pain,et al.  A new computational framework for multi‐scale ocean modelling based on adapting unstructured meshes , 2008 .

[34]  Patrick E. Farrell,et al.  Conservative interpolation between volume meshes by local Galerkin projection , 2011 .