Researcher profile

Barbara Kaltenbacher

Barbara Kaltenbacher contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2026arXiv

Existence, uniqueness, and numerical solutions of the nonlinear periodic Westervelt equation

In this paper, we study the nonlinear periodic Westervelt equation with excitations located within a bounded domain in $\mathbb{R}^d$, where $d \in \{2,3\}$, subject to Robin boundary conditions. This problem is of particular interest for advancing imaging techniques that exploit nonlinearity of the acoustic propagation. We establish the existence and uniqueness of solutions in both the linear and the nonlinear setting, thereby allowing for spatially varying coefficients as relevant in quantitative imaging. Derivation of a multiharmonic formulation enables us to show the generation of higher harmonics (that is, responses at multiples of the fundamental frequency) due nonlinear wave propagation. An iterative scheme for solving the resulting system is proposed that relies on successive resolution of these higher harmonics, and its convergence under smallness conditions on the excitation is proven. Furthermore, we investigate the numerical solution of the resulting system of Helmholtz equations, employing a conforming finite element method for its discretization. Through an implementation of the proposed methodology, we illustrate how acoustic waves propagate in nonlinear media. This study aims to enhance our understanding of ultrasound propagation dynamics, which is essential for obtaining high-quality images from limited in vivo and boundary measurements.

preprint2022arXiv

Determining damping terms in fractional wave equations

This paper deals with the inverse problem of recovering an arbitrary number of fractional damping terms in a wave equation. We develop several approaches on uniqueness and reconstruction, some of them relying on Tauberian theorems on the relation between the asymptotics of solutions in time and Laplace domain. Also the possibility of additionally recovering space dependent coefficients or initial data is discussed. The resulting methods for reconstructing coefficients and fractional orders in these terms are tested numerically. Additionally, we provide an analysis of the forward problem, a multiterm fractional wave equation.

preprint2022arXiv

Time-fractional Moore-Gibson-Thompson equations

In this paper, we consider several time-fractional generalizations of the Jordan-Moore-Gibson-Thompson (JMGT) equations in nonlinear acoustics as well as their linear Moore-Gibson-Thompson (MGT) versions. Following the procedure described in Jordan (2014), these time-fractional acoustic equations are derived from four fractional versions of the Maxwell-Cattaneo law in Compte and Metzler (1997). Additionally to providing well-posedness results for each of them, we also study the respective limits as the fractional order tends to one, leading to the classical third order in time (J)MGT equation.

preprint2021arXiv

All-at-once formulation meets the Bayesian approach: A study of two prototypical linear inverse problems

In this work, the Bayesian approach to inverse problems is formulated in an all-at-once setting. The advantages of the all-at-once formulation are known to include the avoidance of a parameter-to-state map as well as numerical improvements, especially when considering nonlinear problems. In the Bayesian approach, prior knowledge is taken into account with the help of a prior distribution. In addition, the error in the observation equation is formulated by means of a distribution. This method naturally results in a whole posterior distribution for the unknown target, not just point estimates. This allows for further statistical analysis including the computation of credible intervals. We combine the Bayesian setting with the all-at-once formulation, resulting in a novel approach for investigating inverse problems. With this combination we are able to chose a prior not only for the parameter, but also for the state variable, which directly influences the parameter. Furthermore, errors not only in the observation equation, but additionally, in the model can be taken into account. %The aim of this approach is not only to accomplish reasonable reconstructions of the unknown parameter but also to maximize the information gained from measurements through combining it with prior knowledge, obtained either from certain expertise or former investigation in the model. We analyze this approach with the help of two linear standard examples, namely the inverse source problem for the Poisson equation and the backwards heat equation, i.e. a stationary and a time dependent problem. Appropriate function spaces and derivation of adjoint operators are investigated. To assess the degree of ill-posedness, we analyze the singular values of the corresponding all-at-once forward operators. %as well as the convergence of the method. Finally, joint priors are designed and numerically tested.

preprint2021arXiv

Determining kernels in linear viscoelasticity

In this work, we investigate the inverse problem of determining the kernel functions that best describe the mechanical behavior of a complex medium modeled by a general nonlocal viscoelastic wave equation. To this end, we minimize a tracking-type data misfit function under this PDE constraint. We perform the well-posedness analysis of the state and adjoint problems and, using these results, rigorously derive the first-order sensitivities. Numerical experiments in a three-dimensional setting illustrate the method.

preprint2021arXiv

Determining the nonlinearity in an acoustic wave equation

We consider an undetermined coefficient inverse problem for a nonlinear partial differential equation describing high intensity ultrasound propagation as widely used in medical imaging and therapy. The usual nonlinear term in the standard model using the Westervelt equation in pressure formulation is of the form $p p_t$. However, this should be considered as a low order approximation to a more complex physical model where higher order terms will be required. Here we assume a more general case where the form taken is $f(p)\,p_t$ and $f$ is unknown and must be recovered from data measurements. Corresponding to the typical measurement setup, the overposed data consists of time trace observations of the acoustic pressure at a single point or on a one dimensional set $Σ$ representing the receiving transducer array at a fixed time. Additionally to an analysis of well-posedness of the resulting {\sc pde}, we show injectivity of the linearized forward map from $f$ to the overposed data and use this as motivation for several iterative schemes to recover $f$. Numerical simulations will also be shown to illustrate the efficiency of the methods.

preprint2021arXiv

Iterative regularization for constrained minimization formulations of nonlinear inverse problems

In this paper we the formulation of inverse problems as constrained minimization problems and their iterative solution by gradient or Newton type. We carry out a convergence analysis in the sense of regularization methods and discuss applicability to the problem of identifying the spatially varying diffusivity in an elliptic PDE from different sets of observations. Among these is a novel hybrid imaging techology known as impedance acoustic tomography, for which we provide numerical experiments.

preprint2021arXiv

On the identification of the nonlinearity parameter in the Westervelt equation from boundary measurements

We consider an undetermined coefficient inverse problem for a non-\\linear partial differential equation occurring in high intensity ultrasound propagation as used in acoustic tomography. In particular, we investigate the recovery of the nonlinearity coefficient commonly labeled as $B/A$ in the literature which is part of a space dependent coefficient $κ$ in the Westervelt equation governing nonlinear acoustics. Corresponding to the typical measurement setup, the overposed data consists of time trace measurements on some zero or one dimensional set $Σ$ representing the receiving transducer array. After an analysis of the map from $κ$ to the overposed data, we show injectivity of its linearisation and use this as motivation for several iterative schemes to recover $κ$. Numerical simulations will also be shown to illustrate the efficiency of the methods.

preprint2021arXiv

On the simultaneous recovery of the conductivity and the nonlinear reaction term in a parabolic equation

This paper considers the inverse problem of recovering both the unknown, spatially-dependent conductivity $a(x)$ and the nonlinear reaction term $f(u)$ in a reaction-diffusion equation from overposed data. These measurements can consist of: the value of two different solution measurements taken at a later time $T$; time-trace profiles from two solutions; or both final time and time-trace measurements from a single forwards solve data run. We prove both uniqueness results and the convergence of iteration schemes designed to recover these coefficients. The last section of the paper shows numerical reconstructions based on these algorithms.

preprint2021arXiv

On uniqueness and reconstruction of a nonlinear diffusion term in a parabolic equation

The problem of recovering coefficients in a diffusion equation is one of the basic inverse problems. Perhaps the most important term is the one that couples the length and time scales and is often referred to as {\it the\/} diffusion coefficient $a$ in $u_t - \nabla(a\nabla u) = f$. In this paper we seek the unknown $a$ assuming that $a=a(u)$ depends only on the value of the solution at a given point. Such diffusion models are the basic of a wide range of physical phenomena such as nonlinear heat conduction, chemical mixing and population dynamics. We shall look at two types of overposed data in order to effect recovery of $a(u)$: the value of a time trace $u(x_0,t)$ for some fixed point $x_0$ on the boundary of the region $Ω$; or the value of $u$ on an interior curve $Σ$ lying within $Ω$. As examples, these might represent a temperature measurement on the boundary or a census of the population in some subset of $Ω$ taken at a fixed time $T>0$. In the latter case we shall show a uniqueness result that leads to a constructive method for recovery of $a$. Indeed, for both types of measured data we shall show reconstructions based on the iterative algorithms developed in the paper.

preprint2020arXiv

Some application examples of minimization based formulations of inverse problems and their regularization

In this paper we extend a recent idea of formulating and regularizing inverse problems as minimization problems, so without using a forward operator, thus avoiding explicit evaluation of a parameter-to-state map. We do so by rephrasing three application examples in this minimization form, namely (a) electrical impedance tomography with the complete electrode model (b) identification of a nonlinear magnetic permeability from magnetic flux measurements (c) localization of sound sources from microphone array measurements. To establish convergence of the proposed regularization approach for these problems, we first of all extend the existing theory. In particular, we take advantage of the fact that observations are finite dimensional here, so that inversion of the noisy data can to some extent be done separately, using a right inverse of the observation operator. This new approach is actually applicable to a wide range of real world problems.

preprint2020arXiv

The Inverse Problem of Reconstructing Reaction-Diffusion Systems

This paper considers the inverse problem of recovering state-dependent source terms in a reaction-diffusion system from overposed data consisting of the values of the state variables either at a fixed finite time (census-type data) or a time trace of their values at a fixed point on the boundary of the spatial domain. We show both uniqueness results and the convergence of an iteration scheme designed to recover these sources. This leads to a reconstructive method and we shall demonstrate its effectiveness by several illustrative examples.

preprint2019arXiv

On the identification of a nonlinear term in a reaction-diffusion equation

Reaction-diffusion equations are one of the most common partial differential equations used to model physical phenomenon. They arise as the combination of two physical processes: a driving force $f(u)$ that depends on the state variable $u$ and a diffusive mechanism that spreads this effect over a spatial domain. The canonical form is $u_t - \triangle u = f(u)$. Application areas include chemical processes, heat flow models and population dynamics. The direct or forwards problem for such equations is now very well-developed and understood. However, our interest lies in the inverse problem of recovering the reaction term $f(u)$ not just at the level of determining a few parameters in a known functional form, but recovering the complete functional form itself. To achieve this we set up the usual paradigm for the parabolic equation where $u$ is subject to both given initial and boundary data, then prescribe overposed data consisting of the solution at a later time $T$. For example, in the case of a population model this amounts to census data at a fixed time. Our approach will be two-fold.First we will transform the inverse problem into an equivalent nonlinear mapping from which we seek a fixed point. We will be able to prove important features of this map such as a self-mapping property and give conditions under which it is contractive. Second, we consider the direct map from $f$ through the partial differential operator to the overposed data. We will investigate Newton schemes for this case. In recent decades various anomalous processes have been used to generalize classical Brownian diffusion. Amongst the most popular is one that replaces the usual time derivative by a fractional one of order $α\leq 1$. We will also include this model in our analysis. The final section of the paper shows numerical reconstructions that demonstrate the viability of the suggested approaches.