Identification of a dissipation coefficient by a variational method

A generalized inverse problem for the identification of the absorption coefficient for a hyperbolic system is considered. The well-posedness of the problem is examined. It is proved that the regular part of the solution is an L ₂ function, which reduces the inverse problem to minimizing the error functional. The gradient of the functional is determined in explicit form from the adjoint problem, and approximate formulas for its calculation are derived. A regularization algorithm for the solution of the inverse problem is considered. Numerical results obtained for various excitation sources are displayed.     

Existence and uniqueness for a nonlinear inverse reaction‐diffusion problem with a nonlinear source in higher dimensions

This paper analyzes the existence and the uniqueness problem for an n‐dimensional nonlinear inverse reaction‐diffusion problem with a nonlinear source. A transformation is used to obtain a new inverse coefficient problem. Then, a parabolic differential operator L_λ is defined to establish the relation between the solution of L_λ = 0 and the new inverse problem. Following this, it is shown that the inverse problem has at least one solution in the class of admissible coefficients. Furthermore, it is proved that this solution is the unique solution of the undertaken inverse problem. A numerical example is given to illustrate ill‐posedness of the inverse problem. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical method for solving an inverse thermoacoustic problem

In this paper, we consider an inverse problem of determining the initial condition of an initial boundary value problem for the wave equation with some additional information about solving a direct initial boundary value problem. The information is obtained from measurements at the boundary of the solution domain. The purpose of our paper is to construct a numerical algorithm for solving the inverse problem by an iterative method called a method of simple iteration (MSI) and to study the resolution quality of the inverse problem as a function of the number and location of measurement points. Three two-dimensional inverse problem formulations are considered. The results of our numerical calculations are presented. It is shown that the MSI decreases the objective functional at each iteration step. However, due to the ill-posedness of the inverse problem the difference between the exact and approximate solutions decreases up to some fixed number k _min, and then monotonically increases. This shows the regularizing properties of the MSI, and the iteration number can be considered a regularization parameter.     

The Inverse Problem of Simultaneous Determination of the Two Time-Dependent Lower Coefficients in a Nondivergent Parabolic Equation in the Plane

We prove existence and uniqueness theorems for the solution of the inverse problem of simultaneous determination of the t-dependent coefficients of u and u_x in a nondivergent parabolic equation with two independent variables from integral observation of x. Estimates of the maxima of the moduli of these coefficients with constants explicitly expressed in terms of the input data of the problem are given. An example of an inverse problem to which the proved theorems apply is presented.

     

On a backward problem for fractional diffusion equation with Riemann-Liouville derivative

In the present paper, we study the initial inverse problem (backward problem) for a two-dimensional fractional differential equation with Riemann-Liouville derivative. Our model is considered in the random noise of the given data. We show that our problem is not well-posed in the sense of Hadamard. A truncated method is used to construct an approximate function for the solution (called the regularized solution). Furthermore, the error estimate of the regularized solution in L² and H^τ norms is considered and illustrated by numerical example.     

A Tikhonov-based projection iteration for nonlinear Ill-posed problems with sparsity constraints

In this paper, we consider nonlinear inverse problems where the solution is assumed to have a sparse expansion with respect to a preassigned basis or frame. We develop a scheme which allows to minimize a Tikhonov functional where the usual quadratic regularization term is replaced by a one-homogeneous (typically weighted ℓ _p ) penalty on the coefficients (or isometrically transformed coefficients) of such expansions. For (p < 2), the regularized solution will have a sparser expansion with respect to the basis or frame under consideration. The computation of the regularized solution amounts in our setting to a Landweber-fixed-point iteration with a projection applied in each fixed-point iteration step. The performance of the resulting numerical scheme is demonstrated by solving the nonlinear inverse single photon emission computerized tomography (SPECT) problem.     

A Globally Convergent Algorithm for a PDE-Constrained Optimization Problem Arising in Electrical Impedance Tomography

We study the convergence properties of an algorithm for the inverse problem of electrical impedance tomography, which can be reduced to a partial differential equation (PDE) constrained optimization problem. The direct problem consists of the potential equation div(??u) = 0 in a circle, with Neumann condition describing the behavior of the electrostatic potential in a medium with conductivity given by the function ?(x, y). We suppose that at each time a current ψ _i is applied to the boundary of the circle (Neumann's data), and that it is possible to measure the corresponding potential ? _i (Dirichlet data). The inverse problem is to find ?(x, y) given a finite number of Cauchy pairs measurements (? _i , ψ _i ), i = 1,…, N. The problem is formulated as a least squares problem, and the developed algorithm solves the continuous problem using descent iterations in its corresponding finite element approximations. Wolfe's conditions are used to ensure the global convergence of the optimization algorithm for the continuous problem. Although exact data are assumed, measurement errors in data and regularization methods shall be considered in a future work.     

An inverse source problem for a two dimensional time fractional diffusion equation with nonlocal boundary conditions

We consider the inverse source problem for a time fractional diffusion equation. The unknown source term is independent of the time variable, and the problem is considered in two dimensions. A biorthogonal system of functions consisting of two Riesz bases of the space L²[(0,1) × (0,1)], obtained from eigenfunctions and associated functions of the spectral problem and its adjoint problem, is used to represent the solution of the inverse problem. Using the properties of the biorthogonal system of functions, we show the existence and uniqueness of the solution of the inverse problem and its continuous dependence on the data. Copyright © 2012 John Wiley & Sons, Ltd.     

A conformal mapping algorithm for the Bernoulli free boundary value problem

We propose a new numerical method for the solution of the Bernoulli free boundary value problem for harmonic functions in a doubly connected domain D in where an unknown free boundary Γ₀ is determined by prescribed Cauchy data on Γ₀ in addition to a Dirichlet condition on the known boundary Γ₁. Our main idea is to involve the conformal mapping method as proposed and analyzed by Akduman, Haddar, and Kress for the solution of a related inverse boundary value problem. For this, we interpret the free boundary Γ₀ as the unknown boundary in the inverse problem to construct Γ₀ from the Dirichlet condition on Γ₀ and Cauchy data on the known boundary Γ₁. Our method for the Bernoulli problem iterates on the missing normal derivative on Γ₁ by alternating between the application of the conformal mapping method for the inverse problem and solving a mixed Dirichlet–Neumann boundary value problem in D. We present the mathematical foundations of our algorithm and prove a convergence result. Some numerical examples will serve as proof of concept of our approach. Copyright © 2015 John Wiley & Sons, Ltd.     

A Richardson-type iterative approach for identification of delamination boundaries

A direct problem of applied mathematical modelling is to determine the response of a system given the governing partial differential equations, the geometry of interest, the complete boundary and initial conditions, and material properties. When one or more of the conditions for the solution of the direct problem are unknown, an inverse problem can be formulated. One of the methods frequently used for the solution of inverse problems involves finding the values of the unknowns in a mathematical formulation such that the behavior calculated with the model matches the measured response to a degree evaluated in terms of the classical L ₂ norm. Considered in this sense, the inverse problem is equivalent to an ill-posed optimization problem for the estimation of parameters whose solution in the majority of cases is a real mathematical challenge. In this contribution, we report a novel approach that avoids the mathematical difficulties inspired by the ill-posed character of the model. Our method is devoted to the computation of inverse problems furnished by second-order elliptical systems of partial differential equations and falls in the same conceptual line with the method initiated by Kozlov et al. and further extended and algorithmized by Weikl et al. We construct and employ a weak version of the algorithm found by Weikl et al. Proofs for the convergence and regularity of this version are given for the case of a single layer. The computational realization of the algorithm (called briefly AICRA) is applied and numerical results are obtained. Comparison with experiments demonstrates a good significance and representativeness. __________ Translated from Fundamentalnaya i Prikladnaya Matematika, Vol. 12, No. 4, pp. 209–230, 2006.     

Unique solvability of the inverse problem of determination of the leading coefficient in a parabolic equation

We study existence and uniqueness of the solution for the inverse problem of determination of the unknown coefficient ϱ(t) multiplying u _t in a nondivergence parabolic equation. As additional information, the integral of the solution over the domain of space variables with some given weight function is specified. The coefficients of the equation depend both on time and on the space variables.     

A variational method for the solution of biharmonic problems for a rectangular domain

We develop a variational method for the solution of biharmonic problems for a rectangular domain where, at one pair of its opposite sides, the unknown function and its normal derivative take zero values, and, at the other pair, certain inhomogeneous conditions are valid. The cases of semiinfinite and finite domain are considered. The method is based on the minimization of a quadratic functional determining the deviation of the solution from the given inhomogeneous conditions in the norm of L ₂. To solve this variational problem, we apply the expansion of the solution in the systems of complex biharmonic functions (the so-called Papkovich homogeneous solutions), which satisfy identically the given homogeneous conditions at the pair of opposite sides of the rectangle. This representation of the solution is somewhat different from that proposed earlier [V. F. Chekurin, “A variational method for the solution of direct and inverse problems of the theory of elasticity for a semiinfinite strip,” Izv. Ross. Akad. Nauk, Mekh. Tverdogo Tela, No. 2, 58–70 (1999)]. We consider several variants of inhomogeneous boundary conditions arising in the problems of the two-dimensional theory of elasticity. Finally, we give an example of applying the proposed method for the determination of stress distributions in a rectangular area one of whose sides is rigidly fastened and the opposite one is subjected to the action of normal forces. Translated from Matematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 51, No. 1, pp. 88–98, January–March, 2008.     

Application of Rothe’s Method to a Parabolic Inverse Problem with Nonlocal Boundary Condition

We consider an inverse problem for the determination of a purely time-dependent source in a semilinear parabolic equation with a nonlocal boundary condition. An approximation scheme for the solution together with the well-posedness of the problem with the initial value u₀ ∈ H¹(Ω) is presented by means of the Rothe time-discretization method. Further approximation scheme via Rothe’s method is constructed for the problem when u₀ ∈ L²(Ω) and the integral kernel in the nonlocal boundary condition is symmetric.

     

A simultaneous approach to inverse source problem by Green's function

In this paper, we consider an inverse source problem of identification of F(t) function in the linear parabolic equation u_t = u_xx + F(t) and u₀(x) function as the initial condition from the measured final data and local boundary data. Based on the optimal control framework by Green's function, we construct Fréchet derivative of Tikhonov functional. The stability of the minimizer is established from the necessary condition. The CG algorithm based on the Fréchet derivative is applied to the inverse problem, and results are presented for a test example. Copyright © 2014 John Wiley & Sons, Ltd.     

On the Neumann problem for the Helmholtz equation in a plane angle

We prove that the solution of the Neumann problem for the Helmholtz equation in a plane angle Ω with boundary conditions from the space H_−1/2(Γ), where Γ is the boundary of Ω, which is provided by the well‐known Sommerfeld integral, belongs to the Sobolev space H₁(Ω) and depends continuously on the boundary values. To this end, we use another representation of the solution given by the inverse two‐dimensional Fourier transform of an analytic function depending on the Cauchy data of the solution. Copyright © 2000 John Wiley & Sons, Ltd.     

Inverse problem of determining the kernel in an integro-differential equation of parabolic type

We study the inverse problem of determining the multidimensional kernel of the integral term in a parabolic equation of second order. As additional information, the solution of the direct problem is given on the hyperplane x _n = 0. We prove a local existence and uniqueness theorem for the inverse problem.     

Marching schemes for inverse acoustic scattering problems

Summary For the numerical solution of inverse Helmholtz problems the boundary value problem for a Helmholtz equation with spatially variable wave number has to be solved repeatedly. For large wave numbers this is a challenge. In the paper we reformulate the inverse problem as an initial value problem, and describe a marching scheme for the numerical computation that needs only n² log n operations on an n × n grid. We derive stability and error estimates for the marching scheme. We show that the marching solution is close to the low-pass filtered true solution. We present numerical examples that demonstrate the efficacy of the marching scheme.     

On the decidability of the word problem for amalgamated free products of inverse semigroups

We study inverse semigroup amalgams [S ₁,S ₂;U], where S ₁ and S ₂ are finitely presented inverse semigroups with decidable word problem and U is an inverse semigroup with decidable membership problem in S ₁ and S ₂. We use a modified version of Bennett’s work on the structure of Schützenberger graphs of the ℛ-classes of S ₁* _U S ₂ to state sufficient conditions for the amalgamated free products S ₁* _U S ₂ having decidable word problem.     

On the one-dimensional two-phase inverse Stefan problems

New formulations of the inverse nonstationary Stefan problems are considered: (a) forx [0,1] (the inverse problem IP₁; (b) forx [0, (t)] with a degenerate initial condition (the inverse problem IP). Necessary conditions for the existence and uniqueness of a solution to these problems are formulated. On the first phase {x [0, y(t)]{, the solution of the inverse problem is found in the form of a series; on the second phase {x [y(t), 1] orx [y(t), (t)]{, it is found as a sum of heat double-layer potentials. By representing the inverse problem in the form of two connected boundary-value problems for the heat conduction equation in the domains with moving boundaries, it can be reduced to the integral Volterra equations of the second kind. An exact solution of the problem IP is found for the self similar motion of the boundariesx=y(t) andx=(t).Translated from Ukrainskii Matematicheskii Zhurnal, Vol. 45, No. 8, pp. 1058–1065, August, 1993.