BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Barbara Kaltenbacher

Barbara Kaltenbacher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAmath.APNumerical Analysismath.OC

Catalog footprint

What is connected

30works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

30 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Integration based profile likelihood calculation for PDE constrained parameter estimation problems

Partial differential equation (PDE) models are widely used in engineering and natural sciences to describe spatio-temporal processes. The parameters of the considered processes are often unknown and have to be estimated from experimental data. Due to partial observations and measurement noise, these parameter estimates are subject to uncertainty. This uncertainty can be assessed using profile likelihoods, a reliable but computationally intensive approach. In this paper, we introduce an integration based approach for the profile likelihood calculation for inverse problems with PDE constraints. While existing approaches rely on repeated optimization, the proposed approach exploits a dynamical system evolving along the likelihood profile. We derive the dynamical system for the reduced and the full estimation problem and study its properties. To evaluate the proposed method, we compare it with state-of-the-art algorithms for a simple reaction-diffusion model for a cellular patterning process. We observe a good accuracy of the method as well as a significant speed up as compared to established methods. Integration based profile calculation facilitates rigorous uncertainty analysis for computationally demanding parameter estimation problems with PDE constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lavrentiev's regularization method in Hilbert spaces revisited

In this paper, we deal with nonlinear ill-posed problems involving monotone operators and consider Lavrentiev's regularization method. This approach, in contrast to Tikhonov's regularization method, does not make use of the adjoint of the derivative. There are plenty of qualitative and quantitative convergence results in the literature, both in Hilbert and Banach spaces. Our aim here is mainly to contribute to convergence rates results in Hilbert spaces based on some types of error estimates derived under various source conditions and to interpret them in some settings. In particular, we propose and investigate new variational source conditions adapted to these Lavrentiev-type techniques. Another focus of this paper is to exploit the concept of approximate source conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Regularization based on all-at-once formulations for inverse problems

Parameter identification problems typically consist of a model equation, e.g. a (system of) ordinary or partial differential equation(s), and the observation equation. In the conventional reduced setting, the model equation is eliminated via the parameter-to-state map. Alternatively, one might consider both sets of equations (model and observations) as one large system, to which some regularization method is applied. The choice of the formulation (reduced or all-at-once) can make a large difference computationally, depending on which regularization method is used: Whereas almost the same optimality system arises for the reduced and the all-at-once Tikhonov method, the situation is different for iterative methods, especially in the context of nonlinear models. In this paper we will exemplarily provide some convergence results for all-at-once versions of variational, Newton type and gradient based regularization methods. Moreover we will compare the implementation requirements for the respective all-at-one and reduced versions and provide some numerical comparison.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Regularization by Discretization in Banach Spaces

We consider ill-posed linear operator equations with operators acting between Banach spaces. For solution approximation, the methods of choice here are projection methods onto finite dimensional subspaces, thus extending existing results from Hilbert space settings. More precisely, general projection methods, the least squares method and the least error method are analyzed. In order to appropriately choose the dimension of the subspace, we consider a priori and a posteriori choices by the discrepancy principle and by the monotone error rule. Analytical considerations and numerical tests are provided for a collocation method applied to a Volterra integral equation in one space dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The least error method for sparse solution reconstruction

This work deals with a regularization method enforcing solution sparsity of linear ill-posed problems by appropriate discretization in the image space. Namely, we formulate the so called least error method in an $\ell^1$ setting and perform the convergence analysis by choosing the discretization level according to an a priori rule, as well as two a posteriori rules, via the discrepancy principle and the monotone error rule, respectively. Depending on the setting, linear or sublinear convergence rates in the $\ell^1$-norm are obtained under a source condition yielding sparsity of the solution. A part of the study is devoted to analyzing the structure of the approximate solutions and of the involved source elements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An online parameter identification method for time dependent partial differential equations

Online parameter identification is of importance, e.g., for model predictive control. Since the parameters have to be identified simultaneously to the process of the modeled system, dynamical update laws are used for state and parameter estimates. Most of the existing methods for infinite dimensional systems either impose strong assumptions on the model or cannot handle partial observations. Therefore we propose and analyze an online parameter identification method that is less restrictive concerning the underlying model and allows for partial observations and noisy data. The performance of our approach is illustrated by some numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computing quasisolutions of nonlinear inverse problems via efficient minimization of trust region problems

In this paper we present a method for the regularized solution of nonlinear inverse problems, based on Ivanov regularization (also called method of quasi solutions or constrained least squares regularization). This leads to the minimization of a non-convex cost function under a norm constraint, where non-convexity is caused by nonlinearity of the inverse problem. Minimization is done by iterative approximation, using (non-convex) quadratic Taylor expansions of the cost function. This leads to repeated solution of quadratic trust region subproblems with possibly indefinite Hessian. Thus the key step of the method consists in application of an efficient method for solving such quadratic subproblems, developed by Rendl and Wolkowicz [10]. We here present a convergence analysis of the overall method as well as numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On higher regularity for the Westervelt equation with strong nonlinear damping

We show higher interior regularity for the Westervelt equation with strong nonlinear damping term of the $q$-Laplace type. Secondly, we investigate an interface coupling problem for these models, which arise, e.g., in the context of medical applications of high intensity focused ultrasound in the treatment of kidney stones. We show that the solution to the coupled problem exhibits piecewise $H^2$ regularity in space, provided that the gradient of the acoustic pressure is essentially bounded in space and time on the whole domain. This result is of importance in numerical approximations of the present problem, as well as in gradient based algorithms for finding the optimal shape of the focusing acoustic lens in lithotripsy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sensitivity analysis for shape optimization of a focusing acoustic lens in lithotripsy

We are interested in shape sensitivity analysis for an optimization problem arising in medical applications of high intensity focused ultrasound. The goal is to find the optimal shape of a focusing acoustic lens so that the desired acoustic pressure at a kidney stone is achieved. Coupling of the silicone acoustic lens and nonlinearly acoustic fluid region is modeled by the Westervelt equation with nonlinear strong damping and piecewise constant coefficients. We follow the variational approach to calculating the shape derivative of the cost functional which does not require computing the shape derivative of the state variable; however assumptions of certain spatial regularity of the primal and the adjoint state are needed to obtain the derivative, in particular for its strong form according to the Delfour-Hadamard-Zol\' esio Structure Theorem.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A convergence rates result for an iteratively regularized Gauss-Newton-Halley method in Banach space

The use of second order information on the forward operator often comes at a very moderate additional computational price in the context of parameter identification probems for differential equation models. On the other hand the use of general (non-Hilbert) Banach spaces has recently found much interest due to its usefulness in many applications. This motivates us to extend the second order method previously considered by the author in a Hilbert space setting, (see also Hettlich and Rundell 2000) to a Banach space setting and analyze its convergence. We here show rates results for a particular source condition and different exponents in the formulation of Tikhonov regularization in each step. This includes a complementary result on the (first order) iteratively regularized Gauss-Newton method (IRGNM) in case of a one-homogeneous data misfit term, which corresponds to exact penalization. The results clearly show the possible advantages of using second order information, which get most pronounced in this exact penalization case. Numerical simulations for a coefficient identification problem in an elliptic PDE illustrate the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A thermodynamically consistent phenomenological model for ferroelectric and ferroelastic hysteresis

We propose a hysteretic model for electromechanical coupling in piezoelectric materials, with the strain and the electric field as inputs and the stress and the polarization as outputs. This constitutive law satisfies the thermodynamic principles and exhibits good agreement with experimental measurements. Moreover, when it is coupled with the mechanical and electromagnetic balance equations, the resulting PDE system is well-posed under the hypothesis that hysteretic effects take place only in one preferred direction. We prove the existence and uniqueness of its global weak solutions for each initial data with prescribed regularity. One of the tools is a new Lipschitz continuity theorem for the inverse Preisach operator with time dependent coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Absorbing boundary conditions for the Westervelt equation

The focus of this work is on the construction of a family of nonlinear absorbing boundary conditions for the Westervelt equation in one and two space dimensions. The principal ingredient used in the design of such conditions is pseudo-differential calculus. This approach enables to develop high order boundary conditions in a consistent way which are typically more accurate than their low order analogs. Under the hypothesis of small initial data, we establish local well-posedness for the Westervelt equation with the absorbing boundary conditions. The performed numerical experiments illustrate the efficiency of the proposed boundary conditions for different regimes of wave propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Enhanced choice of the parameters in an iteratively regularized Newton-Landweber iteration in Banach space

This paper is a close follow-up of Kaltenbacher and Tomba 2013 and Jin 2012, where Newton-Landweber iterations have been shown to converge either (unconditionally) without rates or (under an additional regularity assumption) with rates. The choice of the parameters in the method were different in each of these two cases. We now found a unified and more general strategy for choosing these parameters that enables both convergence and convergence rates. Moreover, as opposed to the previous one, this choice yields strong convergence as the noise level tends to zero, also in the case of no additional regularity. Additionally, the resulting method appears to be more efficient than the one from Kaltenbacher and Tomba 2013, as our numerical tests show.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Relaxation of regularity for the Westervelt equation by nonlinear damping with application in acoustic-acoustic and elastic-acoustic coupling

In this paper we show local (and partially global) in time existence for the Westervelt equation with several versions of nonlinear damping. This enables us to prove well-posedness with spatially varying $L_\infty$-coefficients, which includes the situation of interface coupling between linear and nonlinear acoustics as well as between linear elasticity and nonlinear acoustics, as relevant, e.g., in high intensity focused ultrasound (HIFU) applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Well-Posedness and asymptotic behavior of solutions for the Blackstock-Crighton-Westervelt equation

We consider a nonlinear fourth order in space partial differential equation arising in the context of the modeling of nonlinear acoustic wave propagation in thermally relaxing viscous fluids. We use the theory of operator semigroups in order to investigate the linearization of the underlying model and see that the underlying semigroup is analytic. This leads to exponential decay results for the linear homogeneous equation. Moreover, we prove local in time well-posedness of the model under the assumption that initial data are sufficiently small by employing a fixed point argument. Global in time well-posedness is obtained by performing energy estimates and using the classical barrier method, again for sufficiently small initial data. Additionally, we provide results concerning exponential decay of solutions of the nonlinear equation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Efficient time integration methods based on operator splitting and application to the Westervelt equation

Efficient time integration methods based on operator splitting are introduced for the Westervelt equation, a nonlinear damped wave equation that arises in nonlinear acoustics as mathematical model for the propagation of sound waves in high intensity ultrasound applications. For the first-order Lie-Trotter splitting method a global error estimate is deduced, confirming that the splitting method remains stable and that the nonstiff convergence order is retained in situations where the problem data are sufficiently regular. Fundamental ingredients in the stability and error analysis are regularity results for the Westervelt equation and related linear evolution equations of hyperbolic and parabolic type. Numerical examples illustrate and complement the theoretical investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Goal oriented adaptivity in the IRGNM for parameter identification in PDEs I:\\ reduced formulation

In this paper we study adaptive discretization of the iteratively regularized Gauss-Newton method IRGNM with an a posteriori (discrepancy principle) choice of the regularization parameter in each Newton step and of the stopping index. We first of all prove convergence and convergence rates under some accuracy requirements formulated in terms of four quantities of interest. Then computation of error estimators for these quantities based on a weighted dual residual method is discussed, which results in an algorithm for adaptive refinement. Finally we extend the results from the Hilbert space setting with quadratic penalty to Banach spaces and general Tikhonov functionals for the regularization of each Newton step.

0 citations0 reviewsNo code linkedOpen access