In Maximum Entropy and Bayesian Methods, G. J. Erickson, J. T. Rychert, and C. R. Smith, eds., pp, 121-135 (Kluwer Academic, Dordrecht, 1998)

Inversion based on computational simulations

Kenneth M. Hanson, Gregory S. Cunningham, and Suhail S. Saquib
Los Alamos National Laboratory


A standard approach to solving inversion problems that involve many parameters uses gradient-based optimization to find the parameters that best match the data. We will discuss enabling techniques that facilitate application of this approach to large-scale computational simulations, which are the only way to investigate many complex physical phenomena. Such simulations may not seem to lend themselves to calculation of the gradient with respect to numerous parameters. However, adjoint differentiation allows one to efficiently compute the gradient of an objective function with respect to all the variables of a simulation. When combined with advanced gradient-based optimization algorithms, adjoint differentiation permits one to solve very large problems of optimization or parameter estimation. These techniques will be illustrated through the simulation of the time-dependent diffusion of infrared light through tissue, which has been used to perform optical tomography. The techniques discussed have a wide range of applicability to modeling including the optimization of models to achieve a desired design goal.

Keywords: simulation, inversion, reconstruction, adjoint differentiation, sensitivity analysis, optimization, model validation

Get full paper (pdf, 228 KB)
Return to publication list
Send e-mail to author at