Angular adaptivity with spherical harmonics for Boltzmann transport

Abstract This paper describes an angular adaptivity algorithm for Boltzmann transport applications which uses Pn and filtered Pn expansions, allowing for different expansion orders across space/energy. Our spatial discretisation is specifically designed to use less memory than competing DG schemes and also gives us direct access to the amount of stabilisation applied at each node. For filtered Pn expansions, we then use our adaptive process in combination with this net amount of stabilisation to compute a spatially dependent filter strength that does not depend on a priori spatial information. This applies heavy filtering only where discontinuities are present, allowing the filtered Pn expansion to retain high-order convergence where possible. Regular and goal-based error metrics are shown and both the adapted Pn and adapted filtered Pn methods show significant reductions in DOFs and runtime. The adapted filtered Pn with our spatially dependent filter shows close to fixed iteration counts and up to high-order is even competitive with P0 discretisations in problems with heavy advection.

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

[2]  Christopher C. Pain,et al.  Goal-based angular adaptivity applied to a wavelet-based discretisation of the neutral particle transport equation , 2015, J. Comput. Phys..

[3]  Mark A. Goffin,et al.  Solving the Boltzmann transport equation with multigrid and adaptive space/angle discretisations , 2015 .

[4]  S. Dargaville,et al.  Scalable angular adaptivity for Boltzmann transport , 2019, J. Comput. Phys..

[5]  Luciano Rezzolla,et al.  A new spherical harmonics scheme for multi-dimensional radiation transport I. Static matter configurations , 2012, J. Comput. Phys..

[6]  Christopher C. Pain,et al.  A Haar wavelet method for angularly discretising the Boltzmann transport equation , 2018, Progress in Nuclear Energy.

[7]  Martin J. Mohlenkamp.  A fast transform for spherical harmonics , 1997 .

[8]  Mats G. Larson,et al.  The Finite Element Method: Theory, Implementation, and Applications , 2013 .

[9]  HyeongKae Park.  Coupled Space-Angle Adaptivity and Goal-Oriented Error Control for Radiation Transport Calculations , 2006 .

[10]  Jon Hill,et al.  Adaptive Haar wavelets for the angular discretisation of spectral wave models , 2016, J. Comput. Phys..

[11]  V. Laboure,et al.  Improved fully-implicit spherical harmonics methods for first and second order forms of the transport equation using Galerkin Finite Element , 2016 .

[12]  T. Hughes,et al.  The variational multiscale method—a paradigm for computational mechanics , 1998 .

[13]  Jean-Luc Guermond,et al.  A robust SN-DG-approximation for radiation transport in optically thick and diffusive regimes , 2012, J. Comput. Phys..

[14]  T. Grasser,et al.  Adaptive variable-order spherical harmonics expansion of the Boltzmann Transport Equation , 2011, 2011 International Conference on Simulation of Semiconductor Processes and Devices.

[15]  R. T. Ackroyd,et al.  Discontinuous finite elements for neutron transport analysis , 1986 .

[16]  Ryan G. McClarren,et al.  Implicit filtered PN for high-energy density thermal radiation transport using discontinuous Galerkin finite elements , 2016, J. Comput. Phys..

[17]  Kerstin Küpper,et al.  Convergence of filtered spherical harmonic equations for radiation transport , 2016 .

[18]  C. C. Pain,et al.  An efficient space-angle subgrid scale discretisation of the neutron transport equation , 2016 .

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

[20]  Cory D. Hauck,et al.  Positivity Limiters for Filtered Spectral Approximations of Linear Kinetic Transport Equations , 2018, J. Sci. Comput..

[21]  R. T. Ackroyd,et al.  Composite finite element solutions for neutron transport , 1988 .

[22]  HyeongKae Park,et al.  Coupled Space-Angle Adaptivity for Radiation Transport Calculations , 2009 .

[23]  Thomas A. Brunner,et al.  Forms of Approximate Radiation Transport , 2002 .

[24]  R. P. Smedley-Stevenson,et al.  Self-Adaptive Spherical Wavelets for Angular Discretizations of the Boltzmann Transport Equation , 2008 .

[25]  A. Adam.  Finite element, adaptive spectral wave modelling , 2016 .

[26]  Adam S. Candy.  Subgrid scale modelling of transport processes , 2008 .

[27]  R. P. Smedley-Stevenson,et al.  Riemann boundary conditions for the Boltzmann transport equation using arbitrary angular approximations , 2011 .

[28]  Eugene Fiume,et al.  Efficient and accurate rotation of finite spherical harmonics expansions , 2012, J. Comput. Phys..

[29]  Christopher C. Pain,et al.  A goal-based angular adaptivity method for thermal radiation modelling in non grey media , 2017 .

[30]  Thomas J. R. Hughes,et al.  A multiscale discontinuous Galerkin method with the computational structure of a continuous Galerkin method , 2006 .

[31]  Ryan G. McClarren,et al.  Robust and accurate filtered spherical harmonics expansions for radiative transfer , 2010, J. Comput. Phys..

[32]  A. S. Shirani,et al.  Hybrid space–angle adaptivity for whole-core particle transport calculations , 2015 .