A Haar wavelet method for angularly discretising the Boltzmann transport equation

Abstract A novel, hierarchical Haar wavelet basis is introduced and used to discretise the angular dimension of the Boltzmann transport equation. This is used in conjunction with a finite element subgrid scale method. This combination is then validated using two steady-state radiation transport problems, namely a 2D dogleg-duct shielding problem and the 2D C5MOX OECD/NEA benchmark. It is shown that the scheme has many similarities to a traditional equal weighted discrete ordinates ( S n ) angular discretisation, but the strong motivation for our hierarchical Haar wavelet method is the potential for adapting in angle in a simple fashion through elimination of redundant wavelets. Initial investigations of this adaptive approach are presented for a shielding and criticality eigenvalue example. It is shown that a 60% reduction in the number of angles needed on most spatial nodes - and rising up to 90% on nodes located in high streaming areas - can be attained without adversely affecting the accuracy of the solution.

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

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

[3]  R. P. Smedley-Stevenson,et al.  Linear and quadratic octahedral wavelets on the sphere for angular discretisations of the Boltzmann transport equation , 2005 .

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

[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]  R. P. Smedley-Stevenson,et al.  The Inner-Element Subgrid Scale Finite Element Method for the Boltzmann Transport Equation , 2010 .

[7]  Mohammad Hossein Heydari,et al.  Legendre wavelets method for solving fractional partial differential equations with Dirichlet boundary conditions , 2014, Appl. Math. Comput..

[8]  Stuart C. Althorpe,et al.  Wavelet-distributed approximating functional method for solving the Navier-Stokes equation , 1998 .

[9]  P. Colella,et al.  An Adaptive Mesh Refinement Algorithm for the Radiative Transport Equation , 1998 .

[10]  Xiaolin Zhou,et al.  Using divergence free wavelets for the numerical solution of the 2-D stationary Navier-Stokes equations , 2005, Appl. Math. Comput..

[11]  Wavelet-based angular dependence analysis of the heterogeneous calculation on MOX fuel lattice , 2009 .

[12]  Stanisław Ryszard Massel,et al.  Wavelet analysis for processing of ocean surface wave records , 2001 .

[13]  Youqi Zheng,et al.  Solution of neutron transport equation using Daubechies’ wavelet expansion in the angular discretization , 2008 .

[14]  Bing Li,et al.  Review: Wavelet-based numerical analysis: A review and classification , 2014 .

[15]  Edward W. Larsen,et al.  Fast iterative methods for discrete-ordinates particle transport calculations , 2002 .

[16]  Oleg V. Vasilyev,et al.  Adaptive wavelet collocation method on the shallow water model , 2014, J. Comput. Phys..

[17]  Fernando T. Pinho,et al.  Adaptive multiresolution approach for solution of hyperbolic PDEs , 2002 .

[18]  S. Saha Ray,et al.  Two-dimensional Haar wavelet Collocation Method for the solution of Stationary Neutron Transport Equation in a homogeneous isotropic medium , 2014 .

[19]  I. Daubechies.  Orthonormal bases of compactly supported wavelets , 1988 .

[20]  P. Vandergheynst,et al.  Wavelets on the 2-sphere: A group-theoretical approach , 1999 .

[21]  Christopher C. Pain,et al.  A POD reduced order model for resolving angular direction in neutron/photon transport problems , 2015, J. Comput. Phys..

[22]  A. Grossmann,et al.  DECOMPOSITION OF HARDY FUNCTIONS INTO SQUARE INTEGRABLE WAVELETS OF CONSTANT SHAPE , 1984 .

[23]  Laurent Jacques,et al.  Wavelets on the sphere: implementation and approximations , 2002 .

[24]  Peter Schröder,et al.  Spherical wavelets: efficiently representing functions on the sphere , 1995, SIGGRAPH.

[25]  Shiyou Yang,et al.  Wavelet-Galerkin method for solving parabolic equations in finite domains , 2001 .

[26]  R. P. Smedley-Stevenson,et al.  A sub‐grid scale finite element agglomeration multigrid method with application to the Boltzmann transport equation , 2012 .

[27]  A. Patra,et al.  A numerical approach based on Haar wavelet operational method to solve neutron point kinetics equation involving imposed reactivity insertions , 2014 .

[28]  Jean C. Ragusa,et al.  A two-mesh adaptive mesh refinement technique for SN neutral-particle transport using a higher-order DGFEM , 2010, J. Comput. Appl. Math..

[29]  B. Kennett,et al.  On a wavelet-based method for the numerical simulation of wave propagation , 2002 .

[30]  Kai Schneider,et al.  Numerical simulation of a mixing layer in an adaptive wavelet basis , 2000 .

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

[32]  Willi Freeden,et al.  Combined Spherical Harmonic and Wavelet Expansion—A Future Concept in Earth's Gravitational Determination , 1997 .

[33]  Hongchun Wu,et al.  An improved three-dimensional wavelet-based method for solving the first-order Boltzmann transport equation , 2009 .

[34]  Richard Kronland-Martinet,et al.  Analysis of Sound Patterns through Wavelet transforms , 1987, Int. J. Pattern Recognit. Artif. Intell..

[35]  G. Chiavassa,et al.  Two adaptive wavelet algorithms for non-linear parabolic partial differential equations ☆ , 2002 .