Nick Schenkels

Nick Schenkels, Wim Vanroose

Many inverse problems can be described by a PDE model with unknown parameters that need to be calibrated based on measurements related to its solution. This can be seen as a constrained minimization problem where one wishes to minimize the mismatch between the observed data and the model predictions, including an extra regularization term, and use the PDE as a constraint. Often, a suitable regularization parameter is determined by solving the problem for a whole range of parameters -- e.g. using the L-curve -- which is computationally very expensive. In this paper we derive two methods that simultaneously solve the inverse problem and determine a suitable value for the regularization parameter. The first one is a direct generalization of the Generalized Arnoldi Tikhonov method for linear inverse problems. The second method is a novel method based on similar ideas, but with a number of advantages for nonlinear problems.

Nick Schenkels, Wim Vanroose

In this paper we derive a Newton type method to solve the non-linear system formed by combining the Tikhonov normal equations and Morozov's discrepancy principle. We prove that by placing a bound on the step size of the Newton iterations the method will always converge to the solution. By projecting the problem onto a low dimensional Krylov subspace and using the method to solve the projected non-linear system we show that we can reduce the computational cost of the method.

Nick Schenkels, Jan Sijbers, Wim van Aarle et al.

Phase contrast tomography is an alternative to classic absorption contrast tomography that leads to higher contrast reconstructions in many applications. We review how phase contrast data can be acquired by using a combination of phase and absorption gratings. Using algebraic reconstruction techniques the object can be reconstructed from the measured data. In order to solve the resulting linear system we propose the Generalized Bidiagonal Tikhonov (GBiT) method, an adaptation of the generalized Arnoldi-Tikhonov method that uses the bidiagonal decomposition of the matrix instead of the Arnoldi decomposition. We also study the effect of the finite difference operator in the model by examining the reconstructions with either a forward difference or a central difference approximation. We validate our conclusions with simulated and experimental data.