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.

A Stability Estimate in the Source Problem for the Radiative Transfer Equation

AN STABILITY ESTIMATE IN 3D SOURCE PROBLEM FOR THE STATIONARY RADIATIVE TRANSFER EQUATION

Stability estimates for the inverse source problem with passive measurements

Stability for a multi-frequency inverse random source problem

Stability for the Inverse Source Problems in Elastic and Electromagnetic Waves

Recovering Point Sources for the Inhomogeneous Helmholtz Equation *

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

Stability in the inverse source problem for elastic and electromagnetic waves with multi-frequencies

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

Inverse Source Problem in Radiative Transport with Nonlinear Anisotropic Sources

An Inverse Random Source Problem for the Helmholtz Equation.

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

Stability for inverse source problems of the stochastic Helmholtz equation with a white noise

Compactness and existence theory for a general class of stationary radiative transfer equations

Conditional Stability for an Inverse Source Problem and an Application to the Estimation of Air Dose Rate of Radioactive Substances by Drone Data

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

A new nonlocal thermodynamical equilibrium radiative transfer method for cool stars

Stability estimates for the 1d multifrequency inverse medium problem

K. Schwarzschild's problem in radiation transfer theory

Stability Estimate for an Inverse Boundary Coefficient Problem in Thermal Imaging