arrow
返回

ITERATIVE METHODS FOR NEUTRON TRANSPORT EIGENVALUE PROBLEMS

delete2011-01-01
delete10
delete
OA
AI
F
Fynn Scheben *
I
Ivan G. Graham
DOI:10.1137/100799022delete
delete原文链接
delete分享
delete收藏
查看原文
摘要

摘要

En 中文
We discuss iterative methods for computing criticality in nuclear reactors. In general this requires the solution of a generalized eigenvalue problem for an unsymmetric integro-differential operator in six independent variables, modeling transport, scattering, and fission, where the dependent variable is the neutron angular flux. In engineering practice this problem is often solved iteratively, using some variant of the inverse power method. Because of the high dimension, matrix representations for the operators are often not available and the inner solves needed for the eigenvalue iteration are implemented by matrix-free inner iterations. This leads to technically complicated inexact iterative methods, for which there appears to be no published rigorous convergence theory. For the monoenergetic homogeneous model case with isotropic scattering and vacuum boundary conditions, we show that, before discretization, the general nonsymmetric eigenproblem for the angular flux is equivalent to a certain related eigenproblem for the scalar flux, involving a symmetric positive definite weakly singular integral operator (in space only). This correspondence to a symmetric problem (in a space of reduced dimension) permits us to give a convergence theory for inexact inverse iteration and related methods. In particular this theory provides rather precise criteria on how accurate the inner solves need to be in order for the whole iterative method to converge. We also give examples of discretizations which have a corresponding symmetric finite-dimensional reduced form. The theory is illustrated with numerical examples for several test problems of physical relevance, using GMRES as the inner solver.
Keyword:
neutron transport
criticality
generalized eigenvalue problem
symmetry
inexact inverse iteration
AI总结

AI总结

对已上传原文的论文进行重点信息的提取,主要内容包括:简要概述、研究摘要、背景介绍、关键亮点、图文解析、展望与总结。

期刊

SIAM Journal on Scientific Computing 封面图
SIAM Journal on Scientific Computing
IF:
2.6
论文数:
5.1K
被引数:
1.8W

机构

U
university of bath
学者数:
1.1W
论文数: 1.3W
被引数: 13