Adjoint functions have been used with forward functions to compute gradients in implicit (iterative) solution methods for inverse problems in optical tomography, geoscience, thermal science and other fields, but only once has this approach been used for inverse solutions to the Boltzmann transport equation. In this paper, this approach is used to develop an inverse method that requires only angle-independent flux measurements, rather than angle-dependent measurements as was done previously. The method is applied to a simplified form of the transport equation that does not include scattering. The resulting procedure uses measured values of gamma-ray fluxes of discrete, characteristic energies to determine interface locations in a multilayer shield. The method was implemented with a Newton-Raphson optimisation algorithm and it worked very well in numerical 1-D spherical test cases. A more sophisticated optimisation method would better exploit the potential of the inverse method.
[1]
Methods of solving ill-posed inverse problems
,
1983
.
[2]
Jeffrey A. Favorite,et al.
Using the Schwinger Variational Functional for the Solution of Inverse Transport Problems
,
2004
.
[3]
A. Klose,et al.
Optical tomography using the time-independent equation of radiative transfer-Part 1: Forward model
,
2002
.
[4]
S. J. Norton,et al.
A general nonlinear inverse transport algorithm using forward and adjoint flux computations
,
1997
.
[5]
C. R. Hagelberg,et al.
Mantle circulation models with variational data assimilation: inferring past mantle flow and structure from plate motion histories and seismic tomography
,
2001
.
[6]
Y. Azmy.
Post-convergence automatic differentiation of iterative schemes
,
1997
.