Abstract:The criticality problem in nuclear engineering asks for the principal eigen-pair of a Boltzmann operator describing neutron transport in a reactor core. Being able to reliably design, and control such reactors requires assessing these quantities within quantifiable accuracy tolerances. In this paper we propose a paradigm that deviates from the common practice of approximately solving the corresponding spectral problem with a fixed, presumably sufficiently fine discretization. Instead, the present approach is based on first contriving iterative schemes, formulated in function space, that are shown to converge at a quantitative rate without assuming any a priori excess regularity properties, and that exploit only properties of the optical parameters in the underlying radiative transfer model. We develop the analytical and numerical tools for approximately realizing each iteration step withing judiciously chosen accuracy tolerances, verified by a posteriori estimates, so as to still warrant quantifiable convergence to the exact eigen-pair. This is carried out in full first for a Newton scheme. Since this is only locally convergent we analyze in addition the convergence of a power iteration in function space to produce sufficiently accurate initial guesses. Here we have to deal with intrinsic difficulties posed by compact but unsymmetric operators preventing standard arguments used in the finite dimensional case. Our main point is that we can avoid any condition on an initial guess to be already in a small neighborhood of the exact solution. We close with a discussion of remaining intrinsic obstructions to a certifiable numerical implementation, mainly related to not knowing the gap between the principal eigenvalue and the next smaller one in modulus.

Discrete Eigenvalues of the Radiative Transfer Equation with Legendre and Chebyshev Polynomials Solutions

Radiative Transfer Equation Solution for Many Scattering Types

Some Numerical Results of the Albedo Problem in Spherical Geometry with Linearly Anisotropic Scattering

Exact Analytical Solution of the One-Dimensional Time-Dependent Radiative Transfer Equation with Linear Scattering

Infinite Series Solution of the Time-Dependent Radiative Transfer Equation in Anisotropically Scattering Media

Exact solutions for radiative transfer with partial frequency redistribution

Numerical analysis of a spherical harmonic discontinuous Galerkin method for scaled radiative transfer equations with isotropic scattering

An adaptive moving mesh discontinuous Galerkin method for the radiative transfer equation

Accuracy Controlled Schemes for the Eigenvalue Problem of the Radiative Transfer Equation

A class of Galerkin schemes for time-dependent radiative transfer

Numerical Simulation of Polarized Light and Temperature with a Refractive Interface

Multilevel Method for Thermal Radiative Transfer Problems with Method of Long Characteristics for the Boltzmann Transport Equation

Spectral Element Method for Vector Radiative Transfer Equation

Equilibrium and Non-Equilibrium diffusion approximation for the radiative transfer equation

The Criticality Problem for the Anisotropic Scattering with Two Discrete Eigenvalues

A discrete-ordinate weak Galerkin method for radiative transfer equation

Reflective Conditions for Radiative Transfer in Integral Form with H-Matrices

Formal solutions for polarized radiative transfer. III. Stiffness and instability

Analytical Solutions of Model Kinetic Radiative Transfer Equations and Energy Equation

Two-Group Radiative Transfer Benchmarks for the Non-Equilibrium Diffusion Model

An Asymptotic-Preserving IMEX Method for Nonlinear Radiative Transfer Equation