On the uniqueness and reconstruction for an inverse problem of the fractional diffusion process

Consider an inverse problem for the time-fractional diffusion equation in one dimensional spatial space. The aim is to determine the initial status and heat flux on the boundary simultaneously from heat measurement data given on the other boundary. Using the Laplace transform and the unique extension technique, the uniqueness for this inverse problem is proven. Then we construct a regularizing scheme for the reconstruction of boundary flux for known initial status. The convergence rate of the regularizing solution is established under some a priori information about the exact solution. Moreover, the initial distribution can also be recovered approximately from our regularizing scheme. Finally we present some numerical examples, which show the validity of the proposed reconstruction scheme.     

An Algorithm for Cavity Reconstruction in Electrical Impedance Tomography

We consider the inverse problem of finding cavities within some object from electrostatic measurements on the boundary. By a cavity we understand any object with a different electrical conductivity from the background material of the body. We give an algorithm for solving this inverse problem based on the output nonlinear least-square formulation and the regularized Newton-type iteration. In particular, we present a number of numerical results to highlight the potential and the limitations of this method.     

An inverse backward problem for degenerate parabolic equations

This work studies the inverse problem of reconstructing an initial value function in the degenerate parabolic equation using the final measurement data. Problems of this type have important applications in the field of financial engineering. Being different from other inverse backward parabolic problems, the mathematical model in our article may be allowed to degenerate at some part of boundaries, which may lead to the corresponding boundary conditions missing. The conditional stability of the solution is obtained using the logarithmic convexity method. A finite difference scheme is constructed to solve the direct problem and the corresponding stability and convergence are proved. The Landweber iteration algorithm is applied to the inverse problem and some typical numerical experiments are also performed in the paper. The numerical results show that the proposed method is stable and the unknown initial value is recovered very well.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1900–1923, 2017     

A numerical way of solving the inverse problem for the wave equation in a medium with local inhomogeneity

The inverse problem of determining the boundary of local inhomogeneity for measuring a field in a bounded receivers location domain in a three-dimensional medium is considered for the wave equation. The problem is reduced to a system of integral equations. An iteration approach to solving the inverse problem is proposed, and the results from numerical experiments are presented.     

Numerical method for determining the inhomogeneity boundary in the Dirichlet problem for Laplace’s equation in a piecewise homogeneous medium

The Dirichlet problem for Laplace’s equation in a two-dimensional domain filled with a piecewise homogeneous medium is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem of determining the inhomogeneity boundary from additional information on the solution of the Dirichlet problem is considered. A numerical method based on the linearization of the nonlinear operator equation for the unknown boundary is proposed for solving the inverse problem. The results of numerical experiments are presented.     

On numerical implementations of a new iterative method with boundary condition splitting for solving the nonstationary stokes problem in a strip with periodicity condition

Based on finite-difference approximations in time and a bilinear finite-element approximation in spatial variables, numerical implementations of a new iterative method with boundary condition splitting are constructed for solving the Dirichlet initial-boundary value problem for the nonstationary Stokes system. The problem is considered in a strip with a periodicity condition along it. At each iteration step of the method, the original problem splits into two much simpler boundary value problems that can be stably numerically approximated. As a result, this approach can be used to construct new effective and stable numerical methods for solving the nonstationary Stokes problem. The velocity and pressure are approximated by identical bilinear finite elements, and there is no need to satisfy the well-known difficult-to-verify Ladyzhenskaya-Brezzi-Babuska condition, as is usually required when the problem is discretized as a whole. Numerical iterative methods are constructed that are first- and second-order accurate in time and second-order accurate in space in the max norm for both velocity and pressure. The numerical methods have fairly high convergence rates corresponding to those of the original iterative method at the differential level (the error decreases approximately 7 times per iteration step). Numerical results are presented that illustrate the capabilities of the methods developed.     

Isospectral flow method for nonnegative inverse eigenvalue problem with prescribed structure

The nonnegative inverse eigenvalue problem is that given a family of complex numbers λ={λ₁,…,λ_n}, find a nonnegative matrix of order n with spectrum λ. This problem is difficult and remains unsolved partially. In this paper, we focus on its generalization that the reconstructed nonnegative matrices should have some prescribed entries. It is easy to see that this new problem will come back to the common nonnegative inverse eigenvalue problem if there is no constraint of the locations of entries. A numerical isospectral flow method which is developed by hybridizing the optimization theory and steepest descent method is used to study the reconstruction. Moreover, an error estimate of the numerical iteration for ordinary differential equations on the matrix manifold is presented. After that, a numerical method for the nonnegative symmetric inverse eigenvalue problem with prescribed entries and its error estimate are considered. Finally, the approaches are verified by the numerical test results.     

A method for finding coefficients of a quasilinear hyperbolic equation

The inverse problem of finding the coefficients q(s) and p(s) in the equation u _tt = a ² u _xx + q(u)u _t ? p(u)u _x is investigated. As overdetermination required in the inverse setting, two additional conditions are set: a boundary condition and a condition with a fixed value of the timelike variable. An iteration method for solving the inverse problem is proposed based on an equivalent system of integral equations of the second kind. A uniqueness theorem and an existence theorem in a small domain are proved for the inverse problem to substantiate the convergence of the algorithm.     

Numerical solution of the inverse electrocardiography problem with the use of the Tikhonov regularization method

The inverse electrocardiography problem related to medical diagnostics is considered in terms of potentials. Within the framework of the quasi-stationary model of the electric field of the heart, the solution of the problem is reduced to the solution of the Cauchy problem for the Laplace equation in R ³. A numerical algorithm based on the Tikhonov regularization method is proposed for the solution of this problem. The Cauchy problem for the Laplace equation is reduced to an operator equation of the first kind, which is solved via minimization of the Tikhonov functional with the regularization parameter chosen according to the discrepancy principle. In addition, an algorithm based on numerical solution of the corresponding Euler equation is proposed for minimization of the Tikhonov functional. The Euler equation is solved using an iteration method that involves solution of mixed boundary value problems for the Laplace equation. An individual mixed problem is solved by means of the method of boundary integral equations of the potential theory. In the study, the inverse electrocardiography problem is solved in region Ω close to the real geometry of the torso and heart.     

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 numerical method for determining the localized initial condition in the FitzHugh-Nagumo and Aliev-Panfilov models

The inverse problem for the FitzHugh-Nagumo and Aliev-Panfilov models describing wave propagation in excitable media is considered. The problem lies in determining a localized initial condition from measurements on the external boundary of a plane region. A numerical method for solving the inverse problem is proposed, and the results from a numerical solution of the inverse problem for regions similar to different sections of a heart are presented.     

Numerical methods for determining the inhomogeneity boundary in a boundary value problem for Laplace’s equation in a piecewise homogeneous medium

A boundary value problem for Laplace’s equation in a bounded two-dimensional domain filled with a piecewise homogeneous medium is considered. The boundary of the inhomogeneity is assumed to be unknown. The inverse problem of determining the inhomogeneity boundary and the solution of the equation given the solution and its normal derivative on the boundary of the domain is discussed. Numerical methods are proposed for solving the inverse problem, and the results of numerical experiments are presented.     

A full multigrid method for nonlinear eigenvalue problems

We introduce a type of full multigrid method for the nonlinear eigenvalue problem. The main idea is to transform the solution of the nonlinear eigenvalue problem into a series of solutions of the corresponding linear boundary value problems on the sequence of finite element spaces and nonlinear eigenvalue problems on the coarsest finite element space. The linearized boundary value problems are solved by some multigrid iterations. Besides the multigrid iteration, all other efficient iteration methods for solving boundary value problems can serve as the linear problem solver. We prove that the computational work of this new scheme is truly optimal, the same as solving the linear corresponding boundary value problem. In this case, this type of iteration scheme certainly improves the overfull efficiency of solving nonlinear eigenvalue problems. Some numerical experiments are presented to validate the efficiency of the new method.     

Inverse problem for the diffusion equation in the case of spherical symmetry

The initial boundary value problem for the diffusion equation is considered in the case of spherical symmetry and an unknown initial condition. Additional information used for determining the unknown initial condition is an external volume potential whose density is the Laplace operator applied to the solution of the initial boundary value problem. The uniqueness of the solution of the inverse problem is studied depending on the parameters entering into the boundary conditions. It is shown that the solution of the inverse problem is either unique or not unique up to a one-dimensional linear subspace.     

Integrability and beyond

An algebraic approach to solving nonlinear functional equations in the Riemann theta functions is stated. By the inverse scattering method and some general methods of the theory of partial differential equations, the solution of the initial boundary value problem for the nonlinear Schrödinger equation is presented. Bibliography:17 titles.     

Initial value problems for creeping flow of Maxwell fluids

We consider the flow of nonlinear Maxwell fluids in the unsteady quasistatic case, where the effect of inertia is neglected. We study the well-posedness of the resulting PDE initial-boundary value problem locally in time. This well-posedness depends on the unique solvability of an elliptic boundary value problem. We first present results for the 3D case with sufficiently small initial data and for a simple shear flow problem with arbitrary initial data; after that we extend our results to some 3D flow problems with large initial data.We solve our problem using an iteration between linear subproblems. The limit of the iteration provides the solution of our original problem.     

An inverse spectral problem for integro-differential Dirac operators with general convolution kernels

ABSTRACT

An integro-differential Dirac system with a convolution kernel consisting of four independent functions is considered. We prove that the kernel is uniquely determined by specifying the spectra of two boundary value problems with one common boundary condition. The proof is based on the reduction of this nonlinear inverse problem to solving some nonlinear integral equation, which we solve globally. On this basis we also obtain a constructive procedure for solving the inverse problem along with necessary and sufficient conditions for its solvability in an appropriate class of kernels.     

Reconstruction of boundary controls in parabolic systems

In the paper, an inverse dynamic problem is considered. It consists in reconstructing a priori unknown boundary controls in dynamical systems described by boundary value problems for partial differential equations of parabolic type. The source information for solving the inverse problem is the results of approximate measurements of the states of the observed system’s motion. The problem is solved in the static case; i.e., to solve it, we use all the measurement data accumulated during some specified observation interval. The problem under consideration is ill-posed. To solve it, we propose the Tikhonov method with a stabilizer containing the sum of the mean-square norm and total time variation of the control. The use of such nondifferentiable stabilizer allows us to obtain more precise results than the approximation of the desired control in the Lebesgue spaces. In particular, this method provides the pointwise and piecewise uniform convergences of regularized approximations and makes possible the numerical reconstruction of the subtle structure of the desired control. In the paper, the subgradient projection method for obtaining a minimizing sequence for the Tikhonov functional is described and substantiated. Also, we demonstrate the two-stage finitedimensional approximation of the problem and present the results of numerical simulation.     

On the determination of the coefficients in the viscoelasticity equations

For the integrodifferential viscoelasticity equations, we study the problem of determining the coefficients of the equations and the kernels occurring in the integral terms of the system of equations. The density of the medium is assumed to be given. We suppose that the inhomogeneity support of the sought functions is included in some compact domain B ₀. We consider a series of inverse problems in which an impulse source is concentrated at the points y of the boundary of B ₀. The point y is the parameter of the problem. The given information about the solution is the trace of the solution to the Cauchy problem with zero initial data. This trace is given on the boundary of B ₀ for all y ∈ ?B ₀ and for a finite time interval. The main result of the article consists in obtaining uniqueness theorems for a solution to the initial inverse problem.     

Uniqueness and nonuniqueness of the solution to the problem of determining the source in the heat equation

An initial–boundary value problem for the two-dimensional heat equation with a source is considered. The source is the sum of two unknown functions of spatial variables multiplied by exponentially decaying functions of time. The inverse problem is stated of determining two unknown functions of spatial variables from additional information on the solution of the initial–boundary value problem, which is a function of time and one of the spatial variables. It is shown that, in the general case, this inverse problem has an infinite set of solutions. It is proved that the solution of the inverse problem is unique in the class of sufficiently smooth compactly supported functions such that the supports of the unknown functions do not intersect. This result is extended to the case of a source involving an arbitrary finite number of unknown functions of spatial variables multiplied by exponentially decaying functions of time.