136 research outputs found
Angular adaptivity with spherical harmonics for Boltzmann transport
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.Comment: arXiv admin note: text overlap with arXiv:1901.0492
Anisotropic Adaptivity and Subgrid Scale Modelling for the Solution of the Neutron Transport Equation with an Emphasis on Shielding Applications
This thesis demonstrates advanced new discretisation and adaptive meshing technologies that improve the accuracy and stability of using finite element discretisations applied to the Boltzmann transport equation (BTE). This equation describes the advective transport of neutral particles such as neutrons and photons within a domain. The BTE is difficult to solve, due to its large phase space (three dimensions of space, two of angle and one each of energy and time) and the presence of non-physical oscillations in many situations. This work explores the use of a finite element method that combines the advantages of the two schemes: the discontinuous and continuous Galerkin methods. The new discretisation uses multiscale (subgrid) finite elements that work locally within each element in the finite element mesh in addition to a global, continuous, formulation. The use of higher order functions that describe the variation of the angular flux over each element is also explored using these subgrid finite element schemes. In addition to the spatial discretisation, methods have also been developed to optimise the finite element mesh in order to reduce resulting errors in the solution over the domain, or locally in situations where there is a goal of specific interest (such as a dose in a detector region).
The chapters of this thesis have been structured to be submitted individually for journal publication, and are arranged as follows. Chapter 1 introduces the reader to motivation behind the research contained within this thesis. Chapter 2 introduces the forms of the BTE that are used within this thesis. Chapter 3 provides the methods that are used, together with examples, of the validation and verification of the software that was developed as a result of this work, the transport code RADIANT. Chapter 4 introduces the inner element subgrid scale finite element discretisation of the BTE that forms the basis of the discretisations within RADIANT and explores its convergence and computational times on a set of benchmark problems. Chapter 5 develops the error metrics that are used to optimise the mesh in order to reduce the discretisation error within a finite element mesh using anisotropic adaptivity that can use elongated elements that accurately resolves computational demanding regions, such as in the presence of shocks. The work of this chapter is then extended in Chapter 6 that forms error metrics for goal based adaptivity to minimise the error in a detector response. Finally, conclusions from this thesis and suggestions for future work that may be explored are discussed in Chapter 7.Open Acces
Scalable angular adaptivity for Boltzmann transport
This paper describes an angular adaptivity algorithm for Boltzmann transport
applications which for the first time shows evidence of
scaling in both runtime and memory usage, where is the number of adapted
angles. This adaptivity uses Haar wavelets, which perform structured
-adaptivity built on top of a hierarchical P FEM discretisation of a 2D
angular domain, allowing different anisotropic angular resolution to be applied
across space/energy. Fixed angular refinement, along with regular and
goal-based error metrics are shown in three example problems taken from
neutronics/radiative transfer applications. We use a spatial discretisation
designed to use less memory than competing alternatives in general applications
and gives us the flexibility to use a matrix-free multgrid method as our
iterative method. This relies on scalable matrix-vector products using Fast
Wavelet Transforms and allows the use of traditional sweep algorithms if
desired
A geometry preserving, conservative, mesh-to-mesh isogeometric interpolation algorithm for spatial adaptivity of the multigroup, second-order even-parity form of the neutron transport equation
In this paper a method is presented for the application of energy-dependent spatial meshes applied to the multigroup, second-order, even-parity form of the neutron transport equation using Isogeometric Analysis (IGA). The computation of the inter-group regenerative source terms is based on conservative interpolation by Galerkin projection. The use of Non-Uniform Rational B-splines (NURBS) from the original computer-aided design (CAD) model allows for efficient implementation and calculation of the spatial projection operations while avoiding the complications of matching different geometric approximations faced by traditional finite element methods (FEM). The rate-of-convergence was verified using the method of manufactured solutions (MMS) and found to preserve the theoretical rates when interpolating between spatial meshes of different refinements. The scheme’s numerical efficiency was then studied using a series of two-energy group pincell test cases where a significant saving in the number of degrees-of-freedom can be found if the energy group with a complex variation in the solution is refined more than an energy group with a simpler solution function. Finally, the method was applied to a heterogeneous, seven-group reactor pincell where the spatial meshes for each energy group were adaptively selected for refinement. It was observed that by refining selected energy groups a reduction in the total number of degrees-of-freedom for the same total L2 error can be obtained
- …