Antti Hannukainen

Antti Hannukainen, Nuutti Hyvönen, Helle Majander et al.

Quantitative photoacoustic tomography is an emerging imaging technique aimed at estimating the distribution of optical parameters inside tissues from photoacoustic images, which are formed by combining optical information and ultrasonic propagation. This optical parameter estimation problem is ill-posed and needs to be approached within the framework of inverse problems. Photoacoustic images are three-dimensional and high-resolution. Furthermore, high-resolution reconstructions of the optical parameters are targeted. Therefore, in order to provide a practical method for quantitative photoacoustic tomography, the inversion algorithm needs to be able to perform successfully with problems of prominent size. In this work, an efficient approach for the inverse problem of quantitative photoacoustic tomography is proposed, assuming an edge-preferring prior for the optical parameters. The method is based on iteratively combining priorconditioned LSQR with a lagged diffusivity step and a linearisation of the measurement model, with the needed multiplications by Jacobians performed in a matrix-free manner. The algorithm is tested with three-dimensional numerical simulations. The results show that the approach can be used to produce accurate and high quality estimates of absorption and diffusion in complex three-dimensional geometries with moderate computation time and cost.

Antti Hannukainen, Nuutti Hyvönen, Lauri Mustonen

This work tackles an inverse boundary value problem for a $p$-Laplace type partial differential equation parametrized by a smoothening parameter $τ\geq 0$. The aim is to numerically test reconstructing a conductivity type coefficient in the equation when Dirichlet boundary values of certain solutions to the corresponding Neumann problem serve as data. The numerical studies are based on a straightforward linearization of the forward map, and they demonstrate that the accuracy of such an approach depends nontrivially on $1 < p < \infty$ and the chosen parametrization for the unknown coefficient. The numerical considerations are complemented by proving that the forward operator, which maps a Hölder continuous conductivity coefficient to the solution of the Neumann problem, is Fréchet differentiable, excluding the degenerate case $τ=0$ that corresponds to the classical (weighted) $p$-Laplace equation.

Antti Hannukainen

In this paper, we analyze the convergence of the preconditioned GMRES method for the first order finite element discretizations of the Helmholtz equation in media with losses. We consider a Laplace preconditioner, an inexact Laplace preconditioner and a two-level preconditioner. Our analysis is based on bounding the field of values of the preconditioned system matrix in the complex plane. The analysis takes the non-normal nature of the linear system naturally into account and allows us to easily consider certain type of inexact Laplace preconditioners via a perturbation argument. For the two-level preconditioner, our convergence analysis takes into account a media, which has not been considered in previous works.

Tom Gustafsson, Antti Hannukainen, Vili Kohonen et al.

We present guidelines for deriving new Nitsche Finite Element Methods to enforce equality and inequality constraints that act on the value of the unknown mechanical quantity. We first formulate the problem as a stabilized finite element method for the saddle point formulation where a Lagrange multiplier enforces the underlying constraint. The Nitsche method is then presented in a general minimization form, suitable for adding constraints to nonlinear finite element methods and allowing straightforward computational implementation with automatic differentation. This extends the method beyond classical boundary condition enforcement. To validate these ideas, we present Nitsche formulations for a range of problems in solid mechanics and give numerical evidence of the convergence rates of the Nitsche method.

A. Hannukainen, J. -C. Mourrat, H. Stoppels

We present an efficient method for the computation of homogenized coefficients of divergence-form operators with random coefficients. The approach is based on a multiscale representation of the homogenized coefficients. We then implement the method numerically using a finite-element method with hierarchical hybrid grids, which is a semi-implicit method allowing for significant gains in memory usage and execution time. Finally, we demonstrate the efficiency of our approach on two- and three-dimensional examples, for piecewise-constant coefficients with corner discontinuities. For moderate ellipticity contrast and for a precision of a few percentage points, our method allows to compute the homogenized coefficients on a laptop computer in a few seconds, in two dimensions, or in a few minutes, in three dimensions.

Negar M. Harandi, Daniel Aalto, Antti Hannukainen et al.

A state-of-the-art 1D acoustic synthesizer has been previously developed, and coupled to speaker-specific biomechanical models of oropharynx in ArtiSynth. As expected, the formant frequencies of the synthesized vowel sounds were shown to be different from those of the recorded audio. Such discrepancy was hypothesized to be due to the simplified geometry of the vocal tract model as well as the one dimensional implementation of Navier-Stokes equations. In this paper, we calculate Helmholtz resonances of our vocal tract geometries using 3D finite element method (FEM), and compare them with the formant frequencies obtained from the 1D method and audio. We hope such comparison helps with clarifying the limitations of our current models and/or speech synthesizer.

Antti Hannukainen

Most finite element methods for solving time-harmonic wave-propagation problems lead to a linear system with a non-normal coefficient matrix. The non-normality is due to boundary conditions and losses. One way to solve these systems is to use a preconditioned iterative method. Detailed mathematical analysis of the convergence properties of these methods is important for developing new and understanding old preconditioners. Due to non-normality, there is currently very little existing literature in this direction. In this paper, we study the convergence of GMRES for such systems by deriving inclusion and exclusion regions for the pseudospectrum of the coefficient matrix. All analysis is done a priori by relating the properties of the weak problem to the coefficient matrix. The inclusion is derived from the stability properties of the problem and the exclusion is established via field of values and boundedness of the weak form. The derived tools are applied to estimate the pseudospectrum of time-harmonic Helmholtz equation with first-order absorbing boundary conditions, with and without a shifted-Laplace preconditioner.