Martin Berggren

2papers

2 Papers

2.2LGJun 3
Generalized TV--$\ell_p$ Structured Priors for Bayesian $T_1$ Mapping

Disi Lin, Martin Berggren, Tommy Löfstedt

We propose an extended family of structured spatial priors that incorporates the total variation (TV) function with $\ell_p$ norms. The prior is proven to be proper and incorporated into a Bayesian regression framework to enable uncertainty quantification in $T_1$ mapping, with posterior inference performed using the No-U-Turn Sampler (NUTS). This TV--$\ell_p$ construction is proven to constitute a well-defined family of prior distributions, and it naturally enforces spatial consistency and smooth variations in the estimated parameter maps. The method was evaluated in comparison to maximum-likelihood estimation and several Bayesian alternative priors based on the uniform, Gamma, and bounded TV priors. The evaluation includes experiments on synthetic brain and cardiac $T_1$ mapping datasets, as well as a real in-vivo breast $T_1$ mapping dataset. The results show that the TV--$\ell_p$ prior yields more concentrated posterior densities, indicating reduced uncertainty. It also consistently achieves lower variance and smaller (negative) bias, leading to more reliable estimates. Overall, embedding a TV-based structured penalty along with $\ell_p$ norms in a prior in a Bayesian model improves spatial coherence in $T_1$ maps and enhances uncertainty quantification, offering a robust approach for $T_1$ mapping with uncertainties.

NAFeb 29, 2016
A Nitsche-type Method for Helmholtz Equation with an Embedded Acoustically Permeable Interface

Esubalewe Lakie Yedeg, Eddie Wadbro, Peter Hansbo et al.

We propose a new finite element method for Helmholtz equation in the situation where an acoustically permeable interface is embedded in the computational domain. A variant of Nitsche's method, different from the standard one, weakly enforces the impedance conditions for transmission through the interface. As opposed to a standard finite-element discretization of the problem, our method seamlessly handles a complex-valued impedance function $Z$ that is allowed to vanish. In the case of a vanishing impedance, the proposed method reduces to the classic Nitsche method to weakly enforce continuity over the interface. We show stability of the method, in terms of a discrete Gårding inequality, for a quite general class of surface impedance functions, provided that possible surface waves are sufficiently resolved by the mesh. Moreover, we prove an a priori error estimate under the assumption that the absolute value of the impedance is bounded away from zero almost everywhere. Numerical experiments illustrate the performance of the method for a number of test cases in 2D and 3D with different interface conditions.